Gene expression profiling is being widely applied in cancer research to identify biomarkers for clinical endpoint prediction. Since RNA-seq provides a powerful tool for transcriptome-based applications beyond the limitations of microarrays, we sought to systematically evaluate the performance of RNA-seq-based and microarray-based classifiers in this MAQC-III/SEQC study for clinical endpoint prediction using neuroblastoma as a model.

A remarkable observation emerging from recent cancer genome analyses is the identification of chromothripsis as a one-off genomic catastrophe, resulting in massive somatic DNA structural rearrangements (SRs). Largely due to lack of suitable model systems, the mechanistic basis of chromothripsis has remained elusive. We developed an integrative method termed "complex alterations after selection and transformation (CAST)," enabling efficient in vitro generation of complex DNA rearrangements including chromothripsis, using cell perturbations coupled with a strong selection barrier followed by massively parallel sequencing. We employed this methodology to characterize catastrophic SR formation processes, their temporal sequence, and their impact on gene expression and cell division. Our in vitro system uncovered a propensity of chromothripsis to occur in cells with damaged telomeres, and in particular in hyperploid cells. Analysis of primary medulloblastoma cancer genomes verified the link between hyperploidy and chromothripsis in vivo. CAST provides the foundation for mechanistic dissection of complex DNA rearrangement processes.

Harmonizing quality control (QC) of large-scale second and third-generation sequencing datasets is key for enabling downstream computational and biological analyses. We present Alfred, an efficient and versatile command-line application that computes multi-sample QC metrics in a read-group aware manner, across a wide variety of sequencing assays and technologies. In addition to standard QC metrics such as GC bias, base composition, insert size and sequencing coverage distributions it supports haplotype-aware and allele-specific feature counting and feature annotation. The versatility of Alfred allows for easy pipeline integration in high-throughput settings, including DNA sequencing facilities and large-scale research initiatives, enabling continuous monitoring of sequence data quality and characteristics across samples. Alfred supports haplo-tagging of BAM/CRAM files to conduct haplotype-resolved analyses in conjunction with a variety of next-generation sequencing based assays. Alfred's companion web application enables interactive exploration of results and comparison to public datasets.

The NUP214-ABL1 fusion is a constitutively activated tyrosine kinase that is significantly associated with overexpression of the TLX1 and TLX3 transcription factors in T cell acute lymphoblastic leukemia (T-ALL). Here we show that NUP214-ABL1 cooperates with TLX1 in driving T-ALL development using a transgenic mouse model and human T-ALL cells. Using integrated ChIP-sequencing, ATAC-sequencing, and RNA-sequencing data, we demonstrate that TLX1 and STAT5, the downstream effector of NUP214-ABL1, co-bind poised enhancer regions, and cooperatively activate the expression of key proto-oncogenes such as MYC and BCL2. Inhibition of STAT5, downregulation of TLX1 or MYC, or interference with enhancer function through BET-inhibitor treatment leads to reduction of target gene expression and induction of leukemia cell death.

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Nonsense-mediated mRNA decay (NMD) is a cellular surveillance pathway that recognizes and degrades mRNAs with premature termination codons (PTCs). The mechanisms underlying translation termination are key to the understanding of RNA surveillance mechanisms such as NMD and crucial for the development of therapeutic strategies for NMD-related diseases. Here, we have used a fully reconstituted in vitro translation system to probe the NMD proteins for interaction with the termination apparatus. We discovered that UPF3B (i) interacts with the release factors, (ii) delays translation termination and (iii) dissociates post-termination ribosomal complexes that are devoid of the nascent peptide. Furthermore, we identified UPF1 and ribosomes as new interaction partners of UPF3B. These previously unknown functions of UPF3B during the early and late phases of translation termination suggest that UPF3B is involved in the crosstalk between the NMD machinery and the PTC-bound ribosome, a central mechanistic step of RNA surveillance.

MYCN is a master regulator controlling many processes necessary for tumor cell survival. Here, we unravel a microRNA network that causes tumor suppressive effects in MYCN-amplified neuroblastoma cells. In profiling studies, histone deacetylase (HDAC) inhibitor treatment most strongly induced miR-183. Enforced miR-183 expression triggered apoptosis, and inhibited anchorage-independent colony formation in vitro and xenograft growth in mice. Furthermore, the mechanism of miR-183 induction was found to contribute to the cell death phenotype induced by HDAC inhibitors. Experiments to identify the HDAC(s) involved in miR-183 transcriptional regulation showed that HDAC2 depletion induced miR-183. HDAC2 overexpression reduced miR-183 levels and counteracted the induction caused by HDAC2 depletion or HDAC inhibitor treatment. MYCN was found to recruit HDAC2 in the same complexes to the miR-183 promoter, and HDAC2 depletion enhanced promoter-associated histone H4 pan-acetylation, suggesting epigenetic changes preceded transcriptional activation. These data reveal miR-183 tumor suppressive properties in neuroblastoma that are jointly repressed by MYCN and HDAC2, and suggest a novel way to bypass MYCN function.

Accurate estimates of mutation rates provide critical information to analyze genome evolution and organism fitness. We used whole-genome DNA sequencing, pulse-field gel electrophoresis, and comparative genome hybridization to determine mutation rates in diploid vegetative and meiotic mutation accumulation lines of Saccharomyces cerevisiae. The vegetative lines underwent only mitotic divisions while the meiotic lines underwent a meiotic cycle every ∼20 vegetative divisions. Similar base substitution rates were estimated for both lines. Given our experimental design, these measures indicated that the meiotic mutation rate is within the range of being equal to zero to being 55-fold higher than the vegetative rate. Mutations detected in vegetative lines were all heterozygous while those in meiotic lines were homozygous. A quantitative analysis of intra-tetrad mating events in the meiotic lines showed that inter-spore mating is primarily responsible for rapidly fixing mutations to homozygosity as well as for removing mutations. We did not observe 1-2 nt insertion/deletion (in-del) mutations in any of the sequenced lines and only one structural variant in a non-telomeric location was found. However, a large number of structural variations in subtelomeric sequences were seen in both vegetative and meiotic lines that did not affect viability. Our results indicate that the diploid yeast nuclear genome is remarkably stable during the vegetative and meiotic cell cycles and support the hypothesis that peripheral regions of chromosomes are more dynamic than gene-rich central sections where structural rearrangements could be deleterious. This work also provides an improved estimate for the mutational load carried by diploid organisms.

Nonsense-mediated mRNA decay (NMD) downmodulates mRNAs that have in-frame premature termination codons and prevents translation of potentially harmful truncated proteins from aberrant mRNAs. Two new approaches have identified physiological NMD substrates, and suggest that NMD functions as a multipurpose tool in the modulation of gene expression.

Double-strand breaks (DSBs) are the most toxic type of DNA lesions. Cells repair these lesions using either end protection- or end resection-coupled mechanisms. To study DSB repair choice, we present the Color Assay Tracing-Repair (CAT-R) to simultaneously quantify DSB repair via end protection and end resection pathways. CAT-R introduces DSBs using CRISPR/Cas9 in a tandem fluorescent reporter, whose repair distinguishes small insertions/deletions from large deletions. We demonstrate CAT-R applications in chemical and genetic screens. First, we evaluate 21 compounds currently in clinical trials which target the DNA damage response. Second, we examine how 417 factors involved in DNA damage response influence the choice between end protection and end resection. Finally, we show that impairing nucleotide excision repair favors error-free repair, providing an alternative way for improving CRISPR/Cas9-based knock-ins. CAT-R is a high-throughput, versatile assay to assess DSB repair choice, which facilitates comprehensive studies of DNA repair and drug efficiency testing.

The occurrence and formation of genomic structural variants (SVs) is known to be influenced by the 3D chromatin architecture, but the extent and magnitude have been challenging to study. Here, we apply Hi-C to study chromatin organization before and after induction of chromothripsis in human cells. We use Hi-C to manually assemble the derivative chromosomes following the occurrence of massive complex rearrangements, which allows us to study the sources of SV formation and their consequences on gene regulation. We observe an action-reaction interplay whereby the 3D chromatin architecture directly impacts the location and formation of SVs. In turn, the SVs reshape the chromatin organization to alter the local topologies, replication timing, and gene regulation in cis We show that SVs have a strong tendency to occur between similar chromatin compartments and replication timing regions. Moreover, we find that SVs frequently occur at 3D loop anchors, that SVs can cause a switch in chromatin compartments and replication timing, and that this is a major source of SV-mediated effects on nearby gene expression changes. Finally, we provide evidence for a general mechanistic bias of the 3D chromatin on SV occurrence using data from more than 2700 patient-derived cancer genomes.

Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

Deciphering translation is of paramount importance for the understanding of many diseases, and antibiotics played a pivotal role in this endeavour. Blasticidin S (BlaS) targets translation by binding to the peptidyl transferase center of the large ribosomal subunit. Using biochemical, structural and cellular approaches, we show here that BlaS inhibits both translation elongation and termination in Mammalia. Bound to mammalian terminating ribosomes, BlaS distorts the 3'CCA tail of the P-site tRNA to a larger extent than previously reported for bacterial ribosomes, thus delaying both, peptide bond formation and peptidyl-tRNA hydrolysis. While BlaS does not inhibit stop codon recognition by the eukaryotic release factor 1 (eRF1), it interferes with eRF1's accommodation into the peptidyl transferase center and subsequent peptide release. In human cells, BlaS inhibits nonsense-mediated mRNA decay and, at subinhibitory concentrations, modulates translation dynamics at premature termination codons leading to enhanced protein production.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Complexity of lung microenvironment and changes in cellular composition during disease make it exceptionally hard to understand molecular mechanisms driving development of chronic lung diseases. Although recent advances in cell type-resolved approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a versatile workflow that allows storage of lung tissue with high viability, permits thorough sample quality check before cell isolation, and befits sequencing-based profiling. We demonstrate that cryopreservation enables isolation of multiple cell types from both healthy and diseased lungs. Basal cells from cryopreserved airways retain their differentiation ability, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for omics profiling of lung cells. Moreover, we obtained high-quality single-cell RNA-sequencing data of cells from cryopreserved human lungs, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our workflow is well suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to viable human tissue.

Chromosome instability (CIN) is associated with poor survival and therapeutic outcome in a number of malignancies. Despite this correlation, CIN can also lead to growth disadvantages. Here, we show that simultaneous overexpression of the mitotic checkpoint protein Mad2 with Kras(G12D) or Her2 in mammary glands of adult mice results in mitotic checkpoint overactivation and a delay in tumor onset. Time-lapse imaging of organotypic cultures and pathologic analysis prior to tumor establishment reveals error-prone mitosis, mitotic arrest, and cell death. Nonetheless, Mad2 expression persists and increases karyotype complexity in Kras tumors. Faced with the selective pressure of oncogene withdrawal, Mad2-positive tumors have a higher frequency of developing persistent subclones that avoid remission and continue to grow.

Mammalian nonsense-mediated mRNA decay (NMD) is a eukaryotic surveillance mechanism that degrades mRNAs containing premature translation termination codons. Phosphorylation of the essential NMD effector UPF1 by the phosphoinositide-3-kinase-like kinase (PIKK) SMG-1 is a key step in NMD and occurs when SMG-1, its two regulatory factors SMG-8 and SMG-9, and UPF1 form a complex at a terminating ribosome. Electron cryo-microscopy of the SMG-1-8-9-UPF1 complex shows the head and arm architecture characteristic of PIKKs and reveals different states of UPF1 docking. UPF1 is recruited to the SMG-1 kinase domain and C-terminal insertion domain, inducing an opening of the head domain that provides access to the active site. SMG-8 and SMG-9 interact with the SMG-1 C-insertion and promote high-affinity UPF1 binding to SMG-1-8-9, as well as decelerated SMG-1 kinase activity and enhanced stringency of phosphorylation site selection. The presence of UPF2 destabilizes the SMG-1-8-9-UPF1 complex leading to substrate release. Our results suggest an intricate molecular network of SMG-8, SMG-9 and the SMG-1 C-insertion domain that governs UPF1 substrate recruitment and phosphorylation by SMG-1 kinase, an event that is central to trigger mRNA decay.

The exon junction complex (EJC) participates in the regulation of many post-transcriptional steps of gene expression. EJCs are deposited on messenger RNAs (mRNAs) during splicing and their core consists of eIF4A3, MLN51, Y14, and MAGOH. Here, we show that two genes encoding MAGOH paralogs (referred to as MAGOH and MAGOHB) are expressed in mammals. In macrophages, the expression of MAGOHB, but not MAGOH mRNA, increases rapidly after LPS stimulation. Both MAGOH proteins interact with other EJC components, incorporate into mRNA-bound EJCs, and activate nonsense-mediated decay. Furthermore, the simultaneous depletion of MAGOH and MAGOHB, but not individual depletions, impair nonsense-mediated decay in human cells. Hence, our results establish that the core composition of mammalian EJCs is more complex than previously recognized.

Synapsins are a family of neuronal phosphoproteins associated with the cytosolic surface of synaptic vesicles. Experimental evidence suggests a role for synapsins in synaptic vesicle clustering and recycling at the presynaptic terminal, as well as in neuronal development and synaptogenesis. Synapsin knock-out (Syn1(-/-) ) mice display an epileptic phenotype and mutations in the SYN1 gene have been identified in individuals affected by epilepsy and/or autism spectrum disorder. We investigated the impact of the c.1067G>A nonsense transition, the first mutation described in a family affected by X-linked syndromic epilepsy, on the expression and functional properties of the synapsin I protein. We found that the presence of a premature termination codon in the human SYN1 transcript renders it susceptible to nonsense-mediated mRNA decay (NMD). Given that the NMD efficiency is highly variable among individuals and cell types, we investigated also the effects of expression of the mutant protein and found that it is expressed at lower levels compared to wild-type synapsin I, forms perinuclear aggregates and is unable to reach presynaptic terminals in mature hippocampal neurons grown in culture. Taken together, these data indicate that in patients carrying the W356× mutation the function of synapsin I is markedly impaired, due to both the strongly decreased translation and the altered function of the NMD-escaped protein, and support the value of Syn1(-/-) mice as an experimental model mimicking the human pathology.

As a consequence of the accumulation of insertion events over evolutionary time, mobile elements now comprise nearly half of the human genome. The Alu, L1, and SVA mobile element families are still duplicating, generating variation between individual genomes. Mobile element insertions (MEI) have been identified as causes for genetic diseases, including hemophilia, neurofibromatosis, and various cancers. Here we present a comprehensive map of 7,380 MEI polymorphisms from the 1000 Genomes Project whole-genome sequencing data of 185 samples in three major populations detected with two detection methods. This catalog enables us to systematically study mutation rates, population segregation, genomic distribution, and functional properties of MEI polymorphisms and to compare MEI to SNP variation from the same individuals. Population allele frequencies of MEI and SNPs are described, broadly, by the same neutral ancestral processes despite vastly different mutation mechanisms and rates, except in coding regions where MEI are virtually absent, presumably due to strong negative selection. A direct comparison of MEI and SNP diversity levels suggests a differential mobile element insertion rate among populations.