Next Article in Journal
Faujasite-Type Zeolite Obtained from Ecuadorian Clay as a Support of ZnTiO3/TiO2 NPs for Cyanide Removal in Aqueous Solutions
Next Article in Special Issue
In Silico Identification and Analysis of Proteins Containing the Phox Homology Phosphoinositide-Binding Domain in Kinetoplastea Protists: Evolutionary Conservation and Uniqueness of Phox-Homology-Domain-Containing Protein Architectures
Previous Article in Journal
Expression Analysis of Lipocalin 2 (LCN2) in Reproductive and Non-Reproductive Tissues of Esr1-Deficient Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Anastasiya A. Kobelyatskaya
1,*,†,
Elena A. Pudova
1,†,
Irina V. Katunina
1,
Anastasiya V. Snezhkina
1,
Maria S. Fedorova
1,
Vladislav S. Pavlov
1,
Anastasiya O. Kotelnikova
1,
Kirill M. Nyushko
2,
Boris Y. Alekseev
2,
George S. Krasnov
1 and
Anna V. Kudryavtseva
1,*
1
Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia
2
National Medical Research Radiological Center, Ministry of Health of the Russian Federation, 125284 Moscow, Russia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(11), 9282; https://doi.org/10.3390/ijms24119282
Submission received: 7 April 2023 / Revised: 12 May 2023 / Accepted: 24 May 2023 / Published: 25 May 2023
(This article belongs to the Special Issue Essential Molecules in Life: Regulation, Defense, and Longevity)

Abstract

:
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.

1. Introduction

Prostate cancer (PCa) is one of the most common types of cancer among men worldwide [1]. The majority of PCa cases are diagnosed as having a localized form, which represents the early malignant process confined to the prostate gland without spreading beyond its borders. In the early stages of localized PCa (LPCa), symptoms may include minor changes in the urinary system or even no symptoms at all. The standard clinical diagnosis of LPCa includes various methods such as prostate palpation, measurement of levels of prostate-specific antigen (PSA), ultrasound/magnetic resonance imaging, and biopsy.
LPCa often has a favorable prognosis, and there is a wide range of available therapeutic approaches (active surveillance, radiation therapy, radical prostatectomy, or their combinations) that can be chosen based on the patient’s risk stratification for biochemical recurrence. One of the main systems for classifying risk groups in PCa is the D’Amico classification system, which defines three risk groups: low, intermediate, and high. This risk group assessment is based on three factors: the Gleason score (GS), PSA level, and stage of the disease. The low-risk group is characterized by a GS ≤ 6, PSA level < 10 ng/mL, and T1–T2a stage; intermediate risk by GS = 7, PSA level 10–20 ng/mL, and/or T2b–T2c stage; and high risk by GS ≥ 8, PSA level > 20 ng/mL, and/or stage ≥ T3a [2].
Patient stratification into risk groups is a widely used approach for assessing the prognosis of PCa. However, in some cases, this method may not be sufficiently informative, potentially leading to incorrect conclusions and inappropriate treatments. One of the limitations of patient stratification into such risk groups is that this method does not consider many other factors that may affect the disease prognosis, such as age, comorbidities, and genetics. Studies have shown that even patients with low D’Amico risk may have a poor prognosis if they have other risk factors that were not considered during stratification [3].
Thus, stratification of patients into risk groups is a useful tool for assessing PCa prognosis, but other risk factors and individual patient characteristics must also be considered.
One of the main factors complicating the diagnosis and treatment of PCa, as well as of other types of cancer, is molecular heterogeneity. This phenomenon is defined by differences in the molecular properties of tumor cells, such as changes in gene expression, mutations, and other factors that affect cell function and behavior [4,5]. Furthermore, the development of aggressive forms of PCa requires only a few driver alterations [6,7,8]. Molecular heterogeneity is also present within LPCa, and this can lead to different prognosis and outcomes [4,5].
It has been shown that 74% of all PCa cases can be attributed to one of a range of molecular subtypes that have been identified on the basis of analysis of somatic mutations, changes in copy number, gene expression, gene fusions, and DNA methylation [9].
There are currently seven major molecular subtypes of PCa that have been identified as part of the Prostate Adenocarcinoma Project of The Cancer Genome Atlas (TCGA) consortium [9]. Four of the seven subtypes are characterized by the presence of fusion transcripts between the TMPRSS2 gene exons and the exons of the ETS family of genes (the erythroblast transformation-specific family of transcription factors): ERG, ETV1, ETV4, and FLI1 (the frequency of occurrence of these subtypes is 46%, 8%, 4%, and 1%, respectively). Three other subtypes are characterized by the presence of point mutations in one of the following genes: SPOP, FOXA1, or IDH1 (the frequency of occurrence of these subtypes is 11%, 3%, and 1%, respectively) [9,10]. Thus, about half of all cases of prostate cancer have a fused TMPRSS2-ERG transcript, which is formed due to an intrachromosomal rearrangement leading to the fusion of two genes: TMPRSS2 and ERG.
Considering the molecular heterogeneity of PCa in particular, the selection of the group of tumors characterized by the TMPRSS2-ERG fusion transcript allows researchers to focus on specific molecular subtypes of PCa and, therefore, form a more homogeneous group of PCa cases. To improve prognoses and to determine the best treatment approach for each patient with LPCa, a more detailed study of tumor molecular heterogeneity is necessary. New biomarkers and genomic technologies can help in this direction and allow for the identification of the biological nature of any given PCa and its prognosis [11].
This study involved RNA-Seq profiling of 58 LPCa samples from intermediate- and high-risk groups and 43 locally advanced PCa (LAPCa) samples from a cohort of Russian patients, with the aim of identifying transcriptional profile differences between these groups, taking into account the TMPRSS2-ERG subtype, and searching for promising genes as prognostic markers.

2. Results

2.1. Differentially Expressed Genes among Risk Groups within Localized PCa

LPCa can be classified into one of three risk groups: low, intermediate, and high. Based on the LPCa cases we analyzed, the low-risk group is quite rare (about 2%), whereas the intermediate one is the most frequent (80–82%), and the high-risk group accounts for 16–18%.
Initially, we compared primary PCa tumors belonging to the intermediate- (n = 47) and high-risk (n = 10) groups. Based on RNA-Seq data obtained for PCas in Russian patients having tumors of the high and intermediate LPCa risk groups, just six (↓COL17A1, ↑FAM83D, ↑GCSAML, ↓IER2, ↑IFI44L, ↓MYH3) differentially expressed (DE) genes were identified (p-value ≤ 0.05 according to the QLF and MW tests, Log2CPM ≥ 3, |Log2FC| ≥ 1, Supplementary Table S1-1). This is definitely too small a number of DE genes; apparently, these risk groups do not have extensive transcriptomic variations among themselves.

2.2. Differentially Expressed Genes between LAPCa and LPCa for High-Risk Group

At this stage of our work, besides the LPCa samples, we included LAPCa cases in the study, in order to find transcriptome differences between these grades of tumor extension.
Including the expression data of LPCa cases, as well as the data obtained in our previous work devoted to LAPCa samples, we performed differential expression analysis (DEA) between the groups of LPCa (n = 10) and LAPCa (n = 43) samples classified as high-risk. As a result, 243 DE genes were identified (p-values ≤ 0.05 according to the QLF and MW tests, Log2CPM ≥ 3, |Log2FC| ≥ 1, Supplementary Table S1-2). Figure 1 shows the expression profiles of the top 50 DE genes.
Based on the identified DE gene profile, we further analyzed the enrichment of biological pathways associated with LAPCa versus LPCa within the high-risk group. As a result of the analysis based on the GSEA algorithm and the KEGG Human 2021 database, we found statistically significant changes in 68 biological processes (FDR ≤ 0.05, Figure 2, Supplementary Table S2-1).
The identified biological pathways include many processes with known involvement in the development and progression of various types of tumors, including “Non-small cell lung cancer”, “Thyroid cancer”, and “Melanoma”, as well as the effects of “Proteoglycans in cancer” and the “Ras signaling pathway”. A complete list of analysis results is presented in Supplementary Table S2-1.

2.3. Differentially Expressed Genes between LAPCa and LPCa within the TMPRSS2-ERG Molecular Subtype

It is known that the frequency of occurrence of the TMPRSS2-ERG subtype varies from 40% to 50% in PCa [9,10]. In addition, most researchers consider TMPRSS2-ERG as a factor involved in increased aggressiveness, propensity for invasion, and metastasis [12,13,14]. However, several studies have demonstrated that TMPRSS2-ERG is either a precursor of a good prognosis or has no association with progression and prognosis at all [15,16]. In the present study, the incidence of the TMPRSS2-ERG subtype was about 42% for LAPCa and 28% for LPCa. Furthermore, we searched for associations of the TMPRSS2-ERG subtype with clinical and pathomorphological criteria. Spearman’s correlation analysis did not show a significant association of the TMPRSS2-ERG subtype with any of the following criteria: age, tumor extension groups (LAPCa, LPCa), risk groups, Gleason Score, ISUP, pT, or preoperative PSA value.
Insomuch as the TMPRSS2-ERG subtype results in a more homogeneous group, we performed DEA between the groups of LPCa (n = 16) and LAPCa (n = 18) samples only within the TMPRSS2-ERG subtype.
As a result, 207 DE genes were identified (p-values ≤ 0.05 according to the QLF and MW tests, Log2CPM ≥ 3, |Log2FC| ≥ 1, Supplementary Table S1-3). Figure 3 shows the expression profiles of the top 50 of these DE genes.
Based on the identified DE genes, we also performed an enrichment analysis of a number of biological pathways. As a result of the analysis, we found statistically significant changes in 16 biological processes (FDR ≤ 0.05, Figure 4, Supplementary Table S2-2).
After clarifying the TMPRSS2-ERG molecular subtype in the LAPCa category, various signaling processes that have a known involvement in PCa progression, such as the “TGF-beta signaling pathway”, were seen to become the most significant. A complete list of analysis results is presented in Supplementary Table S2-2.
We also identified DE genes whose expression was statistically significantly associated with the LAPCa group within the TMPRSS2-ERG subtype: BHLHA15, CIBAR1, CLIC4, CORO1B, CRB3, DNAJB4, DNM3OS, EEF1A1P5, HECTD2, ID4, MFSD3, MIR222HG, OGN, RPLP0P6, SH3BGRL, and ZNF483. The results of the differential expression of these genes are presented in Table 1.
It should be noted that, according to the differential expression of the CIBAR1, CLIC4, DNAJB4, and EEF1A1P5 genes within the TMPRSS2-ERG molecular subtype, we observed the highest correlation of these with the LAPCa group.
We also considered the predictive potential of selected genes by analyzing ROC curves based on a logistic regression algorithm. The results of the analysis are presented in Table 2.
According to the results obtained, the EEF1A1P5 and RPLP0P6 genes had the highest values (AUC > 0.9) of the AUC metric in the test data, both when analyzed in the total set of samples and within the TMPRSS2-ERG molecular subtype. It is also worth noting that an increase in value of more than 0.9 for the AUC metrics after clarification of the molecular subtype was found in the CIBAR1, CLIC4, HECTD2, OGN, and ZNF483 genes.

2.4. Differentially Expressed Genes Associated with ISUP 3 at Intermediate Risk for Localized PCa

Based on the RNA-Seq data obtained for ISUP = 3 (n = 17) and ISUP = 2 (n = 10) LPCa in our Russian patient cohort within the GS = 7, 36 differentially expressed (DE) genes were identified (p-values ≤ 0.05 according to the QLF and MW tests, Log2CPM ≥ 3, |Log2FC| ≥ 1, Supplementary Table S1-4). Figure 5 shows the expression profiles of these DE genes.
Based on the results of our analysis of the enrichment of biological pathways associated with the ISUP 3 group in the LPCa of the intermediate-risk group, we found statistically significant changes in only one cancer-associated biological process, the “PPAR signaling pathway” (NES = 1.93; FDR = 0.02, Figure 6). A complete list of analysis results is presented in Supplementary Table S2-3.

2.5. Validation of the Relative Expression of Genes Associated with the ISUP 3 Group in the Intermediate-Risk LPCa Group of Russian Patients

We also carried out selection and validation of the relative expression of promising genes associated with the ISUP 3 group in our cohort of patients, in order that they might later be considered as additional prognostic markers. Based on the most significant differential expression results, the LPL, MYC, and TWIST1 genes were selected for validation by qPCR. According to the results of this validation of the relative expressions of the genes under consideration, statistically significant results were confirmed (Figure 7).

3. Discussion

In the current work, we performed a comprehensive study of LPCa, considering the risk groups, the degree of differentiation of the tumor cells (ISUP classification), and inclusion in the TMPRSS2-ERG molecular subtype based on RNA-Seq profiling.
As our results demonstrate, only the PCa cases belonging to high- and intermediate-risk groups have insufficiently stable transcriptomic variations between themselves. We found changes in only six genes (↓COL17A1, ↑FAM83D, ↑GCSAML, ↓IER2, ↑IFI44L, ↓MYH3), but little is known about their involvement in PCa. However, when considering two grades of tumor extension (LPCa and LAPCa), we were able to find significant differences in gene expression. It is necessary to mention that when considering the TMPRSS2-ERG molecular subtype (i.e., when only TMPRSS2-ERG-positive cases are included in the analysis), the significance of the changes in the expression of the previously detected genes (CIBAR1, CLIC4, EEF1A1P5, OGN, RPLP0P6, and ZNF483) increased.
The CIBAR1, CLIC4, OGN, and ZNF483 are protein-coding genes, but nothing is known about their association with PCa and cancer in general. EEF1A1P5 and RPLP0P6 are pseudogenes, and there are currently no data on their association with PCa progression. However, there is evidence that EEF1A1P5 gene transcripts and the RPLP0P6 protein are present in exosomes of various tumor cell lines [17,18].
Regarding the HECTD2 gene in PCa, it has been shown that a decrease in the expression of this gene significantly affects androgen-induced and AR-mediated transcription, while suppression of HECTD2 also enhances the growth of LNCaP cells [19]. According to our data, we also observed a decrease in the expression of the HECTD2 gene in the more advanced stage of the high-risk group—LAPCa.
When considering the most significantly enhanced cancer-associated biological processes in the LAPCa category within the TMPRSS2-ERG molecular subtype, such signaling pathways as the “cAMP signaling pathway” and the “TGF-beta signaling pathway” come to the fore. Aberrant signaling in these pathways has been implicated in various types of tumor. Transforming growth factor β (TGF-β) is a key regulator of many biological processes, including metastasis and invasion. This protein binds to intracellular receptors and activates a signaling pathway that is involved in the regulation of cell proliferation and differentiation. Multiple studies have shown that such TGF-β signaling is associated with poor prognosis in PCa [20,21,22]. cAMP signaling can play both a tumor-suppressing and tumor-promoting role, depending on the type of tumor. This cascade can also be used to regulate the growth, migration, invasion, and metabolism of cancer cells, including those in PCa [23,24]. Thus, our results indicate that the identified signaling pathways play an important role in PCa and can potentially be used to assess likely prognoses. Further research is needed to better understand the mechanisms of action of these signaling pathways at different stages of PCa progression and their potential value in developing more effective treatments.
Our study also considered another clinical problem—that based on the molecular heterogeneity of prostate tumors. According to the D’Amico classification for LPCa, patients with PSA from 10 to 20 ng/mL, a Gleason score (GS) of 7 (ISUP 2/3), and cT2b belong to the intermediate-risk group. The standard of care for patients in this group is radical prostatectomy with/without extended pelvic lymphadenectomy or external beam radiation therapy.
This group of patients is one of the most heterogeneous, including both patients with GS (3+4) and GS (4+3), as well as a variety of PSA levels. The main difficulty in the treatment of patients in this risk group is the high probability of disease progression after radical treatment. The recommendations of the American Association of Urology (2017) propose a division of this intermediate-risk group into favorable (GS 3+4, risk group 2 (ISUP 2)) and unfavorable (GS 4+3, risk group 3 (ISUP 3)). Unfortunately, this division does not allow a qualitative change in the approach to the treatment of patients in this category, and the main approach remains radical prostatectomy or radiation therapy with androgen deprivation therapy.
We investigated the features of the transcriptome profile in our LPCa patients in the intermediate-risk group, based on the results of differential gene expression between the ISUP 2 and 3 groups.
From our analysis of biological pathway enrichment, it was found that the ISUP 3 group was characterized by a statistically significant increase in the regulation of the PPAR signaling pathway (NES = 1.93; FDR = 0.015). This signaling cascade is one of the most important mechanisms for regulating lipid and glucose metabolism, as well as cell growth and differentiation. The main participants in this pathway, according to our sample of patients, are the genes ADIPOQ (LogFC = 6.92; MW p-value = 0.04), FABP4 (LogFC = 4.72; MW p-value = 0.03), and LPL (LogFC = 3.45; MW p-value = 0.01).
Based on the identified profiles of the DE genes, we selected the most significant of them in terms of statistical and expression metrics for subsequent qPCR validation: LPL, MYC, and TWIST1. Statistically significant differences between the ISUP 2 and ISUP 3 groups were shown, based on the relative expression of all the selected genes (p-value ≤ 0.05).
The TWIST1 and MYC genes are well-known and important regulators of cancer-associated processes and remain objects of active research in the field of oncology. Various studies have shown a relationship between high levels of expression of these genes and the aggressiveness of oncological diseases, including PCa [25,26,27,28]. Based on our results with Russian patients, we have also demonstrated that increased expression of the TWIST1 and MYC genes is associated with an unfavorable intermediate risk in LPCa.
However, the LPL gene, which encodes the enzyme lipoprotein lipase and plays a key role in the metabolism of fats and carbohydrates, is of particular interest [29]. LPL is an important enzyme that is secreted by extracellular lipolysis and can potentially be supplied by tumor cells or adjacent adipose tissue cells into the tumor microenvironment [30].
The main function of this enzyme is to hydrolyze triglycerides, resulting in the release of fatty acids, essential building blocks of biological membranes, so such dysregulation of fatty acid metabolism is a vital component of lipid metabolism reprogramming in cancer. Tumor cells can use circulating free fatty acids as an energy source through lipolysis, for membrane biosynthesis or in signaling processes [31].
There is also experimental evidence that lipogenesis in tumors, associated with increased expression of the fatty acid synthase gene (FASN), is strongly dependent on the activity and/or expression of important oncogenes and tumor suppressors, including MYC, which cooperates with the sterol regulatory element-binding proteins (SREBPs) and can induce in vitro and in vivo lipogenesis, thus playing an important role in initiating and maintaining oncogenic growth [32].

4. Materials and Methods

4.1. Materials

The PCa samples were obtained from Russian patients who had undergone surgical intervention in the P.A. Hertzen Moscow Oncology Research Center (a branch of the National Medical Research Radiological Center, Ministry of Health of the Russian Federation) between 2015 and 2020. All materials were collected and characterized by the organization’s Pathology Department according to the WHO Classification of Tumours of the Urinary System and Male Genital Organs [33]. Each sample contained a minimum of 70% of tumor cells. Following surgical resection, the tissue samples were immediately frozen and stored in liquid nitrogen. In the current study, we used 58 samples of LPCa (adenocarcinoma) samples obtained from patients who underwent surgical treatment but had not received neoadjuvant therapy (Table 3). Additionally, RNA-Seq data of 43 lymph node-negative LAPCa samples obtained in our previous study were included [34]. The samples have the following characteristics: no regional metastasis (N0 category); negative resection margins; any PSA value; and any Gleason score. Samples with the presence of regional metastasis (pN1) were not included in the study as they are characterized by a specific transcriptomic expression pattern.

4.2. Methods

4.2.1. Total RNA Isolation, Library Preparation, and Next Generation Sequencing (NGS)

Samples of frozen tumor tissues were first homogenized using a MagNA Lyser device (Roche, Basel, Switzerland). Subsequent total RNA isolation was performed using the MagNA Pure Compact RNA Kit (Roche) on the MagNA Pure Compact System (Roche) according to the manufacturer’s protocol. The concentration of isolated total RNA was assessed on a Qubit 4.0 fluorimeter (Thermo Fisher Scientific, Waltham, MA, USA) using the Qubit RNA BR Assay Kit (Thermo Fisher Scientific). The RIN (RNA integrity number) parameter, which characterizes the integrity of RNA, was evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The RIN for all samples studied was no less than 7.
Sample preparation of mRNA libraries was performed using the TruSeq Stranded mRNA Kit (Illumina, San Diego, CA, USA) as described previously [35]. The size of the resulting mRNA library was ~260 bp.
High-throughput sequencing of mRNA libraries was performed on a NextSeq 500 System (Illumina) using NextSeq 500/550 High Output Kit v2.5 (Illumina) in 75 bp single-ended read mode. On average, about 20 million reads were received for each sample.

4.2.2. Bioinformatic Analysis

FastQC (v.0.11.9, Cambridge, UK) (Babraham Bioinformatics-FastQC A Quality Control Tool for High throughput Sequence Data, accessed on 3 March 2021; available online, https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and Trimmomatic (v.0.33, Jülich, Germany) [36] were used for the quality control and trimming of reads, respectively. The STAR splice-aware aligner (v.2.7.1, Cold Spring Harbor, NY, USA) [37] was used to map the reads to the reference genome (GRCh37.p13, GENCODE, Cambridge, UK) (GENCODE—Human Release 19, accessed on 25 April 2021; available online, https://www.gencodegenes.org/human/release_19.html). FeatureCounts (Subread package v.1.6.4, Parkville, VIC, Australia) [38] was used to calculate the read counts per gene.
Differential expression analysis was carried out in the R environment (v.3.6.3, Vienna, Austria) (R: The R Project for Statistical Computing, accessed on 3 March 2023; available online, https://www.r-project.org/) using the edgeR package (v.3.24.3, NSW, Australia) [39], as described previously [34]. The results were considered significant at p-values of the quasi-likelihood F-test (QLF) and the Mann-Whitney U-test (MW) ≤ 0.05.
Gene set enrichment analysis (GSEA) was performed using the GSEApy package in Jupyter Notebook, Python (v.3.6) [40]. Annotation of the results was obtained, based on the KEGG Human 2021. The results were considered significant at FDR ≤ 0.05.
ROC analysis was performed on the basis of the Logistic Regression algorithm, using the scikit-learn library in Jupyter Notebook, Python (v.3.6).

4.2.3. Reverse Transcription and Quantitative PCR (qPCR)

cDNA samples were obtained from the mRNA template using Mint reverse transcriptase and oligo(dT) primer (20 µM) according to the manufacturer’s protocol (Evrogen, Moscow, Russia).
qPCR was performed in three technical replicates in total reaction volume of 10 μL on an Applied Biosystems 7500 instrument (Thermo Fisher Scientific). The TaqMan Gene Expression Assay Hs03063375_ft (Thermo Fisher Scientific) was used to determine the presence of the TMPRSS2-ERG fusion transcript. ROX was used as a reference dye. The PUM1 gene was used as a reference for analysis of relative mRNA expression. The sequences of primers used to validate markers based on mRNA expression are shown in Table 4.
The following process was used for amplification: 95 °C for 15 min; 40 cycles at 95 °C for 15 s; and 60 °C for 60 s. To assess the level of expression, the method of relative measurements (ΔCT) was used and calculations were performed using the ATG program (Analysis of Transcription of Genes) [41]. Visualization and statistical analysis of expression results were performed using the MW test in the R environment (v.3.6.3, Vienna, Austria).

5. Conclusions

We performed RNA-Seq profiling of 58 LPCa and 43 LAPCa tissue samples, obtained as a result of radical prostatectomy on a cohort of Russian patients. Bioinformatics analysis revealed that in the high-risk group, LAPCs showed enrichment of certain biological pathways, both across the entire sample and within the TMPRSS2-ERG molecular subtype. Such enrichment could be further investigated in the search for new potential therapeutic targets in the studied categories of PCa. We also determined which genes showed the greatest significant differences in expression, allowing the LPCa and LAPCa categories in the high-risk group to be distinguished, when taking into account the TMPRSS2-ERG molecular subtype. The highest predictive potential was found for the CIBAR1, CLIC4, EEF1A1P5, OGN, RPLP0P6, and ZNF483 genes.
The study also examined the transcriptomic features of the intermediate-risk group of LPCa, within GS = 7 (ISUP classification—groups 2 and 3). A statistically significant enrichment of the “PPAR signaling pathway” in the ISUP 3 group was shown. Based on the identified transcriptomic profile, the LPL, MYC, and TWIST1 genes were selected as promising additional prognostic markers, the statistical significance of which was confirmed based on qPCR validation.

Supplementary Materials

The supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms24119282/s1.

Author Contributions

Conceptualization, A.V.K. and B.Y.A.; data curation, K.M.N.; formal analysis, V.S.P., A.O.K., and I.V.K.; funding acquisition, G.S.K.; investigation, A.A.K. and E.A.P.; methodology, M.S.F. and E.A.P.; project administration, G.S.K., E.A.P., and A.A.K.; resources, G.S.K.; software, G.S.K.; supervision, A.V.S.; visualization, A.A.K. and E.A.P.; writing—original draft, A.A.K. and E.A.P.; writing—review and editing, E.A.P. and A.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant 18-75-10127.

Institutional Review Board Statement

Ethical review and approval were waived, as the current research was a retrospective study using medical records/biological specimens obtained in previous clinical practices.

Informed Consent Statement

Patient consent was waived due to the study being retrospective and its use of pseudonymized data.

Data Availability Statement

All data generated or analyzed during this study are available (GSE229904). The dataset also includes expression data from our previous studies devoted to the analysis of primary PCa tumors.

Acknowledgments

The authors thank the National Medical Research Center of Radiology for assistance in data interpretation and discussion, and EIMB RAS for providing the computing power and techniques for the data analysis. This work was performed using equipment at the EIMB RAS “Genome” center (http://www.eimb.ru/ru1/ckp/ccu_genome_c.php, accessed on 1 March 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pernar, C.H.; Ebot, E.M.; Wilson, K.M.; Mucci, L.A. The Epidemiology of Prostate Cancer. Cold Spring Harb. Perspect. Med. 2018, 8, a030361. [Google Scholar] [CrossRef] [PubMed]
  2. D’Amico, A.V.; Whittington, R.; Malkowicz, S.B.; Schultz, D.; Blank, K.; Broderick, G.A.; Tomaszewski, J.E.; Renshaw, A.A.; Kaplan, I.; Beard, C.J.; et al. Biochemical outcome after radical prostatectomy, external beam radiation therapy, or interstitial radiation therapy for clinically localized prostate cancer. JAMA 1998, 280, 969–974. [Google Scholar] [CrossRef] [PubMed]
  3. Bostrom, P.J.; Bjartell, A.S.; Catto, J.W.; Eggener, S.E.; Lilja, H.; Loeb, S.; Schalken, J.; Schlomm, T.; Cooperberg, M.R. Genomic Predictors of Outcome in Prostate Cancer. Eur. Urol. 2015, 68, 1033–1044. [Google Scholar] [CrossRef] [PubMed]
  4. Grasso, C.S.; Wu, Y.M.; Robinson, D.R.; Cao, X.; Dhanasekaran, S.M.; Khan, A.P.; Quist, M.J.; Jing, X.; Lonigro, R.J.; Brenner, J.C.; et al. The mutational landscape of lethal castration-resistant prostate cancer. Nature 2012, 487, 239–243. [Google Scholar] [CrossRef]
  5. Taylor, B.S.; Schultz, N.; Hieronymus, H.; Gopalan, A.; Xiao, Y.; Carver, B.S.; Arora, V.K.; Kaushik, P.; Cerami, E.; Reva, B.; et al. Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010, 18, 11–22. [Google Scholar] [CrossRef]
  6. Barbieri, C.E.; Bangma, C.H.; Bjartell, A.; Catto, J.W.; Culig, Z.; Gronberg, H.; Luo, J.; Visakorpi, T.; Rubin, M.A. The mutational landscape of prostate cancer. Eur. Urol. 2013, 64, 567–576. [Google Scholar] [CrossRef]
  7. Haffner, M.C.; Mosbruger, T.; Esopi, D.M.; Fedor, H.; Heaphy, C.M.; Walker, D.A.; Adejola, N.; Gurel, M.; Hicks, J.; Meeker, A.K.; et al. Tracking the clonal origin of lethal prostate cancer. J. Clin. Investig. 2013, 123, 4918–4922. [Google Scholar] [CrossRef]
  8. Kumar, A.; White, T.A.; MacKenzie, A.P.; Clegg, N.; Lee, C.; Dumpit, R.F.; Coleman, I.; Ng, S.B.; Salipante, S.J.; Rieder, M.J.; et al. Exome sequencing identifies a spectrum of mutation frequencies in advanced and lethal prostate cancers. Proc. Natl. Acad. Sci. USA 2011, 108, 17087–17092. [Google Scholar] [CrossRef]
  9. Cancer Genome Atlas Research Network. The Molecular Taxonomy of Primary Prostate Cancer. Cell 2015, 163, 1011–1025. [Google Scholar] [CrossRef]
  10. Arora, K.; Barbieri, C.E. Molecular Subtypes of Prostate Cancer. Curr. Oncol. Rep. 2018, 20, 58. [Google Scholar] [CrossRef]
  11. Kobelyatskaya, A.A.; Kudryavtsev, A.A.; Kudryavtseva, A.V.; Snezhkina, A.V.; Fedorova, M.S.; Kalinin, D.V.; Pavlov, V.S.; Guvatova, Z.G.; Naberezhnev, P.A.; Nyushko, K.M.; et al. ALDH3A2, ODF2, QSOX2, and MicroRNA-503-5p Expression to Forecast Recurrence in TMPRSS2-ERG-Positive Prostate Cancer. Int. J. Mol. Sci. 2022, 23, 11695. [Google Scholar] [CrossRef]
  12. Perner, S.; Mosquera, J.M.; Demichelis, F.; Hofer, M.D.; Paris, P.L.; Simko, J.; Collins, C.; Bismar, T.A.; Chinnaiyan, A.M.; De Marzo, A.M.; et al. TMPRSS2-ERG fusion prostate cancer: An early molecular event associated with invasion. Am. J. Surg. Pathol. 2007, 31, 882–888. [Google Scholar] [CrossRef]
  13. Adamo, P.; Ladomery, M.R. The oncogene ERG: A key factor in prostate cancer. Oncogene 2016, 35, 403–414. [Google Scholar] [CrossRef]
  14. Hagglof, C.; Hammarsten, P.; Stromvall, K.; Egevad, L.; Josefsson, A.; Stattin, P.; Granfors, T.; Bergh, A. TMPRSS2-ERG expression predicts prostate cancer survival and associates with stromal biomarkers. PLoS ONE 2014, 9, e86824. [Google Scholar] [CrossRef]
  15. Kim, S.H.; Kim, S.H.; Joung, J.Y.; Lee, G.K.; Hong, E.K.; Kang, K.M.; Yu, A.; Nam, B.H.; Chung, J.; Seo, H.K.; et al. Overexpression of ERG and Wild-Type PTEN Are Associated with Favorable Clinical Prognosis and Low Biochemical Recurrence in Prostate Cancer. PLoS ONE 2015, 10, e0122498. [Google Scholar] [CrossRef]
  16. Gerke, J.S.; Orth, M.F.; Tolkach, Y.; Romero-Perez, L.; Wehweck, F.S.; Stein, S.; Musa, J.; Knott, M.M.L.; Holting, T.L.B.; Li, J.; et al. Integrative clinical transcriptome analysis reveals TMPRSS2-ERG dependency of prognostic biomarkers in prostate adenocarcinoma. Int. J. Cancer 2020, 146, 2036–2046. [Google Scholar] [CrossRef]
  17. Zimta, A.A.; Sigurjonsson, O.E.; Gulei, D.; Tomuleasa, C. The Malignant Role of Exosomes as Nanocarriers of Rare RNA Species. Int. J. Mol. Sci. 2020, 21, 5866. [Google Scholar] [CrossRef]
  18. Carvalho, A.S.; Baeta, H.; Silva, B.C.; Moraes, M.C.S.; Bodo, C.; Beck, H.C.; Rodriguez, M.S.; Saraswat, M.; Pandey, A.; Matthiesen, R. Extra-cellular vesicles carry proteome of cancer hallmarks. Front. Biosci. 2020, 25, 398–436. [Google Scholar] [CrossRef]
  19. Sun, T.; Wang, X.; He, H.H.; Sweeney, C.J.; Liu, S.X.; Brown, M.; Balk, S.; Lee, G.S.; Kantoff, P.W. MiR-221 promotes the development of androgen independence in prostate cancer cells via downregulation of HECTD2 and RAB1A. Oncogene 2014, 33, 2790–2800. [Google Scholar] [CrossRef]
  20. Mirzaei, S.; Paskeh, M.D.A.; Saghari, Y.; Zarrabi, A.; Hamblin, M.R.; Entezari, M.; Hashemi, M.; Aref, A.R.; Hushmandi, K.; Kumar, A.P.; et al. Transforming growth factor-beta (TGF-beta) in prostate cancer: A dual function mediator? Int. J. Biol. Macromol. 2022, 206, 435–452. [Google Scholar] [CrossRef]
  21. Pudova, E.A.; Krasnov, G.S.; Nyushko, K.M.; Kobelyatskaya, A.A.; Savvateeva, M.V.; Poloznikov, A.A.; Dolotkazin, D.R.; Klimina, K.M.; Guvatova, Z.G.; Simanovsky, S.A.; et al. miRNAs expression signature potentially associated with lymphatic dissemination in locally advanced prostate cancer. BMC Med. Genom. 2020, 13, 129. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, T.; Wang, W.; Shi, G.; Hao, M.; Wang, Y.; Yao, M.; Huang, Y.; Du, L.; Zhang, X.; Ye, D.; et al. Targeting HIC1/TGF-beta axis-shaped prostate cancer microenvironment restrains its progression. Cell. Death Dis. 2022, 13, 624. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, H.; Kong, Q.; Wang, J.; Jiang, Y.; Hua, H. Complex roles of cAMP-PKA-CREB signaling in cancer. Exp. Hematol. Oncol. 2020, 9, 32. [Google Scholar] [CrossRef]
  24. Merkle, D.; Hoffmann, R. Roles of cAMP and cAMP-dependent protein kinase in the progression of prostate cancer: Cross-talk with the androgen receptor. Cell. Signal 2011, 23, 507–515. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, Z.; Rahman, M.A.; Chen, Z.G.; Shin, D.M. Multiple biological functions of Twist1 in various cancers. Oncotarget 2017, 8, 20380–20393. [Google Scholar] [CrossRef]
  26. Gajula, R.P.; Chettiar, S.T.; Williams, R.D.; Thiyagarajan, S.; Kato, Y.; Aziz, K.; Wang, R.; Gandhi, N.; Wild, A.T.; Vesuna, F.; et al. The twist box domain is required for Twist1-induced prostate cancer metastasis. Mol. Cancer Res. 2013, 11, 1387–1400. [Google Scholar] [CrossRef]
  27. Madden, S.K.; de Araujo, A.D.; Gerhardt, M.; Fairlie, D.P.; Mason, J.M. Taking the Myc out of cancer: Toward therapeutic strategies to directly inhibit c-Myc. Mol. Cancer 2021, 20, 3. [Google Scholar] [CrossRef]
  28. Rebello, R.J.; Pearson, R.B.; Hannan, R.D.; Furic, L. Therapeutic Approaches Targeting MYC-Driven Prostate Cancer. Genes 2017, 8, 71. [Google Scholar] [CrossRef]
  29. Pingitore, P.; Lepore, S.M.; Pirazzi, C.; Mancina, R.M.; Motta, B.M.; Valenti, L.; Berge, K.E.; Retterstol, K.; Leren, T.P.; Wiklund, O.; et al. Identification and characterization of two novel mutations in the LPL gene causing type I hyperlipoproteinemia. J. Clin. Lipidol. 2016, 10, 816–823. [Google Scholar] [CrossRef]
  30. Kuemmerle, N.B.; Rysman, E.; Lombardo, P.S.; Flanagan, A.J.; Lipe, B.C.; Wells, W.A.; Pettus, J.R.; Froehlich, H.M.; Memoli, V.A.; Morganelli, P.M.; et al. Lipoprotein lipase links dietary fat to solid tumor cell proliferation. Mol. Cancer Ther. 2011, 10, 427–436. [Google Scholar] [CrossRef]
  31. Wang, W.; Bai, L.; Li, W.; Cui, J. The Lipid Metabolic Landscape of Cancers and New Therapeutic Perspectives. Front. Oncol. 2020, 10, 605154. [Google Scholar] [CrossRef]
  32. Gouw, A.M.; Margulis, K.; Liu, N.S.; Raman, S.J.; Mancuso, A.; Toal, G.G.; Tong, L.; Mosley, A.; Hsieh, A.L.; Sullivan, D.K.; et al. The MYC Oncogene Cooperates with Sterol-Regulated Element-Binding Protein to Regulate Lipogenesis Essential for Neoplastic Growth. Cell. Metab. 2019, 30, 556–572.e55. [Google Scholar] [CrossRef]
  33. Moch, H.; Cubilla, A.L.; Humphrey, P.A.; Reuter, V.E.; Ulbright, T.M. The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs-Part A: Renal, Penile, and Testicular Tumours. Eur. Urol. 2016, 70, 93–105. [Google Scholar] [CrossRef]
  34. Kobelyatskaya, A.A.; Pudova, E.A.; Snezhkina, A.V.; Fedorova, M.S.; Pavlov, V.S.; Guvatova, Z.G.; Savvateeva, M.V.; Melnikova, N.V.; Dmitriev, A.A.; Trofimov, D.Y.; et al. Impact TMPRSS2-ERG Molecular Subtype on Prostate Cancer Recurrence. Life 2021, 11, 588. [Google Scholar] [CrossRef]
  35. Pudova, E.A.; Kobelyatskaya, A.A.; Katunina, I.V.; Snezhkina, A.V.; Fedorova, M.S.; Pavlov, V.S.; Bakhtogarimov, I.R.; Lantsova, M.S.; Kokin, S.P.; Nyushko, K.M.; et al. Lymphatic Dissemination in Prostate Cancer: Features of the Transcriptomic Profile and Prognostic Models. Int. J. Mol. Sci. 2023, 24, 2418. [Google Scholar] [CrossRef]
  36. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  37. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  38. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef]
  39. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
  40. Fang, Z.; Liu, X.; Peltz, G. GSEApy: A comprehensive package for performing gene set enrichment analysis in Python. Bioinformatics 2023, 39, btac757. [Google Scholar] [CrossRef]
  41. Senchenko, V.N.; Krasnov, G.S.; Dmitriev, A.A.; Kudryavtseva, A.V.; Anedchenko, E.A.; Braga, E.A.; Pronina, I.V.; Kondratieva, T.T.; Ivanov, S.V.; Zabarovsky, E.R.; et al. Differential expression of CHL1 gene during development of major human cancers. PLoS ONE 2011, 6, e15612. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Heatmap of top 50 differentially expressed genes between LAPCa and LPCa. LAPCa—locally advanced PCa; LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Figure 1. Heatmap of top 50 differentially expressed genes between LAPCa and LPCa. LAPCa—locally advanced PCa; LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Ijms 24 09282 g001
Figure 2. Dotplot of GSEA results illustrating the KEGG-based biological processes associated with the LAPCa category versus LPCa within the high-risk group. The figure shows statistically significant enriched biological pathways based on coexpressed genes. NES—normalized enrichment score.
Figure 2. Dotplot of GSEA results illustrating the KEGG-based biological processes associated with the LAPCa category versus LPCa within the high-risk group. The figure shows statistically significant enriched biological pathways based on coexpressed genes. NES—normalized enrichment score.
Ijms 24 09282 g002
Figure 3. Heatmap of top 50 differentially expressed genes between LAPCa and LPCa for the TMPRSS2-ERG subtype. LAPCa—locally advanced PCa; LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Figure 3. Heatmap of top 50 differentially expressed genes between LAPCa and LPCa for the TMPRSS2-ERG subtype. LAPCa—locally advanced PCa; LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Ijms 24 09282 g003
Figure 4. Dotplot of GSEA results illustrating KEGG-based biological processes associated with the LAPCa category versus LPCa for the high-risk group within the TMPRSS2-ERG molecular subtype. The figure shows statistically significantly enriched biological pathways based on the coexpression of genes. NES—normalized enrichment score.
Figure 4. Dotplot of GSEA results illustrating KEGG-based biological processes associated with the LAPCa category versus LPCa for the high-risk group within the TMPRSS2-ERG molecular subtype. The figure shows statistically significantly enriched biological pathways based on the coexpression of genes. NES—normalized enrichment score.
Ijms 24 09282 g004
Figure 5. Heatmap of identified differentially expressed genes between ISUP 3 and ISUP 2 for LPCa. LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Figure 5. Heatmap of identified differentially expressed genes between ISUP 3 and ISUP 2 for LPCa. LPCa—localized PCa. Red color indicates upregulated genes, blue—downregulated.
Ijms 24 09282 g005
Figure 6. GSEA enrichment plot in ISUP 3 group versus ISUP 2 group in LPCa intermediate-risk group.
Figure 6. GSEA enrichment plot in ISUP 3 group versus ISUP 2 group in LPCa intermediate-risk group.
Ijms 24 09282 g006
Figure 7. Boxplots of the relative expression of mRNA LPL, MYC, and TWIST1 genes between the ISUP 3 and ISUP 2 groups in the LPCa intermediate-risk group of the Russian patient cohort. ISUP 2 group is marked in red, ISUP 3 group—blue.
Figure 7. Boxplots of the relative expression of mRNA LPL, MYC, and TWIST1 genes between the ISUP 3 and ISUP 2 groups in the LPCa intermediate-risk group of the Russian patient cohort. ISUP 2 group is marked in red, ISUP 3 group—blue.
Ijms 24 09282 g007
Table 1. Differential expression of genes statistically significantly associated with the LAPCa group, including within the TMPRSS2-ERG subtype. LAPCa—locally advanced prostate cancer, LPCa—localized prostate cancer, Log2FC—log2 transformed fold change, rs—Spearman correlation coefficient, grey background—downregulated genes.
Table 1. Differential expression of genes statistically significantly associated with the LAPCa group, including within the TMPRSS2-ERG subtype. LAPCa—locally advanced prostate cancer, LPCa—localized prostate cancer, Log2FC—log2 transformed fold change, rs—Spearman correlation coefficient, grey background—downregulated genes.
Gene Name
(Biotype)
DescriptionLAPCa vs. LPCaLAPCa vs. LPCa (TMPRSS2-ERG)
Log2FCrsLog2FCrs
BHLHA15
(protein coding)
basic helix-loop-helix family member a151.790.571.400.68
CIBAR1
(protein coding)
CBY1 interacting BAR domain containing 1−1.05−0.65−1.20−0.81
CLIC4
(protein coding)
chloride intracellular channel 4−1.02−0.53−1.14−0.81
CORO1B
(protein coding)
coronin 1B1.170.631.040.77
CRB3
(protein coding)
crumbs cell polarity complex component 31.150.581.030.72
DNAJB4
(protein coding)
DnaJ heat shock protein family (Hsp40) member B4−1.13−0.62−1.22−0.80
DNM3OS
(lncRNA)
DNM3 opposite strand/antisense RNA−1.22−0.62−1.58−0.75
EEF1A1P5
(processed pseudogene)
eukaryotic translation elongation factor 1 alpha 1 pseudogene 52.450.772.700.87
HECTD2
(protein coding)
HECT domain E3 ubiquitin protein ligase 2−1.05−0.63−1.04−0.77
ID4
(protein coding)
inhibitor of DNA binding 4, HLH protein−1.34−0.555.85−0.67
MFSD3
(protein coding)
major facilitator superfamily domain containing 31.100.611.090.73
MIR222HG
(lncRNA)
miR222/221 cluster host gene−1.27−0.56−1.27−0.60
OGN
(protein coding)
osteoglycin−1.10−0.54−1.62−0.76
RPLP0P6
(processed pseudogene)
ribosomal protein lateral stalk subunit P0 pseudogene 61.830.701.940.77
SH3BGRL
(protein coding)
SH3 domain binding glutamate rich protein like−1.06−0.63−1.12−0.76
ZNF483
(protein coding)
zinc finger protein 483−1.02−0.58−1.15−0.78
Table 2. Results of ROC analysis based on the logistic regression algorithm for test data based on a cohort of Russian patients with high-risk LPCa and LAPCa, considering the TMPRSS2-ERG molecular subtype. LAPCa—locally advanced prostate cancer, LPCa—localized prostate cancer, AUC—Area Under Curve, CI—confidence interval, grey background—downregulated genes.
Table 2. Results of ROC analysis based on the logistic regression algorithm for test data based on a cohort of Russian patients with high-risk LPCa and LAPCa, considering the TMPRSS2-ERG molecular subtype. LAPCa—locally advanced prostate cancer, LPCa—localized prostate cancer, AUC—Area Under Curve, CI—confidence interval, grey background—downregulated genes.
Gene NameLAPCa vs. LPCaLAPCa vs. LPCa
(TMPRSS2-ERG)
AUC (CI, 95%)Accuracyp-ValueAUC (CI, 95%)Accuracyp-Value
BHLHA150.864
(0.581–0.964)
0.7506.0 × 10−60.867
(0.820–1.000)
0.7081.2 × 10−2
CIBAR10.821
(0.701–1.000)
0.7083.0 × 10−60.979
(0.877–1.000)
0.9178.8 × 10−3
CLIC40.711
(0.611–0.985)
0.7508.1 × 10−50.930
(0.891–1.000)
0.8338.5 × 10−3
CORO1B0.830
(0.701–1.000)
0.7501.0 × 10−60.889
(0.897–1.000)
0.7922.6 × 10−2
CRB30.798
(0.661–1.000)
0.7502.0 × 10−60.734
(0.686–1.000)
0.5426.0 × 10−3
DNAJB40.867
(0.654–0.992)
0.7921.0 × 10−50.824
(0.790–1.000)
0.7088.4 × 10−3
DNM3OS0.867
(0.718–1.000)
0.7921.2 × 10−50.824
(0.807–1.000)
0.7089.7 × 10−3
EEF1A1P50.979
(0.904–1.000)
0.9176.5 × 10−51.000
(0.977–1.000)
0.8752.9 × 10−2
HECTD20.815
(0.691–1.000)
0.7086.0 × 10−60.977
(0.947–1.000)
0.8755.7 × 10−3
ID40.832
(0.486–0.924)
0.7503.9 × 10−50.811
(0.592–1.000)
0.7508.3 × 10−3
MFSD30.837
(0.604–0.992)
0.6672.0 × 10−50.874
(0.754–1.000)
0.7088.5 × 10−3
MIR222HG0.778
(0.512–0.950)
0.7081.3 × 10−40.879
(0.763–1.000)
0.7504.1 × 10−2
OGN0.764
(0.541–0.988)
0.7083.9 × 10−50.963
(0.797–1.000)
0.8751.6 × 10−2
RPLP0P60.930
(0.745–1.000)
0.7925.2 × 10−60.984
(0.764–1.000)
0.9582.3 × 10−2
SH3BGRL0.857
(0.581–0.981)
0.7926.0 × 10−60.818
(0.757–1.000)
0.7929.2 × 10−3
ZNF4830.879
(0.501–0.944)
0.7501.5 × 10−50.914
(0.891–1.000)
0.7087.3 × 10−3
Table 3. Clinical and pathological characteristics of the studied cohort. pT—primary tumor estimation; N—regional lymph nodes; M—distant metastases; ISUP—The International Society of Urological Pathology.
Table 3. Clinical and pathological characteristics of the studied cohort. pT—primary tumor estimation; N—regional lymph nodes; M—distant metastases; ISUP—The International Society of Urological Pathology.
CriterionLPCa, nLAPCa, n
PCa samples5843
Age, years64 (41–78)64 (46–77)
pTpT2a30
pT2b50
pT2c500
pT3a030
pT3b013
pT400
pNpN05843
pN100
cMcM05843
cM100
Gleason score6226
73229
825
913
1000
ISUP grade1226
22017
31212
425
513
PSA, ng/mL11.1 (0.3–30)12.8 (3.4–27.6)
Risk groupLow10
Intermediate470
High1043
Molecular subtype TMPRSS2-ERGYes1618
No4225
Table 4. Primer sequences for assessing the level of mRNA expression.
Table 4. Primer sequences for assessing the level of mRNA expression.
mRNAPrimer Sequence (5′→3′)Product Length, b.p.
LPLF: CAGCCCTACCCTTGTTAGTTATT
R: ACGTTGGAGGATGTGCTATTT
103
MYCF: ATCTCTGGGAGGAATGCTACTA
R: ATCTGCGTGGCTACAGATAAG
95
TWISTF: CGGAGACCTAGATGTCATTGTTT
R: ACGCCCTGTTTCTTTGAATTTG
146
PUM1F: TGGACCATTTCGCCCTTTAG
R: CAGAGAGTTGTTGCCGTAGAA
103
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kobelyatskaya, A.A.; Pudova, E.A.; Katunina, I.V.; Snezhkina, A.V.; Fedorova, M.S.; Pavlov, V.S.; Kotelnikova, A.O.; Nyushko, K.M.; Alekseev, B.Y.; Krasnov, G.S.; et al. Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype. Int. J. Mol. Sci. 2023, 24, 9282. https://doi.org/10.3390/ijms24119282

AMA Style

Kobelyatskaya AA, Pudova EA, Katunina IV, Snezhkina AV, Fedorova MS, Pavlov VS, Kotelnikova AO, Nyushko KM, Alekseev BY, Krasnov GS, et al. Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype. International Journal of Molecular Sciences. 2023; 24(11):9282. https://doi.org/10.3390/ijms24119282

Chicago/Turabian Style

Kobelyatskaya, Anastasiya A., Elena A. Pudova, Irina V. Katunina, Anastasiya V. Snezhkina, Maria S. Fedorova, Vladislav S. Pavlov, Anastasiya O. Kotelnikova, Kirill M. Nyushko, Boris Y. Alekseev, George S. Krasnov, and et al. 2023. "Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype" International Journal of Molecular Sciences 24, no. 11: 9282. https://doi.org/10.3390/ijms24119282

APA Style

Kobelyatskaya, A. A., Pudova, E. A., Katunina, I. V., Snezhkina, A. V., Fedorova, M. S., Pavlov, V. S., Kotelnikova, A. O., Nyushko, K. M., Alekseev, B. Y., Krasnov, G. S., & Kudryavtseva, A. V. (2023). Transcriptome Profiling of Prostate Cancer, Considering Risk Groups and the TMPRSS2-ERG Molecular Subtype. International Journal of Molecular Sciences, 24(11), 9282. https://doi.org/10.3390/ijms24119282

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop