1
|
Zolotovskaia M, Kovalenko M, Pugacheva P, Tkachev V, Simonov A, Sorokin M, Seryakov A, Garazha A, Gaifullin N, Sekacheva M, Zakharova G, Buzdin AA. Algorithmically Reconstructed Molecular Pathways as the New Generation of Prognostic Molecular Biomarkers in Human Solid Cancers. Proteomes 2023; 11:26. [PMID: 37755705 PMCID: PMC10535530 DOI: 10.3390/proteomes11030026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Revised: 08/18/2023] [Accepted: 08/22/2023] [Indexed: 09/28/2023] Open
Abstract
Individual gene expression and molecular pathway activation profiles were shown to be effective biomarkers in many cancers. Here, we used the human interactome model to algorithmically build 7470 molecular pathways centered around individual gene products. We assessed their associations with tumor type and survival in comparison with the previous generation of molecular pathway biomarkers (3022 "classical" pathways) and with the RNA transcripts or proteomic profiles of individual genes, for 8141 and 1117 samples, respectively. For all analytes in RNA and proteomic data, respectively, we found a total of 7441 and 7343 potential biomarker associations for gene-centric pathways, 3020 and 2950 for classical pathways, and 24,349 and 6742 for individual genes. Overall, the percentage of RNA biomarkers was statistically significantly higher for both types of pathways than for individual genes (p < 0.05). In turn, both types of pathways showed comparable performance. The percentage of cancer-type-specific biomarkers was comparable between proteomic and transcriptomic levels, but the proportion of survival biomarkers was dramatically lower for proteomic data. Thus, we conclude that pathway activation level is the advanced type of biomarker for RNA and proteomic data, and momentary algorithmic computer building of pathways is a new credible alternative to time-consuming hypothesis-driven manual pathway curation and reconstruction.
Collapse
Affiliation(s)
- Marianna Zolotovskaia
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
- Omicsway Corp., Walnut, CA 91789, USA
- Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Maks Kovalenko
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
| | - Polina Pugacheva
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
| | | | - Alexander Simonov
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
- Omicsway Corp., Walnut, CA 91789, USA
| | - Maxim Sorokin
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
- Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), 1200 Brussels, Belgium
| | | | | | - Nurshat Gaifullin
- Department of Pathology, Faculty of Medicine, Lomonosov Moscow State University, 119991 Moscow, Russia
| | - Marina Sekacheva
- Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Galina Zakharova
- Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Anton A. Buzdin
- Laboratory for Translational Genomic Bioinformatics, Moscow Institute of Physics and Technology (State University), 141701 Dolgoprudny, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), 1200 Brussels, Belgium
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
- Laboratory of Systems Biology, Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia
| |
Collapse
|
2
|
Kantelhardt S. New Strategies in Diagnosis and Treatments for Brain Tumors. Cancers (Basel) 2023; 15:cancers15112879. [PMID: 37296841 DOI: 10.3390/cancers15112879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 05/17/2023] [Indexed: 06/12/2023] Open
Abstract
In general, cancer is one of the most frequent causes of death [...].
Collapse
Affiliation(s)
- Sven Kantelhardt
- Department of Neurosurgery, Vivantes Hospital im Friedrichshain, Landsberger Allee 49, 10249 Berlin, Germany
| |
Collapse
|
3
|
Tabbal M, Hachim MY, Jan RK, Adrian TE. Using publicly available datasets to identify population-based transcriptomic landscape contributing to the aggressiveness of breast cancer in young women. Front Genet 2023; 13:1039037. [PMID: 36685821 PMCID: PMC9845274 DOI: 10.3389/fgene.2022.1039037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2022] [Accepted: 11/30/2022] [Indexed: 01/05/2023] Open
Abstract
Introduction: Although the risk of breast cancer increases with advancing age, some regions have larger number of young breast cancer patients (≤45 years-old), such as the Middle East, Eastern Asia, and North Africa, with more aggressive and poorly differentiated tumors. We aimed to conduct an in-silico analysis in an attempt to understand the aggressive nature of early-onset breast cancer, and to identify potential drivers of early-onset breast cancer using gene expression profiling datasets in a population-dependent manner. Methods: Functional genomics experiments data were acquired from cBioPortal database for cancer genomics, followed by the stratification of patients based on the age at representation of breast cancer and race. Differential gene expression analysis and gene amplification status analysis were carried out, followed by hub gene, transcription factor, and signalling pathway identification. Results: PAM50 subtype analysis revealed that young patients (≤45 years-old) had four-fold more basal tumors and worst progression-free survival (median of 101 months), compared with the 45-65 years group (median of 168 months). Fourteen genes were amplified in more than 14% of patients with an early-onset breast cancer. Interestingly, FREM2, LINC00332, and LINC00366 were exclusively amplified in younger patients. Gene expression data from three different populations (Asian, White, and African) revealed a unique transcriptomic profile of young patients, which was also reflected on the PAM50 subtype analysis. Our data indicates a higher tendency of young African patients to develop basal tumors, while young Asian patients are more prone to developing Luminal A tumors. Most genes that were found to be upregulated in younger patients are involved in important signaling pathways that promote cancer progression and metastasis, such as MAPK pathway, Reelin pathway and the PI3K/Akt pathway. Conclusion: This study provides strong evidence that the molecular profile of tumors derived from young breast cancer patients of different populations is unique and may explain the aggressiveness of these tumors, stressing the need to conduct population- based multi-omic analyses to identify the potential drivers for tumorigenesis and molecular profiles of young breast cancer patients.
Collapse
|
4
|
Zakharova G, Efimov V, Raevskiy M, Rumiantsev P, Gudkov A, Belogurova-Ovchinnikova O, Sorokin M, Buzdin A. Reclassification of TCGA Diffuse Glioma Profiles Linked to Transcriptomic, Epigenetic, Genomic and Clinical Data, According to the 2021 WHO CNS Tumor Classification. Int J Mol Sci 2022; 24:ijms24010157. [PMID: 36613601 PMCID: PMC9820617 DOI: 10.3390/ijms24010157] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2022] [Revised: 11/25/2022] [Accepted: 12/20/2022] [Indexed: 12/24/2022] Open
Abstract
In 2021, the fifth edition of the WHO classification of tumors of the central nervous system (WHO CNS5) was published. Molecular features of tumors were directly incorporated into the diagnostic decision tree, thus affecting both the typing and staging of the tumor. It has changed the traditional approach, based solely on histopathological classification. The Cancer Genome Atlas project (TCGA) is one of the main sources of molecular information about gliomas, including clinically annotated transcriptomic and genomic profiles. Although TCGA itself has played a pivotal role in developing the WHO CNS5 classification, its proprietary databases still retain outdated diagnoses which frequently appear incorrect and misleading according to the WHO CNS5 standards. We aimed to define the up-to-date annotations for gliomas from TCGA's database that other scientists can use in their research. Based on WHO CNS5 guidelines, we developed an algorithm for the reclassification of TCGA glioma samples by molecular features. We updated tumor type and diagnosis for 828 out of a total of 1122 TCGA glioma cases, after which available transcriptomic and methylation data showed clustering features more consistent with the updated grouping. We also observed better stratification by overall survival for the updated diagnoses, yet WHO grade 3 IDH-mutant oligodendrogliomas and astrocytomas are still indistinguishable. We also detected altered performance in the previous diagnostic transcriptomic molecular biomarkers (expression of SPRY1, CRNDE and FREM2 genes and FREM2 molecular pathway) and prognostic gene signature (FN1, ITGA5, OSMR, and NGFR) after reclassification. Thus, we conclude that further efforts are needed to reconsider glioma molecular biomarkers.
Collapse
Affiliation(s)
- Galina Zakharova
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Victor Efimov
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Mikhail Raevskiy
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | - Pavel Rumiantsev
- Multidisciplinary Medical Center, Group of Clinics, 194044 Saint-Petersburg, Russia
| | - Alexander Gudkov
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
| | | | - Maksim Sorokin
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), 1200 Brussels, Belgium
| | - Anton Buzdin
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119048 Moscow, Russia
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), 1200 Brussels, Belgium
- Correspondence:
| |
Collapse
|
5
|
Molecular markers related to patient outcome in patients with IDH-mutant astrocytomas grade 2 to 4: A systematic review. Eur J Cancer 2022; 175:214-223. [PMID: 36152406 DOI: 10.1016/j.ejca.2022.08.016] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 08/09/2022] [Accepted: 08/10/2022] [Indexed: 11/21/2022]
Abstract
BACKGROUND Grading and classification of IDH-mutant astrocytomas has shifted from solely histology towards histology combined with molecular diagnostics. In this systematic review, we give an overview of all currently known clinically relevant molecular markers within IDH-mutant astrocytomas grade 2 to 4. METHODS A literature search was performed in five electronic databases for English original papers on patient outcome with respect to a molecular marker as determined by DNA/RNA sequencing, micro-arrays, or DNA methylation profiling in IDH-mutant astrocytomas grade 2 to 4. Papers were included if molecular diagnostics were performed on tumour tissue of at least 15 IDH-mutant astrocytoma patients, and if the investigated molecular markers were not limited to the diagnostic markers MGMT, ATRX, TERT, and/or TP53. RESULTS The literature search identified 4508 unique articles, published between August 2012 and December 2021, of which ultimately 44 articles were included. Numerous molecular markers from these papers were significantly correlated to patient outcome. The associations between patient outcome and non-canonical IDH mutations, PI3K mutations, high expression of MSH2, high expression of RAD18, homozygous deletion of CDKN2A/B, amplification of PDGFRA, copy number neutral loss of chromosomal arm 17p, loss of chromosomal arm 19q, the G-CIMP-low DNA methylation cluster, high total CNV, and high tumour mutation burden were confirmed in multiple studies. CONCLUSIONS Multiple genetic and epigenetic markers are associated with survival in IDH-mutant astrocytoma patients. Commonly affected are the RB signalling pathway, the RTK-PI3K-mTOR signalling pathway, genomic stability markers, and (epigenetic) gene regulation.
Collapse
|
6
|
Nayak C, Singh SK. Integrated Transcriptome Profiling Identifies Prognostic Hub Genes as Therapeutic Targets of Glioblastoma: Evidenced by Bioinformatics Analysis. ACS OMEGA 2022; 7:22531-22550. [PMID: 35811900 PMCID: PMC9260928 DOI: 10.1021/acsomega.2c01820] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Accepted: 06/01/2022] [Indexed: 06/15/2023]
Abstract
Glioblastoma (GBM) is the most devastating and frequent type of primary brain tumor with high morbidity and mortality. Despite the use of surgical resection followed by radio- and chemotherapy as standard therapy, the progression of GBM remains dismal with a median overall survival of <15 months. GBM embodies a populace of cancer stem cells (GSCs) that is associated with tumor initiation, invasion, therapeutic resistance, and post-treatment reoccurrence. However, understanding the potential mechanisms of stemness and their candidate biomarkers remains limited. Hence in this investigation, we aimed to illuminate potential candidate hub genes and key pathways associated with the pathogenesis of GSC in the development of GBM. The integrated analysis discovered differentially expressed genes (DEGs) between the brain cancer tissues (GBM and GSC) and normal brain tissues. Multiple approaches, including gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were employed to functionally annotate the DEGs and visualize them through the R program. The significant hub genes were identified through the protein-protein interaction network, Venn diagram analysis, and survival analysis. We observed that the upregulated DEGs were prominently involved in the ECM-receptor interaction pathway. The downregulated genes were mainly associated with the axon guidance pathway. Five significant hub genes (CTNNB1, ITGB1, TNC, EGFR, and SHOX2) were screened out through multiple analyses. GO and KEGG analyses of hub genes uncovered that these genes were primarily enriched in disease-associated pathways such as the inhibition of apoptosis and the DNA damage repair mechanism, activation of the cell cycle, EMT (epithelial-mesenchymal transition), hormone AR (androgen receptor), hormone ER (estrogen receptor), PI3K/AKT (phosphatidylinositol 3-kinase and AKT), RTK (receptor tyrosine kinase), and TSC/mTOR (tuberous sclerosis complex and mammalian target of rapamycin). Consequently, the epigenetic regulatory network disclosed that hub genes played a vital role in the progression of GBM. Finally, candidate drugs were predicted that can be used as possible drugs to treat GBM patients. Overall, our investigation offered five hub genes (CTNNB1, ITGB1, TNC, EGFR, and SHOX2) that could be used as precise diagnostic and prognostic candidate biomarkers of GBM and might be used as personalized therapeutic targets to obstruct gliomagenesis.
Collapse
|
7
|
Zolotovskaia MA, Kovalenko MA, Tkachev VS, Simonov AM, Sorokin MI, Kim E, Kuzmin DV, Karademir-Yilmaz B, Buzdin AA. Next-Generation Grade and Survival Expression Biomarkers of Human Gliomas Based on Algorithmically Reconstructed Molecular Pathways. Int J Mol Sci 2022; 23:ijms23137330. [PMID: 35806337 PMCID: PMC9266372 DOI: 10.3390/ijms23137330] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Revised: 06/24/2022] [Accepted: 06/25/2022] [Indexed: 02/04/2023] Open
Abstract
In gliomas, expression of certain marker genes is strongly associated with survival and tumor type and often exceeds histological assessments. Using a human interactome model, we algorithmically reconstructed 7494 new-type molecular pathways that are centered each on an individual protein. Each single-gene expression and gene-centric pathway activation was tested as a survival and tumor grade biomarker in gliomas and their diagnostic subgroups (IDH mutant or wild type, IDH mutant with 1p/19q co-deletion, MGMT promoter methylated or unmethylated), including the three major molecular subtypes of glioblastoma (proneural, mesenchymal, classical). We used three datasets from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas, which in total include 527 glioblastoma and 1097 low grade glioma profiles. We identified 2724 such gene and 2418 pathway survival biomarkers out of total 17,717 genes and 7494 pathways analyzed. We then assessed tumor grade and molecular subtype biomarkers and with the threshold of AUC > 0.7 identified 1322/982 gene biomarkers and 472/537 pathway biomarkers. This suggests roughly two times greater efficacy of the reconstructed pathway approach compared to gene biomarkers. Thus, we conclude that activation levels of algorithmically reconstructed gene-centric pathways are a potent class of new-generation diagnostic and prognostic biomarkers for gliomas.
Collapse
Affiliation(s)
- Marianna A. Zolotovskaia
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia; (M.A.K.); (A.M.S.); (M.I.S.); (D.V.K.)
- Correspondence: ; Tel.: +7-9165612175
| | - Max A. Kovalenko
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia; (M.A.K.); (A.M.S.); (M.I.S.); (D.V.K.)
| | | | - Alexander M. Simonov
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia; (M.A.K.); (A.M.S.); (M.I.S.); (D.V.K.)
- Omicsway Corp., Walnut, CA 91789, USA;
| | - Maxim I. Sorokin
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia; (M.A.K.); (A.M.S.); (M.I.S.); (D.V.K.)
- Omicsway Corp., Walnut, CA 91789, USA;
- Laboratory of Clinical and Genomic Bioinformatics, I.M. Sechenov First Moscow State Medical University, 119991 Moscow, Russia
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia;
| | - Ella Kim
- Clinic for Neurosurgery, Laboratory of Experimental Neurooncology, Johannes Gutenberg University Medical Centre, Langenbeckstrasse 1, 55124 Mainz, Germany;
| | - Denis V. Kuzmin
- Moscow Institute of Physics and Technology, 141701 Dolgoprudny, Russia; (M.A.K.); (A.M.S.); (M.I.S.); (D.V.K.)
| | - Betul Karademir-Yilmaz
- Department of Biochemistry, School of Medicine/Genetic and Metabolic Diseases Research and Investigation Center (GEMHAM), Marmara University, Istanbul 34854, Turkey;
| | - Anton A. Buzdin
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, 117997 Moscow, Russia;
- World-Class Research Center “Digital Biodesign and Personalized Healthcare”, Sechenov First Moscow State Medical University, 119991 Moscow, Russia
- PathoBiology Group, European Organization for Research and Treatment of Cancer (EORTC), 1200 Brussels, Belgium
| |
Collapse
|
8
|
Zhang C, Zhang Y, Tan G, Mi W, Zhong X, Zhang Y, Zhao Z, Li F, Xu Y, Zhang Y. Prognostic Features of the Tumor Immune Microenvironment in Glioma and Their Clinical Applications: Analysis of Multiple Cohorts. Front Immunol 2022; 13:853074. [PMID: 35677045 PMCID: PMC9168240 DOI: 10.3389/fimmu.2022.853074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 04/25/2022] [Indexed: 11/13/2022] Open
Abstract
Glioma is the most common malignant tumor of the central nervous system. Tumor purity is a source of important prognostic factor for glioma patients, showing the key roles of the microenvironment in glioma prognosis. In this study, we systematically screened functional characterization related to the tumor immune microenvironment and constructed a risk model named Glioma MicroEnvironment Functional Signature (GMEFS) based on eight cohorts. The prognostic value of the GMEFS model was also verified in another two glioma cohorts, glioblastoma (GBM) and low-grade glioma (LGG) cohorts, from The Cancer Genome Atlas (TCGA). Nomograms were established in the training and testing cohorts to validate the clinical use of this model. Furthermore, the relationships between the risk score, intrinsic molecular subtypes, tumor purity, and tumor-infiltrating immune cell abundance were also evaluated. Meanwhile, the performance of the GMEFS model in glioma formation and glioma recurrence was systematically analyzed based on 16 glioma cohorts from the Gene Expression Omnibus (GEO) database. Based on multiple-cohort integrated analysis, risk subpathway signatures were identified, and a drug–subpathway association network was further constructed to explore candidate therapy target regions. Three subpathways derived from Focal adhesion (path: 04510) were identified and contained known targets including platelet derived growth factor receptor alpha (PDGFRA), epidermal growth factor receptor (EGFR), and erb-b2 receptor tyrosine kinase 2 (ERBB2). In conclusion, the novel functional signatures identified in this study could serve as a robust prognostic biomarker, and this study provided a framework to identify candidate therapeutic target regions, which further guide glioma patients’ clinical decision.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Feng Li
- *Correspondence: Yunpeng Zhang, ; Yanjun Xu, ; Feng Li,
| | - Yanjun Xu
- *Correspondence: Yunpeng Zhang, ; Yanjun Xu, ; Feng Li,
| | - Yunpeng Zhang
- *Correspondence: Yunpeng Zhang, ; Yanjun Xu, ; Feng Li,
| |
Collapse
|