Whole-genome sequence analyses of Western Central African Pygmy hunter-gatherers reveal a complex demographic history and identify candidate genes under positive natural selection

Expanding the stdpopsim species catalog, and lessons learned for realistic genome simulations

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

Advances in integrative African genomics

There has been a rapid increase in human genome sequencing in the past two decades, resulting in the identification of millions of previously unknown genetic variants. However, African populations are under-represented in sequencing efforts. Additional sequencing from diverse African populations and the construction of African-specific reference genomes is needed to better characterize the full spectrum of variation in humans. However, sequencing alone is insufficient to address the molecular and cellular mechanisms underlying variable phenotypes and disease risks. Determining functional consequences of genetic variation using multi-omics approaches is a fundamental post-genomic challenge. We discuss approaches to close the knowledge gaps about African genomic diversity and review advances in African integrative genomic studies and their implications for precision medicine.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

The evolution of human skin pigmentation involved the interactions of genetic, environmental, and cultural variables

The primary biological role of human skin pigmentation is as a mediator of penetration of ultraviolet radiation (UVR) into the deep layers of skin and the cutaneous circulation. Since the origin of Homo sapiens, dark, protective constitutive pigmentation and strong tanning abilities have been favored under conditions of high UVR and represent the baseline condition for modern humans. The evolution of partly depigmented skin and variable tanning abilities has occurred multiple times in prehistory, as populations have dispersed into environments with lower and more seasonal UVR regimes, with unique complements of genes and cultural practices. The evolution of extremes of dark pigmentation and depigmentation has been rare and occurred only under conditions of extremely high or low environmental UVR, promoted by positive selection on variant pigmentation genes followed by limited gene flow. Over time, the evolution of human skin pigmentation has been influenced by the nature and course of human dispersals and modifications of cultural practices, which have modified the nature and actions of skin pigmentation genes. Throughout most of prehistory and history, the evolution of human skin pigmentation has been a contingent and non-deterministic process.

Genomic inference of contemporary effective population size in a large island population of collared flycatchers (Ficedula albicollis)

Due to its central importance to many aspects of evolutionary biology and population genetics, the long-term effective population size (N_e ) has been estimated for numerous species and populations. However, estimating contemporary N_e is difficult and in practice this parameter is often unknown. In principle, contemporary N_e can be estimated using either analyses of temporal changes in allele frequencies, or the extent of linkage disequilibrium (LD) between unlinked markers. We applied these approaches to estimate contemporary N_e of a relatively recently founded island population of collared flycatchers (Ficedula albicollis). We sequenced the genomes of 85 birds sampled in 1993 and 2015, and applied several temporal methods to estimate N_e at a few thousand (4000-7000). The approach based on LD provided higher estimates of N_e (20,000-32,000) and was associated with high variance, often resulting in infinite N_e . We conclude that whole-genome sequencing data offers new possibilities to estimate high (>1000) contemporary N_e , but also note that such estimates remain challenging, in particular for LD-based methods for contemporary N_e estimation.

Characterization of structural variation in Tibetans reveals new evidence of high-altitude adaptation and introgression

BACKGROUND

Structural variation (SV) acts as an essential mutational force shaping the evolution and function of the human genome. However, few studies have examined the role of SVs in high-altitude adaptation and little is known of adaptive introgressed SVs in Tibetans so far.

RESULTS

Here, we generate a comprehensive catalog of SVs in a Chinese Tibetan (n = 15) and Han (n = 10) population using nanopore sequencing technology. Among a total of 38,216 unique SVs in the catalog, 27% are sequence-resolved for the first time. We systematically assess the distribution of these SVs across repeat sequences and functional genomic regions. Through genotyping in additional 276 genomes, we identify 69 Tibetan-Han stratified SVs and 80 candidate adaptive genes. We also discover a few adaptive introgressed SV candidates and provide evidence for a deletion of 335 base pairs at 1p36.32.

CONCLUSIONS

Overall, our results highlight the important role of SVs in the evolutionary processes of Tibetans' adaptation to the Qinghai-Tibet Plateau and provide a valuable resource for future high-altitude adaptation studies.

The distinct morphological phenotypes of Southeast Asian aborigines are shaped by novel mechanisms for adaptation to tropical rainforests

Southeast Asian aborigines, the hunter-gatherer populations living in tropical rainforests, exhibit distinct morphological phenotypes, including short stature, dark skin, curly hair and a wide and snub nose. The underlying genetic architecture and evolutionary mechanism of these phenotypes remain a long-term mystery. We conducted whole genome deep sequencing of 81 Cambodian aborigines from eight ethnic groups. Through a genome-wide scan of selective sweeps, we discovered key genes harboring Cambodian-enriched mutations that may contribute to their phenotypes, including two hair morphogenesis genes (TCHH and TCHHL1), one nasal morphology gene (PAX3) and a set of genes (such as ENTPD1-AS1) associated with short stature. The identified new genes and novel mutations suggest an independent origin of the distinct phenotypes in Cambodian aborigines through parallel evolution, refuting the long-standing argument on the common ancestry of these phenotypes among the worldwide rainforest hunter-gatherers. Notably, our discovery reveals that various types of molecular mechanisms, including antisense transcription and epigenetic regulation, contribute to human morphogenesis, providing novel insights into the genetics of human environmental adaptation.

The deep population history in Africa

Africa is the continent with the greatest genetic diversity among humans and the level of diversity is further enhanced by incorporating non-majority groups, which are often understudied. Many of today's minority populations historically practiced foraging lifestyles, which were the only subsistence strategies prior to the rise of agriculture and pastoralism, but only a few groups practicing these strategies remain today. Genomic investigations of Holocene human remains excavated across the African continent show that the genetic landscape was vastly different compared to today's genetic landscape and that many groups that today are population isolate inhabited larger regions in the past. It is becoming clear that there are periods of isolation among groups and geographic areas, but also genetic contact over large distances throughout human history in Africa. Genomic information from minority populations and from prehistoric remains provide an invaluable source of information on the human past, in particular deep human population history, as Holocene large-scale population movements obscure past patterns of population structure. Here we revisit questions on the nature and time of the radiation of early humans in Africa, the extent of gene-flow among human populations as well as introgression from archaic and extinct lineages on the continent.

High Levels of Genetic Diversity within Nilo-Saharan Populations: Implications for Human Adaptation

Africa contains more human genetic variation than any other continent, but the majority of the population-scale analyses of the African peoples have focused on just two of the four major linguistic groups, the Niger-Congo and Afro-Asiatic, leaving the Nilo-Saharan and Khoisan populations under-represented. In order to assess genetic variation and signatures of selection within a Nilo-Saharan population and between the Nilo-Saharan and Niger-Congo and Afro-Asiatic, we sequenced 50 genomes from the Nilo-Saharan Lugbara population of North-West Uganda and 250 genomes from 6 previously unsequenced Niger-Congo populations. We compared these data to data from a further 16 Eurasian and African populations including the Gumuz, another putative Nilo-Saharan population from Ethiopia. Of the 21 million variants identified in the Nilo-Saharan population, 3.57 million (17%) were not represented in dbSNP and included predicted non-synonymous mutations with possible phenotypic effects. We found greater genetic differentiation between the Nilo-Saharan Lugbara and Gumuz populations than between any two Afro-Asiatic or Niger-Congo populations. F3 tests showed that Gumuz contributed a genetic component to most Niger-Congo B populations whereas Lugabara did not. We scanned the genomes of the Lugbara for evidence of selective sweeps. We found selective sweeps at four loci (SLC24A5, SNX13, TYRP1, and UVRAG) associated with skin pigmentation, three of which already have been reported to be under selection. These selective sweeps point toward adaptations to the intense UV radiation of the Sahel.

A Late Pleistocene human humerus from Rusinga Island, Lake Victoria, Kenya

In 2010, a hominin right humerus fragment (KNM-RU 58330) was surface collected in a small gully at Nyamita North in the Late Pleistocene Wasiriya Beds of Rusinga Island, Kenya. A combination of stratigraphic and geochronological evidence suggests the specimen is likely between ∼49 and 36 ka in age. The associated fauna is diverse and dominated by semiarid grassland taxa. The small sample of associated Middle Stone Age artifacts includes Levallois flakes, cores, and retouched points. The 139 mm humeral fragment preserves the shaft from distal to the lesser tubercle to 14 mm below the distal end of the weakly projecting deltoid tuberosity. Key morphological features include a narrow and weakly marked pectoralis major insertion and a distinctive medial bend in the diaphysis at the deltoid insertion. This bend is unusual among recent human humeri but occurs in a few Late Pleistocene humeri. The dimensions of the distal end of the fragment predict a length of 317.9 ± 16.4 mm based on recent samples of African ancestry. A novel method of predicting humeral length from the distance between the middle of the pectoralis major and the bottom of the deltoid insertion predicts a length of 317.3 mm ± 17.6 mm. Cross-sectional geometry at the midshaft shows a relatively high percentage of cortical bone and a moderate degree of flattening of the shaft. The Nyamita humerus is anatomically modern in its morphology and adds to the small sample of hominins from the Late Pleistocene associated with Middle Stone Age artifacts known from East Africa. It may sample a population closely related to the people of the out-of-Africa migration.

Identification of novel genetic variants associated with short stature in a Baka Pygmies population

Human growth is a complex trait determined by genetic factors in combination with external stimuli, including environment, nutrition and hormonal status. In the past, several genome-wide association studies (GWAS) have collectively identified hundreds of genetic variants having a putative effect on determining adult height in different worldwide populations. Theoretically, a valuable approach to better understand the mechanisms of complex traits as adult height is to study a population exhibiting extreme stature phenotypes, such as African Baka Pygmies. After phenotypic characterization, we sequenced the whole exomes of a cohort of Baka Pygmies and their non-Pygmies Bantu neighbors to highlight genetic variants associated with the reduced stature. Whole exome data analysis revealed 29 single nucleotide polymorphisms (SNPs) significantly associated with the reduced height in the Baka group. Among these variants, we focused on SNP rs7629425, located in the 5′-UTR of the Hyaluronidase-2 (HYAL2) gene. The frequency of the alternative allele was significantly increased compared to African and non-African populations. In vitro luciferase assay showed significant differences in transcription modulation by rs7629425 C/T alleles. In conclusion, our results suggested that the HYAL2 gene variants may play a role in the etiology of short stature in Baka Pygmies population.

Adaptive archaic introgression of copy number variants and the discovery of previously unknown human genes

Copy number variants (CNVs) are subject to stronger selective pressure than single-nucleotide variants, but their roles in archaic introgression and adaptation have not been systematically investigated. We show that stratified CNVs are significantly associated with signatures of positive selection in Melanesians and provide evidence for adaptive introgression of large CNVs at chromosomes 16p11.2 and 8p21.3 from Denisovans and Neanderthals, respectively. Using long-read sequence data, we reconstruct the structure and complex evolutionary history of these polymorphisms and show that both encode positively selected genes absent from most human populations. Our results collectively suggest that large CNVs originating in archaic hominins and introgressed into modern humans have played an important role in local population adaptation and represent an insufficiently studied source of large-scale genetic variation.

From Summary Statistics to Gene Trees: Methods for Inferring Positive Selection

Methods to detect signals of natural selection from genomic data have traditionally emphasized the use of simple summary statistics. Here, we review a new generation of methods that consider combinations of conventional summary statistics and/or richer features derived from inferred gene trees and ancestral recombination graphs (ARGs). We also review recent advances in methods for population genetic simulation and ARG reconstruction. Finally, we describe opportunities for future work on a variety of related topics, including the genetics of speciation, estimation of selection coefficients, and inference of selection on polygenic traits. Together, these emerging methods offer promising new directions in the study of natural selection.

The Role of Phylogenetically Conserved Elements in Shaping Patterns of Human Genomic Diversity

Evolutionary genetic studies have shown a positive correlation between levels of nucleotide diversity and either rates of recombination or genetic distance to genes. Both positive-directional and purifying selection have been offered as the source of these correlations via genetic hitchhiking and background selection, respectively. Phylogenetically conserved elements (CEs) are short (∼100 bp), widely distributed (comprising ∼5% of genome), sequences that are often found far from genes. While the function of many CEs is unknown, CEs also are associated with reduced diversity at linked sites. Using high coverage (>80×) whole genome data from two human populations, the Yoruba and the CEU, we perform fine scale evaluations of diversity, rates of recombination, and linkage to genes. We find that the local rate of recombination has a stronger effect on levels of diversity than linkage to genes, and that these effects of recombination persist even in regions far from genes. Our whole genome modeling demonstrates that, rather than recombination or GC-biased gene conversion, selection on sites within or linked to CEs better explains the observed genomic diversity patterns. A major implication is that very few sites in the human genome are predicted to be free of the effects of selection. These sites, which we refer to as the human "neutralome," comprise only 1.2% of the autosomes and 5.1% of the X chromosome. Demographic analysis of the neutralome reveals larger population sizes and lower rates of growth for ancestral human populations than inferred by previous analyses.

Recovering signals of ghost archaic introgression in African populations

While introgression from Neanderthals and Denisovans has been documented in modern humans outside Africa, the contribution of archaic hominins to the genetic variation of present-day Africans remains poorly understood. We provide complementary lines of evidence for archaic introgression into four West African populations. Our analyses of site frequency spectra indicate that these populations derive 2 to 19% of their genetic ancestry from an archaic population that diverged before the split of Neanderthals and modern humans. Using a method that can identify segments of archaic ancestry without the need for reference archaic genomes, we built genome-wide maps of archaic ancestry in the Yoruba and the Mende populations. Analyses of these maps reveal segments of archaic ancestry at high frequency in these populations that represent potential targets of adaptive introgression. Our results reveal the substantial contribution of archaic ancestry in shaping the gene pool of present-day West African populations.

Complex Phenotypes: Mechanisms Underlying Variation in Human Stature

PURPOSE OF REVIEW

The goal of the review is to provide a comprehensive overview of the current understanding of the mechanisms underlying variation in human stature.

RECENT FINDINGS

Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today, it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, the use of an evolutionary perspective has also revealed important mechanisms that have shaped height variation across the planet. This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

Searching for archaic contribution in Africa

Context: Africa's role in the narrative of human evolution is indisputably emphasised in the emergence of Homo sapiens. However, once humans dispersed beyond Africa, the history of those who stayed remains vastly under-studied, lacking the proper attention the birthplace of both modern and archaic humans deserves. The sequencing of Neanderthal and Denisovan genomes has elucidated evidence of admixture between archaic and modern humans outside of Africa, but has not aided efforts in answering whether archaic admixture happened within Africa. Objectives: This article reviews the state of research for archaic introgression in African populations and discusses recent insights into this topic. Methods: Gathering published sources and recently released preprints, this review reports on the different methods developed for detecting archaic introgression. Particularly it discusses how relevant these are when implemented on African populations and what findings these studies have shown so far. Results: Methods for detecting archaic introgression have been predominantly developed and implemented on non-African populations. Recent preprints present new methods considering African populations. While a number of studies using these methods suggest archaic introgression in Africa, without an African archaic genome to validate these results, such findings remain as putative archaic introgression. Conclusion: In light of the caveats with implementing current archaic introgression detection methods in Africa, we recommend future studies to concentrate on unravelling the complicated demographic history of Africa through means of ancient DNA where possible and through more focused efforts to sequence modern DNA from more representative populations across the African continent.

Whole-genome sequence analysis of a Pan African set of samples reveals archaic gene flow from an extinct basal population of modern humans into sub-Saharan populations

Background

Population demography and gene flow among African groups, as well as the putative archaic introgression of ancient hominins, have been poorly explored at the genome level.

Results

Here, we examine 15 African populations covering all major continental linguistic groups, ecosystems, and lifestyles within Africa through analysis of whole-genome sequence data of 21 individuals sequenced at deep coverage. We observe a remarkable correlation among genetic diversity and geographic distance, with the hunter-gatherer groups being more genetically differentiated and having larger effective population sizes throughout most modern-human history. Admixture signals are found between neighbor populations from both hunter-gatherer and agriculturalists groups, whereas North African individuals are closely related to Eurasian populations. Regarding archaic gene flow, we test six complex demographic models that consider recent admixture as well as archaic introgression. We identify the fingerprint of an archaic introgression event in the sub-Saharan populations included in the models (~ 4.0% in Khoisan, ~ 4.3% in Mbuti Pygmies, and ~ 5.8% in Mandenka) from an early divergent and currently extinct ghost modern human lineage.

Conclusion

The present study represents an in-depth genomic analysis of a Pan African set of individuals, which emphasizes their complex relationships and demographic history at population level.

Electronic supplementary material

The online version of this article (10.1186/s13059-019-1684-5) contains supplementary material, which is available to authorized users.

Clarifying distinct models of modern human origins in Africa

Accumulating genomic, fossil and archaeological data from Africa have led to a renewed interest in models of modern human origins. However, such discussions are often discipline-specific, with limited integration of evidence across the different fields. Further, geneticists typically require explicit specification of parameters to test competing demographic models, but these have been poorly outlined for some scenarios. Here, we describe four possible models for the origins of Homo sapiens in Africa based on published literature from paleoanthropology and human genetics. We briefly outline expectations for data patterns under each model, with a special focus on genetic data. Additionally, we present schematics for each model, doing our best to qualitatively describe demographic histories for which genetic parameters can be specifically attached. Finally, it is our hope that this perspective provides context for discussions of human origins in other manuscripts presented in this special issue.

Polygenic adaptation and convergent evolution on growth and cardiac genetic pathways in African and Asian rainforest hunter-gatherers

Different human populations facing similar environmental challenges have sometimes evolved convergent biological adaptations, for example, hypoxia resistance at high altitudes and depigmented skin in northern latitudes on separate continents. The "pygmy" phenotype (small adult body size), characteristic of hunter-gatherer populations inhabiting both African and Asian tropical rainforests, is often highlighted as another case of convergent adaptation in humans. However, the degree to which phenotypic convergence in this polygenic trait is due to convergent versus population-specific genetic changes is unknown. To address this question, we analyzed high-coverage sequence data from the protein-coding portion of the genomes of two pairs of populations: Batwa rainforest hunter-gatherers and neighboring Bakiga agriculturalists from Uganda and Andamanese rainforest hunter-gatherers and Brahmin agriculturalists from India. We observed signatures of convergent positive selection between the rainforest hunter-gatherers across the set of genes with "growth factor binding" functions ([Formula: see text]). Unexpectedly, for the rainforest groups, we also observed convergent and population-specific signatures of positive selection in pathways related to cardiac development (e.g., "cardiac muscle tissue development"; [Formula: see text]). We hypothesize that the growth hormone subresponsiveness likely underlying the adult small body-size phenotype may have led to compensatory changes in cardiac pathways, in which this hormone also plays an essential role. Importantly, in the agriculturalist populations, we did not observe similar patterns of positive selection on sets of genes associated with growth or cardiac development, indicating our results most likely reflect a history of convergent adaptation to the similar ecology of rainforests rather than a more general evolutionary pattern.

The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii)

Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers.

The demographic and adaptive history of central African hunter-gatherers and farmers

Central Africa, a forested region that supports an exceptionally high biodiversity, hosts the world's largest group of hunter-gatherers, who live in close proximity with groups that have adopted agriculture over the past 5000 years. Our understanding of the prehistory of these populations has been dramatically hampered by the almost total absence of fossil remains in this region, a limitation that has recently been circumvented by population genomics approaches. Different studies have estimated that ancestors of rainforest hunter-gatherers and Bantu-speaking farmers separated more than 60 000 years ago, supporting the occurrence of ancient population structure in Africa since the Late Pleistocene. Conversely, the Holocene in central Africa was characterized by large-scale population migrations associated with the emergence of agriculture, and increased genetic interactions between autochthonous rainforest hunter-gatherers and expanding Bantu-speaking farmers. Genomic scans have detected numerous candidate loci for positive selection in these populations, including convergent adaptation for short stature in groups of rainforest hunter-gatherers and local adaptation to endemic malaria in western and central Africans. Furthermore, there is recent increasing evidence that adaptive variation has been acquired by various African populations through admixture, suggesting a previously unappreciated role of intraspecies gene flow in local adaptation. Ancient and modern DNA studies will greatly broaden, and probably challenge, our view on the past history of central Africa, where introgression from yet uncharacterized archaic hominins and long-term adaptation to distinct ecological niches are suspected.

African genetic diversity provides novel insights into evolutionary history and local adaptations

Genetic variation and susceptibility to disease are shaped by human demographic history and adaptation. We can now study the genomes of extant Africans and uncover traces of population migration, admixture, assimilation and selection by applying sophisticated computational algorithms. There are four major ethnolinguistic divisions among present day Africans: Hunter-gatherer populations in southern and central Africa; Nilo-Saharan speakers from north and northeast Africa; Afro-Asiatic speakers from north and east Africa; and Niger-Congo speakers who are the predominant ethnolinguistic group spread across most of sub-Saharan Africa. The enormous ethnolinguistic diversity in sub-Saharan African populations is largely paralleled by extensive genetic diversity and until a decade ago, little was known about detailed origins and divergence of these groups. Results from large-scale population genetic studies, and more recently whole genome sequence data, are unravelling the critical role of events like migration and admixture and environmental factors including diet, infectious diseases and climatic conditions in shaping current population diversity. It is now possible to start providing quantitative estimates of divergence times, population size and dynamic processes that have affected populations and their genetic risk for disease. Finally, the availability of ancient genomes from Africa provides historical insights of unprecedented depth. In this review, we highlight some key interpretations that have emerged from recent African genome studies.

A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers

The evolutionary and biological bases of the Central African "pygmy" phenotype, a characteristic of rainforest hunter-gatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations.

Exome Sequencing Provides Evidence of Polygenic Adaptation to a Fat-Rich Animal Diet in Indigenous Siberian Populations

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

Comparison of Single Genome and Allele Frequency Data Reveals Discordant Demographic Histories

Inference of demographic history from genetic data is a primary goal of population genetics of model and nonmodel organisms. Whole genome-based approaches such as the pairwise/multiple sequentially Markovian coalescent methods use genomic data from one to four individuals to infer the demographic history of an entire population, while site frequency spectrum (SFS)-based methods use the distribution of allele frequencies in a sample to reconstruct the same historical events. Although both methods are extensively used in empirical studies and perform well on data simulated under simple models, there have been only limited comparisons of them in more complex and realistic settings. Here we use published demographic models based on data from three human populations (Yoruba, descendants of northwest-Europeans, and Han Chinese) as an empirical test case to study the behavior of both inference procedures. We find that several of the demographic histories inferred by the whole genome-based methods do not predict the genome-wide distribution of heterozygosity, nor do they predict the empirical SFS. However, using simulated data, we also find that the whole genome methods can reconstruct the complex demographic models inferred by SFS-based methods, suggesting that the discordant patterns of genetic variation are not attributable to a lack of statistical power, but may reflect unmodeled complexities in the underlying demography. More generally, our findings indicate that demographic inference from a small number of genomes, routine in genomic studies of nonmodel organisms, should be interpreted cautiously, as these models cannot recapitulate other summaries of the data.

Tracing the peopling of the world through genomics

Advances in the sequencing and the analysis of the genomes of both modern and ancient peoples have facilitated a number of breakthroughs in our understanding of human evolutionary history. These include the discovery of interbreeding between anatomically modern humans and extinct hominins; the development of an increasingly detailed description of the complex dispersal of modern humans out of Africa and their population expansion worldwide; and the characterization of many of the genetic adaptions of humans to local environmental conditions. Our interpretation of the evolutionary history and adaptation of humans is being transformed by analyses of these new genomic data.

Tracing the peopling of the world through genomics

Advances in the sequencing and the analysis of the genomes of both modern and ancient peoples have facilitated a number of breakthroughs in our understanding of human evolutionary history. These include the discovery of interbreeding between anatomically modern humans and extinct hominins; the development of an increasingly detailed description of the complex dispersal of modern humans out of Africa and their population expansion worldwide; and the characterization of many of the genetic adaptions of humans to local environmental conditions. Our interpretation of the evolutionary history and adaptation of humans is being transformed by analyses of these new genomic data.

Inferences of African evolutionary history from genomic data

Africa is the origin of anatomically modern humans and a continent of linguistic, cultural, environmental, phenotypic, and genetic diversity. However, African populations remain underrepresented in genetic studies, which have largely focused on individuals with European and Asian ancestry. The expansion of high-throughput 'omic' technologies to interrogate multiple tissue types across many biomolecules-DNA, proteins, epigenetic modifications, metabolites, and others-has heralded a new era of investigation into African history. In this review, we summarize how some of these recent advances have been applied to contemporary sub-Saharan African populations to inform studies on human origins and adaptation.

A genomic perspective on the generation and maintenance of genetic diversity in herbivorous insects

Understanding the processes that generate and maintain genetic variation within populations is a central goal in evolutionary biology. Theory predicts that some of this variation is maintained as a consequence of adapting to variable habitats. Studies in herbivorous insects have played a key role in confirming this prediction. Here, we highlight theoretical and conceptual models for the maintenance of genetic diversity in herbivorous insects, empirical genomic studies testing these models, and pressing questions within the realm of evolutionary and functional genomic studies. To address key gaps, we propose an integrative approach combining population genomic scans for adaptation, genome-wide characterization of targets of selection through experimental manipulations, mapping the genetic architecture of traits influencing fitness, and functional studies. We also stress the importance of studying the maintenance of genetic variation across biological scales-from variation within populations to divergence among populations-to form a comprehensive view of adaptation in herbivorous insects.

Hunter-gatherer genomics: evolutionary insights and ethical considerations

Hunting and gathering societies currently comprise only a small proportion of all human populations. However, the geographic and environmental diversity of modern hunter-gatherer groups, their inherent dependence on ecological resources, and their connection to patterns of behavior and subsistence that represent the vast majority of human history provide opportunities for scientific research to deliver major insights into the evolutionary history of our species. We review recent evolutionary genomic studies of hunter-gatherers, focusing especially on those that identify and functionally characterize phenotypic adaptations to local environments. We also call attention to specific ethical issues that scientists conducting hunter-gatherer genomics research ought to consider, including potential social and economic tensions between traditionally mobile hunter-gatherers and the land ownership-based nation-states by which they are governed, and the implications of genomic-based evidence of long-term evolutionary associations with particular habitats.

Model-based analyses of whole-genome data reveal a complex evolutionary history involving archaic introgression in Central African Pygmies

Comparisons of whole-genome sequences from ancient and contemporary samples have pointed to several instances of archaic admixture through interbreeding between the ancestors of modern non-Africans and now extinct hominids such as Neanderthals and Denisovans. One implication of these findings is that some adaptive features in contemporary humans may have entered the population via gene flow with archaic forms in Eurasia. Within Africa, fossil evidence suggests that anatomically modern humans (AMH) and various archaic forms coexisted for much of the last 200,000 yr; however, the absence of ancient DNA in Africa has limited our ability to make a direct comparison between archaic and modern human genomes. Here, we use statistical inference based on high coverage whole-genome data (greater than 60×) from contemporary African Pygmy hunter-gatherers as an alternative means to study the evolutionary history of the genus Homo. Using whole-genome simulations that consider demographic histories that include both isolation and gene flow with neighboring farming populations, our inference method rejects the hypothesis that the ancestors of AMH were genetically isolated in Africa, thus providing the first whole genome-level evidence of African archaic admixture. Our inferences also suggest a complex human evolutionary history in Africa, which involves at least a single admixture event from an unknown archaic population into the ancestors of AMH, likely within the last 30,000 yr.