Common inversion polymorphism at 17q21.31 affects expression of multiple genes in tissue-specific manner

Unraveling the complex role of MAPT-containing H1 and H2 haplotypes in neurodegenerative diseases

A ~ 1 Mb inversion polymorphism exists within the 17q21.31 locus of the human genome as direct (H1) and inverted (H2) haplotype clades. This inversion region demonstrates high linkage disequilibrium, but the frequency of each haplotype differs across ancestries. While the H1 haplotype exists in all populations and shows a normal pattern of genetic variability and recombination, the H2 haplotype is enriched in European ancestry populations, is less frequent in African ancestry populations, and nearly absent in East Asian ancestry populations. H1 is a known risk factor for several neurodegenerative diseases, and has been associated with many other traits, suggesting its importance in cellular phenotypes of the brain and entire body. Conversely, H2 is protective for these diseases, but is associated with predisposition to recurrent microdeletion syndromes and neurodevelopmental disorders such as autism. Many single nucleotide variants and copy number variants define H1/H2 haplotypes and sub-haplotypes, but identifying the causal variant(s) for specific diseases and phenotypes is complex due to the extended linkage equilibrium. In this review, we assess the current knowledge of this inversion region regarding genomic structure, gene expression, cellular phenotypes, and disease association. We discuss recent discoveries and challenges, evaluate gaps in knowledge, and highlight the importance of understanding the effect of the 17q21.31 haplotypes to promote advances in precision medicine and drug discovery for several diseases.

MAPT H2 haplotype and risk of Pick's disease in the Pick's disease International Consortium: a genetic association study

BACKGROUND

Pick's disease is a rare and predominantly sporadic form of frontotemporal dementia that is classified as a primary tauopathy. Pick's disease is pathologically defined by the presence in the frontal and temporal lobes of Pick bodies, composed of hyperphosphorylated, three-repeat tau protein, encoded by the MAPT gene. MAPT has two distinct haplotypes, H1 and H2; the MAPT H1 haplotype is the major genetic risk factor for four-repeat tauopathies (eg, progressive supranuclear palsy and corticobasal degeneration), and the MAPT H2 haplotype is protective for these disorders. The primary aim of this study was to evaluate the association of MAPT H2 with Pick's disease risk, age at onset, and disease duration.

METHODS

In this genetic association study, we used data from the Pick's disease International Consortium, which we established to enable collection of data from individuals with pathologically confirmed Pick's disease worldwide. For this analysis, we collected brain samples from individuals with pathologically confirmed Pick's disease from 35 sites (brainbanks and hospitals) in North America, Europe, and Australia between Jan 1, 2020, and Jan 31, 2023. Neurologically healthy controls were recruited from the Mayo Clinic (FL, USA, or MN, USA between March 1, 1998, and Sept 1, 2019). For the primary analysis, individuals were directly genotyped for the MAPT H1-H2 haplotype-defining variant rs8070723. In a secondary analysis, we genotyped and constructed the six-variant-defined (rs1467967-rs242557-rs3785883-rs2471738-rs8070723-rs7521) MAPT H1 subhaplotypes. Associations of MAPT variants and MAPT haplotypes with Pick's disease risk, age at onset, and disease duration were examined using logistic and linear regression models; odds ratios (ORs) and β coefficients were estimated and correspond to each additional minor allele or each additional copy of the given haplotype.

FINDINGS

We obtained brain samples from 338 people with pathologically confirmed Pick's disease (205 [61%] male and 133 [39%] female; 338 [100%] White) and 1312 neurologically healthy controls (611 [47%] male and 701 [53%] female; 1312 [100%] White). The MAPT H2 haplotype was associated with increased risk of Pick's disease compared with the H1 haplotype (OR 1·35 [95% CI 1·12 to 1·64], p=0·0021). MAPT H2 was not associated with age at onset (β -0·54 [95% CI -1·94 to 0·87], p=0·45) or disease duration (β 0·05 [-0·06 to 0·16], p=0·35). Although not significant after correcting for multiple testing, associations were observed at p less than 0·05: with risk of Pick's disease for the H1f subhaplotype (OR 0·11 [0·01 to 0·99], p=0·049); with age at onset for H1b (β 2·66 [0·63 to 4·70], p=0·011), H1i (β -3·66 [-6·83 to -0·48], p=0·025), and H1u (β -5·25 [-10·42 to -0·07], p=0·048); and with disease duration for H1x (β -0·57 [-1·07 to -0·07], p=0·026).

INTERPRETATION

The Pick's disease International Consortium provides an opportunity to do large studies to enhance our understanding of the pathobiology of Pick's disease. This study shows that, in contrast to the decreased risk of four-repeat tauopathies, the MAPT H2 haplotype is associated with an increased risk of Pick's disease in people of European ancestry. This finding could inform development of isoform-related therapeutics for tauopathies.

FUNDING

Wellcome Trust, Rotha Abraham Trust, Brain Research UK, the Dolby Fund, Dementia Research Institute (Medical Research Council), US National Institutes of Health, and the Mayo Clinic Foundation.

Associations between polygenic risk scores and accelerated brain ageing in smokers

BACKGROUND

Smoking contributes to a variety of neurodegenerative diseases and neurobiological abnormalities, suggesting that smoking is associated with accelerated brain aging. However, the neurobiological mechanisms affected by smoking, and whether they are genetically influenced, remain to be investigated.

METHODS

Using structural magnetic resonance imaging data from the UK Biobank (n = 33 293), a brain age predictor was trained on non-smoking healthy groups and tested on smokers to obtain the BrainAge Gap (BAG). The cumulative effect of multiple common genetic variants associated with smoking was then calculated to acquire a polygenic risk score (PRS). The relationship between PRS, BAG, total gray matter volume (tGMV), and smoking parameters was explored and further genes included in the PRS were annotated to identify potential molecular mechanisms affected by smoking.

RESULTS

The BrainAge in smokers was predicted with very high accuracy (r = 0.725, MAE = 4.16). Smokers had a greater BAG (Cohen's d = 0.074, p < 0.0001) and higher PRS (Cohen's d = 0.63, p < 0.0001) than non-smokers. A higher PRS was associated with increased amount of smoking, mediated by BAG and tGMV. Several neurotransmitters and ion channel pathways were enriched in the group of smoking-related genes involved in addiction, brain synaptic plasticity, and some neurological disorders.

CONCLUSION

By using a simplified single indicator of the entire brain (BAG) in combination with the PRS, this study highlights the greater BAG in smokers and its linkage with genes and smoking behavior, providing insight into the neurobiological underpinnings and potential features of smoking-related aging.

Interindividual variation in human cortical cell type abundance and expression

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and disease.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

Mendelian Randomization and Transcriptome-Wide Association Analysis Identified Genes That Were Pleiotropically Associated with Intraocular Pressure

BACKGROUND

Intraocular pressure (IOP) is a major modifiable risk factor for glaucoma. However, the mechanisms underlying the controlling of IOP remain to be elucidated.

OBJECTIVE

To prioritize genes that are pleiotropically associated with IOP.

METHODS

We adopted a two-sample Mendelian randomization method, named summary-based Mendelian randomization (SMR), to examine the pleiotropic effect of gene expression on IOP. The SMR analyses were based on summarized data from a genome-wide association study (GWAS) on IOP. We conducted separate SMR analyses using Genotype-Tissue Expression (GTEx) and Consortium for the Architecture of Gene Expression (CAGE) expression quantitative trait loci (eQTL) data. Additionally, we performed a transcriptome-wide association study (TWAS) to identify genes whose cis-regulated expression levels were associated with IOP.

RESULTS

We identified 19 and 25 genes showing pleiotropic association with IOP using the GTEx and CAGE eQTL data, respectively. RP11-259G18.3 (P_SMR = 2.66 × 10^-6), KANSL1-AS1 (P_SMR = 2.78 × 10^-6), and RP11-259G18.2 (P_SMR = 2.91 × 10^-6) were the top three genes using the GTEx eQTL data. LRRC37A4 (P_SMR = 1.19 × 10^-5), MGC57346 (P_SMR = 1.19 × 10^-5), and RNF167 (P_SMR = 1.53 × 10^-5) were the top three genes using the CAGE eQTL data. Most of the identified genes were found in or near the 17q21.31 genomic region. Additionally, our TWAS analysis identified 18 significant genes whose expression was associated with IOP. Of these, 12 and 4 were also identified by the SMR analysis using the GTEx and CAGE eQTL data, respectively.

CONCLUSIONS

Our findings suggest that the 17q21.31 genomic region may play a critical role in the regulation of IOP.

Creating the Pick's disease International Consortium: Association study of MAPT H2 haplotype with risk of Pick's disease

Background

Pick's disease (PiD) is a rare and predominantly sporadic form of frontotemporal dementia that is classified as a primary tauopathy. PiD is pathologically defined by argyrophilic inclusion Pick bodies and ballooned neurons in the frontal and temporal brain lobes. PiD is characterised by the presence of Pick bodies which are formed from aggregated, hyperphosphorylated, 3-repeat tau proteins, encoded by the MAPT gene. The MAPT H2 haplotype has consistently been associated with a decreased disease risk of the 4-repeat tauopathies of progressive supranuclear palsy and corticobasal degeneration, however its role in susceptibility to PiD is unclear. The primary aim of this study was to evaluate the association between MAPT H2 and risk of PiD.

Methods

We established the Pick's disease International Consortium (PIC) and collected 338 (60.7% male) pathologically confirmed PiD brains from 39 sites worldwide. 1,312 neurologically healthy clinical controls were recruited from Mayo Clinic Jacksonville, FL (N=881) or Rochester, MN (N=431). For the primary analysis, subjects were directly genotyped for MAPT H1-H2 haplotype-defining variant rs8070723. In secondary analysis, we genotyped and constructed the six-variant MAPT H1 subhaplotypes (rs1467967, rs242557, rs3785883, rs2471738, rs8070723, and rs7521).

Findings

Our primary analysis found that the MAPT H2 haplotype was associated with increased risk of PiD (OR: 1.35, 95% CI: 1.12-1.64 P=0.002). In secondary analysis involving H1 subhaplotypes, a protective association with PiD was observed for the H1f haplotype (0.0% vs. 1.2%, P=0.049), with a similar trend noted for H1b (OR: 0.76, 95% CI: 0.58-1.00, P=0.051). The 4-repeat tauopathy risk haplotype MAPT H1c was not associated with PiD susceptibility (OR: 0.93, 95% CI: 0.70-1.25, P=0.65).

Interpretation

The PIC represents the first opportunity to perform relatively large-scale studies to enhance our understanding of the pathobiology of PiD. This study demonstrates that in contrast to its protective role in 4R tauopathies, the MAPT H2 haplotype is associated with an increased risk of PiD. This finding is critical in directing isoform-related therapeutics for tauopathies.

Different MAPT haplotypes influence expression of total MAPT in postmortem brain tissue

The MAPT gene, encoding the microtubule-associated protein tau on chromosome 17q21.31, is result of an inversion polymorphism, leading to two allelic variants (H1 and H2). Homozygosity for the more common haplotype H1 is associated with an increased risk for several tauopathies, but also for the synucleinopathy Parkinson's disease (PD). In the present study, we aimed to clarify whether the MAPT haplotype influences expression of MAPT and SNCA, encoding the protein α-synuclein (α-syn), on mRNA and protein levels in postmortem brains of PD patients and controls. We also investigated mRNA expression of several other MAPT haplotype-encoded genes. Postmortem tissues from cortex of fusiform gyrus (ctx-fg) and of the cerebellar hemisphere (ctx-cbl) of neuropathologically confirmed PD patients (n = 95) and age- and sex-matched controls (n = 81) were MAPT haplotype genotyped to identify cases homozygous for either H1 or H2. Relative expression of genes was quantified using real-time qPCR; soluble and insoluble protein levels of tau and α-syn were determined by Western blotting. Homozygosity for H1 versus H2 was associated with increased total MAPT mRNA expression in ctx-fg regardless of disease state. Inversely, H2 homozygosity was associated with markedly increased expression of the corresponding antisense MAPT-AS1 in ctx-cbl. PD patients had higher levels of insoluble 0N3R and 1N4R tau isoforms regardless of the MAPT genotype. The increased presence of insoluble α-syn in PD patients in ctx-fg validated the selected postmortem brain tissue. Our findings in this small, but well controlled cohort of PD and controls support a putative biological relevance of tau in PD. However, we did not identify any link between the disease-predisposing H1/H1 associated overexpression of MAPT with PD status. Further studies are required to gain a deeper understanding of the potential regulatory role of MAPT-AS1 and its association to the disease-protective H2/H2 condition in the context of PD.

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.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci

BACKGROUND

Known risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort.

METHODS

Single nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer-related cell types.

RESULTS

We identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types.

CONCLUSIONS

CNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

TWAS Atlas: a curated knowledgebase of transcriptome-wide association studies

Transcriptome-wide association studies (TWASs), as a practical and prevalent approach for detecting the associations between genetically regulated genes and traits, are now leading to a better understanding of the complex mechanisms of genetic variants in regulating various diseases and traits. Despite the ever-increasing TWAS outputs, there is still a lack of databases curating massive public TWAS information and knowledge. To fill this gap, here we present TWAS Atlas (https://ngdc.cncb.ac.cn/twas/), an integrated knowledgebase of TWAS findings manually curated from extensive literature. In the current implementation, TWAS Atlas collects 401,266 high-quality human gene-trait associations from 200 publications, covering 22,247 genes and 257 traits across 135 tissue types. In particular, an interactive knowledge graph of the collected gene-trait associations is constructed together with single nucleotide polymorphism (SNP)-gene associations to build up comprehensive regulatory networks at multi-omics levels. In addition, TWAS Atlas, as a user-friendly web interface, efficiently enables users to browse, search and download all association information, relevant research metadata and annotation information of interest. Taken together, TWAS Atlas is of great value for promoting the utility and availability of TWAS results in explaining the complex genetic basis as well as providing new insights for human health and disease research.

Association Between Brain Structure and Alcohol Use Behaviors in Adults: A Mendelian Randomization and Multiomics Study

IMPORTANCE

Past studies have identified associations between brain macrostructure and alcohol use behaviors. However, identifying directional associations between these phenotypes is difficult due to the limitations of observational studies.

OBJECTIVE

To use mendelian randomization (MR) to identify directional associations between brain structure and alcohol use and elucidate the transcriptomic and cellular underpinnings of identified associations.

DESIGN, SETTING, AND PARTICIPANTS

The main source data comprised summary statistics from population-based and case-control genome-wide association studies (GWAS) of neuroimaging, behavioral, and clinical phenotypes (N = 763 874). Using these data, bidirectional and multivariable MR was performed analyzing associations between brain macrostructure and alcohol use. Downstream transcriptome-wide association studies (TWAS) and cell-type enrichment analyses investigated the biology underlying identified associations. The study approach was data driven and did not test any a priori hypotheses. Data were analyzed August 2021 to May 2022.

MAIN OUTCOMES AND MEASURES

Brain structure phenotypes (global cortical thickness [GCT] and global cortical surface area [GCSA] in 33 709 individuals and left-right subcortical volumes in 19 629 individuals) and alcohol use behaviors (alcoholic drinks per week [DPW] in 537 349 individuals, binge drinking frequency in 143 685 individuals, and alcohol use disorder in 8845 individuals vs 20 657 control individuals [total of 29 502]).

RESULTS

The main bidirectional MR analyses were performed in samples totaling 763 874 individuals, among whom more than 94% were of European ancestry, 52% to 54% were female, and the mean cohort ages were 40 to 63 years. Negative associations were identified between genetically predicted GCT and binge drinking (β, -2.52; 95% CI, -4.13 to -0.91) and DPW (β, -0.88; 95% CI, -1.37 to -0.40) at a false discovery rate (FDR) of 0.05. These associations remained significant in multivariable MR models that accounted for neuropsychiatric phenotypes, substance use, trauma, and neurodegeneration. TWAS of GCT and alcohol use behaviors identified 5 genes at the 17q21.31 locus oppositely associated with GCT and binge drinking or DPW (FDR = 0.05). Cell-type enrichment analyses implicated glutamatergic cortical neurons in alcohol use behaviors.

CONCLUSIONS AND RELEVANCE

The findings in this study show that the associations between GCT and alcohol use may reflect a predispositional influence of GCT and that 17q21.31 genes and glutamatergic cortical neurons may play a role in this association. While replication studies are needed, these findings should enhance the understanding of associations between brain structure and alcohol use.

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

The early-life exposome modulates the effect of polymorphic inversions on DNA methylation

Polymorphic genomic inversions are chromosomal variants with intrinsic variability that play important roles in evolution, environmental adaptation, and complex traits. We investigated the DNA methylation patterns of three common human inversions, at 8p23.1, 16p11.2, and 17q21.31 in 1,009 blood samples from children from the Human Early Life Exposome (HELIX) project and in 39 prenatal heart tissue samples. We found inversion-state specific methylation patterns within and nearby flanking each inversion region in both datasets. Additionally, numerous inversion-exposure interactions on methylation levels were identified from early-life exposome data comprising 64 exposures. For instance, children homozygous at inv-8p23.1 and higher meat intake were more susceptible to TDH hypermethylation (P = 3.8 × 10⁻²²); being the inversion, exposure, and gene known risk factors for adult obesity. Inv-8p23.1 associated hypermethylation of GATA4 was also detected across numerous exposures. Our data suggests that the pleiotropic influence of inversions during development and lifetime could be substantially mediated by allele-specific methylation patterns which can be modulated by the exposome.

Analysis of the relationship between presence of common DNA sequence inversions and DNA methylation patterns suggests a role for environmental exposures (such as food intake) in mediating inversion state-specific methylation patterns.

Discovery of genomic loci of the human cerebral cortex using genetically informed brain atlases

To determine the impact of genetic variants on the brain, we used genetically informed brain atlases in genome-wide association studies of regional cortical surface area and thickness in 39,898 adults and 9136 children. We uncovered 440 genome-wide significant loci in the discovery cohort and 800 from a post hoc combined meta-analysis. Loci in adulthood were largely captured in childhood, showing signatures of negative selection, and were linked to early neurodevelopment and pathways associated with neuropsychiatric risk. Opposing gradations of decreased surface area and increased thickness were associated with common inversion polymorphisms. Inferior frontal regions, encompassing Broca's area, which is important for speech, were enriched for human-specific genomic elements. Thus, a mixed genetic landscape of conserved and human-specific features is concordant with brain hierarchy and morphogenetic gradients.

Dysregulated coordination of MAPT exon 2 and exon 10 splicing underlies different tau pathologies in PSP and AD

Understanding regulation of MAPT splicing is important to the etiology of many nerurodegenerative diseases, including Alzheimer disease (AD) and progressive supranuclear palsy (PSP), in which different tau isoforms accumulate in pathologic inclusions. MAPT, the gene encoding the tau protein, undergoes complex alternative pre-mRNA splicing to generate six isoforms. Tauopathies can be categorized by the presence of tau aggregates containing either 3 (3R) or 4 (4R) microtubule-binding domain repeats (determined by inclusion/exclusion of exon 10), but the role of the N-terminal domain of the protein, determined by inclusion/exclusion of exons 2 and 3 has been less well studied. Using a correlational screen in human brain tissue, we observed coordination of MAPT exons 2 and 10 splicing. Expressions of exon 2 splicing regulators and subsequently exon 2 inclusion are differentially disrupted in PSP and AD brain, resulting in the accumulation of 1N4R isoforms in PSP and 0N isoforms in AD temporal cortex. Furthermore, we identified different N-terminal isoforms of tau present in neurofibrillary tangles, dystrophic neurites and tufted astrocytes, indicating a role for differential N-terminal splicing in the development of disparate tau neuropathologies. We conclude that N-terminal splicing and combinatorial regulation with exon 10 inclusion/exclusion is likely to be important to our understanding of tauopathies.

TIGAR-V2: Efficient TWAS tool with nonparametric Bayesian eQTL weights of 49 tissue types from GTEx V8

Standard transcriptome-wide association study (TWAS) methods first train gene expression prediction models using reference transcriptomic data and then test the association between the predicted genetically regulated gene expression and phenotype of interest. Most existing TWAS tools require cumbersome preparation of genotype input files and extra coding to enable parallel computation. To improve the efficiency of TWAS tools, we developed Transcriptome-Integrated Genetic Association Resource V2 (TIGAR-V2), which directly reads Variant Call Format (VCF) files, enables parallel computation, and reduces up to 90% of computation cost (mainly due to loading genotype data) compared to the original version. TIGAR-V2 can train gene expression imputation models using either nonparametric Bayesian Dirichlet process regression (DPR) or Elastic-Net (as used by PrediXcan), perform TWASs using either individual-level or summary-level genome-wide association study (GWAS) data, and implement both burden and variance-component statistics for gene-based association tests. We trained gene expression prediction models by DPR for 49 tissues using Genotype-Tissue Expression (GTEx) V8 by TIGAR-V2 and illustrated the usefulness of these Bayesian cis-expression quantitative trait locus (eQTL) weights through TWASs of breast and ovarian cancer utilizing public GWAS summary statistics. We identified 88 and 37 risk genes, respectively, for breast and ovarian cancer, most of which are either known or near previously identified GWAS (∼95%) or TWAS (∼40%) risk genes and three novel independent TWAS risk genes with known functions in carcinogenesis. These findings suggest that TWASs can provide biological insight into the transcriptional regulation of complex diseases. The TIGAR-V2 tool, trained Bayesian cis-eQTL weights, and linkage disequilibrium (LD) information from GTEx V8 are publicly available, providing a useful resource for mapping risk genes of complex diseases.

iPS Cell-Based Model for MAPT Haplotype as a Risk Factor for Human Tauopathies Identifies No Major Differences in TAU Expression

The H1 haplotype of the microtubule-associated protein tau (MAPT) gene is a common genetic risk factor for some neurodegenerative diseases such as progressive supranuclear palsy, corticobasal degeneration, and Parkinson's disease. The molecular mechanism causing the increased risk for the named diseases, however, remains unclear. In this paper, we present a valuable tool of eight small molecule neural precursor cell lines (smNPC) homozygous for the MAPT haplotypes (four H1/H1 and four H2/H2 cell lines), which can be used to identify MAPT-dependent phenotypes. The employed differentiation protocol is fast due to overexpression of NEUROGENIN-2 and therefore suitable for high-throughput approaches. A basic characterization of all human cell lines was performed, and their TAU and α-SYNUCLEIN profiles were compared during a differentiation time of 30 days. We could identify higher levels of conformationally altered TAU in cell lines carrying the H2 haplotype. Additionally, we found increased expression levels of α-SYNUCLEIN in H1/H1 cells. With this resource, we aim to fill a gap in neurodegenerative disease modeling with induced pluripotent stem cells (iPSC) for sporadic tauopathies.

Pleiotropy-guided transcriptome imputation from normal and tumor tissues identifies candidate susceptibility genes for breast and ovarian cancer

Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Gene Expression Imputation Across Multiple Tissue Types Provides Insight Into the Genetic Architecture of Frontotemporal Dementia and Its Clinical Subtypes

BACKGROUND

The etiology of frontotemporal dementia (FTD) is poorly understood. To identify genes with predicted expression levels associated with FTD, we integrated summary statistics with external reference gene expression data using a transcriptome-wide association study approach.

METHODS

FUSION software was used to leverage FTD summary statistics (all FTD: n = 2154 cases, n = 4308 controls; behavioral variant FTD: n = 1337 cases, n = 2754 controls; semantic dementia: n = 308 cases, n = 616 controls; progressive nonfluent aphasia: n = 269 cases, n = 538 controls; FTD with motor neuron disease: n = 200 cases, n = 400 controls) from the International FTD-Genomics Consortium with 53 expression quantitative loci tissue type panels (n = 12,205; 5 consortia). Significance was assessed using a 5% false discovery rate threshold.

RESULTS

We identified 73 significant gene-tissue associations for FTD, representing 44 unique genes in 34 tissue types. Most significant findings were derived from dorsolateral prefrontal cortex splicing data (n = 19 genes, 26%). The 17q21.31 inversion locus contained 23 significant associations, representing 6 unique genes. Other top hits included SEC22B (a gene involved in vesicle trafficking), TRGV5, and ZNF302. A single gene finding (RAB38) was observed for behavioral variant FTD. For other clinical subtypes, no significant associations were observed.

CONCLUSIONS

We identified novel candidate genes (e.g., SEC22B) and previously reported risk regions (e.g., 17q21.31) for FTD. Most significant associations were observed in dorsolateral prefrontal cortex splicing data despite the modest sample size of this reference panel. This suggests that our findings are specific to FTD and are likely to be biologically relevant highlights of genes at different FTD risk loci that are contributing to the disease pathology.

An integrative systems-based analysis of substance use: eQTL-informed gene-based tests, gene networks, and biological mechanisms

Genome-wide association studies have identified multiple genetic risk factors underlying susceptibility to substance use, however, the functional genes and biological mechanisms remain poorly understood. The discovery and characterization of risk genes can be facilitated by the integration of genome-wide association data and gene expression data across biologically relevant tissues and/or cell types to identify genes whose expression is altered by DNA sequence variation (expression quantitative trait loci; eQTLs). The integration of gene expression data can be extended to the study of genetic co-expression, under the biologically valid assumption that genes form co-expression networks to influence the manifestation of a disease or trait. Here, we integrate genome-wide association data with gene expression data from 13 brain tissues to identify candidate risk genes for 8 substance use phenotypes. We then test for the enrichment of candidate risk genes within tissue-specific gene co-expression networks to identify modules (or groups) of functionally related genes whose dysregulation is associated with variation in substance use. We identified eight gene modules in brain that were enriched with gene-based association signals for substance use phenotypes. For example, a single module of 40 co-expressed genes was enriched with gene-based associations for drinks per week and biological pathways involved in GABA synthesis, release, reuptake and degradation. Our study demonstrates the utility of eQTL and gene co-expression analysis to uncover novel biological mechanisms for substance use traits.

Genome-wide association study identifies 48 common genetic variants associated with handedness

Handedness has been extensively studied because of its relationship with language and the over-representation of left-handers in some neurodevelopmental disorders. Using data from the UK Biobank, 23andMe and the International Handedness Consortium, we conducted a genome-wide association meta-analysis of handedness (N = 1,766,671). We found 41 loci associated (P < 5 × 10-8) with left-handedness and 7 associated with ambidexterity. Tissue-enrichment analysis implicated the CNS in the aetiology of handedness. Pathways including regulation of microtubules and brain morphology were also highlighted. We found suggestive positive genetic correlations between left-handedness and neuropsychiatric traits, including schizophrenia and bipolar disorder. Furthermore, the genetic correlation between left-handedness and ambidexterity is low (rG = 0.26), which implies that these traits are largely influenced by different genetic mechanisms. Our findings suggest that handedness is highly polygenic and that the genetic variants that predispose to left-handedness may underlie part of the association with some psychiatric disorders.

The genomic basis of mood instability: identification of 46 loci in 363,705 UK Biobank participants, genetic correlation with psychiatric disorders, and association with gene expression and function

Genome-wide association studies (GWAS) of psychiatric phenotypes have tended to focus on categorical diagnoses, but to understand the biology of mental illness it may be more useful to study traits which cut across traditional boundaries. Here, we report the results of a GWAS of mood instability as a trait in a large population cohort (UK Biobank, n = 363,705). We also assess the clinical and biological relevance of the findings, including whether genetic associations show enrichment for nervous system pathways. Forty six unique loci associated with mood instability were identified with a SNP heritability estimate of 9%. Linkage Disequilibrium Score Regression (LDSR) analyses identified genetic correlations with Major Depressive Disorder (MDD), Bipolar Disorder (BD), Schizophrenia, anxiety, and Post Traumatic Stress Disorder (PTSD). Gene-level and gene set analyses identified 244 significant genes and 6 enriched gene sets. Tissue expression analysis of the SNP-level data found enrichment in multiple brain regions, and eQTL analyses highlighted an inversion on chromosome 17 plus two brain-specific eQTLs. In addition, we used a Phenotype Linkage Network (PLN) analysis and community analysis to assess for enrichment of nervous system gene sets using mouse orthologue databases. The PLN analysis found enrichment in nervous system PLNs for a community containing serotonin and melatonin receptors. In summary, this work has identified novel loci, tissues and gene sets contributing to mood instability. These findings may be relevant for the identification of novel trans-diagnostic drug targets and could help to inform future stratified medicine innovations in mental health.

Polymorphic Inversions Underlie the Shared Genetic Susceptibility of Obesity-Related Diseases

The burden of several common diseases including obesity, diabetes, hypertension, asthma, and depression is increasing in most world populations. However, the mechanisms underlying the numerous epidemiological and genetic correlations among these disorders remain largely unknown. We investigated whether common polymorphic inversions underlie the shared genetic influence of these disorders. We performed an inversion association analysis including 21 inversions and 25 obesity-related traits on a total of 408,898 Europeans and validated the results in 67,299 independent individuals. Seven inversions were associated with multiple diseases while inversions at 8p23.1, 16p11.2, and 11q13.2 were strongly associated with the co-occurrence of obesity with other common diseases. Transcriptome analysis across numerous tissues revealed strong candidate genes for obesity-related traits. Analyses in human pancreatic islets indicated the potential mechanism of inversions in the susceptibility of diabetes by disrupting the cis-regulatory effect of SNPs from their target genes. Our data underscore the role of inversions as major genetic contributors to the joint susceptibility to common complex diseases.

Determining the impact of uncharacterized inversions in the human genome by droplet digital PCR

Despite the interest in characterizing genomic variation, the presence of large repeats at the breakpoints hinders the analysis of many structural variants. This is especially problematic for inversions, since there is typically no gain or loss of DNA. Here, we tested novel linkage-based droplet digital PCR (ddPCR) assays to study 20 inversions ranging from 3.1 to 742 kb flanked by inverted repeats (IRs) up to 134 kb long. Of those, we validated 13 inversions predicted by different genome-wide techniques. In addition, we obtained new experimental human population information across 95 African, European, and East Asian individuals for 16 inversions, including four already validated variants without high-throughput genotyping methods. Through comparison with previous data, independent replicates and both inversion breakpoints, we demonstrate that the technique is highly accurate and reproducible. Most studied inversions are widespread across continents, and their frequency is negatively correlated with genetic length. Moreover, all except two show clear signs of being recurrent, and we could better define the factors affecting recurrence levels and estimate the inversion rate across the genome. Finally, the generated genotypes have allowed us to check inversion functional effects, validating gene expression differences reported before for two inversions and finding new candidate associations. Therefore, the developed methodology makes it possible to screen these and other complex genomic variants quickly in a large number of samples for the first time, highlighting the importance of direct genotyping to assess their potential consequences and clinical implications.

Common polymorphic inversions at 17q21.31 and 8p23.1 associate with cancer prognosis

BACKGROUND

Chromosomal inversions are structural genetic variants where a chromosome segment changes its orientation. While sporadic de novo inversions are known genetic risk factors for cancer susceptibility, it is unknown if common polymorphic inversions are also associated with the prognosis of common tumors, as they have been linked to other complex diseases. We studied the association of two well-characterized human inversions at 17q21.31 and 8p23.1 with the prognosis of lung, liver, breast, colorectal, and stomach cancers.

RESULTS

Using data from The Cancer Genome Atlas (TCGA), we observed that inv8p23.1 was associated with overall survival in breast cancer and that inv17q21.31 was associated with overall survival in stomach cancer. In the meta-analysis of two independent studies, inv17q21.31 heterozygosity was significantly associated with colorectal disease-free survival. We found that the association was mediated by the de-methylation of cg08283464 and cg03999934, also linked to lower disease-free survival.

CONCLUSIONS

Our results suggest that chromosomal inversions are important genetic factors of tumor prognosis, likely affecting changes in methylation patterns.

Language impairment in primary progressive aphasia and other neurodegenerative diseases

Primary progressive aphasia (PPA) is a progressive neurodegenerative disease that disrupts the language capacity of an individual by selectively affecting the language network of brain. Although aphasic literature is replete with reports of brain damage responsible for various types of PPA, it does not provide a comprehensive understanding of whether PPA is an independent pathological condition or an atypical syndrome of neurodegenerative diseases (NDD). To address this ambiguity, we provide a detailed description of PPA, its variants and their brain anatomy. Subsequently, we unravel the relationship between PPA and NDDs like Alzheimer's, Parkinson's and Dyslexia. To substantiate the relationship further, we also provide a brief account of their genetic aetiology. In the final section, we offer an exhaustive approach towards the treatment of PPA by combining the existing language the rapies with clinical and pharmacological interventions.

Evolutionary and functional impact of common polymorphic inversions in the human genome

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1–415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.

Inversions are a little-studied type of genomic variation that could contribute to phenotypic traits. Here the authors characterize 45 common polymorphic inversions in human populations and investigate their evolutionary and functional impact.

The intramembrane protease SPPL2c promotes male germ cell development by cleaving phospholamban

Signal peptide peptidase (SPP) and the four homologous SPP-like (SPPL) proteases constitute a family of intramembrane aspartyl proteases with selectivity for type II-oriented transmembrane segments. Here, we analyse the physiological function of the orphan protease SPPL2c, previously considered to represent a non-expressed pseudogene. We demonstrate proteolytic activity of SPPL2c towards selected tail-anchored proteins. Despite shared ER localisation, SPPL2c and SPP exhibit distinct, though partially overlapping substrate spectra and inhibitory profiles, and are organised in different high molecular weight complexes. Interestingly, SPPL2c is specifically expressed in murine and human testis where it is primarily localised in spermatids. In mice, SPPL2c deficiency leads to a partial loss of elongated spermatids and reduced motility of mature spermatozoa, but preserved fertility. However, matings of male and female SPPL2c ^-/- mice exhibit reduced litter sizes. Using proteomics we identify the sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA2)-regulating protein phospholamban (PLN) as a physiological SPPL2c substrate. Accumulation of PLN correlates with a decrease in intracellular Ca²⁺ levels in elongated spermatids that likely contribute to the compromised male germ cell differentiation and function of SPPL2c ^-/- mice.

Expression quantitative trait loci in the developing human brain and their enrichment in neuropsychiatric disorders

BACKGROUND

Genetic influences on gene expression in the human fetal brain plausibly impact upon a variety of postnatal brain-related traits, including susceptibility to neuropsychiatric disorders. However, to date, there have been no studies that have mapped genome-wide expression quantitative trait loci (eQTL) specifically in the human prenatal brain.

RESULTS

We performed deep RNA sequencing and genome-wide genotyping on a unique collection of 120 human brains from the second trimester of gestation to provide the first eQTL dataset derived exclusively from the human fetal brain. We identify high confidence cis-acting eQTL at the individual transcript as well as whole gene level, including many mapping to a common inversion polymorphism on chromosome 17q21. Fetal brain eQTL are enriched among risk variants for postnatal conditions including attention deficit hyperactivity disorder, schizophrenia, and bipolar disorder. We further identify changes in gene expression within the prenatal brain that potentially mediate risk for neuropsychiatric traits, including increased expression of C4A in association with genetic risk for schizophrenia, increased expression of LRRC57 in association with genetic risk for bipolar disorder, and altered expression of multiple genes within the chromosome 17q21 inversion in association with variants influencing the personality trait of neuroticism.

CONCLUSIONS

We have mapped eQTL operating in the human fetal brain, providing evidence that these confer risk to certain neuropsychiatric disorders, and identifying gene expression changes that potentially mediate susceptibility to these conditions.

Early Life Stress and Epigenetics in Late-onset Alzheimer's Dementia: A Systematic Review

Involvement of life stress in Late-Onset Alzheimer's Disease (LOAD) has been evinced in longitudinal cohort epidemiological studies, and endocrinologic evidence suggests involvements of catecholamine and corticosteroid systems in LOAD. Early Life Stress (ELS) rodent models have successfully demonstrated sequelae of maternal separation resulting in LOAD-analogous pathology, thereby supporting a role of insulin receptor signalling pertaining to GSK-3beta facilitated tau hyper-phosphorylation and amyloidogenic processing. Discussed are relevant ELS studies, and findings from three mitogen-activated protein kinase pathways (JNK/SAPK pathway, ERK pathway, p38/MAPK pathway) relevant for mediating environmental stresses. Further considered were the roles of autophagy impairment, neuroinflammation, and brain insulin resistance. For the meta-analytic evaluation, 224 candidate gene loci were extracted from reviews of animal studies of LOAD pathophysiological mechanisms, of which 60 had no positive results in human LOAD association studies. These loci were combined with 89 gene loci confirmed as LOAD risk genes in previous GWAS and WES. Of the 313 risk gene loci evaluated, there were 35 human reports on epigenomic modifications in terms of methylation or histone acetylation. 64 microRNA gene regulation mechanisms were published for the compiled loci. Genomic association studies support close relations of both noradrenergic and glucocorticoid systems with LOAD. For HPA involvement, a CRHR1 haplotype with MAPT was described, but further association of only HSD11B1 with LOAD found; however, association of FKBP1 and NC3R1 polymorphisms was documented in support of stress influence to LOAD. In the brain insulin system, IGF2R, INSR, INSRR, and plasticity regulator ARC, were associated with LOAD. Pertaining to compromised myelin stability in LOAD, relevant associations were found for BIN1, RELN, SORL1, SORCS1, CNP, MAG, and MOG. Regarding epigenetic modifications, both methylation variability and de-acetylation were reported for LOAD. The majority of up-to-date epigenomic findings include reported modifications in the well-known LOAD core pathology loci MAPT, BACE1, APP (with FOS, EGR1), PSEN1, PSEN2, and highlight a central role of BDNF. Pertaining to ELS, relevant loci are FKBP5, EGR1, GSK3B; critical roles of inflammation are indicated by CRP, TNFA, NFKB1 modifications; for cholesterol biosynthesis, DHCR24; for myelin stability BIN1, SORL1, CNP; pertaining to (epi)genetic mechanisms, hTERT, MBD2, DNMT1, MTHFR2. Findings on gene regulation were accumulated for BACE1, MAPK signalling, TLR4, BDNF, insulin signalling, with most reports for miR-132 and miR-27. Unclear in epigenomic studies remains the role of noradrenergic signalling, previously demonstrated by neuropathological findings of childhood nucleus caeruleus degeneration for LOAD tauopathy.

A Transcriptome-Wide Association Study Among 97,898 Women to Identify Candidate Susceptibility Genes for Epithelial Ovarian Cancer Risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10-6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.

Nested Inversion Polymorphisms Predispose Chromosome 22q11.2 to Meiotic Rearrangements

Inversion polymorphisms between low-copy repeats (LCRs) might predispose chromosomes to meiotic non-allelic homologous recombination (NAHR) events and thus lead to genomic disorders. However, for the 22q11.2 deletion syndrome (22q11.2DS), the most common genomic disorder, no such inversions have been uncovered as of yet. Using fiber-FISH, we demonstrate that parents transmitting the de novo 3 Mb LCR22A-D 22q11.2 deletion, the reciprocal duplication, and the smaller 1.5 Mb LCR22A-B 22q11.2 deletion carry inversions of LCR22B-D or LCR22C-D. Hence, the inversions predispose chromosome 22q11.2 to meiotic rearrangements and increase the individual risk for transmitting rearrangements. Interestingly, the inversions are nested or flanking rather than coinciding with the deletion or duplication sizes. This finding raises the possibility that inversions are a prerequisite not only for 22q11.2 rearrangements but also for all NAHR-mediated genomic disorders.

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.

Fitness consequences of polymorphic inversions in the zebra finch genome

Background

Inversion polymorphisms constitute an evolutionary puzzle: they should increase embryo mortality in heterokaryotypic individuals but still they are widespread in some taxa. Some insect species have evolved mechanisms to reduce the cost of embryo mortality but humans have not. In birds, a detailed analysis is missing although intraspecific inversion polymorphisms are regarded as common. In Australian zebra finches (Taeniopygia guttata), two polymorphic inversions are known cytogenetically and we set out to detect these two and potentially additional inversions using genomic tools and study their effects on embryo mortality and other fitness-related and morphological traits.

Results

Using whole-genome SNP data, we screened 948 wild zebra finches for polymorphic inversions and describe four large (12–63 Mb) intraspecific inversion polymorphisms with allele frequencies close to 50 %. Using additional data from 5229 birds and 9764 eggs from wild and three captive zebra finch populations, we show that only the largest inversions increase embryo mortality in heterokaryotypic males, with surprisingly small effect sizes. We test for a heterozygote advantage on other fitness components but find no evidence for heterosis for any of the inversions. Yet, we find strong additive effects on several morphological traits.

Conclusions

The mechanism that has carried the derived inversion haplotypes to such high allele frequencies remains elusive. It appears that selection has effectively minimized the costs associated with inversions in zebra finches. The highly skewed distribution of recombination events towards the chromosome ends in zebra finches and other estrildid species may function to minimize crossovers in the inverted regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1056-3) contains supplementary material, which is available to authorized users.

Gene expression, methylation and neuropathology correlations at progressive supranuclear palsy risk loci

To determine the effects of single nucleotide polymorphisms (SNPs) identified in a genome-wide association study of progressive supranuclear palsy (PSP), we tested their association with brain gene expression, CpG methylation and neuropathology. In 175 autopsied PSP subjects, we performed associations between seven PSP risk variants and temporal cortex levels of 20 genes in-cis, within ±100 kb. Methylation measures were collected using reduced representation bisulfite sequencing in 43 PSP brains. To determine whether SNP/expression associations are due to epigenetic modifications, CpG methylation levels of associated genes were tested against relevant variants. Quantitative neuropathology endophenotypes were tested for SNP associations in 422 PSP subjects. Brain levels of LRRC37A4 and ARL17B were associated with rs8070723; MOBP with rs1768208 and both ARL17A and ARL17B with rs242557. Expression associations for LRRC37A4 and MOBP were available in an additional 100 PSP subjects. Meta-analysis revealed highly significant associations for PSP risk alleles of rs8070723 and rs1768208 with higher LRRC37A4 and MOBP brain levels, respectively. Methylation levels of one CpG in the 3' region of ARL17B associated with rs242557 and rs8070723. Additionally, methylation levels of an intronic ARL17A CpG associated with rs242557 and that of an intronic MOBP CpG with rs1768208. MAPT and MOBP region risk alleles also associated with higher levels of neuropathology. Strongest associations were observed for rs242557/coiled bodies and tufted astrocytes; and for rs1768208/coiled bodies and tau threads. These findings suggest that PSP variants at MAPT and MOBP loci may confer PSP risk via influencing gene expression and tau neuropathology. MOBP, LRRC37A4, ARL17A and ARL17B warrant further assessment as candidate PSP risk genes. Our findings have implications for the mechanism of action of variants at some of the top PSP risk loci.

MAPT haplotype H1G is associated with increased risk of dementia with Lewy bodies

INTRODUCTION

The MAPT H1 haplotype has been associated with several neurodegenerative diseases. We were interested in exploring the role of MAPT haplotypic variation in risk of dementia with Lewy bodies (DLB).

METHOD

We genotyped six MAPT haplotype tagging SNPs and screened 431 clinical DLB cases, 347 pathologically defined high-likelihood DLB cases, and 1049 controls.

RESULT

We performed haplotypic association tests and detected an association with the protective H2 haplotype in our combined series (odds ratio [OR] = 0.75). We fine-mapped the locus and identified a relatively rare haplotype, H1G, that is associated with an increased risk of DLB (OR = 3.30, P = .0017). This association was replicated in our pathologically defined series (OR = 2.26, P = .035).

DISCUSSION

These results support a role for H1 and specifically H1G in susceptibility to DLB. However, the exact functional variant at the locus is still unknown, and additional studies are warranted to fully explain genetic risk of DLB at the MAPT locus.

Population Genomic Analysis Reveals Differential Evolutionary Histories and Patterns of Diversity across Subgenomes and Subpopulations of Brassica napus L

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

Widespread Impact of Chromosomal Inversions on Gene Expression Uncovers Robustness via Phenotypic Buffering

The nonrandom gene organization in eukaryotes plays a significant role in genome evolution and function. Chromosomal structural changes impact meiotic fitness and, in several organisms, are associated with speciation and rapid adaptation to different environments. Small sized chromosomal inversions, encompassing few genes, are pervasive in Saccharomyces “sensu stricto” species, while larger inversions are less common in yeasts compared with higher eukaryotes. To explore the effect of gene order on phenotype, reproductive isolation, and gene expression, we engineered 16 Saccharomyces cerevisiae strains carrying all possible paracentric and pericentric inversions between Ty1 elements, a natural substrate for rearrangements. We found that 4 inversions were lethal, while the other 12 did not show any fitness advantage or disadvantage in rich and minimal media. At meiosis, only a weak negative correlation with fitness was seen with the size of the inverted region. However, significantly lower fertility was seen in heterozygote invertant strains carrying recombination hotspots within the breakpoints. Altered transcription was observed throughout the genome rather than being overrepresented within the inversions. In spite of the large difference in gene expression in the inverted strains, mitotic fitness was not impaired in the majority of the 94 conditions tested, indicating that the robustness of the expression network buffers the deleterious effects of structural changes in several environments. Overall, our results support the notion that transcriptional changes may compensate for Ty-mediated rearrangements resulting in the maintenance of a constant phenotype, and suggest that large inversions in yeast are unlikely to be a selectable trait during vegetative growth.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000–50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

Since the discovery of chromosomal inversions almost 100 years ago, how they are maintained in natural populations has been a highly debated issue. One of the hypotheses is that inversion breakpoints could affect genes and modify gene expression levels, although evidence of this came only from laboratory mutants. In humans, a few inversions have been shown to associate with expression differences, but in all cases the molecular causes have remained elusive. Here, we have carried out a complete characterization of a new human polymorphic inversion and determined that it is specific to East Asian populations. In addition, we demonstrate that it disrupts the ZNF257 gene and, through the translocation of the first exon and regulatory sequences, creates a previously nonexistent fusion transcript, which together are associated to expression changes in several other genes. Finally, we investigate the potential evolutionary and phenotypic consequences of the inversion, and suggest that it is probably deleterious. This is therefore the first example of a natural polymorphic inversion that has position effects and creates a new chimeric gene, contributing to answer an old question in evolutionary biology.

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.

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.

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.

Identifying the Association Between Alzheimer's Disease and Parkinson's Disease Using Genome-Wide Association Studies and Protein-Protein Interaction Network

Alzheimer's disease (AD) and Parkinson's disease (PD) are the first and second most common neurodegenerative diseases in the elderly. Shared clinical and pathological features have been reported. Recent large-scale genome-wide association studies (GWAS) have been conducted and reported a number of AD and PD variants. Until now, the underlying genetic mechanisms for all these newly identified PD variants as well as the association between AD and PD are still unclear exactly. We think that PD variants may contribute to AD and PD by influence on brain gene expression. Here, we conducted a systems analysis using (1) AD and PD variants (P < 5.00E-08) identified by the published GWAS; (2) four brain expression GWAS datasets using expression quantitative trait loci from the cerebellum and temporal cortex; (3) large-scale AD GWAS from the Alzheimer Disease Genetics Consortium (ADGC); (4) a protein-protein interaction network. Our results indicated that PD variants around the 17q21 were associated with gene expression and suggestive AD risk. We also identified significant interaction among AD and PD susceptibility genes. We believe that our findings may explain the underlying genetic mechanisms for newly identified PD variants in PD and AD, as well as the association between AD and PD, which may be very useful for future genetic studies for both diseases.

Evidence for dysregulation of genome-wide recombination in oocytes with nondisjoined chromosomes 21

In oocytes with nondisjoined chromosomes 21 due to a meiosis I (MI) error, recombination is significantly reduced along chromosome 21; several lines of evidence indicate that this contributes to the nondisjunction event. A pilot study found evidence that these oocytes also have reduced recombination genome-wide when compared with controls. This suggests that factors that act globally may be contributing to the reduced recombination on chromosome 21, and hence, the nondisjunction event. To identify the source of these factors, we examined two levels of recombination count regulation in oocytes: (i) regulation at the maternal level that leads to correlation in genome-wide recombination across her oocytes and (ii) regulation at the oocyte level that leads to correlation in recombination count among the chromosomes of an oocyte. We sought to determine whether either of these properties was altered in oocytes with an MI error. As it relates to maternal regulation, we found that both oocytes with an MI error (N = 94) and their siblings (N = 64) had reduced recombination when compared with controls (N = 2723). At the oocyte level, we found that the correlation in recombination count among the chromosomes of an oocyte is reduced in oocytes with MI errors compared with that of their siblings or controls. These results suggest that regulation at the maternal level predisposes MI error oocytes to reduced levels of recombination, but additional oocyte-specific dysregulation contributes to the nondisjunction event.