TERT Expression and Clinical Outcome in Pulmonary Carcinoids.

PURPOSE: The clinical course of pulmonary carcinoids ranges from indolent to fatal disease, suggesting that specific molecular alterations drive progression toward the fully malignant state. A similar spectrum of clinical phenotypes occurs in pediatric neuroblastoma, in which activation of telomerase reverse transcriptase (TERT) is decisive in determining the course of disease. We therefore investigated whether TERT expression defines the clinical fate of patients with pulmonary carcinoid. METHODS: TERT expression was examined by RNA sequencing in a test cohort and a validation cohort of pulmonary carcinoids (n = 88 and n = 105, respectively). A natural TERT expression cutoff was determined in the test cohort on the basis of the distribution of TERT expression, and its prognostic value was assessed by Kaplan-Meier survival estimates and multivariable analyses. Telomerase activity was validated by telomere repeat amplification protocol assay. RESULTS: Similar to neuroblastoma, TERT expression exhibited a bimodal distribution in pulmonary carcinoids, separating tumors into TERT-high and TERT-low subgroups. A natural TERT cutoff discriminated unfavorable from favorable clinical courses with high accuracy both in the test cohort (5-year overall survival [OS], 0.547 ± 0.132 v 1.0; P < .001) and the validation cohort (5-year OS, 0.788 ± 0.063 v 0.913 ± 0.048; P < .001). In line with these findings, telomerase activity was largely absent in TERT-low tumors, whereas it was readily detectable in TERT-high carcinoids. In multivariable analysis considering TERT expression, histology (typical v atypical carcinoid), and stage (≤IIA v ≥IIB), high TERT expression was an independent prognostic marker for poor survival, with a hazard ratio of 5.243 (95% CI, 1.943 to 14.148; P = .001). CONCLUSION: Our data demonstrate that high TERT expression defines clinically aggressive pulmonary carcinoids with fatal outcome, similar to neuroblastoma, indicating that activation of TERT may be a defining feature of lethal cancers.

PRDX6 contributes to selenocysteine metabolism and ferroptosis resistance.

Selenocysteine (Sec) metabolism is crucial for cellular function and ferroptosis prevention and has traditionally been thought to begin with the uptake of the Sec carrier selenoprotein P (SELENOP). Following uptake, Sec released from SELENOP undergoes metabolisation via selenocysteine lyase (SCLY), producing selenide, a substrate used by selenophosphate synthetase 2 (SEPHS2), which provides the essential selenium donor - selenophosphate - for the biosynthesis of the selenocysteine tRNA. Here, we report the discovery of an alternative pathway mediating Sec metabolisation that is independent of SCLY and mediated by peroxiredoxin 6 (PRDX6). Mechanistically, we demonstrate that PRDX6 can readily react with selenide and interact with SEPHS2, potentially acting as a selenium delivery system. Moreover, we demonstrate the presence and functional significance of this alternative route in cancer cells where we reveal a notable association between elevated expression of PRDX6 with a highly aggressive neuroblastoma subtype. Altogether, our study sheds light on a previously unrecognized aspect of Sec metabolism and its implications in ferroptosis, offering new avenues for therapeutic exploitation.

METTL3 drives telomere targeting of TERRA lncRNA through m6A-dependent R-loop formation: a therapeutic target for ALT-positive neuroblastoma.

Telomerase-negative tumors maintain telomere length by alternative lengthening of telomeres (ALT), but the underlying mechanism behind ALT remains poorly understood. A proportion of aggressive neuroblastoma (NB), particularly relapsed tumors, are positive for ALT (ALT+), suggesting that a better dissection of the ALT mechanism could lead to novel therapeutic opportunities. TERRA, a long non-coding RNA (lncRNA) derived from telomere ends, localizes to telomeres in a R-loop-dependent manner and plays a crucial role in telomere maintenance. Here we present evidence that RNA modification at the N6 position of internal adenosine (m6A) in TERRA by the methyltransferase METTL3 is essential for telomere maintenance in ALT+ cells, and the loss of TERRA m6A/METTL3 results in telomere damage. We observed that m6A modification is abundant in R-loop enriched TERRA, and the m6A-mediated recruitment of hnRNPA2B1 to TERRA is critical for R-loop formation. Our findings suggest that m6A drives telomere targeting of TERRA via R-loops, and this m6A-mediated R-loop formation could be a widespread mechanism employed by other chromatin-interacting lncRNAs. Furthermore, treatment of ALT+ NB cells with a METTL3 inhibitor resulted in compromised telomere targeting of TERRA and accumulation of DNA damage at telomeres, indicating that METTL3 inhibition may represent a therapeutic approach for ALT+ NB.

Somatic rearrangements causing oncogenic ectodomain deletions of FGFR1 in squamous cell lung cancer.

The discovery of frequent 8p11-p12 amplifications in squamous cell lung cancer (SQLC) has fueled hopes that FGFR1, located inside this amplicon, might be a therapeutic target. In a clinical trial, only 11% of patients with 8p11 amplification (detected by FISH) responded to FGFR kinase inhibitor treatment. To understand the mechanism of FGFR1 dependency, we performed deep genomic characterization of 52 SQLCs with 8p11-p12 amplification, including 10 tumors obtained from patients who had been treated with FGFR inhibitors. We discovered somatically altered variants of FGFR1 with deletion of exons 1-8 that resulted from intragenic tail-to-tail rearrangements. These ectodomain-deficient FGFR1 variants (ΔEC-FGFR1) were expressed in the affected tumors and were tumorigenic in both in vitro and in vivo models of lung cancer. Mechanistically, breakage-fusion-bridges were the source of 8p11-p12 amplification, resulting from frequent head-to-head and tail-to-tail rearrangements. Generally, tail-to-tail rearrangements within or in close proximity upstream of FGFR1 were associated with FGFR1 dependency. Thus, the genomic events shaping the architecture of the 8p11-p12 amplicon provide a mechanistic explanation for the emergence of FGFR1-driven SQLC. Specifically, we believe that FGFR1 ectodomain-deficient and FGFR1-centered amplifications caused by tail-to-tail rearrangements are a novel somatic genomic event that might be predictive of therapeutically relevant FGFR1 dependency.

LRP8-mediated selenocysteine uptake is a targetable vulnerability in MYCN-amplified neuroblastoma.

Ferroptosis has emerged as an attractive strategy in cancer therapy. Understanding the operational networks regulating ferroptosis may unravel vulnerabilities that could be harnessed for therapeutic benefit. Using CRISPR-activation screens in ferroptosis hypersensitive cells, we identify the selenoprotein P (SELENOP) receptor, LRP8, as a key determinant protecting MYCN-amplified neuroblastoma cells from ferroptosis. Genetic deletion of LRP8 leads to ferroptosis as a result of an insufficient supply of selenocysteine, which is required for the translation of the antiferroptotic selenoprotein GPX4. This dependency is caused by low expression of alternative selenium uptake pathways such as system Xc- . The identification of LRP8 as a specific vulnerability of MYCN-amplified neuroblastoma cells was confirmed in constitutive and inducible LRP8 knockout orthotopic xenografts. These findings disclose a yet-unaccounted mechanism of selective ferroptosis induction that might be explored as a therapeutic strategy for high-risk neuroblastoma and potentially other MYCN-amplified entities.

Genomic ALK alterations in primary and relapsed neuroblastoma.

BACKGROUND: Genomic alterations of the anaplastic lymphoma kinase gene (ALK) occur recurrently in neuroblastoma, a pediatric malignancy of the sympathetic nervous system. However, information on their development over time has remained sparse. METHODS: ALK alterations were assessed in neuroblastomas at diagnosis and/or relapse from a total of 943 patients, covering all stages of disease. Longitudinal information on diagnostic and relapsed samples from individual patients was available in 101 and 102 cases for mutation and amplification status, respectively. RESULTS: At diagnosis, ALK point mutations occurred in 10.5% of all cases, with highest frequencies in stage 4 patients <18 months. At relapse, ALK alteration frequency increased by 70%, both in high-risk and non-high-risk cases. The increase was most likely due to de novo mutations, frequently leading to R1275Q substitutions, which are sensitive to pharmacological ALK inhibition. By contrast, the frequency of ALK amplifications did not change over the course of the disease. ALK amplifications, but not mutations, were associated with poor patient outcome. CONCLUSIONS: The considerably increased frequency of ALK mutations at relapse and their high prevalence in young stage 4 patients suggest surveying the genomic ALK status regularly in these patient cohorts, and to evaluate ALK-targeted treatment also in intermediate-risk patients.

Reliable assessment of telomere maintenance mechanisms in neuroblastoma.

BACKGROUND: Telomere maintenance mechanisms (TMM) are a hallmark of high-risk neuroblastoma, and are conferred by activation of telomerase or alternative lengthening of telomeres (ALT). However, detection of TMM is not yet part of the clinical routine, and consensus on TMM detection, especially on ALT assessment, remains to be achieved. METHODS: Whole genome sequencing (WGS) data of 68 primary neuroblastoma samples were analyzed. Telomere length was calculated from WGS data or by telomere restriction fragment analysis (n = 39). ALT was assessed by C-circle assay (CCA, n = 67) and detection of ALT-associated PML nuclear bodies (APB) by combined fluorescence in situ hybridization and immunofluorescence staining (n = 68). RNA sequencing was performed (n = 64) to determine expression of TERT and telomeric long non-coding RNA (TERRA). Telomerase activity was examined by telomerase repeat amplification protocol (TRAP, n = 15). RESULTS: Tumors were considered as telomerase-positive if they harbored a TERT rearrangement, MYCN amplification or high TERT expression (45.6%, 31/68), and ALT-positive if they were positive for APB and CCA (19.1%, 13/68). If all these markers were absent, tumors were considered TMM-negative (25.0%, 17/68). According to these criteria, the majority of samples were classified unambiguously (89.7%, 61/68). Assessment of additional ALT-associated parameters clarified the TMM status of the remaining seven cases with high likelihood: ALT-positive tumors had higher TERRA expression, longer telomeres, more telomere insertions, a characteristic pattern of telomere variant repeats, and were associated with ATRX mutations. CONCLUSIONS: We here propose a workflow to reliably detect TMM in neuroblastoma. We show that unambiguous classification is feasible following a stepwise approach that determines both, activation of telomerase and ALT. The workflow proposed in this study can be used in clinical routine and provides a framework to systematically and reliably determine telomere maintenance mechanisms for risk stratification and treatment allocation of neuroblastoma patients.

Telomerase-targeting compounds Imetelstat and 6-thio-dG act synergistically with chemotherapy in high-risk neuroblastoma models.

BACKGROUND: The majority of high-risk neuroblastomas harbor telomerase activity, and telomerase-interacting compounds, such as 6-thio-2'-deoxyguanosine (6-thio-dG), have been found to impair the growth of telomerase-positive neuroblastoma cell lines. It has remained unclear, however, how such drugs can be combined with other compounds used in current treatment concepts for neuroblastoma patients. METHODS: Growth-inhibitory effects of varying concentrations of 6-thio-dG in combination with etoposide, doxorubicin or ceritinib were determined in eight telomerase-positive neuroblastoma cell lines with distinct genetic backgrounds. Tumor growth inhibition of subcutaneous xenografts from three different cell lines was assessed upon treatment with 6-thio-dG, the competitive telomerase inhibitor imetelstat, etoposide, or combinations of these compounds. RESULTS: Robust synergistic anti-tumor effects were observed for combinations of 6-thio-dG and etoposide or doxorubicin, but not for 6-thio-dG and ceritinib, in telomerase-positive neuroblastoma cell lines in vitro. Treatment of mouse xenografts with combinations of 6-thio-dG and etoposide significantly attenuated tumor growth and improved mouse survival over etoposide alone in two of three cell line models. Treatment of xenograft tumors by imetelstat monotherapy decreased telomerase activity by roughly 50% and significantly improved survival over control in all three models, whereas treatment with imetelstat plus etoposide led to enhanced survival over etoposide monotherapy in one model. Mechanistically, the synergistic effect was found to be due to both increased apoptosis and cell cycle arrest. CONCLUSION: Our study indicates that telomerase is an actionable target in telomerase-positive neuroblastoma, and demonstrates that combination therapies including telomerase-interacting compounds may improve the efficacy of established cytotoxic drugs. Targeting telomerase may thus represent a therapeutic option in high-risk neuroblastoma patients.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition.

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

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.

Genetic Alterations and Resectability Predict Outcome in Patients with Neuroblastoma Assigned to High-Risk Solely by MYCN Amplification.

BACKGROUND: To identify variables predicting outcome in neuroblastoma patients assigned to the high-risk group solely by the presence of MYCN oncogene amplification (MNA). METHODS: Clinical characteristics, genomic information, and outcome of 190 patients solely assigned to high-risk neuroblastoma by MNA were analyzed and compared to 205 patients with stage 4 neuroblastoma aged ≥18 months with MNA (control group). RESULTS: Event-free survival (EFS) and overall survival (OS) at 10 years were 47% (95%-CI 39-54%) and 56% (95%-CI 49-63%), respectively, which was significantly better than EFS and OS of the control group (EFS 25%, 95%-CI 18-31%, p < 0.001; OS 32% 95%-CI 25-39%, p < 0.001). The presence of RAS-/p53-pathway gene alterations was associated with impaired 10-year EFS and OS (19% vs. 55%, and 19% vs. 67%, respectively; both p < 0.001). In time-dependent multivariable analyses, alterations of RAS-/p53-pathway genes and the extent of the best primary tumor resection were the only independent prognostic variables for OS (p < 0.001 and p = 0.011, respectively). CONCLUSIONS: Neuroblastoma patients attributed to high risk solely by MYCN amplification have generally a more favorable outcome. Mutations of genes of the RAS and/or p53 pathways and incomplete resection are the main risk factors predicting poor outcome.

Defective Interfering Genomes and the Full-Length Viral Genome Trigger RIG-I After Infection With Vesicular Stomatitis Virus in a Replication Dependent Manner.

Replication competent vesicular stomatitis virus (VSV) is the basis of a vaccine against Ebola and VSV strains are developed as oncolytic viruses. Both functions depend on the ability of VSV to induce adequate amounts of interferon-α/ß. It is therefore important to understand how VSV triggers interferon responses. VSV activates innate immunity via retinoic acid-inducible gene I (RIG-I), a sensor for viral RNA. Our results show that VSV needs to replicate for a robust interferon response. Analysis of RIG-I-associated RNA identified a copy-back defective-interfering (DI) genome and full-length viral genomes as main trigger of RIG-I. VSV stocks depleted of DI genomes lost most of their interferon-stimulating activity. The remaining full-length genome and leader-N-read-through sequences, however, still triggered RIG-I. Awareness for DI genomes as trigger of innate immune responses will help to standardize DI genome content and to purposefully deplete or use DI genomes as natural adjuvants in VSV-based therapeutics.

Comprehensive germline-genomic and clinical profiling in 160 unselected children and adolescents with cancer.

In childhood cancer, the frequency of cancer-associated germline variants and their inheritance patterns are not thoroughly investigated. Moreover, the identification of children carrying a genetic predisposition by clinical means remains challenging. In this single-center study, we performed trio whole-exome sequencing and comprehensive clinical evaluation of a prospectively enrolled cohort of 160 children with cancer and their parents. We identified in 11/160 patients a pathogenic germline variant predisposing to cancer and a further eleven patients carried a prioritized VUS with a strong association to the cancerogenesis of the patient. Through clinical screening, 51 patients (31.3%) were identified as suspicious for an underlying cancer predisposition syndrome (CPS), but only in ten of those patients a pathogenic variant could be identified. In contrast, one patient with a classical CPS and ten patients with prioritized VUS were classified as unremarkable in the clinical work-up. Taken together, a monogenetic causative variant was detected in 13.8% of our patients using WES. Nevertheless, the still unclarified clinical suspicious cases emphasize the need to consider other genetic mechanisms including new target genes, structural variants, or polygenic interactions not previously associated with cancer predisposition.

Integrative genomic analysis of pediatric T-cell lymphoblastic lymphoma reveals candidates of clinical significance.

T-cell lymphoblastic lymphoma (T-LBL) is a heterogeneous malignancy of lymphoblasts committed to T-cell lineage. The dismal outcomes (15%-30%) after T-LBL relapse warrant establishing risk-based treatment. To our knowledge, this study presents the first comprehensive, systematic, integrated, genome-wide analysis including relapsed cases that identifies molecular markers of prognostic relevance for T-LBL. NOTCH1 was identified as the putative driver for T-LBL. An activated NOTCH/PI3K-AKT signaling axis and alterations in cell cycle regulators constitute the core oncogenic program for T-LBL. Mutated KMT2D was identified as a prognostic marker. The cumulative incidence of relapse was 47% ± 17% in patients with KMT2D mutations, compared with 14% ± 3% in wild-type KMT2D. Structural analysis of the mutated domains of KMT2D revealed a plausible impact on structure and functional consequences. These findings provide new insights into the pathogenesis of T-LBL, including high translational potential. The ongoing LBL 2018 trial (www.clinicaltrials.gov #NCT04043494) allows for prospective validation and subsequent fine tuning of the stratification criteria for T-LBL risk groups to improve survival of pediatric patients.

Targeted Therapy of TERT-Rearranged Neuroblastoma with BET Bromodomain Inhibitor and Proteasome Inhibitor Combination Therapy.

PURPOSE: TERT gene rearrangement with transcriptional superenhancers leads to TERT overexpression and neuroblastoma. No targeted therapy is available for clinical trials in patients with TERT-rearranged neuroblastoma. EXPERIMENTAL DESIGN: Anticancer agents exerting the best synergistic anticancer effects with BET bromodomain inhibitors were identified by screening an FDA-approved oncology drug library. The synergistic effects of the BET bromodomain inhibitor OTX015 and the proteasome inhibitor carfilzomib were examined by immunoblot and flow cytometry analysis. The anticancer efficacy of OTX015 and carfilzomib combination therapy was investigated in mice xenografted with TERT-rearranged neuroblastoma cell lines or patient-derived xenograft (PDX) tumor cells, and the role of TERT reduction in the anticancer efficacy was examined through rescue experiments in mice. RESULTS: The BET bromodomain protein BRD4 promoted TERT-rearranged neuroblastoma cell proliferation through upregulating TERT expression. Screening of an approved oncology drug library identified the proteasome inhibitor carfilzomib as the agent exerting the best synergistic anticancer effects with BET bromodomain inhibitors including OTX015. OTX015 and carfilzomib synergistically reduced TERT protein expression, induced endoplasmic reticulum stress, and induced TERT-rearranged neuroblastoma cell apoptosis which was blocked by TERT overexpression and endoplasmic reticulum stress antagonists. In mice xenografted with TERT-rearranged neuroblastoma cell lines or PDX tumor cells, OTX015 and carfilzomib synergistically blocked TERT expression, induced tumor cell apoptosis, suppressed tumor progression, and improved mouse survival, which was largely reversed by forced TERT overexpression. CONCLUSIONS: OTX015 and carfilzomib combination therapy is likely to be translated into the first clinical trial of a targeted therapy in patients with TERT-rearranged neuroblastoma.

The ETS transcription factor ETV5 is a target of activated ALK in neuroblastoma contributing to increased tumour aggressiveness.

Neuroblastoma is an aggressive childhood cancer arising from sympatho-adrenergic neuronal progenitors. The low survival rates for high-risk disease point to an urgent need for novel targeted therapeutic approaches. Detailed molecular characterization of the neuroblastoma genomic landscape indicates that ALK-activating mutations are present in 10% of primary tumours. Together with other mutations causing RAS/MAPK pathway activation, ALK mutations are also enriched in relapsed cases and ALK activation was shown to accelerate MYCN-driven tumour formation through hitherto unknown ALK-driven target genes. To gain further insight into how ALK contributes to neuroblastoma aggressiveness, we searched for known oncogenes in our previously reported ALK-driven gene signature. We identified ETV5, a bona fide oncogene in prostate cancer, as robustly upregulated in neuroblastoma cells harbouring ALK mutations, and show high ETV5 levels downstream of the RAS/MAPK axis. Increased ETV5 expression significantly impacted migration, invasion and colony formation in vitro, and ETV5 knockdown reduced proliferation in a murine xenograft model. We also established a gene signature associated with ETV5 knockdown that correlates with poor patient survival. Taken together, our data highlight ETV5 as an intrinsic component of oncogenic ALK-driven signalling through the MAPK axis and propose that ETV5 upregulation in neuroblastoma may contribute to tumour aggressiveness.

Infectious stimuli promote malignant B-cell acute lymphoblastic leukemia in the absence of AID.

The prerequisite to prevent childhood B-cell acute lymphoblastic leukemia (B-ALL) is to decipher its etiology. The current model suggests that infection triggers B-ALL development through induction of activation-induced cytidine deaminase (AID; also known as AICDA) in precursor B-cells. This evidence has been largely acquired through the use of ex vivo functional studies. However, whether this mechanism governs native non-transplant B-ALL development is unknown. Here we show that, surprisingly, AID genetic deletion does not affect B-ALL development in Pax5-haploinsufficient mice prone to B-ALL upon natural infection exposure. We next test the effect of premature AID expression from earliest pro-B-cell stages in B-cell transformation. The generation of AID off-target mutagenic activity in precursor B-cells does not promote B-ALL. Likewise, known drivers of human B-ALL are not preferentially targeted by AID. Overall these results suggest that infections promote B-ALL through AID-independent mechanisms, providing evidence for a new model of childhood B-ALL development.

Author Correction: Integrative analysis identifies lincRNAs up- and downstream of neuroblastoma driver genes.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

JMJD6 is a tumorigenic factor and therapeutic target in neuroblastoma.

Chromosome 17q21-ter is commonly gained in neuroblastoma, but it is unclear which gene in the region is important for tumorigenesis. The JMJD6 gene at 17q21-ter activates gene transcription. Here we show that JMJD6 forms protein complexes with N-Myc and BRD4, and is important for E2F2, N-Myc and c-Myc transcription. Knocking down JMJD6 reduces neuroblastoma cell proliferation and survival in vitro and tumor progression in mice, and high levels of JMJD6 expression in human neuroblastoma tissues independently predict poor patient prognosis. In addition, JMJD6 gene is associated with transcriptional super-enhancers. Combination therapy with the CDK7/super-enhancer inhibitor THZ1 and the histone deacetylase inhibitor panobinostat synergistically reduces JMJD6, E2F2, N-Myc, c-Myc expression, induces apoptosis in vitro and leads to neuroblastoma tumor regression in mice, which are significantly reversed by forced JMJD6 over-expression. Our findings therefore identify JMJD6 as a neuroblastoma tumorigenesis factor, and the combination therapy as a treatment strategy.

Integrative analysis identifies lincRNAs up- and downstream of neuroblastoma driver genes.

Long intergenic non-coding RNAs (lincRNAs) are emerging as integral components of signaling pathways in various cancer types. In neuroblastoma, only a handful of lincRNAs are known as upstream regulators or downstream effectors of oncogenes. Here, we exploit RNA sequencing data of primary neuroblastoma tumors, neuroblast precursor cells, neuroblastoma cell lines and various cellular perturbation model systems to define the neuroblastoma lincRNome and map lincRNAs up- and downstream of neuroblastoma driver genes MYCN, ALK and PHOX2B. Each of these driver genes controls the expression of a particular subset of lincRNAs, several of which are associated with poor survival and are differentially expressed in neuroblastoma tumors compared to neuroblasts. By integrating RNA sequencing data from both primary tumor tissue and cancer cell lines, we demonstrate that several of these lincRNAs are expressed in stromal cells. Deconvolution of primary tumor gene expression data revealed a strong association between stromal cell composition and driver gene status, resulting in differential expression of these lincRNAs. We also explored lincRNAs that putatively act upstream of neuroblastoma driver genes, either as presumed modulators of driver gene activity, or as modulators of effectors regulating driver gene expression. This analysis revealed strong associations between the neuroblastoma lincRNAs MIAT and MEG3 and MYCN and PHOX2B activity or expression. Together, our results provide a comprehensive catalogue of the neuroblastoma lincRNome, highlighting lincRNAs up- and downstream of key neuroblastoma driver genes. This catalogue forms a solid basis for further functional validation of candidate neuroblastoma lincRNAs.