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.

Chromothripsis followed by circular recombination drives oncogene amplification in human cancer.

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. Using whole-genome sequencing data of the pediatric tumor neuroblastoma, we here identified a type of amplification, termed 'seismic amplification', that is characterized by multiple rearrangements and discontinuous copy number levels. Overall, seismic amplifications occurred in 9.9% (274 of 2,756) of cases across 38 cancer types, and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification showed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and overexpression in a substantial fraction of human malignancies.

Quantification of Translocation-Specific ctDNA Provides an Integrating Parameter for Early Assessment of Treatment Response and Risk Stratification in Ewing Sarcoma.

PURPOSE: We evaluated the predictive and prognostic value of circulating tumor DNA (ctDNA) in patients with Ewing sarcoma (EWS) treated in the EWING2008 trial. EXPERIMENTAL DESIGN: Plasma samples from 102 patients with EWS enrolled in the EWING2008 trial were obtained before and during induction chemotherapy. Genomic EWSR1 fusion sequence spanning primers and probes were used for highly specific and sensitive quantification of the levels of ctDNA by digital droplet PCR. ctDNA levels were correlated to established clinical risk factors and outcome parameters. RESULTS: Pretreatment ctDNA copy numbers were correlated with event-free and overall survival. The reduction in ctDNA levels below the detection limit was observed in most cases after only two blocks of vincristine, ifosfamide, doxorubicin, and etoposide (VIDE) induction chemotherapy. The persistence of ctDNA after two VIDE blocks was a strong predictor of poor outcomes. ctDNA levels correlated well with most established clinical risk factors; an inverse correlation was found only for the histologic response to induction therapy. ctDNA levels did not provide simple representations of tumor volume, but integrated information from various tumor characteristics represented an independent EWS tumor marker with predictive and prognostic value. CONCLUSIONS: ctDNA copy number in the plasma of patients with EWS is a quantifiable parameter for early risk stratification and can be used as a dynamic noninvasive biomarker for early prediction of treatment response and outcome of patients.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden.

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

A scalable CRISPR/Cas9-based fluorescent reporter assay to study DNA double-strand break repair choice.

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.

Genome-wide Screens Implicate Loss of Cullin Ring Ligase 3 in Persistent Proliferation and Genome Instability in TP53-Deficient Cells.

TP53 deficiency is the most common alteration in cancer; however, this alone is typically insufficient to drive tumorigenesis. To identify genes promoting tumorigenesis in combination with TP53 deficiency, we perform genome-wide CRISPR-Cas9 knockout screens coupled with proliferation and transformation assays in isogenic cell lines. Loss of several known tumor suppressors enhances cellular proliferation and transformation. Loss of neddylation pathway genes promotes uncontrolled proliferation exclusively in TP53-deficient cells. Combined loss of CUL3 and TP53 activates an oncogenic transcriptional program governed by the nuclear factor κB (NF-κB), AP-1, and transforming growth factor ß (TGF-ß) pathways. This program maintains persistent cellular proliferation, induces partial epithelial to mesenchymal transition, and increases DNA damage, genomic instability, and chromosomal rearrangements. Our findings reveal CUL3 loss as a key event stimulating persistent proliferation in TP53-deficient cells. These findings may be clinically relevant, since TP53-CUL3-deficient cells are highly sensitive to ataxia telangiectasia mutated (ATM) inhibition, exposing a vulnerability that could be exploited for cancer treatment.

Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories.

Identifying the earliest somatic changes in prostate cancer can give important insights into tumor evolution and aids in stratifying high- from low-risk disease. We integrated whole genome, transcriptome and methylome analysis of early-onset prostate cancers (diagnosis ≤55 years). Characterization across 292 prostate cancer genomes revealed age-related genomic alterations and a clock-like enzymatic-driven mutational process contributing to the earliest mutations in prostate cancer patients. Our integrative analysis identified four molecular subgroups, including a particularly aggressive subgroup with recurrent duplications associated with increased expression of ESRP1, which we validate in 12,000 tissue microarray tumors. Finally, we combined the patterns of molecular co-occurrence and risk-based subgroup information to deconvolve the molecular and clinical trajectories of prostate cancer from single patient samples.

Systematic Identification of Determinants for Single-Strand Annealing-Mediated Deletion Formation in Saccharomyces cerevisiae.

To ensure genomic integrity, living organisms have evolved diverse molecular processes for sensing and repairing damaged DNA. If improperly repaired, DNA damage can give rise to different types of mutations, an important class of which are genomic structural variants (SVs). In spite of their importance for phenotypic variation and genome evolution, potential contributors to SV formation in Saccharomyces cerevisiae (budding yeast), a highly tractable model organism, are not fully recognized. Here, we developed and applied a genome-wide assay to identify yeast gene knockout mutants associated with de novo deletion formation, in particular single-strand annealing (SSA)-mediated deletion formation, in a systematic manner. In addition to genes previously linked to genome instability, our approach implicates novel genes involved in chromatin remodeling and meiosis in affecting the rate of SSA-mediated deletion formation in the presence or absence of stress conditions induced by DNA-damaging agents. We closely examined two candidate genes, the chromatin remodeling gene IOC4 and the meiosis-related gene MSH4, which when knocked-out resulted in gene expression alterations affecting genes involved in cell division and chromosome organization, as well as DNA repair and recombination, respectively. Our high-throughput approach facilitates the systematic identification of processes linked to the formation of a major class of genetic variation.

Pan-cancer analysis of somatic copy-number alterations implicates IRS4 and IGF2 in enhancer hijacking.

Extensive prior research focused on somatic copy-number alterations (SCNAs) affecting cancer genes, yet the extent to which recurrent SCNAs exert their influence through rearrangement of cis-regulatory elements (CREs) remains unclear. Here we present a framework for inferring cancer-related gene overexpression resulting from CRE reorganization (e.g., enhancer hijacking) by integrating SCNAs, gene expression data and information on topologically associating domains (TADs). Analysis of 7,416 cancer genomes uncovered several pan-cancer candidate genes, including IRS4, SMARCA1 and TERT. We demonstrate that IRS4 overexpression in lung cancer is associated with recurrent deletions in cis, and we present evidence supporting a tumor-promoting role. We additionally pursued cancer-type-specific analyses and uncovered IGF2 as a target for enhancer hijacking in colorectal cancer. Recurrent tandem duplications intersecting with a TAD boundary mediate de novo formation of a 3D contact domain comprising IGF2 and a lineage-specific super-enhancer, resulting in high-level gene activation. Our framework enables systematic inference of CRE rearrangements mediating dysregulation in cancer.

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.

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.

Immortalization capacity of HPV types is inversely related to chromosomal instability.

High-risk human papillomavirus (hrHPV) types induce immortalization of primary human epithelial cells. Previously we demonstrated that immortalization of human foreskin keratinocytes (HFKs) is HPV type dependent, as reflected by the presence or absence of a crisis period before reaching immortality. This study determined how the immortalization capacity of ten hrHPV types relates to DNA damage induction and overall genomic instability in HFKs.Twenty five cell cultures obtained by transduction of ten hrHPV types (i.e. HPV16/18/31/33/35/45/51/59/66/70 E6E7) in two or three HFK donors each were studied.All hrHPV-transduced HFKs showed an increased number of double strand DNA breaks compared to controls, without exhibiting significant differences between types. However, immortal descendants of HPV-transduced HFKs that underwent a prior crisis period (HPV45/51/59/66/70-transduced HFKs) showed significantly more chromosomal aberrations compared to those without crisis (HPV16/18/31/33/35-transduced HFKs). Notably, the hTERT locus at 5p was exclusively gained in cells with a history of crisis and coincided with increased expression. Chromothripsis was detected in one cell line in which multiple rearrangements within chromosome 8 resulted in a gain of MYC.Together we demonstrated that upon HPV-induced immortalization, the number of chromosomal aberrations is inversely related to the viral immortalization capacity. We propose that hrHPV types with reduced immortalization capacity in vitro, reflected by a crisis period, require more genetic host cell aberrations to facilitate immortalization than types that can immortalize without crisis. This may in part explain the observed differences in HPV-type prevalence in cervical cancers and emphasizes that changes in the host cell genome contribute to HPV-induced carcinogenesis.

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.

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.

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.

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.

Genomics and drug profiling of fatal TCF3-HLF-positive acute lymphoblastic leukemia identifies recurrent mutation patterns and therapeutic options.

TCF3-HLF-positive acute lymphoblastic leukemia (ALL) is currently incurable. Using an integrated approach, we uncovered distinct mutation, gene expression and drug response profiles in TCF3-HLF-positive and treatment-responsive TCF3-PBX1-positive ALL. We identified recurrent intragenic deletions of PAX5 or VPREB1 in constellation with the fusion of TCF3 and HLF. Moreover somatic mutations in the non-translocated allele of TCF3 and a reduction of PAX5 gene dosage in TCF3-HLF ALL suggest cooperation within a restricted genetic context. The enrichment for stem cell and myeloid features in the TCF3-HLF signature may reflect reprogramming by TCF3-HLF of a lymphoid-committed cell of origin toward a hybrid, drug-resistant hematopoietic state. Drug response profiling of matched patient-derived xenografts revealed a distinct profile for TCF3-HLF ALL with resistance to conventional chemotherapeutics but sensitivity to glucocorticoids, anthracyclines and agents in clinical development. Striking on-target sensitivity was achieved with the BCL2-specific inhibitor venetoclax (ABT-199). This integrated approach thus provides alternative treatment options for this deadly disease.

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.

Analysis of deletion breakpoints from 1,092 humans reveals details of mutation mechanisms.

Investigating genomic structural variants at basepair resolution is crucial for understanding their formation mechanisms. We identify and analyse 8,943 deletion breakpoints in 1,092 samples from the 1000 Genomes Project. We find breakpoints have more nearby SNPs and indels than the genomic average, likely a consequence of relaxed selection. By investigating the correlation of breakpoints with DNA methylation, Hi-C interactions, and histone marks and the substitution patterns of nucleotides near them, we find that breakpoints with the signature of non-allelic homologous recombination (NAHR) are associated with open chromatin. We hypothesize that some NAHR deletions occur without DNA replication and cell division, in embryonic and germline cells. In contrast, breakpoints associated with non-homologous (NH) mechanisms often have sequence microinsertions, templated from later replicating genomic sites, spaced at two characteristic distances from the breakpoint. These microinsertions are consistent with template-switching events and suggest a particular spatiotemporal configuration for DNA during the events.