ABSTRACT
Human cytomegalovirus (HCMV) is a leading infectious cause of birth defects and the most common opportunistic infection that causes life-threatening diseases post-transplantation; however, an effective vaccine remains elusive. V160 is a live-attenuated replication defective HCMV vaccine that showed a 42.4% efficacy against primary HCMV infection among seronegative women in a phase 2b clinical trial. Here, we integrated the multicolor flow cytometry, longitudinal T cell receptor (TCR) sequencing, and single-cell RNA/TCR sequencing approaches to characterize the magnitude, phenotype, and functional quality of human T cell responses to V160. We demonstrated that V160 de novo induces IE-1 and pp65 specific durable polyfunctional effector CD8 T cells that are comparable to those induced by natural HCMV infection. We identified a variety of V160-responsive T cell clones which exhibit distinctive "transient" and "durable" expansion kinetics, and revealed a transcriptional signature that marks durable CD8 T cells post-vaccination. Our study enhances the understanding of human T-cell immune responses to V160 vaccination.
ABSTRACT
OBJECTIVE: The accuracy of deep learning models for many disease prediction problems is affected by time-varying covariates, rare incidence, covariate imbalance and delayed diagnosis when using structured electronic health records data. The situation is further exasperated when predicting the risk of one disease on condition of another disease, such as the hepatocellular carcinoma risk among patients with nonalcoholic fatty liver disease due to slow, chronic progression, the scarce of data with both disease conditions and the sex bias of the diseases. The goal of this study is to investigate the extent to which the aforementioned issues influence deep learning performance, and then devised strategies to tackle these challenges. These strategies were applied to improve hepatocellular carcinoma risk prediction among patients with nonalcoholic fatty liver disease. METHODS: We evaluated two representative deep learning models in the task of predicting the occurrence of hepatocellular carcinoma in a cohort of patients with nonalcoholic fatty liver disease (n = 220,838) from a national EHR database. The disease prediction task was carefully formulated as a classification problem while taking censorship and the length of follow-up into consideration. RESULTS: We developed a novel backward masking scheme to deal with the issue of delayed diagnosis which is very common in EHR data analysis and evaluate how the length of longitudinal information after the index date affects disease prediction. We observed that modeling time-varying covariates improved the performance of the algorithms and transfer learning mitigated reduced performance caused by the lack of data. In addition, covariate imbalance, such as sex bias in data impaired performance. Deep learning models trained on one sex and evaluated in the other sex showed reduced performance, indicating the importance of assessing covariate imbalance while preparing data for model training. CONCLUSIONS: The strategies developed in this work can significantly improve the performance of hepatocellular carcinoma risk prediction among patients with nonalcoholic fatty liver disease. Furthermore, our novel strategies can be generalized to apply to other disease risk predictions using structured electronic health records, especially for disease risks on condition of another disease.
Subject(s)
Carcinoma, Hepatocellular , Deep Learning , Liver Neoplasms , Non-alcoholic Fatty Liver Disease , Humans , Carcinoma, Hepatocellular/diagnosis , Carcinoma, Hepatocellular/epidemiology , Non-alcoholic Fatty Liver Disease/complications , Non-alcoholic Fatty Liver Disease/diagnosis , Non-alcoholic Fatty Liver Disease/epidemiology , Liver Neoplasms/diagnosis , Liver Neoplasms/epidemiology , Electronic Health RecordsABSTRACT
SMARCB1 loss has long been observed in many solid tumors. However, there is a need to elucidate targetable pathways driving growth and metastasis in SMARCB1-deficient tumors. Here, we demonstrate that SMARCB1 deficiency, defined as genomic SMARCB1 copy number loss associated with reduced mRNA, drives disease progression in patients with bladder cancer by engaging STAT3. SMARCB1 loss increases the chromatin accessibility of the STAT3 locus in vitro. Orthotopically implanted SMARCB1 knockout (KO) cell lines exhibit increased tumor growth and metastasis. SMARCB1-deficient tumors show an increased IL6/JAK/STAT3 signaling axis in in vivo models and patients. Furthermore, a pSTAT3 selective inhibitor, TTI-101, reduces tumor growth in SMARCB1 KO orthotopic cell line-derived xenografts and a SMARCB1-deficient patient derived xenograft model. We have identified a gene signature generated from SMARCB1 KO tumors that predicts SMARCB1 deficiency in patients. Overall, these findings support the clinical evaluation of STAT3 inhibitors for the treatment of SMARCB1-deficient bladder cancer.
Subject(s)
Interleukin-6 , Urinary Bladder Neoplasms , Humans , Interleukin-6/genetics , Interleukin-6/metabolism , Signal Transduction/genetics , SMARCB1 Protein/genetics , SMARCB1 Protein/metabolism , Urinary Bladder Neoplasms/genetics , Cell Line, Tumor , STAT3 Transcription Factor/genetics , STAT3 Transcription Factor/metabolismABSTRACT
Background: Deep learning models showed great success and potential when applied to many biomedical problems. However, the accuracy of deep learning models for many disease prediction problems is affected by time-varying covariates, rare incidence, and covariate imbalance when using structured electronic health records data. The situation is further exasperated when predicting the risk of one disease on condition of another disease, such as the hepatocellular carcinoma risk among patients with nonalcoholic fatty liver disease due to slow, chronic progression, the scarce of data with both disease conditions and the sex bias of the diseases. Objective: The goal of this study is to investigate the extent to which time-varying covariates, rare incidence, and covariate imbalance influence deep learning performance, and then devised strategies to tackle these challenges. These strategies were applied to improve hepatocellular carcinoma risk prediction among patients with nonalcoholic fatty liver disease. Methods: We evaluated two representative deep learning models in the task of predicting the occurrence of hepatocellular carcinoma in a cohort of patients with nonalcoholic fatty liver disease (n = 220,838) from a national EHR database. The disease prediction task was carefully formulated as a classification problem while taking censorship and the length of follow-up into consideration. Results: We developed a novel backward masking scheme to evaluate how the length of longitudinal information after the index date affects disease prediction. We observed that modeling time-varying covariates improved the performance of the algorithms and transfer learning mitigated reduced performance caused by the lack of data. In addition, covariate imbalance, such as sex bias in data impaired performance. Deep learning models trained on one sex and evaluated in the other sex showed reduced performance, indicating the importance of assessing covariate imbalance while preparing data for model training. Conclusions: Devising proper strategies to address challenges from time-varying covariates, lack of data, and covariate imbalance can be key to counteracting data bias and accurately predicting disease occurrence using deep learning models. The novel strategies developed in this work can significantly improve the performance of hepatocellular carcinoma risk prediction among patients with nonalcoholic fatty liver disease. Furthermore, our novel strategies can be generalized to apply to other disease risk predictions using structured electronic health records, especially for disease risks on condition of another disease.
ABSTRACT
Background: The integration of next-generation sequencing (NGS) comprehensive gene profiling (CGP) into clinical practice is playing an increasingly important role in oncology. Therefore, the HKU-HKSH Multi-disciplinary Molecular Tumour Board (MTB) was established to advance precision oncology in Hong Kong. A multicenter retrospective study investigated the feasibility of the HKU-HKSH MTB in determining genome-guided therapy for treatment-refractory solid cancers in Hong Kong. Methods: Patients who were presented at the HKU-HKSH MTB between August 2018 and June 2022 were included in this study. The primary study endpoints were the proportion of patients who receive MTB-guided therapy based on genomic analysis and overall survival (OS). Secondary endpoints included the proportion of patients with actionable genomic alterations, objective response rate (ORR), and disease control rate (DCR). The Kaplan-Meier method was used in the survival analyses, and hazard ratios were calculated using univariate Cox regression. Findings: 122 patients were reviewed at the HKU-HKSH MTB, and 63% (n = 77) adopted treatment per the MTB recommendations. These patients achieved a significantly longer median OS than those who did not receive MTB-guided therapy (12.7 months vs. 5.2 months, P = 0.0073). Their ORR and DCR were 29% and 65%, respectively. Interpretation: Our study demonstrated that among patients with heavily pre-treated advanced solid cancers, MTB-guided treatment could positively impact survival outcomes, thus illustrating the applicability of NGS CGPs in real-world clinical practice. Funding: The study was supported by the Li Shu Pui Medical Foundation. Dr Aya El Helali was supported by the Li Shu Pui Medical Foundation Fellowship grant from the Li Shu Pui Medical Foundation. Funders had no role in study design, data collection, data analysis, interpretation, or writing of the report.
ABSTRACT
Background: Tumour mutational burden (TMB) has emerged as a predictive marker for responsiveness to immune checkpoint inhibitors (ICI) in multiple tumour types. It can be calculated from somatic mutations detected from whole exome or targeted panel sequencing data. As mutations are unevenly distributed across the cancer genome, the clinical implications from TMB calculated using different genomic regions are not clear. Methods: Pan-cancer data of 10,179 samples were collected from The Cancer Genome Atlas cohort and 6,831 cancer patients with either ICI or non-ICI treatment outcomes were derived from published papers. TMB was calculated as the count of non-synonymous mutations and normalised by the size of genomic regions. Dirichlet method, linear regression and Poisson calibration models are used to unify TMB from different gene panels. Results: We found that panels based on cancer genes usually overestimate TMB compared to whole exome, potentially leading to misclassification of patients to receive ICI. The overestimation is caused by positive selection for mutations in cancer genes and cannot be completely addressed by the removal of mutational hotspots. We compared different approaches to address this discrepancy and developed a generalised statistical model capable of interconverting TMB derived from whole exome and different panel sequencing data, enabling TMB correction for patient stratification for ICI treatment. We show that in a cohort of lung cancer patients treated with ICI, when using a TMB cutoff of 10 mut/Mb, our corrected TMB outperforms the original panel-based TMB. Conclusion: Cancer gene-based panels usually overestimate TMB, and these findings will be valuable for unifying TMB calculations across cancer gene panels in clinical practice.
ABSTRACT
Multicellularity was a watershed development in evolution. However, it also meant that individual cells could escape regulatory mechanisms that restrict proliferation at a severe cost to the organism: cancer. From the standpoint of cellular organization, evolutionary complexity scales to organize different molecules within the intracellular milieu. The recent realization that many biomolecules can "phase-separate" into membraneless organelles, reorganizing cellular biochemistry in space and time, has led to an explosion of research activity in this area. In this review, we explore mechanistic connections between phase separation and cancer-associated processes and emerging examples of how these become deranged in malignancy. SIGNIFICANCE: One of the fundamental functions of phase separation is to rapidly and dynamically respond to environmental perturbations. Importantly, these changes often lead to alterations in cancer-relevant pathways and processes. This review covers recent advances in the field, including emerging principles and mechanisms of phase separation in cancer.
Subject(s)
Neoplasms , Organelles , Humans , Neoplasms/metabolism , Organelles/metabolism , ResearchABSTRACT
Acquired resistance and clonal heterogeneity are critical challenges in cancer treatment, and the lack of effective computational tools hampers the discovery of new treatments to overcome resistance. Using high-throughput transcriptomic databases of compound perturbation profiles, we have developed a bioinformatic strategy for identifying candidate drugs to overcome resistance with combinatorial therapy. We devised this strategy during an investigation into the acquired resistance against PARP inhibitors (PARPi) in a triple-negative inflammatory breast cancer cell line. In this study, we derived multiple PARPi-resistant clones and characterized their transcriptomic adaptations compared to the parental clone. The transcriptomes of the resistant clones showed substantial heterogeneity, highlighting the importance of characterizing multiple clones from the same tumour. Surprisingly, we found that these transcriptomic changes may not actually confer PARPi resistance, but they may nevertheless induce a shared secondary vulnerability. By modeling our data in relation to transcriptomic perturbation profiles of compounds, we uncovered deficiencies in Ras signaling that resulted from transcriptional adaptation to long-term PARPi treatment across multiple resistant clones. Due to these induced deficiencies, we predicted that the resistant clones would be sensitive to pharmacological reinforcement of PARPi-induced transcriptional adaptation. We then experimentally validated this predicted vulnerability that is shared by multiple resistant clones. Our results thus provide a promising paradigm for integrating transcriptomic data with compound perturbation profiles in order to identify drugs that can exploit an induced vulnerability and overcome therapeutic resistance, thus providing another strategy towards precision oncology.
ABSTRACT
MOTIVATION: Cross-sectional analyses of primary cancer genomes have identified regions of recurrent somatic copy-number alteration, many of which result from positive selection during cancer formation and contain driver genes. However, no effective approach exists for identifying genomic loci under significantly different degrees of selection in cancers of different subtypes, anatomic sites or disease stages. RESULTS: CNGPLD is a new tool for performing case-control somatic copy-number analysis that facilitates the discovery of differentially amplified or deleted copy-number aberrations in a case group of cancer compared with a control group of cancer. This tool uses a Gaussian process statistical framework in order to account for the covariance structure of copy-number data along genomic coordinates and to control the false discovery rate at the region level. AVAILABILITY AND IMPLEMENTATION: CNGPLD is freely available at https://bitbucket.org/djhshih/cngpld as an R package. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Genome , Neoplasms , Humans , Cross-Sectional Studies , Genomics , DNA Copy Number Variations , Neoplasms/genetics , Case-Control Studies , SoftwareABSTRACT
Nucleotide excision repair (NER) resolves DNA adducts, such as those caused by ultraviolet light. Deficient NER (dNER) results in a higher mutation rate that can predispose to cancer development and premature ageing phenotypes. Here, we used isogenic dNER model cell lines to establish a gene expression signature that can accurately predict functional NER capacity in both cell lines and patient samples. Critically, none of the identified NER deficient cell lines harbored mutations in any NER genes, suggesting that the prevalence of NER defects may currently be underestimated. Identification of compounds that induce the dNER gene expression signature led to the discovery that NER can be functionally impaired by GSK3 inhibition, leading to synergy when combined with cisplatin treatment. Furthermore, we predicted and validated multiple novel drugs that are synthetically lethal with NER defects using the dNER gene signature as a drug discovery platform. Taken together, our work provides a dynamic predictor of NER function that may be applied for therapeutic stratification as well as development of novel biological insights in human tumors.
Subject(s)
DNA Repair/genetics , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Neoplasms/drug therapy , Neoplasms/genetics , Cell Line, Tumor , Humans , Reproducibility of ResultsABSTRACT
Defect in DNA damage response (DDR) is a common feature of cancer cells, which regulates tumor growth and therapeutic response. Recently, the approval of immune checkpoint blockade (ICB) for tumors with defective mismatch repair has paved the way for investigating the role of other DDR defects in sensitizing cancer to ICB therapy. Despite great progress in understanding DDR pathways, the mechanisms that link DDR defects and ICB response remain incompletely understood. Further, the clinical activity of ICB in patients with DDR defective tumors has not been well described. Here, we discuss recent studies demonstrating that biomarkers in DDR pathways may serve as potential predictors to guide the selection of patients for ICB therapy. A better understanding of the relationship between deficiency in DDR and response to ICB would facilitate efforts in optimizing the efficacy of immunotherapy.
ABSTRACT
Primary cutaneous γδ T cell lymphomas (PCGDTLs) represent a heterogeneous group of uncommon but aggressive cancers. Herein, we perform genome-wide DNA, RNA, and T cell receptor (TCR) sequencing on 29 cutaneous γδ lymphomas. We find that PCGDTLs are not uniformly derived from Vδ2 cells. Instead, the cell-of-origin depends on the tissue compartment from which the lymphomas are derived. Lymphomas arising from the outer layer of skin are derived from Vδ1 cells, the predominant γδ cell in the epidermis and dermis. In contrast, panniculitic lymphomas arise from Vδ2 cells, the predominant γδ T cell in the fat. We also show that TCR chain usage is non-random, suggesting common antigens for Vδ1 and Vδ2 lymphomas respectively. In addition, Vδ1 and Vδ2 PCGDTLs harbor similar genomic landscapes with potentially targetable oncogenic mutations in the JAK/STAT, MAPK, MYC, and chromatin modification pathways. Collectively, these findings suggest a paradigm for classifying, staging, and treating these diseases.
Subject(s)
Lymphoma, T-Cell, Cutaneous/genetics , Lymphoma, T-Cell, Cutaneous/pathology , Receptors, Antigen, T-Cell, gamma-delta/metabolism , Skin Neoplasms/genetics , Skin Neoplasms/pathology , Amino Acid Sequence , Antigens, CD1d/metabolism , Chromatin Assembly and Disassembly , Epitopes/immunology , Genome, Human , HEK293 Cells , Humans , Lymph Nodes/pathology , Models, Biological , Mutation/genetics , Phenotype , Principal Component Analysis , Signal Transduction , Skin/pathology , Transcription, Genetic , Transcriptome/geneticsABSTRACT
Brain metastases from lung adenocarcinoma (BM-LUAD) frequently cause patient mortality. To identify genomic alterations that promote brain metastases, we performed whole-exome sequencing of 73 BM-LUAD cases. Using case-control analyses, we discovered candidate drivers of brain metastasis by identifying genes with more frequent copy-number aberrations in BM-LUAD compared to 503 primary LUADs. We identified three regions with significantly higher amplification frequencies in BM-LUAD, including MYC (12 versus 6%), YAP1 (7 versus 0.8%) and MMP13 (10 versus 0.6%), and significantly more frequent deletions in CDKN2A/B (27 versus 13%). We confirmed that the amplification frequencies of MYC, YAP1 and MMP13 were elevated in an independent cohort of 105 patients with BM-LUAD. Functional assessment in patient-derived xenograft mouse models validated the notion that MYC, YAP1 or MMP13 overexpression increased the incidence of brain metastasis. These results demonstrate that somatic alterations contribute to brain metastases and that genomic sequencing of a sufficient number of metastatic tumors can reveal previously unknown metastatic drivers.
Subject(s)
Adenocarcinoma of Lung/genetics , Brain Neoplasms/genetics , Lung Neoplasms/genetics , Neoplasm Metastasis/genetics , Adenocarcinoma of Lung/pathology , Animals , Brain Neoplasms/pathology , Case-Control Studies , Cell Line , DNA Copy Number Variations/genetics , Female , Genes, myc/genetics , Genomics/methods , HEK293 Cells , Humans , Lung Neoplasms/pathology , Male , Matrix Metalloproteinase 13/genetics , Mice , Mice, Nude , Mutation/genetics , Neoplasm Metastasis/pathology , Transcription Factors/genetics , Exome SequencingABSTRACT
Deficient DNA mismatch repair (dMMR) induces a hypermutator phenotype that can lead to tumorigenesis; however, the functional impact of the high mutation burden resulting from this phenotype remains poorly explored. Here, we demonstrate that dMMR-induced destabilizing mutations lead to proteome instability in dMMR tumors, resulting in an abundance of misfolded protein aggregates. To compensate, dMMR cells utilize a Nedd8-mediated degradation pathway to facilitate clearance of misfolded proteins. Blockade of this Nedd8 clearance pathway with MLN4924 causes accumulation of misfolded protein aggregates, ultimately inducing immunogenic cell death in dMMR cancer cells. To leverage this immunogenic cell death, we combined MLN4924 treatment with PD1 inhibition and found the combination was synergistic, significantly improving efficacy over either treatment alone.
Subject(s)
Cyclopentanes/pharmacology , DNA Mismatch Repair , Endometrial Neoplasms/drug therapy , Proteome/genetics , Pyrimidines/pharmacology , Animals , Cell Line, Tumor , Colorectal Neoplasms/drug therapy , Colorectal Neoplasms/genetics , Colorectal Neoplasms/immunology , Endometrial Neoplasms/genetics , Endometrial Neoplasms/immunology , Female , HCT116 Cells , Humans , Immunotherapy/methods , Mice, Inbred C57BL , Mice, Transgenic , Microsatellite Instability , Mutation , NEDD8 Protein/antagonists & inhibitors , NEDD8 Protein/metabolism , Programmed Cell Death 1 Receptor/antagonists & inhibitors , Programmed Cell Death 1 Receptor/metabolism , Protein Stability , Xenograft Model Antitumor AssaysABSTRACT
The chromatin remodeling factor chromodomain helicase DNA-binding protein 4 (CHD4) is a core component of the nucleosome remodeling and deacetylase (NuRD) complex. Due to its important role in DNA damage repair, CHD4 has been identified as a key determinant in cancer progression, stem cell differentiation, and T cell and B cell development. Accumulating evidence has revealed that CHD4 can function in NuRD dependent and independent manner in response to DNA damage. Mutations of CHD4 have been shown to diminish its functions, which indicates that interpretation of its mutations may provide tangible benefit for patients. The expression of CHD4 play a dual role in sensitizing cancer cells to chemotherapeutic agents, which provides new insights into the contribution of CHD4 to tumor biology and new therapeutic avenues.
ABSTRACT
Posterior fossa ependymoma comprise three distinct molecular variants, termed PF-EPN-A (PFA), PF-EPN-B (PFB), and PF-EPN-SE (subependymoma). Clinically, they are very disparate and PFB tumors are currently being considered for a trial of radiation avoidance. However, to move forward, unraveling the heterogeneity within PFB would be highly desirable. To discern the molecular heterogeneity within PFB, we performed an integrated analysis consisting of DNA methylation profiling, copy-number profiling, gene expression profiling, and clinical correlation across a cohort of 212 primary posterior fossa PFB tumors. Unsupervised spectral clustering and t-SNE analysis of genome-wide methylation data revealed five distinct subtypes of PFB tumors, termed PFB1-5, with distinct demographics, copy-number alterations, and gene expression profiles. All PFB subtypes were distinct from PFA and posterior fossa subependymomas. Of the five subtypes, PFB4 and PFB5 are more discrete, consisting of younger and older patients, respectively, with a strong female-gender enrichment in PFB5 (age: p = 0.011, gender: p = 0.04). Broad copy-number aberrations were common; however, many events such as chromosome 2 loss, 5 gain, and 17 loss were enriched in specific subtypes and 1q gain was enriched in PFB1. Late relapses were common across all five subtypes, but deaths were uncommon and present in only two subtypes (PFB1 and PFB3). Unlike the case in PFA ependymoma, 1q gain was not a robust marker of poor progression-free survival; however, chromosome 13q loss may represent a novel marker for risk stratification across the spectrum of PFB subtypes. Similar to PFA ependymoma, there exists a significant intertumoral heterogeneity within PFB, with distinct molecular subtypes identified. Even when accounting for this heterogeneity, extent of resection remains the strongest predictor of poor outcome. However, this biological heterogeneity must be accounted for in future preclinical modeling and personalized therapies.
Subject(s)
DNA Copy Number Variations/genetics , Ependymoma/classification , Ependymoma/genetics , Infratentorial Neoplasms/classification , Infratentorial Neoplasms/genetics , Adolescent , Adult , Age Factors , Child , Cohort Studies , DNA Methylation/genetics , Ependymoma/pathology , Ependymoma/surgery , Female , Gene Expression Profiling , Humans , Infratentorial Neoplasms/pathology , Infratentorial Neoplasms/surgery , Kaplan-Meier Estimate , Male , Microarray Analysis , Middle Aged , Young AdultABSTRACT
While the preponderance of morbidity and mortality in medulloblastoma patients are due to metastatic disease, most research focuses on the primary tumor due to a dearth of metastatic tissue samples and model systems. Medulloblastoma metastases are found almost exclusively on the leptomeningeal surface of the brain and spinal cord; dissemination is therefore thought to occur through shedding of primary tumor cells into the cerebrospinal fluid followed by distal re-implantation on the leptomeninges. We present evidence for medulloblastoma circulating tumor cells (CTCs) in therapy-naive patients and demonstrate in vivo, through flank xenografting and parabiosis, that medulloblastoma CTCs can spread through the blood to the leptomeningeal space to form leptomeningeal metastases. Medulloblastoma leptomeningeal metastases express high levels of the chemokine CCL2, and expression of CCL2 in medulloblastoma in vivo is sufficient to drive leptomeningeal dissemination. Hematogenous dissemination of medulloblastoma offers a new opportunity to diagnose and treat lethal disseminated medulloblastoma.
Subject(s)
Medulloblastoma/blood supply , Medulloblastoma/pathology , Meningeal Neoplasms/blood supply , Meningeal Neoplasms/secondary , Allografts , Animals , Cell Line, Tumor , Chemokine CCL2/metabolism , Chromosomes, Human, Pair 10/genetics , Female , Humans , Male , Medulloblastoma/genetics , Mice, SCID , Neoplastic Cells, Circulating , ParabiosisABSTRACT
While molecular subgrouping has revolutionized medulloblastoma classification, the extent of heterogeneity within subgroups is unknown. Similarity network fusion (SNF) applied to genome-wide DNA methylation and gene expression data across 763 primary samples identifies very homogeneous clusters of patients, supporting the presence of medulloblastoma subtypes. After integration of somatic copy-number alterations, and clinical features specific to each cluster, we identify 12 different subtypes of medulloblastoma. Integrative analysis using SNF further delineates group 3 from group 4 medulloblastoma, which is not as readily apparent through analyses of individual data types. Two clear subtypes of infants with Sonic Hedgehog medulloblastoma with disparate outcomes and biology are identified. Medulloblastoma subtypes identified through integrative clustering have important implications for stratification of future clinical trials.