Genome sequencing of pediatric medulloblastoma links catastrophic DNA rearrangements with TP53 mutations.

Genomic rearrangements are thought to occur progressively during tumor development. Recent findings, however, suggest an alternative mechanism, involving massive chromosome rearrangements in a one-step catastrophic event termed chromothripsis. We report the whole-genome sequencing-based analysis of a Sonic-Hedgehog medulloblastoma (SHH-MB) brain tumor from a patient with a germline TP53 mutation (Li-Fraumeni syndrome), uncovering massive, complex chromosome rearrangements. Integrating TP53 status with microarray and deep sequencing-based DNA rearrangement data in additional patients reveals a striking association between TP53 mutation and chromothripsis in SHH-MBs. Analysis of additional tumor entities substantiates a link between TP53 mutation and chromothripsis, and indicates a context-specific role for p53 in catastrophic DNA rearrangements. Among these, we observed a strong association between somatic TP53 mutations and chromothripsis in acute myeloid leukemia. These findings connect p53 status and chromothripsis in specific tumor types, providing a genetic basis for understanding particularly aggressive subtypes of cancer.

Somatic structural variant formation is guided by and influences genome architecture.

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.

An integrated map of structural variation in 2,504 human genomes.

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma.

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

Assembly and diploid architecture of an individual human genome via single-molecule technologies.

We present the first comprehensive analysis of a diploid human genome that combines single-molecule sequencing with single-molecule genome maps. Our hybrid assembly markedly improves upon the contiguity observed from traditional shotgun sequencing approaches, with scaffold N50 values approaching 30 Mb, and we identified complex structural variants (SVs) missed by other high-throughput approaches. Furthermore, by combining Illumina short-read data with long reads, we phased both single-nucleotide variants and SVs, generating haplotypes with over 99% consistency with previous trio-based studies. Our work shows that it is now possible to integrate single-molecule and high-throughput sequence data to generate de novo assembled genomes that approach reference quality.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes.

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-ß signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

Mapping copy number variation by population-scale genome sequencing.

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

Identification of ZCCHC8 as fusion partner of ROS1 in a case of congenital glioblastoma multiforme with a t(6;12)(q21;q24.3).

Congenital gliobastoma multiforme (GBM) is rare and little is known about the molecular defects underlying the initiation and progression of this tumor type. We present a case of congenital GBM analyzed by conventional cytogenetics, fluorescence in situ hybridization, array comparative genomic hybridization and next generation sequencing. On cytogenetic analysis we detected a reciprocal translocation t(6;12)(q21;q24.3). By sequencing, the translocation was shown to form a fusion between the 5' region of ZCCHC8 and the 3' region of ROS1. RT-PCR analyses confirmed the existence of an in-frame fusion transcript with ZCCHC8 exons 1-3 joined to ROS1 exons 36-43. In addition to the ZCCHC8-ROS1 fusion, we detected a deletion in the short arm of chromosome 9, including homozygous loss of the CDKN2A/2B locus in 9p21.3 and heterozygous deletion of the HAUS6 gene in 9p22.1. The latter encodes a protein involved in faithful chromosome segregation by regulating microtubule nucleation and its deletion might be associated with the marked subclonal changes of ploidy observed in the tumor. This report adds the ZCCHC8-ROS1 fusion as oncogenic driver in GBM and supports the role of ROS1 activation in the pathogenesis of a subset of GBM. © 2016 Wiley Periodicals, Inc.

MCM3AP and POMP Mutations Cause a DNA-Repair and DNA-Damage-Signaling Defect in an Immunodeficient Child.

Immunodeficiency patients with DNA repair defects exhibit radiosensitivity and proneness to leukemia/lymphoma formation. Though progress has been made in identifying the underlying mutations, in most patients the genetic basis is unknown. Two de novo mutated candidate genes, MCM3AP encoding germinal center-associated nuclear protein (GANP) and POMP encoding proteasome maturation protein (POMP), were identified by whole-exome sequencing (WES) and confirmed by Sanger sequencing in a child with complex phenotype displaying immunodeficiency, genomic instability, skin changes, and myelodysplasia. GANP was previously described to promote B-cell maturation by nuclear targeting of activation-induced cytidine deaminase (AID) and to control AID-dependent hyperrecombination. POMP is required for 20S proteasome assembly and, thus, for efficient NF-κB signaling. Patient-derived cells were characterized by impaired homologous recombination, moderate radio- and cross-linker sensitivity associated with accumulation of damage, impaired DNA damage-induced NF-κB signaling, and reduced nuclear AID levels. Complementation by wild-type (WT)-GANP normalized DNA repair and WT-POMP rescued defective NF-κB signaling. In conclusion, we identified for the first time mutations in MCM3AP and POMP in an immunodeficiency patient. These mutations lead to cooperative effects on DNA recombination and damage signaling. Digenic/polygenic mutations may constitute a novel genetic basis in immunodeficiency patients with DNA repair defects.

Next-generation sequencing-based detection of germline L1-mediated transductions.

BACKGROUND: While active LINE-1 (L1) elements possess the ability to mobilize flanking sequences to different genomic loci through a process termed transduction influencing genomic content and structure, an approach for detecting polymorphic germline non-reference transductions in massively-parallel sequencing data has been lacking. RESULTS: Here we present the computational approach TIGER (Transduction Inference in GERmline genomes), enabling the discovery of non-reference L1-mediated transductions by combining L1 discovery with detection of unique insertion sequences and detailed characterization of insertion sites. We employed TIGER to characterize polymorphic transductions in fifteen genomes from non-human primate species (chimpanzee, orangutan and rhesus macaque), as well as in a human genome. We achieved high accuracy as confirmed by PCR and two single molecule DNA sequencing techniques, and uncovered differences in relative rates of transduction between primate species. CONCLUSIONS: By enabling detection of polymorphic transductions, TIGER makes this form of relevant structural variation amenable for population and personal genome analysis.

Impact of genomic structural variation in Drosophila melanogaster based on population-scale sequencing.

Genomic structural variation (SV) is a major determinant for phenotypic variation. Although it has been extensively studied in humans, the nucleotide resolution structure of SVs within the widely used model organism Drosophila remains unknown. We report a highly accurate, densely validated map of unbalanced SVs comprising 8962 deletions and 916 tandem duplications in 39 lines derived from short-read DNA sequencing in a natural population (the "Drosophila melanogaster Genetic Reference Panel," DGRP). Most SVs (>90%) were inferred at nucleotide resolution, and a large fraction was genotyped across all samples. Comprehensive analyses of SV formation mechanisms using the short-read data revealed an abundance of SVs formed by mobile element and nonhomologous end-joining-mediated rearrangements, and clustering of variants into SV hotspots. We further observed a strong depletion of SVs overlapping genes, which, along with population genetics analyses, suggests that these SVs are often deleterious. We inferred several gene fusion events also highlighting the potential role of SVs in the generation of novel protein products. Expression quantitative trait locus (eQTL) mapping revealed the functional impact of our high-resolution SV map, with quantifiable effects at >100 genic loci. Our map represents a resource for population-level studies of SVs in an important model organism.

A cell-based model system links chromothripsis with hyperploidy.

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.

Primate genome architecture influences structural variation mechanisms and functional consequences.

Although nucleotide resolution maps of genomic structural variants (SVs) have provided insights into the origin and impact of phenotypic diversity in humans, comparable maps in nonhuman primates have thus far been lacking. Using massively parallel DNA sequencing, we constructed fine-resolution genomic structural variation maps in five chimpanzees, five orang-utans, and five rhesus macaques. The SV maps, which are comprised of thousands of deletions, duplications, and mobile element insertions, revealed a high activity of retrotransposition in macaques compared with great apes. By comparison, nonallelic homologous recombination is specifically active in the great apes, which is correlated with architectural differences between the genomes of great apes and macaque. Transcriptome analyses across nonhuman primates and humans revealed effects of species-specific whole-gene duplication on gene expression. We identified 13 gene duplications coinciding with the species-specific gain of tissue-specific gene expression in keeping with a role of gene duplication in the promotion of diversification and the acquisition of unique functions. Differences in the present day activity of SV formation mechanisms that our study revealed may contribute to ongoing diversification and adaptation of great ape and Old World monkey lineages.

Pediatric T-cell lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation.

Relapsed precursor T-cell acute lymphoblastic leukemia is characterized by resistance against chemotherapy and is frequently fatal. We aimed at understanding the molecular mechanisms resulting in relapse of T-cell acute lymphoblastic leukemia and analyzed 13 patients at first diagnosis, remission and relapse by whole exome sequencing, targeted ultra-deep sequencing, multiplex ligation dependent probe amplification and DNA methylation array. Compared to primary T-cell acute lymphoblastic leukemia, in relapse the number of single nucleotide variants and small insertions and deletions approximately doubled from 11.5 to 26. Targeted ultra-deep sequencing sensitively detected subclones that were selected for in relapse. The mutational pattern defined two types of relapses. While both are characterized by selection of subclones and acquisition of novel mutations, 'type 1' relapse derives from the primary leukemia whereas 'type 2' relapse originates from a common pre-leukemic ancestor. Relapse-specific changes included activation of the nucleotidase NT5C2 resulting in resistance to chemotherapy and mutations of epigenetic modulators, exemplified by SUZ12, WHSC1 and SMARCA4. While mutations present in primary leukemia and in relapse were enriched for known drivers of leukemia, relapse-specific changes revealed an association with general cancer-promoting mechanisms. This study thus identifies mechanisms that drive progression of pediatric T-cell acute lymphoblastic leukemia to relapse and may explain the characteristic treatment resistance of this condition.

A common microdeletion affecting a hippocampus- and amygdala-specific isoform of tryptophan hydroxylase 2 is not associated with affective disorders.

OBJECTIVES: Copy number variants (CNVs) have been shown to affect susceptibility for neuropsychiatric disorders. To date, studies implicating the serotonergic system in complex conditions have just focused on single nucleotide polymorphisms (SNPs). We therefore sought to identify novel common genetic copy number polymorphisms affecting genes of the serotonergic system, and to assess their putative role in bipolar affective disorder (BPAD) and major depressive disorder (MDD). METHODS: A selection of 41 genes of the serotonergic system encoding receptors, the serotonin transporter, metabolic enzymes and chaperones were investigated using a paired-end mapping (PEM) approach on next-generation sequencing data from the pilot project of the 1000 Genomes Project. For association testing, 593 patients with MDD, 1,145 patients with BPAD, and 1,738 healthy controls were included in the study. RESULTS: PEM led to the identification of a microdeletion in the gene encoding tryptophan hydroxylase 2 (TPH2), affecting an amygdala- and hippocampus-specific isoform. It was not associated with BPAD or MDD using a case-control association approach. CONCLUSIONS: We did not find evidence for a role of the TPH2 microdeletion in the pathoetiology of affective disorders. Further studies examining its putative role in behavioral traits regulated by the limbic system are warranted.

A comprehensive map of mobile element insertion polymorphisms in humans.

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.

Whole-exome sequencing links caspase recruitment domain 11 (CARD11) inactivation to severe combined immunodeficiency.

BACKGROUND: Primary immunodeficiencies represent model diseases for the mechanistic understanding of the human innate and adaptive immune response. They are clinically highly relevant per se because in patients with severe combined immunodeficiency (SCID), infections caused by opportunistic pathogens are typically life-threatening early in life. OBJECTIVES: We aimed at defining and functionally characterizing a novel form of SCID in an infant of consanguineous parents who presented with life-threatening Pneumocystis jirovecii pneumonia using a comprehensive immunologic and whole-exome genetic diagnostic strategy. METHODS: Analysis of leukocyte subpopulations was performed by using multicolor flow cytometry and was combined with stimulation tests for T-cell function. The search for a disease-causing mutation was performed with diagnostic whole-exome sequencing and systematic variant categorization. Reconstitution assays were used for validating the loss-of-function mutation. RESULTS: The novel entity of SCID was characterized by agammaglobulinemia and profoundly deficient T-cell function despite quantitatively normal T and B lymphocytes. Genetic analysis revealed a single pathogenic homozygous nonsense mutation of the caspase recruitment domain 11 (CARD11) gene. In reconstitution assays we demonstrated that the patient-derived truncated CARD11 protein is defective in antigen receptor signaling and nuclear factor κB activation. CONCLUSION: We show that an inactivating CARD11 mutation links defective nuclear factor κB signaling to a novel cause of autosomal recessive SCID.

DELLY: structural variant discovery by integrated paired-end and split-read analysis.

MOTIVATION: The discovery of genomic structural variants (SVs) at high sensitivity and specificity is an essential requirement for characterizing naturally occurring variation and for understanding pathological somatic rearrangements in personal genome sequencing data. Of particular interest are integrated methods that accurately identify simple and complex rearrangements in heterogeneous sequencing datasets at single-nucleotide resolution, as an optimal basis for investigating the formation mechanisms and functional consequences of SVs. RESULTS: We have developed an SV discovery method, called DELLY, that integrates short insert paired-ends, long-range mate-pairs and split-read alignments to accurately delineate genomic rearrangements at single-nucleotide resolution. DELLY is suitable for detecting copy-number variable deletion and tandem duplication events as well as balanced rearrangements such as inversions or reciprocal translocations. DELLY, thus, enables to ascertain the full spectrum of genomic rearrangements, including complex events. On simulated data, DELLY compares favorably to other SV prediction methods across a wide range of sequencing parameters. On real data, DELLY reliably uncovers SVs from the 1000 Genomes Project and cancer genomes, and validation experiments of randomly selected deletion loci show a high specificity. AVAILABILITY: DELLY is available at www.korbel.embl.de/software.html CONTACT: tobias.rausch@embl.de.

Fine mapping of a QTL on chromosome 13 for submaximal exercise capacity training response: the HERITAGE Family Study.

Although regular exercise improves submaximal aerobic capacity, there is large variability in its response to exercise training. While this variation is thought to be partly due to genetic differences, relatively little is known about the causal genes. Submaximal aerobic capacity traits in the current report include the responses of oxygen consumption (ΔVO(2)60), power output (ΔWORK60), and cardiac output (ΔQ60) at 60% of VO2max to a standardized 20-week endurance exercise training program. Genome-wide linkage analysis in 475 HERITAGE Family Study Caucasians identified a locus on chromosome 13q for ΔVO(2)60 (LOD = 3.11). Follow-up fine mapping involved a dense marker panel of over 1,800 single-nucleotide polymorphisms (SNPs) in a 7.9-Mb region (21.1-29.1 Mb from p-terminus). Single-SNP analyses found 14 SNPs moderately associated with both ΔVO(2)60 at P ≤ 0.005 and the correlated traits of ΔWORK60 and ΔQ60 at P < 0.05. Haplotype analyses provided several strong signals (P < 1.0 × 10(-5)) for ΔVO(2)60. Overall, association analyses narrowed the target region and included potential biological candidate genes (MIPEP and SGCG). Consistent with maximal heritability estimates of 23%, up to 20% of the phenotypic variance in ΔVO(2)60 was accounted for by these SNPs. These results implicate candidate genes on chromosome 13q12 for the ability to improve submaximal exercise capacity in response to regular exercise. Submaximal exercise at 60% of maximal capacity is an exercise intensity that falls well within the range recommended in the Physical Activity Guidelines for Americans and thus has potential public health relevance.