Molecular and clinical heterogeneity within MYC-family amplified medulloblastoma is associated with survival outcomes: A multicenter cohort study.

BACKGROUND: MYC/MYCN are the most frequent oncogene amplifications in medulloblastoma (MB) and its primary biomarkers of high-risk (HR) disease. However, while many patients' MYC(N)-amplified tumors are treatment-refractory, some achieve long-term survival. We therefore investigated clinicobiological heterogeneity within MYC(N)-amplified MB and determined its relevance for improved disease management. METHODS: We characterized the clinical and molecular correlates of MYC- (MYC-MB; n = 64) and MYCN-amplified MBs (MYCN-MB; n = 95), drawn from >1600 diagnostic cases. RESULTS: Most MYC-MBs were molecular group 3 (46/58; 79% assessable) and aged ≥3 years at diagnosis (44/64 [69%]). We identified a "canonical" very high-risk (VHR) MYC-amplified group (n = 51/62; 82%) with dismal survival irrespective of treatment (11% 5-year progression-free survival [PFS]), defined by co-occurrence with ≥1 additional established risk factor(s) (subtotal surgical-resection [STR], metastatic disease, LCA pathology), and commonly group 3/4 subgroup 2 with a high proportion of amplified cells. The majority of remaining noncanonical MYC-MBs survived (i.e. non-group 3/group 3 without other risk features; 11/62 (18%); 61% 5-year PFS). MYCN survival was primarily related to molecular group; MYCN-amplified SHH MB, and group 3/4 MB with additional risk factors, respectively defined VHR and HR groups (VHR, 39% [35/89]; 20% 5-year PFS/HR, 33% [29/89]; 46% 5-year PFS). Twenty-two out of 35 assessable MYCN-amplified SHH tumors harbored TP53 mutations; 9/12 (75%) with data were germline. MYCN-amplified group 3/4 MB with no other risk factors (28%; 25/89) had 70% 5-year PFS. CONCLUSIONS: MYC(N)-amplified MB displays significant clinicobiological heterogeneity. Diagnostics incorporating molecular groups, subgroups, and clinical factors enable their risk assessment. VHR "canonical" MYC tumors are essentially incurable and SHH-MYCN-amplified MBs fare extremely poorly (20% survival at 5 years); both require urgent development of alternative treatment strategies. Conventional risk-adapted therapies are appropriate for more responsive groups, such as noncanonical MYC and non-SHH-MYCN MB.

Poglut2/3 double knockout in mice results in neonatal lethality with reduced levels of fibrillin in lung tissues.

Fibrillin microfibrils play a critical role in the formation of elastic fibers, tissue/organ development, and cardiopulmonary function. These microfibrils not only provide structural support and flexibility to tissues, but they also regulate growth factor signaling through a plethora of microfibril-binding proteins in the extracellular space. Mutations in fibrillins are associated with human diseases affecting cardiovascular, pulmonary, skeletal, and ocular systems. Fibrillins consist of up to 47 epidermal growth factor-like repeats, of which more than half are modified by protein O-glucosyltransferase 2 (POGLUT2) and/or POGLUT3. Loss of these modifications reduces secretion of N-terminal fibrillin constructs overexpressed in vitro. Here, we investigated the role of POGLUT2 and POGLUT3 in vivo using a Poglut2/3 double knockout (DKO) mouse model. Blocking O-glucosylation caused neonatal death with skeletal, pulmonary, and eye defects reminiscent of fibrillin/elastin mutations. Proteomic analyses of DKO dermal fibroblast medium and extracellular matrix provided evidence that fibrillins were more sensitive to loss of O-glucose compared to other POGLUT2/3 substrates. This conclusion was supported by immunofluorescent analyses of late gestation DKO lungs where FBN levels were reduced and microfibrils appeared fragmented in the pulmonary arteries and veins, bronchioles, and developing saccules. Defects in fibrillin microfibrils likely contributed to impaired elastic fiber formation and histological changes observed in DKO lung blood vessels, bronchioles, and saccules. Collectively, these results highlight the importance of POGLUT2/3-mediated O-glucosylation in vivo and open the possibility that O-glucose modifications on fibrillin influence microfibril assembly and or protein interactions in the ECM environment.

MRI-derived brain iron, grey matter volume, and risk of dementia and Parkinson's disease: Observational and genetic analysis in the UK Biobank cohort.

BACKGROUND: Iron overload is observed in neurodegenerative diseases, especially Alzheimer's disease (AD) and Parkinson's disease (PD). Homozygotes for the iron-overload (haemochromatosis) causing HFE p.C282Y variant have increased risk of dementia and PD. Whether brain iron deposition is causal or secondary to the neurodegenerative processes in the general population is unclear. METHODS: We analysed 39,533 UK Biobank participants of European genetic ancestry with brain MRI data. We studied brain iron estimated by R2* and quantitative susceptibility mapping (QSM) in 8 subcortical regions: accumbens, amygdala, caudate, hippocampus, pallidum, putamen, substantia nigra, and thalamus. We performed genome-wide associations studies (GWAS) and used Mendelian Randomization (MR) methods to estimate the causal effect of brain iron on grey matter volume, and risk of AD, non-AD and PD. We also used MR to test whether genetic liability to AD or PD causally increased brain iron (R2* and QSM). FINDINGS: In GWAS of R2* and QSM we replicated 83% of previously reported genetic loci and identified 174 further loci across all eight brain regions. Higher genetically predicted brain iron, using both R2* and QSM, was associated with lower grey matter volumes in the caudate, putamen and thalamus (e.g., Beta-putamenQSM: -0.37, p = 2*10-46). Higher genetically predicted thalamus R2* was associated with increased risk of non-AD dementia (OR 1.36(1.16;1.60), p = 2*10-4) but not AD (p > 0.05). In males, genetically predicted putamen R2* increased non-AD dementia risk, but not in females. Higher genetically predicted iron in the caudate, putamen, and substantia nigra was associated with an increased risk of PD (Odds Ratio QSM âˆ¼ substantia-nigra 1.21(1.07;1.37), p = 0.003). Genetic liability to AD or PD was not associated with R2* or QSM in the dementia or PD-associated regions. INTERPRETATION: Our genetic analysis supports a causal effect of higher iron deposition in specific subcortical brain regions for Parkinson's disease, grey matter volume, and non-Alzheimer's dementia.

Sortase-Modified Cholera Toxoids Show Specific Golgi Localization.

Cholera toxoid is an established tool for use in cellular tracing in neuroscience and cell biology. We use a sortase labeling approach to generate site-specific N-terminally modified variants of both the A2-B5 heterohexamer and B5 pentamer forms of the toxoid. Both forms of the toxoid are endocytosed by GM1-positive mammalian cells, and while the heterohexameric toxoid was principally localized in the ER, the B5 pentamer showed an unexpectedly specific localization in the medial/trans-Golgi. This study suggests a future role for specifically labeled cholera toxoids in live-cell imaging beyond their current applications in neuronal tracing and labeling of lipid rafts in fixed cells.

Metabolite profiles of medulloblastoma for rapid and non-invasive detection of molecular disease groups.

BACKGROUND: The malignant childhood brain tumour, medulloblastoma, is classified clinically into molecular groups which guide therapy. DNA-methylation profiling is the current classification 'gold-standard', typically delivered 3-4 weeks post-surgery. Pre-surgery non-invasive diagnostics thus offer significant potential to improve early diagnosis and clinical management. Here, we determine tumour metabolite profiles of the four medulloblastoma groups, assess their diagnostic utility using tumour tissue and potential for non-invasive diagnosis using in vivo magnetic resonance spectroscopy (MRS). METHODS: Metabolite profiles were acquired by high-resolution magic-angle spinning NMR spectroscopy (MAS) from 86 medulloblastomas (from 59 male and 27 female patients), previously classified by DNA-methylation array (WNT (n = 9), SHH (n = 22), Group3 (n = 21), Group4 (n = 34)); RNA-seq data was available for sixty. Unsupervised class-discovery was performed and a support vector machine (SVM) constructed to assess diagnostic performance. The SVM classifier was adapted to use only metabolites (n = 10) routinely quantified from in vivo MRS data, and re-tested. Glutamate was assessed as a predictor of overall survival. FINDINGS: Group-specific metabolite profiles were identified; tumours clustered with good concordance to their reference molecular group (93%). GABA was only detected in WNT, taurine was low in SHH and lipids were high in Group3. The tissue-based metabolite SVM classifier had a cross-validated accuracy of 89% (100% for WNT) and, adapted to use metabolites routinely quantified in vivo, gave a combined classification accuracy of 90% for SHH, Group3 and Group4. Glutamate predicted survival after incorporating known risk-factors (HR = 3.39, 95% CI 1.4-8.1, p = 0.025). INTERPRETATION: Tissue metabolite profiles characterise medulloblastoma molecular groups. Their combination with machine learning can aid rapid diagnosis from tissue and potentially in vivo. Specific metabolites provide important information; GABA identifying WNT and glutamate conferring poor prognosis. FUNDING: Children with Cancer UK, Cancer Research UK, Children's Cancer North and a Newcastle University PhD studentship.

Iron and risk of dementia: Mendelian randomisation analysis in UK Biobank.

BACKGROUND: Brain iron deposition is common in dementia, but whether serum iron is a causal risk factor is unknown. We aimed to determine whether genetic predisposition to higher serum iron status biomarkers increased risk of dementia and atrophy of grey matter. METHODS: We analysed UK Biobank participants clustered into European (N=451284), African (N=7477) and South Asian (N=9570) groups by genetic similarity to the 1000 genomes project. Using Mendelian randomisation methods, we estimated the association between genetically predicted serum iron (transferrin saturation [TSAT] and ferritin), grey matter volume and genetic liability to clinically defined dementia (including Alzheimer's disease [AD], non-AD dementia, and vascular dementia) from hospital and primary care records. We also performed time-to-event (competing risks) analysis of the TSAT polygenic score on risk of clinically defined non-AD dementia. RESULTS: In Europeans, higher genetically predicted TSAT increased genetic liability to dementia (Odds Ratio [OR]: 1.15, 95% Confidence Intervals [CI] 1.04 to 1.26, p=0.0051), non-AD dementia (OR: 1.27, 95% CI 1.12 to 1.45, p=0.00018) and vascular dementia (OR: 1.37, 95% CI 1.12 to 1.69, p=0.0023), but not AD (OR: 1.00, 95% CI 0.86 to 1.15, p=0.97). Higher TSAT was also associated with increased risk of non-AD dementia in participants of African, but not South Asian groups. In survival analysis using a TSAT polygenic score, the effect was independent of apolipoprotein-E ε4 genotype (with adjustment subdistribution Hazard Ratio: 1.74, 95% CI 1.33 to 2.28, p=0.00006). Genetically predicted TSAT was associated with lower grey matter volume in caudate, putamen and thalamus, and not in other areas of interest. DISCUSSION: Genetic evidence supports a causal relationship between higher TSAT and risk of clinically defined non-AD and vascular dementia, in European and African groups. This association appears to be independent of apolipoprotein-E ε4.

Quantitative N- or C-Terminal Labelling of Proteins with Unactivated Peptides by Use of Sortases and a d-Aminopeptidase.

Quantitative and selective labelling of proteins is widely used in both academic and industrial laboratories, and catalytic labelling of proteins using transpeptidases, such as sortases, has proved to be a popular strategy for such selective modification. A major challenge for this class of enzymes is that the majority of procedures require an excess of the labelling reagent or, alternatively, activated substrates rather than simple commercially sourced peptides. We report the use of a coupled enzyme strategy which enables quantitative N- and C-terminal labelling of proteins using unactivated labelling peptides. The use of an aminopeptidase in conjunction with a transpeptidase allows sequence-specific degradation of the peptide by-product, shifting the equilibrium to favor product formation, which greatly enhances the reaction efficiency. Subsequent optimisation of the reaction allows N-terminal labelling of proteins using essentially equimolar ratios of peptide label to protein and C-terminal labelling with only a small excess. Minimizing the amount of substrate required for quantitative labelling has the potential to improve industrial processes and facilitate the use of transpeptidation as a method for protein labelling.

Modelling how plant cell-cycle progression leads to cell size regulation.

Populations of cells typically maintain a consistent size, despite cell division rarely being precisely symmetrical. Therefore, cells must possess a mechanism of "size control", whereby the cell volume at birth affects cell-cycle progression. While size control mechanisms have been elucidated in a number of other organisms, it is not yet clear how this mechanism functions in plants. Here, we present a mathematical model of the key interactions in the plant cell cycle. Model simulations reveal that the network of interactions exhibits limit-cycle solutions, with biological switches underpinning both the G1/S and G2/M cell-cycle transitions. Embedding this network model within growing cells, we test hypotheses as to how cell-cycle progression can depend on cell size. We investigate two different mechanisms at both the G1/S and G2/M transitions: (i) differential expression of cell-cycle activator and inhibitor proteins (with synthesis of inhibitor proteins being independent of cell size), and (ii) equal inheritance of inhibitor proteins after cell division. The model demonstrates that both these mechanisms can lead to larger daughter cells progressing through the cell cycle more rapidly, and can thus contribute to cell-size control. To test how these features enable size homeostasis over multiple generations, we then simulated these mechanisms in a cell-population model with multiple rounds of cell division. These simulations suggested that integration of size-control mechanisms at both G1/S and G2/M provides long-term cell-size homeostasis. We concluded that while both size independence and equal inheritance of inhibitor proteins can reduce variations in cell size across individual cell-cycle phases, combining size-control mechanisms at both G1/S and G2/M is essential to maintain size homeostasis over multiple generations. Thus, our study reveals how features of the cell-cycle network enable cell-cycle progression to depend on cell size, and provides a mechanistic understanding of how plant cell populations maintain consistent size over generations.

Deriving a continuum score for group 3 and 4 medulloblastoma tumor samples analyzed via RNA-sequencing or DNA methylation microarray.

Here, we present a protocol for deriving a continuum score for group 3 and 4 medulloblastoma tumor samples analyzed via RNA-sequencing or DNA methylation microarray. We describe steps for utilizing NMF-defined group 3/group 4 metagenes to calculate a continuum score between 0 and 1 that can be projected onto new sample data analyzed via RNA-sequencing. We then detail procedures for reverse engineering a continuum score for samples analyzed via DNA methylation microarray using a random forest classifier.

Toward machine-assisted tuning avoiding the underestimation of uncertainty in climate change projections.

Documenting the uncertainty of climate change projections is a fundamental objective of the inter-comparison exercises organized to feed into the Intergovernmental Panel on Climate Change (IPCC) reports. Usually, each modeling center contributes to these exercises with one or two configurations of its climate model, corresponding to a particular choice of "free parameter" values, resulting from a long and often tedious "model tuning" phase. How much uncertainty is omitted by this selection and how might readers of IPCC reports and users of climate projections be misled by its omission? We show here how recent machine learning approaches can transform the way climate model tuning is approached, opening the way to a simultaneous acceleration of model improvement and parametric uncertainty quantification. We show how an automatic selection of model configurations defined by different values of free parameters can produce different "warming worlds," all consistent with present-day observations of the climate system.

The proteomic landscape of soft tissue sarcomas.

Soft tissue sarcomas (STS) are rare and diverse mesenchymal cancers with limited treatment options. Here we undertake comprehensive proteomic profiling of tumour specimens from 321 STS patients representing 11 histological subtypes. Within leiomyosarcomas, we identify three proteomic subtypes with distinct myogenesis and immune features, anatomical site distribution and survival outcomes. Characterisation of undifferentiated pleomorphic sarcomas and dedifferentiated liposarcomas with low infiltrating CD3 + T-lymphocyte levels nominates the complement cascade as a candidate immunotherapeutic target. Comparative analysis of proteomic and transcriptomic profiles highlights the proteomic-specific features for optimal risk stratification in angiosarcomas. Finally, we define functional signatures termed Sarcoma Proteomic Modules which transcend histological subtype classification and show that a vesicle transport protein signature is an independent prognostic factor for distant metastasis. Our study highlights the utility of proteomics for identifying molecular subgroups with implications for risk stratification and therapy selection and provides a rich resource for future sarcoma research.

Differences in Longevity and Temperature-Driven Extrinsic Incubation Period Correlate with Varying Dengue Risk in the Arizona-Sonora Desert Region.

Dengue transmission is determined by a complex set of interactions between the environment, Aedes aegypti mosquitoes, dengue viruses, and humans. Emergence in new geographic areas can be unpredictable, with some regions having established mosquito populations for decades without locally acquired transmission. Key factors such as mosquito longevity, temperature-driven extrinsic incubation period (EIP), and vector-human contact can strongly influence the potential for disease transmission. To assess how these factors interact at the edge of the geographical range of dengue virus transmission, we conducted mosquito sampling in multiple urban areas located throughout the Arizona-Sonora desert region during the summer rainy seasons from 2013 to 2015. Mosquito population age structure, reflecting mosquito survivorship, was measured using a combination of parity analysis and relative gene expression of an age-related gene, SCP-1. Bloodmeal analysis was conducted on field collected blood-fed mosquitoes. Site-specific temperature was used to estimate the EIP, and this predicted EIP combined with mosquito age were combined to estimate the abundance of "potential" vectors (i.e., mosquitoes old enough to survive the EIP). Comparisons were made across cities by month and year. The dengue endemic cities Hermosillo and Ciudad Obregon, both in the state of Sonora, Mexico, had higher abundance of potential vectors than non-endemic Nogales, Sonora, Mexico. Interestingly, Tucson, Arizona consistently had a higher estimated abundance of potential vectors than dengue endemic regions of Sonora, Mexico. There were no observed city-level differences in species composition of blood meals. Combined, these data offer insights into the critical factors required for dengue transmission at the ecological edge of the mosquito's range. However, further research is needed to integrate an understanding of how social and additional environmental factors constrain and enhance dengue transmission in emerging regions.

Timing is everything: A connection between medulloblastoma prognosis and foetal cerebellar development.

The childhood brain tumour medulloblastoma is typically classified into multiple discrete molecular subgroups with characteristic DNA methylation and expression patterns. Several of these subgroups are used as, or proposed to be, an effective basis for treatment stratification. Here, we highlight the close connection between the findings described in a recent series of studies which, together, strongly imply a continuous association between survival outcome, the transcriptional profile of a Group3/Group4 (i.e. non-WNT/non-SHH) medulloblastoma and the specific point during early foetal cerebellar development at which initial pathogenic disruption took place. This has important implications for future efforts to model the disease by incorporating driving molecular features into their specific developmental context. This further suggests that instead of relying upon discrete DNA methylation subgroups, using expression biomarkers as the basis of a continuous risk predictor may produce a more effective risk stratification of patients with Group3/Group4 medulloblastoma.

Molecular characterisation defines clinically-actionable heterogeneity within Group 4 medulloblastoma and improves disease risk-stratification.

Group 4 tumours (MBGrp4) represent the majority of non-WNT/non-SHH medulloblastomas. Their clinical course is poorly predicted by current risk-factors. MBGrp4 molecular substructures have been identified (e.g. subgroups/cytogenetics/mutations), however their inter-relationships and potential to improve clinical sub-classification and risk-stratification remain undefined. We comprehensively characterised the paediatric MBGrp4 molecular landscape and determined its utility to improve clinical management. A clinically-annotated discovery cohort (n = 362 MBGrp4) was assembled from UK-CCLG institutions and SIOP-UKCCSG-PNET3, HIT-SIOP-PNET4 and PNET HR + 5 clinical trials. Molecular profiling was undertaken, integrating driver mutations, second-generation non-WNT/non-SHH subgroups (1-8) and whole-chromosome aberrations (WCAs). Survival models were derived for patients ≥ 3 years of age who received contemporary multi-modal therapies (n = 323). We first independently derived and validated a favourable-risk WCA group (WCA-FR) characterised by ≥ 2 features from chromosome 7 gain, 8 loss, and 11 loss. Remaining patients were high-risk (WCA-HR). Subgroups 6 and 7 were enriched for WCA-FR (p < 0·0001) and aneuploidy. Subgroup 8 was defined by predominantly balanced genomes with isolated isochromosome 17q (p < 0·0001). While no mutations were associated with outcome and overall mutational burden was low, WCA-HR harboured recurrent chromatin remodelling mutations (p = 0·007). Integration of methylation and WCA groups improved risk-stratification models and outperformed established prognostication schemes. Our MBGrp4 risk-stratification scheme defines: favourable-risk (non-metastatic disease and (i) subgroup 7 or (ii) WCA-FR (21% of patients, 5-year PFS 97%)), very-high-risk (metastatic disease with WCA-HR (36%, 5-year PFS 49%)) and high-risk (remaining patients; 43%, 5-year PFS 67%). These findings validated in an independent MBGrp4 cohort (n = 668). Importantly, our findings demonstrate that previously established disease-wide risk-features (i.e. LCA histology and MYC(N) amplification) have little prognostic relevance in MBGrp4 disease. Novel validated survival models, integrating clinical features, methylation and WCA groups, improve outcome prediction and re-define risk-status for ~ 80% of MBGrp4. Our MBGrp4 favourable-risk group has MBWNT-like excellent outcomes, thereby doubling the proportion of medulloblastoma patients who could benefit from therapy de-escalation approaches, aimed at reducing treatment induced late-effects while sustaining survival outcomes. Novel approaches are urgently required for the very-high-risk patients.

Fucosylation of HLA-DRB1 regulates CD4+ T cell-mediated anti-melanoma immunity and enhances immunotherapy efficacy.

Immunotherapy efficacy is limited in melanoma, and combinations of immunotherapies with other modalities have yielded limited improvements but also adverse events requiring cessation of treatment. In addition to ineffective patient stratification, efficacy is impaired by paucity of intratumoral immune cells (itICs); thus, effective strategies to safely increase itICs are needed. We report that dietary administration of L-fucose induces fucosylation and cell surface enrichment of the major histocompatibility complex (MHC)-II protein HLA-DRB1 in melanoma cells, triggering CD4+ T cell-mediated increases in itICs and anti-tumor immunity, enhancing immune checkpoint blockade responses. Melanoma fucosylation and fucosylated HLA-DRB1 associate with intratumoral T cell abundance and anti-programmed cell death protein 1 (PD1) responder status in patient melanoma specimens, suggesting the potential use of melanoma fucosylation as a strategy for stratifying patients for immunotherapies. Our findings demonstrate that fucosylation is a key mediator of anti-tumor immunity and, importantly, suggest that L-fucose is a powerful agent for safely increasing itICs and immunotherapy efficacy in melanoma.

Membrane Fusion Mediated by Non-covalent Binding of Re-engineered Cholera Toxin Assemblies to Glycolipids.

Membrane fusion is essential for the transport of macromolecules and viruses across membranes. While glycan-binding proteins (lectins) often initiate cellular adhesion, subsequent fusion events require additional protein machinery. No mechanism for membrane fusion arising from simply a protein binding to membrane glycolipids has been described thus far. Herein, we report that a biotinylated protein derived from cholera toxin becomes a fusogenic lectin upon cross-linking with streptavidin. This novel reengineered protein brings about hemifusion and fusion of vesicles as demonstrated by mixing of fluorescently labeled lipids between vesicles as well as content mixing of liposomes filled with fluorescently labeled dextran. Exclusion of the complex at vesicle-vesicle interfaces could also be observed, indicating the formation of hemifusion diaphragms. Discovery of this fusogenic lectin complex demonstrates that new emergent properties can arise from simple changes in protein architecture and provides insights into new mechanisms of lipid-driven fusion.

Acute-Onset Footdrop Caused by Intraneural Ganglion Cyst of the Common Peroneal Nerve: The Effects of Extraneural Pressure Gradients on Cyst Propagation.

Ganglion cysts are relatively common entities, but intraneural ganglia within peripheral nerves are rare and poorly understood. We present a case of a 51-year-old man who presented with acute left dropfoot. Initial magnetic resonance imaging (MRI) was misinterpreted as common peroneal neuritis consistent with a traction injury corroborated by the patient's history. However, after surgical decompression and external neurolysis were performed, the patient's symptoms worsened. Repeated MRI revealed an intraneural ganglion cyst of the common peroneal nerve with connection to the superior tibiofibular joint by means of its anterior recurrent branch that was evident retrospectively on preoperative MRI. It is crucial to carefully inspect atypical cases to further recognize and appreciate the dynamic aspect of this disease or "roller-coaster" phenomenon. Intraneural ganglion cysts rely heavily on intraneural and extraneural pressure gradients for propagation, which can be drawn from the expanded work of the unifying articular theory. This report emphasizes the importance of understanding the pathoanatomical and hydraulic factors to appropriately identify and treat intraneural ganglion cysts. Increased recognition of this pathologic entity as a differential diagnosis for acute onset dropfoot is also highlighted.

Medulloblastoma group 3 and 4 tumors comprise a clinically and biologically significant expression continuum reflecting human cerebellar development.

Medulloblastoma is currently subclassified into distinct DNA methylation subgroups/subtypes with particular clinico-molecular features. Using RNA sequencing (RNA-seq) in large, well-annotated cohorts of medulloblastoma, we show that transcriptionally group 3 and group 4 medulloblastomas exist as intermediates on a bipolar continuum between archetypal group 3 and group 4 entities. Continuum position is prognostic, reflecting a propensity for specific DNA copy-number changes, and specific switches in isoform/enhancer usage and RNA editing. Examining single-cell RNA-seq (scRNA-seq) profiles, we show that intratumoral transcriptional heterogeneity along the continuum is limited in a subtype-dependent manner. By integrating with a human scRNA-seq reference atlas, we show that this continuum is mirrored by an equivalent continuum of transcriptional cell types in early fetal cerebellar development. We identify distinct developmental niches for all four major subgroups and link each to a common developmental antecedent. Our findings show a transcriptional continuum arising from oncogenic disruption of highly specific fetal cerebellar cell types, linked to almost every aspect of group 3/group 4 molecular biology and clinico-pathology.

Single-cell DNA sequencing identifies risk-associated clonal complexity and evolutionary trajectories in childhood medulloblastoma development.

We reconstructed the natural history and temporal evolution of the most common childhood brain malignancy, medulloblastoma, by single-cell whole-genome sequencing (sc-WGS) of tumours representing its major molecular sub-classes and clinical risk groups. Favourable-risk disease sub-types assessed (MBWNT and infant desmoplastic/nodular MBSHH) typically comprised a single clone with no evidence of further evolution. In contrast, highest risk sub-classes (MYC-amplified MBGroup3 and TP53-mutated MBSHH) were most clonally diverse and displayed gradual evolutionary trajectories. Clinically adopted biomarkers (e.g. chromosome 6/17 aberrations; CTNNB1/TP53 mutations) were typically early-clonal/initiating events, exploitable as targets for early-disease detection; in analyses of spatially distinct tumour regions, a single biopsy was sufficient to assess their status. Importantly, sc-WGS revealed novel events which arise later and/or sub-clonally and more commonly display spatial diversity; their clinical significance and role in disease evolution post-diagnosis now require establishment. These findings reveal diverse modes of tumour initiation and evolution in the major medulloblastoma sub-classes, with pathogenic relevance and clinical potential.

Epigenetic regulator genes direct lineage switching in MLL/AF4 leukemia.

The fusion gene MLL/AF4 defines a high-risk subtype of pro-B acute lymphoblastic leukemia. Relapse can be associated with a lineage switch from acute lymphoblastic to acute myeloid leukemia, resulting in poor clinical outcomes caused by resistance to chemotherapies and immunotherapies. In this study, the myeloid relapses shared oncogene fusion breakpoints with their matched lymphoid presentations and originated from various differentiation stages from immature progenitors through to committed B-cell precursors. Lineage switching is linked to substantial changes in chromatin accessibility and rewiring of transcriptional programs, including alternative splicing. These findings indicate that the execution and maintenance of lymphoid lineage differentiation is impaired. The relapsed myeloid phenotype is recurrently associated with the altered expression, splicing, or mutation of chromatin modifiers, including CHD4 coding for the ATPase/helicase of the nucleosome remodelling and deacetylation complex. Perturbation of CHD4 alone or in combination with other mutated epigenetic modifiers induces myeloid gene expression in MLL/AF4+ cell models, indicating that lineage switching in MLL/AF4 leukemia is driven and maintained by disrupted epigenetic regulation.