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.

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.

Ligand binding modes from low resolution GPCR models and mutagenesis: chicken bitter taste receptor as a test-case.

Bitter taste is one of the basic taste modalities, warning against consuming potential poisons. Bitter compounds activate members of the bitter taste receptor (Tas2r) subfamily of G protein-coupled receptors (GPCRs). The number of functional Tas2rs is species-dependent. Chickens represent an intriguing minimalistic model, because they detect the bitter taste of structurally different molecules with merely three bitter taste receptor subtypes. We investigated the binding modes of several known agonists of a representative chicken bitter taste receptor, ggTas2r1. Because of low sequence similarity between ggTas2r1 and crystallized GPCRs (~10% identity, ~30% similarity at most), the combination of computational approaches with site-directed mutagenesis was used to characterize the agonist-bound conformation of ggTas2r1 binding site between TMs 3, 5, 6 and 7. We found that the ligand interactions with N93 in TM3 and/or N247 in TM5, combined with hydrophobic contacts, are typically involved in agonist recognition. Next, the ggTas2r1 structural model was successfully used to identify three quinine analogues (epiquinidine, ethylhydrocupreine, quinidine) as new ggTas2r1 agonists. The integrated approach validated here may be applicable to additional cases where the sequence identity of the GPCR of interest and the existing experimental structures is low.