Inversion variants in the human genome: role in disease and genome architecture

A new era of mutation rate analyses: Concepts and methods

The mutation rate is a pivotal biological characteristic, intricately governed by natural selection and historically garnering considerable attention. Recent advances in high-throughput sequencing and analytical methodologies have profoundly transformed our understanding in this domain, ushering in an unprecedented era of mutation rate research. This paper aims to provide a comprehensive overview of the key concepts and methodologies frequently employed in the study of mutation rates. It examines various types of mutations, explores the evolutionary dynamics and associated theories, and synthesizes both classical and contemporary hypotheses. Furthermore, this review comprehensively explores recent advances in understanding germline and somatic mutations in animals and offers an overview of experimental methodologies, mutational patterns, molecular mechanisms, and driving forces influencing variations in mutation rates across species and tissues. Finally, it proposes several potential research directions and pressing questions for future investigations.

Single variant, yet "double trouble": TSC and KBG syndrome because of a large de novo inversion

Despite the advances in high-throughput sequencing, many rare disease patients remain undiagnosed. In particular, the patients with well-defined clinical phenotypes and established clinical diagnosis, yet missing or partial genetic diagnosis, may hold a clue to more complex genetic mechanisms of a disease that could be missed by available clinical tests. Here, we report a patient with a clinical diagnosis of Tuberous sclerosis, combined with unusual secondary features, but negative clinical tests including TSC1 and TSC2 Short-read whole-genome sequencing combined with advanced bioinformatics analyses were successful in uncovering a de novo pericentric 87-Mb inversion with breakpoints in TSC2 and ANKRD11, which explains the TSC clinical diagnosis, and confirms a second underlying monogenic disorder, KBG syndrome. Our findings illustrate how complex variants, such as large inversions, may be missed by clinical tests and further highlight the importance of well-defined clinical diagnoses in uncovering complex molecular mechanisms of a disease, such as complex variants and "double trouble" effects.

Genome sequencing detects a balanced pericentric inversion with breakpoints that impact the DMD and upstream region of POU3F4 genes

We describe a family with two maternal half-brothers both of whom presented with muscular dystrophy, autism spectrum disorder, developmental delay, and sensorineural hearing loss. The elder brother had onset of features at ~3 months of age, followed by clinical confirmation of muscular dystrophy at 3 years. Skeletal biopsy staining at 4.7 years showed an absence of dystrophin protein which prompted extensive molecular testing over 4 years that included gene panels, targeted single-gene assays, arrays, and karyotyping, all of which failed to identify a clinically significant variant in the DMD gene. At 10 years of age, clinical whole-genome sequencing (cWGS) was performed, which revealed a novel hemizygous ~50.7 Mb balanced pericentric inversion on chromosome X that disrupts the DMD gene in both siblings, consistent with the muscular dystrophy phenotype. This inversion also impacts the upstream regulatory region of POU3F4, structural rearrangements which are known to cause hearing loss. The unaffected mother is a heterozygous carrier for the pericentric inversion. This finding illustrates the ability of cWGS to detect a wide breadth of disease-causing genomic variations including large genomic rearrangements.

Higher chromosomal abnormality rate in blastocysts from a subset of patients with pericentric inversion (Inv) 1 variant

The purpose of this study was to evaluate the incidence of unbalanced chromosome rearrangement in blastocyst-stage embryos from carriers of pericentric inversion of chromosome 1 (PEI-1). A total of 98 embryos from 22 PEI-1 carriers were tested for unbalanced rearrangements, originating from inversion carriers, and overall aneuploidy. Logistic regression analysis indicated that the ratio of inverted segment size to chromosome length was a statistically significant risk factor for unbalanced chromosome rearrangement from PEI-1 carriers (p = 0.003). The optimal cut-off values to predict the risk of unbalanced chromosome rearrangement was 36%, with the incidence being 2.0% in the <36% group and 32.7% in the ≥36% group. The unbalanced embryo rate was 24.4% in male carriers compared to 12.3% in female carriers. Inter-chromosomal effect analysis was performed using 98 blastocysts from PEI-1 carriers and 116 blastocysts from age-matched controls. PEI-1 carriers had similar sporadic aneuploidy rates compared to those of age-matched controls at 32.7 vs. 31.9%, respectively. In conclusion, the risk of unbalanced chromosome rearrangement is affected by inverted segment size in PEI-1 carriers.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.

High level of complexity and global diversity of the 3q29 locus revealed by optical mapping and long-read sequencing

BACKGROUND

High sequence identity between segmental duplications (SDs) can facilitate copy number variants (CNVs) via non-allelic homologous recombination (NAHR). These CNVs are one of the fundamental causes of genomic disorders such as the 3q29 deletion syndrome (del3q29S). There are 21 protein-coding genes lost or gained as a result of such recurrent 1.6-Mbp deletions or duplications, respectively, in the 3q29 locus. While NAHR plays a role in CNV occurrence, the factors that increase the risk of NAHR at this particular locus are not well understood.

METHODS

We employed an optical genome mapping technique to characterize the 3q29 locus in 161 unaffected individuals, 16 probands with del3q29S and their parents, and 2 probands with the 3q29 duplication syndrome (dup3q29S). Long-read sequencing-based haplotype resolved de novo assemblies from 44 unaffected individuals, and 1 trio was used for orthogonal validation of haplotypes and deletion breakpoints.

RESULTS

In total, we discovered 34 haplotypes, of which 19 were novel haplotypes. Among these 19 novel haplotypes, 18 were detected in unaffected individuals, while 1 novel haplotype was detected on the parent-of-origin chromosome of a proband with the del3q29S. Phased assemblies from 44 unaffected individuals enabled the orthogonal validation of 20 haplotypes. In 89% (16/18) of the probands, breakpoints were confined to paralogous copies of a 20-kbp segment within the 3q29 SDs. In one del3q29S proband, the breakpoint was confined to a 374-bp region using long-read sequencing. Furthermore, we categorized del3q29S cases into three classes and dup3q29S cases into two classes based on breakpoints. Finally, we found no evidence of inversions in parent-of-origin chromosomes.

CONCLUSIONS

We have generated the most comprehensive haplotype map for the 3q29 locus using unaffected individuals, probands with del3q29S or dup3q29S, and available parents, and also determined the deletion breakpoint to be within a 374-bp region in one proband with del3q29S. These results should provide a better understanding of the underlying genetic architecture that contributes to the etiology of del3q29S and dup3q29S.

Cytogenetic Investigation of Infertile Patients in Hungary: A 10-Year Retrospective Study

Chromosome abnormalities play a crucial role in reproductive failure. The presence of numerical or structural aberrations may induce recurrent pregnancy loss or primary infertility. The main purpose of our study was to determine the types and frequency of chromosomal aberrations in infertile patients and to compare the frequency of structural aberrations to a control group. Karyotyping was performed in 1489 men and 780 women diagnosed with reproductive failure between 2010 and 2020. The control group included 869 male and 1160 female patients having cytogenetic evaluations for reasons other than infertility. Sex chromosomal aberrations were detected in 33/1489 (2.22%) infertile men and 3/780 (0.38%) infertile women. Structural abnormalities (e.g., translocation, inversion) were observed in 89/1489 (5.98%) infertile men and 58/780 (7.44%) infertile women. The control population showed structural chromosomal abnormalities in 27/869 (3.11%) men and 39/1160 (3.36%) women. There were significant differences in the prevalence of single-cell translocations between infertile individuals (males: 3.5%; females: 3.46%) and control patients (males: 0.46%; females: 0.7%). In summary, this is the first report of cytogenetic alterations in infertile patients in Hungary. The types of chromosomal abnormalities were comparable to previously published data. The prevalence of less-studied single-cell translocations was significantly higher in infertile patients than in the control population, supporting an earlier suggestion that these aberrations may be causally related to infertility.

Genomic architecture and functional effects of potential human inversion supergenes

Supergenes are involved in adaptation in multiple organisms, but they are little known in humans. Genomic inversions are the most common mechanism of supergene generation and maintenance. Here, we review the information about two large inversions that are the best examples of potential human supergenes. In addition, we do an integrative analysis of the newest data to understand better their functional effects and underlying genetic changes. We have found that the highly divergent haplotypes of the 17q21.31 inversion of approximately 1.5 Mb have multiple phenotypic associations, with consistent effects in brain-related traits, red and white blood cells, lung function, male and female characteristics and disease risk. By combining gene expression and nucleotide variation data, we also analysed the molecular differences between haplotypes, including gene duplications, amino acid substitutions and regulatory changes, and identify CRHR1, KANLS1 and MAPT as good candidates to be responsible for these phenotypes. The situation is more complex for the 8p23.1 inversion, where there is no clear genetic differentiation. However, the inversion is associated with several related phenotypes and gene expression differences that could be linked to haplotypes specific of one orientation. Our work, therefore, contributes to the characterization of both exceptional variants and illustrates the important role of inversions. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Correction of a Factor VIII genomic inversion with designer-recombinases

Despite advances in nuclease-based genome editing technologies, correcting human disease-causing genomic inversions remains a challenge. Here, we describe the potential use of a recombinase-based system to correct the 140 kb inversion of the F8 gene frequently found in patients diagnosed with severe Hemophilia A. Employing substrate-linked directed molecular evolution, we develop a coupled heterodimeric recombinase system (RecF8) achieving 30% inversion of the target sequence in human tissue culture cells. Transient RecF8 treatment of endothelial cells, differentiated from patient-derived induced pluripotent stem cells (iPSCs) of a hemophilic donor, results in 12% correction of the inversion and restores Factor VIII mRNA expression. In this work, we present designer-recombinases as an efficient and specific means towards treatment of monogenic diseases caused by large gene inversions.

Correction of disease-causing large genomic inversions remains challenging. Here, the authors developed a dual designer-recombinase system (RecF8) that efficiently corrects a 140 kb inversion frequently found in patients with severe Hemophilia A.

Characteristics and possible mechanisms of formation of microinversions distinguishing human and chimpanzee genomes

Genomic inversions come in various sizes. While long inversions are relatively easy to identify by aligning high-quality genome sequences, unambiguous identification of microinversions is more problematic. Here, using a set of extra stringent criteria to distinguish microinversions from other mutational events, we describe microinversions that occurred after the divergence of humans and chimpanzees. In total, we found 59 definite microinversions that range from 17 to 33 nucleotides in length. In majority of them, human genome sequences matched exactly the reverse-complemented chimpanzee genome sequences, implying that the inverted DNA segment was copied precisely. All these microinversions were flanked by perfect or nearly perfect inverted repeats pointing to their key role in their formation. Template switching at inverted repeats during DNA replication was previously discussed as a possible mechanism for the microinversion formation. However, many of definite microinversions found by us cannot be easily explained via template switching owing to the combination of the short length and imperfect nature of their flanking inverted repeats. We propose a novel, alternative mechanism that involves repair of a double-stranded break within the inverting segment via microhomology-mediated break-induced replication, which can consistently explain all definite microinversion events.

Fine Breakpoint Mapping by Genome Sequencing Reveals the First Large X Inversion Disrupting the NHS Gene in a Patient with Syndromic Cataracts

Inversions are structural variants that are generally balanced. However, they could lead to gene disruptions or have positional effects leading to diseases. Mutations in the NHS gene cause Nance-Horan syndrome, an X-linked disorder characterised by congenital cataracts and dental anomalies. Here, we aimed to characterise a balanced pericentric inversion X(p22q27), maternally inherited, in a child with syndromic bilateral cataracts by breakpoint mapping using whole-genome sequencing (WGS). 30× Illumina paired-end WGS was performed in the proband, and breakpoints were confirmed by Sanger sequencing. EdU assays and FISH analysis were used to assess skewed X-inactivation patterns. RNA expression of involved genes in the breakpoint boundaries was evaluated by droplet-digital PCR. We defined the breakpoint position of the inversion at Xp22.13, with a 15 bp deletion, disrupting the unusually large intron 1 of the canonical NHS isoform, and also perturbing topologically-associated domains (TADs). Moreover, a microhomology region of 5 bp was found on both sides. RNA analysis confirmed null and reduced NHS expression in the proband and his unaffected mother, respectively. In conclusion, we report the first chromosomal inversion disrupting NHS, fine-mapped by WGS. Our data expand the clinical spectrum and the pathogenic mechanisms underlying the NHS defects.

Levodopa-responsive dystonia caused by biallelic PRKN exon inversion invisible to exome sequencing

Biallelic pathogenic variants in PRKN (PARK2), encoding the E3 ubiquitin ligase parkin, lead to early-onset Parkinson's disease. Structural variants, including duplications or deletions, are common in PRKN due to their location within the fragile site FRA6E. These variants are readily detectable by copy number variation analysis. We studied four siblings with levodopa-responsive dystonia by exome sequencing followed by genome sequencing. Affected individuals developed juvenile levodopa-responsive dystonia with subsequent appearance of parkinsonism and motor fluctuations that improved by subthalamic stimulation. Exome sequencing and copy number variation analysis were not diagnostic, yet revealed a shared homozygous block including PRKN. Genome sequencing revealed an inversion within PRKN, with intronic breakpoints flanking exon 5. Breakpoint junction analysis implicated non-homologous end joining and possibly replicative mechanisms as the repair pathways involved. Analysis of cDNA indicated skipping of exon 5 (84 bp) that was replaced by 93 bp of retained intronic sequence, preserving the reading frame yet altering a significant number of residues. Balanced copy number inversions in PRKN are associated with a severe phenotype. Such structural variants, undetected by exome analysis and by copy number variation analysis, should be considered in the relevant clinical setting. These findings raise the possibility that PRKN structural variants are more common than currently estimated.

Chromoanagenesis Event Underlies a de novo Pericentric and Multiple Paracentric Inversions in a Single Chromosome Causing Coffin-Siris Syndrome

Chromoanagenesis is a descriptive term that encompasses classes of catastrophic mutagenic processes that generate localized and complex chromosome rearrangements in both somatic and germline genomes. Herein, we describe a 5-year-old female presenting with a constellation of clinical features consistent with a clinical diagnosis of Coffin–Siris syndrome 1 (CSS1). Initial G-banded karyotyping detected a 90-Mb pericentric and a 47-Mb paracentric inversion on a single chromosome. Subsequent analysis of short-read whole-genome sequencing data and genomic optical mapping revealed additional inversions, all clustered on chromosome 6, one of them disrupting ARID1B for which haploinsufficiency leads to the CSS1 disease trait (MIM:135900). The aggregate structural variant data show that the resolved, the resolved derivative chromosome architecture presents four de novo inversions, one pericentric and three paracentric, involving six breakpoint junctions in what appears to be a shuffling of genomic material on this chromosome. Each junction was resolved to nucleotide-level resolution with mutational signatures suggestive of non-homologous end joining. The disruption of the gene ARID1B is shown to occur between the fourth and fifth exon of the canonical transcript with subsequent qPCR studies confirming a decrease in ARID1B expression in the patient versus healthy controls. Deciphering the underlying genomic architecture of chromosomal rearrangements and complex structural variants may require multiple technologies and can be critical to elucidating the molecular etiology of a patient’s clinical phenotype or resolving unsolved Mendelian disease cases.

InvertypeR: Bayesian inversion genotyping with Strand-seq data

Background

Single cell Strand-seq is a unique tool for the discovery and phasing of genomic inversions. Conventional methods to discover inversions with Strand-seq data are blind to known inversion locations, limiting their statistical power for the detection of inversions smaller than 10 Kb. Moreover, the methods rely on manual inspection to separate false and true positives.

Results

Here we describe “InvertypeR”, a method based on a Bayesian binomial model that genotypes inversions using fixed genomic coordinates. We validated InvertypeR by re-genotyping inversions reported for three trios by the Human Genome Structural Variation Consortium. Although 6.3% of the family inversion genotypes in the original study showed Mendelian discordance, this was reduced to 0.5% using InvertypeR. By applying InvertypeR to published inversion coordinates and predicted inversion hotspots (n = 3701), as well as coordinates from conventional inversion discovery, we furthermore genotyped 66 inversions not previously reported for the three trios.

Conclusions

InvertypeR discovers, genotypes, and phases inversions without relying on manual inspection. For greater accessibility, results are presented as phased chromosome ideograms with inversions linked to Strand-seq data in the genome browser. InvertypeR increases the power of Strand-seq for studies on the role of inversions in phenotypic variation, genome instability, and human disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07892-9.

Chromosomal Aberrations in 224 Couples with Recurrent Pregnancy Loss

Background

Recurrent pregnancy loss (RPL) is a major reproductive health issue, affecting 2%-5% of couples. Genetic factors, mainly chromosomal abnormalities, are the most common cause of early miscarriage accounting for 50%-60% of first trimester abortion.

Aim

To estimate the prevalence and nature of chromosomal anomalies in couples with recurrent miscarriage.

Patients and Methods

This study included 224 couples with a history of 2 or more abortions. Both partners were karyotyped as part of the primary investigation. Cytogenetic analysis was carried out using the standard method.

Results

A total of 224 couples with a history of two or more recurrent abortions were enrolled in this study. Chromosomal abnormalities were detected in 26 couples (11.6%) and 28 individuals (6.25%). We found a structural chromosome abnormality in 17/28 patients (60.7%); 12 patients had a reciprocal translocation (42.9%) including one patient with an additional inversion of the Y chromosome, 4 (14.3%) had a Robertsonian translocation, and one patient (3.6%) carried a paracentric inversion of chromosome 2. Numerical chromosome aberrations were detected in 5 patients; three patients (10.7%) with sex chromosome abnormalities and two (7.1%) with a marker chromosome. Six patients (21.4%) showed a heteromorphic variant involving chromosome 9.

Conclusion

The prevalence of chromosomal abnormalities in couples with RPL is within the range reported worldwide. Cytogenetic analysis should become an integral part of the investigations of couples with at least two pregnancy losses of undetermined etiology.

Directional Genomic Hybridization (dGH) for Detection of Intrachromosomal Rearrangements

Fluorescence in situ Hybridization (FISH) techniques, including whole chromosome painting (WCP), spectral karyotyping (SKY), and multicolor FISH (mFISH), are used extensively to characterize and enumerate inter-chromosomal rearrangements (e.g., translocations). Directional genomic hybridization (dGH) is a relatively new cytogenomics-based methodology that combines the strand-specific strategy of Chromosome Orientation-FISH (CO-FISH) with bioinformatics-driven design of single-stranded DNA probe sets that are unique and of like orientation. Such a strategy produces directional probe sets that hybridize to one-and only one-chromatid of prepared (single-stranded) metaphase chromosomes, thereby facilitating high-resolution visualization of intra-chromosomal rearrangements, specifically inversions, and greatly improving our ability to detect such otherwise cryptic structural variants within the genome. In addition to its usefulness in the study of various disease states, including cancer, relevant applications of dGH include monitoring cytogenetic damage caused by exposure to clastogenic agents (e.g., ionizing radiation). dGH can be applied as a discovery tool to globally assess the integrity of the genome, but it can also be used in a more targeted fashion to interrogate fine structural changes at the kilobase level. Consequently, dGH is capable of providing significant mechanistic insight and information not easily obtainable by other approaches.

De Deyn P, C. van Swieten J, M. van Duijn C, van der Zee J, Sleegers K, Van Broeckhoven C. Loss of DPP6 in neurodegenerative dementia: a genetic player in the dysfunction of neuronal excitability

Emerging evidence suggested a converging mechanism in neurodegenerative brain diseases (NBD) involving early neuronal network dysfunctions and alterations in the homeostasis of neuronal firing as culprits of neurodegeneration. In this study, we used paired-end short-read and direct long-read whole genome sequencing to investigate an unresolved autosomal dominant dementia family significantly linked to 7q36. We identified and validated a chromosomal inversion of ca. 4 Mb, segregating on the disease haplotype and disrupting the coding sequence of dipeptidyl-peptidase 6 gene (DPP6). DPP6 resequencing identified significantly more rare variants-nonsense, frameshift, and missense-in early-onset Alzheimer's disease (EOAD, p value = 0.03, OR = 2.21 95% CI 1.05-4.82) and frontotemporal dementia (FTD, p = 0.006, OR = 2.59, 95% CI 1.28-5.49) patient cohorts. DPP6 is a type II transmembrane protein with a highly structured extracellular domain and is mainly expressed in brain, where it binds to the potassium channel Kv4.2 enhancing its expression, regulating its gating properties and controlling the dendritic excitability of hippocampal neurons. Using in vitro modeling, we showed that the missense variants found in patients destabilize DPP6 and reduce its membrane expression (p < 0.001 and p < 0.0001) leading to a loss of protein. Reduced DPP6 and/or Kv4.2 expression was also detected in brain tissue of missense variant carriers. Loss of DPP6 is known to cause neuronal hyperexcitability and behavioral alterations in Dpp6-KO mice. Taken together, the results of our genomic, genetic, expression and modeling analyses, provided direct evidence supporting the involvement of DPP6 loss in dementia. We propose that loss of function variants have a higher penetrance and disease impact, whereas the missense variants have a variable risk contribution to disease that can vary from high to low penetrance. Our findings of DPP6, as novel gene in dementia, strengthen the involvement of neuronal hyperexcitability and alteration in the homeostasis of neuronal firing as a disease mechanism to further investigate.

A comparison of the association between large haplotype blocks under selection and the presence/absence of inversions

Inversions may contribute to ecological adaptation and phenotypic diversity, and with the advent of "second" and "third" generation sequencing technologies, the ability to detect inversion polymorphisms has been enhanced dramatically. A key molecular consequence of an inversion is the suppression of recombination allowing independent accumulation of genetic changes between alleles over time. This may lead to the development of divergent haplotype blocks maintained by balancing selection. Thus, divergent haplotype blocks are often considered as indicating the presence of an inversion. In this paper, we first review the features of a 7.7 Mb inversion causing the Rose-comb phenotype in chicken, as a model for how inversions evolve and directly affect phenotypes. Second, we compare the genetic basis for alternative mating strategies in ruff and timing of reproduction in Atlantic herring, both associated with divergent haplotype blocks. Alternative male mating strategies in ruff are associated with a 4.5 Mb inversion that occurred about 4 million years ago. In fact, the ruff inversion shares some striking features with the Rose-comb inversion such as disruption of a gene at one of the inversion breakpoints and generation of a new allele by recombination between the inverted and noninverted alleles. In contrast, inversions do not appear to be a major reason for the fairly large haplotype blocks (range 10-200 kb) associated with ecological adaptation in the herring. Thus, it is important to note that divergent haplotypes may also be maintained by natural selection in the absence of structural variation.

Infertility patients with chromosome inversions are not susceptible to an inter-chromosomal effect

PURPOSE

The aim of this study was to evaluate the incidence of an inter-chromosomal effect (ICE) in blastocyst-stage embryos from carriers of balanced chromosome inversions.

METHODS

Infertility patients (n = 52) with balanced inversions (n = 66 cycles), and maternal age-matched controls that concurrently cycled (n = 66), consented to an IVF cycle with preimplantation genetic testing for aneuploidy (PGT-A). Blastocyst-stage embryos underwent trophectoderm biopsy for PGT-A with only euploid blastocysts transferred in a subsequent frozen embryo transfer. Subtypes of inversions were included in aggregate: paracentric/pericentric, polymorphic/non-polymorphic, male/female carriers, and varying inversion sizes.

RESULTS

The incidence of aneuploidy was not significantly higher for the inversion patients compared to the controls (inversion = 48.8% vs. control = 47.2% ns). Following euploid blastocyst transfer, there were excellent live birth outcomes.

CONCLUSIONS

Carriers of balanced chromosome inversions did not exhibit higher aneuploidy rates for chromosomes that were not involved in the inversion compared to maternal age-matched controls, signifying the absence of an inter-chromosomal effect for this data set. These results provide the largest investigation of blastocyst embryos regarding the debated existence of an ICE resulting from the presence of an inversion during meiosis. However, further studies are warranted to investigate an ICE among inversions subtypes that were outside the scope of this study.

The Effect of Common Inversion Polymorphisms In(2L)t and In(3R)Mo on Patterns of Transcriptional Variation in Drosophila melanogaster

Chromosomal inversions are a ubiquitous feature of genetic variation. Theoretical models describe several mechanisms by which inversions can drive adaptation and be maintained as polymorphisms. While inversions have been shown previously to be under selection, or contain genetic variation under selection, the specific phenotypic consequences of inversions leading to their maintenance remain unclear. Here we use genomic sequence and expression data from the Drosophila Genetic Reference Panel (DGRP) to explore the effects of two cosmopolitan inversions, In(2L)t and In(3R)Mo, on patterns of transcriptional variation. We demonstrate that each inversion has a significant effect on transcript abundance for hundreds of genes across the genome. Inversion-affected loci (IAL) appear both within inversions as well as on unlinked chromosomes. Importantly, IAL do not appear to be influenced by the previously reported genome-wide expression correlation structure. We found that five genes involved with sterol uptake, four of which are Niemann-Pick Type 2 orthologs, are upregulated in flies with In(3R)Mo but do not have SNPs in linkage disequilibrium (LD) with the inversion. We speculate that this upregulation is driven by genetic variation in mod(mdg4) that is in LD with In(3R)Mo. We find that there is little evidence for a regional or position effect of inversions on gene expression at the chromosomal level, but do find evidence for the distal breakpoint of In(3R)Mo interrupting one gene and possibly disassociating the two flanking genes from regulatory elements.

Detailed analysis of inversions predicted between two human genomes: errors, real polymorphisms, and their origin and population distribution

The growing catalogue of structural variants in humans often overlooks inversions as one of the most difficult types of variation to study, even though they affect phenotypic traits in diverse organisms. Here, we have analysed in detail 90 inversions predicted from the comparison of two independently assembled human genomes: the reference genome (NCBI36/HG18) and HuRef. Surprisingly, we found that two thirds of these predictions (62) represent errors either in assembly comparison or in one of the assemblies, including 27 misassembled regions in HG18. Next, we validated 22 of the remaining 28 potential polymorphic inversions using different PCR techniques and characterized their breakpoints and ancestral state. In addition, we determined experimentally the derived allele frequency in Europeans for 17 inversions (DAF = 0.01-0.80), as well as the distribution in 14 worldwide populations for 12 of them based on the 1000 Genomes Project data. Among the validated inversions, nine have inverted repeats (IRs) at their breakpoints, and two show nucleotide variation patterns consistent with a recurrent origin. Conversely, inversions without IRs have a unique origin and almost all of them show deletions or insertions at the breakpoints in the derived allele mediated by microhomology sequences, which highlights the importance of mechanisms like FoSTeS/MMBIR in the generation of complex rearrangements in the human genome. Finally, we found several inversions located within genes and at least one candidate to be positively selected in Africa. Thus, our study emphasizes the importance of careful analysis and validation of large-scale genomic predictions to extract reliable biological conclusions.

Single-cell template strand sequencing by Strand-seq enables the characterization of individual homologs

The ability to distinguish between genome sequences of homologous chromosomes in single cells is important for studies of copy-neutral genomic rearrangements (such as inversions and translocations), building chromosome-length haplotypes, refining genome assemblies, mapping sister chromatid exchange events and exploring cellular heterogeneity. Strand-seq is a single-cell sequencing technology that resolves the individual homologs within a cell by restricting sequence analysis to the DNA template strands used during DNA replication. This protocol, which takes up to 4 d to complete, relies on the directionality of DNA, in which each single strand of a DNA molecule is distinguished based on its 5'-3' orientation. Culturing cells in a thymidine analog for one round of cell division labels nascent DNA strands, allowing for their selective removal during genomic library construction. To preserve directionality of template strands, genomic preamplification is bypassed and labeled nascent strands are nicked and not amplified during library preparation. Each single-cell library is multiplexed for pooling and sequencing, and the resulting sequence data are aligned, mapping to either the minus or plus strand of the reference genome, to assign template strand states for each chromosome in the cell. The major adaptations to conventional single-cell sequencing protocols include harvesting of daughter cells after a single round of BrdU incorporation, bypassing of whole-genome amplification, and removal of the BrdU⁺ strand during Strand-seq library preparation. By sequencing just template strands, the structure and identity of each homolog are preserved.

Characterizing polymorphic inversions in human genomes by single-cell sequencing

Identifying genomic features that differ between individuals and cells can help uncover the functional variants that drive phenotypes and disease susceptibilities. For this, single-cell studies are paramount, as it becomes increasingly clear that the contribution of rare but functional cellular subpopulations is important for disease prognosis, management, and progression. Until now, studying these associations has been challenged by our inability to map structural rearrangements accurately and comprehensively. To overcome this, we coupled single-cell sequencing of DNA template strands (Strand-seq) with custom analysis software to rapidly discover, map, and genotype genomic rearrangements at high resolution. This allowed us to explore the distribution and frequency of inversions in a heterogeneous cell population, identify several polymorphic domains in complex regions of the genome, and locate rare alleles in the reference assembly. We then mapped the entire genomic complement of inversions within two unrelated individuals to characterize their distinct inversion profiles and built a nonredundant global reference of structural rearrangements in the human genome. The work described here provides a powerful new framework to study structural variation and genomic heterogeneity in single-cell samples, whether from individuals for population studies or tissue types for biomarker discovery.

Genome Editing of Structural Variations: Modeling and Gene Correction

The analysis of chromosomal structural variations (SVs), such as inversions and translocations, was made possible by the completion of the human genome project and the development of genome-wide sequencing technologies. SVs contribute to genetic diversity and evolution, although some SVs can cause diseases such as hemophilia A in humans. Genome engineering technology using programmable nucleases (e.g., ZFNs, TALENs, and CRISPR/Cas9) has been rapidly developed, enabling precise and efficient genome editing for SV research. Here, we review advances in modeling and gene correction of SVs, focusing on inversion, translocation, and nucleotide repeat expansion.

Identifying micro-inversions using high-throughput sequencing reads

BACKGROUND

The identification of inversions of DNA segments shorter than read length (e.g., 100 bp), defined as micro-inversions (MIs), remains challenging for next-generation sequencing reads. It is acknowledged that MIs are important genomic variation and may play roles in causing genetic disease. However, current alignment methods are generally insensitive to detect MIs. Here we develop a novel tool, MID (Micro-Inversion Detector), to identify MIs in human genomes using next-generation sequencing reads.

RESULTS

The algorithm of MID is designed based on a dynamic programming path-finding approach. What makes MID different from other variant detection tools is that MID can handle small MIs and multiple breakpoints within an unmapped read. Moreover, MID improves reliability in low coverage data by integrating multiple samples. Our evaluation demonstrated that MID outperforms Gustaf, which can currently detect inversions from 30 bp to 500 bp.

CONCLUSIONS

To our knowledge, MID is the first method that can efficiently and reliably identify MIs from unmapped short next-generation sequencing reads. MID is reliable on low coverage data, which is suitable for large-scale projects such as the 1000 Genomes Project (1KGP). MID identified previously unknown MIs from the 1KGP that overlap with genes and regulatory elements in the human genome. We also identified MIs in cancer cell lines from Cancer Cell Line Encyclopedia (CCLE). Therefore our tool is expected to be useful to improve the study of MIs as a type of genetic variant in the human genome. The source code can be downloaded from: http://cqb.pku.edu.cn/ZhuLab/MID .

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.

Chromosomal polymorphism in mammals: an evolutionary perspective

Although chromosome rearrangements (CRs) are central to studies of genome evolution, our understanding of the evolutionary consequences of the early stages of karyotypic differentiation (i.e. polymorphism), especially the non-meiotic impacts, is surprisingly limited. We review the available data on chromosomal polymorphisms in mammals so as to identify taxa that hold promise for developing a more comprehensive understanding of chromosomal change. In doing so, we address several key questions: (i) to what extent are mammalian karyotypes polymorphic, and what types of rearrangements are principally involved? (ii) Are some mammalian lineages more prone to chromosomal polymorphism than others? More specifically, do (karyotypically) polymorphic mammalian species belong to lineages that are also characterized by past, extensive karyotype repatterning? (iii) How long can chromosomal polymorphisms persist in mammals? We discuss the evolutionary implications of these questions and propose several research avenues that may shed light on the role of chromosome change in the diversification of mammalian populations and species.

Exome sequencing and whole genome sequencing for the detection of copy number variation

Many laboratories now use genomic microarrays as their first-tier diagnostic test for copy number variation (CNV) detection. In addition, whole exome sequencing is increasingly being offered as a diagnostic test for heterogeneous disorders. Although mostly used for the detection of point mutations and small insertion-deletions, exome sequencing can also be used to call CNVs, allowing combined small and large variant analysis. Whole genome sequencing in addition to these advantages also offers the potential to characterize CNVs to unprecedented levels of accuracy, providing position and orientation information. In this review, we discuss the clinical potential of CNV identification in whole exome sequencing and whole genome sequencing data and the implications this has on diagnostic laboratories.

A decade of structural variants: description, history and methods to detect structural variation

In the past decade, the view on genomic structural variation (SV) has been changed completely. SVs, previously considered rare events, are now recognized as the largest source of interindividual genetic variation affecting more bases than single nucleotide polymorphisms, variable number of tandem repeats and other small genetic variants. They have also been shown to play a role in phenotypic variation and in disease. In this review, the authors will provide an introduction to SV; a short historical perspective on the research of this source of genomic variation; a description of the types of structural variants, and on how they may have arisen; and an overview on methods of detecting structural variants, focusing on the analysis of high-throughput sequencing data.

Following the footprints of polymorphic inversions on SNP data: from detection to association tests

Inversion polymorphisms have important phenotypic and evolutionary consequences in humans. Two different methodologies have been used to infer inversions from SNP dense data, enabling the use of large cohorts for their study. One approach relies on the differences in linkage disequilibrium across breakpoints; the other one captures the internal haplotype groups that tag the inversion status of chromosomes. In this article, we assessed the convergence of the two methods in the detection of 20 human inversions that have been reported in the literature. The methods converged in four inversions including inv-8p23, for which we studied its association with low-BMI in American children. Using a novel haplotype tagging method with control on inversion ancestry, we computed the frequency of inv-8p23 in two American cohorts and observed inversion haplotype admixture. Accounting for haplotype ancestry, we found that the European inverted allele in children carries a recessive risk of underweight, validated in an independent Spanish cohort (combined: OR= 2.00, P = 0.001). While the footprints of inversions on SNP data are complex, we show that systematic analyses, such as convergence of different methods and controlling for ancestry, can reveal the contribution of inversions to the ancestral composition of populations and to the heritability of human disease.

Making the difference: integrating structural variation detection tools

From prokaryotes to eukaryotes, phenotypic variation, adaptation and speciation has been associated with structural variation between genomes of individuals within the same species. Many computer algorithms detecting such variations (callers) have recently been developed, spurred by the advent of the next-generation sequencing technology. Such callers mainly exploit split-read mapping or paired-end read mapping. However, as different callers are geared towards different types of structural variation, there is still no single caller that can be considered a community standard; instead, increasingly the various callers are combined in integrated pipelines. In this article, we review a wide range of callers, discuss challenges in the integration step and present a survey of pipelines used in population genomics studies. Based on our findings, we provide general recommendations on how to set-up such pipelines. Finally, we present an outlook on future challenges in structural variation detection.

A 15 Mb large paracentric chromosome 21 inversion identified in Czech population through a pair of flanking duplications

Background

Inversions are balanced structural chromosome rearrangements, which can influence gene expression and the risk of unbalanced chromosome constitution in offspring. Many examples of inversion polymorphisms exist in human, affecting both heterochromatic regions and euchromatin.

Results

We describe a novel, 15 Mb long paracentric inversion, inv(21)(q21.1q22.11), affecting more than a third of human 21q. Despite of its length, the inversion cannot be detected using karyotyping due to similar band patterns on the normal and inverted chromosomes, and is therefore likely to escape attention. Its identification was aided by the repeated observation of the same pair of 150 kb long duplications present in cis on chromosome 21 in three Czech families subjected to microarray analysis. The finding prompted us to hypothesise that this co-occurrence of two remote duplications could be associated with an inversion of the intervening segment, and this speculation turned out to be right. The inversion was confirmed in a series of FISH experiments which also showed that the second copy of each of the duplications was always located at the opposite end of the inversion. The presence of the same pair of duplications in additional individuals reported in public databases indicates that the inversion may also be present in other populations. Three out of the total of about 4000 chromosomes 21 examined in our sample carried the duplications and were inverted, corresponding to carrier frequency of about 1/660. Although the breakpoints affect protein-coding genes, the occurrence of the inversion in normal parents and siblings of our patients and the occurrence of the duplications in unaffected controls in databases indicate that this rare variant is rather non-pathogenic. The inverted segment carried an identical shared haplotype in the three families studied. The haplotypes, however, diverged very rapidly in the flanking regions, possibly pointing to an ancient founder event at the origin of the inversion.

Conclusions

The identification of inv(21)(q21.1q22.11) supports the notion that paracentric inversions are the most common form of chromosomal variation and that some of them may still remain undetected.

Genome-wide association tests of inversions with application to psoriasis

Although inversions have occasionally been found to be associated with disease susceptibility through interrupting a gene or its regulatory region, or by increasing the risk for deleterious secondary rearrangements, no association study has been specifically conducted for risks associated with inversions, mainly because existing approaches to detecting and genotyping inversions do not readily scale to a large number of samples. Based on our recently proposed approach to identifying and genotyping inversions using principal components analysis (PCA), we herein develop a method of detecting association between inversions and disease in a genome-wide fashion. Our method uses genotype data for single nucleotide polymorphisms (SNPs), and is thus cost-efficient and computationally fast. For an inversion polymorphism, local PCA around the inversion region is performed to infer the inversion genotypes of all samples. For many inversions, we found that some of the SNPs inside an inversion region are fixed in the two lineages of different orientations and thus can serve as surrogate markers. Our method can be applied to case-control and quantitative trait association studies to identify inversions that may interrupt a gene or the connection between a gene and its regulatory agents. Our method also offers a new venue to identify inversions that are responsible for disease-causing secondary rearrangements. We illustrated our proposed approach to case-control data for psoriasis and identified novel associations with a few inversion polymorphisms.

Validation and genotyping of multiple human polymorphic inversions mediated by inverted repeats reveals a high degree of recurrence

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.

InvFEST, a database integrating information of polymorphic inversions in the human genome

The newest genomic advances have uncovered an unprecedented degree of structural variation throughout genomes, with great amounts of data accumulating rapidly. Here we introduce InvFEST (http://invfestdb.uab.cat), a database combining multiple sources of information to generate a complete catalogue of non-redundant human polymorphic inversions. Due to the complexity of this type of changes and the underlying high false-positive discovery rate, it is necessary to integrate all the available data to get a reliable estimate of the real number of inversions. InvFEST automatically merges predictions into different inversions, refines the breakpoint locations, and finds associations with genes and segmental duplications. In addition, it includes data on experimental validation, population frequency, functional effects and evolutionary history. All this information is readily accessible through a complete and user-friendly web report for each inversion. In its current version, InvFEST combines information from 34 different studies and contains 1092 candidate inversions, which are categorized based on internal scores and manual curation. Therefore, InvFEST aims to represent the most reliable set of human inversions and become a central repository to share information, guide future studies and contribute to the analysis of the functional and evolutionary impact of inversions on the human genome.

On the structural plasticity of the human genome: chromosomal inversions revisited

With the aid of novel and powerful molecular biology techniques, recent years have witnessed a dramatic increase in the number of studies reporting the involvement of complex structural variants in several genomic disorders. In fact, with the discovery of Copy Number Variants (CNVs) and other forms of unbalanced structural variation, much attention has been directed to the detection and characterization of such rearrangements, as well as the identification of the mechanisms involved in their formation. However, it has long been appreciated that chromosomes can undergo other forms of structural changes - balanced rearrangements - that do not involve quantitative variation of genetic material. Indeed, a particular subtype of balanced rearrangement – inversions – was recently found to be far more common than had been predicted from traditional cytogenetics. Chromosomal inversions alter the orientation of a specific genomic sequence and, unless involving breaks in coding or regulatory regions (and, disregarding complex trans effects, in their close vicinity), appear to be phenotypically silent. Such a surprising finding, which is difficult to reconcile with the classical interpretation of inversions as a mechanism causing subfertility (and ultimately reproductive isolation), motivated a new series of theoretical and empirical studies dedicated to understand their role in human genome evolution and to explore their possible association to complex genetic disorders. With this review, we attempt to describe the latest methodological improvements to inversions detection at a genome wide level, while exploring some of the possible implications of inversion rearrangements on the evolution of the human genome.

On the power and the systematic biases of the detection of chromosomal inversions by paired-end genome sequencing

One of the most used techniques to study structural variation at a genome level is paired-end mapping (PEM). PEM has the advantage of being able to detect balanced events, such as inversions and translocations. However, inversions are still quite difficult to predict reliably, especially from high-throughput sequencing data. We simulated realistic PEM experiments with different combinations of read and library fragment lengths, including sequencing errors and meaningful base-qualities, to quantify and track down the origin of false positives and negatives along sequencing, mapping, and downstream analysis. We show that PEM is very appropriate to detect a wide range of inversions, even with low coverage data. However, ≥% of inversions located between segmental duplications are expected to go undetected by the most common sequencing strategies. In general, longer DNA libraries improve the detectability of inversions far better than increments of the coverage depth or the read length. Finally, we review the performance of three algorithms to detect inversions--SVDetect, GRIAL, and VariationHunter--, identify common pitfalls, and reveal important differences in their breakpoint precisions. These results stress the importance of the sequencing strategy for the detection of structural variants, especially inversions, and offer guidelines for the design of future genome sequencing projects.

Directional genomic hybridization for chromosomal inversion discovery and detection

Chromosomal rearrangements are a source of structural variation within the genome that figure prominently in human disease, where the importance of translocations and deletions is well recognized. In principle, inversions—reversals in the orientation of DNA sequences within a chromosome—should have similar detrimental potential. However, the study of inversions has been hampered by traditional approaches used for their detection, which are not particularly robust. Even with significant advances in whole genome approaches, changes in the absolute orientation of DNA remain difficult to detect routinely. Consequently, our understanding of inversions is still surprisingly limited, as is our appreciation for their frequency and involvement in human disease. Here, we introduce the directional genomic hybridization methodology of chromatid painting—a whole new way of looking at structural features of the genome—that can be employed with high resolution on a cell-by-cell basis, and demonstrate its basic capabilities for genome-wide discovery and targeted detection of inversions. Bioinformatics enabled development of sequence- and strand-specific directional probe sets, which when coupled with single-stranded hybridization, greatly improved the resolution and ease of inversion detection. We highlight examples of the far-ranging applicability of this cytogenomics-based approach, which include confirmation of the alignment of the human genome database and evidence that individuals themselves share similar sequence directionality, as well as use in comparative and evolutionary studies for any species whose genome has been sequenced. In addition to applications related to basic mechanistic studies, the information obtainable with strand-specific hybridization strategies may ultimately enable novel gene discovery, thereby benefitting the diagnosis and treatment of a variety of human disease states and disorders including cancer, autism, and idiopathic infertility.

Genome-wide identification of structural variants in genes encoding drug targets: possible implications for individualized drug therapy

OBJECTIVE

The objective of the present study was to identify structural variants of drug target-encoding genes on a genome-wide scale. We also aimed at identifying drugs that are potentially amenable for individualization of treatments based on knowledge about structural variation in the genes encoding their targets.

METHODS

Information about human drug targets of all therapeutic drugs and nutraceuticals approved by the Food and Drug Administration and with an Anatomical Therapeutic Chemical (ATC) code, namely, 876, was obtained from the DrugBank and applied to interrogate the Database of Genomic Variants.

RESULTS

We identified 1721 structural variants, which affected 495 of 1005 different genes encoding drug targets. About 70% of the Food and Drug Administration-approved drugs had targets subjected to structural variation, in particular copy number variation. The distribution of drugs with targets encoded by copy number variable genes differed between ATC groups with pronounced over-representation in ATC group N, that is, nervous system drugs (P=3.75e-5). Several narrow therapeutic index drugs with targets subjected to genomic structural variation were identified. Information about the frequencies of the structural variants and drug consumption allowed us to identify potential candidates for implementation in pharmacogenetic testing and individualized drug treatments.

CONCLUSION

Structural variants in pharmacodynamic genes may play a larger role in determining individual variation in drug responses than currently believed. Inclusion of such genes in pharmacogenetic testing holds promising prospects as they may have large effect sizes.

Statistical challenges associated with detecting copy number variations with next-generation sequencing

MOTIVATION

Analysing next-generation sequencing (NGS) data for copy number variations (CNVs) detection is a relatively new and challenging field, with no accepted standard protocols or quality control measures so far. There are by now several algorithms developed for each of the four broad methods for CNV detection using NGS, namely the depth of coverage (DOC), read-pair, split-read and assembly-based methods. However, because of the complexity of the genome and the short read lengths from NGS technology, there are still many challenges associated with the analysis of NGS data for CNVs, no matter which method or algorithm is used.

RESULTS

In this review, we describe and discuss areas of potential biases in CNV detection for each of the four methods. In particular, we focus on issues pertaining to (i) mappability, (ii) GC-content bias, (iii) quality control measures of reads and (iv) difficulty in identifying duplications. To gain insights to some of the issues discussed, we also download real data from the 1000 Genomes Project and analyse its DOC data. We show examples of how reads in repeated regions can affect CNV detection, demonstrate current GC-correction algorithms, investigate sensitivity of DOC algorithm before and after quality control of reads and discuss reasons for which duplications are harder to detect than deletions.

Recent advances in targeted genome engineering in mammalian systems

Targeted genome engineering enables researchers to disrupt, insert, or replace a genomic sequence precisely at a predetermined locus. One well-established technology to edit a mammalian genome is known as gene targeting, which is based on the homologous recombination (HR) mechanism. However, the low HR frequency in mammalian cells (except for mice) prevents its wide application. To address this limitation, a custom-designed nuclease is used to introduce a site-specific DNA double-strand break (DSB) on the chromosome and the subsequent repair of the DSB by the HR mechanism or the non-homologous end joining mechanism results in efficient targeted genome modifications. Engineered homing endonucleases (also called meganucleases), zinc finger nucleases, and transcription activator-like effector nucleases represent the three major classes of custom-designed nucleases that have been successfully applied in many different organisms for targeted genome engineering. This article reviews the recent developments of these genome engineering tools and highlights a few representative applications in mammalian systems. Recent advances in gene delivery strategies of these custom-designed nucleases are also briefly discussed.

Investigation of inversion polymorphisms in the human genome using principal components analysis

Despite the significant advances made over the last few years in mapping inversions with the advent of paired-end sequencing approaches, our understanding of the prevalence and spectrum of inversions in the human genome has lagged behind other types of structural variants, mainly due to the lack of a cost-efficient method applicable to large-scale samples. We propose a novel method based on principal components analysis (PCA) to characterize inversion polymorphisms using high-density SNP genotype data. Our method applies to non-recurrent inversions for which recombination between the inverted and non-inverted segments in inversion heterozygotes is suppressed due to the loss of unbalanced gametes. Inside such an inversion region, an effect similar to population substructure is thus created: two distinct “populations” of inversion homozygotes of different orientations and their 1∶1 admixture, namely the inversion heterozygotes. This kind of substructure can be readily detected by performing PCA locally in the inversion regions. Using simulations, we demonstrated that the proposed method can be used to detect and genotype inversion polymorphisms using unphased genotype data. We applied our method to the phase III HapMap data and inferred the inversion genotypes of known inversion polymorphisms at 8p23.1 and 17q21.31. These inversion genotypes were validated by comparing with literature results and by checking Mendelian consistency using the family data whenever available. Based on the PCA-approach, we also performed a preliminary genome-wide scan for inversions using the HapMap data, which resulted in 2040 candidate inversions, 169 of which overlapped with previously reported inversions. Our method can be readily applied to the abundant SNP data, and is expected to play an important role in developing human genome maps of inversions and exploring associations between inversions and susceptibility of diseases.

Human genetics and genomics a decade after the release of the draft sequence of the human genome

Substantial progress has been made in human genetics and genomics research over the past ten years since the publication of the draft sequence of the human genome in 2001. Findings emanating directly from the Human Genome Project, together with those from follow-on studies, have had an enormous impact on our understanding of the architecture and function of the human genome. Major developments have been made in cataloguing genetic variation, the International HapMap Project, and with respect to advances in genotyping technologies. These developments are vital for the emergence of genome-wide association studies in the investigation of complex diseases and traits. In parallel, the advent of high-throughput sequencing technologies has ushered in the 'personal genome sequencing' era for both normal and cancer genomes, and made possible large-scale genome sequencing studies such as the 1000 Genomes Project and the International Cancer Genome Consortium. The high-throughput sequencing and sequence-capture technologies are also providing new opportunities to study Mendelian disorders through exome sequencing and whole-genome sequencing. This paper reviews these major developments in human genetics and genomics over the past decade.

The origin, global distribution, and functional impact of the human 8p23 inversion polymorphism

Genomic inversions are an increasingly recognized source of genetic variation. However, a lack of reliable high-throughput genotyping assays for these structures has precluded a full understanding of an inversion's phylogenetic, phenotypic, and population genetic properties. We characterize these properties for one of the largest polymorphic inversions in man (the ∼4.5-Mb 8p23.1 inversion), a structure that encompasses numerous signals of natural selection and disease association. We developed and validated a flexible bioinformatics tool that utilizes SNP data to enable accurate, high-throughput genotyping of the 8p23.1 inversion. This tool was applied retrospectively to diverse genome-wide data sets, revealing significant population stratification that largely follows a clinal “serial founder effect” distribution model. Phylogenetic analyses establish the inversion's ancestral origin within the Homo lineage, indicating that 8p23.1 inversion has occurred independently in the Pan lineage. The human inversion breakpoint was localized to an inverted pair of human endogenous retrovirus elements within the large, flanking low-copy repeats; experimental validation of this breakpoint confirmed these elements as the likely intermediary substrates that sponsored inversion formation. In five data sets, mRNA levels of disease-associated genes were robustly associated with inversion genotype. Moreover, a haplotype associated with systemic lupus erythematosus was restricted to the derived inversion state. We conclude that the 8p23.1 inversion is an evolutionarily dynamic structure that can now be accommodated into the understanding of human genetic and phenotypic diversity.

High rates of de novo 15q11q13 inversions in human spermatozoa

Low-Copy Repeats predispose the 15q11-q13 region to non-allelic homologous recombination. We have already demonstrated that a significant percentage of Prader-Willi syndrome (PWS) fathers have an increased susceptibility to generate 15q11q13 deletions in spermatozoa, suggesting the participation of intrachromatid exchanges. This work has been focused on assessing the incidence of de novo 15q11q13 inversions in spermatozoa of control donors and PWS fathers in order to determine the basal rates of inversions and to confirm the intrachromatid mechanism as the main cause of 15q11q13 anomalies.

Semen samples from 10 control donors and 16 PWS fathers were processed and analyzed by triple-color FISH. Three differentially labeled BAC-clones were used: one proximal and two distal of the 15q11-q13 region. Signal associations allowed the discrimination between normal and inverted haplotypes, which were confirmed by laser-scanning confocal microscopy.

Two types of inversions were detected which correspond to the segments involved in Class I and II PWS deletions. No significant differences were observed in the mean frequencies of inversions between controls and PWS fathers (3.59% ± 0.46 and 9.51% ± 0.87 vs 3.06% ± 0.33 and 10.07% ± 0.74). Individual comparisons showed significant increases of inversions in four PWS fathers (P < 0.05) previously reported as patients with increases of 15q11q13 deletions.

Results suggest that the incidence of heterozygous inversion carriers in the general population could reach significant values. This situation could have important implications, as they have been described as predisposing haplotypes for genomic disorders. As a whole, results confirm the high instability of the 15q11-q13 region, which is prone to different types of de novo reorganizations by intrachromatid NAHR.

Human copy number variation and complex genetic disease

Copy number variants (CNVs) play an important role in human disease and population diversity. Advancements in technology have allowed for the analysis of CNVs in thousands of individuals with disease in addition to thousands of controls. These studies have identified rare CNVs associated with neuropsychiatric diseases such as autism, schizophrenia, and intellectual disability. In addition, copy number polymorphisms (CNPs) are present at higher frequencies in the population, show high diversity in copy number, sequence, and structure, and have been associated with multiple phenotypes, primarily related to immune or environmental response. However, the landscape of copy number variation still remains largely unexplored, especially for smaller CNVs and those embedded within complex regions of the human genome. An integrated approach including characterization of single nucleotide variants and CNVs in a large number of individuals with disease and normal genomes holds the promise of thoroughly elucidating the genetic basis of human disease and diversity.