Mechanisms underlying structural variant formation in genomic disorders

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.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes

BACKGROUND

African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages.

RESULTS

We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage.

CONCLUSION

This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Integrated Analysis of Germline and Tumor DNA Identifies New Candidate Genes Involved in Familial Colorectal Cancer

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson’s two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

Mitotic CDK Promotes Replisome Disassembly, Fork Breakage, and Complex DNA Rearrangements

DNA replication errors generate complex chromosomal rearrangements and thereby contribute to tumorigenesis and other human diseases. One mechanism that triggers these errors is mitotic entry before the completion of DNA replication. To address how mitosis might affect DNA replication, we used Xenopus egg extracts. When mitotic CDK (Cyclin B1-CDK1) is used to drive interphase egg extracts into a mitotic state, the replicative CMG (CDC45/MCM2-7/GINS) helicase undergoes ubiquitylation on its MCM7 subunit, dependent on the E3 ubiquitin ligase TRAIP. Whether replisomes have stalled or undergone termination, CMG ubiquitylation is followed by its extraction from chromatin by the CDC48/p97 ATPase. TRAIP-dependent CMG unloading during mitosis is also seen in C. elegans early embryos. At stalled forks, CMG removal results in fork breakage and end joining events involving deletions and templated insertions. Our results identify a mitotic pathway of global replisome disassembly that can trigger replication fork collapse and DNA rearrangements.

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.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

Variational Inference for Coupled Hidden Markov Models Applied to the Joint Detection of Copy Number Variations

Hidden Markov models provide a natural statistical framework for the detection of the copy number variations (CNV) in genomics. In this context, we define a hidden Markov process that underlies all individuals jointly in order to detect and to classify genomics regions in different states (typically, deletion, normal or amplification). Structural variations from different individuals may be dependent. It is the case in agronomy where varietal selection program exists and species share a common phylogenetic past. We propose to take into account these dependencies inthe HMM model. When dealing with a large number of series, maximum likelihood inference (performed classically using the EM algorithm) becomes intractable. We thus propose an approximate inference algorithm based on a variational approach (VEM), implemented in the CHMM R package. A simulation study is performed to assess the performance of the proposed method and an application to the detection of structural variations in plant genomes is presented.

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.

Factors affecting template switch recombination associated with restarted DNA replication

Homologous recombination helps ensure the timely completion of genome duplication by restarting collapsed replication forks. However, this beneficial function is not without risk as replication restarted by homologous recombination is prone to template switching (TS) that can generate deleterious genome rearrangements associated with diseases such as cancer. Previously we established an assay for studying TS in Schizosaccharomyces pombe (Nguyen et al., 2015). Here, we show that TS is detected up to 75 kb downstream of a collapsed replication fork and can be triggered by head-on collision between the restarted fork and RNA Polymerase III transcription. The Pif1 DNA helicase, Pfh1, promotes efficient restart and also suppresses TS. A further three conserved helicases (Fbh1, Rqh1 and Srs2) strongly suppress TS, but there is no change in TS frequency in cells lacking Fml1 or Mus81. We discuss how these factors likely influence TS.

Heterochromatin suppresses gross chromosomal rearrangements at centromeres by repressing Tfs1/TFIIS-dependent transcription

Heterochromatin, characterized by histone H3 lysine 9 (H3K9) methylation, assembles on repetitive regions including centromeres. Although centromeric heterochromatin is important for correct segregation of chromosomes, its exact role in maintaining centromere integrity remains elusive. Here, we found in fission yeast that heterochromatin suppresses gross chromosomal rearrangements (GCRs) at centromeres. Mutations in Clr4/Suv39 methyltransferase increased the formation of isochromosomes, whose breakpoints were located in centromere repeats. H3K9A and H3K9R mutations also increased GCRs, suggesting that Clr4 suppresses centromeric GCRs via H3K9 methylation. HP1 homologs Swi6 and Chp2 and the RNAi component Chp1 were the chromodomain proteins essential for full suppression of GCRs. Remarkably, mutations in RNA polymerase II (RNAPII) or Tfs1/TFIIS, the transcription factor that facilitates restart of RNAPII after backtracking, specifically bypassed the requirement of Clr4 for suppressing GCRs. These results demonstrate that heterochromatin suppresses GCRs by repressing Tfs1-dependent transcription of centromere repeats.

Accelerated rates of large-scale mutations in the presence of copper and nickel

Mutation rate variation has been under intense investigation for decades. Despite these efforts, little is known about the extent to which environmental stressors accelerate mutation rates and influence the genetic load of populations. Moreover, most studies on stressors have focused on unicellular organisms and point mutations rather than large-scale deletions and duplications (copy number variations [CNVs]). We estimated mutation rates in Daphnia pulex exposed to low levels of environmental stressors as well as the effect of selection on de novo mutations. We conducted a mutation accumulation (MA) experiment in which selection was minimized, coupled with an experiment in which a population was propagated under competitive conditions in a benign environment. After an average of 103 generations of MA propagation, we sequenced 60 genomes and found significantly accelerated rates of deletions and duplications in MA lines exposed to ecologically relevant concentrations of metals. Whereas control lines had gene deletion and duplication rates comparable to other multicellular eukaryotes (1.8 × 10-6 per gene per generation), the presence of nickel and copper increased these rates fourfold. The realized mutation rate under selection was reduced to 0.4× that of control MA lines, providing evidence that CNVs contribute to mutational load. Our CNV breakpoint analysis revealed that nonhomologous recombination associated with regions of DNA fragility is the primary source of CNVs, plausibly linking metal-induced DNA strand breaks with higher CNV rates. Our findings suggest that environmental stress, in particular multiple stressors, can have profound effects on large-scale mutation rates and mutational load of multicellular organisms.

Pelizaeus-Merzbacher Disease: Molecular and Cellular Pathologies and Associated Phenotypes

Pelizaeus-Merzbacher disease (PMD) represents a group of disorders known as hypomyelinating leukodystrophies, which are characterized by abnormal development and maintenance of myelin in the central nervous system. PMD is caused by different types of mutations in the proteolipid protein 1 (PLP1) gene, which encodes a major myelin membrane lipoprotein. These mutations in the PLP1 gene result in distinct cellular and molecular pathologies and a spectrum of clinical phenotypes. In this chapter, I discuss the historical aspects and current understanding of the mechanisms underlying how different PLP1 mutations disrupt the normal process of myelination and result in PMD and other disorders.

DNA Damage and Repair in Human Reproductive Cells

The fundamental underlying paradigm of sexual reproduction is the production of male and female gametes of sufficient genetic difference and quality that, following syngamy, they result in embryos with genomic potential to allow for future adaptive change and the ability to respond to selective pressure. The fusion of dissimilar gametes resulting in the formation of a normal and viable embryo is known as anisogamy, and is concomitant with precise structural, physiological, and molecular control of gamete function for species survival. However, along the reproductive life cycle of all organisms, both male and female gametes can be exposed to an array of “stressors” that may adversely affect the composition and biological integrity of their proteins, lipids and nucleic acids, that may consequently compromise their capacity to produce normal embryos. The aim of this review is to highlight gamete genome organization, differences in the chronology of gamete production between the male and female, the inherent DNA protective mechanisms in these reproductive cells, the aetiology of DNA damage in germ cells, and the remarkable DNA repair mechanisms, pre- and post-syngamy, that function to maintain genome integrity.

The copy number variation landscape of congenital anomalies of the kidney and urinary tract

Congenital anomalies of the kidney and urinary tract (CAKUT) are a major cause of pediatric kidney failure. We performed a genome-wide analysis of copy number variants (CNVs) in 2,824 cases and 21,498 controls. Affected individuals carried a significant burden of rare exonic (i.e. affecting coding regions) CNVs and were enriched for known genomic disorders (GD). Kidney anomaly (KA) cases were most enriched for exonic CNVs, encompassing GD-CNVs and novel deletions; obstructive uropathy (OU) had a lower CNV burden and an intermediate prevalence of GD-CNVs; vesicoureteral reflux (VUR) had the fewest GD-CNVs but was enriched for novel exonic CNVs, particularly duplications. Six loci (1q21, 4p16.1-p16.3, 16p11.2, 16p13.11, 17q12, and 22q11.2) accounted for 65% of patients with GD-CNVs. Deletions at 17q12, 4p16.1-p16.3, and 22q11.2 were specific for KA; the 16p11.2 locus showed extensive pleiotropy. Using a multidisciplinary approach, we identified TBX6 as a driver for the CAKUT subphenotypes in the 16p11.2 microdeletion syndrome.

Detection and visualization of complex structural variants from long reads

Background

With applications in cancer, drug metabolism, and disease etiology, understanding structural variation in the human genome is critical in advancing the thrusts of individualized medicine. However, structural variants (SVs) remain challenging to detect with high sensitivity using short read sequencing technologies. This problem is exacerbated when considering complex SVs comprised of multiple overlapping or nested rearrangements. Longer reads, such as those from Pacific Biosciences platforms, often span multiple breakpoints of such events, and thus provide a way to unravel small-scale complexities in SVs with higher confidence.

Results

We present CORGi (COmplex Rearrangement detection with Graph-search), a method for the detection and visualization of complex local genomic rearrangements. This method leverages the ability of long reads to span multiple breakpoints to untangle SVs that appear very complicated with respect to a reference genome. We validated our approach against both simulated long reads, and real data from two long read sequencing technologies. We demonstrate the ability of our method to identify breakpoints inserted in synthetic data with high accuracy, and the ability to detect and plot SVs from NA12878 germline, achieving 88.4% concordance between the two sets of sequence data. The patterns of complexity we find in many NA12878 SVs match known mechanisms associated with DNA replication and structural variant formation, and highlight the ability of our method to automatically label complex SVs with an intuitive combination of adjacent or overlapping reference transformations.

Conclusions

CORGi is a method for interrogating genomic regions suspected to contain local rearrangements using long reads. Using pairwise alignments and graph search CORGi produces labels and visualizations for local SVs of arbitrary complexity.

Common fragile site instability in normal cells: Lessons and perspectives

Mechanisms and events related to common fragile site (CFS) instability are well known in cancer cells. Here, we argue that normal cells remain an important experimental model to address questions related to CFS instability in the absence of alterations in cell cycle and DNA damage repair pathways, which are common features acquired in cancer. Furthermore, a major gap of knowledge concerns the stability of CFSs during gametogenesis. CFS instability in meiotic or postmeiotic stages of the germ cell line could generate chromosome deletions or large rearrangements. This in turn can lead to the functional loss of the several CFS-associated genes with tumor suppressor function. Our hypothesis is that such mutations can potentially result in genetic predisposition to develop cancer. Indirect evidence for CFS instability in human germ cells has been provided by genomic investigations in family pedigrees associated with genetic disease. The issue of CFS instability in the germ cell line should represent one of the future efforts, and may take advantage of the existence of sequence and functional conservation of CFSs between rodents and humans.

Keeping ribosomal DNA intact: a repeating challenge

More than half of the human genome consists of repetitive sequences, with the ribosomal DNA (rDNA) representing two of the largest repeats. Repetitive rDNA sequences may form a threat to genomic integrity and cellular homeostasis due to the challenging aspects of their transcription, replication, and repair. Predisposition to cancer, premature aging, and neurological impairment in ataxia-telangiectasia and Bloom syndrome, for instance, coincide with increased cellular rDNA repeat instability. However, the mechanisms by which rDNA instability contributes to these hereditary syndromes and tumorigenesis remain unknown. Here, we review how cells govern rDNA stability and how rDNA break repair influences expansion and contraction of repeat length, a process likely associated with human disease. Recent advancements in CRISPR-based genome engineering may help to explain how cells keep their rDNA intact in the near future.

Endogenous DNA Double-Strand Breaks during DNA Transactions: Emerging Insights and Methods for Genome-Wide Profiling

DNA double-strand breaks (DSBs) jeopardize genome integrity and can-when repaired unfaithfully-give rise to structural rearrangements associated with cancer. Exogenous agents such as ionizing radiation or chemotherapy can invoke DSBs, but a vast amount of breakage arises during vital endogenous DNA transactions, such as replication and transcription. Additionally, chromatin looping involved in 3D genome organization and gene regulation is increasingly recognized as a possible contributor to DSB events. In this review, we first discuss insights into the mechanisms of endogenous DSB formation, showcasing the trade-off between essential DNA transactions and the intrinsic challenges that these processes impose on genomic integrity. In the second part, we highlight emerging methods for genome-wide profiling of DSBs, and discuss future directions of research that will help advance our understanding of genome-wide DSB formation and repair.

Mitotic antipairing of homologous and sex chromosomes via spatial restriction of two haploid sets

Mitotic recombination must be prevented to maintain genetic stability across daughter cells, but the underlying mechanism remains elusive. We report that mammalian cells impede homologous chromosome pairing during mitosis by keeping the two haploid chromosome sets apart, positioning them to either side of a meridional plane defined by the centrosomes. Chromosome oscillation analysis revealed collective genome behavior of noninteracting chromosome sets. Male translocation mice with a maternal-derived supernumerary chromosome display the tracer chromosome exclusively to the haploid set containing the X chromosome. This haploid set-based antipairing motif is shared by multiple cell types, is doubled in tetraploid cells, and is lost in carcinoma cells. The data provide a model of nuclear polarity through the antipairing of homologous chromosomes during mitosis.

Pairing homologous chromosomes is required for recombination. However, in nonmeiotic stages it can lead to detrimental consequences, such as allelic misregulation and genome instability, and is rare in human somatic cells. How mitotic recombination is prevented—and how genetic stability is maintained across daughter cells—is a fundamental, unanswered question. Here, we report that both human and mouse cells impede homologous chromosome pairing by keeping two haploid chromosome sets apart throughout mitosis. Four-dimensional analysis of chromosomes during cell division revealed that a haploid chromosome set resides on either side of a meridional plane, crossing two centrosomes. Simultaneous tracking of chromosome oscillation and the spindle axis, using fluorescent CENP-A and centrin1, respectively, demonstrates collective genome behavior/segregation of two haploid sets throughout mitosis. Using 3D chromosome imaging of a translocation mouse with a supernumerary chromosome, we found that this maternally derived chromosome is positioned by parental origin. These data, taken together, support the identity of haploid sets by parental origin. This haploid set-based antipairing motif is shared by multiple cell types, doubles in tetraploid cells, and is lost in a carcinoma cell line. The data support a mechanism of nuclear polarity that sequesters two haploid sets along a subcellular axis. This topological segregation of haploid sets revisits an old model/paradigm and provides implications for maintaining mitotic fidelity.

Complex structural variants in Mendelian disorders: identification and breakpoint resolution using short- and long-read genome sequencing

Background

Studies have shown that complex structural variants (cxSVs) contribute to human genomic variation and can cause Mendelian disease. We aimed to identify cxSVs relevant to Mendelian disease using short-read whole-genome sequencing (WGS), resolve the precise variant configuration and investigate possible mechanisms of cxSV formation.

Methods

We performed short-read WGS and analysis of breakpoint junctions to identify cxSVs in a cohort of 1324 undiagnosed rare disease patients. Long-read WGS and gene expression analysis were used to resolve one case.

Results

We identified three pathogenic cxSVs: a de novo duplication-inversion-inversion-deletion affecting ARID1B, a de novo deletion-inversion-duplication affecting HNRNPU and a homozygous deletion-inversion-deletion affecting CEP78. Additionally, a de novo duplication-inversion-duplication overlapping CDKL5 was resolved by long-read WGS demonstrating the presence of both a disrupted and an intact copy of CDKL5 on the same allele, and gene expression analysis showed both parental alleles of CDKL5 were expressed. Breakpoint analysis in all the cxSVs revealed both microhomology and longer repetitive elements.

Conclusions

Our results corroborate that cxSVs cause Mendelian disease, and we recommend their consideration during clinical investigations. We show that resolution of breakpoints can be critical to interpret pathogenicity and present evidence of replication-based mechanisms in cxSV formation.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0606-6) contains supplementary material, which is available to authorized users.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1's proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102-104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18-21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

Small supernumerary marker chromosomes: A legacy of trisomy rescue?

We studied by a whole genomic approach and trios genotyping, 12 de novo, nonrecurrent small supernumerary marker chromosomes (sSMC), detected as mosaics during pre- or postnatal diagnosis and associated with increased maternal age. Four sSMCs contained pericentromeric portions only, whereas eight had additional non-contiguous portions of the same chromosome, assembled together in a disordered fashion by repair-based mechanisms in a chromothriptic event. Maternal hetero/isodisomy was detected with a paternal origin of the sSMC in some cases, whereas in others two maternal alleles in the sSMC region and biparental haplotypes of the homologs were detected. In other cases, the homologs were biparental while the sSMC had the same haplotype of the maternally inherited chromosome. These findings strongly suggest that most sSMCs are the result of a multiple-step mechanism, initiated by maternal meiotic nondisjunction followed by postzygotic anaphase lagging of the supernumerary chromosome and its subsequent chromothripsis.

Breaking point: the genesis and impact of structural variation in tumours

Somatic structural variants undoubtedly play important roles in driving tumourigenesis. This is evident despite the substantial technical challenges that remain in accurately detecting structural variants and their breakpoints in tumours and in spite of our incomplete understanding of the impact of structural variants on cellular function. Developments in these areas of research contribute to the ongoing discovery of structural variation with a clear impact on the evolution of the tumour and on the clinical importance to the patient. Recent large whole genome sequencing studies have reinforced our impression of each tumour as a unique combination of mutations but paradoxically have also discovered similar genome-wide patterns of single-nucleotide and structural variation between tumours. Statistical methods have been developed to deconvolute mutation patterns, or signatures, that recur across samples, providing information about the mutagens and repair processes that may be active in a given tumour. These signatures can guide treatment by, for example, highlighting vulnerabilities in a particular tumour to a particular chemotherapy. Thus, although the complete reconstruction of the full evolutionary trajectory of a tumour genome remains currently out of reach, valuable data are already emerging to improve the treatment of cancer.

Ten things you should know about transposable elements

Transposable elements (TEs) are major components of eukaryotic genomes. However, the extent of their impact on genome evolution, function, and disease remain a matter of intense interrogation. The rise of genomics and large-scale functional assays has shed new light on the multi-faceted activities of TEs and implies that they should no longer be marginalized. Here, we introduce the fundamental properties of TEs and their complex interactions with their cellular environment, which are crucial to understanding their impact and manifold consequences for organismal biology. While we draw examples primarily from mammalian systems, the core concepts outlined here are relevant to a broad range of organisms.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease

Purpose

RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease.

Methods

Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members.

Results

Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry.

Conclusion

This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.

Copy number variations and founder effect underlying complete IL-10Rβ deficiency in Portuguese kindreds

Mutations in interleukin-10 receptor (IL-10R) genes are one cause of very early-onset inflammatory bowel disease with perianal lesions, which can be cured by hematopoietic stem cell transplantation. Using a functional test, which assesses responsiveness of peripheral monocytes to IL-10, we identified three unrelated Portuguese patients carrying two novel IL-10RB mutations. In the three patients, sequencing of genomic DNA identified the same large deletion of exon 3 which precluded protein expression. This mutation was homozygous in two patients born from consanguineous families and heterozygous in the third patient born from unrelated parents. Microsatellite analysis of the IL10RB genomic region revealed a common haplotype in the three Portuguese families pointing to a founder deletion inherited from a common ancestor 400 years ago. In the third patient, surface expression of IL-10R was normal but signaling in response to IL-10 was impaired. Complementary DNA sequencing and next-generation sequencing of IL10RB locus with custom-made probes revealed a ≈ 6 Kb duplication encompassing the exon 6 which leads to a frameshift mutation and a loss of the TYK2-interacting Box 2 motif. Altogether, we describe two novel copy number variations in IL10RB, one with founder effect and one preserving cell surface expression but abolishing signaling.

The dark side of centromeres: types, causes and consequences of structural abnormalities implicating centromeric DNA

Centromeres are the chromosomal domains required to ensure faithful transmission of the genome during cell division. They have a central role in preventing aneuploidy, by orchestrating the assembly of several components required for chromosome separation. However, centromeres also adopt a complex structure that makes them susceptible to being sites of chromosome rearrangements. Therefore, preservation of centromere integrity is a difficult, but important task for the cell. In this review, we discuss how centromeres could potentially be a source of genome instability and how centromere aberrations and rearrangements are linked with human diseases such as cancer.

Deletion of exon 4 in LAMA2 is the most frequent mutation in Chinese patients with laminin α2-related muscular dystrophy

Although recessive mutations in LAMA2 are already known to cause laminin α2-related muscular dystrophy, a rare neuromuscular disorder, large deletions or duplications within this gene are not well-characterized. In this study, we applied next-generation sequencing-based copy number variation profiling in 114 individuals clinically diagnosed with laminin α2-related muscular dystrophy, including 96 who harboured LAMA2 mutations and 34 who harboured intragenic rearrangements. In total, we detected 18 distinct LAMA2 copy number variations that have been reported only among Chinese, 10 of which are novel. The frequency of CNVs in the cohort was 19.3%. Deletion of exon 4 was detected in 10 alleles of eight patients, accounting for 27% of all copy number variations. These patients are Han Chinese and were found to have the same haplotype and sequence at the breakpoint junction, suggesting that exon 4 deletion is a founder mutation in Chinese Han and a mutation hotspot. Moreover, the data highlight our approach, a modified next-generation sequencing assay, as a robust and sensitive tool to detect LAMA2 variants; the assay identifies 85.7% of breakpoint junctions directly alongside sequence information. The method can be applied to clinical samples to determine causal variants underlying various Mendelian disorders.

Xq26.3 Duplication in a Boy With Motor Delay and Low Muscle Tone Refines the X-Linked Acrogigantism Genetic Locus

We describe a 4-year-old boy with developmental delay who was found to carry by clinical grade (CG) molecular cytogenetics (MCs) a chromosome Xq26 microduplication. The report prompted a referral of the patient for possible X-linked acrogigantism (X-LAG), a well-defined condition (MIM300942) due to chromosomal microduplication of a nearby region. The patient was evaluated clinically and investigated for endocrine abnormalities related to X-LAG and not only did he not have acrogigantism, but his growth parameters and other hormones were all normal. We then performed high definition MCs and the duplication copy number variant (CNV) was confirmed to precisely map outside the X-LAG critical region and definitely did not harbor the X-LAG candidate gene, GPR101. The patient's phenotype resembled that of other patients with Xq26 CNVs. The case is instructive for the need for high definition MCs when CG MCs' results are inconsistent with the patient's phenotype. It is also useful for further supporting the contention that GPR101 is the gene responsible for X-LAG.

Early adolescent Rai1 reactivation reverses transcriptional and social interaction deficits in a mouse model of Smith-Magenis syndrome

Losing one copy of the RAI1 gene causes Smith–Magenis syndrome (SMS), a neurodevelopmental disorder. Using a newly generated SMS mouse model, this study demonstrates that restoring the Rai1 gene dose in an early postnatal window could repair gene expression and social interaction deficits in this SMS model. The SMS mouse model also showed a reduced density of dendritic spines, anatomical correlates of excitatory synapses, in the prefrontal cortex. Artificial activation of prefrontal cortex neurons partially alleviated the behavioral deficits. These findings suggest that, similar to Rett syndrome, SMS is caused by disruption of a chromatin-modifying gene with reversible developmental phenotypes, highlighting the potential treatment windows in childhood or adolescence.

Haploinsufficiency of Retinoic Acid Induced 1 (RAI1) causes Smith–Magenis syndrome (SMS), a syndromic autism spectrum disorder associated with craniofacial abnormalities, intellectual disability, and behavioral problems. There is currently no cure for SMS. Here, we generated a genetic mouse model to determine the reversibility of SMS-like neurobehavioral phenotypes in Rai1 heterozygous mice. We show that normalizing the Rai1 level 3–4 wk after birth corrected the expression of genes related to neural developmental pathways and fully reversed a social interaction deficit caused by Rai1 haploinsufficiency. In contrast, Rai1 reactivation 7–8 wk after birth was not beneficial. We also demonstrated that the correct Rai1 dose is required in both excitatory and inhibitory neurons for proper social interactions. Finally, we found that Rai1 heterozygous mice exhibited a reduction of dendritic spines in the medial prefrontal cortex (mPFC) and that optogenetic activation of mPFC neurons in adults improved the social interaction deficit of Rai1 heterozygous mice. Together, these results suggest the existence of a postnatal temporal window during which restoring Rai1 can improve the transcriptional and social behavioral deficits in a mouse model of SMS. It is possible that circuit-level interventions would be beneficial beyond this critical window.

Characterization of copy number variants in a large multibreed population of beef and dairy cattle using high-density single nucleotide polymorphism genotype data

Copy number variants (CNVs) are a form of genomic variation that changes the structure of the genome through deletion or duplication of stretches of DNA. The objective of the present study was to characterize CNVs in a large multibreed population of beef and dairy bulls. The CNVs were called on the autosomes of 5,551 cattle from 22 different beef and dairy breeds, using 2 freely available software suites, QuantiSNP and PennCNV. All CNVs were classified into either deletions or duplications. The median concordance between PennCNV and QuantiSNP, per animal, was 18.5% for deletions and 0% for duplications. The low concordance rate between PennCNV and QuantiSNP indicated that neither algorithm, by itself, could identify all CNVs in the population. In total, PennCNV and QuantiSNP collectively identified 747,129 deletions and 432,523 duplications; 80.2% of all duplications and 69.1% of all deletions were present only once in the population. Only 0.154% of all CNVs identified were present in more than 50 animals in the population. The distribution of the percentage of the autosomes that were composed of deletions, per animal, was positively skewed, as was the distribution for the percentage of the autosomes that were composed of duplications, per animal. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of deletions were 0.019%, 0.037%, and 0.201%, respectively. The first quartile, median, and third quartile of the distribution of the percentage of the autosomes that were composed of duplications were 0.013%, 0.028%, and 0.076%, respectively. The distributions of the number of deletions and duplications per animal were both positively skewed. The interquartile range for the number of deletions per animal in the population was between 16 and 117, whereas for duplications it was between 8 and 23. Per animal, there tended to be twice as many deletions as duplications. The distribution of the length of deletions was positively skewed, as was the distribution of the length of duplications. The interquartile range for the length of deletions in the population was between 25 and 101 kb, and for duplications the interquartile range was between 46 and 235 kb. Per animal, duplications tended to be twice as long as deletions. This study provides a description of the characteristics and distribution of CNVs in a large multibreed population of beef and dairy cattle.

Chromosomal microarray analysis on uncultured chorionic villus sampling can be complicated by confined placental mosaicism for aneuploidy and microdeletions

OBJECTIVE

This study aims to establish the incidence and implications of confined placental mosaicism (CPM) in the context of prenatal chromosomal microarray analysis (CMA).

METHODS

We retrospectively reviewed prenatal array data on 1382 consecutive chorionic villus sampling (CVS) specimens spanning the past 6 years, focusing on those for which whole CVS biopsy (both cytotrophoblast and mesenchymal cells) was used for CMA and cultured cells (primarily mesenchyme) was also analyzed or amniotic fluid (AF)/newborn blood was used for confirmation, to determine the frequency of mosaic abnormal findings that were the result of CPM.

RESULTS

Out of a total of 1382 consecutive CVS cases, we identified 42 (42/1382 = 3.0%) cases with abnormal array findings suggestive of mosaicism. Among them, 10 cases were unequivocally interpreted as CPM based on a normal AF/newborn blood confirmatory result. In addition, another 10 cases were interpreted as provisional CPM based on normal results on cultured cells. Notably, 40% (8/20) of the cases revealed complex findings, including multiple mosaic aneuploidies, mosaic submicroscopic copy number variation (CNV), and mosaic aneuploidy plus mosaic CNV.

CONCLUSION

Abnormal CMA results from CVS specimens should be interpreted with caution when mosaicism is evident or suspected. Furthermore, confirmatory testing on amniotic fluid, which contains cells derived from the fetus, is recommended in these cases.

Preserving replication fork integrity and competence via the homologous recombination pathway

Flaws in the DNA replication process have emerged as a leading driver of genome instability in human diseases. Alteration to replication fork progression is a defining feature of replication stress and the consequent failure to maintain fork integrity and complete genome duplication within a single round of S-phase compromises genetic integrity. This includes increased mutation rates, small and large scale genomic rearrangement and deleterious consequences for the subsequent mitosis that result in the transmission of additional DNA damage to the daughter cells. Therefore, preserving fork integrity and replication competence is an important aspect of how cells respond to replication stress and avoid genetic change. Homologous recombination is a pivotal pathway in the maintenance of genome integrity in the face of replication stress. Here we review our recent understanding of the mechanisms by which homologous recombination acts to protect, restart and repair replication forks. We discuss the dynamics of these genetically distinct functions and their contribution to faithful mitoticsegregation.

LINE- and Alu-containing genomic instability hotspot at 16q24.1 associated with recurrent and nonrecurrent CNV deletions causative for ACDMPV

Transposable elements modify human genome by inserting into new loci or by mediating homology-, microhomology-, or homeology-driven DNA recombination or repair, resulting in genomic structural variation. Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal neonatal developmental lung disorder caused by point mutations or copy-number variant (CNV) deletions of FOXF1 or its distant tissue-specific enhancer. Eighty-five percent of 45 ACDMPV-causative CNV deletions, of which junctions have been sequenced, had at least one of their two breakpoints located in a retrotransposon, with more than half of them being Alu elements. We describe a novel ∼35 kb-large genomic instability hotspot at 16q24.1, involving two evolutionarily young LINE-1 (L1) elements, L1PA2 and L1PA3, flanking AluY, two AluSx, AluSx1, and AluJr elements. The occurrence of L1s at this location coincided with the branching out of the Homo-Pan-Gorilla clade, and was preceded by the insertion of AluSx, AluSx1, and AluJr. Our data show that, in addition to mediating recurrent CNVs, L1 and Alu retrotransposons can predispose the human genome to formation of variably sized CNVs, both of clinical and evolutionary relevance. Nonetheless, epigenetic or other genomic features of this locus might also contribute to its increased instability.

Human copy number variants are enriched in regions of low mappability

Copy number variants (CNVs) are known to affect a large portion of the human genome and have been implicated in many diseases. Although whole-genome sequencing (WGS) can help identify CNVs, most analytical methods suffer from limited sensitivity and specificity, especially in regions of low mappability. To address this, we use PopSV, a CNV caller that relies on multiple samples to control for technical variation. We demonstrate that our calls are stable across different types of repeat-rich regions and validate the accuracy of our predictions using orthogonal approaches. Applying PopSV to 640 human genomes, we find that low-mappability regions are approximately 5 times more likely to harbor germline CNVs, in stark contrast to the nearly uniform distribution observed for somatic CNVs in 95 cancer genomes. In addition to known enrichments in segmental duplication and near centromeres and telomeres, we also report that CNVs are enriched in specific types of satellite and in some of the most recent families of transposable elements. Finally, using this comprehensive approach, we identify 3455 regions with recurrent CNVs that were missing from existing catalogs. In particular, we identify 347 genes with a novel exonic CNV in low-mappability regions, including 29 genes previously associated with disease.

Patterns and mechanisms of structural variations in human cancer

Next-generation sequencing technology has enabled the comprehensive detection of genomic alterations in human somatic cells, including point mutations, chromosomal rearrangements, and structural variations (SVs). Using sophisticated bioinformatics algorithms, unbiased catalogs of SVs are emerging from thousands of human cancer genomes for the first time. Via careful examination of SV breakpoints at single-nucleotide resolution as well as local DNA copy number changes, diverse patterns of genomic rearrangements are being revealed. These "SV signatures" provide deep insight into the mutational processes that have shaped genome changes in human somatic cells. This review summarizes the characteristics of recently identified complex SVs, including chromothripsis, chromoplexy, microhomology-mediated breakage-induced replication (MMBIR), and others, to provide a holistic snapshot of the current knowledge on genomic rearrangements in somatic cells.

The impact of transposable elements in adaptive evolution

The growing knowledge about the influence of transposable elements (TEs) on (a) long-term genome and transcriptome evolution; (b) genomic, transcriptomic and epigenetic variation within populations; and (c) patterns of somatic genetic differences in individuals continues to spur the interest of evolutionary biologists in the role of TEs in adaptive evolution. As TEs can trigger a broad range of molecular variation in a population with potentially severe fitness and phenotypic consequences for individuals, different mechanisms evolved to keep TE activity in check, allowing for a dynamic interplay between the host, its TEs and the environment in evolution. Here, we review evidence for adaptive phenotypic changes associated with TEs and the basic molecular mechanisms by which the underlying genetic changes arise: (a) domestication, (b) exaptation, (c) host gene regulation, (d) TE-mediated formation of intronless gene copies-so-called retrogenes and (e) overall increased genome plasticity. Furthermore, we review and discuss how the stress-dependent incapacitation of defence mechanisms against the activity of TEs might facilitate adaptive responses to environmental challenges and how such mechanisms might be particularly relevant in species frequently facing novel environments, such as invasive, pathogenic or parasitic species.

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.

Mechanisms of genetic instability caused by (CGG)n repeats in an experimental mammalian system

We describe a new experimental system to study genome instability caused by fragile X (CGG)_n repeats in mammalian cells. It is based on a selectable cassette carrying the HyTK gene under the control of the FMR1 promoter with (CGG)_n repeats in its 5′-UTR, which was integrated into the unique RL5 site in murine erythroid leukemia cells. Carrier-size (CGG)_n repeats dramatically elevate the frequency of the reporter’s inactivation making cells ganciclovir-resistant. These resistant clones have a unique mutational signature: a change in the repeat length concurrent with mutagenesis in the reporter gene. Inactivation of genes implicated in break-induced replication including POLD3, POLD4, RAD52, RAD51 and SMARCAL1, reduced the frequency of ganciclovir-resistant clones to the baseline level that was observed in the absence of (CGG)_n repeats. We propose that replication fork collapse at carrier-size (CGG)_n repeats can trigger break-induced replication, which result in simultaneous repeat length changes and mutagenesis at a distance.

Chromatin loop anchors are associated with genome instability in cancer and recombination hotspots in the germline

Background

Chromatin loops form a basic unit of interphase nuclear organization, with chromatin loop anchor points providing contacts between regulatory regions and promoters. However, the mutational landscape at these anchor points remains under-studied. Here, we describe the unusual patterns of somatic mutations and germline variation associated with loop anchor points and explore the underlying features influencing these patterns.

Results

Analyses of whole genome sequencing datasets reveal that anchor points are strongly depleted for single nucleotide variants (SNVs) in tumours. Despite low SNV rates in their genomic neighbourhood, anchor points emerge as sites of evolutionary innovation, showing enrichment for structural variant (SV) breakpoints and a peak of SNVs at focal CTCF sites within the anchor points. Both CTCF-bound and non-CTCF anchor points harbour an excess of SV breakpoints in multiple tumour types and are prone to double-strand breaks in cell lines. Common fragile sites, which are hotspots for genome instability, also show elevated numbers of intersecting loop anchor points. Recurrently disrupted anchor points are enriched for genes with functions in cell cycle transitions and regions associated with predisposition to cancer. We also discover a novel class of CTCF-bound anchor points which overlap meiotic recombination hotspots and are enriched for the core PRDM9 binding motif, suggesting that the anchor points have been foci for diversity generated during recent human evolution.

Conclusions

We suggest that the unusual chromatin environment at loop anchor points underlies the elevated rates of variation observed, marking them as sites of regulatory importance but also genomic fragility.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1483-4) contains supplementary material, which is available to authorized users.

The 3D genome: From fundamental principles to disease and cancer

In higher eukaryotes, the three-dimensional (3D) organization of the genome is intimately related to numerous key biological functions including gene expression, DNA repair and DNA replication regulations. Alteration of this 3D organization is detrimental to the organism and can give rise to a broad range of diseases such as cancers. Here, we review recent advances in the field. We first describe how the genome is packed in 3D to form chromosome territories, compartments and domains. We also give an overview of the recent techniques that allow to map the genome in 3D up to the kilobase resolution. We then discuss potential mechanisms by which genome misfolding can affect proper gene expression by distal enhancers, and how the 3D genome influences the formation of genomic rearrangements.

The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be 'induced' by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10%) and pathogenic/likely pathogenic CNVs (overlap 99.69%) was significantly higher than the random distribution (p < 1 × 10^-6 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

A Variety of Alu-Mediated Copy Number Variations Can Underlie IL-12Rβ1 Deficiency

PURPOSE

Inborn errors of IFN-γ immunity underlie Mendelian susceptibility to mycobacterial disease (MSMD). Autosomal recessive complete IL-12Rβ1 deficiency is the most frequent genetic etiology of MSMD. Only two of the 84 known mutations are copy number variations (CNVs), identified in two of the 213 IL-12Rβ1-deficient patients and two of the 164 kindreds reported. These two CNVs are large deletions found in the heterozygous or homozygous state. We searched for novel families with IL-12Rβ1 deficiency due to CNVs.

METHODS

We studied six MSMD patients from five unrelated kindreds displaying adverse reactions to BCG vaccination. Three of the patients also presented systemic salmonellosis, two had mucocutaneous candidiasis, and one had disseminated histoplasmosis. We searched for CNVs and other variations by IL12RB1-targeted next-generation sequencing (NGS).

RESULTS

We identified six new IL-12Rβ1-deficient patients with a complete loss of IL-12Rβ1 expression on phytohemagglutinin-activated T cells and/or EBV-transformed B cells. The cells of these patients did not respond to IL-12 and IL-23. Five different CNVs encompassing IL12RB1 (four deletions and one duplication) were identified in these patients by NGS coverage analysis, either in the homozygous state (n = 1) or in trans (n = 4) with a single-nucleotide variation (n = 3) or a small indel (n = 1). Seven of the nine mutations are novel. Interestingly, four of the five CNVs were predicted to be driven by nearby Alu elements, as well as the two previously reported large deletions. The IL12RB1 locus is actually enriched in Alu elements (44.7%), when compared with the rest of the genome (10.5%).

CONCLUSION

The IL12RB1 locus is Alu-enriched and therefore prone to rearrangements at various positions. CNVs should be considered in the genetic diagnosis of IL-12Rβ1 deficiency.

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.

Accurate detection of complex structural variations using single-molecule sequencing

Structural variations are the greatest source of genetic variation, but they remain poorly understood because of technological limitations. Single-molecule long-read sequencing has the potential to dramatically advance the field, although high error rates are a challenge with existing methods. Addressing this need, we introduce open-source methods for long-read alignment (NGMLR; https://github.com/philres/ngmlr ) and structural variant identification (Sniffles; https://github.com/fritzsedlazeck/Sniffles ) that provide unprecedented sensitivity and precision for variant detection, even in repeat-rich regions and for complex nested events that can have substantial effects on human health. In several long-read datasets, including healthy and cancerous human genomes, we discovered thousands of novel variants and categorized systematic errors in short-read approaches. NGMLR and Sniffles can automatically filter false events and operate on low-coverage data, thereby reducing the high costs that have hindered the application of long reads in clinical and research settings.

γ-glutamyl transpeptidase deficiency caused by a large homozygous intragenic deletion in GGT1

γ-Glutamyl transpeptidase deficiency (glutathionuria, OMIM 231950) is a rare disease, with only six patients reported in the literature, although this condition has probably been underdiagnosed due the difficulty to routinely analyze glutathione in clinical samples and to the fact that no genetic defect has been coupled to the disease so far. We report two siblings with mild psychomotor developmental delay and mild neurological symptoms, who presented a markedly increased excretion of glutathione in urine and a very low γ-glutamyl transpeptidase activity in serum. Whole-genome sequencing revealed the presence of a 16.9 kb homozygous deletion in GGT1, one of the genes encoding enzymes with γ-glutamyl transpeptidase activity in the human genome. Close analysis revealed the presence of a 13 bp insertion at the deletion junction. This is the first report of a genetic variant as the cause of glutathionuria. In addition, genetic characterization of the patients' parents and a healthy sibling has provided direct genetic evidence regarding the autosomal recessive nature of this disease.