Genomic architecture and functional effects of potential human inversion supergenes

MAPT haplotype-associated transcriptomic changes in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a neurodegenerative movement and cognitive disorder characterized by abnormal accumulation of the microtubule-associated protein tau in the brain. Biochemically, inclusions in PSP are enriched for tau proteoforms with four microtubule-binding domain repeats (4R), an isoform that arises from alternative tau pre-mRNA splicing. While preferential aggregation and reduced degradation of 4R tau protein is thought to play a role in inclusion formation and toxicity, an alternative hypothesis is that altered expression of tau mRNA isoforms plays a causal role. This stems from the observation that PSP is associated with common variation in the tau gene (MAPT) at the 17q21.31 locus which contains low copy number repeats flanking a large recurrent genomic inversion. The complex genomic structural changes at the locus give rise to two dominant haplotypes, termed H1 and H2, that have the potential to markedly influence gene expression. Here, we explored haplotype-dependent differences in gene expression using a bulk RNA-seq dataset derived from human post-mortem brain tissue from PSP (n = 84) and controls (n = 77) using a rigorous computational pipeline, including alternative pre-mRNA splicing. We found 3579 differentially expressed genes in the temporal cortex and 10,011 in the cerebellum. We also found 7214 differential splicing events in the temporal cortex and 18,802 in the cerebellum. In the cerebellum, total tau mRNA levels and the proportion of transcripts encoding 4R tau were significantly increased in PSP compared to controls. In the temporal cortex, the proportion of reads that expressed 4R tau was increased in cases compared to controls. 4R tau mRNA levels were significantly associated with the H1 haplotype in the temporal cortex. Further, we observed a marked haplotype-dependent difference in KANSL1 expression that was strongly associated with H1 in both brain regions. These findings support the hypothesis that sporadic PSP is associated with haplotype-dependent increases in 4R tau mRNA that might play a causal role in this disorder.

Context-specific eQTLs reveal causal genes underlying shared genetic architecture of critically ill COVID-19 and idiopathic pulmonary fibrosis

Most genetic variants identified through genome-wide association studies (GWAS) are suspected to be regulatory in nature, but only a small fraction colocalize with expression quantitative trait loci (eQTLs, variants associated with expression of a gene). Therefore, it is hypothesized but largely untested that integration of disease GWAS with context-specific eQTLs will reveal the underlying genes driving disease associations. We used colocalization and transcriptomic analyses to identify shared genetic variants and likely causal genes associated with critically ill COVID-19 and idiopathic pulmonary fibrosis. We first identified five genome-wide significant variants associated with both diseases. Four of the variants did not demonstrate clear colocalization between GWAS and healthy lung eQTL signals. Instead, two of the four variants colocalized only in cell-type and disease-specific eQTL datasets. These analyses pointed to higher ATP11A expression from the C allele of rs12585036, in monocytes and in lung tissue from primarily smokers, which increased risk of IPF and decreased risk of critically ill COVID-19. We also found lower DPP9 expression (and higher methylation at a specific CpG) from the G allele of rs12610495, acting in fibroblasts and in IPF lungs, and increased risk of IPF and critically ill COVID-19. We further found differential expression of the identified causal genes in diseased lungs when compared to non-diseased lungs, specifically in epithelial and immune cell types. These findings highlight the power of integrating GWAS, context-specific eQTLs, and transcriptomics of diseased tissue to harness human genetic variation to identify causal genes and where they function during multiple diseases.

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.

Discovering genetic mechanisms underlying the co-occurrence of Parkinson's disease and non-motor traits

Understanding the biological mechanisms that underlie the non-motor symptoms of Parkinson's disease (PD) requires comprehensive frameworks that unravel the complex interplay of genetic risk factors. Here, we used a disease-agnostic brain cortex gene regulatory network integrated with Mendelian Randomization analyses that identified 19 genes whose changes in expression were causally linked to PD. We further used the network to identify genes that are regulated by PD-associated genome-wide association study (GWAS) SNPs. Extended protein interaction networks derived from PD-risk genes and PD-associated SNPs identified convergent impacts on biological pathways and phenotypes, connecting PD with established co-occurring traits, including non-motor symptoms. These findings hold promise for therapeutic development. In conclusion, while distinct sets of genes likely influence PD risk and outcomes, the existence of genes in common and intersecting pathways associated with other traits suggests that they may contribute to both increased PD risk and symptom heterogeneity observed in people with Parkinson's.

Impact of population structure in the estimation of recent historical effective population size by the software GONE

BACKGROUND

Effective population size (Ne) is a crucial parameter in conservation genetics and animal breeding. A recent method, implemented by the software GONE, has been shown to be rather accurate in estimating recent historical changes in Ne from a single sample of individuals. However, GONE estimations assume that the population being studied has remained isolated for a period of time, that is, without migration or confluence of other populations. If this occurs, the estimates of Ne can be heavily biased. In this paper, we evaluate the impact of migration and admixture on the estimates of historical Ne provided by GONE through a series of computer simulations considering several scenarios: (a) the mixture of two or more ancestral populations; (b) subpopulations that continuously exchange individuals through migration; (c) populations receiving migrants from a large source; and (d) populations with balanced systems of chromosomal inversions, which also generate genetic structure.

RESULTS

Our results indicate that the estimates of historical Ne provided by GONE may be substantially biased when there has been a recent mixture of populations that were previously separated for a long period of time. Similarly, biases may occur when the rate of continued migration between populations is low, or when chromosomal inversions are present at high frequencies. However, some biases due to population structuring can be eliminated by conducting population structure analyses and restricting the estimation to the differentiated groups. In addition, disregarding the genomic regions that are involved in inversions can also remove biases in the estimates of Ne.

CONCLUSIONS

Different kinds of deviations from isolation and panmixia of the populations can generate biases in the recent historical estimates of Ne. Therefore, estimation of past demography could benefit from performing population structure analyses beforehand, by mitigating the impact of these biases on historical Ne estimates.

The shared genetic architecture and evolution of human language and musical rhythm

Rhythm and language-related traits are phenotypically correlated, but their genetic overlap is largely unknown. Here, we leveraged two large-scale genome-wide association studies performed to shed light on the shared genetics of rhythm (N=606,825) and dyslexia (N=1,138,870). Our results reveal an intricate shared genetic and neurobiological architecture, and lay groundwork for resolving longstanding debates about the potential co-evolution of human language and musical traits.

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.

Genomic hotspots of chromosome rearrangements explain conserved synteny despite high rates of chromosome evolution in a holocentric lineage

Holocentric organisms, unlike typical monocentric organisms, have kinetochore activity distributed along almost the whole length of the chromosome. Because of this, chromosome rearrangements through fission and fusion are more likely to become fixed in holocentric species, which may account for the extraordinary rates of chromosome evolution that many holocentric lineages exhibit. Long blocks of genome synteny have been reported in animals with holocentric chromosomes despite high rates of chromosome rearrangements. Nothing is known from plants, however, despite the fact that holocentricity appears to have played a key role in the diversification of one of the largest angiosperm genera, Carex (Cyperaceae). In the current study, we compared genomes of Carex species and a distantly related Cyperaceae species to characterize conserved and rearranged genome regions. Our analyses span divergence times ranging between 2 and 50 million years. We also compared a C. scoparia chromosome-level genome assembly with a linkage map of the same species to study rearrangements at a population level and suppression of recombination patterns. We found longer genome synteny blocks than expected under a null model of random rearrangement breakpoints, even between very distantly related species. We also found repetitive DNA to be non-randomly associated with holocentromeres and rearranged regions of the genome. The evidence of conserved synteny in sedges despite high rates of chromosome fission and fusion suggests that conserved genomic hotspots of chromosome evolution related to repetitive DNA shape the evolution of recombination, gene order and crossability in sedges. This finding may help explain why sedges are able to maintain species cohesion even in the face of high interspecific chromosome rearrangements.

Effects of urban living environments on mental health in adults

Urban-living individuals are exposed to many environmental factors that may combine and interact to influence mental health. While individual factors of an urban environment have been investigated in isolation, no attempt has been made to model how complex, real-life exposure to living in the city relates to brain and mental health, and how this is moderated by genetic factors. Using the data of 156,075 participants from the UK Biobank, we carried out sparse canonical correlation analyses to investigate the relationships between urban environments and psychiatric symptoms. We found an environmental profile of social deprivation, air pollution, street network and urban land-use density that was positively correlated with an affective symptom group (r = 0.22, Pperm < 0.001), mediated by brain volume differences consistent with reward processing, and moderated by genes enriched for stress response, including CRHR1, explaining 2.01% of the variance in brain volume differences. Protective factors such as greenness and generous destination accessibility were negatively correlated with an anxiety symptom group (r = 0.10, Pperm < 0.001), mediated by brain regions necessary for emotion regulation and moderated by EXD3, explaining 1.65% of the variance. The third urban environmental profile was correlated with an emotional instability symptom group (r = 0.03, Pperm < 0.001). Our findings suggest that different environmental profiles of urban living may influence specific psychiatric symptom groups through distinct neurobiological pathways.

Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes

The human Major Histocompatibility Complex (MHC) or Human Leukocyte Antigen (HLA) super-locus is a highly polymorphic genomic region that encodes more than 140 coding genes including the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. In recent years, MHC genotyping and haplotyping using Sanger sequencing and next-generation sequencing (NGS) methods have produced many hundreds of genomic sequences of the HLA super-locus for comparative studies of the genetic architecture and diversity between the same and different haplotypes. In this special issue on 'The Current Landscape of HLA Genomics and Genetics', we provide a short review of some of the recent analytical developments used to investigate the SNP polymorphisms, structural variants (indels), transcription and haplotypes of the HLA super-locus. This review highlights the importance of using reference cell-lines, population studies, and NGS methods to improve and update our understanding of the mechanisms, architectural structures and combinations of human MHC genomic alleles (SNPs and indels) that better define and characterise haplotypes and their association with various phenotypes and diseases.

Genomic architecture of supergenes: connecting form and function

Supergenes are tightly linked sets of loci that are inherited together and control complex phenotypes. While classical supergenes-governing traits such as wing patterns in Heliconius butterflies or heterostyly in Primula-have been studied since the Modern Synthesis, we still understand very little about how they evolve and persist in nature. The genetic architecture of supergenes is a critical factor affecting their evolutionary fate, as it can change key parameters such as recombination rate and effective population size, potentially redirecting molecular evolution of the supergene in addition to the surrounding genomic region. To understand supergene evolution, we must link genomic architecture with evolutionary patterns and processes. This is now becoming possible with recent advances in sequencing technology and powerful forward computer simulations. The present theme issue brings together theoretical and empirical papers, as well as opinion and synthesis papers, which showcase the architectural diversity of supergenes and connect this to critical processes in supergene evolution, such as polymorphism maintenance and mutation accumulation. Here, we summarize those insights to highlight new ideas and methods that illuminate the path forward for the study of supergenes in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.