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.

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.

Multiple migrations to the Philippines during the last 50,000 years.

Island Southeast Asia has recently produced several surprises regarding human history, but the region's complex demography remains poorly understood. Here, we report ∼2.3 million genotypes from 1,028 individuals representing 115 indigenous Philippine populations and genome-sequence data from two ∼8,000-y-old individuals from Liangdao in the Taiwan Strait. We show that the Philippine islands were populated by at least five waves of human migration: initially by Northern and Southern Negritos (distantly related to Australian and Papuan groups), followed by Manobo, Sama, Papuan, and Cordilleran-related populations. The ancestors of Cordillerans diverged from indigenous peoples of Taiwan at least ∼8,000 y ago, prior to the arrival of paddy field rice agriculture in the Philippines ∼2,500 y ago, where some of their descendants remain to be the least admixed East Asian groups carrying an ancestry shared by all Austronesian-speaking populations. These observations contradict an exclusive "out-of-Taiwan" model of farming-language-people dispersal within the last four millennia for the Philippines and Island Southeast Asia. Sama-related ethnic groups of southwestern Philippines additionally experienced some minimal South Asian gene flow starting ∼1,000 y ago. Lastly, only a few lowlanders, accounting for <1% of all individuals, presented a low level of West Eurasian admixture, indicating a limited genetic legacy of Spanish colonization in the Philippines. Altogether, our findings reveal a multilayered history of the Philippines, which served as a crucial gateway for the movement of people that ultimately changed the genetic landscape of the Asia-Pacific region.

Demographic and Selection Histories of Populations Across the Sahel/Savannah Belt.

The Sahel/Savannah belt harbors diverse populations with different demographic histories and different subsistence patterns. However, populations from this large African region are notably under-represented in genomic research. To investigate the population structure and adaptation history of populations from the Sahel/Savannah space, we generated dense genome-wide genotype data of 327 individuals-comprising 14 ethnolinguistic groups, including 10 previously unsampled populations. Our results highlight fine-scale population structure and complex patterns of admixture, particularly in Fulani groups and Arabic-speaking populations. Among all studied Sahelian populations, only the Rashaayda Arabic-speaking population from eastern Sudan shows a lack of gene flow from African groups, which is consistent with the short history of this population in the African continent. They are recent migrants from Saudi Arabia with evidence of strong genetic isolation during the last few generations and a strong demographic bottleneck. This population also presents a strong selection signal in a genomic region around the CNR1 gene associated with substance dependence and chronic stress. In Western Sahelian populations, signatures of selection were detected in several other genetic regions, including pathways associated with lactase persistence, immune response, and malaria resistance. Taken together, these findings refine our current knowledge of genetic diversity, population structure, migration, admixture and adaptation of human populations in the Sahel/Savannah belt and contribute to our understanding of human history and health.

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.

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.

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.

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.

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

The southern African indigenous Khoe-San populations harbor the most divergent lineages of all living peoples. Exploring their genomes is key to understanding deep human history. We sequenced 25 full genomes from five Khoe-San populations, revealing many novel variants, that 25% of variants are unique to the Khoe-San, and that the Khoe-San group harbors the greatest level of diversity across the globe. In line with previous studies, we found several gene regions with extreme values in genome-wide scans for selection, potentially caused by natural selection in the lineage leading to Homo sapiens and more recent in time. These gene regions included immunity-, sperm-, brain-, diet-, and muscle-related genes. When accounting for recent admixture, all Khoe-San groups display genetic diversity approaching the levels in other African groups and a reduction in effective population size starting around 100,000 years ago. Hence, all human groups show a reduction in effective population size commencing around the time of the Out-of-Africa migrations, which coincides with changes in the paleoclimate records, changes that potentially impacted all humans at the time.

Tales of Human Migration, Admixture, and Selection in Africa.

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

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.

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.

Sahelian pastoralism from the perspective of variants associated with lactase persistence.

OBJECTIVES: Archeological evidence shows that first nomadic pastoralists came to the African Sahel from northeastern Sahara, where milking is reported by ~7.5 ka. A second wave of pastoralists arrived with the expansion of Arabic tribes in 7th-14th century CE. All Sahelian pastoralists depend on milk production but genetic diversity underlying their lactase persistence (LP) is poorly understood. MATERIALS AND METHODS: We investigated SNP variants associated with LP in 1,241 individuals from 29 mostly pastoralist populations in the Sahel. Then, we analyzed six SNPs in the neighboring fragment (419 kb) in the Fulani and Tuareg with the -13910*T mutation, reconstructed haplotypes, and calculated expansion age and growth rate of this variant. RESULTS: Our results reveal a geographic localization of two different LP variants in the Sahel: -13910*T west of Lake Chad (Fulani and Tuareg pastoralists) and -13915*G east of there (mostly Arabic-speaking pastoralists). We show that -13910*T has a more diversified haplotype background among the Fulani than among the Tuareg and that the age estimate for expansion of this variant among the Fulani (~8.5 ka) corresponds to introduction of cattle to the area. CONCLUSIONS: This is the first study showing that the "Eurasian" LP allele -13910*T is widespread both in northern Europe and in the Sahel; however, it is limited to pastoralists in the Sahel. Since the Fulani haplotype with -13910*T is shared with contemporary Eurasians, its origin could be in a region encompassing the Near East and northeastern Africa in a population ancestral to both Saharan pastoralists and European farmers.

Northeast African genomic variation shaped by the continuity of indigenous groups and Eurasian migrations.

Northeast Africa has a long history of human habitation, with fossil-finds from the earliest anatomically modern humans, and housing ancient civilizations. The region is also the gate-way out of Africa, as well as a portal for migration into Africa from Eurasia via the Middle East and the Arabian Peninsula. We investigate the population history of northeast Africa by genotyping ~3.9 million SNPs in 221 individuals from 18 populations sampled in Sudan and South Sudan and combine this data with published genome-wide data from surrounding areas. We find a strong genetic divide between the populations from the northeastern parts of the region (Nubians, central Arab populations, and the Beja) and populations towards the west and south (Nilotes, Darfur and Kordofan populations). This differentiation is mainly caused by a large Eurasian ancestry component of the northeast populations likely driven by migration of Middle Eastern groups followed by admixture that affected the local populations in a north-to-south succession of events. Genetic evidence points to an early admixture event in the Nubians, concurrent with historical contact between North Sudanese and Arab groups. We estimate the admixture in current-day Sudanese Arab populations to about 700 years ago, coinciding with the fall of Dongola in 1315/1316 AD, a wave of admixture that reached the Darfurian/Kordofanian populations some 400-200 years ago. In contrast to the northeastern populations, the current-day Nilotic populations from the south of the region display little or no admixture from Eurasian groups indicating long-term isolation and population continuity in these areas of northeast Africa.

Population history and genetic adaptation of the Fulani nomads: inferences from genome-wide data and the lactase persistence trait.

BACKGROUND: Human population history in the Holocene was profoundly impacted by changes in lifestyle following the invention and adoption of food-production practices. These changes triggered significant increases in population sizes and expansions over large distances. Here we investigate the population history of the Fulani, a pastoral population extending throughout the African Sahel/Savannah belt. RESULTS: Based on genome-wide analyses we propose that ancestors of the Fulani population experienced admixture between a West African group and a group carrying both European and North African ancestries. This admixture was likely coupled with newly adopted herding practices, as it resulted in signatures of genetic adaptation in contemporary Fulani genomes, including the control element of the LCT gene enabling carriers to digest lactose throughout their lives. The lactase persistence (LP) trait in the Fulani is conferred by the presence of the allele T-13910, which is also present at high frequencies in Europe. We establish that the T-13910 LP allele in Fulani individuals analysed in this study lies on a European haplotype background thus excluding parallel convergent evolution. We furthermore directly link the T-13910 haplotype with the Lactase Persistence phenotype through a Genome Wide Association study (GWAS) and identify another genomic region in the vicinity of the SPRY2 gene associated with glycaemic measurements after lactose intake. CONCLUSIONS: Our findings suggest that Eurasian admixture and the European LP allele was introduced into the Fulani through contact with a North African population/s. We furthermore confirm the link between the lactose digestion phenotype in the Fulani to the MCM6/LCT locus by reporting the first GWAS of the lactase persistence trait. We also explored other signals of recent adaptation in the Fulani and identified additional candidates for selection to adapt to herding life-styles.

Patterns of variation in cis-regulatory regions: examining evidence of purifying selection.

BACKGROUND: With only 2 % of the human genome consisting of protein coding genes, functionality across the rest of the genome has been the subject of much debate. This has gained further impetus in recent years due to a rapidly growing catalogue of genomic elements, based primarily on biochemical signatures (e.g. the ENCODE project). While the assessment of functionality is a complex task, the presence of selection acting on a genomic region is a strong indicator of importance. In this study, we apply population genetic methods to investigate signals overlaying several classes of regulatory elements. RESULTS: We disentangle signals of purifying selection acting directly on regulatory elements from the confounding factors of demography and purifying selection linked to e.g. nearby protein coding regions. We confirm the importance of regulatory regions proximal to coding sequence, while also finding differential levels of selection at distal regions. We note differences in purifying selection among transcription factor families. Signals of constraint at some genomic classes were also strongly dependent on their physical location relative to coding sequence. In addition, levels of selection efficacy across genomic classes differed between African and non-African populations. CONCLUSIONS: In order to assign a valid signal of selection to a particular class of genomic sequence, we show that it is crucial to isolate the signal by accounting for the effects of demography and linked-purifying selection. Our study highlights the intricate interplay of factors affecting signals of selection on functional elements.

Adaptation to infectious disease exposure in indigenous Southern African populations.

Genetic analyses can provide information about human evolutionary history that cannot always be gleaned from other sources. We evaluated evidence of selective pressure due to introduced infectious diseases in the genomes of two indigenous southern African San groups-the ‡Khomani who had abundant contact with other people migrating into the region and the more isolated Ju|'hoansi. We used a dual approach to test for increased selection on immune genes compared with the rest of the genome in these groups. First, we calculated summary values of statistics that measure genomic signatures of adaptation to contrast selection signatures in immune genes and all genes. Second, we located regions of the genome with extreme values of three selection statistics and examined these regions for enrichment of immune genes. We found stronger and more abundant signals of selection in immune genes in the ‡Khomani than in the Ju|'hoansi. We confirm this finding within each population to avoid effects of different demographic histories of the two populations. We identified eight immune genes that have potentially been targets of strong selection in the ‡Khomani, whereas in the Ju|'hoansi, no immune genes were found in the genomic regions with the strongest signals of selection. We suggest that the more abundant signatures of selection at immune genes in the ‡Khomani could be explained by their more frequent contact with immigrant groups, which likely led to increased exposure and adaptation to introduced infectious diseases.

Human adaptation to arsenic-rich environments.

Adaptation drives genomic changes; however, evidence of specific adaptations in humans remains limited. We found that inhabitants of the northern Argentinean Andes, an arid region where elevated arsenic concentrations in available drinking water is common, have unique arsenic metabolism, with efficient methylation and excretion of the major metabolite dimethylated arsenic and a less excretion of the highly toxic monomethylated metabolite. We genotyped women from this population for 4,301,332 single nucleotide polymorphisms (SNPs) and found a strong association between the AS3MT (arsenic [+3 oxidation state] methyltransferase) gene and mono- and dimethylated arsenic in urine, suggesting that AS3MT functions as the major gene for arsenic metabolism in humans. We found strong genetic differentiation around AS3MT in the Argentinean Andes population, compared with a highly related Peruvian population (FST = 0.014) from a region with much less environmental arsenic. Also, 13 of the 100 SNPs with the highest genome-wide Locus-Specific Branch Length occurred near AS3MT. In addition, our examination of extended haplotype homozygosity indicated a selective sweep of the Argentinean Andes population, in contrast to Peruvian and Colombian populations. Our data show that adaptation to tolerate the environmental stressor arsenic has likely driven an increase in the frequencies of protective variants of AS3MT, providing the first evidence of human adaptation to a toxic chemical.

The disappearing San of southeastern Africa and their genetic affinities.

Southern Africa was likely exclusively inhabited by San hunter-gatherers before ~2000 years ago. Around that time, East African groups assimilated with local San groups and gave rise to the Khoekhoe herders. Subsequently, Bantu-speaking farmers, arriving from the north (~1800 years ago), assimilated and displaced San and Khoekhoe groups, a process that intensified with the arrival of European colonists ~350 years ago. In contrast to the western parts of southern Africa, where several Khoe-San groups still live today, the eastern parts are largely populated by Bantu speakers and individuals of non-African descent. Only a few scattered groups with oral traditions of Khoe-San ancestry remain. Advances in genetic research open up new ways to understand the population history of southeastern Africa. We investigate the genomic variation of the remaining individuals from two South African groups with oral histories connecting them to eastern San groups, i.e., the San from Lake Chrissie and the Duma San of the uKhahlamba-Drakensberg. Using ~2.2 million genetic markers, combined with comparative published data sets, we show that the Lake Chrissie San have genetic ancestry from both Khoe-San (likely the ||Xegwi San) and Bantu speakers. Specifically, we found that the Lake Chrissie San are closely related to the current southern San groups (i.e., the Karretjie people). Duma San individuals, on the other hand, were genetically similar to southeastern Bantu speakers from South Africa. This study illustrates how genetic tools can be used to assess hypotheses about the ancestry of people who seemingly lost their historic roots, only recalling a vague oral tradition of their origin.