A Hematogenous Route for Medulloblastoma Leptomeningeal Metastases.

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.

Developing deep learning-based strategies to predict the risk of hepatocellular carcinoma among patients with nonalcoholic fatty liver disease from electronic health records.

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.

CNGPLD: case-control copy-number analysis using Gaussian process latent difference.

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.

Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

A Gene Expression Signature to Predict Nucleotide Excision Repair Defects and Novel Therapeutic Approaches.

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.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma.

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

Heterogeneity within the PF-EPN-B ependymoma subgroup.

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.

Clonal selection drives genetic divergence of metastatic medulloblastoma.

Medulloblastoma, the most common malignant paediatric brain tumour, arises in the cerebellum and disseminates through the cerebrospinal fluid in the leptomeningeal space to coat the brain and spinal cord. Dissemination, a marker of poor prognosis, is found in up to 40% of children at diagnosis and in most children at the time of recurrence. Affected children therefore are treated with radiation to the entire developing brain and spinal cord, followed by high-dose chemotherapy, with the ensuing deleterious effects on the developing nervous system. The mechanisms of dissemination through the cerebrospinal fluid are poorly studied, and medulloblastoma metastases have been assumed to be biologically similar to the primary tumour. Here we show that in both mouse and human medulloblastoma, the metastases from an individual are extremely similar to each other but are divergent from the matched primary tumour. Clonal genetic events in the metastases can be demonstrated in a restricted subclone of the primary tumour, suggesting that only rare cells within the primary tumour have the ability to metastasize. Failure to account for the bicompartmental nature of metastatic medulloblastoma could be a major barrier to the development of effective targeted therapies.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes.

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-ß signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

Transcriptional signature of durable effector T cells elicited by a replication defective HCMV vaccine.

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.

The IL6/JAK/STAT3 signaling axis is a therapeutic vulnerability in SMARCB1-deficient bladder cancer.

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.

TERT promoter mutations are highly recurrent in SHH subgroup medulloblastoma.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

Developing deep learning-based strategies to predict the risk of hepatocellular carcinoma among patients with nonalcoholic fatty liver disease from electronic health records.

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.

Tumour mutational burden is overestimated by target cancer gene panels.

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.

The impact of the multi-disciplinary molecular tumour board and integrative next generation sequencing on clinical outcomes in advanced solid tumours.

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.

Rapid, reliable, and reproducible molecular sub-grouping of clinical medulloblastoma samples.

The diagnosis of medulloblastoma likely encompasses several distinct entities, with recent evidence for the existence of at least four unique molecular subgroups that exhibit distinct genetic, transcriptional, demographic, and clinical features. Assignment of molecular subgroup through routine profiling of high-quality RNA on expression microarrays is likely impractical in the clinical setting. The planning and execution of medulloblastoma clinical trials that stratify by subgroup, or which are targeted to a specific subgroup requires technologies that can be economically, rapidly, reliably, and reproducibly applied to formalin-fixed paraffin embedded (FFPE) specimens. In the current study, we have developed an assay that accurately measures the expression level of 22 medulloblastoma subgroup-specific signature genes (CodeSet) using nanoString nCounter Technology. Comparison of the nanoString assay with Affymetrix expression array data on a training series of 101 medulloblastomas of known subgroup demonstrated a high concordance (Pearson correlation r = 0.86). The assay was validated on a second set of 130 non-overlapping medulloblastomas of known subgroup, correctly assigning 98% (127/130) of tumors to the appropriate subgroup. Reproducibility was demonstrated by repeating the assay in three independent laboratories in Canada, the United States, and Switzerland. Finally, the nanoString assay could confidently predict subgroup in 88% of recent FFPE cases, of which 100% had accurate subgroup assignment. We present an assay based on nanoString technology that is capable of rapidly, reliably, and reproducibly assigning clinical FFPE medulloblastoma samples to their molecular subgroup, and which is highly suited for future medulloblastoma clinical trials.

Biomolecular Condensation: A New Phase in Cancer Research.

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.

Exploiting induced vulnerability to overcome PARPi resistance and clonal heterogeneity in BRCA mutant triple-negative inflammatory breast cancer.

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.

Role of DNA repair defects in predicting immunotherapy response.

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.