Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype.

Molecular heterogeneity in prostate cancer (PCa) is one of the key reasons underlying the differing likelihoods of recurrence after surgical treatment in individual patients of the same clinical category. In this study, we performed RNA-Seq profiling of 58 localized PCa and 43 locally advanced PCa tissue samples obtained as a result of radical prostatectomy on a cohort of Russian patients. Based on bioinformatics analysis, we examined features of the transcriptome profiles within the high-risk group, including within the most commonly represented molecular subtype, TMPRSS2-ERG. The most significantly affected biological processes in the samples were also identified, so that they may be further studied in the search for new potential therapeutic targets for the categories of PCa under consideration. The highest predictive potential was found with the EEF1A1P5, RPLP0P6, ZNF483, CIBAR1, HECTD2, OGN, and CLIC4 genes. We also reviewed the main transcriptome changes in the groups at intermediate risk of PCa-Gleason Score 7 (groups 2 and 3 according to the ISUP classification)-on the basis of which the LPL, MYC, and TWIST1 genes were identified as promising additional prognostic markers, the statistical significance of which was confirmed using qPCR validation.

Lymphatic Dissemination in Prostate Cancer: Features of the Transcriptomic Profile and Prognostic Models.

Radical prostatectomy is the gold standard treatment for prostate cancer (PCa); however, it does not always completely cure PCa, and patients often experience a recurrence of the disease. In addition, the clinical and pathological parameters used to assess the prognosis and choose further tactics for treating a patient are insufficiently informative and need to be supplemented with new markers. In this study, we performed RNA-Seq of PCa tissue samples, aimed at identifying potential prognostic markers at the level of gene expression and miRNAs associated with one of the key signs of cancer aggressiveness-lymphatic dissemination. The relative expression of candidate markers was validated by quantitative PCR, including an independent sample of patients based on archival material. Statistically significant results, derived from an independent set of samples, were confirmed for miR-148a-3p and miR-615-3p, as well as for the CST2, OCLN, and PCAT4 genes. Considering the obtained validation data, we also analyzed the predictive value of models based on various combinations of identified markers using algorithms based on machine learning. The highest predictive potential was shown for the "CST2 + OCLN + pT" model (AUC = 0.863) based on the CatBoost Classifier algorithm.

Docetaxel Resistance in Castration-Resistant Prostate Cancer: Transcriptomic Determinants and the Effect of Inhibiting Wnt/ß-Catenin Signaling by XAV939.

Castration-resistant prostate cancer (CRPC) is a common form of prostate cancer in which docetaxel-based chemotherapy is used as the first line. The present study is devoted to the analysis of transcriptome profiles of tumor cells in the development of resistance to docetaxel as well as to the assessment of the combined effect with the XAV939 tankyrase inhibitor on maintaining the sensitivity of tumor cells to chemotherapy. RNA-Seq was performed for experimental PC3 cell lines as well as for plasma exosome samples from patients with CRPC. We have identified key biological processes and identified a signature based on the expression of 17 mRNA isoforms associated with the development of docetaxel resistance in PC3 cells. Transcripts were found in exosome samples, the increased expression of which was associated with the onset of progression of CRPC during therapy. The suppression of pathways associated with the participation of cellular microtubules has also been shown when cells are treated with docetaxel in the presence of XAV939. These results highlight the importance of further research into XAV939 as a therapeutic agent in the treatment of CRPC; moreover, we have proposed a number of mRNA isoforms with high predictive potential, which can be considered as promising markers of response to docetaxel.

ALDH3A2, ODF2, QSOX2, and MicroRNA-503-5p Expression to Forecast Recurrence in TMPRSS2-ERG-Positive Prostate Cancer.

Following radical surgery, patients may suffer a relapse. It is important to identify such patients so that therapy tactics can be modified appropriately. Existing stratification schemes do not display the probability of recurrence with enough precision since locally advanced prostate cancer (PCa) is classified as high-risk but is not ranked in greater detail. Between 40 and 50% of PCa cases belong to the TMPRSS2-ERG subtype that is a sufficiently homogeneous group for high-precision prognostic marker search to be possible. This study includes two independent cohorts and is based on high throughput sequencing and qPCR data. As a result, we have been able to suggest a perspective-trained model involving a deep neural network based on both qPCR data for mRNA and miRNA and clinicopathological criteria that can be used for recurrence risk forecasts in patients with TMPRSS2-ERG-positive, locally advanced PCa (the model uses ALDH3A2 + ODF2 + QSOX2 + hsa-miR-503-5p + ISUP + pT, with an AUC = 0.944). In addition to the prognostic model's use of identified differentially expressed genes and miRNAs, miRNA-target pairs were found that correlate with the prognosis and can be presented as an interactome network.

HK3 overexpression associated with epithelial-mesenchymal transition in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is a common cancer worldwide. The main cause of death in CRC includes tumor progression and metastasis. At molecular level, these processes may be triggered by epithelial-mesenchymal transition (EMT) and necessitates specific alterations in cell metabolism. Although several EMT-related metabolic changes have been described in CRC, the mechanism is still poorly understood. RESULTS: Using CrossHub software, we analyzed RNA-Seq expression profile data of CRC derived from The Cancer Genome Atlas (TCGA) project. Correlation analysis between the change in the expression of genes involved in glycolysis and EMT was performed. We obtained the set of genes with significant correlation coefficients, which included 21 EMT-related genes and a single glycolytic gene, HK3. The mRNA level of these genes was measured in 78 paired colorectal cancer samples by quantitative polymerase chain reaction (qPCR). Upregulation of HK3 and deregulation of 11 genes (COL1A1, TWIST1, NFATC1, GLIPR2, SFPR1, FLNA, GREM1, SFRP2, ZEB2, SPP1, and RARRES1) involved in EMT were found. The results of correlation study showed that the expression of HK3 demonstrated a strong correlation with 7 of the 21 examined genes (ZEB2, GREM1, TGFB3, TGFB1, SNAI2, TWIST1, and COL1A1) in CRC. CONCLUSIONS: Upregulation of HK3 is associated with EMT in CRC and may be a crucial metabolic adaptation for rapid proliferation, survival, and metastases of CRC cells.

Upregulation of NETO2 gene in colorectal cancer.

BACKGROUND: Neuropilin and tolloid-like 2 (NETO2) is a single-pass transmembrane protein that has been shown primarily implicated in neuron-specific processes. Upregulation of NETO2 gene was also detected in several cancer types. In colorectal cancer (CRC), it was associated with tumor progression, invasion, and metastasis, and seems to be involved in epithelial-mesenchymal transition (EMT). However, the mechanism of NETO2 action is still poorly understood. RESULTS: We have revealed significant increase in the expression of NETO2 gene and deregulation of eight EMT-related genes in CRC. Four of them were upregulated (TWIST1, SNAIL1, LEF1, and FOXA2); the mRNA levels of other genes (FOXA1, BMP2, BMP5, and SMAD7) were decreased. Expression of NETO2 gene was weakly correlated with that of genes involved in the EMT process. CONCLUSIONS: We found considerable NETO2 upregulation, but no significant correlation between the expression of NETO2 and EMT-related genes in CRC. Thus, NETO2 may be involved in CRC progression, but is not directly associated with EMT.

Differential expression of alternatively spliced transcripts related to energy metabolism in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is one of the most common malignant tumors worldwide. CRC molecular pathogenesis is heterogeneous and may be followed by mutations in oncogenes and tumor suppressor genes, chromosomal and microsatellite instability, alternative splicing alterations, hypermethylation of CpG islands, oxidative stress, impairment of different signaling pathways and energy metabolism. In the present work, we have studied the alterations of alternative splicing patterns of genes related to energy metabolism in CRC. RESULTS: Using CrossHub software, we analyzed The Cancer Genome Atlas (TCGA) RNA-Seq datasets derived from colon tumor and matched normal tissues. The expression of 1014 alternative mRNA isoforms involved in cell energy metabolism was examined. We found 7 genes with differentially expressed alternative transcripts whereas overall expression of these genes was not significantly altered in CRC. A set of 8 differentially expressed transcripts of interest has been validated by qPCR. These eight isoforms encoded by OGDH, COL6A3, ICAM1, PHPT1, PPP2R5D, SLC29A1, and TRIB3 genes were up-regulated in colorectal tumors, and this is in concordance with the bioinformatics data. The alternative transcript NM_057167 of COL6A3 was also strongly up-regulated in breast, lung, prostate, and kidney tumors. Alternative transcript of SLC29A1 (NM_001078177) was up-regulated only in CRC samples, but not in the other tested tumor types. CONCLUSIONS: We identified tumor-specific expression of alternative spliced transcripts of seven genes involved in energy metabolism in CRC. Our results bring new knowledge on alternative splicing in colorectal cancer and suggest a set of mRNA isoforms that could be used for cancer diagnosis and development of treatment methods.

Effect of lentivirus-mediated shRNA inactivation of HK1, HK2, and HK3 genes in colorectal cancer and melanoma cells.

BACKGROUND: The switch from oxidative phosphorylation to glycolysis in proliferating cancer cells, even under aerobic conditions, has been shown first in 1926 by Otto Warburg. Today this phenomenon is known as the "Warburg effect" and recognized as a hallmark of cancer. The metabolic shift to glycolysis is associated with the alterations in signaling pathways involved in energy metabolism, including glucose uptake and fermentation, and regulation of mitochondrial functions. Hexokinases (HKs), which catalyze the first step of glycolysis, have been identified to play a role in tumorigenesis of human colorectal cancer (CRC) and melanoma. However, the mechanism of action of HKs in the promotion of tumor growth remains unclear. RESULTS: The purpose of the present study was to investigate the effect of silencing of hexokinase genes (HK1, HK2, and HK3) in colorectal cancer (HT-29, SW 480, HCT-15, RKO, and HCT 116) and melanoma (MDA-MB-435S and SK-MEL-28) cell lines using short hairpin RNA (shRNA) lentiviral vectors. shRNA lentiviral plasmid vectors pLSLP-HK1, pLSLP-HK2, and pLSLP-HK3 were constructed and then transfected separately or co-transfected into the cells. HK2 inactivation was associated with increased expression of HK1 in colorectal cancer cell lines pointing to the compensation effect. Simultaneous attenuation of HK1 and HK2 levels led to decreased cell viability. Co-transfection with shRNA vectors against HK1, HK2, and HK3 mRNAs resulted in a rapid cell death via apoptosis. CONCLUSIONS: We have demonstrated that simultaneous inactivation of HK1 and HK2 was sufficient to decrease proliferation and viability of melanoma and colorectal cancer cells. Our results suggest that HK1 and HK2 could be the key therapeutic targets for reducing aerobic glycolysis in examined cancers.

Dynamic Profiling of Exosomal microRNAs in Blood Plasma of Patients with Castration-Resistant Prostate Cancer.

Prostate cancer is one of the most common and socially significant cancers among men. The aim of this study was to identify significant changes in the expression of exosomal miRNAs associated with an increase in the level of prostate specific antigen in castration-resistant prostate cancer during therapy and to evaluate them as potential prognostic markers for this category of disease. High-throughput miRNA sequencing was performed on 49 blood plasma samples taken from 11 Russian patients with castration-resistant cancer during therapy. Bioinformatic analysis of the obtained miRNA-seq data was carried out. Additionally, miRNA-seq data from the PRJNA562276 project were analyzed to identify exosomal miRNAs associated with castration-resistant prostate cancer. We found 34 differentially expressed miRNAs associated with the progression of castration-resistant prostate cancer during therapy in Russian patients. It was also shown that hsa-miRNA-148a-3p expression can serve as a potential prognostic marker. We found the exosomal miRNA expression signature associated with castration-resistant prostate cancer progression, in particular on the Russian patient cohort. Many of these miRNAs are well-known players in either oncogenic transformation or tumor suppression. Further experimental studies with extended sampling are required to validate these results.

Impact TMPRSS2-ERG Molecular Subtype on Prostate Cancer Recurrence.

Currently, seven molecular subtypes of prostate cancer (PCa) are known, the most common of which being the subtype characterized by the presence of the TMPRSS2-ERG fusion transcript. While there is a considerable amount of work devoted to the influence of this transcript on the prognosis of the disease, data on its role in the progression and prognosis of PCa remain controversial. The present study is devoted to the analysis of the association between the TMPRSS2-ERG transcript and the biochemical recurrence of PCa. The study included two cohorts: the RNA-Seq sample of Russian patients with PCa (n = 72) and the TCGA-PRAD data (n = 203). The results of the analysis of the association between the TMPRSS2-ERG transcript and biochemical recurrence were contradictory. The differential expression analysis (biochemical recurrence cases versus biochemical recurrence-free) and the gene set enrichment analysis revealed a list of genes involved in major cellular pathways. The GNL3, QSOX2, SSPO, and SYS1 genes were selected as predictors of the potential prognostic model (AUC = 1.000 for a cohort of Russian patients with PCa and AUC = 0.779 for a TCGA-PRAD cohort).

Gene Expression Changes and Associated Pathways Involved in the Progression of Prostate Cancer Advanced Stages.

Prostate cancer (PC) is one of the most common cancers among men worldwide, and advanced PCs, such as locally advanced PC (LAPC) and castration-resistant PC (CRPC), present the greatest challenges in clinical management. Current indicators have limited capacity to predict the disease course; therefore, better prognostic markers are greatly needed. In this study, we performed a bioinformatic analysis of The Cancer Genome Atlas (TCGA) datasets, including RNA-Seq data from the prostate adenocarcinoma (PRAD; n = 55) and West Coast Dream Team - metastatic CRPC (WCDT-MCRPC; n = 84) projects, to evaluate the transcriptome changes associated with progression-free survival (PFS) for LAPC and CRPC, respectively. We identified the genes whose expression was positively/negatively correlated with PFS. In LAPC, the genes with the most significant negative correlations were ZC2HC1A, SQLE, and KIF11, and the genes with the most significant positive correlations were SOD3, LRRC26, MIR22HG, MEG3, and MIR29B2CHG. In CRPC, the most significant positive correlations were found for BET1, CTAGE5, IFNGR1, and GIMAP6, and the most significant negative correlations were found for CLPB, PRPF19, ZNF610, MPST, and LINC02001. In addition, we performed a gene network interaction analysis using STRINGdb, which revealed a significant relationship between genes predominantly involved in the cell cycle and characterized by upregulated expression in early recurrence. Based on the results, we propose several genes that can be used as potential prognostic markers.

miRNAs expression signature potentially associated with lymphatic dissemination in locally advanced prostate cancer.

BACKGROUND: Prostate cancer is one of the most common and socially significant cancers among men. The aim of our study was to reveal changes in miRNA expression profiles associated with lymphatic dissemination in prostate cancer and to identify the most prominent miRNAs as potential prognostic markers for future studies. METHODS: High-throughput miRNA sequencing was performed for 44 prostate cancer specimens taken from Russian patients, with and without lymphatic dissemination (N1 - 20 samples; N0 - 24 samples). RESULTS: We found at least 18 microRNAs with differential expression between N0 and N1 sample groups: miR-182-5p, miR-183-5p, miR-96-5p, miR-25-3p, miR-93-5p, miR-7-5p, miR-615-3p, miR-10b, miR-1248 (N1-miRs; elevated expression in N1 cohort; p < 0.05); miR-1271-5p, miR-184, miR-222-3p, miR-221-5p, miR-221-3p, miR-455-3p, miR-143-5p, miR-181c-3p and miR-455-5p (N0-miRs; elevated expression in N0; p < 0.05). The expression levels of N1-miRs were highly correlated between each other (the same is applied for N0-miRs) and the expression levels of N0-miRs and N1-miRs were anti-correlated. The tumor samples can be divided into two groups depending on the expression ratio between N0-miRs and N1-miRs. CONCLUSIONS: We found the miRNA expression signature associated with lymphatic dissemination, in particular on the Russian patient cohort. Many of these miRNAs are well-known players in either oncogenic transformation or tumor suppression. Further experimental studies with extended sampling are required to validate these results.

Bioinformatic identification of differentially expressed genes associated with prognosis of locally advanced lymph node-positive prostate cancer.

Prostate cancer (PCa) is one of the primary causes of cancer-related mortality in men worldwide. Patients with locally advanced PCa with metastases in regional lymph nodes are usually marked as a high-risk group. One of the chief concerns for this group is to make an informed decision about the necessity of conducting adjuvant androgen deprivation therapy after radical surgical treatment. During the oncogenic transformation and progression of the disease, the expression of many genes is altered. Some of these genes can serve as markers for diagnosis, predicting the prognosis or effectiveness of drug therapy, as well as possible therapeutic targets. We undertook bioinformatic analysis of the RNA-seq data deposited in The Cancer Genome Atlas consortium database to identify possible prognostic markers. We compared the groups with favorable and unfavorable prognosis for the cohort of patients with PCa showing lymph node metastasis (pT2N1M0, pT3N1M0, and pT4N1M0) and for the most common molecular type carrying the fusion transcript TMPRSS2-ERG. For the entire cohort, we revealed at least six potential markers (IDO1, UGT2B15, IFNG, MUC6, CXCL11, and GBP1). Most of these genes are involved in the positive regulation of immune response. For the TMPRSS2-ERG subtype, we also identified six genes, the expression of which may be associated with prognosis: TOB1, GALNT7, INAFM1, APELA, RAC3, and NNMT. The identified genes, after additional studies and validation in the extended cohort, could serve as a prognostic marker of locally advanced lymph node-positive PCa.

Differentially Expressed Genes Associated With Prognosis in Locally Advanced Lymph Node-Negative Prostate Cancer.

Older age is one of the main risk factors for cancer development. The incidence of prostate cancer, as a multifactorial disease, also depends upon demographic factors, race, and genetic predisposition. Prostate cancer most frequently occurs in men over 60 years of age, indicating a clear association between older age and disease onset. Carcinogenesis is followed by the deregulation of many genes, and some of these changes could serve as biomarkers for diagnosis, prognosis, prediction of drug therapy efficacy, as well as possible therapeutic targets. We have performed a bioinformatic analysis of a The Cancer Genome Atlas (TCGA) data and RNA-Seq profiling of a Russian patient cohort to reveal prognostic markers of locally advanced lymph node-negative prostate cancer (lymph node-negative LAPC). We also aimed to identify markers of the most common molecular subtype of prostate cancer carrying a fusion transcript TMPRSS2-ERG. We have found several genes that were differently expressed between the favorable and unfavorable prognosis groups and involved in the enriched KEGG pathways based on the TCGA (B4GALNT4, PTK6, and CHAT) and Russian patient cohort data (AKR1C1 and AKR1C3). Additionally, we revealed such genes for the TMPRSS2-ERG prostate cancer molecular subtype (B4GALNT4, ASRGL1, MYBPC1, RGS11, SLC6A14, GALNT13, and ST6GALNAC1). Obtained results contribute to a better understanding of the molecular mechanisms behind prostate cancer progression and could be used for further development of the LAPC prognosis marker panel.