Successful RNAi applications depend on strategies allowing robust and persistent expression of minimal gene silencing triggers without perturbing endogenous gene expression. Here, we propose a novel avenue which is integration of a promoterless shmiRNA, i.e. a shRNA embedded in a micro-RNA (miRNA) scaffold, into an engineered genomic miRNA locus. For proof-of-concept, we used TALE or CRISPR/Cas9 nucleases to site-specifically integrate an anti-hepatitis C virus (HCV) shmiRNA into the liver-specific miR-122/hcr locus in hepatoma cells, with the aim to obtain cellular clones that are genetically protected against HCV infection. Using reporter assays, Northern blotting and qRT-PCR, we confirmed anti-HCV shmiRNA expression as well as miR-122 integrity and functionality in selected cellular progeny. Moreover, we employed a comprehensive battery of PCR, cDNA/miRNA profiling and whole genome sequencing analyses to validate targeted integration of a single shmiRNA molecule at the expected position, and to rule out deleterious effects on the genomes or transcriptomes of the engineered cells. Importantly, a subgenomic HCV replicon and a full-length reporter virus, but not a Dengue virus control, were significantly impaired in the modified cells. Our original combination of DNA engineering and RNAi expression technologies benefits numerous applications, from miRNA, genome and transgenesis research, to human gene therapy.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

Medulloblastoma (MB) is the most common malignant brain tumor in children. Patients whose tumors exhibit overexpression or amplification of the MYC oncogene (c-MYC) usually have an extremely poor prognosis, but there are no animal models of this subtype of the disease. Here, we show that cerebellar stem cells expressing Myc and mutant Trp53 (p53) generate aggressive tumors following orthotopic transplantation. These tumors consist of large, pleiomorphic cells and resemble human MYC-driven MB at a molecular level. Notably, antagonists of PI3K/mTOR signaling, but not Hedgehog signaling, inhibit growth of tumor cells. These findings suggest that cerebellar stem cells can give rise to MYC-driven MB and identify a novel model that can be used to test therapies for this devastating disease.

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.

Genomic characterization has begun to redefine diagnostic classifications of cancers. However, it remains a challenge to infer disease phenotypes from genomic alterations alone. To help realize the promise of genomics, we have performed a quantitative proteomics investigation using Stable Isotope Labeling by Amino Acids in Cell Culture (SILAC) and 41 tissue samples spanning the 4 genomically based subgroups of medulloblastoma and control cerebellum. We have identified and quantitated thousands of proteins across these groups and find that we are able to recapitulate the genomic subgroups based upon subgroup restricted and differentially abundant proteins while also identifying subgroup specific protein isoforms. Integrating our proteomic measurements with genomic data, we calculate a poor correlation between mRNA and protein abundance. Using EPIC 850 k methylation array data on the same tissues, we also investigate the influence of copy number alterations and DNA methylation on the proteome in an attempt to characterize the impact of these genetic features on the proteome. Reciprocally, we are able to use the proteome to identify which genomic alterations result in altered protein abundance and thus are most likely to impact biology. Finally, we are able to assemble protein-based pathways yielding potential avenues for clinical intervention. From these, we validate the EIF4F cap-dependent translation pathway as a novel druggable pathway in medulloblastoma. Thus, quantitative proteomics complements genomic platforms to yield a more complete understanding of functional tumor biology and identify novel therapeutic targets for medulloblastoma.

Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology.

Musashi1 (Msi1) is a highly conserved RNA-binding protein that is required during the development of the nervous system. Msi1 has been characterized as a stem cell marker, controlling the balance between self-renewal and differentiation, and has also been implicated in tumorigenesis, being highly expressed in multiple tumor types. We analyzed Msi1 expression in a large cohort of medulloblastoma samples and found that Msi1 is highly expressed in tumor tissue compared with normal cerebellum. Notably, high Msi1 expression levels proved to be a sign of poor prognosis. Msi1 expression was determined to be particularly high in molecular subgroups 3 and 4 of medulloblastoma. We determined that Msi1 is required for tumorigenesis because inhibition of Msi1 expression by small-interfering RNAs reduced the growth of Daoy medulloblastoma cells in xenografts. To characterize the participation of Msi1 in medulloblastoma, we conducted different high-throughput analyses. Ribonucleoprotein immunoprecipitation followed by microarray analysis (RIP-chip) was used to identify mRNA species preferentially associated with Msi1 protein in Daoy cells. We also used cluster analysis to identify genes with similar or opposite expression patterns to Msi1 in our medulloblastoma cohort. A network study identified RAC1, CTGF, SDCBP, SRC, PRL, and SHC1 as major nodes of an Msi1-associated network. Our results suggest that Msi1 functions as a regulator of multiple processes in medulloblastoma formation and could become an important therapeutic target.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

Mutation is a fundamental process in tumorigenesis. However, the degree to which the rate of somatic mutation varies across the human genome and the mechanistic basis underlying this variation remain to be fully elucidated. Here, we performed a cross-cancer comparison of 402 whole genomes comprising a diverse set of childhood and adult tumors, including both solid and hematopoietic malignancies. Surprisingly, we found that the inactive X chromosome of many female cancer genomes accumulates on average twice and up to four times as many somatic mutations per megabase, as compared to the individual autosomes. Whole-genome sequencing of clonally expanded hematopoietic stem/progenitor cells (HSPCs) from healthy individuals and a premalignant myelodysplastic syndrome (MDS) sample revealed no X chromosome hypermutation. Our data suggest that hypermutation of the inactive X chromosome is an early and frequent feature of tumorigenesis resulting from DNA replication stress in aberrantly proliferating cells.

Massive parallel sequencing (MPS) technologies have paved the way into new areas of research including individualized medicine. However, sequencing of trace amounts of DNA or RNA still remains a major challenge, especially for degraded nucleic acids like circulating DNA. This together with high cost and time requirements impedes many important applications of MPS in medicine and fundamental science. We have established a fast, cheap and highly efficient protocol called 'Capture and Amplification by Tailing and Switching' (CATS) to directly generate ready-to-sequence libraries for MPS from nanogram and picogram quantities of both DNA and RNA. Furthermore, those DNA libraries are strand-specific, can be prepared within 2-3 h and do not require preliminary sample amplification steps. To exemplify the capacity of the technique, we have generated and sequenced DNA libraries from hundred-picogram amounts of circulating nucleic acids isolated from human blood plasma, one nanogram of mRNA-enriched total RNA from cultured cells and few nanograms of bisulfite-converted DNA. The approach for DNA library preparation from minimal and fragmented input described here will find broad application in diverse research areas such as translational medicine including therapy monitoring, prediction, prognosis and early detection of various human disorders and will permit high-throughput DNA sequencing from previously inaccessible material such as minute forensic and archeological samples.

Recurrent medulloblastoma is a therapeutic challenge because it is almost always fatal. Studies have confirmed that medulloblastoma consists of at least four distinct subgroups. We sought to delineate subgroup-specific differences in medulloblastoma recurrence patterns.

Two recurrent mutations, K27M and G34R/V, within histone variant H3.3 were recently identified in ∼50% of pHGGs. Both mutations define clinically and biologically distinct subgroups of pHGGs. Here, we provide further insight about the dominant-negative effect of K27M mutant H3.3, leading to a global reduction of the repressive histone mark H3K27me3. We demonstrate that this is caused by aberrant recruitment of the PRC2 complex to K27M mutant H3.3 and enzymatic inhibition of the H3K27me3-establishing methyltransferase EZH2. By performing chromatin immunoprecipitation followed by next-generation sequencing and whole-genome bisulfite sequencing in primary pHGGs, we show that reduced H3K27me3 levels and DNA hypomethylation act in concert to activate gene expression in K27M mutant pHGGs.

The most aggressive of four medulloblastoma (MB) subgroups are cMyc-driven group 3 (G3) tumors, some of which overexpress EZH2, the histone H3K27 mono-, di-, and trimethylase of polycomb-repressive complex 2. Ezh2 has a context-dependent role in different cancers as an oncogene or tumor suppressor and retards tumor progression in a mouse model of G3 MB. Engineered deletions of Ezh2 in G3 MBs by gene editing nucleases accelerated tumorigenesis, whereas Ezh2 re-expression reversed attendant histone modifications and slowed tumor progression. Candidate oncogenic drivers suppressed by Ezh2 included Gfi1, a proto-oncogene frequently activated in human G3 MBs. Gfi1 disruption antagonized the tumor-promoting effects of Ezh2 loss; conversely, Gfi1 overexpression collaborated with Myc to bypass effects of Trp53 inactivation in driving MB progression in primary cerebellar neuronal progenitors. Although negative regulation of Gfi1 by Ezh2 may restrain MB development, Gfi1 activation can bypass these effects.

A hierarchically organized cell compartment drives colorectal cancer (CRC) progression. Genetic barcoding allows monitoring of the clonal output of tumorigenic cells without prospective isolation. In this study, we asked whether tumor clone-initiating cells (TcICs) were genetically heterogeneous and whether differences in self-renewal and activation reflected differential kinetics among individual subclones or functional hierarchies within subclones. Monitoring genomic subclone kinetics in three patient tumors and corresponding serial xenografts and spheroids by high-coverage whole-genome sequencing, clustering of genetic aberrations, subclone combinatorics, and mutational signature analysis revealed at least two to four genetic subclones per sample. Long-term growth in serial xenografts and spheroids was driven by multiple genomic subclones with profoundly differing growth dynamics and hence different quantitative contributions over time. Strikingly, genetic barcoding demonstrated stable functional heterogeneity of CRC TcICs during serial xenografting despite near-complete changes in genomic subclone contribution. This demonstrates that functional heterogeneity is, at least frequently, present within genomic subclones and independent of mutational subclone differences.

Classification and variable selection play an important role in knowledge discovery in high-dimensional data. Although Support Vector Machine (SVM) algorithms are among the most powerful classification and prediction methods with a wide range of scientific applications, the SVM does not include automatic feature selection and therefore a number of feature selection procedures have been developed. Regularisation approaches extend SVM to a feature selection method in a flexible way using penalty functions like LASSO, SCAD and Elastic Net.We propose a novel penalty function for SVM classification tasks, Elastic SCAD, a combination of SCAD and ridge penalties which overcomes the limitations of each penalty alone.Since SVM models are extremely sensitive to the choice of tuning parameters, we adopted an interval search algorithm, which in comparison to a fixed grid search finds rapidly and more precisely a global optimal solution.

The emergence of high-throughput, next-generation sequencing technologies has dramatically altered the way we assess genomes in population genetics and in cancer genomics. Currently, there are four commonly used whole-genome sequencing platforms on the market: Illumina's HiSeq2000, Life Technologies' SOLiD 4 and its completely redesigned 5500xl SOLiD, and Complete Genomics' technology. A number of earlier studies have compared a subset of those sequencing platforms or compared those platforms with Sanger sequencing, which is prohibitively expensive for whole genome studies. Here we present a detailed comparison of the performance of all currently available whole genome sequencing platforms, especially regarding their ability to call SNVs and to evenly cover the genome and specific genomic regions. Unlike earlier studies, we base our comparison on four different samples, allowing us to assess the between-sample variation of the platforms. We find a pronounced GC bias in GC-rich regions for Life Technologies' platforms, with Complete Genomics performing best here, while we see the least bias in GC-poor regions for HiSeq2000 and 5500xl. HiSeq2000 gives the most uniform coverage and displays the least sample-to-sample variation. In contrast, Complete Genomics exhibits by far the smallest fraction of bases not covered, while the SOLiD platforms reveal remarkable shortcomings, especially in covering CpG islands. When comparing the performance of the four platforms for calling SNPs, HiSeq2000 and Complete Genomics achieve the highest sensitivity, while the SOLiD platforms show the lowest false positive rate. Finally, we find that integrating sequencing data from different platforms offers the potential to combine the strengths of different technologies. In summary, our results detail the strengths and weaknesses of all four whole-genome sequencing platforms. It indicates application areas that call for a specific sequencing platform and disallow other platforms. This helps to identify the proper sequencing platform for whole genome studies with different application scopes.

Gene knock down by RNAi is a highly effective approach to silence gene expression in experimental as well as therapeutic settings. However, this widely used methodology entails serious pitfalls, especially concerning specificity of the RNAi molecules.

Chronic lymphocytic leukemia is a malignancy of mature B cells that strongly depend on microenvironmental factors, and their deprivation has been identified as a promising treatment approach for this incurable disease. Cytokine array screening of 247 chronic lymphocytic leukemia serum samples revealed elevated levels of tumor necrosis factor (TNF) receptor-1 which were associated with poor clinical outcome. We detected a microenvironment-induced expression of TNF receptor-1 in chronic lymphocytic leukemia cells in vitro, and an aberrantly high expression of this receptor in the proliferation centers of patients' lymph nodes. Stimulation of TNF receptor-1 with TNF-α enhanced nuclear factor κ-light-chain-enhancer of activated B cells (NFκB) activity and viability of chronic lymphocytic leukemia cells, which was inhibited by wogonin. The therapeutic effects of wogonin were analyzed in mice after adoptive transfer of Eμ-T-cell leukemia 1 (TCL1) leukemic cells. Wogonin treatment prevented leukemia development when given early after transplantation. The treatment of full-blown leukemia resulted in the loss of the TNF receptor-1 on chronic lymphocytic leukemia cells and their mobilization to blood. Targeting TNF receptor-1 signaling is therefore proposed for the treatment of chronic lymphocytic leukemia.

Over the past decade, wingless-activated (WNT) medulloblastoma has been identified as a candidate for therapy de-escalation based on excellent survival; however, a paucity of relapses has precluded additional analyses of markers of relapse. To address this gap in knowledge, an international cohort of 93 molecularly confirmed WNT MB was assembled, where 5-year progression-free survival is 0.84 (95%, 0.763-0.925) with 15 relapsed individuals identified. Maintenance chemotherapy is identified as a strong predictor of relapse, with individuals receiving high doses of cyclophosphamide or ifosphamide having only one very late molecularly confirmed relapse (p = 0.032). The anatomical location of recurrence is metastatic in 12 of 15 relapses, with 8 of 12 metastatic relapses in the lateral ventricles. Maintenance chemotherapy, specifically cumulative cyclophosphamide doses, is a significant predictor of relapse across WNT MB. Future efforts to de-escalate therapy need to carefully consider not only the radiation dose but also the chemotherapy regimen and the propensity for metastatic relapses.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .