Highly punctuated patterns of population structure on the X chromosome and implications for African evolutionary history

No evidence that ACE2 or TMPRSS2 drive population disparity in COVID risks

Early in the SARS-CoV2 pandemic, in this journal, Hou et al. (BMC Med 18:216, 2020) interpreted public genotype data, run through functional prediction tools, as suggesting that members of particular human populations carry potentially COVID-risk-increasing variants in genes ACE2 and TMPRSS2 far more often than do members of other populations. Beyond resting on predictions rather than clinical outcomes, and focusing on variants too rare to typify population members even jointly, their claim mistook a well known artifact (that large samples reveal more of a population's variants than do small samples) as if showing real and congruent population differences for the two genes, rather than lopsided population sampling in their shared source data. We explain that artifact, and contrast it with empirical findings, now ample, that other loci shape personal COVID risks far more significantly than do ACE2 and TMPRSS2-and that variation in ACE2 and TMPRSS2 per se unlikely exacerbates any net population disparity in the effects of such more risk-informative loci.

Joint identification of sex and sex-linked scaffolds in non-model organisms using low depth sequencing data

Being able to assign sex to individuals and identify autosomal and sex-linked scaffolds are essential in most population genomic analyses. Non-model organisms often have genome assemblies at scaffold-level and lack characterization of sex-linked scaffolds. Previous methods to identify sex and sex-linked scaffolds have relied on synteny between the non-model organism and a closely related species or prior knowledge about the sex of the samples to identify sex-linked scaffolds. In the latter case, the difference in depth of coverage between the autosomes and the sex chromosomes are used. Here, we present "sex assignment through coverage" (SATC), a method to assign sex to samples and identify sex-linked scaffolds from next generation sequencing (NGS) data. The method works for species with a homogametic/heterogametic sex determination system and only requires a scaffold-level reference assembly and sampling of both sexes with whole genome sequencing (WGS) data. We use the sequencing depth distribution across scaffolds to jointly identify: (i) male and female individuals, and (ii) sex-linked scaffolds. This is achieved through projecting the scaffold depths into a low-dimensional space using principal component analysis (PCA) and subsequent Gaussian mixture clustering. We demonstrate the applicability of our method using data from five mammal species and a bird species complex. The method is freely available at https://github.com/popgenDK/SATC as R code and a graphical user interface (GUI).

Common X-Chromosome Variants Are Associated with Parkinson Disease Risk

OBJECTIVE

The objective of this study was to identify genetic variants on the X-chromosome associated with Parkinson disease (PD) risk.

METHODS

We performed an X-chromosome-wide association study (XWAS) of PD risk by meta-analyzing results from sex-stratified analyses. To avoid spurious associations, we designed a specific harmonization pipeline for the X-chromosome and focused on a European ancestry sample. We included 11,142 cases, 280,164 controls, and 5,379 proxy cases, based on parental history of PD. Additionally, we tested the association of significant variants with (1) PD risk in an independent replication with 1,561 cases and 2,465 controls and (2) putamen volume in 33,360 individuals from the UK Biobank.

RESULTS

In the discovery meta-analysis, we identified rs7066890 (odds ratio [OR] = 1.10, 95% confidence interval [CI] = 1.06-1.14, p = 2.2 × 10-9 ), intron of GPM6B, and rs28602900 (OR = 1.10, 95% CI = 1.07-1.14, p = 1.6 × 10-8 ) in a high gene density region including RPL10, ATP6A1, FAM50A, and PLXNA3. The rs28602900 association with PD was replicated (OR = 1.16, 95% CI = 1.03-1.30, p = 0.016) and shown to colocalize with a significant expression quantitative locus (eQTL) regulating RPL10 expression in the putamen and other brain tissues in the Genotype-Tissue Expression Project. Additionally, the rs28602900 locus was found to be associated with reduced brain putamen volume. No results reached genome-wide significance in the sex-stratified analyses.

INTERPRETATION

We report the first XWAS of PD and identify 2 genome-wide significant loci. The rs28602900 association was replicated in an independent PD dataset and showed concordant effects in its association with putamen volume. Critically, rs26802900 is a significant eQTL of RPL10. These results support a role for ribosomal proteins in PD pathogenesis and show that the X-chromosome contributes to PD genetic risk. ANN NEUROL 2021;90:22-34.

Quantifying the relative contributions of the X chromosome, autosomes, and mitochondrial genome to local adaptation

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X-linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline-end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly-a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.

IBD Sharing between Africans, Neandertals, and Denisovans

Interbreeding between ancestors of humans and other hominins outside of Africa has been studied intensively, while their common history within Africa still lacks proper attention. However, shedding light on human evolution in this time period about which little is known, is essential for understanding subsequent events outside of Africa. We investigate the genetic relationships of humans, Neandertals, and Denisovans by identifying very short DNA segments in the 1000 Genomes Phase 3 data that these hominins share identical by descent (IBD). By focusing on low frequency and rare variants, we identify very short IBD segments with high confidence. These segments reveal events from a very distant past because shorter IBD segments are presumably older than longer ones. We extracted two types of very old IBD segments that are not only shared among humans, but also with Neandertals and/or Denisovans. The first type contains longer segments that are found primarily in Asians and Europeans where more segments are found in South Asians than in East Asians for both Neandertal and Denisovan. These longer segments indicate complex admixture events outside of Africa. The second type consists of shorter segments that are shared mainly by Africans and therefore may indicate events involving ancestors of humans and other ancient hominins within Africa. Our results from the autosomes are further supported by an analysis of chromosome X, on which segments that are shared by Africans and match the Neandertal and/or Denisovan genome were even more prominent. Our results indicate that interbreeding with other hominins was a common feature of human evolution starting already long before ancestors of modern humans left Africa.

The fire ant social chromosome supergene variant Sb shows low diversity but high divergence from SB

Variation in social behaviour is common, yet little is known about the genetic architectures underpinning its evolution. A rare exception is in the fire ant Solenopsis invicta: Alternative variants of a supergene region determine whether a colony will have exactly one or up to dozens of queens. The two variants of this region are carried by a pair of 'social chromosomes', SB and Sb, which resemble a pair of sex chromosomes. Recombination is suppressed between the two chromosomes in the supergene region. While the X-like SB can recombine with itself in SB/SB queens, recombination is effectively absent in the Y-like Sb because Sb/Sb queens die before reproducing. Here, we analyse whole-genome sequences of eight haploid SB males and eight haploid Sb males. We find extensive SB-Sb differentiation throughout the >19-Mb-long supergene region. We find no evidence of 'evolutionary strata' with different levels of divergence comparable to those reported in several sex chromosomes. A high proportion of substitutions between the SB and Sb haplotypes are nonsynonymous, suggesting inefficacy of purifying selection in Sb sequences, similar to that for Y-linked sequences in XY systems. Finally, we show that the Sb haplotype of the supergene region has 635-fold less nucleotide diversity than the rest of the genome. We discuss how this reduction could be due to a recent selective sweep affecting Sb specifically or associated with a population bottleneck during the invasion of North America by the sampled population.

The Relative Contributions of the X Chromosome and Autosomes to Local Adaptation

Models of sex chromosome and autosome evolution yield key predictions about the genomic basis of adaptive divergence, and such models have been important in guiding empirical research in comparative genomics and studies of speciation. In addition to the adaptive differentiation that occurs between species over time, selection also favors genetic divergence across geographic space, with subpopulations of single species evolving conspicuous differences in traits involved in adaptation to local environmental conditions. The potential contribution of sex chromosomes (the X or Z) to local adaptation remains unclear, as we currently lack theory that directly links spatial variation in selection to local adaptation of X-linked and autosomal genes. Here, we develop population genetic models that explicitly consider the effects of genetic dominance, effective population size, and sex-specific migration and selection on the relative contributions of X-linked and autosomal genes to local adaptation. We show that X-linked genes should nearly always disproportionately contribute to local adaptation in the presence of gene flow. We also show that considerations of dominance and effective population size-which play pivotal roles in the theory of faster-X adaptation between species-have surprisingly little influence on the relative contribution of the X chromosome to local adaptation. Instead, sex-biased migration is the primary mediator of the strength of spatial large-X effects. Our results yield novel predictions about the role of sex chromosomes in local adaptation. We outline empirical approaches in evolutionary quantitative genetics and genomics that could build upon this new theory.

XWAS: A Software Toolset for Genetic Data Analysis and Association Studies of the X Chromosome

XWAS is a new software suite for the analysis of the X chromosome in association studies and similar genetic studies. The X chromosome plays an important role in human disease and traits of many species, especially those with sexually dimorphic characteristics. Special attention needs to be given to its analysis due to the unique inheritance pattern, which leads to analytical complications that have resulted in the majority of genome-wide association studies (GWAS) either not considering X or mishandling it with toolsets that had been designed for non-sex chromosomes. We hence developed XWAS to fill the need for tools that are specially designed for analysis of X. Following extensive, stringent, and X-specific quality control, XWAS offers an array of statistical tests of association, including: 1) the standard test between a SNP (single nucleotide polymorphism) and disease risk, including after first stratifying individuals by sex, 2) a test for a differential effect of a SNP on disease between males and females, 3) motivated by X-inactivation, a test for higher variance of a trait in heterozygous females as compared with homozygous females, and 4) for all tests, a version that allows for combining evidence from all SNPs across a gene. We applied the toolset analysis pipeline to 16 GWAS datasets of immune-related disorders and 7 risk factors of coronary artery disease, and discovered several new X-linked genetic associations. XWAS will provide the tools and incentive for others to incorporate the X chromosome into GWAS and similar studies in any species with an XX/XY system, hence enabling discoveries of novel loci implicated in many diseases and in their sexual dimorphism.

Beyond 2/3 and 1/3: The Complex Signatures of Sex-Biased Admixture on the X Chromosome

Sex-biased demography, in which parameters governing migration and population size differ between females and males, has been studied through comparisons of X chromosomes, which are inherited sex-specifically, and autosomes, which are not. A common form of sex bias in humans is sex-biased admixture, in which at least one of the source populations differs in its proportions of females and males contributing to an admixed population. Studies of sex-biased admixture often examine the mean ancestry for markers on the X chromosome in relation to the autosomes. A simple framework noting that in a population with equally many females and males, two-thirds of X chromosomes appear in females, suggests that the mean X-chromosomal admixture fraction is a linear combination of female and male admixture parameters, with coefficients 2/3 and 1/3, respectively. Extending a mechanistic admixture model to accommodate the X chromosome, we demonstrate that this prediction is not generally true in admixture models, although it holds in the limit for an admixture process occurring as a single event. For a model with constant ongoing admixture, we determine the mean X-chromosomal admixture, comparing admixture on female and male X chromosomes to corresponding autosomal values. Surprisingly, in reanalyzing African-American genetic data to estimate sex-specific contributions from African and European sources, we find that the range of contributions compatible with the excess African ancestry on the X chromosome compared to autosomes has a wide spread, permitting scenarios either without male-biased contributions from Europe or without female-biased contributions from Africa.

Autosomal admixture levels are informative about sex bias in admixed populations

Sex-biased admixture has been observed in a wide variety of admixed populations. Genetic variation in sex chromosomes and functions of quantities computed from sex chromosomes and autosomes have often been examined to infer patterns of sex-biased admixture, typically using statistical approaches that do not mechanistically model the complexity of a sex-specific history of admixture. Here, expanding on a model of Verdu and Rosenberg (2011) that did not include sex specificity, we develop a model that mechanistically examines sex-specific admixture histories. Under the model, multiple source populations contribute to an admixed population, potentially with their male and female contributions varying over time. In an admixed population descended from two source groups, we derive the moments of the distribution of the autosomal admixture fraction from a specific source population as a function of sex-specific introgression parameters and time. Considering admixture processes that are constant in time, we demonstrate that surprisingly, although the mean autosomal admixture fraction from a specific source population does not reveal a sex bias in the admixture history, the variance of autosomal admixture is informative about sex bias. Specifically, the long-term variance decreases as the sex bias from a contributing source population increases. This result can be viewed as analogous to the reduction in effective population size for populations with an unequal number of breeding males and females. Our approach suggests that it may be possible to use the effect of sex-biased admixture on autosomal DNA to assist with methods for inference of the history of complex sex-biased admixture processes.

A tale of two haplotypes: the EDA2R/AR Intergenic region is the most divergent genomic segment between Africans and East Asians in the human genome

Single nucleotide polymorphisms (SNPs) with large allele frequency differences between human populations are relatively rare. The longest run of SNPs with an allele frequency difference of one between the Yoruba of Nigeria and the Han Chinese is found on the long arm of the X chromosome in the intergenic region separating the EDA2R and AR genes. It has been proposed that the unusual allele frequency distributions of these SNPs are the result of a selective sweep affecting African populations that occurred after the out-of-Africa migration. To investigate the evolutionary history of the EDA2R/AR intergenic region, we characterized the haplotype structure of 52 of its highly differentiated SNPs. Using a publicly available data set of 3,000 X chromosomes from 65 human populations, we found that nearly all human X chromosomes carry one of two modal haplotypes for these 52 SNPs. The predominance of two highly divergent haplotypes at this locus was confirmed by use of a subset of individuals sequenced to high coverage. The first of these haplotypes, the α-haplotype is at high frequencies in most of the African populations surveyed and likely arose before the separation of African populations into distinct genetic entities. The second, the β-haplotype, is frequent or fixed in all non-African populations and likely arose in East Africa before the out-of-Africa migration. We also observed a small group or rare haplotypes with no clear relationship to the α- and β-haplotypes. These haplotypes occur at relatively high frequencies in African hunter-gatherer populations, such as the San and Mbuti Pygmies. Our analysis indicates that these haplotypes are part of a pool of diverse, ancestral haplotypes that have now been almost entirely replaced by the α- and β-haplotypes. We suggest that the rise of the α- and β-haplotypes was the result of the demographic forces that human populations experienced during the formation of modern African populations and the out-of-Africa migration. However, we also present evidence that this region is the target of selection in the form of positive selection on the α- and β-haplotypes and of purifying the selection against α/β recombinants.

Contrasting X-linked and autosomal diversity across 14 human populations

Contrasting the genetic diversity of the human X chromosome (X) and autosomes has facilitated understanding historical differences between males and females and the influence of natural selection. Previous studies based on smaller data sets have left questions regarding how empirical patterns extend to additional populations and which forces can explain them. Here, we address these questions by analyzing the ratio of X-to-autosomal (X/A) nucleotide diversity with the complete genomes of 569 females from 14 populations. Results show that X/A diversity is similar within each continental group but notably lower in European (EUR) and East Asian (ASN) populations than in African (AFR) populations. X/A diversity increases in all populations with increasing distance from genes, highlighting the stronger impact of diversity-reducing selection on X than on the autosomes. However, relative X/A diversity (between two populations) is invariant with distance from genes, suggesting that selection does not drive the relative reduction in X/A diversity in non-Africans (0.842 ± 0.012 for EUR-to-AFR and 0.820 ± 0.032 for ASN-to-AFR comparisons). Finally, an array of models with varying population bottlenecks, expansions, and migration from the latest studies of human demographic history account for about half of the observed reduction in relative X/A diversity from the expected value of 1. They predict values between 0.91 and 0.94 for EUR-to-AFR comparisons and between 0.91 and 0.92 for ASN-to-AFR comparisons. Further reductions can be predicted by more extreme demographic events in excess of those captured by the latest studies but, in the absence of these, also by historical sex-biased demographic events or other processes.

X-linked MTMR8 diversity and evolutionary history of sub-Saharan populations

The genetic diversity within an 11 kb segment of the MTMR8 gene in a sample of 111 sub-Saharan and 49 non-African X chromosomes was investigated to assess the early evolutionary history of sub-Saharan Africans and the out-of-Africa expansion. The analyses revealed a complex genetic structure of the Africans that contributed to the emergence of modern humans. We observed partitioning of two thirds of old lineages among southern, west/central and east African populations indicating ancient population stratification predating the out of Africa migration. Age estimates of these lineages, older than coalescence times of uniparentally inherited markers, raise the question whether contemporary humans originated from a single population or as an amalgamation of different populations separated by years of independent evolution, thus suggesting a greater antiquity of our species than generally assumed. While the oldest sub-Saharan lineages, ∼500 thousand years, are found among Khoe-San from southern-Africa, a distinct haplotype found among Biaka is likely due to admixture from an even older population. An East African population that gave rise to non-Africans underwent a selective sweep affecting the subcentromeric region where MTMR8 is located. This and similar sweeps in four other regions of the X chromosome, documented in the literature, effectively reduced genetic diversity of non-African chromosomes and therefore may have exacerbated the effect of the demographic bottleneck usually ascribed to the out of Africa migration. Our data is suggestive, however, that a bottleneck, occurred in Africa before range expansion.

HMOX1 gene promoter alleles and high HO-1 levels are associated with severe malaria in Gambian children

Heme oxygenase 1 (HO-1) is an essential enzyme induced by heme and multiple stimuli associated with critical illness. In humans, polymorphisms in the HMOX1 gene promoter may influence the magnitude of HO-1 expression. In many diseases including murine malaria, HO-1 induction produces protective anti-inflammatory effects, but observations from patients suggest these may be limited to a narrow range of HO-1 induction, prompting us to investigate the role of HO-1 in malaria infection. In 307 Gambian children with either severe or uncomplicated P. falciparum malaria, we characterized the associations of HMOX1 promoter polymorphisms, HMOX1 mRNA inducibility, HO-1 protein levels in leucocytes (flow cytometry), and plasma (ELISA) with disease severity. The (GT)(n) repeat polymorphism in the HMOX1 promoter was associated with HMOX1 mRNA expression in white blood cells in vitro, and with severe disease and death, while high HO-1 levels were associated with severe disease. Neutrophils were the main HO-1-expressing cells in peripheral blood, and HMOX1 mRNA expression was upregulated by heme-moieties of lysed erythrocytes. We provide mechanistic evidence that induction of HMOX1 expression in neutrophils potentiates the respiratory burst, and propose this may be part of the causal pathway explaining the association between short (GT)(n) repeats and increased disease severity in malaria and other critical illnesses. Our findings suggest a genetic predisposition to higher levels of HO-1 is associated with severe illness, and enhances the neutrophil burst leading to oxidative damage of endothelial cells. These add important information to the discussion about possible therapeutic manipulation of HO-1 in critically ill patients.

Estimators of the human effective sex ratio detect sex biases on different timescales

Determining historical sex ratios throughout human evolution can provide insight into patterns of genomic variation, the structure and composition of ancient populations, and the cultural factors that influence the sex ratio (e.g., sex-specific migration rates). Although numerous studies have suggested that unequal sex ratios have existed in human evolutionary history, a coherent picture of sex-biased processes has yet to emerge. For example, two recent studies compared human X chromosome to autosomal variation to make inferences about historical sex ratios but reached seemingly contradictory conclusions, with one study finding evidence for a male bias and the other study identifying a female bias. Here, we show that a large part of this discrepancy can be explained by methodological differences. Specifically, through reanalysis of empirical data, derivation of explicit analytical formulae, and extensive simulations we demonstrate that two estimators of the effective sex ratio based on population structure and nucleotide diversity preferentially detect biases that have occurred on different timescales. Our results clarify apparently contradictory evidence on the role of sex-biased processes in human evolutionary history and show that extant patterns of human genomic variation are consistent with both a recent male bias and an earlier, persistent female bias.

Genetic admixture and population substructure in Guanacaste Costa Rica

The population of Costa Rica (CR) represents an admixture of major continental populations. An investigation of the CR population structure would provide an important foundation for mapping genetic variants underlying common diseases and traits. We conducted an analysis of 1,301 women from the Guanacaste region of CR using 27,904 single nucleotide polymorphisms (SNPs) genotyped on a custom Illumina InfiniumII iSelect chip. The program STRUCTURE was used to compare the CR Guanacaste sample with four continental reference samples, including HapMap Europeans (CEU), East Asians (JPT+CHB), West African Yoruba (YRI), as well as Native Americans (NA) from the Illumina iControl database. Our results show that the CR Guanacaste sample comprises a three-way admixture estimated to be 43% European, 38% Native American and 15% West African. An estimated 4% residual Asian ancestry may be within the error range. Results from principal components analysis reveal a correlation between genetic and geographic distance. The magnitude of linkage disequilibrium (LD) measured by the number of tagging SNPs required to cover the same region in the genome in the CR Guanacaste sample appeared to be weaker than that observed in CEU, JPT+CHB and NA reference samples but stronger than that of the HapMap YRI sample. Based on the clustering pattern observed in both STRUCTURE and principal components analysis, two subpopulations were identified that differ by approximately 20% in LD block size averaged over all LD blocks identified by Haploview. We also show in a simulated association study conducted within the two subpopulations, that the failure to account for population stratification (PS) could lead to a noticeable inflation in the false positive rate. However, we further demonstrate that existing PS adjustment approaches can reduce the inflation to an acceptable level for gene discovery.

Genetic structure in African populations: implications for human demographic history

The continent of Africa is the source of all anatomically modern humans that dispersed across the planet during the past 100,000 years. As such, African populations are characterized by high genetic diversity and low levels of linkage disequilibrium (LD) among loci, as compared to populations from other continents. African populations also possess a number of genetic adaptations that have evolved in response to the diverse climates, diets, geographic environments, and infectious agents that characterize the African continent. Recently, Tishkoff et al. (2009) performed a genome-wide analysis of substructure based on DNA from 2432 Africans from 121 geographically diverse populations. The authors analyzed patterns of variation at 1327 nuclear microsatellite and insertion/deletion markers and identified 14 ancestral population clusters that correlate well with self-described ethnicity and shared cultural or linguistic properties. The results suggest that African populations may have maintained a large and subdivided population structure throughout much of their evolutionary history. In this chapter, we synthesize recent work documenting evidence of African population structure and discuss the implications for inferences about evolutionary history in both African populations and anatomically modern humans as a whole.