Low-copy repeats at the human VIPR2 gene predispose to recurrent and nonrecurrent rearrangements

Inverted triplications formed by iterative template switches generate structural variant diversity at genomic disorder loci

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Tandem and inverted duplications in haemophilia A: Breakpoint characterisation provides insight into possible rearrangement mechanisms

INTRODUCTION

Approximately half of patients with severe haemophilia A are caused by structural variants in the F8 gene. Unlike inversions or deletions directly impairing the integrity of F8, some duplications do not completely disrupt the open reading frame or even retain an intact F8 copy. Currently, only a few duplication breakpoints were precisely characterized, and the corresponding rearrangement mechanisms and clinical outcomes remain to be further investigated.

AIM

Establishing an effective strategy for breakpoint characterization of duplications and revealing their rearrangement mechanisms.

METHODS

AccuCopy is used for the detection of duplications, long-distance PCR for the characterization of tandem duplications, genome walking technique and whole genome sequencing for the characterization of inverted duplications.

RESULTS

Four F8 duplication rearrangements were successfully characterized at the nucleotide level: one tandem duplication (exons 7-11) and three inverted duplications (exons 7-22, exons 2-26, and exons 15-22). Two shared features of inverted duplication were found after carefully analysing our results and breakpoint information in the literature: 1, an inverted fragment was inserted into the original chromosome via two junctions; 2, one junction is mediated by a pair of inverted repetitive elements, while the other consists of two breakpoints with microhomology.

CONCLUSION

Similar breakpoint features motivated us to propose a DNA replication-based model to explain the formation of duplication rearrangements. Based on our model, we further divide the inverted duplications into three basic types: type I with a DEL-NOR/INV-DUP pattern, type II with a DUP-NOR/INV-DUP pattern and type III with a DUP-TRP/INV-DUP pattern.

Overexpression of VIPR2 in mice results in microencephaly with paradoxical increased white matter volume

Large scale studies in populations of European and Han Chinese ancestry found a series of rare gain-of-function microduplications in VIPR2, encoding VPAC2, a receptor that binds vasoactive intestinal peptide and pituitary adenylate cyclase-activating polypeptide with high affinity, that were associated with an up to 13-fold increased risk for schizophrenia. To address how VPAC2 receptor overactivity might affect brain development, we used a well-characterized Nestin-Cre mouse strain and a knock-in approach to overexpress human VPAC2 in the central nervous system. Mice that overexpressed VPAC2 were found to exhibit a significant reduction in brain weight. Magnetic resonance imaging analysis confirmed a decrease in brain size, a specific reduction in the hippocampus grey matter volume and a paradoxical increase in whole-brain white matter volume. Sex-specific changes in behavior such as impaired prepulse inhibition and contextual fear memory were observed in VPAC2 overexpressing mice. The data indicate that the VPAC2 receptor may play a critical role in brain morphogenesis and suggest that overactive VPAC2 signaling during development plays a mechanistic role in some forms of schizophrenia.

Inversions on human chromosomes

Human chromosome inversions are types of balanced structural variations, making them difficult to analyze. Thanks to PEM (paired-end sequencing and mapping), there has been tremendous progress in studying inversions. Inversions play an important role as an evolutionary factor, contributing to the formation of gonosomes, speciation of chimpanzees and humans, and inv17q21.3 or inv8p23.1 exhibit the features of natural selection. Both inversions have been related to pathogenic phenotype by directly affecting a gene structure (e.g., inv5p15.1q14.1), regulating gene expression (e.g., inv7q21.3q35) and by predisposing to other secondary arrangements (e.g., inv7q11.23). A polymorphism of human inversions is documented by the InvFEST database (a database that stores information about clinical predictions, validations, frequency of inversions, etc.), but only a small fraction of these inversions is validated, and a detailed analysis is complicated by the frequent location of breakpoints within regions of repetitive sequences.

Cruciform DNA Structures Act as Legible Templates for Accelerating Homologous Recombination in Transgenic Animals

Inverted repeat (IR) DNA sequences compose cruciform structures. Some genetic disorders are the result of genome inversion or translocation by cruciform DNA structures. The present study examined whether exogenous DNA integration into the chromosomes of transgenic animals was related to cruciform DNA structures. Large imperfect cruciform structures were frequently predicted around predestinated transgene integration sites in host genomes of microinjection-based transgenic (Tg) animals (αLA-LPH Tg goat, Akr1A1eGFP/eGFP Tg mouse, and NFκB-Luc Tg mouse) or CRISPR/Cas9 gene-editing (GE) animals (αLA-AP1 GE mouse). Transgene cassettes were imperfectly matched with their predestinated sequences. According to the analyzed data, we proposed a putative model in which the flexible cruciform DNA structures acted as a legible template for DNA integration into linear DNAs or double-strand break (DSB) alleles. To demonstrate this model, artificial inverted repeat knock-in (KI) reporter plasmids were created to analyze the KI rate using the CRISPR/Cas9 system in NIH3T3 cells. Notably, the KI rate of the 5′ homologous arm inverted repeat donor plasmid (5′IR) with the ROSA gRNA group (31.5%) was significantly higher than the knock-in reporter donor plasmid (KIR) with the ROSA gRNA group (21.3%, p < 0.05). However, the KI rate of the 3′ inverted terminal repeat/inverted repeat donor plasmid (3′ITRIR) group was not different from the KIR group (23.0% vs. 22.0%). These results demonstrated that the legibility of the sequence with the cruciform DNA existing in the transgene promoted homologous recombination (HR) with a higher KI rate. Our findings suggest that flexible cruciform DNAs folded by IR sequences improve the legibility and accelerate DNA 3′-overhang integration into the host genome via homologous recombination machinery.

Mutations in G Protein-Coupled Receptors: Mechanisms, Pathophysiology and Potential Therapeutic Approaches

There are approximately 800 annotated G protein-coupled receptor (GPCR) genes, making these membrane receptors members of the most abundant gene family in the human genome. Besides being involved in manifold physiologic functions and serving as important pharmacotherapeutic targets, mutations in 55 GPCR genes cause about 66 inherited monogenic diseases in humans. Alterations of nine GPCR genes are causatively involved in inherited digenic diseases. In addition to classic gain- and loss-of-function variants, other aspects, such as biased signaling, trans-signaling, ectopic expression, allele variants of GPCRs, pseudogenes, gene fusion, and gene dosage, contribute to the repertoire of GPCR dysfunctions. However, the spectrum of alterations and GPCR involvement is probably much larger because an additional 91 GPCR genes contain homozygous or hemizygous loss-of-function mutations in human individuals with currently unidentified phenotypes. This review highlights the complexity of genomic alteration of GPCR genes as well as their functional consequences and discusses derived therapeutic approaches. SIGNIFICANCE STATEMENT: With the advent of new transgenic and sequencing technologies, the number of monogenic diseases related to G protein-coupled receptor (GPCR) mutants has significantly increased, and our understanding of the functional impact of certain kinds of mutations has substantially improved. Besides the classical gain- and loss-of-function alterations, additional aspects, such as biased signaling, trans-signaling, ectopic expression, allele variants of GPCRs, uniparental disomy, pseudogenes, gene fusion, and gene dosage, need to be elaborated in light of GPCR dysfunctions and possible therapeutic strategies.

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.

Mechanisms underlying structural variant formation in genomic disorders

With the recent burst of technological developments in genomics, and the clinical implementation of genome-wide assays, our understanding of the molecular basis of genomic disorders, specifically the contribution of structural variation to disease burden, is evolving quickly. Ongoing studies have revealed a ubiquitous role for genome architecture in the formation of structural variants at a given locus, both in DNA recombination-based processes and in replication-based processes. These reports showcase the influence of repeat sequences on genomic stability and structural variant complexity and also highlight the tremendous plasticity and dynamic nature of our genome in evolution, health and disease susceptibility.

Methylomic analysis of salivary DNA in childhood ADHD identifies altered DNA methylation in VIPR2

BACKGROUND

Peripheral epigenetic marks hold promise for understanding psychiatric illness and may represent fingerprints of gene-environment interactions. We conducted an initial examination of CpG methylation variation in children with or without attention-deficit/hyperactivity disorder (ADHD).

METHODS

Children age 7-12 were recruited, screened, evaluated and assigned to ADHD or non-ADHD groups by defined research criteria. Two independent age-matched samples were examined, a discovery set (n = 92, all boys, half control, half ADHD) and a confirmation set (n = 20, half ADHD, all boys). 5-methylcytosine levels were quantified in salivary DNA using the Illumina 450 K HumanMethylation array. Genes for which multiple probes were nominally significant and had a beta difference of at least 2% were evaluated for biological relevance and prioritized for confirmation and sequence validation. Gene pathways were explored and described.

RESULTS

Two genes met the criteria for confirmation testing, VIPR2 and MYT1L; both had multiple probes meeting cutoffs and strong biological relevance. Probes on VIPR2 passed FDR correction in the confirmation set and were confirmed through bisulfite sequencing. Enrichment analysis suggested involvement of gene sets or pathways related to inflammatory processes and modulation of monoamine and cholinergic neurotransmission.

CONCLUSIONS

Although it is unknown to what extent CpG methylation seen in peripheral tissue reflect transcriptomic changes in the brain, these initial results indicate that peripheral DNA methylation markers in ADHD may be promising and suggest targeted hypotheses for future study in larger samples.

Human inversions and their functional consequences

Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies.

Disclosing the Hidden Structure and Underlying Mutational Mechanism of a Novel Type of Duplication CNV Responsible for Hereditary Multiple Osteochondromas

The additional mutational complexity associated with copy number variation (CNV) can provide important clues as to the underlying mechanisms of CNV formation. Correct annotation of the additional mutational complexity is, however, a prerequisite for establishing the mutational mechanism. We illustrate this point through the characterization of a novel ∼230 kb EXT1 duplication CNV causing autosomal dominant hereditary multiple osteochondromas. Whole-genome sequencing initially identified the CNV as having a 22-bp insertion at the breakpoint junction and, unprecedentedly, multiple breakpoint-flanking micromutations on both sides of the duplication. Further investigation revealed that this genomic rearrangement had a duplication-inverted triplication-duplication structure, the inverted triplication being a 41-bp sequence synthesized from a nearby template. This permitted the identification of the sequence determinants of both the initiation (an inverted Alu repeat) and termination (a triplex-forming sequence) of break-induced replication and suggested a possible model for the repair of replication-associated double-strand breaks.

Absence of heterozygosity due to template switching during replicative rearrangements

We investigated complex genomic rearrangements (CGRs) consisting of triplication copy-number variants (CNVs) that were accompanied by extended regions of copy-number-neutral absence of heterozygosity (AOH) in subjects with multiple congenital abnormalities. Molecular analyses provided observational evidence that in humans, post-zygotically generated CGRs can lead to regional uniparental disomy (UPD) due to template switches between homologs versus sister chromatids by using microhomology to prime DNA replication-a prediction of the replicative repair model, MMBIR. Our findings suggest that replication-based mechanisms might underlie the formation of diverse types of genomic alterations (CGRs and AOH) implicated in constitutional disorders.

Complex genomic rearrangements at the PLP1 locus include triplication and quadruplication

Inverted repeats (IRs) can facilitate structural variation as crucibles of genomic rearrangement. Complex duplication-inverted triplication-duplication (DUP-TRP/INV-DUP) rearrangements that contain breakpoint junctions within IRs have been recently associated with both MECP2 duplication syndrome (MIM#300260) and Pelizaeus-Merzbacher disease (PMD, MIM#312080). We investigated 17 unrelated PMD subjects with copy number gains at the PLP1 locus including triplication and quadruplication of specific genomic intervals-16/17 were found to have a DUP-TRP/INV-DUP rearrangement product. An IR distal to PLP1 facilitates DUP-TRP/INV-DUP formation as well as an inversion structural variation found frequently amongst normal individuals. We show that a homology-or homeology-driven replicative mechanism of DNA repair can apparently mediate template switches within stretches of microhomology. Moreover, we provide evidence that quadruplication and potentially higher order amplification of a genomic interval can occur in a manner consistent with rolling circle amplification as predicted by the microhomology-mediated break induced replication (MMBIR) model.

Definition and refinement of the 7q36.3 duplication region associated with schizophrenia

Using a very high-resolution oligonucleotide array for copy number variant (CNV) screening of samples comprising schizophrenic patients, we detected a novel CNV within the critical region (NCBI36/hg18, Chr7: 158,630,410-158,719,410) previously shown to be associated with schizophrenia. We investigated the association between the novel CNV identified in the current study and schizophrenia. Three independent samples were used: (1) Screening set, 300 Japanese schizophrenic patients (53.28 ± 14.66 years); (2) Confirmation set, 531 schizophrenic patients (46.03 ± 12.15 years); and (3) 711 healthy controls (47.12 ± 11.03 years). All subjects enrolled in the study were Japanese. Chromosomal position was determined using fluorescence in situ hybridization. We identified a novel duplication within the region associated with schizophrenia identified on 7q36.3 that is adjacent to VIPR2 and is not associated with schizophrenia. In the Japanese population, the 35-kb region that harbors the common, novel CNV should be excluded from the region associated with schizophrenia on 7q36.3.

Search for missing schizophrenia genes will require a new developmental neurogenomic perspective

Even the most powerful experimental designs in search of genetic causes of schizophrenia have not met the desired goal. It is imperative to review the reasons for such an outcome and to formulate novel strategies for the future direction of this research in the new era of individual genomes. Here, we will review aspects of neurodevelopmental hypothesis of schizophrenia in the light of novel genomic and epigenomic insights. Specifically, we will argue for the involvement of de novo mutations and epigenetic modifications during neurodevelopment that may result in schizophrenia. Our conclusion is that the successful elucidation of hereditary mechanisms in neuropsychiatric disorders must begin with attention to discrete endophenotypes; consideration of ontogeny, forethought of genome structure including temporal and spatial patterns of (epi) mutations and the use of judicious techniques that go beyond association studies.

Targeting of vasoactive intestinal peptide receptor 2, VPAC2, a secretin family G-protein coupled receptor, to primary cilia

Primary cilia protrude from the cell surface of many cell types in the human body and function as cellular antennae via ciliary membrane localized receptors. Neurons and glial cells in the brain possess primary cilia, and the malfunction of primary cilia may contribute to neurological deficits present in many cilia-associated disorders. Several rhodopsin family G-protein coupled receptors (GPCRs) are specifically localized to a subset of neuronal primary cilia. However, whether other family GPCRs target to neuronal cilia and whether glial primary cilia harbor any GPCRs are not known. We conducted a screening of GPCRs to determine their ability to target to primary cilia, and identified a secretin family member, Vasoactive Intestinal Receptor 2 (VPAC2), as a novel ciliary GPCR. Here, we show that endogenous VPAC2 targets to primary cilia in various brain regions, including the suprachiasmatic nuclei and the thalamus. Surprisingly, VPAC2 not only localizes to neuronal cilia but also to glial cilia. In addition, we show that VPAC2's C-terminus is both necessary and sufficient for its ciliary targeting and we define a novel ciliary targeting signal: the tetrapeptide RDYR motif in the C-terminus of VPAC2. Furthermore, we demonstrate that VPAC2 ciliary targeting is dependent on Tubby, the BBSome (a complex of Bardet-Biedl syndrome proteins) and the BBSome targeting factor, Arl6.