UBQLN4 Represses Homologous Recombination and Is Overexpressed in Aggressive Tumors.

Genomic instability can be a hallmark of both human genetic disease and cancer. We identify a deleterious UBQLN4 mutation in families with an autosomal recessive syndrome reminiscent of genome instability disorders. UBQLN4 deficiency leads to increased sensitivity to genotoxic stress and delayed DNA double-strand break (DSB) repair. The proteasomal shuttle factor UBQLN4 is phosphorylated by ATM and interacts with ubiquitylated MRE11 to mediate early steps of homologous recombination-mediated DSB repair (HRR). Loss of UBQLN4 leads to chromatin retention of MRE11, promoting non-physiological HRR activity in vitro and in vivo. Conversely, UBQLN4 overexpression represses HRR and favors non-homologous end joining. Moreover, we find UBQLN4 overexpressed in aggressive tumors. In line with an HRR defect in these tumors, UBQLN4 overexpression is associated with PARP1 inhibitor sensitivity. UBQLN4 therefore curtails HRR activity through removal of MRE11 from damaged chromatin and thus offers a therapeutic window for PARP1 inhibitor treatment in UBQLN4-overexpressing tumors.

Mutations in ALK signaling pathways conferring resistance to ALK inhibitor treatment lead to collateral vulnerabilities in neuroblastoma cells.

BACKGROUND: Development of resistance to targeted therapies has tempered initial optimism that precision oncology would improve poor outcomes for cancer patients. Resistance mechanisms, however, can also confer new resistance-specific vulnerabilities, termed collateral sensitivities. Here we investigated anaplastic lymphoma kinase (ALK) inhibitor resistance in neuroblastoma, a childhood cancer frequently affected by activating ALK alterations. METHODS: Genome-wide forward genetic CRISPR-Cas9 based screens were performed to identify genes associated with ALK inhibitor resistance in neuroblastoma cell lines. Furthermore, the neuroblastoma cell line NBLW-R was rendered resistant by continuous exposure to ALK inhibitors. Genes identified to be associated with ALK inhibitor resistance were further investigated by generating suitable cell line models. In addition, tumor and liquid biopsy samples of four patients with ALK-mutated neuroblastomas before ALK inhibitor treatment and during tumor progression under treatment were genomically profiled. RESULTS: Both genome-wide CRISPR-Cas9-based screens and preclinical spontaneous ALKi resistance models identified NF1 loss and activating NRASQ61K mutations to confer resistance to chemically diverse ALKi. Moreover, human neuroblastomas recurrently developed de novo loss of NF1 and activating RAS mutations after ALKi treatment, leading to therapy resistance. Pathway-specific perturbations confirmed that NF1 loss and activating RAS mutations lead to RAS-MAPK signaling even in the presence of ALKi. Intriguingly, NF1 loss rendered neuroblastoma cells hypersensitive to MEK inhibition. CONCLUSIONS: Our results provide a clinically relevant mechanistic model of ALKi resistance in neuroblastoma and highlight new clinically actionable collateral sensitivities in resistant cells.

Neuroblastoma signalling models unveil combination therapies targeting feedback-mediated resistance.

Very high risk neuroblastoma is characterised by increased MAPK signalling, and targeting MAPK signalling is a promising therapeutic strategy. We used a deeply characterised panel of neuroblastoma cell lines and found that the sensitivity to MEK inhibitors varied drastically between these cell lines. By generating quantitative perturbation data and mathematical modelling, we determined potential resistance mechanisms. We found that negative feedbacks within MAPK signalling and via the IGF receptor mediate re-activation of MAPK signalling upon treatment in resistant cell lines. By using cell-line specific models, we predict that combinations of MEK inhibitors with RAF or IGFR inhibitors can overcome resistance, and tested these predictions experimentally. In addition, phospho-proteomic profiling confirmed the cell-specific feedback effects and synergy of MEK and IGFR targeted treatment. Our study shows that a quantitative understanding of signalling and feedback mechanisms facilitated by models can help to develop and optimise therapeutic strategies. Our findings should be considered for the planning of future clinical trials introducing MEKi in the treatment of neuroblastoma.

Telomerase activation by genomic rearrangements in high-risk neuroblastoma.

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system. Roughly half of these tumours regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavourable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive. Here we have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n = 39; low-risk, n = 17) and discovered recurrent genomic rearrangements affecting a chromosomal region at 5p15.33 proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in a mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumour type. In an extended case series (n = 217), TERT rearrangements defined a subgroup of high-risk tumours with particularly poor outcome. Despite a large structural diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumours, TERT expression was also elevated in MYCN-amplified tumours, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.33 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodelling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cell lines bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodelling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the centre of transformation in a large fraction of these tumours.

Chromosome 17/17q gain and unaltered profiles in high resolution array-CGH are prognostically informative in neuroblastoma.

The prognostic relevance of chromosome 17 gain in neuroblastoma is still discussed. This investigation specifies the frequency, type, size, and transcriptional relevance in a large patient cohort. Primary tumor material of 202 patients was analyzed using high-resolution oligonucleotide array-based comparative genomic hybridization (aCGH) and correlated with clinical and survival data. A subset (n = 145) was correlated for differentially expressed genes (DEG) by microarray analysis. Chromosome 17 aCGH analysis showed numerical gain in 94/202 patients (47%), partial gain in 93/202 patients (46%), and no gain in 15/202 patients (7%). The frequency of partial gain was higher in stage 4 neuroblastoma (stage 1 15%; stage 2 12%; stage 3 16%; stage 4S 7%; and stage 4 50%). Overall survival (OS) was superior in patients with numerical gain compared with patients with partial gain or no gain (5-y-OS: 0.95 ± 0.02 vs. 0.63 ± 0.05 vs. 0.60 ± 0.13; P < 0.001). Gene expression analysis demonstrated 95/130 DEGs between tumors with numerical or partial chromosome/no gain. Only one DEG (CCKBR) was detected comparing tumors with partial gain and those with no gain. In patients with partial gain, the distribution of breakpoints did not correlate with stage and 11q status, but with MYCN amplification and 1p status. The "best" breakpoints in cases with partial 17q gain were at 42.5 Mb for event-free and 26.6 Mb for OS. Numerical gain of chromosome 17 is associated with a better prognosis than partial and no gain. The group of tumors with partial gain was similar to the group without gain with respect to stage distribution, outcome, and gene expression profile.

Smad4 binds Hoxa9 in the cytoplasm and protects primitive hematopoietic cells against nuclear activation by Hoxa9 and leukemia transformation.

We studied leukemic stem cells (LSCs) in a Smad4(-/-) mouse model of acute myelogenous leukemia (AML) induced either by the HOXA9 gene or by the fusion oncogene NUP98-HOXA9. Although Hoxa9-Smad4 complexes accumulate in the cytoplasm of normal hematopoietic stem cells and progenitor cells (HSPCs) transduced with these oncogenes, there is no cytoplasmic stabilization of HOXA9 in Smad4(-/-) HSPCs, and as a consequence increased levels of Hoxa9 is observed in the nucleus leading to increased immortalization in vitro. Loss of Smad4 accelerates the development of leukemia in vivo because of an increase in transformation of HSPCs. Therefore, the cytoplasmic binding of Hoxa9 by Smad4 is a mechanism to protect Hoxa9-induced transformation of normal HSPCs. Because Smad4 is a potent tumor suppressor involved in growth control, we developed a strategy to modify the subcellular distribution of Smad4. We successfully disrupted the interaction between Hoxa9 and Smad4 to activate the TGF-ß pathway and apoptosis, leading to a loss of LSCs. Together, these findings reveal a major role for Smad4 in the negative regulation of leukemia initiation and maintenance induced by HOXA9/NUP98-HOXA9 and provide strong evidence that antagonizing Smad4 stabilization by these oncoproteins might be a promising novel therapeutic approach in leukemia.

Mutational topography reflects clinical neuroblastoma heterogeneity.

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

Defining the structure of the general stress regulon of Bacillus subtilis using targeted microarray analysis and random forest classification.

The structure of the SigB-dependent general stress regulon of Bacillus subtilis has previously been characterized by proteomics approaches as well as DNA array-based expression studies. However, comparing the SigB targets published in three previous major transcriptional profiling studies it is obvious that although each of them identified well above 100 target genes, only 67 were identified in all three studies. These substantial differences can likely be attributed to the different strains, growth conditions, microarray platforms and experimental setups used in the studies. In order to gain a better understanding of the structure of this important regulon, a targeted DNA microarray analysis covering most of the known SigB-inducing conditions was performed, and the changes in expression kinetics of 252 potential members of the SigB regulon and appropriate control genes were recorded. Transcriptional data for the B. subtilis wild-type strain 168 and its isogenic sigB mutant BSM29 were analysed using random forest, a machine learning algorithm, by incorporating the knowledge from previous studies. This analysis revealed a strictly SigB-dependent expression pattern for 166 genes following ethanol, butanol, osmotic and oxidative stress, low-temperature growth and heat shock, as well as limitation of oxygen or glucose. Kinetic analysis of the data for the wild-type strain identified 30 additional members of the SigB regulon, which were also subject to control by additional transcriptional regulators, thus displaying atypical SigB-independent induction patterns in the mutant strain under some of the conditions tested. For 19 of these 30 SigB regulon members, published reports support control by secondary regulators along with SigB. Thus, this microarray-based study assigns a total of 196 genes to the SigB-dependent general stress regulon of B. subtilis.

Prospectively isolated CD133/CD24-positive ependymal cells from the adult spinal cord and lateral ventricle wall differ in their long-term in vitro self-renewal and in vivo gene expression.

In contrast to ependymal cells located above the subventricular zone (SVZ) of the adult lateral ventricle wall (LVW), adult spinal cord (SC) ependymal cells possess certain neural stem cell characteristics. The molecular basis of this difference is unknown. In this study, antibodies against multiple cell surface markers were applied to isolate pure populations of SC and LVW ependymal cells, which allowed a direct comparison of their in vitro behavior and in vivo gene expression profile. Isolated CD133(+)/CD24(+)/CD45(-)/CD34(-) ependymal cells from the SC displayed in vitro self-renewal and differentiation capacity, whereas those from the LVW did not. SC ependymal cells showed a higher expression of several genes involved in cell division, cell cycle regulation, and chromosome stability, which is consistent with a long-term self-renewal capacity, and shared certain transcripts with neural stem cells of the embryonic forebrain. They also expressed several retinoic acid (RA)-regulated genes and responded to RA exposure. LVW ependymal cells showed higher transcript levels of many genes regulated by transforming growth factor-ß family members. Among them were Dlx2, Id2, Hey1, which together with Foxg1 could explain their potential to turn into neuroblasts under certain environmental conditions.

Chromothripsis followed by circular recombination drives oncogene amplification in human cancer.

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. Using whole-genome sequencing data of the pediatric tumor neuroblastoma, we here identified a type of amplification, termed 'seismic amplification', that is characterized by multiple rearrangements and discontinuous copy number levels. Overall, seismic amplifications occurred in 9.9% (274 of 2,756) of cases across 38 cancer types, and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification showed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and overexpression in a substantial fraction of human malignancies.

Spatial and temporal intratumour heterogeneity has potential consequences for single biopsy-based neuroblastoma treatment decisions.

Intratumour heterogeneity is a major cause of treatment failure in cancer. We present in-depth analyses combining transcriptomic and genomic profiling with ultra-deep targeted sequencing of multiregional biopsies in 10 patients with neuroblastoma, a devastating childhood tumour. We observe high spatial and temporal heterogeneity in somatic mutations and somatic copy-number alterations which are reflected on the transcriptomic level. Mutations in some druggable target genes including ALK and FGFR1 are heterogeneous at diagnosis and/or relapse, raising the issue whether current target prioritization and molecular risk stratification procedures in single biopsies are sufficiently reliable for therapy decisions. The genetic heterogeneity in gene mutations and chromosome aberrations observed in deep analyses from patient courses suggest clonal evolution before treatment and under treatment pressure, and support early emergence of metastatic clones and ongoing chromosomal instability during disease evolution. We report continuous clonal evolution on mutational and copy number levels in neuroblastoma, and detail its implications for therapy selection, risk stratification and therapy resistance.

Neuroblastoma Risk Assessment and Treatment Stratification with Hybrid Capture-Based Panel Sequencing.

For many years, the risk-based therapy stratification of children with neuroblastoma has relied on clinical and molecular covariates. In recent years, genome analysis has revealed further alterations defining risk, tumor biology, and therapeutic targets. The implementation of a robust and scalable method for analyzing traditional and new molecular markers in routine diagnostics is an urgent clinical need. Here, we investigated targeted panel sequencing as a diagnostic approach to analyze all relevant genomic neuroblastoma risk markers in one assay. Our "neuroblastoma hybrid capture sequencing panel" (NB-HCSP) assay employs a technology for the high-coverage sequencing (>1000×) of 55 selected genes and neuroblastoma-relevant genomic regions, which allows for the detection of single nucleotide changes, structural rearrangements, and copy number alterations. We validated our assay by analyzing 15 neuroblastoma cell lines and a cohort of 20 neuroblastomas, for which reference routine diagnostic data and genome sequencing data were available. We observed a high concordance for risk markers identified by the NB-HSCP assay, clinical routine diagnostics, and genome sequencing. Subsequently, we demonstrated clinical applicability of the NB-HCSP assay by analyzing routine clinical samples. We conclude that the NB-HCSP assay may be implemented into routine diagnostics as a single assay that covers all essential covariates for initial neuroblastoma classification, extended risk stratification, and targeted therapy selection.

Enhancer hijacking determines extrachromosomal circular MYCN amplicon architecture in neuroblastoma.

MYCN amplification drives one in six cases of neuroblastoma. The supernumerary gene copies are commonly found on highly rearranged, extrachromosomal circular DNA (ecDNA). The exact amplicon structure has not been described thus far and the functional relevance of its rearrangements is unknown. Here, we analyze the MYCN amplicon structure using short-read and Nanopore sequencing and its chromatin landscape using ChIP-seq, ATAC-seq and Hi-C. This reveals two distinct classes of amplicons which explain the regulatory requirements for MYCN overexpression. The first class always co-amplifies a proximal enhancer driven by the noradrenergic core regulatory circuit (CRC). The second class of MYCN amplicons is characterized by high structural complexity, lacks key local enhancers, and instead contains distal chromosomal fragments harboring CRC-driven enhancers. Thus, ectopic enhancer hijacking can compensate for the loss of local gene regulatory elements and explains a large component of the structural diversity observed in MYCN amplification.

Publisher Correction: Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

Extrachromosomal circularization of DNA is an important genomic feature in cancer. However, the structure, composition and genome-wide frequency of extrachromosomal circular DNA have not yet been profiled extensively. Here, we combine genomic and transcriptomic approaches to describe the landscape of extrachromosomal circular DNA in neuroblastoma, a tumor arising in childhood from primitive cells of the sympathetic nervous system. Our analysis identifies and characterizes a wide catalog of somatically acquired and undescribed extrachromosomal circular DNAs. Moreover, we find that extrachromosomal circular DNAs are an unanticipated major source of somatic rearrangements, contributing to oncogenic remodeling through chimeric circularization and reintegration of circular DNA into the linear genome. Cancer-causing lesions can emerge out of circle-derived rearrangements and are associated with adverse clinical outcome. It is highly probable that circle-derived rearrangements represent an ongoing mutagenic process. Thus, extrachromosomal circular DNAs represent a multihit mutagenic process, with important functional and clinical implications for the origins of genomic remodeling in cancer.

CD133 is not present on neurogenic astrocytes in the adult subventricular zone, but on embryonic neural stem cells, ependymal cells, and glioblastoma cells.

Human brain tumor stem cells have been enriched using antibodies against the surface protein CD133. An antibody recognizing CD133 also served to isolate normal neural stem cells from fetal human brain, suggesting a possible lineage relationship between normal neural and brain tumor stem cells. Whether CD133-positive brain tumor stem cells can be derived from CD133-positive neural stem or progenitor cells still requires direct experimental evidence, and an important step toward such investigations is the identification and characterization of normal CD133-presenting cells in neurogenic regions of the embryonic and adult brain. Here, we present evidence that CD133 is a marker for embryonic neural stem cells, an intermediate radial glial/ependymal cell type in the early postnatal stage, and for ependymal cells in the adult brain, but not for neurogenic astrocytes in the adult subventricular zone. Our findings suggest two principal possibilities for the origin of brain tumor stem cells: a derivation from CD133-expressing cells, which are normally not present in the adult brain (embryonic neural stem cells and an early postnatal intermediate radial glial/ependymal cell type), or from CD133-positive ependymal cells in the adult brain, which are, however, generally regarded as postmitotic. Alternatively, brain tumor stem cells could be derived from proliferative but CD133-negative neurogenic astrocytes in the adult brain. In the latter case, brain tumor development would involve the production of CD133.

Telomerase Is a Prognostic Marker of Poor Outcome and a Therapeutic Target in Neuroblastoma.

PURPOSE: Telomere maintenance is a hallmark of high-risk neuroblastoma; however, the contribution of telomerase and alternative lengthening of telomeres (ALT) to clinical phenotypes has remained unclear. We aimed to determine the clinical relevance of telomerase activation versus ALT as biomarkers in pretreatment neuroblastoma and to assess the potential value of telomerase as a therapeutic target. MATERIALS AND METHODS: The genomic status of TERT and MYCN was assessed in 457 pretreatment neuroblastomas by fluorescence in situ hybridization. ALT was examined in 273 of 457 tumors by detection of ALT-associated promyelocytic leukemia nuclear bodies, and TERT expression was determined by microarrays in 223 of these. Cytotoxic effects of telomerase-interacting compounds were analyzed in neuroblastoma cell lines in vitro and in vivo. RESULTS: We detected TERT rearrangements in 46 of 457 cases (10.1%), MYCN amplification in 93 of 457 cases (20.4%), and elevated TERT expression in tumors lacking TERT or MYCN alterations in 10 of 223 cases (4.5%). ALT activation was found in 49 of 273 cases (17.9%). All these alterations occurred almost mutually exclusively and were associated with unfavorable prognostic variables and adverse outcome. The presence of activated telomerase (ie, TERT rearrangements, MYCN amplification, or high TERT expression without these alterations) was associated with poorest overall survival and was an independent prognostic marker in multivariable analyses. We also found that the telomerase-interacting compound 6-thio-2'-deoxyguanosine effectively inhibited viability and proliferation of neuroblastoma cells bearing activated telomerase. Similarly, tumor growth was strongly impaired upon 6-thio-2'-deoxyguanosine treatment in telomerase-positive neuroblastoma xenografts in mice. CONCLUSION: Our data suggest telomerase activation and ALT define distinct neuroblastoma subgroups with adverse outcome and that telomerase may represent a promising therapeutic target in many high-risk neuroblastomas.

The mutational landscape of MYCN, Lin28b and ALKF1174L driven murine neuroblastoma mimics human disease.

Genetically engineered mouse models have proven to be essential tools for unraveling fundamental aspects of cancer biology and for testing novel therapeutic strategies. To optimally serve these goals, it is essential that the mouse model faithfully recapitulates the human disease. Recently, novel mouse models for neuroblastoma have been developed. Here, we report on the further genomic characterization through exome sequencing and DNA copy number analysis of four of the currently available murine neuroblastoma model systems (ALK, Th-MYCN, Dbh-MYCN and Lin28b). The murine tumors revealed a low number of genomic alterations - in keeping with human neuroblastoma - and a positive correlation of the number of genetic lesions with the time to onset of tumor formation was observed. Gene copy number alterations are the hallmark of both murine and human disease and frequently affect syntenic genomic regions. Despite low mutational load, the genes mutated in murine disease were found to be enriched for genes mutated in human disease. Taken together, our study further supports the validity of the tested mouse models for mechanistic and preclinical studies of human neuroblastoma.

Sense-Antisense lncRNA Pair Encoded by Locus 6p22.3 Determines Neuroblastoma Susceptibility via the USP36-CHD7-SOX9 Regulatory Axis.

Trait-associated loci often map to genomic regions encoding long noncoding RNAs (lncRNAs), but the role of these lncRNAs in disease etiology is largely unexplored. We show that a pair of sense/antisense lncRNA (6p22lncRNAs) encoded by CASC15 and NBAT1 located at the neuroblastoma (NB) risk-associated 6p22.3 locus are tumor suppressors and show reduced expression in high-risk NBs. Loss of functional synergy between 6p22lncRNAs results in an undifferentiated state that is maintained by a gene-regulatory network, including SOX9 located on 17q, a region frequently gained in NB. 6p22lncRNAs regulate SOX9 expression by controlling CHD7 stability via modulating the cellular localization of USP36, encoded by another 17q gene. This regulatory nexus between 6p22.3 and 17q regions may lead to potential NB treatment strategies.