Mechanisms underlying structural variant formation in genomic disorders

Implications of germline copy-number variations in psychiatric disorders: review of large-scale genetic studies

Copy number variants (CNVs), defined as genome sequences of ≥50 bp that differ in copy number from that in a reference genome, are a common form of structural variation. Germline CNVs account for some of the missing heritability that single nucleotide polymorphisms could not account for. Recent technological advances have had a huge impact on CNV research. Microarray technology enables relatively low-cost, high-throughput, genome-wide measurements, and short-read sequencing technology enables the detection of short CNVs that cannot be detected by microarrays. As a result, large-scale genetic studies have been able to identify a variety of common and rare germline CNVs and their associations with diseases. Rare germline CNVs have been reported to be associated with neuropsychiatric disorders. In this review, we focused on germline CNVs and briefly described their functional characteristics, formation mechanisms, detection methods, related databases, and the latest findings. Finally, we introduced recent large-scale genetic studies to assess associations of CNVs with diseases, especially psychiatric disorders, and discussed the use of CNV-based animal models to investigate the molecular and cellular mechanisms underlying these disorders. The development and implementation of improved detection methods, such as long-read single-molecule sequencing, are expected to provide additional insight into the molecular basis of psychiatric disorders and other complex diseases, thus facilitating basic and clinical research on CNVs.

A high-quality chromosome-level genome assembly reveals genetics for important traits in eggplant

Eggplant (Solanum melongena L.) is an economically important vegetable crop in the Solanaceae family, with extensive diversity among landraces and close relatives. Here, we report a high-quality reference genome for the eggplant inbred line HQ-1315 (S. melongena-HQ) using a combination of Illumina, Nanopore and 10X genomics sequencing technologies and Hi-C technology for genome assembly. The assembled genome has a total size of ~1.17 Gb and 12 chromosomes, with a contig N50 of 5.26 Mb, consisting of 36,582 protein-coding genes. Repetitive sequences comprise 70.09% (811.14 Mb) of the eggplant genome, most of which are long terminal repeat (LTR) retrotransposons (65.80%), followed by long interspersed nuclear elements (LINEs, 1.54%) and DNA transposons (0.85%). The S. melongena-HQ eggplant genome carries a total of 563 accession-specific gene families containing 1009 genes. In total, 73 expanded gene families (892 genes) and 34 contraction gene families (114 genes) were functionally annotated. Comparative analysis of different eggplant genomes identified three types of variations, including single-nucleotide polymorphisms (SNPs), insertions/deletions (indels) and structural variants (SVs). Asymmetric SV accumulation was found in potential regulatory regions of protein-coding genes among the different eggplant genomes. Furthermore, we performed QTL-seq for eggplant fruit length using the S. melongena-HQ reference genome and detected a QTL interval of 71.29-78.26 Mb on chromosome E03. The gene Smechr0301963, which belongs to the SUN gene family, is predicted to be a key candidate gene for eggplant fruit length regulation. Moreover, we anchored a total of 210 linkage markers associated with 71 traits to the eggplant chromosomes and finally obtained 26 QTL hotspots. The eggplant HQ-1315 genome assembly can be accessed at http://eggplant-hq.cn. In conclusion, the eggplant genome presented herein provides a global view of genomic divergence at the whole-genome level and powerful tools for the identification of candidate genes for important traits in eggplant.

Mosaicism in Human Health and Disease

Mosaicism refers to the occurrence of two or more genomes in an individual derived from a single zygote. Germline mosaicism is a mutation that is limited to the gonads and can be transmitted to offspring. Somatic mosaicism is a postzygotic mutation that occurs in the soma, and it may occur at any developmental stage or in adult tissues. Mosaic variation may be classified in six ways: (a) germline or somatic origin, (b) class of DNA mutation (ranging in scale from single base pairs to multiple chromosomes), (c) developmental context, (d) body location(s), (e) functional consequence (including deleterious, neutral, or advantageous), and (f) additional sources of mosaicism, including mitochondrial heteroplasmy, exogenous DNA sources such as vectors, and epigenetic changes such as imprinting and X-chromosome inactivation. Technological advances, including single-cell and other next-generation sequencing, have facilitated improved sensitivity and specificity to detect mosaicism in a variety of biological contexts.

Inverted duplicate DNA sequences increase translocation rates through sequencing nanopores resulting in reduced base calling accuracy

Inverted duplicated DNA sequences are a common feature of structural variants (SVs) and copy number variants (CNVs). Analysis of CNVs containing inverted duplicated DNA sequences using nanopore sequencing identified recurrent aberrant behavior characterized by low confidence, incorrect and missed base calls. Inverted duplicate DNA sequences in both yeast and human samples were observed to have systematic elevation in the electrical current detected at the nanopore, increased translocation rates and decreased sampling rates. The coincidence of inverted duplicated DNA sequences with dramatically reduced sequencing accuracy and an increased translocation rate suggests that secondary DNA structures may interfere with the dynamics of transit of the DNA through the nanopore.

Replication stress at microsatellites causes DNA double-strand breaks and break-induced replication

Short tandemly repeated DNA sequences, termed microsatellites, are abundant in the human genome. These microsatellites exhibit length instability and susceptibility to DNA double-strand breaks (DSBs) due to their tendency to form stable non-B DNA structures. Replication-dependent microsatellite DSBs are linked to genome instability signatures in human developmental diseases and cancers. To probe the causes and consequences of microsatellite DSBs, we designed a dual-fluorescence reporter system to detect DSBs at expanded (CTG/CAG) _n and polypurine/polypyrimidine (Pu/Py) mirror repeat structures alongside the c-myc replication origin integrated at a single ectopic chromosomal site. Restriction cleavage near the (CTG/CAG)₁₀₀ microsatellite leads to homology-directed single-strand annealing between flanking AluY elements and reporter gene deletion that can be detected by flow cytometry. However, in the absence of restriction cleavage, endogenous and exogenous replication stressors induce DSBs at the (CTG/CAG)₁₀₀ and Pu/Py microsatellites. DSBs map to a narrow region at the downstream edge of the (CTG)₁₀₀ lagging-strand template. (CTG/CAG) _n chromosome fragility is repeat length-dependent, whereas instability at the (Pu/Py) microsatellites depends on replication polarity. Strikingly, restriction-generated DSBs and replication-dependent DSBs are not repaired by the same mechanism. Knockdown of DNA damage response proteins increases (Rad18, polymerase (Pol) η, Pol κ) or decreases (Mus81) the sensitivity of the (CTG/CAG)₁₀₀ microsatellites to replication stress. Replication stress and DSBs at the ectopic (CTG/CAG)₁₀₀ microsatellite lead to break-induced replication and high-frequency mutagenesis at a flanking thymidine kinase gene. Our results show that non-B structure-prone microsatellites are susceptible to replication-dependent DSBs that cause genome instability.

Intragenic Deletions of GNAS in Pseudohypoparathyroidism Type 1A Identify a New Region Affecting Methylation of Exon A/B

CONTEXT

Pseudohypoparathyroidism type 1A (PHP1A) and pseudopseudohypoparathyroidism (PPHP) are caused by inactivating mutations in the exons of GNAS that encode the alpha-subunit of the stimulatory G protein (Gsα). In some cases abnormal methylation of exon A/B of GNAS, a hallmark of PHP1B, has been reported.

OBJECTIVE

To identify the underlying genetic basis for PHP1A/PPHP in patients in whom molecular defects were not detected by GNAS sequencing and microarray-based analysis of copy number variations.

METHODS

Whole genome sequencing (WGS) and pyrosequencing of differentially methylated regions (DMRs) of GNAS using genomic deoxyribonucleic acid from affected patients.

RESULTS

We identified 2 novel heterozygous GNAS deletions: a 6.4 kb deletion that includes exon 2 of GNAS in the first proband that was associated with normal methylation (57%) of exon A/B DMR, and a 1438 bp deletion in a second PHP1A patient that encompasses the promoter region and 5' untranslated region of Gsα transcripts, which was inherited from his mother with PPHP. This deletion was associated with reduced methylation (32%) of exon A/B DMR.

CONCLUSIONS

WGS can detect exonic and intronic mutations, including deletions that are too small to be identified by microarray analysis, and therefore is more sensitive than other techniques for molecular analysis of PHP1A/PPHP. One of the deletions we identified led to reduced methylation of exon A/B DMR, further refining a region needed for normal imprinting of this DMR. We propose that deletion of this region can explain why some PHP1A patients have reduced of methylation of the exon A/B DMR.

Dosage-sensitive genes in autism spectrum disorders: From neurobiology to therapy

Autism spectrum disorders (ASDs) are a group of heterogenous neurodevelopmental disorders affecting 1 in 59 children. Syndromic ASDs are commonly associated with chromosomal rearrangements or dosage imbalance involving a single gene. Many of these genes are dosage-sensitive and regulate transcription, protein homeostasis, and synaptic function in the brain. Despite vastly different molecular perturbations, syndromic ASDs share core symptoms including social dysfunction and repetitive behavior. However, each ASD subtype has a unique pathogenic mechanism and combination of comorbidities that require individual attention. We have learned a great deal about how these dosage-sensitive genes control brain development and behaviors from genetically-engineered mice. Here we describe the clinical features of eight monogenic neurodevelopmental disorders caused by dosage imbalance of four genes, as well as recent advances in using genetic mouse models to understand their pathogenic mechanisms and develop intervention strategies. We propose that applying newly developed quantitative molecular and neuroscience technologies will advance our understanding of the unique neurobiology of each disorder and enable the development of personalized therapy.

Detection of simple and complex de novo mutations with multiple reference sequences

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read–derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

XRCC1 promotes replication restart, nascent fork degradation and mutagenic DNA repair in BRCA2-deficient cells

Homologous recombination/end joining (HR/HEJ)-deficient cancers with BRCA mutations utilize alternative DNA double-strand break repair pathways, particularly alternative non-homologous end joining or microhomology-mediated end joining (alt-EJ/MMEJ) during S and G2 cell cycle phases. Depletion of alt-EJ factors, including XRCC1, PARP1 and POLQ, is synthetically lethal with BRCA2 deficiency; yet, XRCC1 roles in HR-deficient cancers and replication stress are enigmatic. Here, we show that after replication stress, XRCC1 forms an active repair complex with POLQ and MRE11 that supports alt-EJ activity in vitro. BRCA2 limits XRCC1 recruitment and repair complex formation to suppress alt-EJ at stalled forks. Without BRCA2 fork protection, XRCC1 enables cells to complete DNA replication at the expense of increased genome instability by promoting MRE11-dependent fork resection and restart. High XRCC1 and MRE11 gene expression negatively impacts Kaplan–Meier survival curves and hazard ratios for HR-deficient breast cancer patients in The Cancer Genome Atlas. The additive effects of depleting both BRCA2 and XRCC1 indicate distinct pathways for replication restart. Our collective data show that XRCC1-mediated processing contributes to replication fork degradation, replication restart and chromosome aberrations in BRCA2-deficient cells, uncovering new roles of XRCC1 and microhomology-mediated repair mechanisms in HR-deficient cancers, with implications for chemotherapeutic strategies targeting POLQ and PARP activities.

Disease-modifying effects of an SCAF4 structural variant in a predominantly SOD1 ALS cohort

OBJECTIVE

To test the hypothesis that rs573116164 will have disease-modifying effects in patients with superoxide dismutase 1 (SOD1) familial amyotrophic lateral sclerosis (fALS), we characterized rs573116164 within a cohort of 190 patients with fALS and 560 healthy age-matched controls to assess the variant for association with various measures of disease.

METHODS

Using a previously described bioinformatics evaluation algorithm, a polymorphic short structural variant associated with SOD1 was identified according to its theoretical effect on gene expression. An 12-18 poly-T repeat (rs573116164) within the 3' untranslated region of serine and arginine rich proteins-related carboxy terminal domain associated factor 4 (SCAF4), a gene that is adjacent to SOD1, was assessed for disease association and influence on survival and age at onset in an fALS cohort using PCR, Sanger sequencing, and capillary separation techniques for allele detection.

RESULTS

In a North American cohort of predominantly SOD1 fALS patients (n =190) and age-matched healthy controls (n = 560), we showed that carriage of an 18T SCAF4 allele was associated with disease within this cohort (odds ratio [OR] 6.6; 95% confidence interval [CI] 3.9-11.2; p = 4.0e-11), but also within non-SOD1 cases (n = 27; OR 5.3; 95% CI 1.9-14.5; p = 0.0014). This finding suggests genetically SOD1-independent effects of SCAF4 on fALS susceptibility. Furthermore, carriage of an 18T allele was associated with a 26-month reduction in survival time (95% CI 6.6-40.8; p = 0.014), but did not affect age at onset of disease.

CONCLUSIONS

The findings in this fALS cohort suggest that rs573116164 could have SOD1-independent and broader relevance in ALS, warranting further investigation in other fALS and sporadic ALS cohorts, as well as studies of functional effects of the 18T variant on gene expression.

Integrated sequencing and array comparative genomic hybridization in familial Parkinson disease

Objective

To determine how single nucleotide variants (SNVs) and copy number variants (CNVs) contribute to molecular diagnosis in familial Parkinson disease (PD), we integrated exome sequencing (ES) and genome-wide array-based comparative genomic hybridization (aCGH) and further probed CNV structure to reveal mutational mechanisms.

Methods

We performed ES on 110 subjects with PD and a positive family history; 99 subjects were also evaluated using genome-wide aCGH. We interrogated ES and aCGH data for pathogenic SNVs and CNVs at Mendelian PD gene loci. We confirmed SNVs via Sanger sequencing and further characterized CNVs with custom-designed high-density aCGH, droplet digital PCR, and breakpoint sequencing.

Results

Using ES, we discovered individuals with known pathogenic SNVs in GBA (p.Glu365Lys, p.Thr408Met, p.Asn409Ser, and p.Leu483Pro) and LRRK2 (p.Arg1441Gly and p.Gly2019Ser). Two subjects were each double heterozygotes for variants in GBA and LRRK2. Based on aCGH, we additionally discovered cases with an SNCA duplication and heterozygous intragenic GBA deletion. Five additional subjects harbored both SNVs (p.Asn52Metfs*29, p.Thr240Met, p.Pro437Leu, and p.Trp453*) and likely disrupting CNVs at the PRKN locus, consistent with compound heterozygosity. In nearly all cases, breakpoint sequencing revealed microhomology, a mutational signature consistent with CNV formation due to DNA replication errors.

Conclusions

Integrated ES and aCGH yielded a genetic diagnosis in 19.3% of our familial PD cohort. Our analyses highlight potential mechanisms for SNCA and PRKN CNV formation, uncover multilocus pathogenic variation, and identify novel SNVs and CNVs for further investigation as potential PD risk alleles.

Optical mapping of the 22q11.2DS region reveals complex repeat structures and preferred locations for non-allelic homologous recombination (NAHR)

The most prevalent microdeletion in humans occurs at 22q11.2, a region rich in chromosome-specific low copy repeats (LCR22s). The structure of this region has defied elucidation due to its size, regional complexity, and haplotype diversity, and is not well represented in the human genome reference. Most individuals with 22q11.2 deletion syndrome (22q11.2DS) carry a de novo hemizygous deletion of ~ 3 Mbp occurring by non-allelic homologous recombination (NAHR) mediated by LCR22s. In this study, optical mapping has been used to elucidate LCR22 structure and variation in 88 individuals in thirty 22q11.2DS families to uncover potential risk factors for germline rearrangements leading to 22q11.2DS offspring. Families were optically mapped to characterize LCR22 structures, NAHR locations, and genomic signatures associated with the deletion. Bioinformatics analyses revealed clear delineations between LCR22 structures in normal and deletion-containing haplotypes. Despite no explicit whole-haplotype predisposing configurations being identified, all NAHR events contain a segmental duplication encompassing FAM230 gene members suggesting preferred recombination sequences. Analysis of deletion breakpoints indicates that preferred recombinations occur between FAM230 and specific segmental duplication orientations within LCR22A and LCR22D, ultimately leading to NAHR. This work represents the most comprehensive analysis of 22q11.2DS NAHR events demonstrating completely contiguous LCR22 structures surrounding and within deletion breakpoints.

Fatal perinatal mitochondrial cardiac failure caused by recurrent de novo duplications in the ATAD3 locus

Background

In about half of all patients with a suspected monogenic disease, genomic investigations fail to identify the diagnosis. A contributing factor is the difficulty with repetitive regions of the genome, such as those generated by segmental duplications. The ATAD3 locus is one such region, in which recessive deletions and dominant duplications have recently been reported to cause lethal perinatal mitochondrial diseases characterized by pontocerebellar hypoplasia or cardiomyopathy, respectively.

Methods

Whole exome, whole genome and long-read DNA sequencing techniques combined with studies of RNA and quantitative proteomics were used to investigate 17 subjects from 16 unrelated families with suspected mitochondrial disease.

Findings

We report six different de novo duplications in the ATAD3 gene locus causing a distinctive presentation including lethal perinatal cardiomyopathy, persistent hyperlactacidemia, and frequently corneal clouding or cataracts and encephalopathy. The recurrent 68 Kb ATAD3 duplications are identifiable from genome and exome sequencing but usually missed by microarrays. The ATAD3 duplications result in the formation of identical chimeric ATAD3A/ATAD3C proteins, altered ATAD3 complexes and a striking reduction in mitochondrial oxidative phosphorylation complex I and its activity in heart tissue.

Conclusions

ATAD3 duplications appear to act in a dominant-negative manner and the de novo inheritance infers a low recurrence risk for families, unlike most pediatric mitochondrial diseases. More than 350 genes underlie mitochondrial diseases. In our experience the ATAD3 locus is now one of the five most common causes of nuclear-encoded pediatric mitochondrial disease but the repetitive nature of the locus means ATAD3 diagnoses may be frequently missed by current genomic strategies.

Funding

Australian NHMRC, US Department of Defense, Japanese AMED and JSPS agencies, Australian Genomics Health Alliance and Australian Mito Foundation.

Genomic Diagnosis for Pediatric Disorders: Revolution and Evolution

Powerful, recent advances in technologies to analyze the genome have had a profound impact on the practice of medical genetics, both in the laboratory and in the clinic. Increasing utilization of genome-wide testing such as chromosomal microarray analysis and exome sequencing have lead a shift toward a "genotype-first" approach. Numerous techniques are now available to diagnose a particular syndrome or phenotype, and while traditional techniques remain efficient tools in certain situations, higher-throughput technologies have become the de facto laboratory tool for diagnosis of most conditions. However, selecting the right assay or technology is challenging, and the wrong choice may lead to prolonged time to diagnosis, or even a missed diagnosis. In this review, we will discuss current core technologies for the diagnosis of classic genetic disorders to shed light on the benefits and disadvantages of these strategies, including diagnostic efficiency, variant interpretation, and secondary findings. Finally, we review upcoming technologies posed to impart further changes in the field of genetic diagnostics as we move toward "genome-first" practice.

Whole Genome Low-Coverage Sequencing Concurrently Detecting Copy Number Variations and Their Underlying Complex Chromosomal Rearrangements by Systematic Breakpoint Mapping in Intellectual Deficiency/Developmental Delay Patients

Simple copy number variations (CNVs) detected by chromosomal microarray (CMA) can result from complex structural changes. Therefore, it is necessary to characterize potential structural changes that cause pathogenic CNVs. We applied whole-genome low-coverage sequencing (WGLCS) to concurrently detect pathogenic CNVs and their associated chromosomal rearrangements in 15 patients. All the patients had an average of 2–3 pathogenic CNVs involving 1–2 chromosomes. WGLCS identified all the 34 pathogenic CNVs found by microarray. By identifying chimeric read pairs, WGLCS mapped 70 breakpoints in these patients, of which 47 were finely mapped at the nucleotide level and confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. In 15 patients, structural rearrangements were defined at molecular level in 13 patients. In 13 patients, WGLCS reveal no additional results in two patients. In another 11 patients, WGLCS revealed new breakpoints or finely mapped the genes disrupted by breakpoints or 1–6 bp microhomology and/or short insertion (4–70 bp) in the breakpoints junctions. However, structural changes in the other two patients still remained unclear after WGLCS was performed. The structural alteration identified in the 13 patients could be divided into the following categories: (1) interstitial inverted duplication with concomitant terminal deletion (inv dup del) (P1,P4,P9,P11); (2) the product of pericentric inversion (P5); (3) ring chromosome (P8); (4) interstitial duplication and/or triplication (P6, P7); and (5) +der(22)t(11;22) (P2,P15); (6) complex structural rearrangements (P3,P12,P14). WGLCS displayed the ability to discover CNVs and define breakpoints and its disrupted genes and its surrounding sequences in one experiment at base-pair-resolution, which help us to learn more about the mechanisms of formation of observed genomic rearrangements, and in which DNA replicative/repair mechanism might contribute to the formation of complex rearrangements in 11 patients. Clear karyotype at molecular level could help provide an accurate evaluation of recurrent risk and guide prenatal diagnosis or reproductive planning.

Improved structural variant interpretation for hereditary cancer susceptibility using long-read sequencing

Purpose

Structural variants (SVs) may be an underestimated cause of hereditary cancer syndromes given the current limitations of short-read next-generation sequencing. Here we investigated the utility of long-read sequencing in resolving germline SVs in cancer susceptibility genes detected through short-read genome sequencing.

Methods

Known or suspected deleterious germline SVs were identified using Illumina genome sequencing across a cohort of 669 advanced cancer patients with paired tumor genome and transcriptome sequencing. Candidate SVs were subsequently assessed by Oxford Nanopore long-read sequencing.

Results

Nanopore sequencing confirmed eight simple pathogenic or likely pathogenic SVs, resolving three additional variants whose impact could not be fully elucidated through short-read sequencing. A recurrent sequencing artifact on chromosome 16p13 and one complex rearrangement on chromosome 5q35 were subsequently classified as likely benign, obviating the need for further clinical assessment. Variant configuration was further resolved in one case with a complex pathogenic rearrangement affecting TSC2.

Conclusion

Our findings demonstrate that long-read sequencing can improve the validation, resolution, and classification of germline SVs. This has important implications for return of results, cascade carrier testing, cancer screening, and prophylactic interventions.

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.

Discovery and quality analysis of a comprehensive set of structural variants and short tandem repeats

Structural variants (SVs) and short tandem repeats (STRs) are important sources of genetic diversity but are not routinely analyzed in genetic studies because they are difficult to accurately identify and genotype. Because SVs and STRs range in size and type, it is necessary to apply multiple algorithms that incorporate different types of evidence from sequencing data and employ complex filtering strategies to discover a comprehensive set of high-quality and reproducible variants. Here we assemble a set of 719 deep whole genome sequencing (WGS) samples (mean 42×) from 477 distinct individuals which we use to discover and genotype a wide spectrum of SV and STR variants using five algorithms. We use 177 unique pairs of genetic replicates to identify factors that affect variant call reproducibility and develop a systematic filtering strategy to create of one of the most complete and well characterized maps of SVs and STRs to date.

SVIM: structural variant identification using mapped long reads

Motivation

Structural variants are defined as genomic variants larger than 50 bp. They have been shown to affect more bases in any given genome than single-nucleotide polymorphisms or small insertions and deletions. Additionally, they have great impact on human phenotype and diversity and have been linked to numerous diseases. Due to their size and association with repeats, they are difficult to detect by shotgun sequencing, especially when based on short reads. Long read, single-molecule sequencing technologies like those offered by Pacific Biosciences or Oxford Nanopore Technologies produce reads with a length of several thousand base pairs. Despite the higher error rate and sequencing cost, long-read sequencing offers many advantages for the detection of structural variants. Yet, available software tools still do not fully exploit the possibilities.

Results

We present SVIM, a tool for the sensitive detection and precise characterization of structural variants from long-read data. SVIM consists of three components for the collection, clustering and combination of structural variant signatures from read alignments. It discriminates five different variant classes including similar types, such as tandem and interspersed duplications and novel element insertions. SVIM is unique in its capability of extracting both the genomic origin and destination of duplications. It compares favorably with existing tools in evaluations on simulated data and real datasets from Pacific Biosciences and Nanopore sequencing machines.

Availability and implementation

The source code and executables of SVIM are available on Github: github.com/eldariont/svim. SVIM has been implemented in Python 3 and published on bioconda and the Python Package Index.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of 5 Gene Signatures in Survival Prediction for Patients with Lung Squamous Cell Carcinoma Based on Integrated Multiomics Data Analysis

Background

Lung squamous cell carcinoma (LSCC) is a frequently diagnosed cancer worldwide, and it has a poor prognosis. The current study is aimed at developing the prediction of LSCC prognosis by integrating multiomics data including transcriptome, copy number variation data, and mutation data analysis, so as to predict patients' survival and discover new therapeutic targets.

Methods

RNASeq, SNP, CNV data, and LSCC patients' clinical follow-up information were downloaded from The Cancer Genome Atlas (TCGA), and the samples were randomly divided into two groups, namely, the training set and the validation set. In the training set, the genes related to prognosis and those with different copy numbers or with different SNPs were integrated to extract features using random forests, and finally, robust biomarkers were screened. In addition, a gene-related prognostic model was established and further verified in the test set and GEO validation set.

Results

We obtained a total of 804 prognostic-related genes and 535 copy amplification genes, 621 copy deletions genes, and 388 significantly mutated genes in genomic variants; noticeably, these genomic variant genes were found closely related to tumor development. A total of 51 candidate genes were obtained by integrating genomic variants and prognostic genes, and 5 characteristic genes (HIST1H2BH, SERPIND1, COL22A1, LCE3C, and ADAMTS17) were screened through random forest feature selection; we found that many of those genes had been reported to be related to LSCC progression. Cox regression analysis was performed to establish 5-gene signature that could serve as an independent prognostic factor for LSCC patients and can stratify risk samples in training set, test set, and external validation set (p < 0.01), and the 5-year survival areas under the curve (AUC) of both training set and validation set were > 0.67.

Conclusion

In the current study, 5 gene signatures were constructed as novel prognostic markers to predict the survival of LSCC patients. The present findings provide new diagnostic and prognostic biomarkers and therapeutic targets for LSCC treatment.

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.

Parental somatic mosaicism for CNV deletions - A need for more sensitive and precise detection methods in clinical diagnostics settings

To further assess the scale and level of parental somatic mosaicism, we queried the CMA database at Baylor Genetics. We selected 50 unrelated families where clinically relevant apparent de novo CNV-deletions were found in the affected probands. Parental blood samples screening using deletion junction-specific PCR revealed four parents with somatic mosaicism. Droplet digital PCR (ddPCR), qPCR, and amplicon-based next-generation sequencing (NGS) were applied to validate these findings. Using ddPCR levels of mosaicism ranged from undetectable to 18.5%. Amplicon-based NGS and qPCR for the father with undetectable mosaicism was able to detect mosaicism at 0.39%. In one mother, ddPCR analysis revealed 15.6%, 10.6%, 8.2%, and undetectable levels of mosaicism in her blood, buccal cells, saliva, and urine samples, respectively. Our data suggest that more sensitive and precise methods, e.g. CNV junction-specific LR-PCR, ddPCR, or qPCR may allow for a more refined assessment of the potential disease recurrence risk for an identified variant.

DNA replication machinery prevents Rad52-dependent single-strand annealing that leads to gross chromosomal rearrangements at centromeres

Homologous recombination between repetitive sequences can lead to gross chromosomal rearrangements (GCRs). At fission yeast centromeres, Rad51-dependent conservative recombination predominantly occurs between inverted repeats, thereby suppressing formation of isochromosomes whose arms are mirror images. However, it is unclear how GCRs occur in the absence of Rad51 and how GCRs are prevented at centromeres. Here, we show that homology-mediated GCRs occur through Rad52-dependent single-strand annealing (SSA). The rad52-R45K mutation, which impairs SSA activity of Rad52 protein, dramatically reduces isochromosome formation in rad51 deletion cells. A ring-like complex Msh2-Msh3 and a structure-specific endonuclease Mus81 function in the Rad52-dependent GCR pathway. Remarkably, mutations in replication fork components, including DNA polymerase α and Swi1/Tof1/Timeless, change the balance between Rad51-dependent recombination and Rad52-dependent SSA at centromeres, increasing Rad52-dependent SSA that forms isochromosomes. Our results uncover a role of DNA replication machinery in the recombination pathway choice that prevents Rad52-dependent GCRs at centromeres.

Duplication with structural modification through extrachromosomal circular and lariat DNA in the human genome

Duplication plays an important role in creating drastic changes in genome evolution. In addition to well-known tandem duplication, duplication can occur such that a duplicated DNA fragment is inserted at another location in the genome. Here, we report several genomic regions in the human genome that could be best explained by two types of insertion-based duplication mechanisms, where a duplicated DNA fragment was modified structurally and then inserted into the genome. In one process, the DNA fragment is turned into an extrachromosomal circular DNA, cut somewhere in the circle, and reintegrated into another location in the genome. And in the other, the DNA fragment forms a “lariat structure” with a “knot”, the strand is swapped at the knot, and is then reintegrated into the genome. Our results suggest that insertion-based duplication may not be a simple process; it may involve a complicated procedures such as structural modification before reintegration. However, the molecular mechanism has yet to be fully understood.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Identification of Structural Variation in Chimpanzees Using Optical Mapping and Nanopore Sequencing

Recent efforts to comprehensively characterize great ape genetic diversity using short-read sequencing and single-nucleotide variants have led to important discoveries related to selection within species, demographic history, and lineage-specific traits. Structural variants (SVs), including deletions and inversions, comprise a larger proportion of genetic differences between and within species, making them an important yet understudied source of trait divergence. Here, we used a combination of long-read and -range sequencing approaches to characterize the structural variant landscape of two additional Pan troglodytes verus individuals, one of whom carries 13% admixture from Pan troglodytes troglodytes. We performed optical mapping of both individuals followed by nanopore sequencing of one individual. Filtering for larger variants (>10 kbp) and combined with genotyping of SVs using short-read data from the Great Ape Genome Project, we identified 425 deletions and 59 inversions, of which 88 and 36, respectively, were novel. Compared with gene expression in humans, we found a significant enrichment of chimpanzee genes with differential expression in lymphoblastoid cell lines and induced pluripotent stem cells, both within deletions and near inversion breakpoints. We examined chromatin-conformation maps from human and chimpanzee using these same cell types and observed alterations in genomic interactions at SV breakpoints. Finally, we focused on 56 genes impacted by SVs in >90% of chimpanzees and absent in humans and gorillas, which may contribute to chimpanzee-specific features. Sequencing a greater set of individuals from diverse subspecies will be critical to establish the complete landscape of genetic variation in chimpanzees.

Whole-genome sequencing reveals complex chromosome rearrangement disrupting NIPBL in infant with Cornelia de Lange syndrome

Clinical laboratory diagnostic evaluation of the genomes of children with suspected genetic disorders, including chromosomal microarray and exome sequencing, cannot detect copy number neutral genomic rearrangements such as inversions, balanced translocations, and complex chromosomal rearrangements (CCRs). We describe an infant with a clinical diagnosis of Cornelia de Lange syndrome (CdLS) in whom chromosome analysis revealed a de novo complex balanced translocation, 46,XY,t(5;7;6)(q11.2;q32;q13)dn. Subsequent molecular characterization by whole-genome sequencing (WGS) identified 23 breakpoints, delineating segments derived from four chromosomes (5;6;7;21) in ancestral or inverted orientation. One of the breakpoints disrupted a known CdLS gene, NIPBL. Further investigation revealed paternal origin of the CCR allele, clustering of the breakpoint junctions, and molecular repair signatures suggestive of a single catastrophic event. Notably, very short DNA segments (25 and 41 bp) were included in the reassembled chromosomes, lending additional support that the DNA repair machinery can detect and repair such segments. Interestingly, there was an independent paternally derived miniscule complex rearrangement, possibly predisposing to subsequent genomic instability. In conclusion, we report a CCR causing a monogenic Mendelian disorder, urging WGS analysis of similar unsolved cases with suspected Mendelian disorders. Breakpoint analysis allowed for identification of the underlying molecular diagnosis and implicated chromoanagenesis in CCR formation.

Structural variant identification and characterization

Structural variant (SV) differences between human genomes can cause germline and mosaic disease as well as inter-individual variation. De-regulation of accurate DNA repair and genomic surveillance mechanisms results in a large number of SVs in cancer. Analysis of the DNA sequences at SV breakpoints can help identify pathways of mutagenesis and regions of the genome that are more susceptible to rearrangement. Large-scale SV analyses have been enabled by high-throughput genome-level sequencing on humans in the past decade. These studies have shed light on the mechanisms and prevalence of complex genomic rearrangements. Recent advancements in both sequencing and other mapping technologies as well as calling algorithms for detection of genomic rearrangements have helped propel SV detection into population-scale studies, and have begun to elucidate previously inaccessible regions of the genome. Here, we discuss the genomic organization of simple and complex SVs, the molecular mechanisms of their formation, and various ways to detect them. We also introduce methods for characterizing SVs and their consequences on human genomes.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Copy-number variation contributes 9% of pathogenicity in the inherited retinal degenerations

PURPOSE

Current sequencing strategies can genetically solve 55-60% of inherited retinal degeneration (IRD) cases, despite recent progress in sequencing. This can partially be attributed to elusive pathogenic variants (PVs) in known IRD genes, including copy-number variations (CNVs), which have been shown as major contributors to unsolved IRD cases.

METHODS

Five hundred IRD patients were analyzed with targeted next-generation sequencing (NGS). The NGS data were used to detect CNVs with ExomeDepth and gCNV and the results were compared with CNV detection with a single-nucleotide polymorphism (SNP) array. Likely causal CNV predictions were validated by quantitative polymerase chain reaction (qPCR).

RESULTS

Likely disease-causing single-nucleotide variants (SNVs) and small indels were found in 55.6% of subjects. PVs in USH2A (11.6%), RPGR (4%), and EYS (4%) were the most common. Likely causal CNVs were found in an additional 8.8% of patients. Of the three CNV detection methods, gCNV showed the highest accuracy. Approximately 30% of unsolved subjects had a single likely PV in a recessive IRD gene.

CONCLUSION

CNV detection using NGS-based algorithms is a reliable method that greatly increases the genetic diagnostic rate of IRDs. Experimentally validating CNVs helps estimate the rate at which IRDs might be solved by a CNV plus a more elusive variant.

Replication-Based Rearrangements Are a Common Mechanism for SNCA Duplication in Parkinson's Disease

BACKGROUND

SNCA multiplication is a genomic cause of familial PD, showing dosage-dependent toxicity. Until now, nonallelic homologous recombination was suggested as the mechanism of SNCA duplication, based on various types of repetitive elements found in the spanning region of the breakpoints. However, the sequence at the breakpoint was analyzed only for 1 case.

OBJECTIVES

We have analyzed the breakpoint sequences of 6 patients with PD who had duplicated SNCA using whole-genome sequencing data to elucidate the mechanism of SNCA duplication.

METHODS

Six patient samples with SNCA duplication underwent whole-genome sequencing. The duplicated regions were defined with nucleotide-resolution breakpoints, which were confirmed by junction polymerase chain reaction and Sanger sequencing. The search for potential non-B DNA-forming sequences and stem-loop structure predictions was conducted.

RESULTS

Duplicated regions ranged from the smallest region of 718.3 kb to the largest one of 4,162 kb. Repetitive elements were found at 8 of the 12 breakpoint sequences on each side of the junction, but none of the pairs shared overt homologies. Five of these six junctions had microhomologies (2-4 bp) at the breakpoint, and a short stretch of sequences was inserted in 3 cases. All except one junction were located within or next to stem-loop structures.

CONCLUSION

Our study has determined that homologous recombination mechanisms involving repetitive elements are not the main cause of the duplication of SNCA. The presence of microhomology at the junctions and their position within stem-loop structures suggest that replication-based rearrangements may be a common mechanism for SNCA amplification. © 2020 International Parkinson and Movement Disorder Society.

Comparison of Arachis monticola with Diploid and Cultivated Tetraploid Genomes Reveals Asymmetric Subgenome Evolution and Improvement of Peanut

Like many important crops, peanut is a polyploid that underwent polyploidization, evolution, and domestication. The wild allotetraploid peanut species Arachis monticola (A. monticola) is an important and unique link from the wild diploid species to cultivated tetraploid species in the Arachis lineage. However, little is known about A. monticola and its role in the evolution and domestication of this important crop. A fully annotated sequence of ≈2.6 Gb A. monticola genome and comparative genomics of the Arachis species is reported. Genomic reconstruction of 17 wild diploids from AA, BB, EE, KK, and CC groups and 30 tetraploids demonstrates a monophyletic origin of A and B subgenomes in allotetraploid peanuts. The wild and cultivated tetraploids undergo asymmetric subgenome evolution, including homoeologous exchanges, homoeolog expression bias, and structural variation (SV), leading to subgenome functional divergence during peanut domestication. Significantly, SV-associated homoeologs tend to show expression bias and correlation with pod size increase from diploids to wild and cultivated tetraploids. Moreover, genomic analysis of disease resistance genes shows the unique alleles present in the wild peanut can be introduced into breeding programs to improve some resistance traits in the cultivated peanuts. These genomic resources are valuable for studying polyploid genome evolution, domestication, and improvement of peanut production and resistance.

The Genomic Complexity of a Large Inversion in Great Tits

Chromosome inversions have clear effects on genome evolution and have been associated with speciation, adaptation, and the evolution of the sex chromosomes. In birds, these inversions may play an important role in hybridization of species and disassortative mating. We identified a large (≈64 Mb) inversion polymorphism in the great tit (Parus major) that encompasses almost 1,000 genes and more than 90% of Chromosome 1A. The inversion occurs at a low frequency in a set of over 2,300 genotyped great tits in the Netherlands with only 5% of the birds being heterozygous for the inversion. In an additional analysis of 29 resequenced birds from across Europe, we found two heterozygotes. The likely inversion breakpoints show considerable genomic complexity, including multiple copy number variable segments. We identified different haplotypes for the inversion, which differ in the degree of recombination in the center of the chromosome. Overall, this remarkable genetic variant is widespread among distinct great tit populations and future studies of the inversion haplotype, including how it affects the fitness of carriers, may help to understand the mechanisms that maintain it.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3′ branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

Genetic Contribution to Congenital Heart Disease (CHD)

Congenital heart defects (CHD) are the most common congenital problems in neonates. The basis for CHD is multifactorial, involving genetic and environmental components. The elucidation of genetic components remains difficult because it is a genetically heterogeneous disease. Currently, the major identified genetic causes include chromosomal abnormalities, large subchromosomal deletions/duplications, and point mutations. However, much more remains to be unraveled. An important insight from the research on the genetics of CHD is that any change at the genetic level that alters the dosage of genes required in any process during heart development results in a developmental defect. The use of conventional gene identification (linkage analysis and direct targeted sequencing) methods followed by the rapid advancements in high-throughput technologies (copy number variant platforms, SNP arrays, and next-generation sequencing) has identified an extensive list of genetic causes. However, the most common presentation of CHD is in the form of sporadic cases. Therefore, it is important to identify their underlying genetic cause. In this review, we revisit the causal genetic factors of CHD and discuss the clinical implications of research in the field.

Visualization tools for human structural variations identified by whole-genome sequencing

Visualizing structural variations (SVs) is a critical step for finding associations between SVs and human traits or diseases. Given that there are many sequencing platforms used for SV identification and given that how best to visualize SVs together with other data, such as read alignments and annotations, depends on research goals, there are dozens of SV visualization tools designed for different research goals and sequencing platforms. Here, we provide a comprehensive survey of over 30 SV visualization tools to help users choose which tools to use. This review targets users who wish to visualize a set of SVs identified from the massively parallel sequencing reads of an individual human genome. We first categorize the ways in which SV visualization tools display SVs into ten major categories, which we denote as view modules. View modules allow readers to understand the features of each SV visualization tool quickly. Next, we introduce the features of individual SV visualization tools from several aspects, including whether SV views are integrated with annotations, whether long-read alignment is displayed, whether underlying data structures are graph-based, the type of SVs shown, whether auditing is possible, whether bird's eye view is available, sequencing platforms, and the number of samples. We hope that this review will serve as a guide for readers on the currently available SV visualization tools and lead to the development of new SV visualization tools in the near future.

A Return to the Origin of the EMGS: Rejuvenating the Quest for Human Germ Cell Mutagens and Determining the Risk to Future Generations

Fifty years ago, the Environmental Mutagen Society (now Environmental Mutagenesis and Genomics Society) was founded with a laser-focus on germ cell mutagenesis and the protection of "our most vital assets"-the sperm and egg genomes. Yet, five decades on, despite the fact that many agents have been demonstrated to induce inherited changes in the offspring of exposed laboratory rodents, there is no consensus on whether human germ cell mutagens exist. We argue that it is time to reevaluate the available data and conclude that we already have evidence for the existence of environmental exposures that impact human germ cells. What is missing are definite data to demonstrate a significant increase in de novo mutations in the offspring of exposed parents. We believe that with over two decades of research advancing knowledge and technologies in genomics, we are at the cusp of generating data to conclusively show that environmental exposures cause heritable de novo changes in the human offspring. We call on the research community to harness our technologies, synergize our efforts, and return to our Founders' original focus. The next 50 years must involve collaborative work between clinicians, epidemiologists, genetic toxicologists, genomics experts and bioinformaticians to precisely define how environmental exposures impact germ cell genomes. It is time for the research and regulatory communities to prepare to interpret the coming outpouring of data and develop a framework for managing, communicating and mitigating the risk of exposure to human germ cell mutagens. Environ. Mol. Mutagen. 61:42-54, 2020. © 2019 Her Majesty the Queen in Right of Canada.

The Driver of Extreme Human-Specific Olduvai Repeat Expansion Remains Highly Active in the Human Genome

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

The Fml1-MHF complex suppresses inter-fork strand annealing in fission yeast

Previously we reported that a process called inter-fork strand annealing (IFSA) causes genomic deletions during the termination of DNA replication when an active replication fork converges on a collapsed fork (Morrow et al., 2017). We also identified the FANCM-related DNA helicase Fml1 as a potential suppressor of IFSA. Here, we confirm that Fml1 does indeed suppress IFSA, and show that this function depends on its catalytic activity and ability to interact with Mhf1-Mhf2 via its C-terminal domain. Finally, a plausible mechanism of IFSA suppression is demonstrated by the finding that Fml1 can catalyse regressed fork restoration in vitro.

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.

LinkedSV for detection of mosaic structural variants from linked-read exome and genome sequencing data

Linked-read sequencing provides long-range information on short-read sequencing data by barcoding reads originating from the same DNA molecule, and can improve detection and breakpoint identification for structural variants (SVs). Here we present LinkedSV for SV detection on linked-read sequencing data. LinkedSV considers barcode overlapping and enriched fragment endpoints as signals to detect large SVs, while it leverages read depth, paired-end signals and local assembly to detect small SVs. Benchmarking studies demonstrate that LinkedSV outperforms existing tools, especially on exome data and on somatic SVs with low variant allele frequencies. We demonstrate clinical cases where LinkedSV identifies disease-causal SVs from linked-read exome sequencing data missed by conventional exome sequencing, and show examples where LinkedSV identifies SVs missed by high-coverage long-read sequencing. In summary, LinkedSV can detect SVs missed by conventional short-read and long-read sequencing approaches, and may resolve negative cases from clinical genome/exome sequencing studies.

Prioritization of genes driving congenital phenotypes of patients with de novo genomic structural variants

BACKGROUND

Genomic structural variants (SVs) can affect many genes and regulatory elements. Therefore, the molecular mechanisms driving the phenotypes of patients carrying de novo SVs are frequently unknown.

METHODS

We applied a combination of systematic experimental and bioinformatic methods to improve the molecular diagnosis of 39 patients with multiple congenital abnormalities and/or intellectual disability harboring apparent de novo SVs, most with an inconclusive diagnosis after regular genetic testing.

RESULTS

In 7 of these cases (18%), whole-genome sequencing analysis revealed disease-relevant complexities of the SVs missed in routine microarray-based analyses. We developed a computational tool to predict the effects on genes directly affected by SVs and on genes indirectly affected likely due to the changes in chromatin organization and impact on regulatory mechanisms. By combining these functional predictions with extensive phenotype information, candidate driver genes were identified in 16/39 (41%) patients. In 8 cases, evidence was found for the involvement of multiple candidate drivers contributing to different parts of the phenotypes. Subsequently, we applied this computational method to two cohorts containing a total of 379 patients with previously detected and classified de novo SVs and identified candidate driver genes in 189 cases (50%), including 40 cases whose SVs were previously not classified as pathogenic. Pathogenic position effects were predicted in 28% of all studied cases with balanced SVs and in 11% of the cases with copy number variants.

CONCLUSIONS

These results demonstrate an integrated computational and experimental approach to predict driver genes based on analyses of WGS data with phenotype association and chromatin organization datasets. These analyses nominate new pathogenic loci and have strong potential to improve the molecular diagnosis of patients with de novo SVs.

Comparative Analyses of Copy-Number Variation in Autism Spectrum Disorder and Schizophrenia Reveal Etiological Overlap and Biological Insights

Compelling evidence in Caucasian populations suggests a role for copy-number variations (CNVs) in autism spectrum disorder (ASD) and schizophrenia (SCZ). We analyzed 1,108 ASD cases, 2,458 SCZ cases, and 2,095 controls in a Japanese population and confirmed an increased burden of rare exonic CNVs in both disorders. Clinically significant (or pathogenic) CNVs, including those at 29 loci common to both disorders, were found in about 8% of ASD and SCZ cases, which was significantly higher than in controls. Phenotypic analysis revealed an association between clinically significant CNVs and intellectual disability. Gene set analysis showed significant overlap of biological pathways in both disorders including oxidative stress response, lipid metabolism/modification, and genomic integrity. Finally, based on bioinformatics analysis, we identified multiple disease-relevant genes in eight well-known ASD/SCZ-associated CNV loci (e.g., 22q11.2, 3q29). Our findings suggest an etiological overlap of ASD and SCZ and provide biological insights into these disorders.

High-throughput long paired-end sequencing of a Fosmid library by PacBio

BACKGROUND

Large insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.

RESULTS

A Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.

CONCLUSIONS

We developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.

Mismatch recognition and subsequent processing have distinct effects on mitotic recombination intermediates and outcomes in yeast

The post-replicative mismatch repair (MMR) system has anti-recombination activity that limits interactions between diverged sequences by recognizing mismatches in strand-exchange intermediates. In contrast to their equivalent roles during replication-error repair, mismatch recognition is more important for anti-recombination than subsequent mismatch processing. To obtain insight into this difference, ectopic substrates with 2% sequence divergence were used to examine mitotic recombination outcome (crossover or noncrossover; CO and NCO, respectively) and to infer molecular intermediates formed during double-strand break repair in Saccharomyces cerevisiae. Experiments were performed in an MMR-proficient strain, a strain with compromised mismatch-recognition activity (msh6Δ) and a strain that retained mismatch-recognition activity but was unable to process mismatches (mlh1Δ). While the loss of either mismatch binding or processing elevated the NCO frequency to a similar extent, CO events increased only when mismatch binding was compromised. The molecular features of NCOs, however, were altered in fundamentally different ways depending on whether mismatch binding or processing was eliminated. These data suggest a model in which mismatch recognition reverses strand-exchange intermediates prior to the initiation of end extension, while subsequent mismatch processing that is linked to end extension specifically destroys NCO intermediates that contain conflicting strand-discrimination signals for mismatch removal.

Structural variant calling: the long and the short of it

Recent research into structural variants (SVs) has established their importance to medicine and molecular biology, elucidating their role in various diseases, regulation of gene expression, ethnic diversity, and large-scale chromosome evolution-giving rise to the differences within populations and among species. Nevertheless, characterizing SVs and determining the optimal approach for a given experimental design remains a computational and scientific challenge. Multiple approaches have emerged to target various SV classes, zygosities, and size ranges. Here, we review these approaches with respect to their ability to infer SVs across the full spectrum of large, complex variations and present computational methods for each approach.

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.

Evaluating nanopore sequencing data processing pipelines for structural variation identification

BACKGROUND

Structural variations (SVs) account for about 1% of the differences among human genomes and play a significant role in phenotypic variation and disease susceptibility. The emerging nanopore sequencing technology can generate long sequence reads and can potentially provide accurate SV identification. However, the tools for aligning long-read data and detecting SVs have not been thoroughly evaluated.

RESULTS

Using four nanopore datasets, including both empirical and simulated reads, we evaluate four alignment tools and three SV detection tools. We also evaluate the impact of sequencing depth on SV detection. Finally, we develop a machine learning approach to integrate call sets from multiple pipelines. Overall SV callers' performance varies depending on the SV types. For an initial data assessment, we recommend using aligner minimap2 in combination with SV caller Sniffles because of their speed and relatively balanced performance. For detailed analysis, we recommend incorporating information from multiple call sets to improve the SV call performance.

CONCLUSIONS

We present a workflow for evaluating aligners and SV callers for nanopore sequencing data and approaches for integrating multiple call sets. Our results indicate that additional optimizations are needed to improve SV detection accuracy and sensitivity, and an integrated call set can provide enhanced performance. The nanopore technology is improving, and the sequencing community is likely to grow accordingly. In turn, better benchmark call sets will be available to more accurately assess the performance of available tools and facilitate further tool development.

Chromatin organization modulates the origin of heritable structural variations in human genome

Structural variations (SVs) in the human genome originate from different mechanisms related to DNA repair, replication errors, and retrotransposition. Our analyses of 26 927 SVs from the 1000 Genomes Project revealed differential distributions and consequences of SVs of different origin, e.g. deletions from non-allelic homologous recombination (NAHR) are more prone to disrupt chromatin organization while processed pseudogenes can create accessible chromatin. Spontaneous double stranded breaks (DSBs) are the best predictor of enrichment of NAHR deletions in open chromatin. This evidence, along with strong physical interaction of NAHR breakpoints belonging to the same deletion suggests that majority of NAHR deletions are non-meiotic i.e. originate from errors during homology directed repair (HDR) of spontaneous DSBs. In turn, the origin of the spontaneous DSBs is associated with transcription factor binding in accessible chromatin revealing the vulnerability of functional, open chromatin. The chromatin itself is enriched with repeats, particularly fixed Alu elements that provide the homology required to maintain stability via HDR. Through co-localization of fixed Alus and NAHR deletions in open chromatin we hypothesize that old Alu expansion had a stabilizing role on the human genome.