Khoe-San Genomes Reveal Unique Variation and Confirm the Deepest Population Divergence in Homo sapiens

Look past the cooperative eye hypothesis: reconsidering the evolution of human eye appearance

The external appearance of the human eye has been prominently linked to the evolution of complex sociocognitive functions in our species. The cooperative eye hypothesis (CEH) proposes that human eyeballs, with their weakly expressed conjunctival and scleral pigmentation, are uniquely conspicuous and evolved under selective pressures to behave cooperatively, therefore signalling attentiveness to conspecifics. Non-human primates are instead assumed to display less-salient eye morphologies that help mask their gaze to facilitate competitive, rather than cooperative actions. Here, we argue that the CEH, although continuing to be influential, lacks robust empirical support. Over the past two decades, multidisciplinary research has undermined its original rationale and central premises: human eye pigmentation does not uniquely stand out among primates, it is not uniform at species level and the available evidence does not conclusively suggest that it facilitates gaze following to notable extents. Hence, the CEH currently provides a theoretical framework that risks confusing, rather than informing, inferences about the evolution of human external eye appearance and its selective drivers. In a call to move past it, we review alternative hypotheses with the potential to elucidate the emergence of the human ocular phenotype from the considerable spectrum of diversity found within the primate order.

Uncharted territory: the role of mitochondrial DNA variation in macrophage-mediated host susceptibility to tuberculosis

Mitochondria form an integral, yet frequently underappreciated, part of the immune response to Mycobacterium tuberculosis (M.tb), particularly within macrophages. Despite growing recognition for their role in infection and immunity, studies investigating how mitochondrial DNA (mtDNA) variation influences host susceptibility to tuberculosis (TB) are limited. Notably, there are no studies in African-based populations, although Africans possess unparalleled human genetic diversity, including the earliest diverged mitochondrial haplogroups, and a high TB burden. This underrepresentation limits the discovery of novel ancestry-specific genetic loci associated with TB. In this review article, we describe the unique characteristics of mtDNA, highlight key mitochondrial functions relevant to macrophage responses during M.tb infection, and summarise published studies that investigate the role of host mtDNA variation in TB susceptibility. We further advocate for the inclusion of African populations in future studies to identify novel TB susceptibility genetic risk loci and expand the current knowledgebase on host TB susceptibility.

The Precision and Power of Population Branch Statistics in Identifying the Genomic Signatures of Local Adaptation

Population branch statistics, which estimate the degree of genetic differentiation along a focal population's lineage, have been used as an alternative to FST-based genome-wide scans for identifying loci associated with local selective sweeps. Beyond the population branch statistic (PBS), the normalized PBSn1 adjusts focal branch length with respect to outgroup branch lengths at the same locus, whereas population branch excess (PBE) incorporates median branch lengths at other loci. PBSn1 and PBE were proposed to be more specific to local selective sweeps as opposed to geographically ubiquitous selection. However, the accuracy and statistical power of branch statistics have not been systematically assessed. To do so, we simulate genomes in representative large and small populations with varying proportions of sites evolving under genetic drift or (approximated) background selection, with local selective sweeps or geographically parallel selective sweeps. We then assess the probability that local selective sweep loci are correctly identified as outliers by FST and by each of the branch statistics. We find that branch statistics consistently outperform FST at identifying local sweeps. Particularly when parallel sweeps are introduced, PBSn1 and PBE correctly identify local sweeps among their top outliers more frequently than PBS. Additionally, we evaluate versions of these statistics based on maximal site differentiation within a window, finding that site-based PBE and PBSn1 are particularly effective at identifying local soft sweeps. These results validate the greater specificity of the rescaled branch statistics PBE and PBSn1 to detect population-specific positive selection, supporting their use in genomic studies focused on local adaptation.

Population genomics of sika deer reveals recent speciation and genetic selective signatures during evolution and domestication

BACKGROUND

Population genomic analysis can reconstruct the phylogenetic relationship and demographic history, and identify genomic selective signatures of a species. To date, fundamental aspects of population genomic analyses, such as intraspecies taxonomy, evolutionary history, and adaptive evolution, of sika deer have not been systematically investigated. Furthermore, accumulating lines of evidences have illustrated that incorrect species delimitation will mislead conservation decisions, and even lead to irreversible mistakes in threatened species.

RESULTS

In this study, we resequenced 81 wild and 71 domesticated sika deer representing 10 main geographic populations and two farms to clarify the species delimitation, demographic and divergence histories, and adaptive evolution of this species. First, our analyses of whole genomes, Y chromosomes and mitochondrial genomes revealed substantial genetic differentiation between the continental and Japanese lineages of sika deer, representing two phylogenetically distinct species. Second, sika deer in Japan were inferred to have experienced a "divergence-mixing-isolation" evolutionary scenario. Third, we identified four candidate genes (XKR4, NPAS3, CTNNA3, and CNTNAP5) possibly involved in body size regulation of sika deer by selective sweep analysis. Furthermore, we also detected two candidate genes (NRP2 and EDIL3) that may be associated with an important economic trait (antler weight) were under selection during the process of domestication.

CONCLUSION

Population genomic analyses revealed that the continental and Japanese lineages represent distinct phylogenetic species. Moreover, our results provide insights into the genetic selection signatures related to body size differences and a valuable genomic resource for future genetic studies and genomics-informed breeding of sika deer.

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

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

Linguistic capacity was present in the Homo sapiens population 135 thousand years ago

Recent genome-level studies on the divergence of early Homo sapiens, based on single nucleotide polymorphisms, suggest that the initial population division within H. sapiens from the original stem occurred approximately 135 thousand years ago. Given that this and all subsequent divisions led to populations with full linguistic capacity, it is reasonable to assume that the potential for language must have been present at the latest by around 135 thousand years ago, before the first division occurred. Had linguistic capacity developed later, we would expect to find some modern human populations without language, or with some fundamentally different mode of communication. Neither is the case. While current evidence does not tell us exactly when language itself appeared, the genomic studies do allow a fairly accurate estimate of the time by which linguistic capacity must have been present in the modern human lineage. Based on the lower boundary of 135 thousand years ago for language, we propose that language may have triggered the widespread appearance of modern human behavior approximately 100 thousand years ago.

Frontal sinus size in South African Later Stone Age Holocene Khoe-San

Frontal size variation is comparatively poorly sampled among sub-Saharan African populations. This study assessed frontal sinus size in a sample of Khoe-San skeletal remains from South African Later Stone Age contexts. Volumes were determined from CT scans of 102 adult crania; individual sex could be estimated in 82 cases. Sinus volume is not sexually dimorphic in this sample. The lack of frontal sinus aplasia is concordant with the low incidences recorded for other sub-Saharan African and most other global populations save those that inhabit high latitudes. There is considerable variation in frontal sinus size among global populations, and the Khoe-San possess among the smallest. The Khoe-San have rather diminutive sinuses compared to sub-Saharan Bantu-speaking populations but resemble a northern African (Sudanese) population. Genetic studies indicate the earliest population divergence within Homo sapiens to have been between the Khoe-San and all other living groups, and that this likely occurred in Africa during the span of Marine Isotope Stages 8-6. There is scant information on frontal sinus development among Late Quaternary African fossils that are likely either closely related or attributable to Homo sapiens. Among these, the MIS 3 cranium from Hofmeyr, South Africa, exhibits distinct Khoe-San cranial affinities and despite its large size has a very small frontal sinus. This raises the possibility that the small frontal sinuses of the Holocene South African Khoe-San might be a feature retained from an earlier MIS 3 population.

Deep history of cultural and linguistic evolution among Central African hunter-gatherers

Human evolutionary history in Central Africa reflects a deep history of population connectivity. However, Central African hunter-gatherers (CAHGs) currently speak languages acquired from their neighbouring farmers. Hence it remains unclear which aspects of CAHG cultural diversity results from long-term evolution preceding agriculture and which reflect borrowing from farmers. On the basis of musical instruments, foraging tools, specialized vocabulary and genome-wide data from ten CAHG populations, we reveal evidence of large-scale cultural interconnectivity among CAHGs before and after the Bantu expansion. We also show that the distribution of hunter-gatherer musical instruments correlates with the oldest genomic segments in our sample predating farming. Music-related words are widely shared between western and eastern groups and likely precede the borrowing of Bantu languages. In contrast, subsistence tools are less frequently exchanged and may result from adaptation to local ecologies. We conclude that CAHG material culture and specialized lexicon reflect a long evolutionary history in Central Africa.

Changes in weight status among "Ju/'hoansi" San hunter-gatherers between 1968/69 and 1987-The effects of nutritional transition, sex and age

BACKGROUND

Changes in lifestyle and dietary habits that hunter-gatherer populations have undergone in recent decades have often led to rising obesity rates with disastrous consequences for their health.

OBJECTIVES

The associations between dietary habits and weight status were studied in 238 "Ju/'hoansi" San (93 women and 145 men) aged between 18 and 65 years in northern Namibia in 1987. Weight status was estimated based on the World Health Organization body mass index (BMI) categories, and dietary habits were recorded using food recall methods. Anthropometrics and weight status were compared with those of a sex- and age-matched sample of "Ju/'hoansi" San people collected by Nancy Howell in 1968/69.

RESULTS

Body weight had increased significantly among "Ju/'hoansi" San people from 1968/69 to 1987. The number of underweight people decreased from 1968/69 to 1987. In 1987, most participants (60.9%) were of normal weight. Overweight was found in 1.3% of the women, but not among men. No participants were obese. Less than 4% of the women and less than 2% of the men consumed exclusively traditional hunter-gatherer food. Westernized food products were significantly (p < 0.001) more common among men and younger people. Dietary patterns were significantly associated with weight status. The less traditional the diet, the higher the BMI (p < 0.001).

CONCLUSION

The transition to domestic agricultural and westernized foods was positively associated with increasing BMI. Overweight, however, was still an extremely rare condition in this population in 1987.

Quantifying hominin morphological diversity at the end of the middle Pleistocene: Implications for the origin of Homo sapiens

OBJECTIVES

The Middle Pleistocene (MP) saw the emergence of new species of hominins: Homo sapiens in Africa, H. neanderthalensis, and possibly Denisovans in Eurasia, whose most recent common ancestor is thought to have lived in Africa around 600 ka ago. However, hominin remains from this period present a wide range of morphological variation making it difficult to securely determine their taxonomic attribution and their phylogenetic position within the Homo genus. This study proposes to reconsider the phenetic relationships between MP hominin fossils in order to clarify evolutionary trends and contacts between the populations they represent.

MATERIALS AND METHODS

We used a Geometric Morphometrics approach to quantify the morphological variation of the calvarium of controversial MP specimens from Africa and Eurasia by using a comparative sample that can be divided into 5 groups: H. ergaster, H. erectus, H. neanderthalensis, and H. sapiens, as well as individuals from current modern human populations. We performed a Generalized Procrustes Analysis, a Principal Component Analysis, and Multinomial Principal Component Logistic Regressions to determine the phenetic affinities of the controversial Middle Pleistocene specimens with the other groups.

RESULTS

MP African and Eurasian specimens represent several populations, some of which show strong affinities with H. neanderthalensis in Europe or H. sapiens in Africa, others presenting multiple affinities.

DISCUSSION

These MP populations might have contributed to the emergence of these two species in different proportions. This study proposes a new framework for the human evolutionary history during the MP.

The Precision and Power of Population Branch Statistics in Identifying the Genomic Signatures of Local Adaptation

Population branch statistics, which estimate the branch lengths of focal populations with respect to two outgroups, have been used as an alternative to FST-based genome-wide scans for identifying loci associated with local selective sweeps. In addition to the original population branch statistic (PBS), there are subsequently proposed branch rescalings: normalized population branch statistic (PBSn1), which adjusts focal branch length with respect to outgroup branch lengths at the same locus, and population branch excess (PBE), which also incorporates median branch lengths at other loci. PBSn1 and PBE have been proposed to be less sensitive to allele frequency divergence generated by background selection or geographically ubiquitous positive selection rather than local selective sweeps. However, the accuracy and statistical power of branch statistics have not been systematically assessed. To do so, we simulate genomes in representative large and small populations with varying proportions of sites evolving under genetic drift or background selection (approximated using variable N e ), local selective sweeps, and geographically parallel selective sweeps. We then assess the probability that local selective sweep loci are correctly identified as outliers by FST and by each of the branch statistics. We find that branch statistics consistently outperform FST at identifying local sweeps. When background selection and/or parallel sweeps are introduced, PBSn1 and especially PBE correctly identify local sweeps among their top outliers at a higher frequency than PBS. These results validate the greater specificity of rescaled branch statistics such as PBE to detect population-specific positive selection, supporting their use in genomic studies focused on local adaptation.

Ghosts of extinct apes: genomic insights into African hominid evolution

We are accustomed to regular announcements of new hominin fossils. There are now some 6000 hominin fossils, and up to 31 species. However, where are the announcements of African ape fossils? The answer is that there are almost none. Our knowledge of African ape evolution is based entirely on genomic analyses, which show that extant diversity is very young. This contrasts with the extensive and deep diversity of hominins known from fossils. Does this difference point to low and late diversification of ape lineages, or high rates of extinction? The comparative evolutionary dynamics of African hominids are central to interpreting living ape adaptations, as well as understanding the patterns of hominin evolution and the nature of the last common ancestor.

The genetic changes that shaped Neandertals, Denisovans, and modern humans

Modern human ancestors diverged from the ancestors of Neandertals and Denisovans about 600,000 years ago. Until about 40,000 years ago, these three groups existed in parallel, occasionally met, and exchanged genes. A critical question is why modern humans, and not the other two groups, survived, became numerous, and developed complex cultures. Here, we discuss genetic differences among the groups and some of their functional consequences. As more present-day genome sequences become available from diverse groups, we predict that very few, if any, differences will distinguish all modern humans from all Neandertals and Denisovans. We propose that the genetic basis of what constitutes a modern human is best thought of as a combination of genetic features, where perhaps none of them is present in each and every present-day individual.

Hunter-gatherer genetics research: Importance and avenues

Major developments in the field of genetics in the past few decades have revolutionised notions of what it means to be human. Although currently only a few populations around the world practise a hunting and gathering lifestyle, this mode of subsistence has characterised members of our species since its very origins and allowed us to migrate across the planet. Therefore, the geographical distribution of hunter-gatherer populations, dependence on local ecosystems and connections to past populations and neighbouring groups have provided unique insights into our evolutionary origins. However, given the vulnerable status of hunter-gatherers worldwide, the development of the field of anthropological genetics requires that we reevaluate how we conduct research with these communities. Here, we review how the inclusion of hunter-gatherer populations in genetics studies has advanced our understanding of human origins, ancient population migrations and interactions as well as phenotypic adaptations and adaptability to different environments, and the important scientific and medical applications of these advancements. At the same time, we highlight the necessity to address yet unresolved questions and identify areas in which the field may benefit from improvements.

Entheseal changes and activity patterns in southern African hunter-gatherer/herders from the Holocene

OBJECTIVES

Activity patterns and lifeways in southern African hunter-gatherer/herders (sAHGH) during the Holocene were dynamic, with subsistence activities and mobility varying through space and time. In this study, spatial and temporal variations in entheseal changes (ECs) are assessed as physical activity markers in sAHGH from the Holocene.

METHODS

The Coimbra method was used to assess fibrocartilaginous ECs in the upper and lower limbs of 118 sAHGH from the Holocene. Descriptive statistics and generalized estimating equations were used to explore the association between ECs, sex, age, ecological biomes, and temporality.

RESULTS

A total of 118 individuals were sampled, comprising 67 males, and 42 females, mostly from the fynbos (59/118), forest (30/118) and succulent karoo biomes (16/118). ECs were identified in 94% of the sample. Interobserver scoring suggests our findings are likely to underrepresent the extent of EC score differences in sAHGH. Findings indicate a complex pattern of physical activity in sAHGH with differences attributable to regional ecology rather than age or sex. More prominent ECs were identified in individuals from the forest biome, relative to the fynbos and succulent karoo biomes. These were consistent with resource search and processing costs, and terrain differences in these ecozones. ECs were only detected temporally relative to the infiltration of pastoralism (at 2000 BP) when the data were stratified by ecological biome.

DISCUSSION

This study provides evidence that regional ecology plays a central role in driving activity patterns regardless of social or cultural organization. Ecological biomes provide a gradient along which the temporal impact of resource limitations on human biology, activity patterns and sociocultural behaviors can be studied.

The genetic legacy of the expansion of Bantu-speaking peoples in Africa

The expansion of people speaking Bantu languages is the most dramatic demographic event in Late Holocene Africa and fundamentally reshaped the linguistic, cultural and biological landscape of the continent1-7. With a comprehensive genomic dataset, including newly generated data of modern-day and ancient DNA from previously unsampled regions in Africa, we contribute insights into this expansion that started 6,000-4,000 years ago in western Africa. We genotyped 1,763 participants, including 1,526 Bantu speakers from 147 populations across 14 African countries, and generated whole-genome sequences from 12 Late Iron Age individuals8. We show that genetic diversity amongst Bantu-speaking populations declines with distance from western Africa, with current-day Zambia and the Democratic Republic of Congo as possible crossroads of interaction. Using spatially explicit methods9 and correlating genetic, linguistic and geographical data, we provide cross-disciplinary support for a serial-founder migration model. We further show that Bantu speakers received significant gene flow from local groups in regions they expanded into. Our genetic dataset provides an exhaustive modern-day African comparative dataset for ancient DNA studies10 and will be important to a wide range of disciplines from science and humanities, as well as to the medical sector studying human genetic variation and health in African and African-descendant populations.

Whole genome sequencing reveals population diversity and variation in HIV-1 specific host genes

HIV infection continues to be a major global public health issue. The population heterogeneity in susceptibility or resistance to HIV-1 and progression upon infection is attributable to, among other factors, host genetic variation. Therefore, identifying population-specific variation and genetic modifiers of HIV infectivity can catapult the invention of effective strategies against HIV-1 in African populations. Here, we investigated whole genome sequences of 390 unrelated HIV-positive and -negative individuals from Botswana. We report 27.7 million single nucleotide variations (SNVs) in the complete genomes of Botswana nationals, of which 2.8 million were missing in public databases. Our population structure analysis revealed a largely homogenous structure in the Botswana population. Admixture analysis showed elevated components shared between the Botswana population and the Niger-Congo (65.9%), Khoe-San (32.9%), and Europeans (1.1%) ancestries in the population of Botswana. Statistical significance of the mutational burden of deleterious and loss-of-function variants per gene against a null model was estimated. The most deleterious variants were enriched in five genes: ACTRT2 (the Actin Related Protein T2), HOXD12 (homeobox D12), ABCB5 (ATP binding cassette subfamily B member 5), ATP8B4 (ATPase phospholipid transporting 8B4) and ABCC12 (ATP Binding Cassette Subfamily C Member 12). These genes are enriched in the glycolysis and gluconeogenesis (p < 2.84e-6) pathways and therefore, may contribute to the emerging field of immunometabolism in which therapy against HIV-1 infection is being evaluated. Published transcriptomic evidence supports the role of the glycolysis/gluconeogenesis pathways in the regulation of susceptibility to HIV, and that cumulative effects of genetic modifiers in glycolysis/gluconeogenesis pathways may potentially have effects on the expression and clinical variability of HIV-1. Identified genes and pathways provide novel avenues for other interventions, with the potential for informing the design of new therapeutics.

Genome-wide variation in the Angolan Namib Desert reveals unique pre-Bantu ancestry

Ancient DNA studies reveal the genetic structure of Africa before the expansion of Bantu-speaking agriculturalists; however, the impact of now extinct hunter-gatherer and herder societies on the genetic makeup of present-day African groups remains elusive. Here, we uncover the genetic legacy of pre-Bantu populations from the Angolan Namib Desert, where we located small-scale groups associated with enigmatic forager traditions, as well as the last speakers of the Khoe-Kwadi family's Kwadi branch. By applying an ancestry decomposition approach to genome-wide data from these and other African populations, we reconstructed the fine-scale histories of contact emerging from the migration of Khoe-Kwadi-speaking pastoralists and identified a deeply divergent ancestry, which is exclusively shared between groups from the Angolan Namib and adjacent areas of Namibia. The unique genetic heritage of the Namib peoples shows how modern DNA research targeting understudied regions of high ethnolinguistic diversity can complement ancient DNA studies in probing the deep genetic structure of the African continent.

Genetic continuity, isolation, and gene flow in Stone Age Central and Eastern Europe

The genomic landscape of Stone Age Europe was shaped by multiple migratory waves and population replacements, but different regions do not all show similar patterns. To refine our understanding of the population dynamics before and after the dawn of the Neolithic, we generated and analyzed genomic sequence data from human remains of 56 individuals from the Mesolithic, Neolithic, and Eneolithic across Central and Eastern Europe. We found that Mesolithic European populations formed a geographically widespread isolation-by-distance zone ranging from Central Europe to Siberia, which was already established 10,000 years ago. We found contrasting patterns of population continuity during the Neolithic transition: people around the lower Dnipro Valley region, Ukraine, showed continuity over 4000 years, from the Mesolithic to the end of the Neolithic, in contrast to almost all other parts of Europe where population turnover drove this cultural change, including vast areas of Central Europe and around the Danube River.

Performance and accuracy evaluation of reference panels for genotype imputation in sub-Saharan African populations

Based on evaluations of imputation performed on a genotype dataset consisting of about 11,000 sub-Saharan African (SSA) participants, we show Trans-Omics for Precision Medicine (TOPMed) and the African Genome Resource (AGR) to be currently the best panels for imputing SSA datasets. We report notable differences in the number of single-nucleotide polymorphisms (SNPs) that are imputed by different panels in datasets from East, West, and South Africa. Comparisons with a subset of 95 SSA high-coverage whole-genome sequences (WGSs) show that despite being about 20-fold smaller, the AGR imputed dataset has higher concordance with the WGSs. Moreover, the level of concordance between imputed and WGS datasets was strongly influenced by the extent of Khoe-San ancestry in a genome, highlighting the need for integration of not only geographically but also ancestrally diverse WGS data in reference panels for further improvement in imputation of SSA datasets. Approaches that integrate imputed data from different panels could also lead to better imputation.

Hypercementosis in Late Pleistocene Homo sapiens fossils from Klasies River Main Site, South Africa

OBJECTIVE

To examine early Homo sapiens fossils from the Late Pleistocene site of Klasies River Main Site, South Africa for evidence of hypercementosis. The specimens represent seven adult individuals dated to between 119,000 and 58,000 years ago. These observations are contextualized in relation to the incidences of hypercementosis among recent human populations and fossil human samples and the potential etiologies of hypercementosis.

DESIGN

The fossils were investigated utilizing micro-CT and nano-CT scanning to visualize and measure cementum apposition on permanent incisor, premolar and molar roots. Cementum thickness was measured at mid-root level, and the volume of the cementum sleeve was calculated for the two fossil specimens that display marked hypercementosis.

RESULTS

Two of the fossils display no evidence of cementum hypertrophy. Three exhibit moderate cementum thickening, barely attaining the quantitative threshold for hypercementosis. Two evince marked hypercementosis. One of the Klasies specimens with marked hypercementosis is judged to be an older individual with periapical abscessing. The second specimen is a younger adult, and seemingly similar in age to other Klasies fossils that exhibit only minimal cementum apposition. However, this second specimen exhibits dento-alveolar ankylosis of the premolar and molars.

CONCLUSIONS

These two fossils from Klasies River Main Site provide the earliest manifestation of hypercementosis in Homo sapiens.

Evolutionary Genetics and Admixture in African Populations

As the ancestral homeland of our species, Africa contains elevated levels of genetic diversity and substantial population structure. Importantly, African genomes are heterogeneous: They contain mixtures of multiple ancestries, each of which have experienced different evolutionary histories. In this review, we view population genetics through the lens of admixture, highlighting how multiple demographic events have shaped African genomes. Each of these historical vignettes paints a recurring picture of population divergence followed by secondary contact. First, we give a brief overview of genetic variation in Africa and examine deep population structure within Africa, including the evidence of ancient introgression from archaic "ghost" populations. Second, we describe the genetic legacies of admixture events that have occurred during the past 10,000 years. This includes gene flow between different click-speaking Khoe-San populations, the stepwise spread of pastoralism from eastern to southern Africa, multiple migrations of Bantu speakers across the continent, as well as admixture from the Middle East and Europe into the Sahel region and North Africa. Furthermore, the genomic signatures of more recent admixture can be found in the Cape Peninsula and throughout the African diaspora. Third, we highlight how natural selection has shaped patterns of genetic variation across the continent, noting that gene flow provides a potent source of adaptive variation and that selective pressures vary across Africa. Finally, we explore the biomedical implications of population structure in Africa on health and disease and call for more ethically conducted studies of genetic variation in Africa.

Rickettsia felis DNA recovered from a child who lived in southern Africa 2000 years ago

The Stone Age record of South Africa provides some of the earliest evidence for the biological and cultural origins of Homo sapiens. While there is extensive genomic evidence for the selection of polymorphisms in response to pathogen-pressure in sub-Saharan Africa, e.g., the sickle cell trait which provides protection against malaria, there is inadequate direct human genomic evidence for ancient human-pathogen infection in the region. Here, we analysed shotgun metagenome libraries derived from the sequencing of a Later Stone Age hunter-gatherer child who lived near Ballito Bay, South Africa, c. 2000 years ago. This resulted in the identification of ancient DNA sequence reads homologous to Rickettsia felis, the causative agent of typhus-like flea-borne rickettsioses, and the reconstruction of an ancient R. felis genome.

The recombination landscape of the Khoe-San likely represents the upper limits of recombination divergence in humans

BACKGROUND

Recombination maps are  important resources for epidemiological and evolutionary analyses; however, there are currently no recombination maps representing any African population outside of those with West African ancestry. We infer the demographic history for the Nama, an indigenous Khoe-San population of southern Africa, and derive a novel, population-specific recombination map from the whole genome sequencing of 54 Nama individuals. We hypothesise that there are no publicly available recombination maps representative of the Nama, considering the deep population divergence and subsequent isolation of the Khoe-San from other African groups.

RESULTS

We show that the recombination landscape of the Nama does not cluster with any continental groups with publicly available representative recombination maps. Finally, we use selection scans as an example of how fine-scale differences between the Nama recombination map and the combined Phase II HapMap recombination map can impact the outcome of selection scans.

CONCLUSIONS

Fine-scale differences in recombination can meaningfully alter the results of a selection scan. The recombination map we infer likely represents an upper bound on the extent of divergence we expect to see for a recombination map in humans and would be of interest to any researcher that wants to test the sensitivity of population genetic or GWAS analysis to recombination map input.

Temporal mapping of derived high-frequency gene variants supports the mosaic nature of the evolution of Homo sapiens

Large-scale estimations of the time of emergence of variants are essential to examine hypotheses concerning human evolution with precision. Using an open repository of genetic variant age estimations, we offer here a temporal evaluation of various evolutionarily relevant datasets, such as Homo sapiens-specific variants, high-frequency variants found in genetic windows under positive selection, introgressed variants from extinct human species, as well as putative regulatory variants specific to various brain regions. We find a recurrent bimodal distribution of high-frequency variants, but also evidence for specific enrichments of gene categories in distinct time windows, pointing to different periods of phenotypic changes, resulting in a mosaic. With a temporal classification of genetic mutations in hand, we then applied a machine learning tool to predict what genes have changed more in certain time windows, and which tissues these genes may have impacted more. Overall, we provide a fine-grained temporal mapping of derived variants in Homo sapiens that helps to illuminate the intricate evolutionary history of our species.

Population interconnectivity over the past 120,000 years explains distribution and diversity of Central African hunter-gatherers

We combined ethnographic, archaeological, genetic, and paleoclimatic data to model the dynamics of Central African hunter-gatherer populations over the past 120,000 years. We show, against common assumptions, that their distribution and density are explained by changing environments rather than by a displacement following recent farming expansions, and that they have maintained large population sizes and genetic diversity, despite fluctuations in niche availability. Our results provide insights into the evolution of genetic and cultural diversity in Homo sapiens.

The evolutionary history of African hunter-gatherers holds key insights into modern human diversity. Here, we combine ethnographic and genetic data on Central African hunter-gatherers (CAHG) to show that their current distribution and density are explained by ecology rather than by a displacement to marginal habitats due to recent farming expansions, as commonly assumed. We also estimate the range of hunter-gatherer presence across Central Africa over the past 120,000 years using paleoclimatic reconstructions, which were statistically validated by our newly compiled dataset of dated archaeological sites. Finally, we show that genomic estimates of divergence times between CAHG groups match our ecological estimates of periods favoring population splits, and that recoveries of connectivity would have facilitated subsequent gene flow. Our results reveal that CAHG stem from a deep history of partially connected populations. This form of sociality allowed the coexistence of relatively large effective population sizes and local differentiation, with important implications for the evolution of genetic and cultural diversity in Homo sapiens.

Climate effects on archaic human habitats and species successions

It has long been believed that climate shifts during the last 2 million years had a pivotal role in the evolution of our genus Homo^1–3. However, given the limited number of representative palaeo-climate datasets from regions of anthropological interest, it has remained challenging to quantify this linkage. Here, we use an unprecedented transient Pleistocene coupled general circulation model simulation in combination with an extensive compilation of fossil and archaeological records to study the spatiotemporal habitat suitability for five hominin species over the past 2 million years. We show that astronomically forced changes in temperature, rainfall and terrestrial net primary production had a major impact on the observed distributions of these species. During the Early Pleistocene, hominins settled primarily in environments with weak orbital-scale climate variability. This behaviour changed substantially after the mid-Pleistocene transition, when archaic humans became global wanderers who adapted to a wide range of spatial climatic gradients. Analysis of the simulated hominin habitat overlap from approximately 300–400 thousand years ago further suggests that antiphased climate disruptions in southern Africa and Eurasia contributed to the evolutionary transformation of Homo heidelbergensis populations into Homo sapiens and Neanderthals, respectively. Our robust numerical simulations of climate-induced habitat changes provide a framework to test hypotheses on our human origin.

A new model simulation of climate change during the past 2 million years indicates that the appearances and disappearances of hominin species correlate with long-term climatic anomalies.

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.

Ostrich eggshell beads reveal 50,000-year-old social network in Africa

Humans evolved in a patchwork of semi-connected populations across Africa1,2; understanding when and how these groups connected is critical to interpreting our present-day biological and cultural diversity. Genetic analyses reveal that eastern and southern African lineages diverged sometime in the Pleistocene epoch, approximately 350-70 thousand years ago (ka)3,4; however, little is known about the exact timing of these interactions, the cultural context of these exchanges or the mechanisms that drove their separation. Here we compare ostrich eggshell bead variations between eastern and southern Africa to explore population dynamics over the past 50,000 years. We found that ostrich eggshell bead technology probably originated in eastern Africa and spread southward approximately 50-33 ka via a regional network. This connection breaks down approximately 33 ka, with populations remaining isolated until herders entered southern Africa after 2 ka. The timing of this disconnection broadly corresponds with the southward shift of the Intertropical Convergence Zone, which caused periodic flooding of the Zambezi River catchment (an area that connects eastern and southern Africa). This suggests that climate exerted some influence in shaping human social contact. Our study implies a later regional divergence than predicted by genetic analyses, identifies an approximately 3,000-kilometre stylistic connection and offers important new insights into the social dimension of ancient interactions.

Male-biased migration from East Africa introduced pastoralism into southern Africa

BACKGROUND

Hunter-gatherer lifestyles dominated the southern African landscape up to ~ 2000 years ago, when herding and farming groups started to arrive in the area. First, herding and livestock, likely of East African origin, appeared in southern Africa, preceding the arrival of the large-scale Bantu-speaking agro-pastoralist expansion that introduced West African-related genetic ancestry into the area. Present-day Khoekhoe-speaking Namaqua (or Nama in short) pastoralists show high proportions of East African admixture, linking the East African ancestry with Khoekhoe herders. Most other historical Khoekhoe populations have, however, disappeared over the last few centuries and their contribution to the genetic structure of present-day populations is not well understood. In our study, we analyzed genome-wide autosomal and full mitochondrial data from a population who trace their ancestry to the Khoekhoe-speaking Hessequa herders from the southern Cape region of what is now South Africa.

RESULTS

We generated genome-wide data from 162 individuals and mitochondrial DNA data of a subset of 87 individuals, sampled in the Western Cape Province, South Africa, where the Hessequa population once lived. Using available comparative data from Khoe-speaking and related groups, we aligned genetic date estimates and admixture proportions to the archaeological proposed dates and routes for the arrival of the East African pastoralists in southern Africa. We identified several Afro-Asiatic-speaking pastoralist groups from Ethiopia and Tanzania who share high affinities with the East African ancestry present in southern Africa. We also found that the East African pastoralist expansion was heavily male-biased, akin to a pastoralist migration previously observed on the genetic level in ancient Europe, by which Pontic-Caspian Steppe pastoralist groups represented by the Yamnaya culture spread across the Eurasian continent during the late Neolithic/Bronze Age.

CONCLUSION

We propose that pastoralism in southern Africa arrived through male-biased migration of an East African Afro-Asiatic-related group(s) who introduced new subsistence and livestock practices to local southern African hunter-gatherers. Our results add to the understanding of historical human migration and mobility in Africa, connected to the spread of food-producing and livestock practices.

Revisiting the out of Africa event with a deep-learning approach

Anatomically modern humans evolved around 300 thousand years ago in Africa. They started to appear in the fossil record outside of Africa as early as 100 thousand years ago, although other hominins existed throughout Eurasia much earlier. Recently, several studies argued in favor of a single out of Africa event for modern humans on the basis of whole-genome sequence analyses. However, the single out of Africa model is in contrast with some of the findings from fossil records, which support two out of Africa events, and uniparental data, which propose a back to Africa movement. Here, we used a deep-learning approach coupled with approximate Bayesian computation and sequential Monte Carlo to revisit these hypotheses from the whole-genome sequence perspective. Our results support the back to Africa model over other alternatives. We estimated that there are two sequential separations between Africa and out of African populations happening around 60-90 thousand years ago and separated by 13-15 thousand years. One of the populations resulting from the more recent split has replaced the older West African population to a large extent, while the other one has founded the out of Africa populations.

Variation in Middle Stone Age mandibular molar enamel-dentine junction topography at Klasies River Main Site assessed by diffeomorphic surface matching

The morphology and variability of the Middle Stone Age (MSA) hominin fossils from Klasies River Main Site have been the focus of investigation for more than four decades. The mandibular remains have figured prominently in discussions relating to robusticity, size dimorphism, and symphyseal morphology. Variation in corpus size between the robust SAM-AP 6223 and the diminutive SAM-AP 6225 mandibles is particularly impressive, and the difference between the buccolingual diameters of their M₂s significantly exceeds recent human sample variation. SAM-AP 6223 and SAM-AP 6225 are the only Klasies specimens with homologous teeth (M₂ and M₃) that permit comparisons of crown morphology. While the differences in dental trait expression at the outer enamel surfaces of these molars are slight, diffeomorphic surface analyses of their underlying enamel-dentine junction (EDJ) topographies reveal differences that are well beyond the means of pairwise differences among comparative samples of Later Stone Age (LSA) Khoesan and recent African homologues. The EDJs of both SAM-AP 6225 molars and the SAM-AP 6223 M₃ fall outside the envelopes that define the morphospace of these two samples. Although the radiocarbon dated LSA individuals examined here differ by a maximum of some 7000 years, and the two Klasies jaws may differ by perhaps as much as 18,000 years, it is difficult to ascribe their differences to time alone. With reference to the morphoscopic traits by which the SAM-AP 6223 and SAM-AP 6225 EDJs differ, the most striking is the expression of the protoconid cingulum. This is very weakly developed on the SAM-AP 6223 molars and distinct in SAM-AP 6225. As such, this diminutive fossil exhibits a more pronounced manifestation of what is likely a plesiomorphic feature, thus adding to the morphological mosaicism that is evident in the Klasies hominin assemblage. Several possible explanations for the variation and mosaicism in this MSA sample are discussed.

Comparison of sequencing data processing pipelines and application to underrepresented African human populations

Background

Population genetic studies of humans make increasing use of high-throughput sequencing in order to capture diversity in an unbiased way. There is an abundance of sequencing technologies, bioinformatic tools and the available genomes are increasing in number. Studies have evaluated and compared some of these technologies and tools, such as the Genome Analysis Toolkit (GATK) and its “Best Practices” bioinformatic pipelines. However, studies often focus on a few genomes of Eurasian origin in order to detect technical issues. We instead surveyed the use of the GATK tools and established a pipeline for processing high coverage full genomes from a diverse set of populations, including Sub-Saharan African groups, in order to reveal challenges from human diversity and stratification.

Results

We surveyed 29 studies using high-throughput sequencing data, and compared their strategies for data pre-processing and variant calling. We found that processing of data is very variable across studies and that the GATK “Best Practices” are seldom followed strictly. We then compared three versions of a GATK pipeline, differing in the inclusion of an indel realignment step and with a modification of the base quality score recalibration step. We applied the pipelines on a diverse set of 28 individuals. We compared the pipelines in terms of count of called variants and overlap of the callsets. We found that the pipelines resulted in similar callsets, in particular after callset filtering. We also ran one of the pipelines on a larger dataset of 179 individuals. We noted that including more individuals at the joint genotyping step resulted in different counts of variants. At the individual level, we observed that the average genome coverage was correlated to the number of variants called.

Conclusions

We conclude that applying the GATK “Best Practices” pipeline, including their recommended reference datasets, to underrepresented populations does not lead to a decrease in the number of called variants compared to alternative pipelines. We recommend to aim for coverage of > 30X if identifying most variants is important, and to work with large sample sizes at the variant calling stage, also for underrepresented individuals and populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04407-x.

Homo sapiens origins and evolution in the Kalahari Basin, southern Africa

The Kalahari Basin, southern Africa preserves a rich archeological record of human origins and evolution spanning the Early, Middle and Late Pleistocene. Since the 1930s, several stratified and dated archeological sites have been identified and investigated, together with numerous open-air localities that provide landscape-scale perspectives. However, next to recent discoveries from nearby coastal regions, the Kalahari Basin has remained peripheral to debates about the origins of Homo sapiens. Though the interior region of southern Africa is generally considered to be less suitable for hunter-gatherer occupation than coastal and near-coastal regions, especially during glacial periods, the archeological record documents human presence in the Kalahari Basin from the Early Pleistocene onwards, and the region is not abandoned during glacial phases. Furthermore, many significant behavioral innovations have an early origin in the Kalahari Basin, which adds support to poly-centric, pan-African models for the emergence of our species.

African Global Representation in Biomedical Sciences

African populations are diverse in their ethnicity, language, culture, and genetics. Although plagued by high disease burdens, until recently the continent has largely been excluded from biomedical studies. Along with limitations in research and clinical infrastructure, human capacity, and funding, this omission has resulted in an underrepresentation of African data and disadvantaged African scientists. This review interrogates the relative abundance of biomedical data from Africa, primarily in genomics and other omics. The visibility of African science through publications is also discussed. A challenge encountered in this review is the relative lack of annotation of data on their geographical or population origin, with African countries represented as a single group. In addition to the abovementioned limitations,the global representation of African data may also be attributed to the hesitation to deposit data in public repositories. Whatever the reason, the disparity should be addressed, as African data have enormous value for scientists in Africa and globally.

Life in Deserts: The Genetic Basis of Mammalian Desert Adaptation

Deserts are among the harshest environments on Earth. The multiple ages of different deserts and their global distribution provide a unique opportunity to study repeated adaptation at different timescales. Here, we summarize recent genomic research on the genetic mechanisms underlying desert adaptations in mammals. Several studies on different desert mammals show large overlap in functional classes of genes and pathways, consistent with the complexity and variety of phenotypes associated with desert adaptation to water and food scarcity and extreme temperatures. However, studies of desert adaptation are also challenged by a lack of accurate genotype-phenotype-environment maps. We encourage development of systems that facilitate functional analyses, but also acknowledge the need for more studies on a wider variety of desert mammals.

Bantu-speaker migration and admixture in southern Africa

The presence of Early and Middle Stone Age human remains and associated archeological artifacts from various sites scattered across southern Africa, suggests this geographic region to be one of the first abodes of anatomically modern humans. Although the presence of hunter-gatherer cultures in this region dates back to deep times, the peopling of southern Africa has largely been reshaped by three major sets of migrations over the last 2000 years. These migrations have led to a confluence of four distinct ancestries (San hunter-gatherer, East-African pastoralist, Bantu-speaker farmer and Eurasian) in populations from this region. In this review, we have summarized the recent insights into the refinement of timelines and routes of the migration of Bantu-speaking populations to southern Africa and their admixture with resident southern African Khoe-San populations. We highlight two recent studies providing evidence for the emergence of fine-scale population structure within some South-Eastern Bantu-speaker groups. We also accentuate whole genome sequencing studies (current and ancient) that have both enhanced our understanding of the peopling of southern Africa and demonstrated a huge potential for novel variant discovery in populations from this region. Finally, we identify some of the major gaps and inconsistencies in our understanding and emphasize the importance of more systematic studies of southern African populations from diverse ethnolinguistic groups and geographic locations.

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.

Later Stone Age human hair from Vaalkrans Shelter, Cape Floristic Region of South Africa, reveals genetic affinity to Khoe groups

Previous studies show that the indigenous people of the southern Cape of South Africa were dramatically impacted by the arrival of European colonists starting ~400 years ago and their descendants are today mixed with Europeans and Asians. To gain insight on the occupants of the Vaalkrans Shelter located at the southernmost tip of Africa, we investigated the genetic make-up of an individual who lived there about 200 years ago. We further contextualize the genetic ancestry of this individual among prehistoric and current groups. From a hair sample excavated at the shelter, which was indirectly dated to about 200 years old, we sequenced the genome (1.01 times coverage) of a Later Stone Age individual. We analyzed the Vaalkrans genome together with genetic data from 10 ancient (pre-colonial) individuals from southern Africa spanning the last 2000 years. We show that the individual from Vaalkrans was a man who traced ~80% of his ancestry to local southern San hunter-gatherers and ~20% to a mixed East African-Eurasian source. This genetic make-up is similar to modern-day Khoekhoe individuals from the Northern Cape Province (South Africa) and Namibia, but in the southern Cape, the Vaalkrans man's descendants have likely been assimilated into mixed-ancestry "Coloured" groups. The Vaalkrans man's genome reveals that Khoekhoe pastoralist groups/individuals lived in the southern Cape as late as 200 years ago, without mixing with non-African colonists or Bantu-speaking farmers. Our findings are also consistent with the model of a Holocene pastoralist migration, originating in Eastern Africa, shaping the genomic landscape of historic and current southern African populations.

Sub-Saharan African information potential to unveil adaptations to infectious disease

Sub-Saharan Africa is the most promising region of the world to conduct high-throughput studies to unveil adaptations to infectious diseases due to several reasons, namely, the longest evolving time-depth in the Homo sapiens phylogenetic tree (at least two-third older than any other worldwide region); the continuous burden of infectious diseases (still number one in health/life threat); and the coexistence of populations practising diverse subsistence modes (nomadic or seminomadic hunter-gatherers and agropastoralists, and sedentary agriculturalists, small urban and megacity groups). In this review, we will present the most up-to-date results that shed light on three main hypotheses related with this adaptation. One is the hypothesis of coevolution between host and pathogen, given enough time for the establishment of this highly dynamic relationship. The second hypothesis enunciates that the agricultural transition was responsible for the increase of the infectious disease burden, due to the huge expansion of the sedentary human population and the cohabitation with domesticates as main reservoirs of pathogens. The third hypothesis states that the boosting of our immune system against pathogens by past selection may have resulted in maladaptation of the developed hygienic societies, leading to an increase of allergic, inflammatory and autoimmune disorders. Further work will enlighten the biological mechanisms behind these main adaptations, which can be insightful for translation into diagnosis, prognosis and treatment interventions.

Closing the Gaps in Genomic Research

Despite Africa's central role in the origin of our species, our knowledge of the genomic diversity in Africa is remarkably sparse. A recent publication by Choudhury et al. underscores the scientific imperative for a broader characterisation of African genomic diversity to better understand demographic history and improve global human health.

The genomic prehistory of peoples speaking Khoisan languages

Peoples speaking so-called Khoisan languages-that is, indigenous languages of southern Africa that do not belong to the Bantu family-are culturally and linguistically diverse. They comprise herders, hunter-gatherers as well as groups of mixed modes of subsistence, and their languages are classified into three distinct language families. This cultural and linguistic variation is mirrored by extensive genetic diversity. We here review the recent genomics literature and discuss the genetic evidence for a formerly wider geographic spread of peoples with Khoisan-related ancestry, for the deep divergence among populations speaking Khoisan languages overlaid by more recent gene flow among these groups and for the impact of admixture with immigrant food-producers in their prehistory.

Regional patterns of diachronic technological change in the Howiesons Poort of southern Africa

The Howiesons Poort (HP) of southern Africa plays an important role in models on the early behavioral evolution of Homo sapiens. The HP is often portrayed as a coherent MSA industry characterized by early complex material culture. Recent work has emphasized parallel technological change through time across southern Africa potentially driven by ecological adaptations or demographic change. Here we examine patterns of diachronic variation within the HP and evaluate potential causal factors behind these changes. We test previous temporal assessments of the technocomplex at the local and regional level based on high-resolution quantitative data on HP lithic assemblages from Sibudu (KwaZulu-Natal) and comparisons with other southern African sites. At Sibudu, consistent unidirectional change in lithic technology characterizes the HP sequence. The results show a gradual reduction in typical HP markers such as the proportion of blades, backed pieces, and HP cores, as well as declining size of blades and backed artifacts. Quantitative comparisons with seven HP sites in South Africa suggest that lithic technology varies between regions over time instead of following similar changes. Concerning hypotheses of causal drivers, directional changes in lithic technology at Sibudu covary with shifting hunting patterns towards larger-sized bovids and a gradual opening of the vegetation. In contrast, variation in lithic technology shows little association with site use, mobility patterns or demographic expansions. Unlike at Sibudu, diachronic changes at other HP sites such as Diepkloof, Klasies River and Klipdrift appear to be associated with aspects of mobility, technological organization and site use. The regional diachronic patterns in the HP partly follow paleoclimatic zones, which could imply different ecological adaptations and distinct connection networks over time. Divergent and at times decoupled changes in lithic traits across sites precludes monocausal explanations for the entire HP, supporting more complex models for the observed technological trajectories.

Y-Chromosome Variation in Southern African Khoe-San Populations Based on Whole-Genome Sequences

Although the human Y chromosome has effectively shown utility in uncovering facets of human evolution and population histories, the ascertainment bias present in early Y-chromosome variant data sets limited the accuracy of diversity and TMRCA estimates obtained from them. The advent of next-generation sequencing, however, has removed this bias and allowed for the discovery of thousands of new variants for use in improving the Y-chromosome phylogeny and computing estimates that are more accurate. Here, we describe the high-coverage sequencing of the whole Y chromosome in a data set of 19 male Khoe-San individuals in comparison with existing whole Y-chromosome sequence data. Due to the increased resolution, we potentially resolve the source of haplogroup B-P70 in the Khoe-San, and reconcile recently published haplogroup A-M51 data with the most recent version of the ISOGG Y-chromosome phylogeny. Our results also improve the positioning of tentatively placed new branches of the ISOGG Y-chromosome phylogeny. The distribution of major Y-chromosome haplogroups in the Khoe-San and other African groups coincide with the emerging picture of African demographic history; with E-M2 linked to the agriculturalist Bantu expansion, E-M35 linked to pastoralist eastern African migrations, B-M112 linked to earlier east-south gene flow, A-M14 linked to shared ancestry with central African rainforest hunter-gatherers, and A-M51 potentially unique to the Khoe-San.