Diverse African genomes reveal selection on ancient modern human introgressions in Neanderthals

A structured coalescent model reveals deep ancestral structure shared by all modern humans

Understanding the history of admixture events and population size changes leading to modern humans is central to human evolutionary genetics. Here we introduce a coalescence-based hidden Markov model, cobraa, that explicitly represents an ancestral population split and rejoin, and demonstrate its application on simulated and real data across multiple species. Using cobraa, we present evidence for an extended period of structure in the history of all modern humans, in which two ancestral populations that diverged ~1.5 million years ago came together in an admixture event ~300 thousand years ago, in a ratio of ~80:20%. Immediately after their divergence, we detect a strong bottleneck in the major ancestral population. We inferred regions of the present-day genome derived from each ancestral population, finding that material from the minority correlates strongly with distance to coding sequence, suggesting it was deleterious against the majority background. Moreover, we found a strong correlation between regions of majority ancestry and human-Neanderthal or human-Denisovan divergence, suggesting the majority population was also ancestral to those archaic humans.

Uncovering the genetic architecture and evolutionary roots of androgenetic alopecia in African men

Androgenetic alopecia is a highly heritable trait. However, much of our understanding about the genetics of male-pattern baldness comes from individuals of European descent. Here, we examined a dataset comprising 2,136 men from Ghana, Nigeria, Senegal, and South Africa that were genotyped using the Men of African Descent and Carcinoma of the Prostate Array. We first tested how genetic predictions of baldness generalize from Europe to Africa and found that polygenic scores from European genome-wide association studies (GWASs) yielded area under the curve statistics that ranged from 0.513 to 0.546, indicating that genetic predictions of baldness generalized poorly from European to African populations. Subsequently, we conducted an African GWAS of androgenetic alopecia, focusing on self-reported baldness patterns at age 45. After correcting for age at recruitment, population structure, and study site, we identified 266 moderately significant associations, 51 of which were independent (p < 10-5, r2 < 0.2). Most baldness associations were autosomal, and the X chromosome does not seem to have a large impact on baldness in African men. Although Neanderthal alleles have previously been associated with skin and hair phenotypes, within the limits of statistical power, we did not find evidence that continental differences in the genetic architecture of baldness are due to Neanderthal introgression. While most loci that are associated with androgenetic alopecia do not have large integrative haplotype scores or fixation index statistics, multiple baldness-associated SNPs near the EDA2R and AR genes have large allele frequency differences between continents. Collectively, our findings illustrate how population genetic differences contribute to the limited portability of polygenic predictions across ancestries.

Archaic introgression and the distribution of shared variation under stabilizing selection

Many phenotypic traits are under stabilizing selection, which maintains a population's mean phenotypic value near some optimum. The dynamics of traits and trait architectures under stabilizing selection have been extensively studied for single populations at steady state. However, natural populations are seldom at steady state and are often structured in some way. Admixture and introgression events may be common, including over human evolutionary history. Because stabilizing selection results in selection against the minor allele at a trait-affecting locus, alleles from the minor parental ancestry will be selected against after admixture. We show that the site-frequency spectrum can be used to model the genetic architecture of such traits, allowing for the study of trait architecture dynamics in complex multi-population settings. We use a simple deterministic two-locus model to predict the reduction of introgressed ancestry around trait-contributing loci. From this and individual-based simulations, we show that introgressed-ancestry is depleted around such loci. When introgression between two diverged populations occurs in both directions, as has been inferred between humans and Neanderthals, the locations of such regions with depleted introgressed ancestry will tend to be shared across populations. We argue that stabilizing selection for shared phenotypic optima may explain recent observations in which regions of depleted human-introgressed ancestry in the Neanderthal genome overlap with Neanderthal-ancestry deserts in humans.

Archaic hominin admixture and its consequences for modern humans

As anatomically modern humans dispersed out of Africa, they encountered and mated with now extinct hominins, including Neanderthals and Denisovans. It is now well established that all non-African individuals derive approximately 2% of their genome from Neanderthal ancestors and individuals of Melanesian and Australian aboriginal ancestry inherited an additional 2%-5% of their genomes from Denisovan ancestors. Attention has started to shift from documenting amounts of archaic admixture and identifying introgressed segments to understanding their molecular, phenotypic, and evolutionary consequences and refining models of human history. Here, we review recent insights into admixture between modern and archaic humans, emphasizing methodological innovations and the functional and phenotypic effects Neanderthal and Denisovan sequences have in contemporary individuals.

Thematic Layers of Genomic Susceptibility for Conservation Monitoring

Population genomics has great potential to inform applied conservation management and associated policy. However, the bioinformatic analyses and interpretation of population genomic datasets can be daunting and difficult to convey to nonspecialists, including on-the-ground conservationists that work with many state, federal and international agencies. We think that individual population genomic metrics of interest can be interpolated and ultimately distilled into thematic GIS layers that represent spatiotemporal genomic potential (or conversely, susceptibility) in conservation monitoring. As examples relevant to ongoing conservation efforts, we use introgressive hybridisation and individual heterozygosity to illustrate a conceptual approach for mapping population genomic susceptibility. The general framework of thematic layers could be extended to integrate key genomic metrics (e.g., runs of homozygosity and genomic load) that are relevant to many conservation efforts.

Detecting adaptive changes in gene copy number distribution accompanying the human out-of-Africa expansion

Genes with multiple copies are likely to be maintained by stabilizing selection, which puts a bound to unlimited expansion of copy number. We designed a model in which copy number variation is generated by unequal recombination, which fits well with several genes surveyed in three human populations. Based on this theoretical model and computer simulations, we were interested in determining whether the gene copy number distribution in the derived European and Asian populations can be explained by a purely demographic scenario or whether shifts in the distribution are signatures of adaptation. Although the copy number distribution in most of the analyzed gene clusters can be explained by a bottleneck, such as in the out-of-Africa expansion of Homo sapiens 60-10 kyrs ago, we identified several candidate genes, such as AMY1A and PGA3, whose copy numbers are likely to differ among African, Asian, and European populations.

The evolutionary fate of Neanderthal DNA in 30,780 admixed genomes with recent African-like ancestry

Following introgression, Neanderthal DNA was initially purged from non-African genomes, but the evolutionary fate of remaining introgressed DNA has not been explored yet. To fill this gap, we analyzed 30,780 admixed genomes with African-like ancestry from the All of Us research program, in which Neanderthal alleles encountered novel genetic backgrounds during the last 15 generations. Observed amounts of Neanderthal DNA approximately match expectations based on ancestry proportions, suggesting neutral evolution. Nevertheless, we identified genomic regions that have significantly less or more Neanderthal ancestry than expected and are associated with spermatogenesis, innate immunity, and other biological processes. We also identified three novel introgression desert-like regions in recently admixed genomes, whose genetic features are compatible with hybrid incompatibilities and intrinsic negative selection. Overall, we find that much of the remaining Neanderthal DNA in human genomes is not under strong selection, and complex evolutionary dynamics have shaped introgression landscapes in our species.

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

Ancient and recent origins of shared polymorphisms in yeast

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Uncovering the genetic architecture and evolutionary roots of androgenetic alopecia in African men

Androgenetic alopecia is a highly heritable trait. However, much of our understanding about the genetics of male pattern baldness comes from individuals of European descent. Here, we examined a novel dataset comprising 2,136 men from Ghana, Nigeria, Senegal, and South Africa that were genotyped using a custom array. We first tested how genetic predictions of baldness generalize from Europe to Africa, finding that polygenic scores from European GWAS yielded AUC statistics that ranged from 0.513 to 0.546, indicating that genetic predictions of baldness in African populations performed notably worse than in European populations. Subsequently, we conducted the first African GWAS of androgenetic alopecia, focusing on self-reported baldness patterns at age 45. After correcting for present age, population structure, and study site, we identified 266 moderately significant associations, 51 of which were independent (p-value < 10-5, r2 < 0.2). Most baldness associations were autosomal, and the X chromosomes does not appear to have a large impact on baldness in African men. Finally, we examined the evolutionary causes of continental differences in genetic architecture. Although Neanderthal alleles have previously been associated with skin and hair phenotypes, we did not find evidence that European-ascertained baldness hits were enriched for signatures of ancient introgression. Most loci that are associated with androgenetic alopecia are evolving neutrally. However, multiple baldness-associated SNPs near the EDA2R and AR genes have large allele frequency differences between continents. Collectively, our findings illustrate how evolutionary history contributes to the limited portability of genetic predictions across ancestries.

Human evolution: Neanderthal footprints in African genomes

Human and Neanderthal populations met and mixed on multiple occasions over evolutionary time, resulting in the exchange of genetic material. New genomic analyses of diverse African populations reveal a history of bidirectional gene flow and selection acting on introgressed alleles.