Reconstructing Prehistoric African Population Structure

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.

Along the Indian Ocean Coast: Genomic Variation in Mozambique Provides New Insights into the Bantu Expansion

The Bantu expansion, which started in West Central Africa around 5,000 BP, constitutes a major migratory movement involving the joint spread of peoples and languages across sub-Saharan Africa. Despite the rich linguistic and archaeological evidence available, the genetic relationships between different Bantu-speaking populations and the migratory routes they followed during various phases of the expansion remain poorly understood. Here, we analyze the genetic profiles of southwestern and southeastern Bantu-speaking peoples located at the edges of the Bantu expansion by generating genome-wide data for 200 individuals from 12 Mozambican and 3 Angolan populations using ∼1.9 million autosomal single nucleotide polymorphisms. Incorporating a wide range of available genetic data, our analyses confirm previous results favoring a “late split” between West and East Bantu speakers, following a joint passage through the rainforest. In addition, we find that Bantu speakers from eastern Africa display genetic substructure, with Mozambican populations forming a gradient of relatedness along a North–South cline stretching from the coastal border between Kenya and Tanzania to South Africa. This gradient is further associated with a southward increase in genetic homogeneity, and involved minimum admixture with resident populations. Together, our results provide the first genetic evidence in support of a rapid North–South dispersal of Bantu peoples along the Indian Ocean Coast, as inferred from the distribution and antiquity of Early Iron Age assemblages associated with the Kwale archaeological tradition.

Evolution of Human Brain Size-Associated NOTCH2NL Genes Proceeds toward Reduced Protein Levels

Ever since the availability of genomes from Neanderthals, Denisovans, and ancient humans, the field of evolutionary genomics has been searching for protein-coding variants that may hold clues to how our species evolved over the last ∼600,000 years. In this study, we identify such variants in the human-specific NOTCH2NL gene family, which were recently identified as possible contributors to the evolutionary expansion of the human brain. We find evidence for the existence of unique protein-coding NOTCH2NL variants in Neanderthals and Denisovans which could affect their ability to activate Notch signaling. Furthermore, in the Neanderthal and Denisovan genomes, we find unusual NOTCH2NL configurations, not found in any of the modern human genomes analyzed. Finally, genetic analysis of archaic and modern humans reveals ongoing adaptive evolution of modern human NOTCH2NL genes, identifying three structural variants acting complementary to drive our genome to produce a lower dosage of NOTCH2NL protein. Because copy-number variations of the 1q21.1 locus, encompassing NOTCH2NL genes, are associated with severe neurological disorders, this seemingly contradicting drive toward low levels of NOTCH2NL protein indicates that the optimal dosage of NOTCH2NL may have not yet been settled in the human population.

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

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

ContamLD: estimation of ancient nuclear DNA contamination using breakdown of linkage disequilibrium

We report a method called ContamLD for estimating autosomal ancient DNA (aDNA) contamination by measuring the breakdown of linkage disequilibrium in a sequenced individual due to the introduction of contaminant DNA. ContamLD leverages the idea that contaminants should have haplotypes uncorrelated to those of the studied individual. Using simulated data, we confirm that ContamLD accurately infers contamination rates with low standard errors: for example, less than 1.5% standard error in cases with less than 10% contamination and 500,000 sequences covering SNPs. This method is optimized for application to aDNA, taking advantage of characteristic aDNA damage patterns to provide calibrated contamination estimates, and is available at https://github.com/nathan-nakatsuka/ContamLD .

A southern African archaeological database of organic containers and materials, 800 cal BC to cal AD 1500: Possible implications for the transition from foraging to livestock-keeping

Although evidence of organic materials has consistently been reported in the archaeology of southern Africa little attention has been given to how this evidence, so slight in comparison to pottery and lithics, might be used to understand the transition from foraging to livestock-keeping in southern African Archaeology. We have compiled a geo-referenced, radiocarbon database of these organic, material culture remains, with particular reference to containers made of ostrich eggshell, wood, gourd, tortoise shell, twine, and leather over a 2300-year period to capture the periods before and after the appearance of livestock. We have mapped the organic materials for the period 800 cal BC to cal AD 1500 and explored the subsistence base of those who used them. This distribution is compared to that of pottery and livestock remains–conventionally the two archaeological markers of pastoralists. The paper interrogates what this might add to the vexed question of how the practice of livestock-keeping and pottery-making spread into and through the region (the hunter-herder debate). Our analysis suggests that ostrich eggshell containers can be used as a proxy for hunter-gatherers. By comparing areas of bead manufacture with those that have evidence only of bead use, we show the areas to which items may have travelled, along already established hunter-gatherer exchange networks. Our results suggest that hunter-gatherers widely and quickly adopted pottery across southern Africa in a process of cultural diffusion and local innovation, and that this was possibly the main mechanism for the dispersal of livestock at 2100 years ago.

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.

Methods for detecting introgressed archaic sequences

Analysis of genome sequences from archaic and modern humans have revealed multiple episodes of admixture between highly-diverged population groups. Statistical methods that attempt to localize DNA segments introduced by these events offer a powerful tool to investigate recent human evolution. We review recent advances in methods for detecting introgressed sequences.

Ancient genomes reveal complex patterns of population movement, interaction, and replacement in sub-Saharan Africa

Africa hosts the greatest human genetic diversity globally, but legacies of ancient population interactions and dispersals across the continent remain understudied. Here, we report genome-wide data from 20 ancient sub-Saharan African individuals, including the first reported ancient DNA from the DRC, Uganda, and Botswana. These data demonstrate the contraction of diverse, once contiguous hunter-gatherer populations, and suggest the resistance to interaction with incoming pastoralists of delayed-return foragers in aquatic environments. We refine models for the spread of food producers into eastern and southern Africa, demonstrating more complex trajectories of admixture than previously suggested. In Botswana, we show that Bantu ancestry post-dates admixture between pastoralists and foragers, suggesting an earlier spread of pastoralism than farming to southern Africa. Our findings demonstrate how processes of migration and admixture have markedly reshaped the genetic map of sub-Saharan Africa in the past few millennia and highlight the utility of combined archaeological and archaeogenetic approaches.

Human adaptation over the past 40,000 years

Over the past few years several methodological and data-driven advances have greatly improved our ability to robustly detect genomic signatures of selection in humans. New methods applied to large samples of present-day genomes provide increased power, while ancient DNA allows precise estimation of timing and tempo. However, despite these advances, we are still limited in our ability to translate these signatures into understanding about which traits were actually under selection, and why. Combining information from different populations and timescales may allow interpretation of selective sweeps. Other modes of selection have proved more difficult to detect. In particular, despite strong evidence of the polygenicity of most human traits, evidence for polygenic selection is weak, and its importance in recent human evolution remains unclear. Balancing selection and archaic introgression seem important for the maintenance of potentially adaptive immune diversity, but perhaps less so for other traits.

Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

The surgeon's role on chemical investigations of the composition of urinary stones

The chemical analysis of an urolith is often interpreted as "stone's composition". However, it must be taken into consideration, that in most cases, only a fragment of the stone has been sent to the laboratory. In some recurrent patients, stone compositions either vary considerably between episodes or the analytical result obtained from the stone fragment does not fit with the data of e.g. current 24 h-urinalysis or urinary pH-records. The question arises, whether this outcome may be the result of an improper stone sampling scheme. On a simple layered 2D-stone model composed of two mineral phases it is shown, how the choice of a stone fragment process may influence the result of "stone composition". Depending on the initial position of fragment within the whole stone, the respective calculated analyses can relevantly differ from the whole stone composition as well as strongly between two fragments. Even under the simplified conditions of a 2D-2-component-model "grown" under defined conditions, the differences between the analyses of the different specimens taken from a stone are in part remarkable. The more it can be argued that these differences increase if a real 3D-urolith is investigated. Further sampling biases may evolve and increase the problem of proper sampling:, e.g., if an urolith's more resistant parts remain intact while ESWL or laser-based stone fragmentation ("dusting"), the weak parts became fully disintegrated and removed from the body as fine-grained sludge-the stone's fine fraction is lost although its composition may carry important information on the stone's pathogenesis. Consequently, a "stone analysis" only obtained from the harder remains reveals an incomplete result, a fact that in principle limits its clinical interpretation. Choice of stone fragment is crucial. The extent of the uncertainty of an analysis resulting from potential selection biases should not be underestimated. Thus, sampling should be considered as an important part of the processes of quality assurance and management. Errors made at this early stage of diagnosis finding will affect the analytical result and thus influence the clarification of the underlying pathomechanism. This can lead to an improper metaphylactic strategy potentially causing recurrent stone formation which otherwise would have been prevented. A decision scheme for analysis of urinary stones removed using endoscopic methods is suggested.

Insights into human genetic variation and population history from 929 diverse genomes

Genome sequences from diverse human groups are needed to understand the structure of genetic variation in our species and the history of, and relationships between, different populations. We present 929 high-coverage genome sequences from 54 diverse human populations, 26 of which are physically phased using linked-read sequencing. Analyses of these genomes reveal an excess of previously undocumented common genetic variation private to southern Africa, central Africa, Oceania, and the Americas, but an absence of such variants fixed between major geographical regions. We also find deep and gradual population separations within Africa, contrasting population size histories between hunter-gatherer and agriculturalist groups in the past 10,000 years, and a contrast between single Neanderthal but multiple Denisovan source populations contributing to present-day human populations.

Trajectories of cultural innovation from the Middle to Later Stone Age in Eastern Africa: Personal ornaments, bone artifacts, and ocher from Panga ya Saidi, Kenya

African Middle Stone Age (MSA) populations used pigments, manufactured and wore personal ornaments, made abstract engravings, and produced fully shaped bone tools. However, ongoing research across Africa reveals variability in the emergence of cultural innovations in the MSA and their subsequent development through the Later Stone Age (LSA). When present, it appears that cultural innovations manifest regional variability, suggestive of distinct cultural traditions. In eastern Africa, several Late Pleistocene sites have produced evidence for novel activities, but the chronologies of key behavioral innovations remain unclear. The 3 m deep, well-dated, Panga ya Saidi sequence in eastern Kenya, encompassing 19 layers covering a time span of 78 kyr beginning in late Marine Isotope Stage 5, is the only known African site recording the interplay between cultural and ecological diversity in a coastal forested environment. Excavations have yielded worked and incised bones, ostrich eggshell beads (OES), beads made from seashells, worked and engraved ocher pieces, fragments of coral, and a belemnite fossil. Here, we provide, for the first time, a detailed analysis of this material. This includes a taphonomic, archeozoological, technological, and functional study of bone artifacts; a technological and morphometric analysis of personal ornaments; and a technological and geochemical analysis of ocher pieces. The interpretation of the results stemming from the analysis of OES beads is guided by an ethnoarcheological perspective and field observations. We demonstrate that key cultural innovations on the eastern African coast are evident by 67 ka and exhibit remarkable diversity through the LSA and Iron Age. We suggest the cultural trajectories evident at Panga ya Saidi were shaped by both regional traditions and cultural/demic diffusion.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Tracking human population structure through time from whole genome sequences

The genetic diversity of humans, like many species, has been shaped by a complex pattern of population separations followed by isolation and subsequent admixture. This pattern, reaching at least as far back as the appearance of our species in the paleontological record, has left its traces in our genomes. Reconstructing a population's history from these traces is a challenging problem. Here we present a novel approach based on the Multiple Sequentially Markovian Coalescent (MSMC) to analyze the separation history between populations. Our approach, called MSMC-IM, uses an improved implementation of the MSMC (MSMC2) to estimate coalescence rates within and across pairs of populations, and then fits a continuous Isolation-Migration model to these rates to obtain a time-dependent estimate of gene flow. We show, using simulations, that our method can identify complex demographic scenarios involving post-split admixture or archaic introgression. We apply MSMC-IM to whole genome sequences from 15 worldwide populations, tracking the process of human genetic diversification. We detect traces of extremely deep ancestry between some African populations, with around 1% of ancestry dating to divergences older than a million years ago.

Human auditory ossicles as an alternative optimal source of ancient DNA

DNA recovery from ancient human remains has revolutionized our ability to reconstruct the genetic landscape of the past. Ancient DNA research has benefited from the identification of skeletal elements, such as the cochlear part of the osseous inner ear, that provides optimal contexts for DNA preservation; however, the rich genetic information obtained from the cochlea must be counterbalanced against the loss of morphological information caused by its sampling. Motivated by similarities in developmental processes and histological properties between the cochlea and auditory ossicles, we evaluate the ossicles as an alternative source of ancient DNA. We show that ossicles perform comparably to the cochlea in terms of DNA recovery, finding no substantial reduction in data quantity and minimal differences in data quality across preservation conditions. Ossicles can be sampled from intact skulls or disarticulated petrous bones without damage to surrounding bone, and we argue that they should be used when available to reduce damage to human remains. Our results identify another optimal skeletal element for ancient DNA analysis and add to a growing toolkit of sampling methods that help to better preserve skeletal remains for future research while maximizing the likelihood that ancient DNA analysis will produce useable results.

Direct Evidence of an Increasing Mutational Load in Humans

The extent to which selection has shaped present-day human populations has attracted intense scrutiny, and examples of local adaptations abound. However, the evolutionary trajectory of alleles that, today, are deleterious has received much less attention. To address this question, the genomes of 2,062 individuals, including 1,179 ancient humans, were reanalyzed to assess how frequencies of risk alleles and their homozygosity changed through space and time in Europe over the past 45,000 years. Although the overall deleterious homozygosity has consistently decreased, risk alleles have steadily increased in frequency over that period of time. Those that increased most are associated with diseases such as asthma, Crohn disease, diabetes, and obesity, which are highly prevalent in present-day populations. These findings may not run against the existence of local adaptations but highlight the limitations imposed by drift and population dynamics on the strength of selection in purging deleterious mutations from human populations.

A western route of prehistoric human migration from Africa into the Iberian Peninsula

Being at the western fringe of Europe, Iberia had a peculiar prehistory and a complex pattern of Neolithization. A few studies, all based on modern populations, reported the presence of DNA of likely African origin in this region, generally concluding it was the result of recent gene flow, probably during the Islamic period. Here, we provide evidence of much older gene flow from Africa to Iberia by sequencing whole genomes from four human remains from northern Portugal and southern Spain dated around 4000 years BP (from the Middle Neolithic to the Bronze Age). We found one of them to carry an unequivocal sub-Saharan mitogenome of most probably West or West-Central African origin, to our knowledge never reported before in prehistoric remains outside Africa. Our analyses of ancient nuclear genomes show small but significant levels of sub-Saharan African affinity in several ancient Iberian samples, which indicates that what we detected was not an occasional individual phenomenon, but an admixture event recognizable at the population level. We interpret this result as evidence of an early migration process from Africa into the Iberian Peninsula through a western route, possibly across the Strait of Gibraltar.

Genomic landscape of the signals of positive natural selection in populations of Northern Eurasia: A view from Northern Russia

Natural selection of beneficial genetic variants played a critical role in human adaptation to a wide range of environmental conditions. Northern Eurasia, despite its severe climate, is home to lots of ethnically diverse populations. The genetic variants associated with the survival of these populations have hardly been analyzed. We searched for the genomic signatures of positive selection in (1) the genome-wide microarray data of 432 people from eight different northern Russian populations and (2) the whole-genome sequences of 250 people from Northern Eurasia from a public repository through testing the extended haplotype homozigosity (EHH) and direct comparison of allele frequency, respectively. The 20 loci with the strongest selection signals were characterized in detail. Among the top EHH hits were the NRG3 and NBEA genes, which are involved in the development and functioning of the neural system, the PTPRM gene, which mediates cell-cell interactions and adhesion, and a region on chromosome 4 (chr4:28.7-28.9 Mb) that contained several loci affiliated with different classes of non-coding RNAs (RN7SL101P, MIR4275, MESTP3, and LINC02364). NBEA and the region on chromosome 4 were novel selection targets that were identified for the first time in Western Siberian populations. Cross-population comparisons of EHH profiles suggested a particular role for the chr4:28.7-28.9 Mb region in the local adaptation of Western Siberians. The strongest selection signal identified in Siberian sequenced genomes was formed by six SNPs on chromosome 11 (chr11:124.9-125.2 Mb). This region included well-known genes SLC37A2 and PKNOX2. SLC37A2 is most-highly expressed in the gut. Its expression is regulated by vitamin D, which is often deficient in northern regions. The PKNOX2 gene is a transcription factor of the homeobox family that is expressed in the brain and many other tissues. This gene is associated with alcohol addiction, which is widespread in many Northern Eurasian populations.

Recovering signals of ghost archaic introgression in African populations

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

Identifying and Interpreting Apparent Neanderthal Ancestry in African Individuals

Admixture has played a prominent role in shaping patterns of human genomic variation, including gene flow with now-extinct hominins like Neanderthals and Denisovans. Here, we describe a novel probabilistic method called IBDmix to identify introgressed hominin sequences, which, unlike existing approaches, does not use a modern reference population. We applied IBDmix to 2,504 individuals from geographically diverse populations to identify and analyze Neanderthal sequences segregating in modern humans. Strikingly, we find that African individuals carry a stronger signal of Neanderthal ancestry than previously thought. We show that this can be explained by genuine Neanderthal ancestry due to migrations back to Africa, predominately from ancestral Europeans, and gene flow into Neanderthals from an early dispersing group of humans out of Africa. Our results refine our understanding of Neanderthal ancestry in African and non-African populations and demonstrate that remnants of Neanderthal genomes survive in every modern human population studied to date.

Genetic Affinities among Southern Africa Hunter-Gatherers and the Impact of Admixing Farmer and Herder Populations

Southern African indigenous groups, traditionally hunter-gatherers (San) and herders (Khoekhoe), are commonly referred to as “Khoe-San” populations and have a long history in southern Africa. Their ancestors were largely isolated up until ∼2,000 years ago before the arrival of pastoralists and farmers in southern Africa. Assessing relationships among regional Khoe-San groups has been challenging due to admixture with immigrant populations that obscure past population affinities and gene flow among these autochthonous communities. We re-evaluate a combined genome-wide data set of previously published southern Africa Khoe-San populations in conjunction with novel data from Khoe-San individuals collected in Xade (Central Kalahari Game Reserve, Botswana) prior to their resettlement outside the reserve. After excluding regions in the genome that trace their ancestry to recent migrant groups, the genetic diversity of 20 Khoe-San groups fitted an isolation-by-distance model. Even though isolation-by-distance explained most genetic affinities between the different autochthonous groups, additional signals of contact between Khoe-San groups could be detected. For instance, we found stronger genetic affinities, than what would be explained by isolation-by-distance gene flow, between the two geographically separated Khoe-San groups, who speak branches of the Kx’a-language family (ǂHoan and Ju). We also scanned the genome-wide data for signals of adaptive gene flow from farmers/herders into Khoe-San groups and identified a number of genomic regions potentially introduced by the arrival of the new groups. This study provides a comprehensive picture of affinities among Khoe-San groups, prior to the arrival of recent migrants, and found that these affinities are primarily determined by the geographic landscape.

A different view on fine-scale population structure in Western African populations

Due to its long genetic evolutionary history, Africans exhibit more genetic variation than any other population in the world. Their genetic diversity further lends itself to subdivisions of Africans into groups of individuals with a genetic similarity of varying degrees of granularity. It remains challenging to detect fine-scale structure in a computationally efficient and meaningful way. In this paper, we present a proof-of-concept of a novel fine-scale population structure detection tool with Western African samples. These samples consist of 1396 individuals from 25 ethnic groups (two groups are African American descendants). The strategy is based on a recently developed tool called IPCAPS. IPCAPS, or Iterative Pruning to CApture Population Structure, is a genetic divisive clustering strategy that enhances iterative pruning PCA, is robust to outliers and does not require a priori computation of haplotypes. Our strategy identified in total 12 groups and 6 groups were revealed as fine-scale structure detected in the samples from Cameroon, Gambia, Mali, Southwest USA, and Barbados. Our finding helped to explain evolutionary processes in the analyzed West African samples and raise awareness for fine-scale structure resolution when conducting genome-wide association and interaction studies.

Ancient West African foragers in the context of African population history

Our knowledge of ancient human population structure in sub-Saharan Africa, particularly prior to the advent of food production, remains limited. Here we report genome-wide DNA data from four children-two of whom were buried approximately 8,000 years ago and two 3,000 years ago-from Shum Laka (Cameroon), one of the earliest known archaeological sites within the probable homeland of the Bantu language group1-11. One individual carried the deeply divergent Y chromosome haplogroup A00, which today is found almost exclusively in the same region12,13. However, the genome-wide ancestry profiles of all four individuals are most similar to those of present-day hunter-gatherers from western Central Africa, which implies that populations in western Cameroon today-as well as speakers of Bantu languages from across the continent-are not descended substantially from the population represented by these four people. We infer an Africa-wide phylogeny that features widespread admixture and three prominent radiations, including one that gave rise to at least four major lineages deep in the history of modern humans.

HOPS: automated detection and authentication of pathogen DNA in archaeological remains

High-throughput DNA sequencing enables large-scale metagenomic analyses of complex biological systems. Such analyses are not restricted to present-day samples and can also be applied to molecular data from archaeological remains. Investigations of ancient microbes can provide valuable information on past bacterial commensals and pathogens, but their molecular detection remains a challenge. Here, we present HOPS (Heuristic Operations for Pathogen Screening), an automated bacterial screening pipeline for ancient DNA sequences that provides detailed information on species identification and authenticity. HOPS is a versatile tool for high-throughput screening of DNA from archaeological material to identify candidates for genome-level analyses.

West Asian sources of the Eurasian component in Ethiopians: a reassessment

The presence of genomic signatures of Eurasian origin in contemporary Ethiopians has been reported by several authors and estimated to have arrived in the area from 3000 years ago. Several studies reported plausible source populations for such a signature, using haplotype based methods on modern data or single-site methods on modern or ancient data. These studies did not reach a consensus and suggested an Anatolian or Sardinia-like proxy, broadly Levantine or Neolithic Levantine as possible sources. We demonstrate, however, that the deeply divergent, autochthonous African component which accounts for ~50% of most contemporary Ethiopian genomes, affects the overall allele frequency spectrum to an extent that makes it hard to control for it and, at once, to discern between subtly different, yet important, Eurasian sources (such as Anatolian or Levant Neolithic ones). Here we re-assess pattern of allele sharing between the Eurasian component of Ethiopians (here called “NAF” for Non African) and ancient and modern proxies. Our results unveil a genomic legacy that may connect the Eurasian genetic component of contemporary Ethiopians with Sea People and with population movements that affected the Mediterranean area and the Levant after the fall of the Minoan civilization.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.

Ostrich eggshell bead diameter in the Holocene: Regional variation with the spread of herding in eastern and southern Africa

Despite their ubiquity in Holocene African archaeological assemblages, ostrich eggshell (OES) beads are rarely studied in detail. An exception is in southern Africa, where there is a proposed relationship between OES bead diameter and the arrival of herding ~2000 years before present. In 1987, Leon Jacobson first observed that beads from forager sites in Namibia tended to be smaller than those associated with herder sites. Studies examining bead size around the Western Cape have generally confirmed Jacobson's findings, though the driving forces of the diameter change remain unknown. Since this time, diameter has become an informal way of distinguishing forager and herder assemblages in southern Africa, but no large-scale studies of OES bead variation have been undertaken. Here we present an expanded analysis of Holocene OES bead diameters from southern, and for the first time, eastern Africa. Results reveal distinct patterns in OES bead size over time, reflecting different local dynamics associated with the spread of herding. In southern Africa, OES diameters display low variability and smaller absolute size through time. While larger beads begin to appear <2000 years ago, most beads in our study remained smaller. In contrast, eastern African OES bead diameters are consistently larger over the last 10,000 years and show no appreciable size change with the introduction of herding. Notably, larger beads thought to be associated with herders in southern Africa fall within the range of eastern African beads, indicating a potential connection between these regions in the Late Holocene consistent with genetic findings. Regional differences in bead size are subtle, on the order of millimeters, yet offer a potentially important line of evidence for investigating the spread of herding in sub-Saharan Africa. In order to understand the meaning of these changes, we encourage additional studies of OES bead assemblages and urge researchers to report individual bead diameters, rather than averages by level.

Archaic hominin introgression into modern human genomes

Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?

Human origins in a southern African palaeo-wetland and first migrations

Anatomically modern humans originated in Africa around 200 thousand years ago (ka)^1-4. Although some of the oldest skeletal remains suggest an eastern African origin², southern Africa is home to contemporary populations that represent the earliest branch of human genetic phylogeny^5,6. Here we generate, to our knowledge, the largest resource for the poorly represented and deepest-rooting maternal L0 mitochondrial DNA branch (198 new mitogenomes for a total of 1,217 mitogenomes) from contemporary southern Africans and show the geographical isolation of L0d1'2, L0k and L0g KhoeSan descendants south of the Zambezi river in Africa. By establishing mitogenomic timelines, frequencies and dispersals, we show that the L0 lineage emerged within the residual Makgadikgadi-Okavango palaeo-wetland of southern Africa⁷, approximately 200 ka (95% confidence interval, 240-165 ka). Genetic divergence points to a sustained 70,000-year-long existence of the L0 lineage before an out-of-homeland northeast-southwest dispersal between 130 and 110 ka. Palaeo-climate proxy and model data suggest that increased humidity opened green corridors, first to the northeast then to the southwest. Subsequent drying of the homeland corresponds to a sustained effective population size (L0k), whereas wet-dry cycles and probable adaptation to marine foraging allowed the southwestern migrants to achieve population growth (L0d1'2), as supported by extensive south-coastal archaeological evidence^8-10. Taken together, we propose a southern African origin of anatomically modern humans with sustained homeland occupation before the first migrations of people that appear to have been driven by regional climate changes.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

Identifying loci under positive selection in complex population histories

Detailed modeling of a species' history is of prime importance for understanding how natural selection operates over time. Most methods designed to detect positive selection along sequenced genomes, however, use simplified representations of past histories as null models of genetic drift. Here, we present the first method that can detect signatures of strong local adaptation across the genome using arbitrarily complex admixture graphs, which are typically used to describe the history of past divergence and admixture events among any number of populations. The method-called graph-aware retrieval of selective sweeps (GRoSS)-has good power to detect loci in the genome with strong evidence for past selective sweeps and can also identify which branch of the graph was most affected by the sweep. As evidence of its utility, we apply the method to bovine, codfish, and human population genomic data containing panels of multiple populations related in complex ways. We find new candidate genes for important adaptive functions, including immunity and metabolism in understudied human populations, as well as muscle mass, milk production, and tameness in specific bovine breeds. We are also able to pinpoint the emergence of large regions of differentiation owing to inversions in the history of Atlantic codfish.

Ancient DNA from the skeletons of Roopkund Lake reveals Mediterranean migrants in India

Situated at over 5,000 meters above sea level in the Himalayan Mountains, Roopkund Lake is home to the scattered skeletal remains of several hundred individuals of unknown origin. We report genome-wide ancient DNA for 38 skeletons from Roopkund Lake, and find that they cluster into three distinct groups. A group of 23 individuals have ancestry that falls within the range of variation of present-day South Asians. A further 14 have ancestry typical of the eastern Mediterranean. We also identify one individual with Southeast Asian-related ancestry. Radiocarbon dating indicates that these remains were not deposited simultaneously. Instead, all of the individuals with South Asian-related ancestry date to ~800 CE (but with evidence of being deposited in more than one event), while all other individuals date to ~1800 CE. These differences are also reflected in stable isotope measurements, which reveal a distinct dietary profile for the two main groups.

Remains of several hundred humans are scattered around Roopkund Lake, situated over 5,000 meters above sea level in the Himalayan Mountains. Here the authors analyze genome-wide data from 38 skeletons and find 3 clusters with different ancestries and dates, showing the people were desposited in multiple catastrophic events.

The genomic history of the Iberian Peninsula over the past 8000 years

We assembled genome-wide data from 271 ancient Iberians, of whom 176 are from the largely unsampled period after 2000 BCE, thereby providing a high-resolution time transect of the Iberian Peninsula. We document high genetic substructure between northwestern and southeastern hunter-gatherers before the spread of farming. We reveal sporadic contacts between Iberia and North Africa by ~2500 BCE and, by ~2000 BCE, the replacement of 40% of Iberia's ancestry and nearly 100% of its Y-chromosomes by people with Steppe ancestry. We show that, in the Iron Age, Steppe ancestry had spread not only into Indo-European-speaking regions but also into non-Indo-European-speaking ones, and we reveal that present-day Basques are best described as a typical Iron Age population without the admixture events that later affected the rest of Iberia. Additionally, we document how, beginning at least in the Roman period, the ancestry of the peninsula was transformed by gene flow from North Africa and the eastern Mediterranean.

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.

Diet and adult age-at-death among mobile foragers: A synthesis of bioarcheological methods

OBJECTIVES

The research explores whether the combined study of cortical bone histology, bone morphology, and dietary stable isotopes can expand insights into past human health and adaptations, particularly dietary sufficiency and life span.

MATERIALS AND METHODS

Midthoracic rib cortices from 54 South African Late Holocene adult skeletons (28 M, 24 F, two sex undetermined) are assessed by transmitted-light microscopy for cross-sectional area measurements, osteon area (On.Ar), osteon population density, and presence/absence of secondary osteon variants. Values for δ¹³ C_{bone collagen} , δ¹⁵ N_{bone collagen} , ¹⁴ C dates, Southwestern and Southern Cape geographic regions, body size measures, estimated ages-at-death from both morphological and histological methods are integrated into analyses, which include Spearman correlations, χ² tests and Kruskal-Wallis ANOVAs.

RESULTS

There is reduced On.Ar variability with higher δ¹⁵ N (r = -.41, p = .005); rib %cortical area and δ¹⁵ N are negatively correlated in the Southern Cape group (r = -.60, p = .03). Osteon variants are more common in older adults; histological ages at death are significantly older than those determined from gross morphology.

DISCUSSION

We found bone tissue relationships with measures of diet composition, but indicators of dietary adequacy remain elusive. Relationships of tissue quality and isotopes suggest that some Southern Cape adults lived long lives. Osteon variants are associated with age-at-death; some association with diet remains possible. Gross morphological methods appear to underestimate adult ages-at-death, at least among small-bodied adults.

Genome-Wide Characterization of Arabian Peninsula Populations: Shedding Light on the History of a Fundamental Bridge between Continents

The Arabian Peninsula (AP) was an important crossroad between Africa, Asia, and Europe, being the cradle of the structure defining these main human population groups, and a continuing path for their admixture. The screening of 741,000 variants in 420 Arabians and 80 Iranians allowed us to quantify the dominant sub-Saharan African admixture in the west of the peninsula, whereas South Asian and Levantine/European influence was stronger in the east, leading to a rift between western and eastern sides of the Peninsula. Dating of the admixture events indicated that Indian Ocean slave trade and Islamization periods were important moments in the genetic makeup of the region. The western-eastern axis was also observable in terms of positive selection of diversity conferring lactose tolerance, with the West AP developing local adaptation and the East AP acquiring the derived allele selected in European populations and existing in South Asia. African selected malaria resistance through the DARC gene was enriched in all Arabian genomes, especially in the western part. Clear European influences associated with skin and eye color were equally frequent across the Peninsula.

PopNetD3-A Network-Based Web Resource for Exploring Population Structure

We present PopNetD3, a web tool that provides an integrated approach for the network-based visualization of population structure based on the PopNet clustering framework. Users first submit a tab-delimited file that defines diversity of SNPs across the genome which is subsequently processed by the PopNet backend to define patterns of conservation at the chromosome level. The resulting population structure is visualized through a dedicated D3-based tool, allowing users to interactively examine chromosomal regions predicted to share ancestry. We illustrate the capabilities of PopNetD3 through an analysis of 16 strains of Neisseria gonorrhoeae. PopNetD3 is capable of processing population data sets consisting of hundreds of individuals and is publicly available online at: http://compsysbio.org/popnetd3 Last Accessed: May 17, 2019.

Genomics of disease risk in globally diverse populations

Risk of disease is multifactorial and can be shaped by socio-economic, demographic, cultural, environmental and genetic factors. Our understanding of the genetic determinants of disease risk has greatly advanced with the advent of genome-wide association studies (GWAS), which detect associations between genetic variants and complex traits or diseases by comparing populations of cases and controls. However, much of this discovery has occurred through GWAS of individuals of European ancestry, with limited representation of other populations, including from Africa, The Americas, Asia and Oceania. Population demography, genetic drift and adaptation to environments over thousands of years have led globally to the diversification of populations. This global genomic diversity can provide new opportunities for discovery and translation into therapies, as well as a better understanding of population disease risk. Large-scale multi-ethnic and representative biobanks and population health resources provide unprecedented opportunities to understand the genetic determinants of disease on a global scale.

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.

Models of archaic admixture and recent history from two-locus statistics

We learn about population history and underlying evolutionary biology through patterns of genetic polymorphism. Many approaches to reconstruct evolutionary histories focus on a limited number of informative statistics describing distributions of allele frequencies or patterns of linkage disequilibrium. We show that many commonly used statistics are part of a broad family of two-locus moments whose expectation can be computed jointly and rapidly under a wide range of scenarios, including complex multi-population demographies with continuous migration and admixture events. A full inspection of these statistics reveals that widely used models of human history fail to predict simple patterns of linkage disequilibrium. To jointly capture the information contained in classical and novel statistics, we implemented a tractable likelihood-based inference framework for demographic history. Using this approach, we show that human evolutionary models that include archaic admixture in Africa, Asia, and Europe provide a much better description of patterns of genetic diversity across the human genome. We estimate that an unidentified, deeply diverged population admixed with modern humans within Africa both before and after the split of African and Eurasian populations, contributing 4 - 8% genetic ancestry to individuals in world-wide populations.

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.