Southern African ancient genomes estimate modern human divergence to 350,000 to 260,000 years ago

Deciphering African late middle Pleistocene hominin diversity and the origin of our species

The origin of Homo sapiens remains a matter of debate. The extent and geographic patterning of morphological diversity among Late Middle Pleistocene (LMP) African hominins is largely unknown, thus precluding the definition of boundaries of variability in early H. sapiens and the interpretation of individual fossils. Here we use a phylogenetic modelling method to predict possible morphologies of a last common ancestor of all modern humans, which we compare to LMP African fossils (KNM-ES 11693, Florisbad, Irhoud 1, Omo II, and LH18). Our results support a complex process for the evolution of H. sapiens, with the recognition of different, geographically localised, populations and lineages in Africa - not all of which contributed to our species' origin. Based on the available fossils, H. sapiens appears to have originated from the coalescence of South and, possibly, East-African source populations, while North-African fossils may represent a population which introgressed into Neandertals during the LMP.

A Rare Deep-Rooting D0 African Y-Chromosomal Haplogroup and Its Implications for the Expansion of Modern Humans Out of Africa

Present-day humans outside Africa descend mainly from a single expansion out ∼50,000-70,000 years ago, but many details of this expansion remain unclear, including the history of the male-specific Y chromosome at this time. Here, we reinvestigate a rare deep-rooting African Y-chromosomal lineage by sequencing the whole genomes of three Nigerian men described in 2003 as carrying haplogroup DE* Y chromosomes, and analyzing them in the context of a calibrated worldwide Y-chromosomal phylogeny. We confirm that these three chromosomes do represent a deep-rooting DE lineage, branching close to the DE bifurcation, but place them on the D branch as an outgroup to all other known D chromosomes, and designate the new lineage D0. We consider three models for the expansion of Y lineages out of Africa ∼50,000-100,000 years ago, incorporating migration back to Africa where necessary to explain present-day Y-lineage distributions. Considering both the Y-chromosomal phylogenetic structure incorporating the D0 lineage, and published evidence for modern humans outside Africa, the most favored model involves an origin of the DE lineage within Africa with D0 and E remaining there, and migration out of the three lineages (C, D, and FT) that now form the vast majority of non-African Y chromosomes. The exit took place 50,300-81,000 years ago (latest date for FT lineage expansion outside Africa - earliest date for the D/D0 lineage split inside Africa), and most likely 50,300-59,400 years ago (considering Neanderthal admixture). This work resolves a long-running debate about Y-chromosomal out-of-Africa/back-to-Africa migrations, and provides insights into the out-of-Africa expansion more generally.

The presence and impact of reference bias on population genomic studies of prehistoric human populations

Haploid high quality reference genomes are an important resource in genomic research projects. A consequence is that DNA fragments carrying the reference allele will be more likely to map successfully, or receive higher quality scores. This reference bias can have effects on downstream population genomic analysis when heterozygous sites are falsely considered homozygous for the reference allele. In palaeogenomic studies of human populations, mapping against the human reference genome is used to identify endogenous human sequences. Ancient DNA studies usually operate with low sequencing coverages and fragmentation of DNA molecules causes a large proportion of the sequenced fragments to be shorter than 50 bp-reducing the amount of accepted mismatches, and increasing the probability of multiple matching sites in the genome. These ancient DNA specific properties are potentially exacerbating the impact of reference bias on downstream analyses, especially since most studies of ancient human populations use pseudo-haploid data, i.e. they randomly sample only one sequencing read per site. We show that reference bias is pervasive in published ancient DNA sequence data of prehistoric humans with some differences between individual genomic regions. We illustrate that the strength of reference bias is negatively correlated with fragment length. Most genomic regions we investigated show little to no mapping bias but even a small proportion of sites with bias can impact analyses of those particular loci or slightly skew genome-wide estimates. Therefore, reference bias has the potential to cause minor but significant differences in the results of downstream analyses such as population allele sharing, heterozygosity estimates and estimates of archaic ancestry. These spurious results highlight how important it is to be aware of these technical artifacts and that we need strategies to mitigate the effect. Therefore, we suggest some post-mapping filtering strategies to resolve reference bias which help to reduce its impact substantially.

A catalog of single nucleotide changes distinguishing modern humans from archaic hominins

Throughout the past decade, studying ancient genomes has provided unique insights into human prehistory, and differences between modern humans and other branches like Neanderthals can enrich our understanding of the molecular basis of unique modern human traits. Modern human variation and the interactions between different hominin lineages are now well studied, making it reasonable to go beyond fixed genetic changes and explore changes that are observed at high frequency in present-day humans. Here, we identify 571 genes with non-synonymous changes at high frequency. We suggest that molecular mechanisms in cell division and networks affecting cellular features of neurons were prominently modified by these changes. Complex phenotypes in brain growth trajectory and cognitive traits are likely influenced by these networks and other non-coding changes presented here. We propose that at least some of these changes contributed to uniquely human traits, and should be prioritized for experimental validation.

Ancient DNA reveals a multistep spread of the first herders into sub-Saharan Africa

How food production first entered eastern Africa ~5000 years ago and the extent to which people moved with livestock is unclear. We present genome-wide data from 41 individuals associated with Later Stone Age, Pastoral Neolithic (PN), and Iron Age contexts in what are now Kenya and Tanzania to examine the genetic impacts of the spreads of herding and farming. Our results support a multiphase model in which admixture between northeastern African-related peoples and eastern African foragers formed multiple pastoralist groups, including a genetically homogeneous PN cluster. Additional admixture with northeastern and western African-related groups occurred by the Iron Age. These findings support several movements of food producers while rejecting models of minimal admixture with foragers and of genetic differentiation between makers of distinct PN artifacts.

Assessing the importance of cultural diffusion in the Bantu spread into southeastern Africa

The subsistence of Neolithic populations is based on agriculture, whereas that of previous populations was based on hunting and gathering. Neolithic spreads due to dispersal of populations are called demic, and those due to the incorporation of hunter-gatherers are called cultural. It is well-known that, after agriculture appeared in West Africa, it spread across most of subequatorial Africa. It has been proposed that this spread took place alongside with that of Bantu languages. In eastern and southeastern Africa, it is also linked to the Early Iron Age. From the beginning of the last millennium BC, cereal agriculture spread rapidly from the Great Lakes area eastwards to the East African coast, and southwards to northeastern South Africa. Here we show that the southwards spread took place substantially more rapidly (1.50–2.27 km/y) than the eastwards spread (0.59–1.27 km/y). Such a faster southwards spread could be the result of a stronger cultural effect. To assess this possibility, we compare these observed ranges to those obtained from a demic-cultural wave-of-advance model. We find that both spreads were driven by demic diffusion, in agreement with most archaeological, linguistic and genetic results. Nonetheless, the southwards spread seems to have indeed a stronger cultural component, which could lead support to the hypothesis that, at the southern areas, the interaction with pastoralist people may have played a significant role.

Dental evolutionary rates and its implications for the Neanderthal-modern human divergence

The origin of Neanderthal and modern human lineages is a matter of intense debate. DNA analyses have generally indicated that both lineages diverged during the middle period of the Middle Pleistocene, an inferred time that has strongly influenced interpretations of the hominin fossil record. This divergence time, however, is not compatible with the anatomical and genetic Neanderthal affinities observed in Middle Pleistocene hominins from Sima de los Huesos (Spain), which are dated to 430 thousand years (ka) ago. Drawing on quantitative analyses of dental evolutionary rates and Bayesian analyses of hominin phylogenetic relationships, I show that any divergence time between Neanderthals and modern humans younger than 800 ka ago would have entailed unexpectedly rapid dental evolution in early Neanderthals from Sima de los Huesos. These results support a pre-800 ka last common ancestor for Neanderthals and modern humans unless hitherto unexplained mechanisms sped up dental evolution in early Neanderthals.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations

BACKGROUND

Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations.

RESULTS

Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations.

CONCLUSIONS

We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.

Out of Africa by spontaneous migration waves

Hominin evolution is characterized by progressive regional differentiation, as well as migration waves, leading to anatomically modern humans that are assumed to have emerged in Africa and spread over the whole world. Why or whether Africa was the source region of modern humans and what caused their spread remains subject of ongoing debate. We present a spatially explicit, stochastic numerical model that includes ongoing mutations, demic diffusion, assortative mating and migration waves. Diffusion and assortative mating alone result in a structured population with relatively homogeneous regions bound by sharp clines. The addition of migration waves results in a power-law distribution of wave areas: for every large wave, many more small waves are expected to occur. This suggests that one or more out-of-Africa migrations would probably have been accompanied by numerous smaller migration waves across the world. The migration waves are considered "spontaneous", as the current model excludes environmental or other extrinsic factors. Large waves preferentially emanate from the central areas of large, compact inhabited areas. During the Pleistocene, Africa was the largest such area most of the time, making Africa the statistically most likely origin of anatomically modern humans, without a need to invoke additional environmental or ecological drivers.

Enamel pearls: Their occurrence in recent human populations and earliest manifestation in the modern human lineage

OBJECTIVE

To document and describe the occurrence of an enamel pearl on the distal root surface of the maxillary M3 of the fossil hominin specimen from Florisbad, South Africa that is dated to ca. 259,000 years B.P., and is an early representative of Homo sapiens or as a member of the evolutionary line that was directly ancestral to modern humans.

DESIGN

The molar was examined macroscopically and by micro-computed tomography (μCT) to enable accurate measurement and visualization of the structure of the enamel pearl.

RESULTS

The single pearl has a diameter of 0.97 mm; it is a Type 2 "composite" pearl comprising an enamel cap and dentine core without pulp chamber involvement. The size of the Florisbad pearl falls within or just below the size ranges of this anomaly in modern human samples. Type 2 pearls are most commonly encountered in recent human populations, and the location of the pearl on the distal root surface of the Florisbad M³ is consistent with its most frequent location in recent humans. Pearls in recent human populations affect between 0.2-4.8% of individuals, and 1.7-6.8% of permanent molars. Pearls have been documented in several prehistoric human dentitions, and all examples are less than 4000 years old.

CONCLUSIONS

Enamel pearls have been associated with periodontal disease, but it is not possible to relate its presence to the advanced periodontal inflammation and alveolar bone loss in the Florisbad fossil. Florisbad presents the earliest evidence of this anomaly in the fossil record pertaining to modern humans.

A dispersal of Homo sapiens from southern to eastern Africa immediately preceded the out-of-Africa migration

Africa was the birth-place of Homo sapiens and has the earliest evidence for symbolic behaviour and complex technologies. The best-attested early flowering of these distinctive features was in a glacial refuge zone on the southern coast 100–70 ka, with fewer indications in eastern Africa until after 70 ka. Yet it was eastern Africa, not the south, that witnessed the first major demographic expansion, ~70–60 ka, which led to the peopling of the rest of the world. One possible explanation is that important cultural traits were transmitted from south to east at this time. Here we identify a mitochondrial signal of such a dispersal soon after ~70 ka – the only time in the last 200,000 years that humid climate conditions encompassed southern and tropical Africa. This dispersal immediately preceded the out-of-Africa expansions, potentially providing the trigger for these expansions by transmitting significant cultural elements from the southern African refuge.

Genome-wide data from the Bubi of Bioko Island clarifies the Atlantic fringe of the Bantu dispersal

BACKGROUND

Bioko is one of the few islands that exist around Africa, the most genetically diverse continent on the planet. The native Bantu-speaking inhabitants of Bioko, the Bubi, are believed to have colonized the island about 2000 years ago. Here, we sequenced the genome of thirteen Bubi individuals at high coverage and analysed their sequences in comparison to mainland populations from the Gulf of Guinea.

RESULTS

We found that, genetically, the closest mainland population to the Bubi are Bantu-speaking groups from Angola instead the geographically closer groups from Cameroon. The Bubi possess a lower proportion of rainforest hunter-gatherer (RHG) ancestry than most other Bantu-speaking groups. However, their RHG component most likely came from the same source and could have reached them by gene flow from the mainland after island settlement. By studying identity by descent (IBD) genomic blocks and runs of homozygosity (ROHs), we found evidence for a significant level of genetic isolation among the Bubi, isolation that can be attributed to the island effect. Additionally, as this population is known to have one of the highest malaria incidence rates in the world we analysed their genome for malaria-resistant alleles. However, we were unable to detect any specific selective sweeps related to this disease.

CONCLUSIONS

By describing their dispersal to the Atlantic islands, the genomic characterization of the Bubi contributes to the understanding of the margins of the massive Bantu migration that shaped all Sub-Saharan African populations.

Ancestral mitochondrial N lineage from the Neolithic 'green' Sahara

Because Africa's climate hampers DNA preservation, knowledge of its genetic variability is mainly restricted to modern samples, even though population genetics dynamics and back-migrations from Eurasia may have modified haplotype frequencies, masking ancient genetic scenarios. Thanks to improved methodologies, ancient genetic data for the African continent are now increasingly available, starting to fill in the gap. Here we present newly obtained mitochondrial genomes from two ~7000-year-old individuals from Takarkori rockshelter, Libya, representing the earliest and first genetic data for the Sahara region. These individuals carry a novel mutation motif linked to the haplogroup N root. Our result demonstrates the presence of an ancestral lineage of the N haplogroup in the Holocene "Green Sahara", associated to a Middle Pastoral (Neolithic) context.

The people behind the samples: Biographical features of Past Hunter-Gatherers from KwaZulu-Natal who yielded aDNA

PURPOSE

Skeletons sampled for ancient human DNA analysis are sometimes complete enough to provide information about the lives of the people they represent. We focus on three Later Stone Age skeletons, ca. 2000 B.P., from coastal KwaZulu-Natal, South Africa, whose ancient genomes have been sequenced (Schlebusch et al., 2017).

METHODS

Bioarchaeological approaches are integrated with aDNA information.

RESULTS

All skeletons are male. Dental development shows that the boy, with prominent cribra orbitalia, died at age 6-7 years. Two men show cranial and spinal trauma, extensive tooth wear, plus mild cribra orbitalia in one.

CONCLUSIONS

Dental wear and trauma of the adults are consistent with hunter-gatherer lives. Even partial aDNA evidence contributes to sex determination. Parasitic infection such as schistosomiasis is the best-fit cause for the child's anemia in this case.

CONTRIBUTION TO KNOWLEDGE

The convergence of genomic and bioarchaeological approaches expands our knowledge of the past lives of a boy and two men whose lives as hunter-gatherers included episodes of trauma and disease.

LIMITATIONS

The skeletons are incomplete, in variable condition, and from poorly characterized local cultural contexts.

SUGGESTIONS FOR FURTHER RESEARCH

Thorough osteobiographic analysis should accompany paleogenomic investigations. Such disciplinary collaboration enriches our understanding of the human past.

The Emergence of Genomic Research in Africa and New Frameworks for Equity in Biomedical Research

Individuals with African ancestry have the greatest genomic diversity in the world, yet they have been underrepresented in genomic research. To advance our understanding of human biology and our ability to trace human history, we must include more samples from Africans in genomic research. Additionally, inclusion of more samples from participants of recent African descent is imperative to provide equitable health care as genomics is increasingly used for diagnosis, treatment, and to understand disease risk. The Human Heredity and Health in Africa initiative (H3Africa) seeks to expand the number of Africans included in genomic research and to do so by expanding the research capacity on the continent. In this article, we discuss how H3Africa is endeavoring to achieve these goals while promoting equitable research collaborations.

Limits of long-term selection against Neandertal introgression

Several studies have suggested that introgressed Neandertal DNA was subjected to negative selection in modern humans. A striking observation in support of this is an apparent monotonic decline in Neandertal ancestry observed in modern humans in Europe over the past 45,000 years. Here, we show that this decline is an artifact likely caused by gene flow between modern human populations, which is not taken into account by statistics previously used to estimate Neandertal ancestry. When we apply a statistic that avoids assumptions about modern human demography by taking advantage of two high-coverage Neandertal genomes, we find no evidence for a change in Neandertal ancestry in Europe over the past 45,000 years. We use whole-genome simulations of selection and introgression to investigate a wide range of model parameters and find that negative selection is not expected to cause a significant long-term decline in genome-wide Neandertal ancestry. Nevertheless, these models recapitulate previously observed signals of selection against Neandertal alleles, in particular the depletion of Neandertal ancestry in conserved genomic regions. Surprisingly, we find that this depletion is strongest in regulatory and conserved noncoding regions and in the most conserved portion of protein-coding sequences.

Approximate Bayesian computation with deep learning supports a third archaic introgression in Asia and Oceania

Since anatomically modern humans dispersed Out of Africa, the evolutionary history of Eurasian populations has been marked by introgressions from presently extinct hominins. Some of these introgressions have been identified using sequenced ancient genomes (Neanderthal and Denisova). Other introgressions have been proposed for still unidentified groups using the genetic diversity present in current human populations. We built a demographic model based on deep learning in an Approximate Bayesian Computation framework to infer the evolutionary history of Eurasian populations including past introgression events in Out of Africa populations fitting the current genetic evidence. In addition to the reported Neanderthal and Denisovan introgressions, our results support a third introgression in all Asian and Oceanian populations from an archaic population. This population is either related to the Neanderthal-Denisova clade or diverged early from the Denisova lineage. We propose the use of deep learning methods for clarifying situations with high complexity in evolutionary genomics.

Rapid evolution of a skin-lightening allele in southern African KhoeSan

Skin pigmentation is under strong directional selection in northern European and Asian populations. The indigenous KhoeSan populations of far southern Africa have lighter skin than other sub-Saharan African populations, potentially reflecting local adaptation to a region of Africa with reduced UV radiation. Here, we demonstrate that a canonical Eurasian skin pigmentation gene, SLC24A5, was introduced to southern Africa via recent migration and experienced strong adaptive evolution in the KhoeSan. To reconstruct the evolution of skin pigmentation, we collected phenotypes from over 400 ≠Khomani San and Nama individuals and high-throughput sequenced candidate pigmentation genes. The derived causal allele in SLC24A5, p.Ala111Thr, significantly lightens basal skin pigmentation in the KhoeSan and explains 8 to 15% of phenotypic variance in these populations. The frequency of this allele (33 to 53%) is far greater than expected from colonial period European gene flow; however, the most common derived haplotype is identical among European, eastern African, and KhoeSan individuals. Using four-population demographic simulations with selection, we show that the allele was introduced into the KhoeSan only 2,000 y ago via a back-to-Africa migration and then experienced a selective sweep (s = 0.04 to 0.05 in ≠Khomani and Nama). The SLC24A5 locus is both a rare example of intense, ongoing adaptation in very recent human history, as well as an adaptive gene flow at a pigmentation locus in humans.

Pliocene and Eocene provide best analogs for near-future climates

As the world warms due to rising greenhouse gas concentrations, the Earth system moves toward climate states without societal precedent, challenging adaptation. Past Earth system states offer possible model systems for the warming world of the coming decades. These include the climate states of the Early Eocene (ca. 50 Ma), the Mid-Pliocene (3.3-3.0 Ma), the Last Interglacial (129-116 ka), the Mid-Holocene (6 ka), preindustrial (ca. 1850 CE), and the 20th century. Here, we quantitatively assess the similarity of future projected climate states to these six geohistorical benchmarks using simulations from the Hadley Centre Coupled Model Version 3 (HadCM3), the Goddard Institute for Space Studies Model E2-R (GISS), and the Community Climate System Model, Versions 3 and 4 (CCSM) Earth system models. Under the Representative Concentration Pathway 8.5 (RCP8.5) emission scenario, by 2030 CE, future climates most closely resemble Mid-Pliocene climates, and by 2150 CE, they most closely resemble Eocene climates. Under RCP4.5, climate stabilizes at Pliocene-like conditions by 2040 CE. Pliocene-like and Eocene-like climates emerge first in continental interiors and then expand outward. Geologically novel climates are uncommon in RCP4.5 (<1%) but reach 8.7% of the globe under RCP8.5, characterized by high temperatures and precipitation. Hence, RCP4.5 is roughly equivalent to stabilizing at Pliocene-like climates, while unmitigated emission trajectories, such as RCP8.5, are similar to reversing millions of years of long-term cooling on the scale of a few human generations. Both the emergence of geologically novel climates and the rapid reversion to Eocene-like climates may be outside the range of evolutionary adaptive capacity.

The evolutionary history of Southern Africa

The genomic variability of Southern African groups is characterized by an exceptional degree of diversity, which is the result of long-term local evolutionary history, migrations and gene-flow. Over the last few years several investigations have identified and described signatures related to these processes, revealing how ancient and more recent events have shaped the structure and ancestry composition of local populations. Here we discuss recent insights into the genetic history of the Southernmost part of the African continent provided by the analysis of modern and ancient genomes. Future work is expected to clarify the population dynamics associated with the emergence of Homo sapiens across Africa and the details of the process of dispersion and admixture associated with the arrival of Bantu-speaking groups in the region.

FADS1 and the Timing of Human Adaptation to Agriculture

Variation at the FADS1/FADS2 gene cluster is functionally associated with differences in lipid metabolism and is often hypothesized to reflect adaptation to an agricultural diet. Here, we test the evidence for this relationship using both modern and ancient DNA data. We show that almost all the inhabitants of Europe carried the ancestral allele until the derived allele was introduced ∼8,500 years ago by Early Neolithic farming populations. However, we also show that it was not under strong selection in these populations. We find that this allele, and other proposed agricultural adaptations at LCT/MCM6 and SLC22A4, were not strongly selected until much later, perhaps as late as the Bronze Age. Similarly, increased copy number variation at the salivary amylase gene AMY1 is not linked to the development of agriculture although, in this case, the putative adaptation precedes the agricultural transition. Our analysis shows that selection at the FADS locus was not tightly linked to the initial introduction of agriculture and the Neolithic transition. Further, it suggests that the strongest signals of recent human adaptation in Europe did not coincide with the Neolithic transition but with more recent changes in environment, diet, or efficiency of selection due to increases in effective population size.

The role of matrilineality in shaping patterns of Y chromosome and mtDNA sequence variation in southwestern Angola

Southwestern Angola is a region characterized by contact between indigenous foragers and incoming food-producers, involving genetic and cultural exchanges between peoples speaking Kx'a, Khoe-Kwadi, and Bantu languages. Although present-day Bantu speakers share a patrilocal residence pattern and matrilineal principle of clan and group membership, a highly stratified social setting divides dominant pastoralists from marginalized groups that subsist on alternative strategies and have previously been thought to have pre-Bantu origins. Here, we compare new high-resolution sequence data from 2.3 Mb of the male-specific region of the Y chromosome (MSY) from 170 individuals with previously reported mitochondrial DNA (mtDNA) genomes, to investigate the population history of seven representative southwestern Angolan groups (Himba, Kuvale, Kwisi, Kwepe, Twa, Tjimba, !Xun), and to study the causes and consequences of sex-biased processes in their genetic variation. We found no clear link between the formerly Kwadi-speaking Kwepe and pre-Bantu eastern African migrants, and no pre-Bantu MSY lineages among Bantu-speaking groups, except for small amounts of "Khoisan" introgression. We therefore propose that irrespective of their subsistence strategies, all Bantu-speaking groups of the area share a male Bantu origin. Additionally, we show that in Bantu-speaking groups, the levels of among-group and between-group variation are higher for mtDNA than for MSY. These results, together with our previous demonstration that the matriclanic systems of southwestern Angolan Bantu groups are genealogically consistent, suggest that matrilineality strongly enhances both female population sizes and interpopulation mtDNA variation.

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.

Origin and age of the causative mutations in KLC2, IMPA1, MED25 and WNT7A unravelled through Brazilian admixed populations

The mutation age and local ancestry of chromosomal segments harbouring mutations associated with autosomal recessive (AR) disorders in Brazilian admixed populations remain unknown; additionally, inbreeding levels for these affected individuals continue to be estimated based on genealogical information. Here, we calculated inbreeding levels using a runs of homozygosity approach, mutation age and local ancestry to infer the origin of each chromosomal segments containing disorder-causing mutations in KLC2, IMPA1, MED25 and WNT7A. Genotyped data were generated from 18 patients affected by AR diseases and combined to the 1000 genome project (1KGP) and Simons genome diversity project (SGDP) databases to infer local ancestry. We found a major European contribution for mutated haplotypes with recent mutation age and inbreeding values found only in Native American and Middle East individuals. These results contribute to identifying the origin of and to understanding how these diseases are maintained and spread in Brazilian and world populations.

Quantifying and reducing spurious alignments for the analysis of ultra-short ancient DNA sequences

BACKGROUND

The study of ancient DNA is hampered by degradation, resulting in short DNA fragments. Advances in laboratory methods have made it possible to retrieve short DNA fragments, thereby improving access to DNA preserved in highly degraded, ancient material. However, such material contains large amounts of microbial contamination in addition to DNA fragments from the ancient organism. The resulting mixture of sequences constitutes a challenge for computational analysis, since microbial sequences are hard to distinguish from the ancient sequences of interest, especially when they are short.

RESULTS

Here, we develop a method to quantify spurious alignments based on the presence or absence of rare variants. We find that spurious alignments are enriched for mismatches and insertion/deletion differences and lack substitution patterns typical of ancient DNA. The impact of spurious alignments can be reduced by filtering on these features and by imposing a sample-specific minimum length cutoff. We apply this approach to sequences from four ~ 430,000-year-old Sima de los Huesos hominin remains, which contain particularly short DNA fragments, and increase the amount of usable sequence data by 17-150%. This allows us to place a third specimen from the site on the Neandertal lineage.

CONCLUSIONS

Our method maximizes the sequence data amenable to genetic analysis from highly degraded ancient material and avoids pitfalls that are associated with the analysis of ultra-short DNA sequences.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.

Selecting among Alternative Scenarios of Human Evolution by Simulated Genetic Gradients

Selecting among alternative scenarios of human evolution is nowadays a common methodology to investigate the history of our species. This strategy is usually based on computer simulations of genetic data under different evolutionary scenarios, followed by a fitting of the simulated data with the real data. A recent trend in the investigation of ancestral evolutionary processes of modern humans is the application of genetic gradients as a measure of fitting, since evolutionary processes such as range expansions, range contractions, and population admixture (among others) can lead to different genetic gradients. In addition, this strategy allows the analysis of the genetic causes of the observed genetic gradients. Here, we review recent findings on the selection among alternative scenarios of human evolution based on simulated genetic gradients, including pros and cons. First, we describe common methodologies to simulate genetic gradients and apply them to select among alternative scenarios of human evolution. Next, we review previous studies on the influence of range expansions, population admixture, last glacial period, and migration with long-distance dispersal on genetic gradients for some regions of the world. Finally, we discuss this analytical approach, including technical limitations, required improvements, and advice. Although here we focus on human evolution, this approach could be extended to study other species.

The Combinatorial Creature: Cortical Phenotypes within and across Lifetimes

The neocortex is one of the most distinctive structures of the mammalian brain, yet also one of the most varied in terms of both size and organization. Multiple processes have contributed to this variability, including evolutionary mechanisms (i.e., alterations in gene sequence) that alter the size, organization, and connections of neocortex, and activity dependent mechanisms that can also modify these same features. Thus, changes to the neocortex can occur over different time-scales, including within a single generation. This combination of genetic and activity dependent mechanisms that create a given cortical phenotype allows the mammalian neocortex to rapidly and flexibly adjust to different body and environmental contexts, and in humans permits culture to impact brain construction.

Genetic structure and sex-biased gene flow in the history of southern African populations

Objectives

We investigated the genetic history of southern African populations with a special focus on their paternal history. We reexamined previous claims that the Y‐chromosome haplogroup E1b1b (E‐M293) was brought to southern Africa by pastoralists from eastern Africa, and investigated patterns of sex‐biased gene flow in southern Africa.

Materials and methods

We analyzed previously published complete mtDNA genome sequences and ∼900 kb of NRY sequences from 23 populations from Namibia, Botswana, and Zambia, as well as haplogroup frequencies from a large sample of southern African populations and 23 newly genotyped Y‐linked STR loci for samples assigned to haplogroup E1b1b.

Results

Our results support an eastern African origin for Y‐chromosome haplogroup E1b1b (E‐M293); however, its current distribution in southern Africa is not strongly associated with pastoralism, suggesting more complex demographic events and/or changes in subsistence practices in this region. The Bantu expansion in southern Africa had a notable genetic impact and was probably a rapid, male‐dominated expansion. Our finding of a significant increase in the intensity of the sex‐biased gene flow from north to south may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

Conclusions

Our study shows that the population history of southern Africa has been complex, with different immigrating groups mixing to different degrees with the autochthonous populations. The Bantu expansion led to heavily sex‐biased admixture as a result of interactions between Khoisan females and Bantu males, with a geographic gradient which may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

Ancient Genomics of Modern Humans: The First Decade

The first decade of ancient genomics has revolutionized the study of human prehistory and evolution. We review new insights based on prehistoric modern human genomes, including greatly increased resolution of the timing and structure of the out-of-Africa expansion, the diversification of present-day non-African populations, and the earliest expansions of those populations into Eurasia and America. Prehistoric genomes now document population transformations on every inhabited continent—in particular the effect of agricultural expansions in Africa, Europe, and Oceania—and record a history of natural selection that shapes present-day phenotypic diversity. Despite these advances, much remains unknown, in particular about the genomic histories of Asia (the most populous continent) and Africa (the continent that contains the most genetic diversity). Ancient genomes from these and other regions, integrated with a growing understanding of the genomic basis of human phenotypic diversity, will be in focus during the next decade of research in the field.

Tales of Human Migration, Admixture, and Selection in Africa

In the last three decades, genetic studies have played an increasingly important role in exploring human history. They have helped to conclusively establish that anatomically modern humans first appeared in Africa roughly 250,000–350,000 years before present and subsequently migrated to other parts of the world. The history of humans in Africa is complex and includes demographic events that influenced patterns of genetic variation across the continent. Through genetic studies, it has become evident that deep African population history is captured by relationships among African hunter–gatherers, as the world's deepest population divergences occur among these groups, and that the deepest population divergence dates to 300,000 years before present. However, the spread of pastoralism and agriculture in the last few thousand years has shaped the geographic distribution of present-day Africans and their genetic diversity. With today's sequencing technologies, we can obtain full genome sequences from diverse sets of extant and prehistoric Africans. The coming years will contribute exciting new insights toward deciphering human evolutionary history in Africa.

In-solution Y-chromosome capture-enrichment on ancient DNA libraries

Background

As most ancient biological samples have low levels of endogenous DNA, it is advantageous to enrich for specific genomic regions prior to sequencing. One approach—in-solution capture-enrichment—retrieves sequences of interest and reduces the fraction of microbial DNA. In this work, we implement a capture-enrichment approach targeting informative regions of the Y chromosome in six human archaeological remains excavated in the Caribbean and dated between 200 and 3000 years BP. We compare the recovery rate of Y-chromosome capture (YCC) alone, whole-genome capture followed by YCC (WGC + YCC) versus non-enriched (pre-capture) libraries.

Results

The six samples show different levels of initial endogenous content, with very low (< 0.05%, 4 samples) or low (0.1–1.54%, 2 samples) percentages of sequenced reads mapping to the human genome. We recover 12–9549 times more targeted unique Y-chromosome sequences after capture, where 0.0–6.2% (WGC + YCC) and 0.0–23.5% (YCC) of the sequence reads were on-target, compared to 0.0–0.00003% pre-capture. In samples with endogenous DNA content greater than 0.1%, we found that WGC followed by YCC (WGC + YCC) yields lower enrichment due to the loss of complexity in consecutive capture experiments, whereas in samples with lower endogenous content, the libraries’ initial low complexity leads to minor proportions of Y-chromosome reads. Finally, increasing recovery of informative sites enabled us to assign Y-chromosome haplogroups to some of the archeological remains and gain insights about their paternal lineages and origins.

Conclusions

We present to our knowledge the first in-solution capture-enrichment method targeting the human Y-chromosome in aDNA sequencing libraries. YCC and WGC + YCC enrichments lead to an increase in the amount of Y-DNA sequences, as compared to libraries not enriched for the Y-chromosome. Our probe design effectively recovers regions of the Y-chromosome bearing phylogenetically informative sites, allowing us to identify paternal lineages with less sequencing than needed for pre-capture libraries. Finally, we recommend considering the endogenous content in the experimental design and avoiding consecutive rounds of capture, as clonality increases considerably with each round.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4945-x) contains supplementary material, which is available to authorized users.

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.

Middle Stone Age human teeth from Magubike rockshelter, Iringa Region, Tanzania

In 2006, six isolated hominin teeth were excavated from Middle Stone Age (MSA) deposits at the Magubike rockshelter in southern Tanzania. They comprise two central incisors, one lateral incisor, one canine, one third premolar, and one fourth premolar. All are fully developed and come from the maxilla. None of the teeth are duplicated, so they may represent a single individual. While there is some evidence of post-depositional alteration, the morphology of these teeth clearly shares features with anatomically modern Homo sapiens. Both metric and non-metric traits are compared to those from other African and non-African dental remains. The degree of biological relatedness between eastern and southern African Stone Age hunter-gatherers has long been a subject of interest, and several characteristics of the Magubike teeth resemble those of the San of southern Africa. Another notable feature is that the three incisors are marked on the labial crown by scratches that are much coarser than microwear striations. These non-masticatory scratches on the Magubike teeth suggest that the use of the front teeth as tools included regularly repeated activities undertaken throughout the life of the individual. The exact age of these teeth is not clear as ESR and radiocarbon dates on associated snail shells give varying results, but a conservative estimate of their minimum age is 45,000 years.

Did Our Species Evolve in Subdivided Populations across Africa, and Why Does It Matter?

We challenge the view that our species, Homo sapiens, evolved within a single population and/or region of Africa. The chronology and physical diversity of Pleistocene human fossils suggest that morphologically varied populations pertaining to the H. sapiens clade lived throughout Africa. Similarly, the African archaeological record demonstrates the polycentric origin and persistence of regionally distinct Pleistocene material culture in a variety of paleoecological settings. Genetic studies also indicate that present-day population structure within Africa extends to deep times, paralleling a paleoenvironmental record of shifting and fractured habitable zones. We argue that these fields support an emerging view of a highly structured African prehistory that should be considered in human evolutionary inferences, prompting new interpretations, questions, and interdisciplinary research directions.

The view that Homo sapiens evolved from a single region/population within Africa has been given primacy in studies of human evolution.

However, developments across multiple fields show that relevant data are no longer consistent with this view.

We argue instead that Homo sapiens evolved within a set of interlinked groups living across Africa, whose connectivity changed through time.

Genetic models therefore need to incorporate a more complex view of ancient migration and divergence in Africa.

We summarize this new framework emphasizing population structure, outline how this changes our understanding of human evolution, and identify new research directions.

Who were the Nataruk people? Mandibular morphology among late Pleistocene and early Holocene fisher-forager populations of West Turkana (Kenya)

Africa is the birthplace of the species Homo sapiens, and Africans today are genetically more diverse than other populations of the world. However, the processes that underpinned the evolution of African populations remain largely obscure. Only a handful of late Pleistocene African fossils (∼50-12 Ka) are known, while the more numerous sites with human fossils of early Holocene age are patchily distributed. In particular, late Pleistocene and early Holocene human diversity in Eastern Africa remains little studied, precluding any analysis of the potential factors that shaped human diversity in the region, and more broadly throughout the continent. These periods include the Last Glacial Maximum (LGM), a moment of extreme aridity in Africa that caused the fragmentation of population ranges and localised extinctions, as well as the 'African Humid Period', a moment of abrupt climate change and enhanced connectivity throughout Africa. East Africa, with its range of environments, may have acted as a refugium during the LGM, and may have played a critical biogeographic role during the heterogene`ous environmental recovery that followed. This environmental context raises a number of questions about the relationships among early Holocene African populations, and about the role played by East Africa in shaping late hunter-gatherer biological diversity. Here, we describe eight mandibles from Nataruk, an early Holocene site (∼10 Ka) in West Turkana, offering the opportunity of exploring population diversity in Africa at the height of the 'African Humid Period'. We use 3D geometric morphometric techniques to analyze the phenotypic variation of a large mandibular sample. Our results show that (i) the Nataruk mandibles are most similar to other African hunter-fisher-gatherer populations, especially to the fossils from Lothagam, another West Turkana locality, and to other early Holocene fossils from the Central Rift Valley (Kenya); and (ii) a phylogenetic connection may have existed between these Eastern African populations and some Nile Valley and Maghrebian groups, who lived at a time when a Green Sahara may have allowed substantial contact, and potential gene flow, across a vast expanse of Northern and Eastern Africa.

Cultural hitchhiking and competition between patrilineal kin groups explain the post-Neolithic Y-chromosome bottleneck

In human populations, changes in genetic variation are driven not only by genetic processes, but can also arise from cultural or social changes. An abrupt population bottleneck specific to human males has been inferred across several Old World (Africa, Europe, Asia) populations 5000–7000 BP. Here, bringing together anthropological theory, recent population genomic studies and mathematical models, we propose a sociocultural hypothesis, involving the formation of patrilineal kin groups and intergroup competition among these groups. Our analysis shows that this sociocultural hypothesis can explain the inference of a population bottleneck. We also show that our hypothesis is consistent with current findings from the archaeogenetics of Old World Eurasia, and is important for conceptions of cultural and social evolution in prehistory.

A population bottleneck 5000-7000 years ago in human males, but not females, has been inferred across several African, European and Asian populations. Here, Zeng and colleagues synthesize theory and data to suggest that competition among patrilineal kin groups produced the bottleneck pattern.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2⁺ HLA-C. The relatively high frequency of C2⁺ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

Pleistocene North African genomes link Near Eastern and sub-Saharan African human populations

North Africa is a key region for understanding human history, but the genetic history of its people is largely unknown. We present genomic data from seven 15,000-year-old modern humans, attributed to the Iberomaurusian culture, from Morocco. We find a genetic affinity with early Holocene Near Easterners, best represented by Levantine Natufians, suggesting a pre-agricultural connection between Africa and the Near East. We do not find evidence for gene flow from Paleolithic Europeans to Late Pleistocene North Africans. The Taforalt individuals derive one-third of their ancestry from sub-Saharan Africans, best approximated by a mixture of genetic components preserved in present-day West and East Africans. Thus, we provide direct evidence for genetic interactions between modern humans across Africa and Eurasia in the Pleistocene.

Biological Sexing of a 4000-Year-Old Egyptian Mummy Head to Assess the Potential of Nuclear DNA Recovery from the Most Damaged and Limited Forensic Specimens

High throughput sequencing (HTS) has been used for a number of years in the field of paleogenomics to facilitate the recovery of small DNA fragments from ancient specimens. Recently, these techniques have also been applied in forensics, where they have been used for the recovery of mitochondrial DNA sequences from samples where traditional PCR-based assays fail because of the very short length of endogenous DNA molecules. Here, we describe the biological sexing of a ~4000-year-old Egyptian mummy using shotgun sequencing and two established methods of biological sex determination (R_X and R_Y), by way of mitochondrial genome analysis as a means of sequence data authentication. This particular case of historical interest increases the potential utility of HTS techniques for forensic purposes by demonstrating that data from the more discriminatory nuclear genome can be recovered from the most damaged specimens, even in cases where mitochondrial DNA cannot be recovered with current PCR-based forensic technologies. Although additional work remains to be done before nuclear DNA recovered via these methods can be used routinely in operational casework for individual identification purposes, these results indicate substantial promise for the retrieval of probative individually identifying DNA data from the most limited and degraded forensic specimens.