Genomic and functional overlap between somatic and germline chromosomal rearrangements

Integrative omics approaches to advance rare disease diagnostics

Over the past decade high-throughput DNA sequencing approaches, namely whole exome and whole genome sequencing became a standard procedure in Mendelian disease diagnostics. Implementation of these technologies greatly facilitated diagnostics and shifted the analysis paradigm from variant identification to prioritisation and evaluation. The diagnostic rates vary widely depending on the cohort size, heterogeneity and disease and range from around 30% to 50% leaving the majority of patients undiagnosed. Advances in omics technologies and computational analysis provide an opportunity to increase these unfavourable rates by providing evidence for disease-causing variant validation and prioritisation. This review aims to provide an overview of the current application of several omics technologies including RNA-sequencing, proteomics, metabolomics and DNA-methylation profiling for diagnostics of rare genetic diseases in general and inborn errors of metabolism in particular.

Detection of aberrant gene expression events in RNA sequencing data

RNA sequencing (RNA-seq) has emerged as a powerful approach to discover disease-causing gene regulatory defects in individuals affected by genetically undiagnosed rare disorders. Pioneering studies have shown that RNA-seq could increase the diagnosis rates over DNA sequencing alone by 8-36%, depending on the disease entity and tissue probed. To accelerate adoption of RNA-seq by human genetics centers, detailed analysis protocols are now needed. We present a step-by-step protocol that details how to robustly detect aberrant expression levels, aberrant splicing and mono-allelic expression in RNA-seq data using dedicated statistical methods. We describe how to generate and assess quality control plots and interpret the analysis results. The protocol is based on the detection of RNA outliers pipeline (DROP), a modular computational workflow that integrates all the analysis steps, can leverage parallel computing infrastructures and generates browsable web page reports.

Copy number variations and constitutional chromothripsis (Review)

Both copy number variations (CNVs) and chromothripsis are phenomena that involve complex genomic rearrangements. Chromothripsis results in CNVs and other structural changes. CNVs are frequently observed in the human genome. Studies on CNVs have been increasing exponentially; the Database of Genomic Variants shows an increase in the number of data published on structural variations added to the database in the last 15 years. CNVs may be a result of replicative and non-replicative mechanisms, and are hypothesized to serve important roles in human health and disease. Chromothripsis is a phenomena of chromosomal rearrangement following chromosomal breaks at multiple locations and involves impaired DNA repair. In 2011, Stephens et al coined the term chromothripsis for this type of fragmenting event. Several proposed mechanisms have been suggested to underlie chromothripsis, such as p53 inactivation, micronuclei formation, abortive apoptosis and telomere fusions in telomere crisis. Chromothripsis gives rise to normal or abnormal phenotypes. In this review, constitutional chromothripsis, which may coexist with multiple de novo CNVs are described and discussed. This reviews aims to summarize recent advances in our understanding of CNVs and chromothripsis, and describe the effects of these phenomena on human health and birth defects.

SMaSH: Sample matching using SNPs in humans

BACKGROUND

Inadvertent sample swaps are a real threat to data quality in any medium to large scale omics studies. While matches between samples from the same individual can in principle be identified from a few well characterized single nucleotide polymorphisms (SNPs), omics data types often only provide low to moderate coverage, thus requiring integration of evidence from a large number of SNPs to determine if two samples derive from the same individual or not.

METHODS

We select about six thousand SNPs in the human genome and develop a Bayesian framework that is able to robustly identify sample matches between next generation sequencing data sets.

RESULTS

We validate our approach on a variety of data sets. Most importantly, we show that our approach can establish identity between different omics data types such as Exome, RNA-Seq, and MethylCap-Seq. We demonstrate how identity detection degrades with sample quality and read coverage, but show that twenty million reads of a fairly low quality RNA-Seq sample are still sufficient for reliable sample identification.

CONCLUSION

Our tool, SMASH, is able to identify sample mismatches in next generation sequencing data sets between different sequencing modalities and for low quality sequencing data.

Investigation of somatic CNVs in brains of synucleinopathy cases using targeted SNCA analysis and single cell sequencing

Synucleinopathies are mostly sporadic neurodegenerative disorders of partly unexplained aetiology, and include Parkinson's disease (PD) and multiple system atrophy (MSA). We have further investigated our recent finding of somatic SNCA (α-synuclein) copy number variants (CNVs, specifically gains) in synucleinopathies, using Fluorescent in-situ Hybridisation for SNCA, and single-cell whole genome sequencing for the first time in a synucleinopathy. In the cingulate cortex, mosaicism levels for SNCA gains were higher in MSA and PD than controls in neurons (> 2% in both diseases), and for MSA also in non-neurons. In MSA substantia nigra (SN), we noted SNCA gains in > 3% of dopaminergic (DA) neurons (identified by neuromelanin) and neuromelanin-negative cells, including olig2-positive oligodendroglia. Cells with CNVs were more likely to have α-synuclein inclusions, in a pattern corresponding to cell categories mostly relevant to the disease: DA neurons in Lewy-body cases, and other cells in the striatonigral degeneration-dominant MSA variant (MSA-SND). Higher mosaicism levels in SN neuromelanin-negative cells may correlate with younger onset in typical MSA-SND, and in cingulate neurons with younger death in PD. Larger sample sizes will, however, be required to confirm these putative findings. We obtained genome-wide somatic CNV profiles from 169 cells from the substantia nigra of two MSA cases, and pons and putamen of one. These showed somatic CNVs in ~ 30% of cells, with clonality and origins in segmental duplications for some. CNVs had distinct profiles based on cell type, with neurons having a mix of gains and losses, and other cells having almost exclusively gains, although control data sets will be required to determine possible disease relevance. We propose that somatic SNCA CNVs may contribute to the aetiology and pathogenesis of synucleinopathies, and that genome-wide somatic CNVs in MSA brain merit further study.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease

BACKGROUND

RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.

METHODS

We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.

RESULTS

We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.

CONCLUSIONS

The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.

On the Complexity of Mechanisms and Consequences of Chromothripsis: An Update

In the present review, we focus on the phenomenon of chromothripsis, a new type of complex chromosomal rearrangements. We discuss the challenges of chromothripsis detection and its distinction from other chromoanagenesis events. Along with already known causes and mechanisms, we introduce aberrant epigenetic regulation as a possible pathway to chromothripsis. We address the issue of chromothripsis characteristics in cancers and benign tumours, as well as chromothripsis inheritance in cases of its occurrence in germ cells, zygotes and early embryos. Summarising the presented data on different phenotypic effect of chromothripsis, we assume that its consequences are most likely determined not by the chromosome shattering and reassembly themselves, but by the genome regions involved in the rearrangement.

Evolutionary Loss of Genomic Proximity to Conserved Noncoding Elements Impacted the Gene Expression Dynamics During Mammalian Brain Development

Conserved noncoding elements (CNEs) have a significant regulatory influence on their neighboring genes. Loss of proximity to CNEs through genomic rearrangements can, therefore, impact the transcriptional states of the cognate genes. Yet, the evolutionary implications of such chromosomal alterations have not been studied. Through genome-wide analysis of CNEs and the cognate genes of representative species from five different mammalian orders, we observed a significant loss of genes' linear proximity to CNEs in the rat lineage. The CNEs and the genes losing proximity had a significant association with fetal, but not postnatal, brain development as assessed through ontology terms, developmental gene expression, chromatin marks, and genetic mutations. The loss of proximity to CNEs correlated with the independent evolutionary loss of fetus-specific upregulation of nearby genes in the rat brain. DNA breakpoints implicated in brain abnormalities of germline origin had significant representation between a CNE and the gene that exhibited loss of proximity, signifying the underlying developmental tolerance of genomic rearrangements that allowed the evolutionary splits of CNEs and the cognate genes in the rodent lineage. Our observations highlighted a nontrivial impact of chromosomal rearrangements in shaping the evolutionary dynamics of mammalian brain development and might explain the loss of brain traits, like cerebral folding of the cortex, in the rodent lineage.

Homologous Recombination and the Formation of Complex Genomic Rearrangements

The maintenance of genome integrity involves multiple independent DNA damage avoidance and repair mechanisms. However, the origin and pathways of the focal chromosomal reshuffling phenomena collectively referred to as chromothripsis remain mechanistically obscure. We discuss here the role, mechanisms, and regulation of homologous recombination (HR) in the formation of simple and complex chromosomal rearrangements. We emphasize features of the recently characterized multi-invasion (MI)-induced rearrangement (MIR) pathway which uniquely amplifies the initial DNA damage. HR intermediates and cellular contexts that endanger genomic stability are discussed as well as the emerging roles of various classes of nucleases in the formation of genome rearrangements. Long-read sequencing and improved mapping of repeats should enable better appreciation of the significance of recombination in generating genomic rearrangements.

Genetic diagnosis of a Chinese multiple endocrine neoplasia type 2A family through whole genome sequencing

Approximately 98% of patients with multiple endocrine neoplasia type 2A (MEN 2A) have an identifiable RET mutation. Prophylactic or early total thyroidectomy or pheochromocytoma/parathyroid removal in patients can be preventative or curative and has become standard management. The general strategy for RET screening on family members at risk is to sequence the most commonly affected exons and, if negative, to extend sequencing to additional exons. However, different families with MEN 2A due to the same RET mutation often have significant variability in the clinical exhibition of disease and aggressiveness of the MTC, which implies additional genetic loci exsit beyond RET coding region. Whole genome sequencing (WGS) greatly expands the breadth of screening from genes associated with a particular disease to the whole genome and, potentially, all the information that the genome contains about diseases or traits. This is presumably due to additive effect of disease modifying factors. In this study, we performed WGS on a typical Chinese MEN 2A proband and identified the pathogenic RET p.C634R mutation. We also identified several neutral variants within RET and pheochromocytoma-related genes. Moreover, we found several interesting structural variants including genetic deletions (RSPO1, OVCH2 and AP3S1, etc.) and fusion transcripts (FSIP1-BAZ2A, etc.).

Translocation Breakpoints Preferentially Occur in Euchromatin and Acrocentric Chromosomes

Chromosomal translocations drive the development of many hematological and some solid cancers. Several factors have been identified to explain the non-random occurrence of translocation breakpoints in the genome. These include chromatin density, gene density and CCCTC-binding factor (CTCF)/cohesin binding site density. However, such factors are at least partially interdependent. Using 13,844 and 1563 karyotypes from human blood and solid cancers, respectively, our multiple regression analysis only identified chromatin density as the primary statistically significant predictor. Specifically, translocation breakpoints preferentially occur in open chromatin. Also, blood and solid tumors show markedly distinct translocation signatures. Strikingly, translocation breakpoints occur significantly more frequently in acrocentric chromosomes than in non-acrocentric chromosomes. Thus, translocations are probably often generated around nucleoli in the inner nucleoplasm, away from the nuclear envelope. Importantly, our findings remain true both in multivariate analyses and after removal of highly recurrent translocations. Finally, we applied pairwise probabilistic co-occurrence modeling. In addition to well-known highly prevalent translocations, such as those resulting in BCR-ABL1 (BCR-ABL) and RUNX1-RUNX1T1 (AML1-ETO) fusion genes, we identified significantly underrepresented translocations with putative fusion genes, which are probably subject to strong negative selection during tumor evolution. Taken together, our findings provide novel insights into the generation and selection of translocations during cancer development.

Distinct distributions of genomic features of the 5’ and 3’ partners of coding somatic cancer gene fusions: arising mechanisms and functional implications

The genomic features and arising mechanisms of coding cancer somatic gene fusions (CSGFs) largely remain elusive. In this study, we show the gene origin stratification pattern of CSGF partners that fusion partners in human cancers are significantly enriched for genes with the gene age ofEuteleostomes and with the gene family age of Bilateria. GC skew (a measurement of G, C nucleotide content bias, (G-C)/(G+C)) is a useful measurement to indicate the DNA leading strand, lagging strand, replication origin, and replication terminal and DNA-RNA R-loop formation. We find that GC skew bias at the 5 prime (5′) but not the 3 prime (3’) partners of CSGFs, coincident with the polarity feature of gene expression breadth that the 5’ partners are more ubiquitous while the 3’ fusion partners are more tissue specific in general. We reveal distinct length and composition distributions of 5’ and 3’ of CSGFs, including sequence features corresponded to the 5’ untranslated regions (UTRs), 3’ UTRs, and the N-terminal sequences of the encoded proteins. Oncogenic somatic gene fusions are most enriched for the 5’ and 3’ genes’ somatic amplification alongside a substantial proportion of other types of combinations. At the function level, 5’ partners of CSGFs appear more likely to be tumour suppressor genes while many 3’ partners appear to be proto-oncogene. Such distinct polarities of CSGFs at the evolutionary, structural, genomic and functional levels indicate the heterogeneous arsing mechanisms of CSGFs including R-loops and suggest potential novel targeted therapeutics specific to CSGF functional categories.

Molecular dissection of germline chromothripsis in a developmental context using patient-derived iPS cells

Background

Germline chromothripsis causes complex genomic rearrangements that are likely to affect multiple genes and their regulatory contexts. The contribution of individual rearrangements and affected genes to the phenotypes of patients with complex germline genomic rearrangements is generally unknown.

Methods

To dissect the impact of germline chromothripsis in a relevant developmental context, we performed trio-based RNA expression analysis on blood cells, induced pluripotent stem cells (iPSCs), and iPSC-derived neuronal cells from a patient with de novo germline chromothripsis and both healthy parents. In addition, Hi-C and 4C-seq experiments were performed to determine the effects of the genomic rearrangements on transcription regulation of genes in the proximity of the breakpoint junctions.

Results

Sixty-seven genes are located within 1 Mb of the complex chromothripsis rearrangements involving 17 breakpoints on four chromosomes. We find that three of these genes (FOXP1, DPYD, and TWIST1) are both associated with developmental disorders and differentially expressed in the patient. Interestingly, the effect on TWIST1 expression was exclusively detectable in the patient’s iPSC-derived neuronal cells, stressing the need for studying developmental disorders in the biologically relevant context. Chromosome conformation capture analyses show that TWIST1 lost genomic interactions with several enhancers due to the chromothripsis event, which likely led to deregulation of TWIST1 expression and contributed to the patient’s craniosynostosis phenotype.

Conclusions

We demonstrate that a combination of patient-derived iPSC differentiation and trio-based molecular profiling is a powerful approach to improve the interpretation of pathogenic complex genomic rearrangements. Here we have applied this approach to identify misexpression of TWIST1, FOXP1, and DPYD as key contributors to the complex congenital phenotype resulting from germline chromothripsis rearrangements.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0399-z) contains supplementary material, which is available to authorized users.

Rare congenital chromosomal aberration dic(X;Y)(p22.33;p11.32) in a patient with primary myelofibrosis

BACKGROUND

Constitutional translocations between sex chromosomes are rather rare in humans with breakpoints at Xp11 and Yq11 as the most frequent. Breakpoints on the short arm of the Y chromosome form one subgroup of t(X;Y), giving rise to a derived chromosome with the centromeres of both the X and Y chromosomes, dic(X;Y). Here, we report a rare congenital chromosomal aberration, 46,X,dic(X;Y)(p22.33;p11.32)[20]/45,X[10], in an adult male.

CASE PRESENTATION

Primary myelofibrosis, a malignant haematological disease, was diagnosed in a 63-year-old man following liver transplantation after hepatocellular carcinoma. By the analysis of the bone marrow sample, the karyotype 46,X,dic(X;Y)(p22.33;p11.32) was detected in all the mitoses analysed and verified with multicolour fluorescence in situ hybridization (mFISH). A cytogenetic examination of stimulated peripheral blood cells revealed the constitutional karyotype 46,X,dic(X;Y)(p22.33;p11.32)[20]/45,X[10]. The cell line 45,X was confirmed with FISH in 35 % of interphase nuclei. The SRY locus was present on the dicentric chromosome. A CGH/SNP array (Illumina) revealed a gain of 153,7 Mbp of the X chromosome and a 803-kbp microdeletion (including the SHOX gene), which were also confirmed with FISH. SHOX encodes a transcriptional factor that regulates the growth of the long bones. The deletion of the SHOX gene together with the Madelung deformity of the forearm and the short stature of the proband led to a diagnosis of Léri-Weill dyschondrosteosis (LWD). The gain of almost the whole X chromosome (153,7 Mbp) was considered a variant of Klinefelter syndrome (KS). The levels of gonadotropins and testosterone were consistent with gonadal dysfunction. A malformation of the right external ear was detected.

CONCLUSIONS

We have reported a structural aberration of the sex chromosomes, dic(X;Y)(p22.33;p11.32). The related genomic imbalance is associated with two known hereditary syndromes, LWD and a KS variant, identified in our proband at an advanced age. Because the breakpoints did not involve cancer genes, we inferred that the two malignancies in the proband were not caused by this abnormality. The possible influence of SHOX haploinsufficiency on the growth regulation of auricular chondrocytes is discussed.

The genomic characteristics and cellular origin of chromothripsis

Human genomes are continuously subjected to mutations, which can drive genetic diseases and cancer. An intriguing recent finding has been the discovery of chromothripsis, a spectacular and complex form of chromosome rearrangement that can occur in the genomes of cancer cells and patients with congenital diseases. Chromothripsis has been described in a large array of human cancers and various types of chromothripsis have appeared, which differ in complexity and genomic hallmarks. From the combined genomic data a consensus hypothesis has been inferred, involving aberrant DNA replication and chromosome shattering as the underlying processes explaining chromothripsis. In addition, recent work has established several cellular models that recapitulate chromothripsis under defined experimental conditions. One of these models indicates that chromothripsis can originate from DNA damage in micronuclei, providing an elegant explanation for the restriction of chromothriptic rearrangements to a single chromosome. Alternatively, chromothripsis can be caused by telomere crisis, a process that involves formation of dicentric chromosomes and chromatin bridges. Here, we summarize the genomic features of chromothripsis and we discuss experimental approaches that allow dissection of the chromothripsis process.

Mechanisms of Origin, Phenotypic Effects and Diagnostic Implications of Complex Chromosome Rearrangements

Complex chromosome rearrangements (CCRs) are currently defined as structural genome variations that involve more than 2 chromosome breaks and result in exchanges of chromosomal segments. They are thought to be extremely rare, but their detection rate is rising because of improvements in molecular cytogenetic technology. Their population frequency is also underestimated, since many CCRs may not elicit a phenotypic effect. CCRs may be the result of fork stalling and template switching, microhomology-mediated break-induced repair, breakage-fusion-bridge cycles, or chromothripsis. Patients with chromosomal instability syndromes show elevated rates of CCRs due to impaired DNA double-strand break responses during meiosis. Therefore, the putative functions of the proteins encoded by ATM, BLM, WRN, ATR, MRE11, NBS1, and RAD51 in preventing CCRs are discussed. CCRs may exert a pathogenic effect by either (1) gene dosage-dependent mechanisms, e.g. haploinsufficiency, (2) mechanisms based on disruption of the genomic architecture, such that genes, parts of genes or regulatory elements are truncated, fused or relocated and thus their interactions disturbed - these mechanisms will predominantly affect gene expression - or (3) mixed mutation mechanisms in which a CCR on one chromosome is combined with a different type of mutation on the other chromosome. Such inferred mechanisms of pathogenicity need corroboration by mRNA sequencing. Also, future studies with in vitro models, such as inducible pluripotent stem cells from patients with CCRs, and transgenic model organisms should substantiate current inferences regarding putative pathogenic effects of CCRs. The ramifications of the growing body of information on CCRs for clinical and experimental genetics and future treatment modalities are briefly illustrated with 2 cases, one of which suggests KDM4C (JMJD2C) as a novel candidate gene for mental retardation.

Human Structural Variation: Mechanisms of Chromosome Rearrangements

Chromosome structural variation (SV) is a normal part of variation in the human genome, but some classes of SV can cause neurodevelopmental disorders. Analysis of the DNA sequence at SV breakpoints can reveal mutational mechanisms and risk factors for chromosome rearrangement. Large-scale SV breakpoint studies have become possible recently owing to advances in next-generation sequencing (NGS) including whole-genome sequencing (WGS). These findings have shed light on complex forms of SV such as triplications, inverted duplications, insertional translocations, and chromothripsis. Sequence-level breakpoint data resolve SV structure and determine how genes are disrupted, fused, and/or misregulated by breakpoints. Recent improvements in breakpoint sequencing have also revealed non-allelic homologous recombination (NAHR) between paralogous long interspersed nuclear element (LINE) or human endogenous retrovirus (HERV) repeats as a cause of deletions, duplications, and translocations. This review covers the genomic organization of simple and complex constitutional SVs, as well as the molecular mechanisms of their formation.

Unbalanced translocations arise from diverse mutational mechanisms including chromothripsis

Unbalanced translocations are a relatively common type of copy number variation and a major contributor to neurodevelopmental disorders. We analyzed the breakpoints of 57 unique unbalanced translocations to investigate the mechanisms of how they form. Fifty-one are simple unbalanced translocations between two different chromosome ends, and six rearrangements have more than three breakpoints involving two to five chromosomes. Sequencing 37 breakpoint junctions revealed that simple translocations have between 0 and 4 base pairs (bp) of microhomology (n = 26), short inserted sequences (n = 8), or paralogous repeats (n = 3) at the junctions, indicating that translocations do not arise primarily from nonallelic homologous recombination but instead form most often via nonhomologous end joining or microhomology-mediated break-induced replication. Three simple translocations fuse genes that are predicted to produce in-frame transcripts of SIRPG-WWOX, SMOC2-PROX1, and PIEZO2-MTA1, which may lead to gain of function. Three complex translocations have inversions, insertions, and multiple breakpoint junctions between only two chromosomes. Whole-genome sequencing and fluorescence in situ hybridization analysis of two de novo translocations revealed at least 18 and 33 breakpoints involving five different chromosomes. Breakpoint sequencing of one maternally inherited translocation involving four chromosomes uncovered multiple breakpoints with inversions and insertions. All of these breakpoint junctions had 0–4 bp of microhomology consistent with chromothripsis, and both de novo events occurred on paternal alleles. Together with other studies, these data suggest that germline chromothripsis arises in the paternal genome and may be transmitted maternally. Breakpoint sequencing of our large collection of chromosome rearrangements provides a comprehensive analysis of the molecular mechanisms behind translocation formation.