Identification of potential key genes and high-frequency mutant genes in prostate cancer by using RNA-Seq data

  • Authors:
    • Ze Zhang
    • He Wu
    • Hong Zhou
    • Yunhe Gu
    • Yufeng Bai
    • Shiliang Yu
    • Ruihua An
    • Jiping Qi
  • View Affiliations

  • Published online on: January 24, 2018     https://doi.org/10.3892/ol.2018.7846
  • Pages: 4550-4556
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

The aim of the present study was to identify potential key genes and single nucleotide variations (SNVs) in prostate cancer. RNA sequencing (RNA‑seq) data, GSE22260, were downloaded from the Gene Expression Omnibus database, including 4 prostate cancer samples and 4 normal tissues samples. RNA‑Seq reads were processed using Tophat and differentially-expressed genes (DEGs) were identified using the Cufflinks package. Gene Ontology enrichment analysis of DEGs was performed. Subsequently, Seqpos was used to identify the potential upstream regulatory elements of DEGs. SNV was analyzed using Genome Analysis Toolkit. In addition, the frequency and risk‑level of mutant genes were calculated using VarioWatch. A total of 150 upregulated and 211 downregulated DEGs were selected and 25 upregulated and 17 downregulated potential upstream regulatory elements were identified, respectively. The SNV annotations of somatic mutations revealed that 65% were base transition and 35% were base transversion. At frequencies ≥2, a total of 17 mutation sites were identified. The mutation site with the highest frequency was located in the folate hydrolase 1B (FOLH1B) gene. Furthermore, 20 high‑risk mutant genes with high frequency were identified using VarioWatch, including ribosomal protein S4 Y‑linked 2 (RPS4Y2), polycystin 1 transient receptor potential channel interacting (PKD1) and FOLH1B. In addition, kallikrein 1 (KLK1) and PKD1 are known tumor suppressor genes. The potential regulatory elements and high‑frequency mutant genes (RPS4Y2, KLK1, PKD1 and FOLH1B) may have key functions in prostate cancer. The results of the present study may provide novel information for the understanding of prostate cancer development.

Introduction

Prostate cancer is the sixth leading cause of cancer-associated mortality in males worldwide, and is the second leading cause among males in the United States (1). The morbidity rates of prostate cancer vary widely across the world and are less frequent in South and East Asia, compared with that in Europe and United States (2). In spite of the high incidence, the etiology of prostate cancer remains unknown (3). Furthermore, early prostate cancer is typically asymptomatic, with ~66% of patients diagnosed with prostate cancer exhibiting no symptoms, preventing the early control of prostate cancer (4). The identification of novel therapeutic targets is required to diagnose early prostate cancer.

Established risk factors of prostate cancer include age, race and exhibiting a family history of prostate cancer (5). Genetic factors, which may increase the risk of developing prostate cancer, are associated with ethnicity and family history of the disease (6). A number of different genes have been implicated in the development of prostate cancer (7,8). A previous study validated that prostate cancer is prone to recurrent gene fusions of androgen-regulated genes, including transmembrane protease serine 2 (TMPRSS2), solute carrier family 45 member 3 and N-myc downstream regulated 1, and E26 transformation-specific (ETS) transcription factors, including ERG ETS transcription factor (ERG) and ETS variant 1 (9). In addition, loss of tumor suppressor genes [including Kruppel-like factor 6 (10), BTG anti-proliferation factor 3 (11) and phosphatase and tensin homolog (12)] and activation of oncogenes [including nuclear receptor coactivator 2 (NCOA2) and MYC proto-oncogene bHLH transcription factor] in the carcinogenesis of prostate cancer have been identified (13). Although previous studies have revealed a number of candidate genes, the causes of prostate cancer remain unclear.

Next-generation RNA sequencing (RNA-seq) enables the identification of genes that may be susceptible to prostate cancer (14). On the basis of RNA-seq data, the transcriptome profiles of primary prostate cancer are identified, including gene fusions, long non-coding RNAs, alternative splicing and somatic mutations (9). Ren et al (15) used RNA-seq to profile genetic aberrations in Chinese patients with prostate cancer and identify three recurrent gene fusions, including TMPRSS2-ERG, ubiquitin specific peptidase 9 Y-linked-testis-specific transcripts, Y-linked 15 on chromosome Y and the interchromosomal translocation of CTAGE family member 5 ER export factor-KH RNA binding domain containing signal transduction associated 3. Kannan et al (16) investigated chimeric RNAs expressed in human prostate cancer and obtained 1.3 billion sequence reads, which enabled the identification of 2,369 chimeric RNA candidates for distinguishing of prostate cancer. On the basis of this, Xu et al (17) identified 116 disruptive mutations in 92 genes with high confidence, including a frameshift insertion/deletion in the coding region of the TNF superfamily member 10 gene associated with apoptosis. However, differentially-expressed genes (DEGs) between prostate cancer samples and normal control samples were not investigated, and a number of somatic mutations remain unknown.

In the present study, raw RNA-seq data from the study by Kannan et al (16) was downloaded from the National Center for Biotechnology Information database and analyzed by a number of bioinformatics methods. First, DEGs between prostate cancer samples and normal control samples were identified. Subsequently, functional enrichment analysis of DEGs was performed, to understand the underlying molecular mechanisms of prostate cancer. Furthermore, the upstream regulatory elements of DEGs were identified and analyzed. In addition, single nucleotide variations (SNVs) were determined and the somatic mutation sites were annotated. The risk level of SNVs was assessed for selection of SNVs with high risk and high frequency. These candidate genes may contribute to further the understanding of prostate cancer.

Materials and methods

Raw RNA-seq data

The RNA-seq data GSE22260 (16) were downloaded from the Gene Expression Omnibus (www.ncbi.nlm.nih.gov) database, including 4 prostate cancer samples (GSM554078, GSM554082, GSM554086 and GSM554088) and 4 matched normal tissues samples (GSM554118, GSM554120, GSM554122 and GSM554124), respectively. The Gleason score of the four cancer samples were 7, 7, 7 and 6 (18). The sequencing platform was Illumina Genome Analyzer IIx (Illumina, Inc., San Diego, CA, USA). Reads were generated via a paired-end approach.

Read alignment

All RNA-Seq reads were mapped to the reference human genome (hg19) of the University of California Santa Cruz (Santa Cruz, CA, USA) by using Tophat software (19). Only the reads that mapped to specific genome locations were retained. A maximum of two mismatches in each read were permitted. Other parameters were set up according to the default settings of Tophat.

DEG analysis

On the basis of Refseq gene annotation (20), transcripts were assembled using Cufflinks (version 0.9.3) (21). Subsequently, the gene expression levels of transcripts were calculated using Cuffdiff (part of the Cufflinks package), on the basis of the fragments/kilobase/million reads method (22). The differentially expressed transcripts were identified by calculating the fold change in the transcript and Student's t-tests in Cufflinks were performed to assess the difference. Transcripts with a log2 fold change >2 and P<0.05 were considered DEGs.

Function enrichment analysis of DEGs

Gene Ontology (GO) (23) enrichment analysis was performed for the aforementioned DEGs, on the basis of the Database for Annotation, Visualization and Integrated Discovery database (www.david.niaid.nih.gov) (24). GO categories were classified into biological process (BP), molecular function and cellular component (CC) GO-terms. P<0.05 was set as the threshold criterion.

DEGs which demonstrated regulatory function were selected, labeled and inputted into the tumor suppressor gene (25) database (bioinfo.mc.vanderbilt.edu/TSGene) and the tumor-associated gene (26) database (www.binfo.ncku.edu.tw/TAG/GeneDoc.php) to select cancer-related genes (tumor suppressor genes or oncogenes) for additional analysis.

Upstream regulatory elements of DEGs

In the present study, the upstream region (1.5 kb) of the transcription start site was defined as the promoter region. For the promoter region of up- and downregulated DGEs, motif identification was performed using Seqpos (27) to identify transcription factors. P<0.00001 and the frequency of motif targeted sequences >50% of up- and downregulated DEGs were set as the threshold criteria.

Identification of SNVs

The Genome Analysis Toolkit (28) software (available at www.broadinstitute.org/gatk) was used to identify the SNVs. The coverage of each credible SNV was >5×. The minimum quality score of reliability SNV was 30 and SNVs with a quality score >50 were defined as high reliability SNVs. Based on the single nucleotide polymorphism (SNP) sites documented in dbSNP137 and 1,000 genome databases, the already known SNV callings in tumor tissues were removed. To remove the interference of RNA editing in the transcriptome, SNV calling was optimized by combining with RNA-seq data of normal controls.

Somatic mutation sites annotation and high-risk mutation sites assessment

The VarioWatch (29) software (genepipe.ncgm.sinica.edu.tw/variowatch/main.do) was used to annotate the SNVs in coding sequence and assess the functional impact of gene products, on the basis of the risk assessment of the software. The SNVs with high-risk level of abnormal protein function were selected.

Results

Differential expression analysis and functional enrichment analysis

According to the differential expression analysis of paired RNA-seq data between prostate cancer samples and normal control samples, a total of 150 upregulated and 211 downregulated DEGs were selected (Table I). On the basis of the annotation information of transcription factors, 8 differentially-expressed transcription factors were identified, 4 of which were upregulated [interferon-γ inducible protein 16 (IFI16), neurogenin 3 (NEUROG3), retinoic acid receptor-γ (RARG) and single-minded family bHLH transcription factor 1 (SIM1)] and 4 were downregulated [diencephalon/mesencephalon homeobox 1 (DMBX1), nuclear receptor coactivator 2 (NCOA2), one cut homeobox 2 (ONECUT2) and zing finger protein 83 (ZNF83)] (Table I).

Table I.

Comparison between identified DEGs from prostate cancer samples and normal control samples.

Table I.

Comparison between identified DEGs from prostate cancer samples and normal control samples.

ExpressionDEGsTF countsTF genes
Downregulated2114IFI16, NEUROG3, RARG, SIM1
Upregulated1504DMBX1, NCOA2, ONECUT2, ZNF83

[i] DEGs, differentially-expressed genes; TF, transcription factor.

The most significantly enriched BP GO-terms of upregulated DEGs were chromatin assembly or disassembly (P=0.0016) and nucleosome assembly (P=0.0024) (Table II). In addition, the most significant CC GO-terms identified were nucleosome (P=0.0006) and protein-DNA complex (P=0.0020). The results of the present study suggested that genes associated with chromatin structure stability may be highly expressed in prostate cancer.

Table II.

Ten most significant GO terms of up- and downregulated DEGs in prostate cancer.

Table II.

Ten most significant GO terms of up- and downregulated DEGs in prostate cancer.

A, Upregulated DEGs

CategoryTermCount, nP-value
CCGO:0000786; nucleosome50.0006
BPGO:0006333; chromatin assembly or disassembly60.0016
CCGO:0032993; protein-DNA complex50.0020
BPGO:0006334; nucleosome assembly50.0024
BPGO:0031497; chromatin assembly50.0027
BPGO:0065004; protein-DNA complex assembly50.0032
BPGO:0034728; nucleosome organization50.0034
BPGO:0051223; regulation of protein transport50.0070
BPGO:0006323; DNA packaging50.0077
CCGO:0000785; chromatin60.0083

B, Downregulated DEGs

CategoryTermCount, nP-value

CCGO:0005902; microvillus40.0033
CCGO:0031226; intrinsic to plasma membrane200.0055
BPGO:0009100; glycoprotein metabolic process70.0082
BPGO:0009101; glycoprotein biosynthetic process60.0120
CCGO:0044459; plasma membrane part290.0129
CCGO:0046658; anchored to plasma membrane30.0144
BPGO:0006955; immune response130.0162
CCGO:0005887; integral to plasma membrane180.0196
CCGO:0005886; plasma membrane430.0196
BPGO:0043413; biopolymer glycosylation50.0250

[i] DEGS, differentially-expressed genes; GO, Gene Ontology; CC, cellular component; BP, biological process.

Of the downregulated DEGs, BP GO terms were significantly enriched in glycoprotein metabolic process (P=0.0082), glycoprotein biosynthetic process (P=0.0120) and immune response (P=0.0162); whereas, CC GO terms were significantly enriched in the microvillus (P=0.0033), intrinsic to plasma membrane (P=0.0055) and plasma membrane part (P=0.0129) (Table II).

Identification of the upstream regulatory elements of DEGs

Motif scanning of the upstream regulatory elements of up-regulated DEGs revealed 25 transcription factors that may regulate these DEGs (Table III). Additionally, 17 transcription factors were identified to regulate downregulated DEGs (Table III).

Table III.

Potential upstream regulatory elements of DEGs in prostate cancer.

Table III.

Potential upstream regulatory elements of DEGs in prostate cancer.

ExpressionCount, nCandidate TFs
Upregulated25BARHL1, CD200, CEBPA, CUX1, EN1, ESX1, HNF4A, HOXA5, IRF8, LHX6, MEOX1, NOX1, NR1H4, PAX6, PBX1, PHOX2A, PITX2, PRKRA, RAX, RORA, SIX4, TCF7L2, VAX2, VSX2, ZEB1
Downregulated17CDX1, EPAS1, FOXC1, HMX1, HMX3, HOXA10, HOXA3, HOXC12, HSF1, HSF2, IKZF2, IRF1, IRF2, MSX1, NR3C1, POU2F1, ZEB1

[i] DEGs, differentially-expressed genes; TF, transcription factor.

Differential expression analysis of the potential upstream regulatory elements demonstrated that the expression of H6 family homeobox 3 and caudal type homeobox 1 were markedly decreased (Fig. 1A). In addition, the following genes were activated in different degrees: Retina and anterior neural fold homeobox, PBX homeobox 1, zinc finger E-box binding homeobox 1, paired box 6 and CCAAT/enhancer binding protein-α (Fig. 1B).

Detection of SNVs in prostate cancer

Following the removal of polymerase chain reaction duplicates and data processing of RNA-seq in 4 prostate cancer samples, 317, 1,004, 310 and 373 somatic mutation sites were identified in 4 prostate cancer samples particularly (Fig. 2). In each sample, base transversion (C:G → A:T; A:T → T:A; C:G → G:C) and base transition (A:T → G:C; C:G → T:A) were the most common types of base mutation (Fig. 2A). Only one insertion site was identified in the 4 prostate cancer samples. The proportions of each mutation type in the 4 prostate cancer samples were similar, suggesting that 65% of mutations were base transition and 35% were base transversion (Fig. 2B).

Mutation sites that occurred in >2 prostate cancer samples were defined as a high-frequency site (frequency ≥2). A total of 17 mutation sites were identified as high-frequency sites (Table IV). The T → A base transversion of chromosome 1189431695 was identified in 3 cancer samples and its loci was located in the folate hydrolase 1B (FOLH1B) gene. An additional 16 mutation sites were detected in 2 cancer samples and loci were all located in the coding sequence (CDS) region of genes.

Table IV.

High-frequency mutation sites in prostate cancer.

Table IV.

High-frequency mutation sites in prostate cancer.

Base

ChrPositionReferenceMutatedFrequency, n
Chr1189431695TA3
ChrY22941450AC2
ChrY22921933AG2
ChrX51807602GA2
ChrX43634502TC2
Chr996214432GC2
Chr773245681AG2
Chr629692069CA2
Chr610749928AG2
Chr1951325077GT2
Chr183253963CG2
Chr1821057145CG2
Chr1760631088AG2
Chr162149965GA2
Chr1528421681AG2
Chr189474720TC2
Chr1147580839AG2

[i] Chr, chromosome.

The frequency of high-risk mutant genes was calculated using VarioWatch (Table V). A total of 20 high-risk mutant genes were identified, of which ribosomal protein S4 Y-linked 2 (RPS4Y2), polycystin 1 transient receptor potential channel interacting (PKD1) and FOLH1B were identified in 3 cancer samples. Furthermore, kallikrein 1 (KLK1) and PKD1 were known tumor suppressor genes (30,31). Although the expression level of the KLK1 gene was increased in prostate cancer, the high frequency of mutation located in the CDS region may lead to dysfunction of the tumor suppressor. However, the expression levels of other 19 high-risk genes were not markedly different in prostate cancer, indicating that the abnormal function of these genes, induced by base mutation, may be associated with the development of prostate cancer.

Table V.

Frequency of high-risk genes in prostate cancer.

Table V.

Frequency of high-risk genes in prostate cancer.

High-risk mutated geneGene nameFrequencyExpression levelTumor-associated gene
PPP4R2Protein phosphatase 4, regulatory subunit 22
CLDN4Claudin 42 Upregulateda Tumor-associateda
FOLH1BFolate hydrolase 1B3 Downregulateda Tumor-associateda
KIF27Kinesin family member 272
MAOBMonoamine oxidase B2 Upregulateda Tumor-associateda
DUXADouble homeobox A2
KLK1Kallikrein 12UpregulatedTumor suppressor gene
MYL12AMyosin, light chain 12A, regulatory, non-sarcomeric2
HERC2Hect domain and RLD 22
HLA-EMajor histocompatibility complex, class I, E2 DownregulatedaTumor suppressor genea
HLA-GMajor histocompatibility complex, class I, G2 Downregulateda Tumor-associateda
RIOK3RIO kinase 32
RPS4Y2Ribosomal protein S4, Y-linked 23
PKD1Polycystic kidney disease 1 (autosomal dominant)3 UpregulatedaTumor suppressor gene
USP10Ubiquitin specific peptidase 102
TLK2Tousled-like kinase 22
TMEM14BTransmembrane protein 14D; transmembrane protein 14B2
TUBA1ATubulin, α 1a2
GBP3Guanylate binding protein 32
MAGED4BMelanoma antigen family D, 4B; melanoma antigen family D, 42Tumor-associated gene

a Data provided by the OncoSearch database (http://oncosearch.biopathway.org/#/search).

Discussion

The present study provided a survey of DEGs and mutant genes in human prostate cancer. RNA-seq data between 4 prostate cancer samples and 4 matched normal samples were analyzed. A total 150 upregulated and 211 downregulated DEGs were identified, 4 of which were upregulated transcription factors (DMBX1, NCOA2, ONECUT2 and ZNF83), and 4 were downregulated transcription factors (IFI16, NEUROG3, RARG and SIM1).

NCOA2 has been suggested as an oncogene in primary tumors by increasing androgen receptor signaling, which is known to have a function in early and late-stage prostate cancer (13). Furthermore, the upregulated expression of transcription factor ONECUT2 was revealed in breast and prostate cancer cell lines (32), and increased IFI16 protein in normal human prostate epithelial cells was associated with cellular senescence-associated cell growth arrest (33). NEUROG3 was identified to be expressed in metastatic neuroendocrine prostate cancer cells (34). However, mechanism of abnormal regulation of DMBX1, ZNF83, RARG and SIM1 in prostate cancer remains unknown. For example, the abnormal regulation of ZNF83 has been identified in hepatocellular carcinoma (35) and colorectal cancer (36), but not in prostate cancer. DMBX1 is a paired-class homeodomain transcription factor and may be determined in the brain, stomach and testis in adult normal tissues (37). A previous study focused on the expression pattern of DMBX1 in the development of the neural network (38). Therefore, these aforementioned transcription factors may have important functions in the progress of prostate cancer.

Determining the SNVs in prostate cancer identified 17 mutation sites with a frequency ≥2. The mutation site with the highest frequency was located in FOLH1B. Furthermore, 20 high-risk mutant genes with highest frequency were identified using VarioWatch and included RPS4Y2, PKD1 and FOLH1B. FOLH1B originates from the duplication of folate hydrolase 1 (FOLH1) (39), which is an established biomarker for prostate cancer (40). FOLH1, also known as prostate-specific membrane antigen 1, is used as a diagnostic and prognostic indicator for prostate cancer, and is associated with aggressiveness and metastasis of prostate cancer (41). Of the identified high-risk mutant genes exhibiting the highest frequency, KLK1 and PKD1 are known tumor suppressor genes and MAGE family member D4B is a tumor-associated gene (30,31). PKD1 has been identified to be downregulated in advanced-stage prostate cancer and was present as a protein complex, combined with the androgen receptor, in prostate cancer cells (42). The kallikrein-related peptidases have been identified in a number of types of cancer, including prostate and ovarian, and a combination of the dysregulation of KLK1, 5 and 13 was associated with poorer disease-free survival for prostate cancer (43). The high frequency of mutation located in KLK1, PKD1 and D4B may lead to dysfunction of the tumor suppressor functions and consequently contribute to the progress of prostate cancer.

There were a number of limitations in the present study; for example, the sample size used for the analysis was small. Furthermore, the results of the present study require validation using other RNA-seq data or microarray data of prostate cancer.

Combined with bioinformatics methods, the RNA-seq data were analyzed to determine candidate genes for diagnosing and/or treating of prostate cancer. The identified DEGs (DMBX1, ZNF83, RARG and SIM1) and mutant genes (RPS4Y2, KLK1 and FOLH1B) may have important functions in the progression of prostate cancer. The results of the present study may enable an improved understanding of the molecular mechanisms that underlie prostate cancer pathogenesis.

References

1 

Siegel R, Ward E, Brawley O and Jemal A: Cancer statistics, 2011: The impact of eliminating socioeconomic and racial disparities on premature cancer deaths. CA Cancer J Clin. 61:212–236. 2011. View Article : Google Scholar : PubMed/NCBI

2 

Jemal A, Bray F, Center MM, Ferlay J, Ward E and Forman D: Global cancer statistics. CA Cancer J Clin. 61:69–90. 2011. View Article : Google Scholar : PubMed/NCBI

3 

Hsing AW and Chokkalingam AP: Prostate cancer epidemiology. Front Biosci. 11:1388–1413. 2006. View Article : Google Scholar : PubMed/NCBI

4 

Miller DC, Hafez KS, Stewart A, Montie JE and Wei JT: Prostate carcinoma presentation, diagnosis and staging: An update form the National Cancer Data Base. Cancer. 98:1169–1178. 2003. View Article : Google Scholar : PubMed/NCBI

5 

Leitzmann MF and Rohrmann S: Risk factors for the onset of prostatic cancer: Age, location, and behavioral correlates. Clin Epidemiol. 4:1–11. 2012. View Article : Google Scholar : PubMed/NCBI

6 

Schaid DJ: The complex genetic epidemiology of prostate cancer. Hum Mol Genet 13 Spec No. 1:R103–R121. 2004. View Article : Google Scholar

7 

Mahmoud AM, Yang W and Bosland MC: Soy isoflavones and prostate cancer: A review of molecular mechanism. J Steroid Biochem Mol Biol. 140:116–132. 2014. View Article : Google Scholar : PubMed/NCBI

8 

Watson PA, Arora VK and Sawyers CL: Emerging mechanisms of resistance to androgen receptor inhibitors in prostate cancer. Nat Rev Cancer. 15:701–711. 2015. View Article : Google Scholar : PubMed/NCBI

9 

Pflueger D, Terry S, Sboner A, Habegger L, Esgueva R, Lin PC, Svensson MA, Kitabayashi N, Moss BJ, MacDonald TY, et al: Discovery of non-ETS gene fusions in human prostate cancer using next-generation RNA sequencing. Genome Res. 21:56–67. 2011. View Article : Google Scholar : PubMed/NCBI

10 

Narla G, Heath KE, Reeves HL, Li D, Giono LE, Kimmelman AC, Glucksman MJ, Narla J, Eng FJ, Chan AM, et al: KLF6, a candidate tumor suppressor gene mutated in prostate cancer. Science. 294:2563–2566. 2001. View Article : Google Scholar : PubMed/NCBI

11 

Majid S, Dar AA, Shahryari V, Hirata H, Ahmad A, Saini S, Tanaka Y, Dahiya AV and Dahiya R: Genistein reverses hypermethylation and induces active histone modifications in tumor suppressor gene B-Cell translocation gene 3 in prostate cancer. Cancer. 116:66–76. 2010.PubMed/NCBI

12 

Carver BS, Tran J, Gopalan A, Chen Z, Shaikh S, Carracedo A, Alimonti A, Nardella C, Varmeh S, Scardino PT, et al: Aberrant ERG expression cooperates with loss of PTEN to promote cancer progression in the prostate. Nat Genet. 41:619–624. 2009. View Article : Google Scholar : PubMed/NCBI

13 

Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, Arora VK, Kaushik P, Cerami E, Reva B, et al: Integrative genomic profiling of human prostate cancer. Cancer cell. 18:11–22. 2010. View Article : Google Scholar : PubMed/NCBI

14 

Wang Z, Gerstein M and Snyder M: RNA-Seq: A revolutionary tool for transcriptomics. Nat Rev Genet. 10:57–63. 2009. View Article : Google Scholar : PubMed/NCBI

15 

Ren S, Peng Z, Mao JH, Yu Y, Yin C, Gao X, Cui Z, Zhang J, Yi K, Xu W, et al: RNA-seq analysis of prostate cancer in the Chinese population identifies recurrent gene fusions, cancer-associated long noncoding RNAs and aberrant alternative splicings. Cell Res. 22:806–821. 2012. View Article : Google Scholar : PubMed/NCBI

16 

Kannan K, Wang L, Wang J, Ittmann MM, Li W and Yen L: Recurrent chimeric RNAs enriched in human prostate cancer identified by deep sequencing. Proc Natl Acad Sci USA. 108:9172–9177. 2011. View Article : Google Scholar : PubMed/NCBI

17 

Xu X, Zhu K, Liu F, Wang Y, Shen J, Jin J, Wang Z, Chen L, Li J and Xu M: Identification of somatic mutations in human prostate cancer by RNA-Seq. Gene. 519:343–347. 2013. View Article : Google Scholar : PubMed/NCBI

18 

Makarov DV, Trock BJ, Humphreys EB, Mangold LA, Walsh PC, Epstein JI and Partin AW: Updated nomogram to predict pathological stage of prostate cance given prostate-specific antigen level, clinical stage, and biopsy Gleason score (Partin tables) based on cases from 2000 to 2005. Urology. 69:1095–1101. 2007. View Article : Google Scholar : PubMed/NCBI

19 

Trapnell C, Pachter L and Salzberg SL: TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics. 25:1105–1111. 2009. View Article : Google Scholar : PubMed/NCBI

20 

Pruitt KD, Tatusova T and Maglott DR: NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 35 (Database Issue). D61–D65. 2007. View Article : Google Scholar

21 

Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL and Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 7:562–578. 2012. View Article : Google Scholar : PubMed/NCBI

22 

Beane J, Vick J, Schembri F, Anderlind C, Gower A, Campbell J, Luo L, Zhang XH, Xiao J, Alekseyev YO, et al: Characterizing the impact of smoking and lung cancer on the airway transcriptome using RNA-Seq. Cancer Prev Res (Phila). 4:803–817. 2011. View Article : Google Scholar : PubMed/NCBI

23 

Hulsegge I, Kommadath A and Smits MA: Globaltest and GOEAST: Two different approaches for Gene Ontology analysis. BMC Proc. 3 Suppl 4:S102009. View Article : Google Scholar : PubMed/NCBI

24 

Huang da W, Sherman BT and Lempicki R: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 4:44–57. 2009. View Article : Google Scholar

25 

Zhao M, Sun J and Zhao Z: TSGene: A web resource for tumor suppressor genes. Nucleic Acids Res. 41:(Database Issue). D970–D976. 2013. View Article : Google Scholar : PubMed/NCBI

26 

Chen JS, Hung WS, Chan HH, Tsai SJ and Sun HS: In silico identification of oncogenic potential of fyn-related kinase in hepatocellular carcinoma. Bioinformatics. 29:420–427. 2013. View Article : Google Scholar : PubMed/NCBI

27 

He HH, Meyer CA, Shin H, Bailey ST, Wei G, Wang Q, Zhang Y, Xu K, Ni M, Lupien M, et al: Nucleosome dynamics define transcriptional enhancers. Nat Genet. 42:343–347. 2010. View Article : Google Scholar : PubMed/NCBI

28 

McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M and DePristo MA: The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20:1297–1303. 2010. View Article : Google Scholar : PubMed/NCBI

29 

Cheng YC, Hsiao FC, Yeh EC, Lin WJ, Tang CY, Tseng HC, Wu HT, Liu CK, Chen CC, Chen YT and Yao A: VarioWatch: Providing large-scale and comprehensive annotations on human genomic variants in the next generation sequencing era. Nucleic Acids Res. 40:(Web Server Issue). W76–W81. 2012. View Article : Google Scholar : PubMed/NCBI

30 

Raea F, Bulmera B, Nicolb D and Clementsa J: The human tissue kallikreins (KLKs 1–3) and a novel KLK1 mRNA transcript are expressed in a renal cell carcinoma cDNA library. Immnuopharmacology. 45:83–88. 1999. View Article : Google Scholar

31 

Shabelnik MY, Kovalevska LM, Yurchenko MY, Shlapatska LM, Rzepetsky Y and Sidorenko S: Differential expression of PKD1 and PKD2 in gastric cancer and analysis of PKD1 and PKD2 function in the model system. Exp Oncol. 33:206–211. 2011.PubMed/NCBI

32 

Yamamoto F and Yamamoto M: Scanning copy number and gene expression on the 18q21-qter chromosomal region by the systematic multiplex PCR and reverse transcription-PCR methods. Electrophoresis. 28:1882–1895. 2007. View Article : Google Scholar : PubMed/NCBI

33 

Alimirah F, Chen J, Davis FJ and Choubey D: IFI16 in human prostate cancer. Mol Cancer Res. 5:251–259. 2007. View Article : Google Scholar : PubMed/NCBI

34 

Telesca D, Inoue LY, Neira M, Etzioni R, Gleave M and Nelson C: Differential expression and network inferences through functional data modeling. Biometrics. 65:793–804. 2009. View Article : Google Scholar : PubMed/NCBI

35 

Dong H, Ge X, Shen Y, Chen L, Kong Y, Zhang H, Man X, Tang L, Yuan H, Wang H, et al: Gene expression profile analysis of human hepatocellular carcinoma using SAGE and LongSAGE. BMC Med Genomics. 2:52009. View Article : Google Scholar : PubMed/NCBI

36 

Jovov B, Araujo-Perez F, Sigel CS, Stratford JK, McCoy AN, Yeh JJ and Keku T: Differential gene expression between African American and European American colorectal cancer patients. PLoS One. 7:e301682012. View Article : Google Scholar : PubMed/NCBI

37 

Ohtoshi A, Nishijima I, Justice MJ and Behringer RR: Dmbx1, A novel evolutionarily conserved paired-like homeobox gene expressed in the brain of mouse embryos. Mech Dev. 110:241–244. 2002. View Article : Google Scholar : PubMed/NCBI

38 

Fujimoto W, Shiuchi T, Miki T, Minokoshi Y, Takahashi Y, Takeuchi A, Kimura K, Saito M, Iwanaga T and Seino S: Dmbx1 is essential in agouti-related protein action. Proc Natl Acad Sci USA. 104:15514–15519. 2007. View Article : Google Scholar : PubMed/NCBI

39 

De Lorenzi L, Genualdo V, Gimelli S, Rossi E, Perucatti A, Iannuzzi A, Zannotti M, Malagutti L, Molteni L, Iannuzzi L and Parma P: Genomic analysis of cattle rob (1;29). Chromosome Res. 20:815–823. 2012. View Article : Google Scholar

40 

Zhang T, Song B, Zhu W, Xu X, Gong QQ, Morando C, Dassopoulos T, Newberry RD, Hunt SR and Li E: An Ileal Crohn's disease gene signature based on whole human genome expression profiles of disease unaffected ileal mucosal biopsies. PLoS One. 7:e371392012. View Article : Google Scholar : PubMed/NCBI

41 

Ergün A, Lawrence CA, Kohanski MA, Brennan TA and Collins JJ: A network biology approach to prostate cancer. Mol Syst Biol. 3:822007. View Article : Google Scholar : PubMed/NCBI

42 

Mak P, Jaggi M, Syed V, Chauhan SC, Hassan S, Biswas H and Balaji KC: Protein kinase D1 (PKD1) influences androgen receptor (AR) function in prostate cancer cells. Biochem Biophys Res Commun. 373:618–623. 2008. View Article : Google Scholar : PubMed/NCBI

43 

Girgis AH, Bui A, White NM and Yousef GM: Integrated genomic characterization of the kallikrein gene locus in cancer. Anticancer Res. 32:957–963. 2012.PubMed/NCBI

Related Articles

Journal Cover

April-2018
Volume 15 Issue 4

Print ISSN: 1792-1074
Online ISSN:1792-1082

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Zhang Z, Wu H, Zhou H, Gu Y, Bai Y, Yu S, An R and Qi J: Identification of potential key genes and high-frequency mutant genes in prostate cancer by using RNA-Seq data. Oncol Lett 15: 4550-4556, 2018
APA
Zhang, Z., Wu, H., Zhou, H., Gu, Y., Bai, Y., Yu, S. ... Qi, J. (2018). Identification of potential key genes and high-frequency mutant genes in prostate cancer by using RNA-Seq data. Oncology Letters, 15, 4550-4556. https://doi.org/10.3892/ol.2018.7846
MLA
Zhang, Z., Wu, H., Zhou, H., Gu, Y., Bai, Y., Yu, S., An, R., Qi, J."Identification of potential key genes and high-frequency mutant genes in prostate cancer by using RNA-Seq data". Oncology Letters 15.4 (2018): 4550-4556.
Chicago
Zhang, Z., Wu, H., Zhou, H., Gu, Y., Bai, Y., Yu, S., An, R., Qi, J."Identification of potential key genes and high-frequency mutant genes in prostate cancer by using RNA-Seq data". Oncology Letters 15, no. 4 (2018): 4550-4556. https://doi.org/10.3892/ol.2018.7846