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.

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.

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.

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.

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.

Comprehensive characterization of complex structural variations in cancer by directly comparing genome sequence reads.

The development of high-throughput sequencing technologies has advanced our understanding of cancer. However, characterizing somatic structural variants in tumor genomes is still challenging because current strategies depend on the initial alignment of reads to a reference genome. Here, we describe SMUFIN (somatic mutation finder), a single program that directly compares sequence reads from normal and tumor genomes to accurately identify and characterize a range of somatic sequence variation, from single-nucleotide variants (SNV) to large structural variants at base pair resolution. Performance tests on modeled tumor genomes showed average sensitivity of 92% and 74% for SNVs and structural variants, with specificities of 95% and 91%, respectively. Analyses of aggressive forms of solid and hematological tumors revealed that SMUFIN identifies breakpoints associated with chromothripsis and chromoplexy with high specificity. SMUFIN provides an integrated solution for the accurate, fast and comprehensive characterization of somatic sequence variation in cancer.

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.

Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma.

Pilocytic astrocytoma, the most common childhood brain tumor, is typically associated with mitogen-activated protein kinase (MAPK) pathway alterations. Surgically inaccessible midline tumors are therapeutically challenging, showing sustained tendency for progression and often becoming a chronic disease with substantial morbidities. Here we describe whole-genome sequencing of 96 pilocytic astrocytomas, with matched RNA sequencing (n = 73), conducted by the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. We identified recurrent activating mutations in FGFR1 and PTPN11 and new NTRK2 fusion genes in non-cerebellar tumors. New BRAF-activating changes were also observed. MAPK pathway alterations affected all tumors analyzed, with no other significant mutations identified, indicating that pilocytic astrocytoma is predominantly a single-pathway disease. Notably, we identified the same FGFR1 mutations in a subset of H3F3A-mutated pediatric glioblastoma with additional alterations in the NF1 gene. Our findings thus identify new potential therapeutic targets in distinct subsets of pilocytic astrocytoma and childhood glioblastoma.

Identification of a Ninein (NIN) mutation in a family with spondyloepimetaphyseal dysplasia with joint laxity (leptodactylic type)-like phenotype.

Spondyloepimetaphyseal dysplasia with joint laxity-leptodactylic type (SEMDJL2) is an autosomal dominant skeletal dysplasia which is characterized by midface hypoplasia, short stature, joint laxity with dislocations, genua valga, progressive scoliosis, and slender fingers. Recently, heterozygous missense mutations in KIF22, a gene which encodes a member of the kinesin-like protein family, have been identified in sporadic as well as familial cases of SEMDJL2. In the present study homozygosity mapping and whole-exome sequencing were combined to analyze a consanguineous family with a phenotype resembling SEMDJL2. We identified homozygous missense mutations in the two nearby genes NIN (Ninein) and POLE2 (DNA polymerase epsilon subunit B) which segregate with the disease in the family and were not present in 500 healthy control individuals and in the 1094 control individuals contained within the 1000-genomes database. We present several lines of evidence that mutant Ninein is most likely causative for the SEMDJL2-like phenotype. The centrosomal protein NIN shows a functional relationship with KIF22 and other proteins associated with chromosome congression/movement, centrosomal function, and ciliogenesis, which have been associated with skeletal dysplasias. Moreover, compound heterozygous missense mutations at more N-terminal positions of Ninein have very recently been identified in a family with microcephalic primordial dwarfism. Together with the present report this strongly supports a fundamental role of Ninein in skeletal development.

Integrative genomic analyses reveal an androgen-driven somatic alteration landscape in early-onset prostate cancer.

Early-onset prostate cancer (EO-PCA) represents the earliest clinical manifestation of prostate cancer. To compare the genomic alteration landscapes of EO-PCA with "classical" (elderly-onset) PCA, we performed deep sequencing-based genomics analyses in 11 tumors diagnosed at young age, and pursued comparative assessments with seven elderly-onset PCA genomes. Remarkable age-related differences in structural rearrangement (SR) formation became evident, suggesting distinct disease pathomechanisms. Whereas EO-PCAs harbored a prevalence of balanced SRs, with a specific abundance of androgen-regulated ETS gene fusions including TMPRSS2:ERG, elderly-onset PCAs displayed primarily non-androgen-associated SRs. Data from a validation cohort of > 10,000 patients showed age-dependent androgen receptor levels and a prevalence of SRs affecting androgen-regulated genes, further substantiating the activity of a characteristic "androgen-type" pathomechanism in EO-PCA.

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.