Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Allele-Specific Chromosome Removal after Cas9 Cleavage in Human Embryos.

Correction of disease-causing mutations in human embryos holds the potential to reduce the burden of inherited genetic disorders and improve fertility treatments for couples with disease-causing mutations in lieu of embryo selection. Here, we evaluate repair outcomes of a Cas9-induced double-strand break (DSB) introduced on the paternal chromosome at the EYS locus, which carries a frameshift mutation causing blindness. We show that the most common repair outcome is microhomology-mediated end joining, which occurs during the first cell cycle in the zygote, leading to embryos with non-mosaic restoration of the reading frame. Notably, about half of the breaks remain unrepaired, resulting in an undetectable paternal allele and, after mitosis, loss of one or both chromosomal arms. Correspondingly, Cas9 off-target cleavage results in chromosomal losses and hemizygous indels because of cleavage of both alleles. These results demonstrate the ability to manipulate chromosome content and reveal significant challenges for mutation correction in human embryos.

The genomes of all lungfish inform on genome expansion and tetrapod evolution.

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Structural variants drive context-dependent oncogene activation in cancer.

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer3. However, only a small minority of SVs are associated with altered gene expression4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

p53 directly represses human LINE1 transposons.

p53 is a potent tumor suppressor and commonly mutated in human cancers. Recently, we demonstrated that p53 genes act to restrict retrotransposons in germline tissues of flies and fish but whether this activity is conserved in somatic human cells is not known. Here we show that p53 constitutively restrains human LINE1s by cooperatively engaging sites in the 5'UTR and stimulating local deposition of repressive histone marks at these transposons. Consistent with this, the elimination of p53 or the removal of corresponding binding sites in LINE1s, prompted these retroelements to become hyperactive. Concurrently, p53 loss instigated chromosomal rearrangements linked to LINE sequences and also provoked inflammatory programs that were dependent on reverse transcriptase produced from LINE1s. Taken together, our observations establish that p53 continuously operates at the LINE1 promoter to restrict autonomous copies of these mobile elements in human cells. Our results further suggest that constitutive restriction of these retroelements may help to explain tumor suppression encoded by p53, since erupting LINE1s produced acute oncogenic threats when p53 was absent.

Small RNA-mediated genome rearrangement pathways in ciliates.

Most eukaryotes employ a combination of transcriptional and post-transcriptional silencing mechanisms to suppress transposons, yet ciliates employ a more extreme approach. They separate germline and somatic functions into distinct nuclei, enabling the elimination of transposons from the active somatic genome through diverse small RNA-mediated genome rearrangement pathways during sexual processes.

Patterns of somatic structural variation in human cancer genomes.

A key mutational process in cancer is structural variation, in which rearrangements delete, amplify or reorder genomic segments that range in size from kilobases to whole chromosomes1-7. Here we develop methods to group, classify and describe somatic structural variants, using data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumour types8. Sixteen signatures of structural variation emerged. Deletions have a multimodal size distribution, assort unevenly across tumour types and patients, are enriched in late-replicating regions and correlate with inversions. Tandem duplications also have a multimodal size distribution, but are enriched in early-replicating regions-as are unbalanced translocations. Replication-based mechanisms of rearrangement generate varied chromosomal structures with low-level copy-number gains and frequent inverted rearrangements. One prominent structure consists of 2-7 templates copied from distinct regions of the genome strung together within one locus. Such cycles of templated insertions correlate with tandem duplications, and-in liver cancer-frequently activate the telomerase gene TERT. A wide variety of rearrangement processes are active in cancer, which generate complex configurations of the genome upon which selection can act.

Complex genomic rearrangements: an underestimated cause of rare diseases.

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

The cryptic complex rearrangements involving the DMD gene: etiologic clues about phenotypical differences revealed by optical genome mapping.

BACKGROUND: Deletion or duplication in the DMD gene is one of the most common causes of Duchenne and Becker muscular dystrophy (DMD/BMD). However, the pathogenicity of complex rearrangements involving DMD, especially segmental duplications with unknown breakpoints, is not well understood. This study aimed to evaluate the structure, pattern, and potential impact of rearrangements involving DMD duplication. METHODS: Two families with DMD segmental duplications exhibiting phenotypical differences were recruited. Optical genome mapping (OGM) was used to explore the cryptic pattern of the rearrangements. Breakpoints were validated using long-range polymerase chain reaction combined with next-generation sequencing and Sanger sequencing. RESULTS: A multi-copy duplication involving exons 64-79 of DMD was identified in Family A without obvious clinical symptoms. Family B exhibited typical DMD neuromuscular manifestations and presented a duplication involving exons 10-13 of DMD. The rearrangement in Family A involved complex in-cis tandem repeats shown by OGM but retained a complete copy (reading frame) of DMD inferred from breakpoint validation. A reversed insertion with a segmental repeat was identified in Family B by OGM, which was predicted to disrupt the normal structure and reading frame of DMD after confirming the breakpoints. CONCLUSIONS: Validating breakpoint and rearrangement pattern is crucial for the functional annotation and pathogenic classification of genomic structural variations. OGM provides valuable insights into etiological analysis of DMD/BMD and enhances our understanding for cryptic effects of complex rearrangements.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities.

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Mechanistic origins of diverse genome rearrangements in cancer.

Cancer genomes frequently harbor structural chromosomal rearrangements that disrupt the linear DNA sequence order and copy number. To date, diverse classes of structural variants have been identified across multiple cancer types. These aberrations span a wide spectrum of complexity, ranging from simple translocations to intricate patterns of rearrangements involving multiple chromosomes. Although most somatic rearrangements are acquired gradually throughout tumorigenesis, recent interrogation of cancer genomes have uncovered novel categories of complex rearrangements that arises rapidly through a one-off catastrophic event, including chromothripsis and chromoplexy. Here we review the cellular and molecular mechanisms contributing to the formation of diverse structural rearrangement classes during cancer development. Genotoxic stress from a myriad of extrinsic and intrinsic sources can trigger DNA double-strand breaks that are subjected to DNA repair with potentially mutagenic outcomes. We also highlight how aberrant nuclear structures generated through mitotic cell division errors, such as rupture-prone micronuclei and chromosome bridges, can instigate massive DNA damage and the formation of complex rearrangements in cancer genomes.

A Cre-LoxP-based approach for combinatorial chromosome rearrangements in human HAP1 cells.

Alterations of human karyotype caused by chromosomal rearrangements are often associated with considerable phenotypic effects. Studying molecular mechanisms underlying these effects requires an efficient and scalable experimental model. Here, we propose a Cre-LoxP-based approach for the generation of combinatorial diversity of chromosomal rearrangements. We demonstrate that using the developed system, both intra- and inter-chromosomal rearrangements can be induced in the human haploid HAP1 cells, although the latter is significantly less effective. The obtained genetically modified HAP1 cell line can be used to dissect genomic effects associated with intra-chromosomal structural variations.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae.

BACKGROUND: Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti. METHODS AND RESULTS: Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae. CONCLUSIONS: Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

A rare KMT2A::CBL transcript in an acute monoblastic leukemia patient with an unfavorable outcome.

BACKGROUND: Lysine [K] methyltransferase 2A (KMT2A, previously known as MLL) gene rearrangements are common in acute leukemias of various lineages and are associated with features such as chemotherapy resistance and rapid relapse. KMT2A::CBL is a rare fusion of unknown pathogenesis generated by a unique interstitial deletion of chromosome 11 that has been reported across a wide age range in both acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL) patients. The leukemogenic effect of the KMT2A::CBL rearrangement and its association with clinical prognosis have not been well clarified. METHODS AND RESULTS: We report the case of a 64-year-old female who was diagnosed with acute monoblastic leukemia (M5a) and who acquired the rare KMT2A::CBL fusion. The patient received multiple cycles of therapy but did not achieve remission and eventually succumbed to severe infection and disease progression. Additionally, we characterized the predicted KMT2A-CBL protein structure in this case to reveal the underlying leukemogenic mechanisms and summarized reported cases of hematological malignancies with KMT2A::CBL fusion to investigate the correlation of gene rearrangements with clinical outcomes. CONCLUSIONS: This report provides novel insights into the leukemogenic potential of the KMT2A::CBL rearrangement and the correlation between gene rearrangements and clinical outcomes.

Malignant Perivascular epithelioid cell tumour of the uterus without TFE3 gene rearrangement: a case report.

BACKGROUND: Perivascular epithelioid cell tumours (PEComas) are soft tissue tumours. These neoplasms belong to the family of mesenchymal tumours, which include angiomyolipomas, clear-cell sugar tumours of the lung, and PEComas not otherwise specified (NOS). The probability of a perivascular epithelioid cell tumour (PEComa) occurring in the uterus is low, and the incidence, diagnosis, treatment, and outcomes of such tumours are still unclear. CASE PRESENTATION: A 51-year-old woman presented a 4-year history of natural menopause. An intrauterine mass was detected via ultrasound examination; the mass showed a tendency to increase but caused no symptoms. The levels of tumour markers were within the normal range. Pathological analysis of the curettage revealed perivascular epithelioid differentiation of the endometrial tumour. Consequently, a laparoscopic total hysterectomy with bilateral adnexectomy was performed. No distant metastasis was detected via whole-body positron emission computed tomography (PETCT) after the operation. Fluorescence in situ hybridization (FISH) revealed no TFE3 gene rearrangement. Next-generation sequencing of bone and soft tissue revealed negative TSC1/2 and TP53 expression. No recurrence or metastasis was observed during the 18-month follow-up period. CONCLUSION: PEComa of the gynecologic tract is a rare and challenging entity. Diffuse HMB-45 expression, TSC alterations and TFE3 rearrangement are characteristic of uterine PEComas. Surgical resection is the first choice. Genetic testing is helpful for determining the nature of the mass and for choosing targeted therapy. Further research is needed to establish treatment protocols.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Chromoanagenesis from radiation-induced genome damage in Populus.

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar (Populus sp.) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

R-loops and regulatory changes in chronologically ageing fission yeast cells drive non-random patterns of genome rearrangements.

Aberrant repair of DNA double-strand breaks can recombine distant chromosomal breakpoints. Chromosomal rearrangements compromise genome function and are a hallmark of ageing. Rearrangements are challenging to detect in non-dividing cell populations, because they reflect individually rare, heterogeneous events. The genomic distribution of de novo rearrangements in non-dividing cells, and their dynamics during ageing, remain therefore poorly characterized. Studies of genomic instability during ageing have focussed on mitochondrial DNA, small genetic variants, or proliferating cells. To characterize genome rearrangements during cellular ageing in non-dividing cells, we interrogated a single diagnostic measure, DNA breakpoint junctions, using Schizosaccharomyces pombe as a model system. Aberrant DNA junctions that accumulated with age were associated with microhomology sequences and R-loops. Global hotspots for age-associated breakpoint formation were evident near telomeric genes and linked to remote breakpoints elsewhere in the genome, including the mitochondrial chromosome. Formation of breakpoint junctions at global hotspots was inhibited by the Sir2 histone deacetylase and might be triggered by an age-dependent de-repression of chromatin silencing. An unexpected mechanism of genomic instability may cause more local hotspots: age-associated reduction in an RNA-binding protein triggering R-loops at target loci. This result suggests that biological processes other than transcription or replication can drive genome rearrangements. Notably, we detected similar signatures of genome rearrangements that accumulated in old brain cells of humans. These findings provide insights into the unique patterns and possible mechanisms of genome rearrangements in non-dividing cells, which can be promoted by ageing-related changes in gene-regulatory proteins.

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus.

Supernumerary mini-chromosomes-a unique type of genomic structural variation-have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

Evaluation of genetic risk of apparently balanced chromosomal rearrangement carriers by breakpoint characterization.

PURPOSE: To report genetic characteristics and associated risk of chromosomal breaks due to chromosomal rearrangements in large samples. METHODS: MicroSeq, a technique that combines chromosome microdissection and next-generation sequencing, was used to identify chromosomal breakpoints. Long-range PCR and Sanger sequencing were used to precisely characterize 100 breakpoints in 50 ABCR carriers. RESULTS: In addition to the recurrent regions of balanced rearrangement breaks in 8q24.13, 11q11.23, and 22q11.21 that had been documented, we have discovered a 10-Mb region of 12q24.13-q24.3 that could potentially be a sparse region of balanced rearrangement breaks. We found that 898 breakpoints caused gene disruption and a total of 188 breakpoints interrupted genes recorded in OMIM. The percentage of breakpoints that disrupted autosomal dominant genes recorded in OMIM was 25.53% (48/188). Fifty-four of the precisely characterized breakpoints had 1-8-bp microhomologous sequences. CONCLUSION: Our findings provide a reference for the evaluation of the pathogenicity of mutations in related genes that cause protein truncation in clinical practice. According to the characteristics of breakpoints, non-homologous end joining and microhomology-mediated break-induced replication may be the main mechanism for ABCRs formation.