1
|
Chen L, Wang Y, Cai C, Ding Y, Kim RS, Lipchik C, Gavin PG, Yothers G, Allegra CJ, Petrelli NJ, Suga JM, Hopkins JO, Saito NG, Evans T, Jujjavarapu S, Wolmark N, Lucas PC, Paik S, Sun M, Pogue-Geile KL, Lu X. Machine Learning Predicts Oxaliplatin Benefit in Early Colon Cancer. J Clin Oncol 2024; 42:1520-1530. [PMID: 38315963 PMCID: PMC11095904 DOI: 10.1200/jco.23.01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Revised: 10/12/2023] [Accepted: 11/13/2023] [Indexed: 02/07/2024] Open
Abstract
PURPOSE A combination of fluorouracil, leucovorin, and oxaliplatin (FOLFOX) is the standard for adjuvant therapy of resected early-stage colon cancer (CC). Oxaliplatin leads to lasting and disabling neurotoxicity. Reserving the regimen for patients who benefit from oxaliplatin would maximize efficacy and minimize unnecessary adverse side effects. METHODS We trained a new machine learning model, referred to as the colon oxaliplatin signature (COLOXIS) model, for predicting response to oxaliplatin-containing regimens. We examined whether COLOXIS was predictive of oxaliplatin benefits in the CC adjuvant setting among 1,065 patients treated with 5-fluorouracil plus leucovorin (FULV; n = 421) or FULV + oxaliplatin (FOLFOX; n = 644) from NSABP C-07 and C-08 phase III trials. The COLOXIS model dichotomizes patients into COLOXIS+ (oxaliplatin responder) and COLOXIS- (nonresponder) groups. Eight-year recurrence-free survival was used to evaluate oxaliplatin benefits within each of the groups, and the predictive value of the COLOXIS model was assessed using the P value associated with the interaction term (int P) between the model prediction and the treatment effect. RESULTS Among 1,065 patients, 526 were predicted as COLOXIS+ and 539 as COLOXIS-. The COLOXIS+ prediction was associated with prognosis for FULV-treated patients (hazard ratio [HR], 1.52 [95% CI, 1.07 to 2.15]; P = .017). The model was predictive of oxaliplatin benefits: COLOXIS+ patients benefited from oxaliplatin (HR, 0.65 [95% CI, 0.48 to 0.89]; P = .0065; int P = .03), but COLOXIS- patients did not (COLOXIS- HR, 1.08 [95% CI, 0.77 to 1.52]; P = .65). CONCLUSION The COLOXIS model is predictive of oxaliplatin benefits in the CC adjuvant setting. The results provide evidence supporting a change in CC adjuvant therapy: reserve oxaliplatin only for COLOXIS+ patients, but further investigation is warranted.
Collapse
Affiliation(s)
- Lujia Chen
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA
| | | | - Chunhui Cai
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA
| | - Ying Ding
- NRG Oncology Statistics and Data Management Center, Pittsburgh, PA
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA
| | - Rim S. Kim
- NSABP/NRG Oncology, Pittsburgh, PA
- AstraZeneca, Oncology Translational Medicine, Gaithersburg, MD
| | | | - Patrick G. Gavin
- NSABP/NRG Oncology, Pittsburgh, PA
- AstraZeneca Respiratory and Immunology, Gaithersburg, MD
| | - Greg Yothers
- NRG Oncology Statistics and Data Management Center, Pittsburgh, PA
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA
| | - Carmen J. Allegra
- Department of Medicine, University of Florida Health, Gainesville, FL
| | - Nicholas J. Petrelli
- Helen F. Graham Cancer Center and Research Institute at Christiana Care, Newark, DE
| | - Jennifer Marie Suga
- Kaiser Permanente Oncology Clinical Trials, KP NCI Community Oncology Research Program (NCORP), Vallejo, CA
| | - Judith O. Hopkins
- Novant Health Forsyth Medical Cancer Institute/Southeast Clinical Oncology Research NCORP, Kernersville, NC
| | - Naoyuki G. Saito
- Melvin and Bren Simon Cancer Center, Indiana University School of Medicine, Indianapolis, IN
| | | | | | - Norman Wolmark
- NSABP/NRG Oncology, Pittsburgh, PA
- UPMC Hillman Cancer Center, Pittsburgh, PA
- UPMC Hillman Cancer Center, University of Pittsburgh School of Medicine, Pittsburgh, PA
| | - Peter C. Lucas
- NSABP/NRG Oncology, Pittsburgh, PA
- UPMC Hillman Cancer Center, Pittsburgh, PA
- UPMC Hillman Cancer Center, University of Pittsburgh School of Medicine, Pittsburgh, PA
| | - Soonmyung Paik
- NSABP/NRG Oncology, Pittsburgh, PA
- Yonsei University College of Medicine, Yonsei Biomedical Research Institute, Seoul, Republic of South Korea
| | - Min Sun
- UPMC Hillman Cancer Center, Pittsburgh, PA
- UPMC Hillman Cancer Center, University of Pittsburgh School of Medicine, Pittsburgh, PA
- DeepRx Inc, Pittsburgh, PA
| | | | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA
- DeepRx Inc, Pittsburgh, PA
| |
Collapse
|
2
|
Wei XW, Lu C, Zhang YC, Fan X, Xu CR, Chen ZH, Wang F, Yang XR, Deng JY, Yang MY, Gou Q, Mei SQ, Luo WC, Zhong RW, Zhong WZ, Yang JJ, Zhang XC, Tu HY, Wu YL, Zhou Q. Redox high phenotype mediated by KEAP1/STK11/SMARCA4/NRF2 mutations diminishes tissue-resident memory CD8+ T cells and attenuates the efficacy of immunotherapy in lung adenocarcinoma. Oncoimmunology 2024; 13:2340154. [PMID: 38601319 PMCID: PMC11005803 DOI: 10.1080/2162402x.2024.2340154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 04/03/2024] [Indexed: 04/12/2024] Open
Abstract
Metabolism reprogramming within the tumor microenvironment (TME) can have a profound impact on immune cells. Identifying the association between metabolic phenotypes and immune cells in lung adenocarcinoma (LUAD) may reveal mechanisms of resistance to immune checkpoint inhibitors (ICIs). Metabolic phenotypes were classified by expression of metabolic genes. Somatic mutations and transcriptomic features were compared across the different metabolic phenotypes. The metabolic phenotype of LUAD is predominantly determined by reductase-oxidative activity and is divided into two categories: redoxhigh LUAD and redoxlow LUAD. Genetically, redoxhigh LUAD is mainly driven by mutations in KEAP1, STK11, NRF2, or SMARCA4. These mutations are more prevalent in redoxhigh LUAD (72.5%) compared to redoxlow LUAD (17.4%), whereas EGFR mutations are more common in redoxlow LUAD (19.0% vs. 0.7%). Single-cell RNA profiling of pre-treatment and post-treatment samples from patients receiving neoadjuvant chemoimmunotherapy revealed that tissue-resident memory CD8+ T cells are responders to ICIs. However, these cells are significantly reduced in redoxhigh LUAD. The redoxhigh phenotype is primarily attributed to tumor cells and is positively associated with mTORC1 signaling. LUAD with the redoxhigh phenotype demonstrates a lower response rate (39.1% vs. 70.8%, p = 0.001), shorter progression-free survival (3.3 vs. 14.6 months, p = 0.004), and overall survival (12.1 vs. 31.2 months, p = 0.022) when treated with ICIs. The redoxhigh phenotype in LUAD is predominantly driven by mutations in KEAP1, STK11, NRF2, and SMARCA4. This phenotype diminishes the number of tissue-resident memory CD8+ T cells and attenuates the efficacy of ICIs.
Collapse
Affiliation(s)
- Xue-Wu Wei
- School of Medicine, South China University of Technology, Guangzhou, China
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Chang Lu
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Yi-Chen Zhang
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Xue Fan
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Chong-Rui Xu
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Zhi-Hong Chen
- Guangdong Provincial Key Laboratory of Translational Medicine in Lung Cancer, Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Fen Wang
- Department of Oncology, Shenzhen Key Laboratory of Gastrointestinal Cancer Translational Research, Cancer Institute, Peking University Shenzhen Hospital, Shenzhen-Peking University-Hong Kong University of Science and Technology Medical Center, Shenzhen, China
| | - Xiao-Rong Yang
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Jia-Yi Deng
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Ming-Yi Yang
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Qing Gou
- Department of Interventional Radiology, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Shi-Qi Mei
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Wei-Chi Luo
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Ri-Wei Zhong
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Wen-Zhao Zhong
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Jin-Ji Yang
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Xu-Chao Zhang
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Hai-Yan Tu
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Yi-Long Wu
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| | - Qing Zhou
- School of Medicine, South China University of Technology, Guangzhou, China
- Guangdong Lung Cancer Institute, Guangdong Provincial People’s Hospital (Guangdong Academy of Medical Sciences), Southern Medical University, Guangzhou, China
| |
Collapse
|
3
|
Rahman MA, Cai C, Bo N, McNamara DM, Ding Y, Cooper GF, Lu X, Liu J. An individualized Bayesian method for estimating genomic variants of hypertension. BMC Genomics 2023; 23:863. [PMID: 37936055 PMCID: PMC10631115 DOI: 10.1186/s12864-023-09757-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 10/19/2023] [Indexed: 11/09/2023] Open
Abstract
BACKGROUND Genomic variants of the disease are often discovered nowadays through population-based genome-wide association studies (GWAS). Identifying genomic variations potentially underlying a phenotype, such as hypertension, in an individual is important for designing personalized treatment; however, population-level models, such as GWAS, may not capture all the important, individualized factors well. In addition, GWAS typically requires a large sample size to detect the association of low-frequency genomic variants with sufficient power. Here, we report an individualized Bayesian inference (IBI) algorithm for estimating the genomic variants that influence complex traits, such as hypertension, at the level of an individual (e.g., a patient). By modeling at the level of the individual, IBI seeks to find genomic variants observed in the individual's genome that provide a strong explanation of the phenotype observed in this individual. RESULTS We applied the IBI algorithm to the data from the Framingham Heart Study to explore the genomic influences of hypertension. Among the top-ranking variants identified by IBI and GWAS, there is a significant number of shared variants (intersection); the unique variants identified only by IBI tend to have relatively lower minor allele frequency than those identified by GWAS. In addition, IBI discovered more individualized and diverse variants that explain hypertension patients better than GWAS. Furthermore, IBI found several well-known low-frequency variants as well as genes related to blood pressure that GWAS missed in the same cohort. Finally, IBI identified top-ranked variants that predicted hypertension better than GWAS, according to the area under the ROC curve. CONCLUSIONS The results support IBI as a promising approach for complementing GWAS, especially in detecting low-frequency genomic variants as well as learning personalized genomic variants of clinical traits and disease, such as the complex trait of hypertension, to help advance precision medicine.
Collapse
Affiliation(s)
- Md Asad Rahman
- Department of Engineering Management and Systems Engineering, Missouri University of Science and Technology, Rolla, MO, USA
| | - Chunhui Cai
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Na Bo
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Dennis M McNamara
- Department of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Ying Ding
- Department of Biostatistics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Gregory F Cooper
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA, USA
| | - Jinling Liu
- Department of Engineering Management and Systems Engineering, Missouri University of Science and Technology, Rolla, MO, USA.
- Department of Biological Sciences, Missouri University of Science and Technology, Rolla, MO, USA.
- Department of Epidemiology, College of Public Health and Health Professions and College of Medicine, University of Florida, Gainesville, FL, USA.
| |
Collapse
|
4
|
Young JD, Ren S, Chen L, Lu X. Revealing the Impact of Genomic Alterations on Cancer Cell Signaling with an Interpretable Deep Learning Model. Cancers (Basel) 2023; 15:3857. [PMID: 37568673 PMCID: PMC10416927 DOI: 10.3390/cancers15153857] [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/30/2023] [Revised: 07/24/2023] [Accepted: 07/26/2023] [Indexed: 08/13/2023] Open
Abstract
Cancer is a disease of aberrant cellular signaling resulting from somatic genomic alterations (SGAs). Heterogeneous SGA events in tumors lead to tumor-specific signaling system aberrations. We interpret the cancer signaling system as a causal graphical model, where SGAs affect signaling proteins, propagate their effects through signal transduction, and ultimately change gene expression. To represent such a system, we developed a deep learning model called redundant-input neural network (RINN) with a transparent redundant-input architecture. Our findings demonstrate that by utilizing SGAs as inputs, the RINN can encode their impact on the signaling system and predict gene expression accurately when measured as the area under ROC curves. Moreover, the RINN can discover the shared functional impact (similar embeddings) of SGAs that perturb a common signaling pathway (e.g., PI3K, Nrf2, and TGF). Furthermore, the RINN exhibits the ability to discover known relationships in cellular signaling systems.
Collapse
Affiliation(s)
- Jonathan D. Young
- Intelligent Systems Program, School of Computing and Information, University of Pittsburgh, Pittsburgh, PA 15260, USA;
| | - Shuangxia Ren
- Intelligent Systems Program, School of Computing and Information, University of Pittsburgh, Pittsburgh, PA 15260, USA;
| | - Lujia Chen
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA;
| | - Xinghua Lu
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA;
| |
Collapse
|
5
|
Ren S, Cooper GF, Chen L, Lu X. An interpretable deep learning framework for genome-informed precision oncology. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.11.548534. [PMID: 37503199 PMCID: PMC10369905 DOI: 10.1101/2023.07.11.548534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Cancers result from aberrations in cellular signaling systems, typically resulting from driver somatic genome alterations (SGAs) in individual tumors. Precision oncology requires understanding the cellular state and selecting medications that induce vulnerability in cancer cells under such conditions. To this end, we developed a computational framework consisting of two components: 1) A representation-learning component, which learns a representation of the cellular signaling systems when perturbed by SGAs, using a biologically-motivated and interpretable deep learning model. 2) A drug-response-prediction component, which predicts the response to drugs by leveraging the information of the cellular state of the cancer cells derived by the first component. Our cell-state-oriented framework significantly enhances the accuracy of genome-informed prediction of drug responses in comparison to models that directly use SGAs as inputs. Importantly, our framework enables the prediction of response to chemotherapy agents based on SGAs, thus expanding genome-informed precision oncology beyond molecularly targeted drugs.
Collapse
Affiliation(s)
- Shuangxia Ren
- Intelligent Systems Program, School of Computing and Information, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Gregory F Cooper
- Intelligent Systems Program, School of Computing and Information, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Lujia Chen
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Xinghua Lu
- Intelligent Systems Program, School of Computing and Information, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| |
Collapse
|
6
|
Bai G, Li Y, Ji Y, Peng Y, Yang Z, Zhao L. Treatments and whole exon sequencing of a case with multiple primary lung cancer. J Cardiothorac Surg 2023; 18:58. [PMID: 36732778 PMCID: PMC9896780 DOI: 10.1186/s13019-023-02161-0] [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: 11/05/2022] [Accepted: 01/24/2023] [Indexed: 02/04/2023] Open
Abstract
INTRODUCTION The number of patients with synchronous multiple primary lung cancer (sMPLC) has increased recently. However, diagnosing and selecting the appropriate therapeutic strategy for this type of disease is not simple. CASE PRESENTATION This report presented a case of sMPLC with lymph node metastasis. With no smoking and cancer history, this patient had seven nodules in the right lung and underwent single-portal video-assisted thoracoscopic surgery (VATS). In addition, she received four cycles of chemotherapy after the operation. Whole exon sequencing (WES) was performed in five resected tissue samples (four tumors and one lymph node). We conducted genomic profiling and clone evolution analysis of the five samples. Gene detection helped to confirm that the metastasis lymph node was transferred from one nodule. There was apparent heterogeneity of gene mutations among the five samples of the patient, with only one shared "neurofilament heavy polypeptide" (NEFH) mutation. A dominant substitution of C > T/G > A was found in all the samples. Pyclone model was used to calculate all tissues' cellular prevalence (CP) values, and NEFH mutations were thought to be the ancestral clones. During the follow-up period, residual lesions showed no apparent changes and limited response to chemotherapy. CONCLUSIONS This report showed an essential role in genomic detection and selecting the appropriate treatment of sMPLC. Surgery remains the primary treatment strategy for this type of disease, and the occurrence and development of sMPLC need more in-depth research.
Collapse
Affiliation(s)
- Guangyu Bai
- grid.506261.60000 0001 0706 7839Department of Thoracic Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, 100021 China
| | - Yuan Li
- grid.506261.60000 0001 0706 7839Department of Thoracic Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, 100021 China
| | - Ying Ji
- grid.24696.3f0000 0004 0369 153XDepartment of Thoracic Surgery, Beijing Chao-Yang Hospital, Capital Medical University, Beijing, 100020 China
| | - Yue Peng
- grid.24696.3f0000 0004 0369 153XDepartment of Thoracic Surgery, Beijing Chao-Yang Hospital, Capital Medical University, Beijing, 100020 China
| | - Zhenlin Yang
- grid.506261.60000 0001 0706 7839Department of Thoracic Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, 100021 China
| | - Liang Zhao
- grid.506261.60000 0001 0706 7839Department of Thoracic Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, 100021 China
| |
Collapse
|
7
|
Chen X, Chen L, Kürten CHL, Jabbari F, Vujanovic L, Ding Y, Lu B, Lu K, Kulkarni A, Tabib T, Lafyatis R, Cooper GF, Ferris R, Lu X. An individualized causal framework for learning intercellular communication networks that define microenvironments of individual tumors. PLoS Comput Biol 2022; 18:e1010761. [PMID: 36548438 PMCID: PMC9822106 DOI: 10.1371/journal.pcbi.1010761] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Revised: 01/06/2023] [Accepted: 11/26/2022] [Indexed: 12/24/2022] Open
Abstract
Cells within a tumor microenvironment (TME) dynamically communicate and influence each other's cellular states through an intercellular communication network (ICN). In cancers, intercellular communications underlie immune evasion mechanisms of individual tumors. We developed an individualized causal analysis framework for discovering tumor specific ICNs. Using head and neck squamous cell carcinoma (HNSCC) tumors as a testbed, we first mined single-cell RNA-sequencing data to discover gene expression modules (GEMs) that reflect the states of transcriptomic processes within tumor and stromal single cells. By deconvoluting bulk transcriptomes of HNSCC tumors profiled by The Cancer Genome Atlas (TCGA), we estimated the activation states of these transcriptomic processes in individual tumors. Finally, we applied individualized causal network learning to discover an ICN within each tumor. Our results show that cellular states of cells in TMEs are coordinated through ICNs that enable multi-way communications among epithelial, fibroblast, endothelial, and immune cells. Further analyses of individual ICNs revealed structural patterns that were shared across subsets of tumors, leading to the discovery of 4 different subtypes of networks that underlie disparate TMEs of HNSCC. Patients with distinct TMEs exhibited significantly different clinical outcomes. Our results show that the capability of estimating individual ICNs reveals heterogeneity of ICNs and sheds light on the importance of intercellular communication in impacting disease development and progression.
Collapse
Affiliation(s)
- Xueer Chen
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Causal Discovery, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Lujia Chen
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Causal Discovery, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Cornelius H. L. Kürten
- Department of Otolaryngology, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
- University of Pittsburgh Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Department of Otorhinolaryngology, Head and Neck Surgery, University Hospital Essen, University Duisburg-Essen, Duisburg, Germany
| | - Fattaneh Jabbari
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Causal Discovery, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Lazar Vujanovic
- Department of Otolaryngology, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
- University of Pittsburgh Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Ying Ding
- Department of Biostatistics, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Binfeng Lu
- Department of Immunology, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Kevin Lu
- Williamsville North High School, Williamsville, New York, United States of America
| | - Aditi Kulkarni
- Department of Otolaryngology, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
| | - Tracy Tabib
- Department of Otorhinolaryngology, Head and Neck Surgery, University Hospital Essen, University Duisburg-Essen, Duisburg, Germany
| | - Robert Lafyatis
- Department of Otorhinolaryngology, Head and Neck Surgery, University Hospital Essen, University Duisburg-Essen, Duisburg, Germany
| | - Gregory F. Cooper
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Causal Discovery, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
- University of Pittsburgh Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Robert Ferris
- Department of Otolaryngology, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
- University of Pittsburgh Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
- Center for Causal Discovery, University of Pittsburgh, Pennsylvania, Pittsburgh, United States of America
- University of Pittsburgh Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| |
Collapse
|
8
|
Genome wide DNA methylation analysis identifies novel molecular subgroups and predicts survival in neuroblastoma. Br J Cancer 2022; 127:2006-2015. [PMID: 36175618 PMCID: PMC9681858 DOI: 10.1038/s41416-022-01988-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 08/22/2022] [Accepted: 09/08/2022] [Indexed: 11/09/2022] Open
Abstract
BACKGROUND Neuroblastoma is the most common malignancy in infancy, accounting for 15% of childhood cancer deaths. Outcome for the high-risk disease remains poor. DNA-methylation patterns are significantly altered in all cancer types and can be utilised for disease stratification. METHODS Genome-wide DNA methylation (n = 223), gene expression (n = 130), genetic/clinical data (n = 213), whole-exome sequencing (n = 130) was derived from the TARGET study. Methylation data were derived from HumanMethylation450 BeadChip arrays. t-SNE was used for the segregation of molecular subgroups. A separate validation cohort of 105 cases was studied. RESULTS Five distinct neuroblastoma molecular subgroups were identified, based on genome-wide DNA-methylation patterns, with unique features in each, including three subgroups associated with known prognostic features and two novel subgroups. As expected, Cluster-4 (infant diagnosis) had significantly better 5-year progression-free survival (PFS) than the four other clusters. However, in addition, the molecular subgrouping identified multiple patient subsets with highly increased risk, most notably infant patients that do not map to Cluster-4 (PFS 50% vs 80% for Cluster-4 infants, P = 0.005), and allowed identification of subgroup-specific methylation differences that may reflect important biological differences within neuroblastoma. CONCLUSIONS Methylation-based clustering of neuroblastoma reveals novel molecular subgroups, with distinct molecular/clinical characteristics and identifies a subgroup of higher-risk infant patients.
Collapse
|
9
|
Tao Y, Ma X, Palmer D, Schwartz R, Lu X, Osmanbeyoglu H. Interpretable deep learning for chromatin-informed inference of transcriptional programs driven by somatic alterations across cancers. Nucleic Acids Res 2022; 50:10869-10881. [PMID: 36243974 PMCID: PMC9638905 DOI: 10.1093/nar/gkac881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 09/23/2022] [Accepted: 09/29/2022] [Indexed: 11/14/2022] Open
Abstract
Cancer is a disease of gene dysregulation, where cells acquire somatic and epigenetic alterations that drive aberrant cellular signaling. These alterations adversely impact transcriptional programs and cause profound changes in gene expression. Interpreting somatic alterations within context-specific transcriptional programs will facilitate personalized therapeutic decisions but is a monumental task. Toward this goal, we develop a partially interpretable neural network model called Chromatin-informed Inference of Transcriptional Regulators Using Self-attention mechanism (CITRUS). CITRUS models the impact of somatic alterations on transcription factors and downstream transcriptional programs. Our approach employs a self-attention mechanism to model the contextual impact of somatic alterations. Furthermore, CITRUS uses a layer of hidden nodes to explicitly represent the state of transcription factors (TFs) to learn the relationships between TFs and their target genes based on TF binding motifs in the open chromatin regions of tumor samples. We apply CITRUS to genomic, transcriptomic, and epigenomic data from 17 cancer types profiled by The Cancer Genome Atlas. CITRUS predicts patient-specific TF activities and reveals transcriptional program variations between and within tumor types. We show that CITRUS yields biological insights into delineating TFs associated with somatic alterations in individual tumors. Thus, CITRUS is a promising tool for precision oncology.
Collapse
Affiliation(s)
- Yifeng Tao
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xiaojun Ma
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
- UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA
| | - Drake Palmer
- UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA
| | - Russell Schwartz
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA, USA
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Xinghua Lu
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
- Department of Pharmaceutical Science, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
| | - Hatice Ulku Osmanbeyoglu
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
- UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA, USA
- Department of Bioengineering, School of Engineering, University of Pittsburgh, Pittsburgh, PA, USA
- Department of Biostatistics, School of Public Health, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
10
|
Liu Z, Cai C, Ma X, Liu J, Chen L, Lui VWY, Cooper GF, Lu X. A Novel Bayesian Framework Infers Driver Activation States and Reveals Pathway-Oriented Molecular Subtypes in Head and Neck Cancer. Cancers (Basel) 2022; 14:cancers14194825. [PMID: 36230748 PMCID: PMC9563147 DOI: 10.3390/cancers14194825] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 09/28/2022] [Accepted: 09/30/2022] [Indexed: 02/08/2023] Open
Abstract
Head and neck squamous cell cancer (HNSCC) is an aggressive cancer resulting from heterogeneous causes. To reveal the underlying drivers and signaling mechanisms of different HNSCC tumors, we developed a novel Bayesian framework to identify drivers of individual tumors and infer the states of driver proteins in cellular signaling system in HNSCC tumors. First, we systematically identify causal relationships between somatic genome alterations (SGAs) and differentially expressed genes (DEGs) for each TCGA HNSCC tumor using the tumor-specific causal inference (TCI) model. Then, we generalize the most statistically significant driver SGAs and their regulated DEGs in TCGA HNSCC cohort. Finally, we develop machine learning models that combine genomic and transcriptomic data to infer the protein functional activation states of driver SGAs in tumors, which enable us to represent a tumor in the space of cellular signaling systems. We discovered four mechanism-oriented subtypes of HNSCC, which show distinguished patterns of activation state of HNSCC driver proteins, and importantly, this subtyping is orthogonal to previously reported transcriptomic-based molecular subtyping of HNSCC. Further, our analysis revealed driver proteins that are likely involved in oncogenic processes induced by HPV infection, even though they are not perturbed by genomic alterations in HPV+ tumors.
Collapse
Affiliation(s)
- Zhengping Liu
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
- School of Medicine, Tsinghua University, Beijing 100190, China
| | - Chunhui Cai
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
- Correspondence:
| | - Xiaojun Ma
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
| | - Jinling Liu
- Department of Engineering Management and Systems Engineering, Missouri University of Science and Technology, Rolla, MO 65409, USA
- Department of Biological Sciences, Missouri University of Science and Technology, Rolla, MO 65409, USA
| | - Lujia Chen
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
| | - Vivian Wai Yan Lui
- Georgia Cancer Center, and Department of Medicine, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| | - Gregory F. Cooper
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
- UPMC Hillman Cancer Center, University of Pittsburgh Medical Center, Pittsburgh, PA 15232, USA
| | - Xinghua Lu
- Department of Biomedical Informatics, School of Medicine, University of Pittsburgh, Pittsburgh 15206, PA, USA
- UPMC Hillman Cancer Center, University of Pittsburgh Medical Center, Pittsburgh, PA 15232, USA
| |
Collapse
|
11
|
Dressler L, Bortolomeazzi M, Keddar MR, Misetic H, Sartini G, Acha-Sagredo A, Montorsi L, Wijewardhane N, Repana D, Nulsen J, Goldman J, Pollitt M, Davis P, Strange A, Ambrose K, Ciccarelli FD. Comparative assessment of genes driving cancer and somatic evolution in non-cancer tissues: an update of the Network of Cancer Genes (NCG) resource. Genome Biol 2022; 23:35. [PMID: 35078504 PMCID: PMC8790917 DOI: 10.1186/s13059-022-02607-z] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Accepted: 01/10/2022] [Indexed: 12/30/2022] Open
Abstract
Background Genetic alterations of somatic cells can drive non-malignant clone formation and promote cancer initiation. However, the link between these processes remains unclear and hampers our understanding of tissue homeostasis and cancer development. Results Here, we collect a literature-based repertoire of 3355 well-known or predicted drivers of cancer and non-cancer somatic evolution in 122 cancer types and 12 non-cancer tissues. Mapping the alterations of these genes in 7953 pan-cancer samples reveals that, despite the large size, the known compendium of drivers is still incomplete and biased towards frequently occurring coding mutations. High overlap exists between drivers of cancer and non-cancer somatic evolution, although significant differences emerge in their recurrence. We confirm and expand the unique properties of drivers and identify a core of evolutionarily conserved and essential genes whose germline variation is strongly counter-selected. Somatic alteration in even one of these genes is sufficient to drive clonal expansion but not malignant transformation. Conclusions Our study offers a comprehensive overview of our current understanding of the genetic events initiating clone expansion and cancer revealing significant gaps and biases that still need to be addressed. The compendium of cancer and non-cancer somatic drivers, their literature support, and properties are accessible in the Network of Cancer Genes and Healthy Drivers resource at http://www.network-cancer-genes.org/. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-022-02607-z.
Collapse
|
12
|
Whole-exome sequencing in eccrine porocarcinoma indicates promising therapeutic strategies. Cancer Gene Ther 2022; 29:697-708. [PMID: 34045664 PMCID: PMC9209330 DOI: 10.1038/s41417-021-00347-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Revised: 04/23/2021] [Accepted: 05/10/2021] [Indexed: 02/06/2023]
Abstract
Malignant sweat gland tumours are rare, with the most common form being Eccrine porocarcinoma (EP). To investigate the mutational landscape of EP, we performed whole-exome sequencing (WES) on 14 formalin-fixed paraffin-embedded samples of matched primary EP and healthy surrounding tissue. Mutational profiling revealed a high overall median mutation rate. This was attributed to signatures of mutational processes related to ultraviolet (UV) exposure, APOBEC enzyme dysregulation, and defective homologous double-strand break repair. All of these processes cause genomic instability and are implicated in carcinogenesis. Recurrent driving somatic alterations were detected in the EP candidate drivers TP53, FAT2, CACNA1S, and KMT2D. The analyses also identified copy number alterations and recurrent gains and losses in several chromosomal regions including that containing BRCA2, as well as deleterious alterations in multiple HRR components. In accordance with this reduced or even a complete loss of BRCA2 protein expression was detected in 50% of the investigated EP tumours. Our results implicate crucial oncogenic driver pathways and suggest that defective homologous double-strand break repair and the p53 pathway are involved in EP aetiology. Targeting of the p53 axis and PARP inhibition, and/or immunotherapy may represent promising treatment strategies.
Collapse
|
13
|
Fang Y, Wang Y, Zeng D, Zhi S, Shu T, Huang N, Zheng S, Wu J, Liu Y, Huang G, Xue Y, Bin J, Liao Y, Shi M, Liao W. Comprehensive analyses reveal TKI-induced remodeling of the tumor immune microenvironment in EGFR/ALK-positive non-small-cell lung cancer. Oncoimmunology 2021; 10:1951019. [PMID: 34345533 PMCID: PMC8288040 DOI: 10.1080/2162402x.2021.1951019] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022] Open
Abstract
Tyrosine kinase inhibitors (TKI) play a pivotal role in the treatment of non-small-cell lung cancer (NSCLC) with mutations in epidermal growth factor receptor (EGFR) and rearrangements in anaplastic lymphoma kinase (ALK). However, the influences of TKIs on the tumor immune microenvironment (TIM), especially dynamic changes of responders, have not yet been fully elucidated. Therefore, RNA sequencing and whole-exome sequencing were performed on EGFR/ALK-positive NSCLC samples before and after TKI treatment. In combination with neoantigen and mutational-load estimations, xCell and single-sample gene set enrichment analysis (ssGSEA) were used to assess tumor immune-cell infiltration and activity. Furthermore, weighted-gene correlation network analysis and the bottleneck method were used to identify the hub genes that affected treatment-related immune responses. We found that TKI treatment remodeled the TIM in treatment-responsive samples. Profound increases in the rate of anti-tumor cell infiltration and cytotoxicity was observed following TKI treatment, while antigen presentation was limited in ALK-rearranged samples. However, no significant change in anti-tumor cell infiltration or cytotoxicity was found between pre-treatment and post-progression samples. Subsequently, we found that neurofilament heavy (NEFH) mutations were enriched in samples after TKI treatment and were associated with reduced neutrophil infiltration. The cytotoxicity of EGFR-mutant NSCLCs with co-driver TP53 mutation and ALK-rearranged samples with wild-type TP53 seems to be more easily induced by TKI. Finally, the immune-associated score generated by hub genes was positively correlated with immune infiltration, immune activation, and a favorable prognosis. In conclusion, the dynamic changes in the TIM provide clues to drug selection and timing for TKI-immunotherapy combinations.
Collapse
Affiliation(s)
- Yisheng Fang
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Yuanyuan Wang
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Dongqiang Zeng
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Shimeng Zhi
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Tingting Shu
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Na Huang
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Siting Zheng
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Jianhua Wu
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Yantan Liu
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Genjie Huang
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Yichen Xue
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Jianping Bin
- Department of Cardiology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Yulin Liao
- Department of Cardiology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Min Shi
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| | - Wangjun Liao
- Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou, Guangdong, China
| |
Collapse
|
14
|
Efficient representations of tumor diversity with paired DNA-RNA aberrations. PLoS Comput Biol 2021; 17:e1008944. [PMID: 34115745 PMCID: PMC8221796 DOI: 10.1371/journal.pcbi.1008944] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 06/23/2021] [Accepted: 04/07/2021] [Indexed: 12/13/2022] Open
Abstract
Cancer cells display massive dysregulation of key regulatory pathways due to now well-catalogued mutations and other DNA-related aberrations. Moreover, enormous heterogeneity has been commonly observed in the identity, frequency and location of these aberrations across individuals with the same cancer type or subtype, and this variation naturally propagates to the transcriptome, resulting in myriad types of dysregulated gene expression programs. Many have argued that a more integrative and quantitative analysis of heterogeneity of DNA and RNA molecular profiles may be necessary for designing more systematic explorations of alternative therapies and improving predictive accuracy. We introduce a representation of multi-omics profiles which is sufficiently rich to account for observed heterogeneity and support the construction of quantitative, integrated, metrics of variation. Starting from the network of interactions existing in Reactome, we build a library of “paired DNA-RNA aberrations” that represent prototypical and recurrent patterns of dysregulation in cancer; each two-gene “Source-Target Pair” (STP) consists of a “source” regulatory gene and a “target” gene whose expression is plausibly “controlled” by the source gene. The STP is then “aberrant” in a joint DNA-RNA profile if the source gene is DNA-aberrant (e.g., mutated, deleted, or duplicated), and the downstream target gene is “RNA-aberrant”, meaning its expression level is outside the normal, baseline range. With M STPs, each sample profile has exactly one of the 2M possible configurations. We concentrate on subsets of STPs, and the corresponding reduced configurations, by selecting tissue-dependent minimal coverings, defined as the smallest family of STPs with the property that every sample in the considered population displays at least one aberrant STP within that family. These minimal coverings can be computed with integer programming. Given such a covering, a natural measure of cross-sample diversity is the extent to which the particular aberrant STPs composing a covering vary from sample to sample; this variability is captured by the entropy of the distribution over configurations. We apply this program to data from TCGA for six distinct tumor types (breast, prostate, lung, colon, liver, and kidney cancer). This enables an efficient simplification of the complex landscape observed in cancer populations, resulting in the identification of novel signatures of molecular alterations which are not detected with frequency-based criteria. Estimates of cancer heterogeneity across tumor phenotypes reveals a stable pattern: entropy increases with disease severity. This framework is then well-suited to accommodate the expanding complexity of cancer genomes and epigenomes emerging from large consortia projects. A large variety of genomic and transcriptomic aberrations are observed in cancer cells, and their identity, location, and frequency can be highly indicative of the particular subtype or molecular phenotype, and thereby inform treatment options. However, elucidating this association between sets of aberrations and subtypes of cancer is severely impeded by considerable diversity in the set of aberrations across samples from the same population. Most attempts at analyzing tumor heterogeneity have dealt with either the genome or transcriptome in isolation. Here we present a novel, multi-omics approach for quantifying heterogeneity by determining a small set of paired DNA-RNA aberrations that incorporates potential downstream effects on gene expression. We apply integer programming to identify a small set of paired aberrations such that at least one among them is present in every sample of a given cancer population. The resulting “coverings” are analyzed for six cancer cohorts from the Cancer Genome Atlas, and facilitate introducing an information-theoretic measure of heterogeneity. Our results identify many known facets of tumorigenesis as well as suggest potential novel genes and interactions of interest.
Collapse
|
15
|
The Role of Fibroblast Growth Factor 19 in Hepatocellular Carcinoma. THE AMERICAN JOURNAL OF PATHOLOGY 2021; 191:1180-1192. [PMID: 34000282 DOI: 10.1016/j.ajpath.2021.04.014] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Revised: 04/09/2021] [Accepted: 04/22/2021] [Indexed: 12/12/2022]
Abstract
Hepatocellular carcinoma (HCC) is the fifth most common type of cancer and the third leading cause of cancer-related deaths worldwide. Liver resection or liver transplantation is the most effective therapy for HCC because drugs approved by the US Food and Drug Administration to treat patients with unresectable HCC have an unfavorable overall survival rate. Therefore, the development of biomarkers for early diagnosis and effective therapy strategies are still necessary to improve patient outcomes. Fibroblast growth factor (FGF) 19 was amplified in patients with HCC from various studies, including patients from The Cancer Genome Atlas. FGF19 plays a syngeneic function with other signaling pathways in primary liver cancer development, such as epidermal growth factor receptor, Wnt/β-catenin, the endoplasmic reticulum-related signaling pathway, STAT3/IL-6, RAS, and extracellular signal-regulated protein kinase, among others. The current review presents a comprehensive description of the FGF19 signaling pathway involved in liver cancer development. The use of big data and bioinformatic analysis can provide useful clues for further studies of the FGF19 pathway in HCC, including its application as a biomarker, targeted therapy, and combination therapy strategies.
Collapse
|
16
|
Grzadkowski MR, Holly HD, Somers J, Demir E. Systematic interrogation of mutation groupings reveals divergent downstream expression programs within key cancer genes. BMC Bioinformatics 2021; 22:233. [PMID: 33957863 PMCID: PMC8101181 DOI: 10.1186/s12859-021-04147-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Accepted: 04/22/2021] [Indexed: 12/15/2022] Open
Abstract
Background Genes implicated in tumorigenesis often exhibit diverse sets of genomic variants in the tumor cohorts within which they are frequently mutated. For many genes, neither the transcriptomic effects of these variants nor their relationship to one another in cancer processes have been well-characterized. We sought to identify the downstream expression effects of these mutations and to determine whether this heterogeneity at the genomic level is reflected in a corresponding heterogeneity at the transcriptomic level. Results By applying a novel hierarchical framework for organizing the mutations present in a cohort along with machine learning pipelines trained on samples’ expression profiles we systematically interrogated the signatures associated with combinations of mutations recurrent in cancer. This allowed us to catalogue the mutations with discernible downstream expression effects across a number of tumor cohorts as well as to uncover and characterize over a hundred cases where subsets of a gene’s mutations are clearly divergent in their function from the remaining mutations of the gene. These findings successfully replicated across a number of disease contexts and were found to have clear implications for the delineation of cancer processes and for clinical decisions. Conclusions The results of cataloguing the downstream effects of mutation subgroupings across cancer cohorts underline the importance of incorporating the diversity present within oncogenes in models designed to capture the downstream effects of their mutations. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04147-y.
Collapse
Affiliation(s)
- Michal R Grzadkowski
- Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR, USA.
| | - Hannah D Holly
- Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR, USA
| | - Julia Somers
- Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR, USA
| | - Emek Demir
- Department of Molecular and Medical Genetics, Oregon Health & Science University, Portland, OR, USA
| |
Collapse
|
17
|
Spatial Distribution of Private Gene Mutations in Clear Cell Renal Cell Carcinoma. Cancers (Basel) 2021; 13:cancers13092163. [PMID: 33946379 PMCID: PMC8124666 DOI: 10.3390/cancers13092163] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Revised: 04/02/2021] [Accepted: 04/27/2021] [Indexed: 12/15/2022] Open
Abstract
Simple Summary Tumours consist of multiple groups of similar cells resulting from differing evolutionary trajectories, i.e., subclones. These subclones are prevalent in clear cell renal cell carcinoma (ccRCC). The aim of this study is to determine how similar or dissimilar the subclones in 89 ccRCC tumours are from one another regarding their gene mutations and expression profiles, i.e., the extent of intra-tumour heterogeneity. The implications of these alterations with respect to signalling pathways is also assessed. Deep sequencing allows for the identification of mutations with low-allele frequencies, providing a more comprehensive view of the heterogeneity present in the tumours. With an average of 62% of mutations having been identified in only one of the two biopsies, some of which in turn are found to impact gene expression, the complex makeup of ccRCC tumours is evident, and this can drastically influence treatment outcome. Abstract Intra-tumour heterogeneity is the molecular hallmark of renal cancer, and the molecular tumour composition determines the treatment outcome of renal cancer patients. In renal cancer tumourigenesis, in general, different tumour clones evolve over time. We analysed intra-tumour heterogeneity and subclonal mutation patterns in 178 tumour samples obtained from 89 clear cell renal cell carcinoma patients. In an initial discovery phase, whole-exome and transcriptome sequencing data from paired tumour biopsies from 16 ccRCC patients were used to design a gene panel for follow-up analysis. In this second phase, 826 selected genes were targeted at deep coverage in an extended cohort of 89 patients for a detailed analysis of tumour heterogeneity. On average, we found 22 mutations per patient. Pairwise comparison of the two biopsies from the same tumour revealed that on average, 62% of the mutations in a patient were detected in one of the two samples. In addition to commonly mutated genes (VHL, PBRM1, SETD2 and BAP1), frequent subclonal mutations with low variant allele frequency (<10%) were observed in TP53 and in mucin coding genes MUC6, MUC16, and MUC3A. Of the 89 ccRCC tumours, 87 (~98%) harboured private mutations, occurring in only one of the paired tumour samples. Clonally exclusive pathway pairs were identified using the WES data set from 16 ccRCC patients. Our findings imply that shared and private mutations significantly contribute to the complexity of differential gene expression and pathway interaction and might explain the clonal evolution of different molecular renal cancer subgroups. Multi-regional sequencing is central for the identification of subclones within ccRCC.
Collapse
|
18
|
Lopes MB, Martins EP, Vinga S, Costa BM. The Role of Network Science in Glioblastoma. Cancers (Basel) 2021; 13:1045. [PMID: 33801334 PMCID: PMC7958335 DOI: 10.3390/cancers13051045] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Revised: 02/19/2021] [Accepted: 02/22/2021] [Indexed: 12/13/2022] Open
Abstract
Network science has long been recognized as a well-established discipline across many biological domains. In the particular case of cancer genomics, network discovery is challenged by the multitude of available high-dimensional heterogeneous views of data. Glioblastoma (GBM) is an example of such a complex and heterogeneous disease that can be tackled by network science. Identifying the architecture of molecular GBM networks is essential to understanding the information flow and better informing drug development and pre-clinical studies. Here, we review network-based strategies that have been used in the study of GBM, along with the available software implementations for reproducibility and further testing on newly coming datasets. Promising results have been obtained from both bulk and single-cell GBM data, placing network discovery at the forefront of developing a molecularly-informed-based personalized medicine.
Collapse
Affiliation(s)
- Marta B. Lopes
- Center for Mathematics and Applications (CMA), FCT, UNL, 2829-516 Caparica, Portugal
- NOVA Laboratory for Computer Science and Informatics (NOVA LINCS), FCT, UNL, 2829-516 Caparica, Portugal
| | - Eduarda P. Martins
- Life and Health Sciences Research Institute (ICVS), School of Medicine, University of Minho, Campus de Gualtar, 4710-057 Braga, Portugal; (E.P.M.); (B.M.C.)
- ICVS/3B’s—PT Government Associate Laboratory, 4710-057/4805-017 Braga/Guimarães, Portugal
| | - Susana Vinga
- INESC-ID, Instituto Superior Técnico, Universidade de Lisboa, 1000-029 Lisbon, Portugal;
- IDMEC, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisbon, Portugal
| | - Bruno M. Costa
- Life and Health Sciences Research Institute (ICVS), School of Medicine, University of Minho, Campus de Gualtar, 4710-057 Braga, Portugal; (E.P.M.); (B.M.C.)
- ICVS/3B’s—PT Government Associate Laboratory, 4710-057/4805-017 Braga/Guimarães, Portugal
| |
Collapse
|
19
|
Tao Y, Cai C, Cohen WW, Lu X. From genome to phenome: Predicting multiple cancer phenotypes based on somatic genomic alterations via the genomic impact transformer. PACIFIC SYMPOSIUM ON BIOCOMPUTING. PACIFIC SYMPOSIUM ON BIOCOMPUTING 2020; 25:79-90. [PMID: 31797588 PMCID: PMC6932864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cancers are mainly caused by somatic genomic alterations (SGAs) that perturb cellular signaling systems and eventually activate oncogenic processes. Therefore, understanding the functional impact of SGAs is a fundamental task in cancer biology and precision oncology. Here, we present a deep neural network model with encoder-decoder architecture, referred to as genomic impact transformer (GIT), to infer the functional impact of SGAs on cellular signaling systems through modeling the statistical relationships between SGA events and differentially expressed genes (DEGs) in tumors. The model utilizes a multi-head self-attention mechanism to identify SGAs that likely cause DEGs, or in other words, differentiating potential driver SGAs from passenger ones in a tumor. GIT model learns a vector (gene embedding) as an abstract representation of functional impact for each SGA-affected gene. Given SGAs of a tumor, the model can instantiate the states of the hidden layer, providing an abstract representation (tumor embedding) reflecting characteristics of perturbed molecular/cellular processes in the tumor, which in turn can be used to predict multiple phenotypes. We apply the GIT model to 4,468 tumors profiled by The Cancer Genome Atlas (TCGA) project. The attention mechanism enables the model to better capture the statistical relationship between SGAs and DEGs than conventional methods, and distinguishes cancer drivers from passengers. The learned gene embeddings capture the functional similarity of SGAs perturbing common pathways. The tumor embeddings are shown to be useful for tumor status representation, and phenotype prediction including patient survival time and drug response of cancer cell lines.
Collapse
Affiliation(s)
- Yifeng Tao
- School of Computer Science, Carnegie Mellon University
| | - Chunhui Cai
- Department of Biomedical Informatics, University of Pittsburgh
| | - William W. Cohen
- School of Computer Science, Carnegie Mellon University,To whom correspondence should be addressed. ,
| | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh,Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA, USA,To whom correspondence should be addressed. ,
| |
Collapse
|
20
|
Jabbari F, Villaruz LC, Davis M, Cooper GF. Lung Cancer Survival Prediction Using Instance-Specific Bayesian Networks. Artif Intell Med 2020. [DOI: 10.1007/978-3-030-59137-3_14] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
21
|
Xue Y, Cooper G, Cai C, Lu S, Hu B, Ma X, Lu X. Tumour-specific Causal Inference Discovers Distinct Disease Mechanisms Underlying Cancer Subtypes. Sci Rep 2019; 9:13225. [PMID: 31519988 PMCID: PMC6744493 DOI: 10.1038/s41598-019-48318-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Accepted: 07/31/2019] [Indexed: 01/22/2023] Open
Abstract
Cancer is a disease mainly caused by somatic genome alterations (SGAs) that perturb cellular signalling systems. Furthermore, the combination of pathway aberrations in a tumour defines its disease mechanism, and distinct disease mechanisms underlie the inter-tumour heterogeneity in terms of disease progression and responses to therapies. Discovering common disease mechanisms shared by tumours would provide guidance for precision oncology but remains a challenge. Here, we present a novel computational framework for revealing distinct combinations of aberrant signalling pathways in tumours. Specifically, we applied the tumour-specific causal inference algorithm (TCI) to identify causal relationships between SGAs and differentially expressed genes (DEGs) within tumours from the Cancer Genome Atlas (TCGA) study. Based on these causal inferences, we adopted a network-based method to identify modules of DEGs, such that the member DEGs within a module tend to be co-regulated by a common pathway. Using the expression status of genes in a module as a surrogate measure of the activation status of the corresponding pathways, we divided breast cancers (BRCAs) into five subgroups and glioblastoma multiformes (GBMs) into six subgroups with distinct combinations of pathway aberrations. The patient groups exhibited significantly different survival patterns, indicating that our approach can identify clinically relevant disease subtypes.
Collapse
Affiliation(s)
- Yifan Xue
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States
| | - Gregory Cooper
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States
| | - Chunhui Cai
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States
| | - Songjian Lu
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States
| | - Baoli Hu
- Department of Neurological Surgery, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States.,Paediatric Neurosurgery, UPMC Children's Hospital of Pittsburgh, Pittsburgh, 15213, United States.,Molecular and Cellular Cancer Biology Program, UPMC Hillman Cancer Centre, Pittsburgh, 15232, United States
| | - Xiaojun Ma
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States
| | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh School of Medicine, Pittsburgh, 15260, United States.
| |
Collapse
|