Integrated analysis of miRNA landscape and cellular networking pathways in stage-specific prostate cancer

Dysregulation of miRNAs has been demonstrated in several human malignancies including prostate cancer. Due to tissue limitation and variable disease progression, stage-specific miRNAs changes in prostate cancer is unknown. Using chip-based microarray, we investigated global miRNA expression in human prostate cancer LNCaP, PC3, DU145 and 22Rv1 cells representing early-stage, advanced-stage and castration resistant prostate cancer in comparison with normal prostate epithelial cells. A total of 292 miRNAs were differentially expressed with 125 upregulated and 167 downregulated. These miRNAs were involved in pathways including drug resistance drug-efflux, adipogenesis, epithelial-to-mesenchymal transition, bone metamorphosis, and Th1/Th2 signaling. Regulation of miRNAs were interlinked with upstream regulators such as Argonaut 2 (AGO2), Double-Stranded RNA-Specific Endoribonuclease (DICER1), Sjogren syndrome antigen B (SSB), neurofibromatosis 2 (NF2), and peroxisome proliferator activated receptor alpha (PPARA), activated during stage-specific disease progression. Candidate target genes and pathways dysregulated in stage-specific prostate cancer were identified using CS-miRTar database and confirmed in clinical specimens. Integrative network analysis suggested some genes targeted by miRNAs include miR-17, let7g, miR-146, miR-204, miR-205, miR-221, miR-301 and miR-520 having a major effect on their dysregulation in prostate cancer. MiRNA-microarray analysis further identified miR-130a, miR-181, miR-328, miR146 and miR-200 as a panel of novel miRNAs associated with drug resistance drug-efflux and epithelial-to-mesenchymal transition in prostate cancer. Our findings provide evidence on miRNA dysregulation and its association with key functional components in stage-specific prostate cancer.


Introduction
Prostate cancer, by far, remains the leading cause of cancer morbidity and mortality in males in the United States, and approximately 20% of these men will develop invasive disease during their lifetime [1]. According to an assessment by the American Cancer Society, in the year 2019, a total of 174,650 men will be diagnosed with prostate cancer and 31,620 prostate cancer-related deaths are predicted [2]. Ethnicity, race, age and family history are some non-modifiable risk factors for prostate cancer, whereas environment and lifestyle factors are modifiable risk factors that contribute to its development [3]. Prostate cancer in humans exhibits a distinctive continuum of features including heterogeneity, multi-focality, inconsistent clinical progression including metastasis to bone, and emergence of androgen-independent disease [4]. The growth of early prostate lesions are androgen dependent, progressing to an androgen-refractory stage, ultimately resulting in patients' death [5]. Current clinical examination including non-invasive imaging, pathological findings, and serum markers are inadequate in prognostication of prostate cancer. Most patients diagnosed with advance-stage disease are treated with androgen deprivation therapy [6], which results in tumor shrinkage [7]. About 70-80% of patients initially respond to this therapy, but the tumor eventually becomes resistant with the emergence of castration-resistant prostate cancer (CRPC) [8]. Although primary prostate cancer is manageable with surgery and/or radiation, however no effective treatment is still available for metastatic CRPC [9]. Recent therapeutic approaches for metastatic forms of prostate cancer are palliative rather than curative. Therefore, understanding the molecular abnormalities during various stages of prostate cancer, in particular CRPC, will lead to the development of approaches for early detection, prevention and therapeutic intervention benefiting patients.
MicroRNAs (miRNAs) are a class of small non-coding regulatory RNAs which are endogenously synthesized molecules about 20-22 nucleotides in length. These miRNAs are well conserved among mammalian species and have essential role in the regulation of �60% of human genes [10]. MiRNAs interfere with the targets' expression at the posttranscriptional stage affecting proliferation and differentiation including other cellular pathways. Studies demonstrate that miRNA expression profiles could distinguish between normal and tumor tissues to identify tumors and their subtypes [11][12][13]. MiRNA can be used as biomarker predicting disease recurrence and clinical outcome in some tumor types [14,15]. Therefore, miRNAs have emerged as promising tools for diagnosis, prognosis and tumor response in patients. It is expected that an integrated analysis of miRNAs landscape and cellular networking pathways in stage-specific prostate cancer will help in understanding the molecular basis of pathogenesis and is likely to provide information on early diagnosis/prognosis and targets for therapeutic intervention.
In our previous study, we determined the miRNA complement of normal prostate epithelial and stromal cells to obtain a comprehensive assessment of the differential expression of these molecules as potential regulators of various signaling pathways [16]. In the present study, using a microarray platform, we compared the miRNA profiles of various established human prostate cancer cell lines viz. LNCaP, PC3, DU145, 22Rv1 with that of normal prostate epithelial cells (PrEc). These cell lines are representative of various stages of prostate cancer including (i) early-stage with androgen sensitivity (LNCaP) (ii) advanced-stage with the loss of androgen receptor (AR) function (PC3 and DU145) and (iii) castration resistant prostate cancer harboring repressed AR (22Rv1). We designed the study to understand the pathological characteristics of prostate cancer linked with miRNA expression at a particular disease stage with emphasis on canonical pathways, upstream regulators, target molecules (miRNA: mRNA) and their regulatory genes followed by clinical validation.

Cell culture
Human prostate cancer cell lines LNCaP, PC3, DU145, 22Rv1 and the normal epithelial prostate cells (PrEc) were used in the study. LNCaP cells were established from a supraclavicular lymph node lesion of human prostate adenocarcinoma. These cells possess wild-type AR, are androgen responsive, and secret prostate-specific antigen (PSA) [17]. The PC3 cells were established from the bone metastasis of a grade IV prostatic adenocarcinoma, are androgenindependent and do not express AR or PSA [18]. Similarly, DU145 cells were derived from a brain metastatic site with mutation in AR and are PSA negative [19]. The 22Rv1 cells were derived from a human prostatic carcinoma xenograft after serial propagation in mice with castration-induced regression and relapse of the parental, androgen-dependent CWR22 cells [20]. All cancer cells were maintained in RPMI1640 cell culture medium (Gibco, Gaithersburg, MD) supplemented with 10% fetal bovine serum (Hyclone, GE Healthcare, Pittsburgh, PA), 100 unit's/ml penicillin (Gibco, Grand Island, NY) and 100 μg/ml streptomycin (Gibco, Grand Island, NY). The cells were maintained in an incubator with a humidified atmosphere of 95% air and 5% CO2 at 37˚C. The PrEc cells were cultured in keratinocyte growth medium with supplements including 5 ng ml-1 human recombinant epidermal growth factor and 0.05 mg ml-1 bovine pituitary extract (Life Technologies/Invitrogen, Carlsbad, CA).

RNA isolation
Total RNA was extracted from LNCaP, PC3, DU145, 22Rv1, and PrEc prostate cell lines using RNeasy Mini kit (Qiagen, Germantown, MD) following manufacturer's protocols. DNA was removed by performing on-column DNase digestion with RNase-free DNase (Qiagen, Germantown, MD). The integrity and size distribution of RNA was examined by agarose gel electrophoresis. From each cell line RNA was isolated in three biological replicates and 3 technical replicates. For quality control, the isolated RNAs were analyzed using the Agilent Bio-analyzer 2100, and RNA with RIN numbers �7.2 were used in experiments.

MiRNA-microarray
MiRNA expression levels were profiled using the Affymetrix array from GenoSensor TM (Tempe, AZ, catlog#1201C) microRNA array version 2007. RNA was isolated from LNCaP, DU145, PC-3, 22Rv1, and PrEc cells and labeled according to the Agilent protocol (version 1.0, April 2007). Arrays were scanned with Agilent Microarray G2565, and Agilent Feature Extraction Software (version 9.5) was used to extract the data. The data summary provides gene names, intensity on three redundancies including average intensity of the triplicates, standard deviation, normalization, and fold changes. Signal intensity measured through fluorescent signal of each chemical detected with either 635 nm or 532 nm and the signals were labeled with either F635 or F532. Bacterial sequences/probes such as gnd and fixB were used as negative and QC controls, which are non-homologous to any targets in the reaction.

Data analysis and normalization
Normalization represents the normalization ratio of gene signal intensity of a specific gene versus global gene signal intensity. The normalization value was then used for the fold change calculation. For quantitative normalization and subsequent data processing GenoExplorer software was used. Quantile normalization of the raw data was performed and differentially expressed miRNAs were identified through volcano plot (p value <0.001) filtering. Hierarchical clustering was accomplished using GraphPad Prism 7 (https://www.graphpad.com) to distinguish differential miRNA expression profile among the samples. The miRNA data of all the cell lines were submitted to the GEO public domain portal with accession number GSE119338.

Functional and pathway enrichment analyses
The log2 ratios and P values for each identified RNA in the microarray dataset was inserted into Ingenuity Pathway Analysis (IPA) using the core analysis platform (Ingenuity Systems, Redwood City, CA). The differentially expressed microarray data of prostate cell lines were matched with those in Ingenuity Knowledge Base to generate molecular networks. A list of biological functions overrepresented in the dataset, and determine overrepresented canonical pathways; unmapped RNAs were excluded from further analysis. The p value cutoff of 0.001 was applied. To evaluate the actual over-represented pathways, or to remove the chances of any randomness in data, the more stringent criteria of the Benjamini-Hochberg (B-H) procedure was used.

Analysis of upstream regulators
The association of miRNA with the upstream regulators, especially transcription factors (TFs), was analyzed using the Ingenuity Pathway Analysis (IPA; Qiagen). The complete dataset with information on Gene ID, FDR, and expression ratio of the original transcriptomic dataset and the datasets corrected by miRNA were used. The annotated microarray was used as background in all analyses. The comparison analysis feature in the IPA Ingenuity1 Knowledge Base was used for visualization and analysis which allowed identification of the upstream regulators of differentially expressed genes, prediction of activation Z-score and an overlap P-value for each molecule. The activation Z-score estimates the status of the upstream regulator by using gene expression level of known target gene. A prediction of activation or inhibition (or no prediction) records for the chance that random data generation for significant predictions. Z-scores greater than 2 or smaller than −2 can be considered as significant.

miRNA and mRNA target interaction
The target of miRNAs was identified using mir-Tar (https://bio.tools/mirtarbase) an integrative web server for identifying miRNA target interactions in humans with Gibbs free energy ΔG (kcal/mol) and alignment score >150.

Prediction of miRNA target genes
We used CSmiR-Tar database to identify miRNA targets expressed in tissue-and stage-specific disease progression [21]. It also allowed to identify specific target(s) of miRNA using two filter base. The first filter offers target expression in specific tissue, whereas the second filter is disease specific target. Both criteria were used to identify miRNAs and their specific targets in prostate tissue in three stage-specific manner including early stage, advanced stage and CRPC. In particular, these 3 stages comprise of 6 different sub-stages; high-grade prostate intraepithelial neoplasia (putative prostate cancer lesion), prostate adenocarcinoma, androgen-independent prostate cancer and metastatic disease as advance stage prostate cancer, whereas patients exhibiting recurrence after androgen deprivation therapy as CRPC-hormone refractory prostate cancer. Formula used to calculate average normalized score include: Normalized score = Original score-Minimum score / Maximum score-Minimum score; Average normalized score (ANS) = sum of the normalized score from different databases / Number of databases (n = 3).

Statistical analysis
The overlap p-value identified transcriptional regulators having ability to elucidate observed gene expression changes. The overlap p-value measures statistically significant overlap between the dataset genes and the genes that are regulated by an upstream regulator, calculated using Fisher's exact test with p-values <0.01 considered to be significant.

Identification of differentially expressed miRNA
The differentially expressed miRNAs in all 4 prostate cancer cell lines including LNCaP, PC3, DU145 and 22Rv1 compared with PrEc are represented in the form of heat map S1 Fig. The red color in the panel shows an increase, whereas the blue color exhibited decrease in miRNA expression with p value <0.001. Results obtained from human prostate cancer LNCaP cells, representative of early-stage cancer, revealed a total of 52 differentially expressed miRNAs, of which, 8 miRNAs were upregulated and 44 miRNAs were downregulated S1 Table. In PC3 cells, 116 miRNAs were differentially expressed including 107 upregulated and 9 downregulated S2 Table. In DU145 cells, 66 miRNAs were differentially expressed with upregulation of 1 miRNA (miR-301; 2.25 fold change) and 65 were downregulated S3 Table. These two cell lines represent advanced-stage cancer. In 22Rv1 cells, a representative of castration-resistant prostate cancer (CRPC), 58 miRNAs were differentially expressed including 9 miRNAs upregulated and 49 downregulated with a fold change cut-off �1.5 S4 Table. The data was scrutinized using statistical software graph pad prism (sigma plot) and generated a volcano plot for visualizing differential expression patterns between two different conditions. The statistically significant differentially expressed miRNAs for all 4 cell lines were represented using log2 fold changes (Fig 1).

miRNA-upstream transcriptional regulator analysis
The IPA knowledge database identified several miRNAs which post-transcriptionally modulate the expression of upstream regulators during particular disease stage viz. early, advance and castration-resistant disease progression (Fig 4A-4P). For instance, the expression of miR- The volcano plot exhibits the relationship between fold-change and significance between the two groups, using a scatter plot view. The y-axis is the negative log2 of P values (a higher value indicates greater significance) and the x-axis is the difference in expression between two experimental groups as measured in log2 space. The blue dots in the figure represents individual miRNAs differentially expressed in four prostate cancer cell lines LNCaP (early-stage), PC3, DU145 (advanced-stage) and 22Rv1 (castration resistant prostate cancer; CRPC). The vertical dotted lines correspond to p value < 0.01 and <0.001, respectively. The horizontal dotted line represents a log 2 fold changes. The left side in the panel corresponds to differentially expressed miRNAs with decreased expression, and the right side in the panel showed miRNA with increased expression.
17 was upregulated during early-stage, advance-stage and CRPC (red color) and regulated the expression of AGO2, SSB, DICER1 and PPARA; whereas miR-221 and miR-205 were downregulated (green color) during all stages of disease progression. However, let-7 was differentially

miRNA and gene regulation
The CSmiR-Tar database was applied to identify miRNA targets expressed during stage-specific disease progression. A list of major miRNA regulating stage-specific prostate cancer progression include miR-17, let7g, miR-146, miR-204, miR-205, miR-221, miR-301 and miR-520. Several genes have been identified which were linked to various signaling pathways altered by the above mentioned miRNAs during prostate cancer progression (Fig 6; S2 Fig). Other candidate miRNAs potentially involved in modulating the expression of target genes in prostate cancer are listed in Table 1.

Gene ontology analysis
Gene ontology analysis of the target genes by PANTHER identified their involvement in different clusters of molecular functions, cellular component, biological process and pathways. Significant results obtained through GO analysis prompt us to uncover stage-specific prostate cancer genes among various target genes. Manual screening of target genes with NCBI Gene database and OMIM prostate cancer genes identified a subset of target genes and pathways which overlap in stage-specific disease progression (Fig 7; Table 1).

Clinical validation
Next, the parallel matched miRNA profiles of cell lines were validated using miRNA profiles of human prostate specimens. Data revealed significant correlation with clinical specimens such as the upregulated expression of miR-301, miR-520, miR-27a and miR-214 in early-stage and CRPC. Likewise, the expression of miR-214 and miR-18b was consistently upregulated across all prostate cancer cell lines and human patient specimens. In contrast, expression of miR-205 was consistently downregulated in all cell lines and human patient cohort. The expression of miR-021 was in-parallel downregulated in all prostate cancer cell lines but upregulated in human dataset (Table 2).

Discussion
Dysregulation of miRNAs has been reported in prostate malignancy from early-to advancedstage and castration resistant disease progression. It is known that miRNAs reside within the fragile regions of the genome, and are implicated to play an important role in cancer progression. Using glass chip-based miRNA-microarray we comprehensively evaluated and profiled miRNA expression using four prostate cancer cell lines as representative of early-stage, advanced-stage and CRPC progression, compared with normal prostate epithelial cells as control (Fig 1; S1 Fig). This study focuses on the comprehensive bioinformatics analysis of miRNA profiles connecting canonical pathways, upstream regulators and their regulatory genes in stage-specific prostate cancer. The associated signaling pathways identified in prostate cancer include adipogenesis, bone metamorphosis, drug resistance drug efflux, epithelial-mesenchymal transition, Th1 and Th2 pathways (Fig 3A), and the upstream regulators viz. Argonaute 2 (AGO2), Double-Stranded RNA-Specific Endoribonuclease (DICER1), Sjogren syndrome antigen B (SSB), neurofibromatosis 2 (NF2), and peroxisome proliferator activated receptor alpha (PPARA) (Fig 3B). In addition, miRNA regulatory genes during stage-specific prostate cancer progression suggest the involvement of p53, EGFR-PI3K-Akt, IGF, interleukins, TGFB, VEGF, JAK/STAT, Wnt signaling and their effectors as the most critical genes in prostate cancer via upregulation of growth factor receptors, specifically EGFR, or through PTEN inactivation (Fig 7).
Among differentially expressed miRNAs, most interestingly, the expression of miR-301 was upregulated in early-stage and CRPC progression, and this high expression of miR-301 was consistent in both serum and tumor tissue in prostate cancer patients compared to patients with benign prostate hyperplasia [22]. The CSmiR-Tar database identified erythroblastic oncogene B 2 (ERBB2) as a potential target gene of miR-301, across all cancer stages, with a binding site in the 5'UTR of ERBB2 which is structurally related to EGFR (Fig 5) [23]. However, the role of miR-301 in context to its target gene EGFR has not yet been explored. Moreover, EGFR was predominantly upregulated in DU145 and PC3 cells, compared to an androgen-responsive     25], in which both EGFR and ERBB2 were expressed [26]. EGF/EGFR signaling activates ERK and Akt pathways [27], that markedly promotes cell migration and survival. Taken together, the above data support that dysregulation of ERBB2/EGFR activates downstream pro-oncogenic pathways which include AKT-PI3K-mTOR pathways, critical for cancer cell migration and proliferation. Based on the above information, miR-301 plays a critical role in prostate cancer progression and may either be developed as a prognostic marker or therapeutic target [28]. Another miRNA, miR-146b/146 exhibited higher expression (5.38/4.7 fold) during advanced-stage prostate cancer. Consistent with the analysis performed in 42 patient samples including 15 benign prostatic hyperplasia, 15 prostate cancer, and 9 treatment-resistant prostate cancer, the expression of miR-146b was upregulated during CRPC progression (GSE59156). Similar overexpression pattern was found in human patient cohort [29], and in androgen-refractory prostate cancer DU145 and PC3 cells [29]. With reference to its target interaction, miR-146 binds at 3´UTR coding region of EGFR during early-stage, PTGS2 in advance-stage and EGFR in CRPC progression (Fig 5). Previous reports suggested that miR-146-5p represses EGFR expression through binding to its 3'-untranslated region [30], and that the reconstitution of miR-146b-5p may be useful for the treatment of invasive prostate cancer [31].

hsa-miR-384
The expression of miR-205 was downregulated in all prostate cancer cell lines progressively in stage-specific disease progression (S1-S4 Tables). Similar pattern of downregulation of miR-205 was noted in prostate cancer tissues from Caucasian-American and African-American men, compared to adjacent normal tissue [32]. It is speculated that downregulation of miR-205 is inversely correlated with methylation status of miR-205 promoter and locus (S2 Fig). miR-205 actively targets the 3'UTR of some key genes such as AR, BCL2, ERBB3, PTEN, and VEGF-A that play important roles in cancer cell invasion and metastasis (Fig 6). During  Table 2. List of potential miRNAs and their expression (Up or Down) were matched parallel with human prostate cancer sample cohorts. The expression of miRNAs from human patient's samples were derived from three GEO databases including the recurrent and non-recurrent cases (GSE88958), with primary prostate cancer stem cells (GSE59156), and the miRNA expression profile of prostate cancer compared with normal prostate tissue (GSE60117). advanced-stage and CRPC, miR-205 targets the same set of genes (PTEN and VEGF-A) along with AR, and E2F1 transcription factor. Moreover, downregulation of miR-205 was inversely correlated with BCL-2 upregulation [33], an anti-apoptotic oncoprotein required for prostate cancer progression [34]. Mechanistically, miR-205 targets the ERBB3 receptor and inhibits Akt activation; the expression of ERBB3 was correlated with biochemical recurrence. Conclusively, miR-205 is a novel molecule having potential to serve as a predictive biomarker and therapeutic target in prostate cancer [35]. Another miRNA, miR-221 was downregulated in the miRNA-microarray dataset (S1-S4 Tables). Downregulation of miR-221 was documented in LNCaP cells as early-stage cancer [36,37], and underlines the involvement of miR-221 in the expansion and maintenance of CRPC phenotype [38]. The miR-221-3p targets genes at early stage such as Annexin A1 (ANXA1), Cyclin Dependent Kinase Inhibitor 1B (CDKN1B), and multiple genes during advance-stage cancer such as CDKN1B, PTEN, Ring Finger Protein 20 (RNF20), Estrogen Receptor 1 (ESR1) and AKT3 (Fig 6; Table 1). In CRPC, miR-221 targets Adenomatous polyposis coli (APC), and AR (Fig 6). Moreover, in context to unique/signature disease stage specific miRNAs, upregulation of miR-204 was associated with early stage whereas downregulation of miR-520a was identified as signature for CRPC (Fig 6). The miR-204 acts as tumor suppressor in human cancer by targeting BCL2 at an early-stage, and, thus protecting cancer cells from undergoing apoptosis as a resistance to androgen deprivation therapy [39]. During advancedstage, miR-204 targets Runt Related Transcription Factor 2 (RUNX2), Baculoviral IAP Repeat Containing 2 (BIRC2). RUNX2 is a bone-specific transcriptional regulator, aberrantly expressed in metastatic prostate cancer cells [40] and BIRC2, also known as Inhibitor of Apoptosis (IAP) plays key role in drug resistance and survival of cancer cells through inhibition of apoptosis [41]. During CRPC progression, miR-204 targets SIRT1 gene, which is involved in cancer drug resistance [42] by FOXO1 deacetylation [43,44].
Another miRNA, miR-520a binds to the 3' UTR of target genes including Microtubule Affinity Regulating Kinase 2 (MARK2), Estrogen Receptor 1 (ESR1) at early-stage cancer and Suppressor of Zeste 12 (SUZ12) during advance-stage and ESR1 during CRPC progression of prostate cancer ( Table 1). The role of MARK2 with respect to prostate cancer was not yet been studied, while previous reports showed that ESR1 play protective roles in prostate cancer development [45], though still unknown. Moreover, the role of SUZ12 was pivotal, indeed mutation in SUZ12 suppressed tumor growth in prostate [46].
The miR-17 is an oncomiR overexpressed in all stages of prostate cancer. The miR-17~92 cluster targets 3 0 UTR of tumor suppressor PTEN mRNA [47]. Studies suggest that miR-17~92 cluster promoter binds to E2F1-3 and activated its transcription in several cell types through a negative feedback loop [48]. The miR-17~92 expression is frequently upregulated in prostate cancer cells, resulting in E2F1-3 depletion thereby evading apoptosis. Besides miR-17 target genes altered during early-stage cancer include TSC1 and HEXIM1. In advance-stage miR-17 targets SMARCE1, VIM, AKT3, AGO2 and PLAG1; whereas TCEB1 and PTEN are major targets of miR-17 during CRPC. Studies suggest that miR-17 overexpression leads to downregulation of STAT3 expression, which in turn, suppresses proliferation in human prostate cancer LNCaP cells. STAT3 is a member of STAT protein family and a transcription factor phosphorylated by receptor-associated Janus kinases (JAK), playing a critical role in many biological processes including cellular proliferation. Let-7g miRNAs has been recognized as a tumor suppressor that plays an important role in various cellular processes and loss of its expression is associated with cancer development and progression [49]. Our previous publication has shown that the Let-7 family members are significantly expressed in both prostate epithelial and stromal cells [16]. Members of the let-7g miRNA targets and negatively regulates KRAS and MYC oncogenes, chromatin associated non-histone proteins viz. HMGA1, HMGA2, cell cycle regulator CDK6, and Lin28B that promotes cell proliferation, invasion and metastasis [50]. In fact, replacement of let-7g member's expression could be developed as a valid therapeutic option in certain types of tumor. Conclusively, the expression of above set of miRNAs including miR17, let-7g, miR-146, miR-204, miR-205, miR-301, miR-221 and miR520 remarkably modulate the expression of several target genes which include AR, ANXA1, BCL2, CDKN1B, EGFR, ESR1, MARK2, RUNX2, PTEN, SIRT1, SUZ12, VEGF-A, TSC1, HEXIM1and others in disease specific stage influencing various signaling pathways in proliferation, invasion, angiogenesis, migration and evasion of apoptosis.
miRNA-microarray expression data coupled with the IPA knowledge database identified associated signaling pathways and its interacting miRNAs such as miR-181a, miR130a, and miR-328 as predominant molecules associated with drug resistance drug efflux pathway in prostate cancer. The expression of 130a/b, miR-181a, (downregulated in early-stage, upregulated in advance-stage, and downregulated in CRPC) and miR-328 (upregulated in advancestage). Moreover, the drug resistance-drug efflux pathway showed that miR-130a/b, and miR-181a were associated with member's ABC transporter family including P-glycoprotein (P-gp/ ABCB1), and breast cancer resistant protein-(BCRP) (S2 Fig). Overexpression of ABC transporter proteins is associated with the resistance of a variety of chemotherapeutic drugs [51]. Likewise, miR-146a was differentially upregulated in PC3 cells and was associated with Th1 and Th2 signaling pathway (S2 Fig). Considerable evidence showed that miR-146a was a potent inhibitor of Th1 differentiation and cell proliferation of human T cells, and dysregulation of miR-146a contributed to the pathogenesis [52].
Another microRNA, miR-200b which was significantly downregulated in CRPC (S4 Table), was associated with epithelial-to-mesenchymal transition (EMT) pathway, and directly associated with Jagged1-Notch1 signaling (S2 Fig). Increasing evidences indicate that Notch1 expression was influenced by miR-200, and together they play an essential role in EMT progression. In consequence, EMT promotes cancer metastases and closely correlates with increased stemness and drug resistance. Moreover, the role of miR-200 warrant further investigation particularly during CRPC disease stage. Conclusively, the role of miRNA which includes miR-130a, miR-181, miR146 and miR-200 are considered to be important in context to prostate cancer.
Transcription factors (TFs) and miRNAs are essential regulators that directly or indirectly controls each other expression and effector genes through feedback and feed-forward loops [53]. TF and miRNA co-regulation is predominant in biological systems and any perturbation to this co-regulation alters the expression of upstream regulators. Earlier, inactivation/activation of upstream translational regulators has not received much attention in prostate cancer. In the present study the expression of AGO2, SSB, and NF2 is downregulated in early-stage cancer and upregulated in advance-stage cancer; whereas PPARA was overexpressed and DICER1 was inhibited across all stages of prostate cancer (Fig 3B). Documented evidence showed that higher expression of AGO2 was directly associated with prostate cancer cell proliferation, apoptosis, and cell cycle regulation [46]. Expression of miR-146 was upregulated only during advance-stage cancer, while miR-17 was upregulated, miR-221, miR-205, let-7g was downregulated across all disease stage (S1-S4 Tables). Consequently, AGO2 was inactivated during early-stage and CRPC and activated at advance-stage ( Fig 3B). Indeed, the interaction between miRNA: mRNA showed that miRNAs (miR-17, miR-205, miR-224, miR-let-7g, miR-10a and miR-130a) binds to the 3'UTR region of upstream regulators (AGO2, PPARA, NF2, DICER1 and SSB), which suppresses their posttranscriptional regulation either at the translation step of degrading the transcript by deadenylation. We successfully found target genes and pathways which have some degree of overlap between stage-specific disease progression (Figs 6 and 7; Table 1). However, the above facts and data warrant further functional research and validation either in preclinical model or established cell lines, which may lead to novel findings in context to prostate cancer, and biomarkers either prognosis or diagnosis.
In conclusion, miRNAs have emerged as essential modulators of immunity, cellular physiology and cancer. Our findings across prostate cancer cell lines identified potential candidate miRNAs involved in regulating target genes essential for the development and progression of prostate cancer. Most miRNAs have multiple targets at different stage of prostate cancer, thus it is not surprising that different targets may be involved in cancer cell apoptosis, cell differentiation and proliferation. Along with miRNAs, we have identified some potential upstream regulators which display critical role in prostate carcinogenesis, and further investigation on them may open a new avenue in prostate cancer research. We have identified a panel of consistently dysregulated miRNA in the present study and thus, validation of these promising candidate miRNAs in experimental studies and in a prospective cohort will not only define their exact role in prostate carcinogenesis but may also develop them as potential prognostic and diagnostic markers. We are currently analyzing the functional role of some miRNAs in prostate cancer. It is our understanding that further studies on these differentially expressed miR-NAs will lead to a better understanding of mechanisms mediating the development and progression of prostate cancer. The internetworking relationship with miRNAs and plasma membrane protein P-glycoprotein (Pgp-plasma membrane protein). The miRNA-130a and miR-181a are downregulated in this pathway (green color) and linked with Pg and BCRP (breast cancer resistant protein). miR-133a and miR379 are involved in regulating the expression of MRP2. While miR-298, miR27a, miR331-5p and miR-130a are involved in regulating the expression of poly-glycoprotein (P-gp). Epithelial mesenchymal transition pathway. The miR-200b was downregulated during early-stage prostate cancer and was involved in inhibiting Jagged-2 (JAG2), one of the NOTCH ligands. Adipogenesis pathway. The expression of miR-326 was upregulated during the metastatic stage of prostate cancer, (red color) and involved in regulating the expression of CCAT/enhancer binding protein β (C/EBPβ), directly linked with the nuclear hormone receptor peroxisome proliferator-activated receptor-gamma (PPAR-γ). Bone metamorphosis signaling pathway. The expression of miR-140, miR-145, and miR-155 was upregulated, along with miR-140 and miR-145 were associated with modulation of gene SOX9, and miR-155 was involved in regulating the gene FOXO3A. Th1 pathway. In this pathway, the expression of miR-146a showed a higher level of expression (red color) and modulated the expression of NF-κB signaling. Th1 and Th2 pathway. In this pathway, the expression of miR-146a showed a lower level of expression (green color) and may modulate the expression of NF-κB signaling. (PDF) S1