Golgi-localized Ring Finger Protein 121 is necessary for MYCN-driven neuroblastoma tumorigenesis.

MYCN amplification predicts poor prognosis in childhood neuroblastoma. To identify MYCN oncogenic signal dependencies we performed N-ethyl-N-nitrosourea (ENU) mutagenesis on the germline of neuroblastoma-prone TH-MYCN transgenic mice to generate founders which had lost tumorigenesis. Sequencing of the mutant mouse genomes identified the Ring Finger Protein 121 (RNF121WT) gene mutated to RNFM158R associated with heritable loss of tumorigenicity. While the RNF121WT protein localised predominantly to the cis-Golgi Complex, the RNF121M158R mutation in Helix 4 of its transmembrane domain caused reduced RNF121 protein stability and absent Golgi localisation. RNF121WT expression markedly increased during TH-MYCN tumorigenesis, whereas hemizygous RNF121WT gene deletion reduced TH-MYCN tumorigenicity. The RNF121WT-enhanced growth of MYCN-amplified neuroblastoma cells depended on RNF121WT transmembrane Helix 5. RNF121WT directly bound MYCN protein and enhanced its stability. High RNF121 mRNA expression associated with poor prognosis in human neuroblastoma tissues and another MYC-driven malignancy, laryngeal cancer. RNF121 is thus an essential oncogenic cofactor for MYCN and a target for drug development.

The transcriptional co-repressor Runx1t1 is essential for MYCN-driven neuroblastoma tumorigenesis.

MYCN oncogene amplification is frequently observed in aggressive childhood neuroblastoma. Using an unbiased large-scale mutagenesis screen in neuroblastoma-prone transgenic mice, we identify a single germline point mutation in the transcriptional corepressor Runx1t1, which abolishes MYCN-driven tumorigenesis. This loss-of-function mutation disrupts a highly conserved zinc finger domain within Runx1t1. Deletion of one Runx1t1 allele in an independent Runx1t1 knockout mouse model is also sufficient to prevent MYCN-driven neuroblastoma development, and reverse ganglia hyperplasia, a known pre-requisite for tumorigenesis. Silencing RUNX1T1 in human neuroblastoma cells decreases colony formation in vitro, and inhibits tumor growth in vivo. Moreover, RUNX1T1 knockdown inhibits the viability of PAX3-FOXO1 fusion-driven rhabdomyosarcoma and MYC-driven small cell lung cancer cells. Despite the role of Runx1t1 in MYCN-driven tumorigenesis neither gene directly regulates the other. We show RUNX1T1 forms part of a transcriptional LSD1-CoREST3-HDAC repressive complex recruited by HAND2 to enhancer regions to regulate chromatin accessibility and cell-fate pathway genes.

Precision-guided treatment in high-risk pediatric cancers.

Recent research showed that precision medicine can identify new treatment strategies for patients with childhood cancers. However, it is unclear which patients will benefit most from precision-guided treatment (PGT). Here we report consecutive data from 384 patients with high-risk pediatric cancer (with an expected cure rate of less than 30%) who had at least 18 months of follow-up on the ZERO Childhood Cancer Precision Medicine Program PRecISion Medicine for Children with Cancer (PRISM) trial. A total of 256 (67%) patients received PGT recommendations and 110 (29%) received a recommended treatment. PGT resulted in a 36% objective response rate and improved 2-year progression-free survival compared with standard of care (26% versus 12%; P = 0.049) or targeted agents not guided by molecular findings (26% versus 5.2%; P = 0.003). PGT based on tier 1 evidence, PGT targeting fusions or commenced before disease progression had the greatest clinical benefit. Our data show that PGT informed by comprehensive molecular profiling significantly improves outcomes for children with high-risk cancers. ClinicalTrials.gov registration: NCT03336931.

A novel network-based method identifies a cuproplasia-related pan-cancer gene signature to predict patient outcome.

Copper is a vital micronutrient involved in many biological processes and is an essential component of tumour cell growth and migration. Copper influences tumour growth through a process called cuproplasia, defined as abnormal copper-dependent cell-growth and proliferation. Copper-chelation therapy targeting this process has demonstrated efficacy in several clinical trials against cancer. While the molecular pathways associated with cuproplasia are partially known, genetic heterogeneity across different cancer types has limited the understanding of how cuproplasia impacts patient survival. Utilising RNA-sequencing data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) datasets, we generated gene regulatory networks to identify the critical cuproplasia-related genes across 23 different cancer types. From this, we identified a novel 8-gene cuproplasia-related gene signature associated with pan-cancer survival, and a 6-gene prognostic risk score model in low grade glioma. These findings highlight the use of gene regulatory networks to identify cuproplasia-related gene signatures that could be used to generate risk score models. This can potentially identify patients who could benefit from copper-chelation therapy and identifies novel targeted therapeutic strategies.

Analytical Workflows for Single-Cell Multiomic Data Using the BD Rhapsody Platform.

The conversion of raw sequencing reads to biologically relevant data in high-throughput single-cell RNA sequencing experiments is a complex and involved process. Drawing meaning from thousands of individual cells to provide biological insight requires ensuring not only that the data are of the highest quality but also that the signal can be separated from noise. In this article, we describe a detailed analytical workflow, including six pipelines, that allows high-quality data analysis in single-cell multiomics. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Image analysis Basic Protocol 2: Sequencing quality control and generation of a gene expression matrix Basic Protocol 3: Gene expression matrix data pre-processing and analysis Basic Protocol 4: Advanced analysis Basic Protocol 5: Conversion to flow cytometry standard (FCS) format Basic Protocol 6: Visualization using graphical interfaces.

MYCN and SNRPD3 cooperate to maintain a balance of alternative splicing events that drives neuroblastoma progression.

Many of the pro-tumorigenic functions of the oncogene MYCN are attributed to its regulation of global gene expression programs. Alternative splicing is another important regulator of gene expression and has been implicated in neuroblastoma development, however, the molecular mechanisms remain unknown. We found that MYCN up-regulated the expression of the core spliceosomal protein, SNRPD3, in models of neuroblastoma initiation and progression. High mRNA expression of SNRPD3 in human neuroblastoma tissues was a strong, independent prognostic factor for poor patient outcome. Repression of SNRPD3 expression correlated with loss of colony formation in vitro and reduced tumorigenicity in vivo. The effect of SNRPD3 on cell viability was in part dependent on MYCN as an oncogenic co-factor. RNA-sequencing revealed a global increase in the number of genes being differentially spliced when MYCN was overexpressed. Surprisingly, depletion of SNRPD3 in the presence of overexpressed MYCN further increased differential splicing, particularly of cell cycle regulators, such as BIRC5 and CDK10. MYCN directly bound SNRPD3, and the protein arginine methyltransferase, PRMT5, consequently increasing SNRPD3 methylation. Indeed, the PRMT5 inhibitor, JNJ-64619178, reduced cell viability and SNRPD3 methylation in neuroblastoma cells with high SNRPD3 and MYCN expression. Our findings demonstrate a functional relationship between MYCN and SNRPD3, which maintains the fidelity of MYCN-driven alternative splicing in the narrow range required for neuroblastoma cell growth. SNRPD3 methylation and its protein-protein interface with MYCN represent novel therapeutic targets. Hypothetical model for SNRPD3 as a co-factor for MYCN oncogenesis. SNRPD3 and MYCN participate in a regulatory loop to balance splicing fidelity in neuroblastoma cells. First MYCN transactivates SNRPD3 to lead to high-level expression. Second, SNRPD3 and MYCN form a protein complex involving PRMT5. Third, this leads to balanced alterative splicing (AS) activitiy that is favorable to neuroblastoma. Together this forms as a therapeutic vulnerability where SNRPD3 perturbation or PRMT5 inhibitors are selectively toxic to neuroblastoma by conditionally disturbing splicing activity.

High-Throughput Drug Screening of Primary Tumor Cells Identifies Therapeutic Strategies for Treating Children with High-Risk Cancer.

For one-third of patients with pediatric cancer enrolled in precision medicine programs, molecular profiling does not result in a therapeutic recommendation. To identify potential strategies for treating these high-risk pediatric patients, we performed in vitro screening of 125 patient-derived samples against a library of 126 anticancer drugs. Tumor cell expansion did not influence drug responses, and 82% of the screens on expanded tumor cells were completed while the patients were still under clinical care. High-throughput drug screening (HTS) confirmed known associations between activating genomic alterations in NTRK, BRAF, and ALK and responses to matching targeted drugs. The in vitro results were further validated in patient-derived xenograft models in vivo and were consistent with clinical responses in treated patients. In addition, effective combinations could be predicted by correlating sensitivity profiles between drugs. Furthermore, molecular integration with HTS identified biomarkers of sensitivity to WEE1 and MEK inhibition. Incorporating HTS into precision medicine programs is a powerful tool to accelerate the improved identification of effective biomarker-driven therapeutic strategies for treating high-risk pediatric cancers. SIGNIFICANCE: Integrating HTS with molecular profiling is a powerful tool for expanding precision medicine to support drug treatment recommendations and broaden the therapeutic options available to high-risk pediatric cancers.

Introme accurately predicts the impact of coding and noncoding variants on gene splicing, with clinical applications.

Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .

A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer.

BACKGROUND: Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. METHODS: We combined immunohistochemistry (IHC) with RNA sequencing and whole-genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Furthermore, we explored transcriptional features of immune archetypes and T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions and assessed the common adult biomarkers such as neoantigen load and TMB. RESULTS: A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. CONCLUSIONS: Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

In vivo loss of tumorigenicity in a patient-derived orthotopic xenograft mouse model of ependymoma.

Introduction: Ependymomas (EPN) are the third most common malignant brain cancer in children. Treatment strategies for pediatric EPN have remained unchanged over recent decades, with 10-year survival rates stagnating at just 67% for children aged 0-14 years. Moreover, a proportion of patients who survive treatment often suffer long-term neurological side effects as a result of therapy. It is evident that there is a need for safer, more effective treatments for pediatric EPN patients. There are ten distinct subgroups of EPN, each with their own molecular and prognostic features. To identify and facilitate the testing of new treatments for EPN, in vivo laboratory models representative of the diverse molecular subtypes are required. Here, we describe the establishment of a patient-derived orthotopic xenograft (PDOX) model of posterior fossa A (PFA) EPN, derived from a metastatic cranial lesion. Methods: Patient and PDOX tumors were analyzed using immunohistochemistry, DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. Results: Both patient and PDOX tumors classified as PFA EPN by methylation profiling, and shared similar histological features consistent with this molecular subgroup. RNA sequencing revealed that gene expression patterns were maintained across the primary and metastatic tumors, as well as the PDOX. Copy number profiling revealed gains of chromosomes 7, 8 and 19, and loss of chromosomes 2q and 6q in the PDOX and matched patient tumor. No clinically significant single nucleotide variants were identified, consistent with the low mutation rates observed in PFA EPN. Overexpression of EZHIP RNA and protein, a common feature of PFA EPN, was also observed. Despite the aggressive nature of the tumor in the patient, this PDOX was unable to be maintained past two passages in vivo. Discussion: Others who have successfully developed PDOX models report some of the lowest success rates for EPN compared to other pediatric brain cancer types attempted, with loss of tumorigenicity not uncommon, highlighting the challenges of propagating these tumors in the laboratory. Here, we discuss our collective experiences with PFA EPN PDOX model generation and propose potential approaches to improve future success in establishing preclinical EPN models.

Inhibitors of the Oncogenic PA2G4-MYCN Protein-Protein Interface.

MYCN is a major oncogenic driver for neuroblastoma tumorigenesis, yet there are no direct MYCN inhibitors. We have previously identified PA2G4 as a direct protein-binding partner of MYCN and drive neuroblastoma tumorigenesis. A small molecule known to bind PA2G4, WS6, significantly decreased tumorigenicity in TH-MYCN neuroblastoma mice, along with the inhibition of PA2G4 and MYCN interactions. Here, we identified a number of novel WS6 analogues, with 80% structural similarity, and used surface plasmon resonance assays to determine their binding affinity. Analogues #5333 and #5338 showed direct binding towards human recombinant PA2G4. Importantly, #5333 and #5338 demonstrated a 70-fold lower toxicity for normal human myofibroblasts compared to WS6. Structure-activity relationship analysis showed that a 2,3 dimethylphenol was the most suitable substituent at the R1 position. Replacing the trifluoromethyl group on the phenyl ring at the R2 position, with a bromine or hydrogen atom, increased the difference between efficacy against neuroblastoma cells and normal myofibroblast toxicity. The WS6 analogues inhibited neuroblastoma cell phenotype in vitro, in part through effects on apoptosis, while their anti-cancer effects required both PA2G4 and MYCN expression. Collectively, chemical inhibition of PA2G4-MYCN binding by WS6 analogues represents a first-in-class drug discovery which may have implications for other MYCN-driven cancers.

Histone H3-wild type diffuse midline gliomas with H3K27me3 loss are a distinct entity with exclusive EGFR or ACVR1 mutation and differential methylation of homeobox genes.

Diffuse midline gliomas (DMG) harbouring H3K27M mutation are paediatric tumours with a dismal outcome. Recently, a new subtype of midline gliomas has been described with similar features to DMG, including loss of H3K27 trimethylation, but lacking the canonical H3K27M mutation (H3-WT). Here, we report a cohort of five H3-WT tumours profiled by whole-genome sequencing, RNA sequencing and DNA methylation profiling and combine their analysis with previously published cases. We show that these tumours have recurrent and mutually exclusive mutations in either ACVR1 or EGFR and are characterised by high expression of EZHIP associated to its promoter hypomethylation. Affected patients share a similar poor prognosis as patients with H3K27M DMG. Global molecular analysis of H3-WT and H3K27M DMG reveal distinct transcriptome and methylome profiles including differential methylation of homeobox genes involved in development and cellular differentiation. Patients have distinct clinical features, with a trend demonstrating ACVR1 mutations occurring in H3-WT tumours at an older age. This in-depth exploration of H3-WT tumours further characterises this novel DMG, H3K27-altered sub-group, characterised by a specific immunohistochemistry profile with H3K27me3 loss, wild-type H3K27M and positive EZHIP. It also gives new insights into the possible mechanism and pathway regulation in these tumours, potentially opening new therapeutic avenues for these tumours which have no known effective treatment. This study has been retrospectively registered on clinicaltrial.gov on 8 November 2017 under the registration number NCT03336931 ( https://clinicaltrials.gov/ct2/show/NCT03336931 ).

Diagnostic classification of childhood cancer using multiscale transcriptomics.

The causes of pediatric cancers' distinctiveness compared to adult-onset tumors of the same type are not completely clear and not fully explained by their genomes. In this study, we used an optimized multilevel RNA clustering approach to derive molecular definitions for most childhood cancers. Applying this method to 13,313 transcriptomes, we constructed a pediatric cancer atlas to explore age-associated changes. Tumor entities were sometimes unexpectedly grouped due to common lineages, drivers or stemness profiles. Some established entities were divided into subgroups that predicted outcome better than current diagnostic approaches. These definitions account for inter-tumoral and intra-tumoral heterogeneity and have the potential of enabling reproducible, quantifiable diagnostics. As a whole, childhood tumors had more transcriptional diversity than adult tumors, maintaining greater expression flexibility. To apply these insights, we designed an ensemble convolutional neural network classifier. We show that this tool was able to match or clarify the diagnosis for 85% of childhood tumors in a prospective cohort. If further validated, this framework could be extended to derive molecular definitions for all cancer types.

Identification of Novel CSF-Derived miRNAs in Treated Paediatric Onset Spinal Muscular Atrophy: An Exploratory Study.

The availability of disease modifying therapies for spinal muscular atrophy (SMA) have created an urgent need to identify clinically meaningful biomarkers that provide insight into disease progression and therapeutic response. microRNAs (miRNA) have been shown to be involved in the pathogenesis of SMA and have the potential to provide insight within the field of SMA. miRNA-sequencing was utilized to identify differential miRNA expression in the cerebrospinal fluid (CSF) in six children with SMA treated with nusinersen in this exploratory study. Fourteen differentially expressed miRNAs were significantly altered in CSF from baseline to follow-up during treatment with nusinersen. The greatest magnitude of change was noted in miR-7-5p, miR-15a-5p, miR-15b-3p/5p, miR-126-5p, miR-128-2-5p and miR-130a-3p which encompassed a spectrum of functions predominantly in neurogenesis, neuronal differentiation and growth. The dominant signaling pathways identified in this study were the mammalian target of rapamycin and the mitogen-activated protein kinase signaling pathways. This study identified multiple miRNAs that were involved in the complex interplay between neurodevelopment and neurodegeneration.

Combined thioredoxin reductase and glutaminase inhibition exerts synergistic anti-tumor activity in MYC-high high-grade serous ovarian carcinoma.

Approximately 50%-55% of high-grade serous ovarian carcinoma (HGSOC) patients have MYC oncogenic pathway activation. Because MYC is not directly targetable, we have analyzed molecular pathways enriched in MYC-high HGSOC tumors to identify potential therapeutic targets. Here, we report that MYC-high HGSOC tumors show enrichment in genes controlled by NRF2, an antioxidant signaling pathway, along with increased thioredoxin redox activity. Treatment of MYC-high HGSOC tumors cells with US Food and Drug Administration (FDA)-approved thioredoxin reductase 1 (TrxR1) inhibitor auranofin resulted in significant growth suppression and apoptosis in MYC-high HGSOC cells in vitro and also significantly reduced tumor growth in an MYC-high HGSOC patient-derived tumor xenograft. We found that auranofin treatment inhibited glycolysis in MYC-high cells via oxidation-induced GAPDH inhibition. Interestingly, in response to auranofin-induced glycolysis inhibition, MYC-high HGSOC cells switched to glutamine metabolism for survival. Depletion of glutamine with either glutamine starvation or glutaminase (GLS1) inhibitor CB-839 exerted synergistic anti-tumor activity with auranofin in HGSOC cells and OVCAR-8 cell line xenograft. These findings suggest that applying a combined therapy of GLS1 inhibitor and TrxR1 inhibitor could effectively treat MYC-high HGSOC patients.

NONO enhances mRNA processing of super-enhancer-associated GATA2 and HAND2 genes in neuroblastoma.

High-risk neuroblastoma patients have poor survival rates and require better therapeutic options. High expression of a multifunctional DNA and RNA-binding protein, NONO, in neuroblastoma is associated with poor patient outcome; however, there is little understanding of the mechanism of NONO-dependent oncogenic gene regulatory activity in neuroblastoma. Here, we used cell imaging, biochemical and genome-wide molecular analysis to reveal complex NONO-dependent regulation of gene expression. NONO forms RNA- and DNA-tethered condensates throughout the nucleus and undergoes phase separation in vitro, modulated by nucleic acid binding. CLIP analyses show that NONO mainly binds to the 5' end of pre-mRNAs and modulates pre-mRNA processing, dependent on its RNA-binding activity. NONO regulates super-enhancer-associated genes, including HAND2 and GATA2. Abrogating NONO RNA binding, or phase separation activity, results in decreased expression of HAND2 and GATA2. Thus, future development of agents that target RNA-binding activity of NONO may have therapeutic potential in this cancer context.

Inhibition of mitochondrial translocase SLC25A5 and histone deacetylation is an effective combination therapy in neuroblastoma.

The mitochondrion is a gatekeeper of apoptotic processes, and mediates drug resistance to several chemotherapy agents used to treat cancer. Neuroblastoma is a common solid cancer in young children with poor clinical outcomes following conventional chemotherapy. We sought druggable mitochondrial protein targets in neuroblastoma cells. Among mitochondria-associated gene targets, we found that high expression of the mitochondrial adenine nucleotide translocase 2 (SLC25A5/ANT2), was a strong predictor of poor neuroblastoma patient prognosis and contributed to a more malignant phenotype in pre-clinical models. Inhibiting this transporter with PENAO reduced cell viability in a panel of neuroblastoma cell lines in a TP53-status-dependant manner. We identified the histone deacetylase inhibitor, suberanilohydroxamic acid (SAHA), as the most effective drug in clinical use against mutant TP53 neuroblastoma cells. SAHA and PENAO synergistically reduced cell viability, and induced apoptosis, in neuroblastoma cells independent of TP53-status. The SAHA and PENAO drug combination significantly delayed tumour progression in pre-clinical neuroblastoma mouse models, suggesting that these clinically advanced inhibitors may be effective in treating the disease.

Advances in CAR T cell immunotherapy for paediatric brain tumours.

Brain tumours are the most common solid tumour in children and the leading cause of cancer related death in children. Current treatments include surgery, chemotherapy and radiotherapy. The need for aggressive treatment means many survivors are left with permanent severe disability, physical, intellectual and social. Recent progress in immunotherapy, including genetically engineered T cells with chimeric antigen receptors (CARs) for treating cancer, may provide new avenues to improved outcomes for patients with paediatric brain cancer. In this review we discuss advances in CAR T cell immunotherapy, the major CAR T cell targets that are in clinical and pre-clinical development with a focus on paediatric brain tumours, the paediatric brain tumour microenvironment and strategies used to improve CAR T cell therapy for paediatric tumours.

RaScALL: Rapid (Ra) screening (Sc) of RNA-seq data for prognostically significant genomic alterations in acute lymphoblastic leukaemia (ALL).

RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, integrating RNA-seq in a clinical setting requires rapid detection and accurate reporting of clinically relevant alterations. Here we present RaScALL, an implementation of the k-mer based variant detection tool km, capable of identifying more than 100 prognostically significant lesions observed in ALL, including gene fusions, single nucleotide variants and focal gene deletions. We compared genomic alterations detected by RaScALL and those reported by alignment-based de novo variant detection tools in a study cohort of 180 Australian patient samples. Results were validated using 100 patient samples from a published North American cohort. RaScALL demonstrated a high degree of accuracy for reporting subtype defining genomic alterations. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% of reported cases in the study cohort (n = 164) and 95% of samples (n = 63) in the validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n = 12) and IKZF1 p.N159Y (n = 4). Intragenic IKZF1 deletions resulting in aberrant transcript isoforms were also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches. The application of RaScALL enables rapid identification and reporting of previously identified genomic alterations of known clinical relevance.