The human SMAD9 (GCC) repeat links to natural selection and late-onset neurocognitive disorders.

INTRODUCTION: Whereas (GCC)-repeats are overrepresented in genic regions, and mutation hotspots, they are largely unexplored with regard to their link with natural selection. Across numerous primate species and tissues, SMAD9 (SMAD Family Member 9) reaches highest level of expression in the human brain. This gene contains a (GCC)-repeat in the interval between + 1 and + 60 of the transcription start site, which is in the high-ranking (GCC)-repeats with respect to length. METHODS: Here we sequenced this (GCC)-repeat in 396 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 181) and controls (N = 215). RESULTS: We detected two predominantly abundant alleles of 7 and 9 repeats, forming 96.2% of the allele pool. The (GCC)7/(GCC)9 ratio was in the reverse order in the NCD group versus controls (p = 0.005), resulting from excess of (GCC)7 in the NCD group (p = 0.003) and (GCC)9 in the controls (p = 0.01). Five genotypes, predominantly consisting of (GCC)7 and lacking (GCC)9 were detected in the NCD group only (p = 0.008). The patients harboring those genotypes received the diagnoses of Alzheimer's disease (AD) and vascular dementia (VD). Five genotypes consisting of (GCC)9 and lacking (GCC)7 were detected in the control group only (p = 0.002). The group-specific genotypes formed approximately 4% of the genotype pool in the human samples studied. CONCLUSION: We propose natural selection and a novel locus for late-onset AD and VD at the SMAD9 (GCC)-repeat in humans.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

Transcription factor RFX7 governs a tumor suppressor network in response to p53 and stress.

Despite its prominence, the mechanisms through which the tumor suppressor p53 regulates most genes remain unclear. Recently, the regulatory factor X 7 (RFX7) emerged as a suppressor of lymphoid neoplasms, but its regulation and target genes mediating tumor suppression remain unknown. Here, we identify a novel p53-RFX7 signaling axis. Integrative analysis of the RFX7 DNA binding landscape and the RFX7-regulated transcriptome in three distinct cell systems reveals that RFX7 directly controls multiple established tumor suppressors, including PDCD4, PIK3IP1, MXD4, and PNRC1, across cell types and is the missing link for their activation in response to p53 and stress. RFX7 target gene expression correlates with cell differentiation and better prognosis in numerous cancer types. Interestingly, we find that RFX7 sensitizes cells to Doxorubicin by promoting apoptosis. Together, our work establishes RFX7's role as a ubiquitous regulator of cell growth and fate determination and a key node in the p53 transcriptional program.

T-cell prolymphocytic leukemia is associated with deregulation of oncogenic microRNAs on transcriptional and epigenetic level.

Deregulation of micro(mi)-RNAs is a common mechanism in tumorigenesis. We investigated the expression of 2083 miRNAs in T-cell prolymphocytic leukemia (T-PLL). Compared to physiologic CD4+ and CD8+ T-cell subsets, 111 miRNAs were differentially expressed in T-PLL. Of these, 33 belonged to miRNA gene clusters linked to cancer. Genomic variants affecting miRNAs were infrequent with the notable exception of copy number aberrations. Remarkably, we found strong upregulation of the miR-200c/-141 cluster in T-PLL to be associated with DNA hypomethylation and active promoter marks. Our findings suggest that copy number aberrations and epigenetic changes could contribute to miRNA deregulation in T-PLL.

Multiomics reveal unique signatures of human epiploic adipose tissue related to systemic insulin resistance.

OBJECTIVE: Human white adipose tissue (AT) is a metabolically active organ with distinct depot-specific functions. Despite their locations close to the gastrointestinal tract, mesenteric AT and epiploic AT (epiAT) have only scarcely been investigated. Here, we aim to characterise these ATs in-depth and estimate their contribution to alterations in whole-body metabolism. DESIGN: Mesenteric, epiploic, omental and abdominal subcutaneous ATs were collected from 70 patients with obesity undergoing Roux-en-Y gastric bypass surgery. The metabolically well-characterised cohort included nine subjects with insulin sensitive (IS) obesity, whose AT samples were analysed in a multiomics approach, including methylome, transcriptome and proteome along with samples from subjects with insulin resistance (IR) matched for age, sex and body mass index (n=9). Findings implying differences between AT depots in these subgroups were validated in the entire cohort (n=70) by quantitative real-time PCR. RESULTS: While mesenteric AT exhibited signatures similar to those found in the omental depot, epiAT was distinct from all other studied fat depots. Multiomics allowed clear discrimination between the IS and IR states in all tissues. The highest discriminatory power between IS and IR was seen in epiAT, where profound differences in the regulation of developmental, metabolic and inflammatory pathways were observed. Gene expression levels of key molecules involved in AT function, metabolic homeostasis and inflammation revealed significant depot-specific differences with epiAT showing the highest expression levels. CONCLUSION: Multi-omics epiAT signatures reflect systemic IR and obesity subphenotypes distinct from other fat depots. Our data suggest a previously unrecognised role of human epiploic fat in the context of obesity, impaired insulin sensitivity and related diseases.

Effects of lifestyle interventions on epigenetic signatures of liver fat: Central randomized controlled trial.

BACKGROUND AND AIMS: In the CENTRAL trial context, we found diverse liver fat dynamics in response to different dietary interventions. Epigenetic mechanisms may contribute to the intraindividual variation. Moreover, genetic factors are involved in developing nonalcoholic fatty-liver disease (NAFLD), a disease reflected by an increase in intrahepatic fat (IHF). In this exploratory analysis, we primarily aimed to examine the effect of lifestyle interventions on DNA-methylation of NAFLD related genes associated with IHF. METHODS: For 120 participants from the CENTRAL trial, an 18-month regimen of either low-fat (LF) or Mediterranean-low carbohydrate (MED/LC) diets, with or without physical activity (PA+/PA-), was instructed. Magnetic resonance imaging was used to measure IHF%, which was analysed for association with CpG specific DNA-methylation levels of 41 selected candidate genes. Single-nucleotide polymorphisms known to be associated with NAFLD within the studied genes were genotyped by TaqMan assays. RESULTS: At baseline, participants (92% men; body mass index = 30.2 kg/m2 ) had mean IHF of 10.7% (59% NAFLD). Baseline-IHF% was inversely correlated with DNA-methylation at individual CpGs within AC074286.1, CRACR2A, A2MP1, FARP1 (P < .05 for all multivariate models). FARP1 rs9584805 showed association with IHF, with the prevalence of NAFLD and baseline methylation level of the CpG site (cg00071727) associated with IHF%. Following 18-month lifestyle intervention, differential DNA-methylation patterns were observed between diets at cg14335324 annotated to A2MP1 (P = .04, LF vs. MED/LC), and differential DNA-methylation between PA groups within AC074286.1, CRACR2A, and FARP1 CpGs (P < .05 for all, PA-vs. PA+). CONCLUSIONS: This study suggests epigenetic markers for IHF and potential epigenetic remodeling after long-term lifestyle interventions.

Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes.

With the rapid increase of sequenced metazoan mitochondrial genomes, a detailed manual annotation is becoming more and more infeasible. While it is easy to identify the approximate location of protein-coding genes within mitogenomes, the peculiar processing of mitochondrial transcripts, however, makes the determination of precise gene boundaries a surprisingly difficult problem. We have analyzed the properties of annotated start and stop codon positions in detail, and use the inferred patterns to devise a new method for predicting gene boundaries in de novo annotations. Our method benefits from empirically observed prevalances of start/stop codons and gene lengths, and considers the dependence of these features on variations of genetic codes. Albeit not being perfect, our new approach yields a drastic improvement in the accuracy of gene boundaries and upgrades the mitochondrial genome annotation server MITOS to an even more sophisticated tool for fully automatic annotation of metazoan mitochondrial genomes.

DREAM and RB cooperate to induce gene repression and cell-cycle arrest in response to p53 activation.

Most human cancers acquire mutations causing defects in the p53 signaling pathway. The tumor suppressor p53 becomes activated in response to genotoxic stress and is essential for arresting the cell cycle to facilitate DNA repair or to initiate apoptosis. p53-induced cell cycle-arrest is mediated by expression of the CDK inhibitor p21WAF1/Cip1, which prevents phosphorylation and inactivation of the pocket proteins RB, p130, and p107. In a hypophosphorylated state, pocket proteins bind to E2F factors forming RB-E2F and DREAM transcriptional repressor complexes. Here, we analyze the influence of RB and DREAM on p53-induced gene repression and cell-cycle arrest. We show that abrogation of DREAM function by knockout of the DREAM component LIN37 results in a reduced repression of cell-cycle genes. We identify the genes repressed by the p53-DREAM pathway and describe a set of genes that is downregulated by p53 independent of LIN37/DREAM. Most strikingly, p53-dependent repression of cell-cycle genes is completely abrogated in LIN37-/-;RB-/- cells leading to a loss of the G1/S checkpoint. Taken together, we show that DREAM and RB are key factors in the p53 signaling pathway to downregulate a large number of cell-cycle genes and to arrest the cell cycle at the G1/S transition.

Reconstruction of rearranged T-cell receptor loci by whole genome and transcriptome sequencing gives insights into the initial steps of T-cell prolymphocytic leukemia.

T-cell prolymphocytic leukemia (T-PLL) is an aggressive tumor with leukemic presentation of mature T-lymphocytes. Here, we aimed at characterizing the initial events in the molecular pathogenesis of T-PLL and particularly, at determining the point in T-cell differentiation when the hallmark oncogenic events, that is, inv(14)(q11q32)/t(14;14)(q11;q32) and t(X;14)(q28;q11) occur. To this end, we mined whole genome and transcriptome sequencing data of 17 and 11 T-PLL cases, respectively. Mapping of the 14q32.1 locus breakpoints identified only TCL1A, which was moreover significantly overexpressed in T-PLL as compared to benign CD4+ and CD8+ T-cells, as the only common oncogenic target of aberrations. In cases with t(14;14), the breakpoints mapped telomeric and in cases with inv(14) centromeric or in the 3'-untranslated region of TCL1A. Regarding the T-cell receptor alpha (TRA) locus-TCL1A breakpoint junctions, all 17 breakpoints involved recombination signal sequences and 15 junctions contained nontemplated (N-) nucleotides. All T-PLL cases studied carried in-frame TRA rearrangements on the intact allele, which skewed significantly toward usage of distal/central TRAV/TRAJ gene segments as compared to the illegitimate TRA rearrangements. Our findings suggest that the oncogenic TRA-TCL1A/MTCP1 rearrangements in T-PLL occur during opening of the TRA locus, that is, during the progression from CD4+ immature single positive to early double positive thymocyte stage, just before physiologic TCL1A expression is silenced. The cell carrying such an oncogenic event continues maturation and rearranges the second TRA allele to achieve a functional T-cell receptor. Thereafter, it switches off RAG and DNTT expression in line with the mature T-cell phenotype at presentation of T-PLL.

Identification of RNA 3´ ends and termination sites in Haloferax volcanii.

Archaeal genomes are densely packed; thus, correct transcription termination is an important factor for orchestrated gene expression. A systematic analysis of RNA 3´ termini, to identify transcription termination sites (TTS) using RNAseq data has hitherto only been performed in two archaea, Methanosarcina mazei and Sulfolobus acidocaldarius. In this study, only regions directly downstream of annotated genes were analysed, and thus, only part of the genome had been investigated. Here, we developed a novel algorithm (Internal Enrichment-Peak Calling) that allows an unbiased, genome-wide identification of RNA 3´ termini independent of annotation. In an RNA fraction enriched for primary transcripts by terminator exonuclease (TEX) treatment we identified 1,543 RNA 3´ termini. Approximately half of these were located in intergenic regions, and the remainder were found in coding regions. A strong sequence signature consistent with known termination events at intergenic loci indicates a clear enrichment for native TTS among them. Using these data we determined distinct putative termination motifs for intergenic (a T stretch) and coding regions (AGATC). In vivo reporter gene tests of selected TTS confirmed termination at these sites, which exemplify the different motifs. For several genes, more than one termination site was detected, resulting in transcripts with different lengths of the 3´ untranslated region (3´ UTR).

metilene: fast and sensitive calling of differentially methylated regions from bisulfite sequencing data.

The detection of differentially methylated regions (DMRs) is a necessary prerequisite for characterizing different epigenetic states. We present a novel program, metilene, to identify DMRs within whole-genome and targeted data with unrivaled specificity and sensitivity. A binary segmentation algorithm combined with a two-dimensional statistical test allows the detection of DMRs in large methylation experiments with multiple groups of samples in minutes rather than days using off-the-shelf hardware. metilene outperforms other state-of-the-art tools for low coverage data and can estimate missing data. Hence, metilene is a versatile tool to study the effect of epigenetic modifications in differentiation/development, tumorigenesis, and systems biology on a global, genome-wide level. Whether in the framework of international consortia with dozens of samples per group, or even without biological replicates, it produces highly significant and reliable results.

DIEGO: detection of differential alternative splicing using Aitchison's geometry.

Motivation: Alternative splicing is a biological process of fundamental importance in most eukaryotes. It plays a pivotal role in cell differentiation and gene regulation and has been associated with a number of different diseases. The widespread availability of RNA-Sequencing capacities allows an ever closer investigation of differentially expressed isoforms. However, most tools for differential alternative splicing (DAS) analysis do not take split reads, i.e. the most direct evidence for a splice event, into account. Here, we present DIEGO, a compositional data analysis method able to detect DAS between two sets of RNA-Seq samples based on split reads. Results: The python tool DIEGO works without isoform annotations and is fast enough to analyze large experiments while being robust and accurate. We provide python and perl parsers for common formats. Availability and implementation: The software is available at: www.bioinf.uni-leipzig.de/Software/DIEGO. Contact: steve@bioinf.uni-leipzig.de. Supplementary information: Supplementary data are available at Bioinformatics online.

Temperature-responsive in vitro RNA structurome of Yersinia pseudotuberculosis.

RNA structures are fundamentally important for RNA function. Dynamic, condition-dependent structural changes are able to modulate gene expression as shown for riboswitches and RNA thermometers. By parallel analysis of RNA structures, we mapped the RNA structurome of Yersinia pseudotuberculosis at three different temperatures. This human pathogen is exquisitely responsive to host body temperature (37 °C), which induces a major metabolic transition. Our analysis profiles the structure of more than 1,750 RNAs at 25 °C, 37 °C, and 42 °C. Average mRNAs tend to be unstructured around the ribosome binding site. We searched for 5'-UTRs that are folded at low temperature and identified novel thermoresponsive RNA structures from diverse gene categories. The regulatory potential of 16 candidates was validated. In summary, we present a dynamic bacterial RNA structurome and find that the expression of virulence-relevant functions in Y. pseudotuberculosis and reprogramming of its metabolism in response to temperature is associated with a restructuring of numerous mRNAs.

MINCR is a MYC-induced lncRNA able to modulate MYC's transcriptional network in Burkitt lymphoma cells.

Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long noncoding RNAs (lncRNAs), on MYC ability to influence the cellular transcriptome. Here, we have intersected RNA-sequencing data from two MYC-inducible cell lines and a cohort of 91 B-cell lymphomas with or without genetic variants resulting in MYC overexpression. We identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them, we focused on a lncRNA that we named MYC-induced long noncoding RNA (MINCR), showing a strong correlation with MYC expression in MYC-positive lymphomas. To understand its cellular role, we performed RNAi and found that MINCR knockdown is associated with an impairment in cell cycle progression. Differential gene expression analysis after RNAi showed a significant enrichment of cell cycle genes among the genes down-regulated after MINCR knockdown. Interestingly, these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of the MYC transcriptional program. Accordingly, MINCR knockdown was associated with a reduction in MYC binding to the promoters of selected cell cycle genes. Finally, we show that down-regulation of Aurora kinases A and B and chromatin licensing and DNA replication factor 1 may explain the reduction in cellular proliferation observed on MINCR knockdown. We, therefore, suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis.

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project "Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing", we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.

The PCBP1 gene encoding poly(rC) binding protein I is recurrently mutated in Burkitt lymphoma.

The genetic hallmark of Burkitt lymphoma is the translocation t(8;14)(q24;q32), or one of its light chain variants, resulting in IG-MYC juxtaposition. However, these translocations alone are insufficient to drive lymphomagenesis, which requires additional genetic changes for malignant transformation. Recent studies of Burkitt lymphoma using next generation sequencing approaches have identified various recurrently mutated genes including ID3, TCF3, CCND3, and TP53. Here, by using similar approaches, we show that PCBP1 is a recurrently mutated gene in Burkitt lymphoma. By whole-genome sequencing, we identified somatic mutations in PCBP1 in 3/17 (18%) Burkitt lymphomas. We confirmed the recurrence of PCBP1 mutations by Sanger sequencing in an independent validation cohort, finding mutations in 3/28 (11%) Burkitt lymphomas and in 6/16 (38%) Burkitt lymphoma cell lines. PCBP1 is an intron-less gene encoding the 356 amino acid poly(rC) binding protein 1, which contains three K-Homology (KH) domains and two nuclear localization signals. The mutations predominantly (10/12, 83%) affect the KH III domain, either by complete domain loss or amino acid changes. Thus, these changes are predicted to alter the various functions of PCBP1, including nuclear trafficking and pre-mRNA splicing. Remarkably, all six primary Burkitt lymphomas with a PCBP1 mutation expressed MUM1/IRF4, which is otherwise detected in around 20-40% of Burkitt lymphomas. We conclude that PCBP1 mutations are recurrent in Burkitt lymphomas and might contribute, in cooperation with other mutations, to its pathogenesis.

Blood methylation pattern reflects epigenetic remodelling in adipose tissue after bariatric surgery.

BACKGROUND: Studies on DNA methylation following bariatric surgery have primarily focused on blood cells, while it is unclear to which extend it may reflect DNA methylation profiles in specific metabolically relevant organs such as adipose tissue. Here, we investigated whether adipose tissue depots specific methylation changes after bariatric surgery are mirrored in blood. METHODS: Using Illumina 850K EPIC technology, we analysed genome-wide DNA methylation in paired blood, subcutaneous and omental visceral AT (SAT/OVAT) samples from nine individuals (N = 6 female) with severe obesity pre- and post-surgery. FINDINGS: The numbers and effect sizes of differentially methylated regions (DMRs) post-bariatric surgery were more pronounced in AT (SAT: 12,865 DMRs from -11.5 to 10.8%; OVAT: 14,632 DMRs from -13.7 to 12.8%) than in blood (9267 DMRs from -8.8 to 7.7%). Cross-tissue DMRs implicated immune-related genes. Among them, 49 regions could be validated with similar methylation changes in blood from independent individuals. Fourteen DMRs correlated with differentially expressed genes in AT post bariatric surgery, including downregulation of PIK3AP1 in both SAT and OVAT. DNA methylation age acceleration was significantly higher in AT compared to blood, but remained unaffected after surgery. INTERPRETATION: Concurrent methylation pattern changes in blood and AT, particularly in immune-related genes, suggest blood DNA methylation mirrors AT's inflammatory state post-bariatric surgery. FUNDING: The funding sources are listed in the Acknowledgments section.

Folding RNA/DNA hybrid duplexes.

MOTIVATION: While there are numerous programs that can predict RNA or DNA secondary structures, a program that predicts RNA/DNA hetero-dimers is still missing. The lack of easy to use tools for predicting their structure may be in part responsible for the small number of reports of biologically relevant RNA/DNA hetero-dimers. RESULTS: We present here an extension to the widely used ViennaRNA Package (Lorenz et al., 2011) for the prediction of the structure of RNA/DNA hetero-dimers. AVAILABILITY: http://www.tbi.univie.ac.at/~ronny/RNA/vrna2.html CONTACT: ronny@tbi.univie.ac.at, berni@bioinf.uni-leipzig.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of two unannotated miRNAs in classic Hodgkin lymphoma cell lines.

MicroRNAs (miRNAs) are small non coding RNAs responsible for posttranscriptional regulation of gene expression. Even though almost 2000 precursors have been described so far, additional miRNAs are still being discovered in normal as well as malignant cells. Alike protein coding genes, miRNAs may acquire oncogenic properties in consequence of altered expression or presence of gain or loss of function mutations. In this study we mined datasets from miRNA expression profiling (miRNA-seq) of 7 classic Hodgkin Lymphoma (cHL) cell lines, 10 non-Hodgkin lymphoma (NHL) cell lines and 56 samples of germinal center derived B-cell lymphomas. Our aim was to discover potential novel cHL oncomiRs not reported in miRBase (release 22.1) and expressed in cHL cell lines but no other B-cell lymphomas. We identified six such miRNA candidates in cHL cell lines and verified the expression of two of them encoded at chr2:212678788-212678849 and chr5:168090507-168090561 (GRCh38). Interestingly, we showed that one of the validated miRNAs (located in an intron of the TENM2 gene) is expressed together with its host gene. TENM2 is characterized by hypomethylation and open chromatin around its TSS in cHL cell lines in contrast to NHL cell lines and germinal centre B-cells respectively. It indicates an epigenetic mechanism responsible for aberrant expression of both, the TENM2 gene and the novel miRNA in cHL cell lines. Despite the GO analysis performed with the input of the in silico predicted novel miRNA target genes did not reveal ontologies typically associated with cHL pathogenesis, it pointed to several interesting candidates involved in i.e. lymphopoiesis. These include the lymphoma related BCL11A gene, the IKZF2 gene involved in lymphocyte development or the transcription initiator GTF2H1.

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.