Nucleotide, cytogenetic and expression impact of the human chromosome 8p23.1 inversion polymorphism

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.

Identifying novel genetic loci associated with polycystic ovary syndrome based on its shared genetic architecture with type 2 diabetes

Genome-wide association studies (GWAS) have identified several common variants associated with polycystic ovary syndrome (PCOS). However, the etiology behind PCOS remains incomplete. Available evidence suggests a potential genetic correlation between PCOS and type 2 diabetes (T2D). The publicly available data may provide an opportunity to enhance the understanding of the PCOS etiology. Here, we quantified the polygenic overlap between PCOS and T2D using summary statistics of PCOS and T2D and then identified the novel genetic variants associated with PCOS behind this phenotypic association. A bivariate causal mixture model (MiXeR model) found a moderate genetic overlap between PCOS and T2D (Dice coefficient = 44.1% and after adjusting for body mass index, 32.1%). The conditional/conjunctional false discovery rate method identified 11 potential risk variants of PCOS conditional on associations with T2D, 9 of which were novel and 6 of which were jointly associated with two phenotypes. The functional annotation of these genetic variants supports a significant role for genes involved in lipid metabolism, immune response, and the insulin signaling pathway. An expression quantitative trait locus functionality analysis successfully repeated that 5 loci were significantly associated with the expression of candidate genes in many tissues, including the whole blood, subcutaneous adipose, adrenal gland, and cerebellum. We found that SCN2A gene is co-localized with PCOS in subcutaneous adipose using GWAS-eQTL co-localization analyses. A total of 11 candidate genes were differentially expressed in multiple tissues of the PCOS samples. These findings provide a new understanding of the shared genetic architecture between PCOS and T2D and the underlying molecular genetic mechanism of PCOS.

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'.

Methionine sulfoxide reductase (Msr) dysfunction in human brain disease

Extensive research has shown that oxidative stress is strongly associated with aging, senescence and several diseases, including neurodegenerative and psychiatric disorders. Oxidative stress is caused by the overproduction of reactive oxygen species (ROS) that can be counteracted by both enzymatic and nonenzymatic antioxidants. One of these antioxidant mechanisms is the widely studied methionine sulfoxide reductase system (Msr). Methionine is one of the most easily oxidized amino acids and Msr can reverse this oxidation and restore protein function, with MsrA and MsrB reducing different stereoisomers. This article focuses on experimental and genetic research performed on Msr and its link to brain diseases. Studies on several model systems as well as genome-wide association studies are compiled to highlight the role of MSRA in schizophrenia, Alzheimer's disease, and Parkinson's disease. Genetic variation of MSRA may also contribute to the risk of psychosis, personality traits, and metabolic factors.

Multiple signals at the extended 8p23 locus are associated with susceptibility to systemic lupus erythematosus

BACKGROUND

A major systemic lupus erythematosus (SLE) susceptibility locus lies within a common inversion polymorphism region (encompassing 3.8 - 4.5  Mb) located at 8p23. Initially implicated genes included FAM167A-BLK and XKR6, of which BLK received major attention due to its known role in B-cell biology. Recently, additional SLE risk carried in non-inverted background was also reported.

OBJECTIVE AND METHODS

In this case -control study, we further investigated the 'extended' 8p23 locus (~ 4  Mb) where we observed multiple SLE signals and assessed these signals for their relation to the inversion affecting this region. The study involved a North American discovery data set (~ 1200  subjects) and a replication data set (> 10 000  subjects) comprising European-descent individuals.

RESULTS

Meta-analysis of 8p23 SNPs, with p < 0.05 in both data sets, identified 51 genome-wide significant SNPs (p < 5.0 × 10^-8). While most of these SNPs were related to previously implicated signals (XKR6-FAM167A-BLK subregion), our results also revealed two 'new' SLE signals, including SGK223-CLDN23-MFHAS1 (6.06 × 10^-9 ≤ meta p ≤ 4.88 × 10^-8) and CTSB (meta p = 4.87 × 10^-8) subregions that are located > 2 Mb upstream and ~ 0.3  Mb downstream from previously reported signals. Functional assessment of relevant SNPs indicated putative cis-effects on the expression of various genes at 8p23. Additional analyses in discovery sample, where the inversion genotypes were inferred, replicated the association of non-inverted status with SLE risk and suggested that a number of SLE risk alleles are predominantly carried in non-inverted background.

CONCLUSIONS

Our results implicate multiple (known+novel) SLE signals/genes at the extended 8p23 locus, beyond previously reported signals/genes, and suggest that this broad locus contributes to SLE risk through the effects of multiple genes/pathways.

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.

Muscle connective tissue controls development of the diaphragm and is a source of congenital diaphragmatic hernias

The diaphragm is an essential mammalian skeletal muscle, and defects in diaphragm development are the cause of congenital diaphragmatic hernias (CDH), a common and often lethal birth defect. The diaphragm is derived from multiple embryonic sources, but how these give rise to the diaphragm is unknown and, despite the identification of many CDH-associated genes, the etiology of CDH is incompletely understood. Using mouse genetics, we show that the pleuroperitoneal folds (PPFs), transient embryonic structures, are the source of the diaphragm’s muscle connective tissue, regulate muscle development, and their striking migration controls diaphragm morphogenesis. Furthermore, Gata4 mosaic mutations in PPF-derived muscle connective tissue fibroblasts result in the development of localized amuscular regions that are biomechanically weaker and more compliant and lead to CDH. Thus the PPFs and muscle connective tissue are critical for diaphragm development and mutations in PPF-derived fibroblasts are a source of CDH.

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.

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.

Nuclear factors: roles related to mitochondrial deafness

Hearing loss (HL) is a common disorder with mitochondrial dysfunction as one of the major causes leading to deafness. Mitochondrial dysfunction may be caused by either mutations in nuclear genes leading to defective nuclear-encoded proteins or mutations in mitochondrial genes leading to defective mitochondrial-encoded products. The specific nuclear genes involved in HL can be classified into two categories depending on whether mitochondrial gene mutations co-exist (modifier genes) or not (deafness-causing genes). TFB1M, MTO1, GTPBP3, and TRMU are modifier genes. A mutation in any of these modifier genes may lead to a deafness phenotype when accompanied by the mitochondrial gene mutation. OPA1, TIMM8A, SMAC/DIABLO, MPV17, PDSS1, BCS1L, SUCLA2, C10ORF2, COX10, PLOG1and RRM2B are deafness-causing genes. A mutation in any of these deafness-causing genes will directly induce variable phenotypic HL.

Congenital diaphragmatic hernia interval on chromosome 8p23.1 characterized by genetics and protein interaction networks

Chromosome 8p23.1 is a common hotspot associated with major congenital malformations, including congenital diaphragmatic hernia (CDH) and cardiac defects. We present findings from high-resolution arrays in patients who carry a loss (n = 18) or a gain (n = 1) of sub-band 8p23.1. We confirm a region involved in both diaphragmatic and heart malformations. Results from a novel CNVConnect algorithm, prioritizing protein-protein interactions between products of genes in the 8p23.1 hotspot and products of previously known CDH causing genes, implicated GATA4, NEIL2, and SOX7 in diaphragmatic defects. Sequence analysis of these genes in 226 chromosomally normal CDH patients, as well as in a small number of deletion 8p23.1 patients, showed rare unreported variants in the coding region; these may be contributing to the diaphragmatic phenotype. We also demonstrated that two of these three genes were expressed in the E11.5-12.5 primordial mouse diaphragm, the developmental stage at which CDH is thought to occur. This combination of bioinformatics and expression studies can be applied to other chromosomal hotspots, as well as private microdeletions or microduplications, to identify causative genes and their interaction networks.

An overview of the genetic structure within the Italian population from genome-wide data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F(ST) statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.

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.

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.

Identification of polymorphic inversions from genotypes

Background

Polymorphic inversions are a source of genetic variability with a direct impact on recombination frequencies. Given the difficulty of their experimental study, computational methods have been developed to infer their existence in a large number of individuals using genome-wide data of nucleotide variation. Methods based on haplotype tagging of known inversions attempt to classify individuals as having a normal or inverted allele. Other methods that measure differences between linkage disequilibrium attempt to identify regions with inversions but unable to classify subjects accurately, an essential requirement for association studies.

Results

We present a novel method to both identify polymorphic inversions from genome-wide genotype data and classify individuals as containing a normal or inverted allele. Our method, a generalization of a published method for haplotype data [1], utilizes linkage between groups of SNPs to partition a set of individuals into normal and inverted subpopulations. We employ a sliding window scan to identify regions likely to have an inversion, and accumulation of evidence from neighboring SNPs is used to accurately determine the inversion status of each subject. Further, our approach detects inversions directly from genotype data, thus increasing its usability to current genome-wide association studies (GWAS).

Conclusions

We demonstrate the accuracy of our method to detect inversions and classify individuals on principled-simulated genotypes, produced by the evolution of an inversion event within a coalescent model [2]. We applied our method to real genotype data from HapMap Phase III to characterize the inversion status of two known inversions within the regions 17q21 and 8p23 across 1184 individuals. Finally, we scan the full genomes of the European Origin (CEU) and Yoruba (YRI) HapMap samples. We find population-based evidence for 9 out of 15 well-established autosomic inversions, and for 52 regions previously predicted by independent experimental methods in ten (9+1) individuals [3,4]. We provide efficient implementations of both genotype and haplotype methods as a unified R package inveRsion.

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.