An Organismal CNV Mutator Phenotype Restricted to Early Human Development

Isodicentric Y Chromosome with Multiple Breakpoints in the Pseudoautosomal Region 1

INTRODUCTION

Isodicentric Y chromosomes are relatively common structural variants of the human genome. The underlying mechanism of isodicentric Y chromosomes with short arm breakpoints [idic(Yq)] remains to be clarified.

CASE PRESENTATION

We encountered a Japanese man with azoospermia and mild short stature. G-banding and array-based comparative genomic hybridization indicated that his karyotype was 45,X/46,X,idic(Y)(qter→p11.32::p11.32→qter) with a ∼1.8 Mb terminal deletion. Whole-genome sequencing suggested that the Y chromosome had four breakpoints in a ∼7 kb region of the pseudoautosomal region 1 (PAR1).

CONCLUSION

This case was assumed to have an idic(Yq) resulting from multiple DNA double-strand breaks in PAR1. This rearrangement may have been facilitated by the PAR1-specific chromatin architecture. The clinical features of the patient can be ascribed to SHOX haploinsufficiency and the presence of a 45,X cell line, although copy-number gains of some Yq genes and the size reduction of PAR1 may also contribute to his spermatogenic failure.

SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemia

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C. Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5 and 37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the hypothesis of a rare phenomenon of hypermutation. Eight-hundred-thirty apparent de novo SNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component having an average VAF of 25% (range: 18.7-37.41%) (n = 446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range 9.55-18.69%) (n = 384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA. We propose that other FA groups should be investigated for genome-wide de novo variants.

SNV/indel hypermutator phenotype in biallelic RAD51C variant - Fanconi anemia

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C . Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5%-37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the rare phenomenon of hypermutation. Eight-hundred-thirty apparent dnSNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component with an average VAF of 25% (range: 18.7-37.41%) (n=446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range: 9.55-18.69%) (n=384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA, and pregnancy loss in general. We propose that other FA groups should be investigated for genome-wide de novo variants.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

Emerging technologies for prenatal diagnosis: The application of whole genome and RNA sequencing

DNA sequencing technologies for clinical genetic testing have been rapidly evolving in recent years, and steadily become more important within the field of prenatal diagnostics. This review aims to give an overview of recent developments and to describe how they have the potential to fill the gaps of the currently clinically implemented methods for prenatal diagnosis of various genetic disorders. It has been shown for postnatal testing that whole genome sequencing provides a set of added benefits compared to exome sequencing, and it is to be expected that this will be the case for prenatal testing as well. RNA-sequencing, already used postnatally, can provide valuable complementary data to DNA-based testing, and aid in variant interpretation. While not ready for clinical implementation, emerging technologies such as long-read and Hi-C sequencing analyses might add to the toolbox for interpreting the expanding genetic data sets generated by genome-wide sequencing. Lastly, we also discuss some more practical implications of introducing these emerging technologies, which generate larger and larger genomic data sets, in the prenatal field.

Genetic and chemotherapeutic influences on germline hypermutation

Mutations in the germline generates all evolutionary genetic variation and is a cause of genetic disease. Parental age is the primary determinant of the number of new germline mutations in an individual's genome1,2. Here we analysed the genome-wide sequences of 21,879 families with rare genetic diseases and identified 12 individuals with a hypermutated genome with between two and seven times more de novo single-nucleotide variants than expected. In most families (9 out of 12), the excess mutations came from the father. Two families had genetic drivers of germline hypermutation, with fathers carrying damaging genetic variation in DNA-repair genes. For five of the families, paternal exposure to chemotherapeutic agents before conception was probably a key driver of hypermutation. Our results suggest that the germline is well protected from mutagenic effects, hypermutation is rare, the number of excess mutations is relatively modest and most individuals with a hypermutated genome will not have a genetic disease.

Cancer Genomic Rearrangements and Copy Number Alterations from Errors in Cell Division

Analysis of cancer genomes has shown that a large fraction of chromosomal changes originate from catastrophic events including whole-genome duplication, chromothripsis, breakage-fusion-bridge cycles, and chromoplexy. Through sophisticated computational analysis of cancer genomes and experimental recapitulation of these catastrophic alterations, we have gained significant insights into the origin, mechanism, and evolutionary dynamics of cancer genome complexity. In this review, we summarize this progress and survey the major unresolved questions, with particular emphasis on the relative contributions of chromosome fragmentation and DNA replication errors to complex chromosomal alterations.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.

Cancer cells are highly susceptible to accumulation of templated insertions linked to MMBIR

Microhomology-mediated break-induced replication (MMBIR) is a DNA repair pathway initiated by polymerase template switching at microhomology, which can produce templated insertions that initiate chromosomal rearrangements leading to neurological and metabolic diseases, and promote complex genomic rearrangements (CGRs) found in cancer. Yet, how often templated insertions accumulate from processes like MMBIR in genomes is poorly understood due to difficulty in directly identifying these events by whole genome sequencing (WGS). Here, by using our newly developed MMBSearch software, we directly detect such templated insertions (MMB-TIs) in human genomes and report substantial differences in frequency and complexity of MMB-TI events between normal and cancer cells. Through analysis of 71 cancer genomes from The Cancer Genome Atlas (TCGA), we observed that MMB-TIs readily accumulate de novo across several cancer types, with particularly high accumulation in some breast and lung cancers. By contrast, MMB-TIs appear only as germline variants in normal human fibroblast cells, and do not accumulate as de novo somatic mutations. Finally, we performed WGS on a lung adenocarcinoma patient case and confirmed MMB-TI-initiated chromosome fusions that disrupted potential tumor suppressors and induced chromothripsis-like CGRs. Based on our findings we propose that MMB-TIs represent a trigger for widespread genomic instability and tumor evolution.

Genome maintenance during embryogenesis

Genome maintenance during embryogenesis is critical, because defects during this period can be perpetuated and thus have a long-term impact on individual's health and longevity. Nevertheless, genome instability is normal during certain aspects of embryonic development, indicating that there is a balance between the exigencies of timely cell proliferation and mutation prevention. In particular, early embryos possess unique cellular and molecular features that underscore the challenge of having an appropriate balance. Here, we discuss genome instability during embryonic development, the mechanisms used in various cell compartments to manage genomic stress and address outstanding questions regarding the balance between genome maintenance mechanisms in key cell types that are important for adulthood and progeny.

Double-strand breaks induce short-scale DNA replication and damage amplification in the fully grown mouse oocytes

Break-induced replication (BIR) is essential for the repair of DNA double-strand breaks (DSBs) with single ends. DSBs-induced microhomology-mediated BIR (mmBIR) and template-switching can increase the risk of complex genome rearrangement. In addition, DSBs can also induce the multi-invasion-mediated DSB amplification. The mmBIR-induced genomic rearrangement has been identified in cancer cells and patients with rare diseases. However, when and how mmBIR is initiated have not been fully and deeply studied. Furthermore, it is not well understood about the conditions for initiation of multi-invasion-mediated DSB amplification. In the G2 phase oocyte of mouse, we identified a type of short-scale BIR (ssBIR) using the DNA replication indicator 5-ethynyl-2'-deoxyuridine (EdU). These ssBIRs could only be induced in the fully grown oocytes but not the growing oocytes. If the DSB oocytes were treated with Rad51 or Chek1/2 inhibitors, both EdU signals and DSB marker γH2A.X foci would decrease. In addition, the DNA polymerase inhibitor Aphidicolin could inhibit the ssBIR and another inhibitor ddATP could reduce the number of γH2A.X foci in the DSB oocytes. In conclusion, our results showed that DNA DSBs in the fully grown oocytes can initiate ssBIR and be amplified by Rad51 or DNA replication.

Somatic copy number variants in neuropsychiatric disorders

Copy number variants (CNVs) have been implicated in neuropsychiatric disorders, with rare-inherited and de novo CNVs (dnCNVs) having large effects on disease liability. Recent studies started exploring a class of dnCNVs that occur post-zygotically, and are therefore present in some but not all cells of the body. Analogous to conditional mutations in animal models, the presence of risk mutations in a fraction of cells has the potential to enlighten how damaging mutations affect cell-type/cell-circuit specific pathologies leading to neuropsychiatric manifestations. Although mosaic CNVs appear to contribute to a modest fraction of risk (0.3-0.5%), expanding our insights about them with more sensitive experimental and statistical methods, has the potential to help clarify mechanisms of neuropsychiatric disease.

Genome diversity and instability in human germ cells and preimplantation embryos

Genome diversity is essential for evolution and is of fundamental importance to human health. Generating genome diversity requires phases of DNA damage and repair that can cause genome instability. Humans have a high incidence of de novo congenital disorders compared to other organisms. Recent access to eggs, sperm and preimplantation embryos is revealing unprecedented rates of genome instability that may result in infertility and de novo mutations that cause genomic imbalance in at least 70% of conceptions. The error type and incidence of de novo mutations differ during developmental stages and are influenced by differences in male and female meiosis. In females, DNA repair is a critical factor that determines fertility and reproductive lifespan. In males, aberrant meiotic recombination causes infertility, embryonic failure and pregnancy loss. Evidence suggest germ cells are remarkably diverse in the type of genome instability that they display and the DNA damage responses they deploy. Additionally, the initial embryonic cell cycles are characterized by a high degree of genome instability that cause congenital disorders and may limit the use of CRISPR-Cas9 for heritable genome editing.

Repair of DNA Breaks by Break-Induced Replication

Double-strand DNA breaks (DSBs) are the most lethal type of DNA damage, making DSB repair critical for cell survival. However, some DSB repair pathways are mutagenic and promote genome rearrangements, leading to genome destabilization. One such pathway is break-induced replication (BIR), which repairs primarily one-ended DSBs, similar to those formed by collapsed replication forks or telomere erosion. BIR is initiated by the invasion of a broken DNA end into a homologous template, synthesizes new DNA within the context of a migrating bubble, and is associated with conservative inheritance of new genetic material. This mode of synthesis is responsible for a high level of genetic instability associated with BIR. Eukaryotic BIR was initially investigated in yeast, but now it is also actively studied in mammalian systems. Additionally, a significant breakthrough has been made regarding the role of microhomology-mediated BIR in the formation of complex genomic rearrangements that underly various human pathologies.

Confined placental mosaicism involving multiple de novo copy number variants associated with fetal growth restriction: A case report

The presence of multiple large (>1 Mb) copy number variants (CNVs) in non‐malignant tissue is rare in human genetics. We present a liveborn male with a birth weight below the first percentile associated with placental mosaicism involving eight 2.4–3.9 Mb de novo duplications. We found that the duplications likely co‐localized to the same cells, were mosaic in the placenta, and impacted maternal and paternal chromosomes. In addition, 27.4 Mb and 240 genes were duplicated in affected cells, including candidate placental genes KISS1 and REN. We ruled out involvement of homologous recombination‐based mechanisms or an altered epigenome in generating the CNVs. This case highlights the diversity of genetic abnormalities in the human placenta and the gaps in our knowledge of how such errors arise.

Chromoanagenesis, the mechanisms of a genomic chaos

Designated under the name of chromoanagenesis, the phenomena of chromothripsis, chromanasynthesis and chromoplexy constitute new types of complex rearrangements, including many genomic alterations localized on a few chromosomal regions, and whose discovery over the last decade has changed our perception about the formation of chromosomal abnormalities and their etiology. Although exhibiting specific features, these new catastrophic mechanisms generally occur within a single cell cycle and their emergence is closely linked to genomic instability. Various non-exclusive exogenous or cellular mechanisms capable of generating chromoanagenesis have been evoked. However, recent experimental data shed light on 2 major processes, which following a defect in the mitotic segregation of chromosomes, can generate a cascade of cellular events leading to chromoanagenesis. These mechanisms are the formation of micronuclei integrating isolated chromosomal material, and the occurrence of chromatin bridges around chromosomal material resulting from telomeric fusions. In both cases, the cellular and molecular mechanisms of fragmentation, repair and transmission of damaged chromosomal material are consistent with the features of chromoanagenesis-related complex chromosomal rearrangements. In this review, we introduce each type of chromoanagenesis, and describe the experimental models that have allowed to validate the existence of chromoanagenesis events and to better understand their cellular mechanisms of formation and transmission, as well as their impact on the stability and the plasticity of the genome.

Efficient detection of copy-number variations using exome data: Batch- and sex-based analyses

Many algorithms to detect copy number variations (CNVs) using exome sequencing (ES) data have been reported and evaluated on their sensitivity and specificity, reproducibility, and precision. However, operational optimization of such algorithms for a better performance has not been fully addressed. ES of 1199 samples including 763 patients with different disease profiles was performed. ES data were analyzed to detect CNVs by both the eXome Hidden Markov Model (XHMM) and modified Nord's method. To efficiently detect rare CNVs, we aimed to decrease sequencing biases by analyzing, at the same time, the data of all unrelated samples sequenced in the same flow cell as a batch, and to eliminate sex effects of X-linked CNVs by analyzing female and male sequences separately. We also applied several filtering steps for more efficient CNV selection. The average number of CNVs detected in one sample was <5. This optimization together with targeted CNV analysis by Nord's method identified pathogenic/likely pathogenic CNVs in 34 patients (4.5%, 34/763). In particular, among 142 patients with epilepsy, the current protocol detected clinically relevant CNVs in 19 (13.4%) patients, whereas the previous protocol identified them in only 14 (9.9%) patients. Thus, this batch-based XHMM analysis efficiently selected rare pathogenic CNVs in genetic diseases.

Characterization of systemic genomic instability in budding yeast

Conventional models of genome evolution are centered around the principle that mutations form independently of each other and build up slowly over time. We characterized the occurrence of bursts of genome-wide loss-of-heterozygosity (LOH) in Saccharomyces cerevisiae, providing support for an additional nonindependent and faster mode of mutation accumulation. We initially characterized a yeast clone isolated for carrying an LOH event at a specific chromosome site, and surprisingly found that it also carried multiple unselected rearrangements elsewhere in its genome. Whole-genome analysis of over 100 additional clones selected for carrying primary LOH tracts revealed that they too contained unselected structural alterations more often than control clones obtained without any selection. We also measured the rates of coincident LOH at two different chromosomes and found that double LOH formed at rates 14- to 150-fold higher than expected if the two underlying single LOH events occurred independently of each other. These results were consistent across different strain backgrounds and in mutants incapable of entering meiosis. Our results indicate that a subset of mitotic cells within a population can experience discrete episodes of systemic genomic instability, when the entire genome becomes vulnerable and multiple chromosomal alterations can form over a narrow time window. They are reminiscent of early reports from the classic yeast genetics literature, as well as recent studies in humans, both in cancer and genomic disorder contexts. The experimental model we describe provides a system to further dissect the fundamental biological processes responsible for punctuated bursts of structural genomic variation.

Postzygotic Somatic Mutations in the Human Brain Expand the Threshold-Liability Model of Schizophrenia

The search for what causes schizophrenia has been onerous. This research has included extensive assessment of a variety of genetic and environmental factors using ever emerging high-resolution technologies and traditional understanding of the biology of the brain. These efforts have identified a large number of schizophrenia-associated genes, some of which are altered by mutational and epi-mutational mechanisms in a threshold liability model of schizophrenia development. The results, however, have limited predictability and the actual cause of the disease remains unknown. This current state asks for conceptualizing the problem differently in light of novel insights into the nature of mutations, the biology of the brain and the fine precision and resolution of emerging technologies. There is mounting evidence that mutations acquired during postzygotic development are more common than germline mutations. Also, the postzygotic somatic mutations including epimutations (PZMs), which often lead to somatic mosaicism, are relatively common in the mammalian brain in comparison to most other tissues and PZMs are more common in patients with neurodevelopmental mental disorders, including schizophrenia. Further, previously inaccessible, detection of PZMs is becoming feasible with the advent of novel technologies that include single-cell genomics and epigenomics and the use of exquisite experimental designs including use of monozygotic twins discordant for the disease. These developments allow us to propose a working hypothesis and expand the threshold liability model of schizophrenia that already encompasses familial genetic, epigenetic and environmental factors to include somatic de novo PZMs. Further, we offer a test for this expanded model using currently available genome sequences and methylome data on monozygotic twins discordant for schizophrenia (MZD) and their parents. The results of this analysis argue that PZMs play a significant role in the development of schizophrenia and explain extensive heterogeneity seen across patients. It also offers the potential to convincingly link PZMs to both nervous system health and disease, an area that has remained challenging to study and relatively under explored.

Clinical genomics and contextualizing genome variation in the diagnostic laboratory

INTRODUCTION

The human genome contains the instructions for the development and biological homeostasis of the human organism and the genetic transmission of traits. Genome variation in human populations is the basis of evolution; individual or personal genomes vary tremendously, making each of us truly unique.

AREAS COVERED

Assaying this individual variation using genomic technologies has many applications in clinical medicine, from elucidating the biology of disease to designing strategies to ameliorate perturbations from homeostasis. Detecting pathogenic rare variation in a genome may provide a molecular diagnosis that can be informative for patient management and family healthcare.

EXPERT OPINION

Despite the increasing clinical use of unbiased genomic testing, including chromosome microarray analysis (CMA) with array comparative genomic hybridization (aCGH) or SNP arrays, clinical exome sequencing (cES), and whole-genome sequencing (WGS), to survey genome-wide for molecular aberrations, clinical acumen paired with an understanding of the limitations of each testing type will be needed to achieve molecular diagnoses. Potential opportunities for improving case solved rates, functionally annotating the majority of genes in the human genome, and further understanding genetic contributions to disease will empower clinical genomics and the precision medicine initiative.

Prenatal Diagnosis of a De Novo Partial Trisomy 6q and Partial Monosomy 18p Associated with Cephalocele: A Case Report

A 28-year-old woman underwent amniocentesis at 18 weeks’ gestation upon detection of increased fetal nuchal fold and parietal cephalocele on the second trimester ultrasound examination. Prenatal microarray showed a de novo unbalanced translocation resulting in a gain in 6q and loss in 18p. A female infant was delivered at 38 weeks’ gestation. At birth, cephalocele and webbed neck were noted as major dysmorphic features. The case presented here shows how a combination of different genetic studies is used to accurately elucidate a chromosomal anomaly in a prenatal setting.

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.

Established and Novel Mechanisms Leading to de novo Genomic Rearrangements in the Human Germline

During gametogenesis, the human genome can acquire various de novo rearrangements. Most constitutional genomic rearrangements are created through 1 of the 4 well-known mechanisms, i.e., nonallelic homologous recombination, erroneous repair after double-strand DNA breaks, replication errors, and retrotransposition. However, recent studies have identified 2 types of extremely complex rearrangements that cannot be simply explained by these mechanisms. The first type consists of chaotic structural changes in 1 or a few chromosomes that result from “chromoanagenesis (an umbrella term that covers chromothripsis, chromoanasynthesis, and chromoplexy).” The other type is large independent rearrangements in multiple chromosomes indicative of “transient multifocal genomic crisis.” Germline chromoanagenesis (chromothripsis) likely occurs predominantly during spermatogenesis or postzygotic embryogenesis, while multifocal genomic crisis appears to be limited to a specific time window during oogenesis and early embryogenesis or during spermatogenesis. This review article introduces the current understanding of the molecular basis of de novo rearrangements in the germline.

A de novo 2.2 Mb recurrent 17q23.1q23.2 deletion unmasks novel putative regulatory non-coding SNVs associated with lethal lung hypoplasia and pulmonary hypertension: a case report

BACKGROUND

Application of whole genome sequencing (WGS) enables identification of non-coding variants that play a phenotype-modifying role and are undetectable by exome sequencing. Recently, non-coding regulatory single nucleotide variants (SNVs) have been reported in patients with lethal lung developmental disorders (LLDDs) or congenital scoliosis with recurrent copy-number variant (CNV) deletions at 17q23.1q23.2 or 16p11.2, respectively.

CASE PRESENTATION

Here, we report a deceased newborn with pulmonary hypertension and pulmonary interstitial emphysema with features suggestive of pulmonary hypoplasia, resulting in respiratory failure and neonatal death soon after birth. Using the array comparative genomic hybridization and WGS, two heterozygous recurrent CNV deletions: ~ 2.2 Mb on 17q23.1q23.2, involving TBX4, and ~ 600 kb on 16p11.2, involving TBX6, that both arose de novo on maternal chromosomes were identified. In the predicted lung-specific enhancer upstream to TBX4, we have detected seven novel putative regulatory non-coding SNVs that were absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies.

CONCLUSIONS

Our findings further support a recently reported model of complex compound inheritance of LLDD in which both non-coding and coding heterozygous TBX4 variants contribute to the lung phenotype. In addition, this is the first report of a patient with combined de novo heterozygous recurrent 17q23.1q23.2 and 16p11.2 CNV deletions.

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.

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.

Distinct patterns of complex rearrangements and a mutational signature of microhomeology are frequently observed in PLP1 copy number gain structural variants

BACKGROUND

We investigated the features of the genomic rearrangements in a cohort of 50 male individuals with proteolipid protein 1 (PLP1) copy number gain events who were ascertained with Pelizaeus-Merzbacher disease (PMD; MIM: 312080). We then compared our new data to previous structural variant mutagenesis studies involving the Xq22 region of the human genome. The aggregate data from 159 sequenced join-points (discontinuous sequences in the reference genome that are joined during the rearrangement process) were studied. Analysis of these data from 150 individuals enabled the spectrum and relative distribution of the underlying genomic mutational signatures to be delineated.

METHODS

Genomic rearrangements in PMD individuals with PLP1 copy number gain events were investigated by high-density customized array or clinical chromosomal microarray analysis and breakpoint junction sequence analysis.

RESULTS

High-density customized array showed that the majority of cases (33/50; ~ 66%) present with single duplications, although complex genomic rearrangements (CGRs) are also frequent (17/50; ~ 34%). Breakpoint mapping to nucleotide resolution revealed further previously unknown structural and sequence complexities, even in single duplications. Meta-analysis of all studied rearrangements that occur at the PLP1 locus showed that single duplications were found in ~ 54% of individuals and that, among all CGR cases, triplication flanked by duplications is the most frequent CGR array CGH pattern observed. Importantly, in ~ 32% of join-points, there is evidence for a mutational signature of microhomeology (highly similar yet imperfect sequence matches).

CONCLUSIONS

These data reveal a high frequency of CGRs at the PLP1 locus and support the assertion that replication-based mechanisms are prominent contributors to the formation of CGRs at Xq22. We propose that microhomeology can facilitate template switching, by stabilizing strand annealing of the primer using W-C base complementarity, and is a mutational signature for replicative repair.

Xq22 deletions and correlation with distinct neurological disease traits in females: Further evidence for a contiguous gene syndrome

Xq22 deletions that encompass PLP1 (Xq22-PLP1-DEL) are notable for variable expressivity of neurological disease traits in females ranging from a mild late-onset form of spastic paraplegia type 2 (MIM# 312920), sometimes associated with skewed X-inactivation, to an early-onset neurological disease trait (EONDT) of severe developmental delay, intellectual disability, and behavioral abnormalities. Size and gene content of Xq22-PLP1-DEL vary and were proposed as potential molecular etiologies underlying variable expressivity in carrier females where two smallest regions of overlap (SROs) were suggested to influence disease. We ascertained a cohort of eight unrelated patients harboring Xq22-PLP1-DEL and performed high-density array comparative genomic hybridization and breakpoint-junction sequencing. Molecular characterization of Xq22-PLP1-DEL from 17 cases (eight herein and nine published) revealed an overrepresentation of breakpoints that reside within repeats (11/17, ~65%) and the clustering of ~47% of proximal breakpoints in a genomic instability hotspot with characteristic non-B DNA density. These findings implicate a potential role for genomic architecture in stimulating the formation of Xq22-PLP1-DEL. The correlation of Xq22-PLP1-DEL gene content with neurological disease trait in female cases enabled refinement of the associated SROs to a single genomic interval containing six genes. Our data support the hypothesis that genes contiguous to PLP1 contribute to EONDT.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

In vitro fertilization does not increase the incidence of de novo copy number alterations in fetal and placental lineages

Although chromosomal instability (CIN) is a common phenomenon in cleavage-stage embryogenesis following in vitro fertilization (IVF)^1-3, its rate in naturally conceived human embryos is unknown. CIN leads to mosaic embryos that contain a combination of genetically normal and abnormal cells, and is significantly higher in in vitro-produced preimplantation embryos as compared to in vivo-conceived preimplantation embryos⁴. Even though embryos with CIN-derived complex aneuploidies may arrest between the cleavage and blastocyst stages of embryogenesis^5,6, a high number of embryos containing abnormal cells can pass this strong selection barrier^7,8. However, neither the prevalence nor extent of CIN during prenatal development and at birth, following IVF treatment, is well understood. Here we profiled the genomic landscape of fetal and placental tissues postpartum from both IVF and naturally conceived children, to investigate the prevalence and persistence of large genetic aberrations that probably arose from IVF-related CIN. We demonstrate that CIN is not preserved at later stages of prenatal development, and that de novo numerical aberrations or large structural DNA imbalances occur at similar rates in IVF and naturally conceived live-born neonates. Our findings affirm that human IVF treatment has no detrimental effect on the chromosomal constitution of fetal and placental lineages.

Exome-wide assessment of the functional impact and pathogenicity of multinucleotide mutations

Approximately 2% of de novo single-nucleotide variants (SNVs) appear as part of clustered mutations that create multinucleotide variants (MNVs). MNVs are an important source of genomic variability as they are more likely to alter an encoded protein than a SNV, which has important implications in disease as well as evolution. Previous studies of MNVs have focused on their mutational origins and have not systematically evaluated their functional impact and contribution to disease. We identified 69,940 MNVs and 91 de novo MNVs in 6688 exome-sequenced parent–offspring trios from the Deciphering Developmental Disorders Study comprising families with severe developmental disorders. We replicated the previously described MNV mutational signatures associated with DNA polymerase zeta, an error-prone translesion polymerase, and the APOBEC family of DNA deaminases. We estimate the simultaneous MNV germline mutation rate to be 1.78 × 10⁻¹⁰ mutations per base pair per generation. We found that most MNVs within a single codon create a missense change that could not have been created by a SNV. MNV-induced missense changes were, on average, more physicochemically divergent, were more depleted in highly constrained genes (pLI ≥ 0.9), and were under stronger purifying selection compared with SNV-induced missense changes. We found that de novo MNVs were significantly enriched in genes previously associated with developmental disorders in affected children. This shows that MNVs can be more damaging than SNVs even when both induce missense changes, and are an important variant type to consider in relation to human disease.

Paternal exposure to benzo(a)pyrene induces genome-wide mutations in mouse offspring

Understanding the effects of environmental exposures on germline mutation rates has been a decades-long pursuit in genetics. We used next-generation sequencing and comparative genomic hybridization arrays to investigate genome-wide mutations in the offspring of male mice exposed to benzo(a)pyrene (BaP), a common environmental pollutant. We demonstrate that offspring developing from sperm exposed during the mitotic or post-mitotic phases of spermatogenesis have significantly more de novo single nucleotide variants (1.8-fold; P < 0.01) than controls. Both phases of spermatogenesis are susceptible to the induction of heritable mutations, although mutations arising from post-fertilization events are more common after post-mitotic exposure. In addition, the mutation spectra in sperm and offspring of BaP-exposed males are consistent. Finally, we report a significant increase in transmitted copy number duplications (P = 0.001) in BaP-exposed sires. Our study demonstrates that germ cell mutagen exposures induce genome-wide mutations in the offspring that may be associated with adverse health outcomes.

Transient multifocal genomic crisis creating chromothriptic and non-chromothriptic rearrangements in prezygotic testicular germ cells

Background

The co-occurrence of multiple de novo copy number variations (CNVs) is a rare phenomenon in the human genome. Recently, an “organismal CNV mutator phenotype” has been reported to result in transient genomic instability introducing multiple de novo CNVs in primary oocytes and early-stage zygotes. These findings opened a new area of human genome research.

Methods

We performed genome-wide copy number analysis for ~ 2100 individuals with various congenital defects. Furthermore, extensive molecular analyses, including synthetic long-read whole-genome sequencing and haplotype-phasing, were carried out for an individual with multiple de novo CNVs.

Results

A boy was found to have de novo rearrangements on five chromosomes. The rearrangements comprised simple duplication and inversion as well as chaotic changes, all of which affected paternally derived chromosomes. Postzygotic genomic instability was ruled out. The duplicated regions on 6q and 13q contained both diallelic and triallelic loci, indicating that the genomic rearrangements were initially created during premeiotic mitosis and subsequently modified by physiological cross-over during meiosis I. Breakpoints of the rearrangements were indicative of non-homologous end joining, replication-based errors, and/or chromothripsis. The mutagenic event was independent of specific local DNA motifs or de novo point mutations, but may be driven by spermatogenesis-specific factors.

Conclusions

These results indicate that during spermatogenesis, a transient multifocal genomic crisis can introduce several chromothriptic and non-chromothriptic changes into the genome. These findings broaden the concept of the “organismal CNV mutator phenotype”. This study provides insights into mechanisms for altering the global chromosomal architecture of human embryos.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0526-3) contains supplementary material, which is available to authorized users.

Preserving Genome Integrity During the Early Embryonic DNA Replication Cycles

During the very early stages of embryonic development chromosome replication occurs under rather challenging conditions, including a very short cell cycle, absence of transcription, a relaxed DNA damage response and, in certain animal species, a highly contracted S-phase. This raises the puzzling question of how the genome can be faithfully replicated in such a peculiar metabolic context. Recent studies have provided new insights into this issue, and unveiled that embryos are prone to accumulate genetic and genomic alterations, most likely due to restricted cellular functions, in particular reduced DNA synthesis quality control. These findings may explain the low rate of successful development in mammals and the occurrence of diseases, such as abnormal developmental features and cancer. In this review, we will discuss recent findings in this field and put forward perspectives to further study this fascinating question.

Ontogenetic and Pathogenetic Views on Somatic Chromosomal Mosaicism

Intercellular karyotypic variability has been a focus of genetic research for more than 50 years. It has been repeatedly shown that chromosome heterogeneity manifesting as chromosomal mosaicism is associated with a variety of human diseases. Due to the ability of changing dynamically throughout the ontogeny, chromosomal mosaicism may mediate genome/chromosome instability and intercellular diversity in health and disease in a bottleneck fashion. However, the ubiquity of negligibly small populations of cells with abnormal karyotypes results in difficulties of the interpretation and detection, which may be nonetheless solved by post-genomic cytogenomic technologies. In the post-genomic era, it has become possible to uncover molecular and cellular pathways to genome/chromosome instability (chromosomal mosaicism or heterogeneity) using advanced whole-genome scanning technologies and bioinformatic tools. Furthermore, the opportunities to determine the effect of chromosomal abnormalities on the cellular phenotype seem to be useful for uncovering the intrinsic consequences of chromosomal mosaicism. Accordingly, a post-genomic review of chromosomal mosaicism in the ontogenetic and pathogenetic contexts appears to be required. Here, we review chromosomal mosaicism in its widest sense and discuss further directions of cyto(post)genomic research dedicated to chromosomal heterogeneity.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

Background

Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.

Methods

Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.

Results

We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.

Conclusions

These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0633-y) contains supplementary material, which is available to authorized users.

Chromosome segregation errors generate a diverse spectrum of simple and complex genomic rearrangements

Cancer genomes are frequently characterized by numerical and structural chromosomal abnormalities. Here we integrated a centromere-specific inactivation approach with selection for a conditionally essential gene, a strategy termed CEN-SELECT, to systematically interrogate the structural landscape of mis-segregated chromosomes. We show that single-chromosome mis-segregation into a micronucleus can directly trigger a broad spectrum of genomic rearrangement types. Cytogenetic profiling revealed that mis-segregated chromosomes exhibit 120-fold-higher susceptibility to developing seven major categories of structural aberrations, including translocations, insertions, deletions, and complex reassembly through chromothripsis coupled to classical non-homologous end joining. Whole-genome sequencing of clonally propagated rearrangements identified random patterns of clustered breakpoints with copy-number alterations resulting in interspersed gene deletions and extrachromosomal DNA amplification events. We conclude that individual chromosome segregation errors during mitotic cell division are sufficient to drive extensive structural variations that recapitulate genomic features commonly associated with human disease.

Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Origins of Intraindividual Genetic Variation in Human Fetuses

BACKGROUND

Intraindividual copy number variation (CNV) origin is largely unknown. They might be due to aging and/or common genome instability at the preimplantation stage while contribution of preimplantation in human intraindividual CNVs occurrence is unknown. To address this question, we investigated mosaicism and its origin in the fetuses of natural conception.

METHODS

We studied normal fetuses following therapeutic abortion due to maternal indications. We analyzed the genome of 22 tissues of each fetus by array comparative genomic hybridization for intraindividual CNVs. Each tissue was studied in 2 microarray experiments; the reciprocal aberrations larger than 40 Kb, identified by comparing tissues of each fetus, were subsequently validated using quantitative polymerase chain reaction.

RESULTS

Through intraindividual comparison, frequency of reciprocal events varied from 2 to 9. According to the distribution pattern of the frequent CNV in derivatives of different germ layers, we found that its origin is early development including preimplantation, whereas CNVs with low frequency have occurred in later stages. Shared CNVs in both fetuses were belonged to thymus and related to the functional role of genes located in these CNVs.

CONCLUSIONS

The origin of some of fetal CNVs is preimplantation stage. Each organ might inherit CNVs with an unpredictable pattern due to the extensive cell mixing/migration in embryonic development.

Predicting human genes susceptible to genomic instability associated with Alu/Alu-mediated rearrangements

Alu elements, the short interspersed element numbering more than 1 million copies per human genome, can mediate the formation of copy number variants (CNVs) between substrate pairs. These Alu/Alu-mediated rearrangements (AAMRs) can result in pathogenic variants that cause diseases. To investigate the impact of AAMR on gene variation and human health, we first characterized Alus that are involved in mediating CNVs (CNV-Alus) and observed that these Alus tend to be evolutionarily younger. We then computationally generated, with the assistance of a supercomputer, a test data set consisting of 78 million Alu pairs and predicted ∼18% of them are potentially susceptible to AAMR. We further determined the relative risk of AAMR in 12,074 OMIM genes using the count of predicted CNV-Alu pairs and experimentally validated the predictions with 89 samples selected by correlating predicted hotspots with a database of CNVs identified by clinical chromosomal microarrays (CMAs) on the genomes of approximately 54,000 subjects. We fine-mapped 47 duplications, 40 deletions, and two complex rearrangements and examined a total of 52 breakpoint junctions of simple CNVs. Overall, 94% of the candidate breakpoints were at least partially Alu mediated. We successfully predicted all (100%) of Alu pairs that mediated deletions (n = 21) and achieved an 87% positive predictive value overall when including AAMR-generated deletions and duplications. We provided a tool, AluAluCNVpredictor, for assessing AAMR hotspots and their role in human disease. These results demonstrate the utility of our predictive model and provide insights into the genomic features and molecular mechanisms underlying AAMR.

Marker chromosome genomic structure and temporal origin implicate a chromoanasynthesis event in a family with pleiotropic psychiatric phenotypes

Small supernumerary marker chromosomes (sSMC) are chromosomal fragments difficult to characterize genomically. Here, we detail a proband with schizoaffective disorder and a mother with bipolar disorder with psychotic features who present with a marker chromosome that segregates with disease. We explored the architecture of this marker and investigated its temporal origin. Array comparative genomic hybridization (aCGH) analysis revealed three duplications and three triplications that spanned the short arm of chromosome 9, suggestive of a chromoanasynthesis-like event. Segregation of marker genotypes, phased using sSMC mosaicism in the mother, provided evidence that it was generated during a germline-level event in the proband's maternal grandmother. Whole-genome sequencing (WGS) was performed to resolve the structure and junctions of the chromosomal fragments, revealing further complexities. While structural variations have been previously associated with neuropsychiatric disorders and marker chromosomes, here we detail the precise architecture, human life-cycle genesis, and propose a DNA replicative/repair mechanism underlying formation.

Integrating CNVs into meta-QTL identified GBP4 as positional candidate for adult cattle stature

Copy number variation (CNV) of DNA sequences, functionally significant but yet fully ascertained, is believed to confer considerable increments in unexplained heritability of quantitative traits. Identification of phenotype-associated CNVs (paCNVs) therefore is a pressing need in CNV studies to speed up their exploitation in cattle breeding programs. Here, we provided a new avenue to achieve this goal that is to project the published CNV data onto meta-quantitative trait loci (meta-QTL) map which connects causal genes with phenotypes. Any CNVs overlapping meta-QTL therefore will be potential paCNVs. This study reported potential paCNVs in Bos taurus autosome 3 (BTA3). Notably, overview indexes and CNVs both highlighted a narrower region (BTA3 54,500,000-55,000,000 bp, named BTA3_INQTL_6) within one constructed meta-QTL. Then, we ascertained guanylate-binding protein 4 (GBP4) among the nine positional candidate genes was significantly associated with adult cattle stature, including body weight (BW, P < 0.05) and withers height (WHT, P < 0.05), fitting GBP4 CNV either with three levels or with six levels in the model. Although higher copy number downregulated the mRNA levels of GBP2 (P < 0.05) and GBP4 (P < 0.05) in 1-Mb window (54.0-55.0 Mb) in muscle and adipose, additional analyses will be needed to clarify the causality behind the ascertained association.