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.

Eurasian back-migration into Northeast Africa was a complex and multifaceted process.

Recent studies have identified Northeast Africa as an important area for human movements during the Holocene. Eurasian populations have moved back into Northeastern Africa and contributed to the genetic composition of its people. By gathering the largest reference dataset to date of Northeast, North, and East African as well as Middle Eastern populations, we give new depth to our knowledge of Northeast African demographic history. By employing local ancestry methods, we isolated the Non-African parts of modern-day Northeast African genomes and identified the best putative source populations. Egyptians and Sudanese Copts bore most similarities to Levantine populations whilst other populations in the region generally had predominantly genetic contributions from the Arabian peninsula rather than Levantine populations for their Non-African genetic component. We also date admixture events and investigated which factors influenced the date of admixture and find that major linguistic families were associated with the date of Eurasian admixture. Taken as a whole we detect complex patterns of admixture and diverse origins of Eurasian admixture in Northeast African populations of today.

Genetic structure and characteristics of Tibetan chickens.

Tibetan chicken is one of the most common and widely distributed highland breeds, and is often used as a model organism for understanding genetic adaptation to extreme environments in Tibet. Despite its apparent geographical diversity and large variations in plumage patterns, the genetic differences within breed were not accounted for in most studies and have not been systematically investigated. In order to reveal and genetically differentiate the current existing TBC sub-populations that might have major implications for genomic research in TBCs, we systematically evaluated the population structure and demography of current TBC populations. Based on 344 whole-genome sequenced birds including 115 Tibetan chickens that were mostly sampled from family-farms across Tibet, we revealed a clear separation of Tibetan chickens into 4 sub-populations that broadly aligns with their geographical distribution. Moreover, population structure, population size dynamics, and the extent of admixture jointly suggest complex demographic histories of these sub-populations, including possible multiple origins, inbreeding, and introgressions. While most of the candidate selected regions found between the TBC sub-populations and Red Jungle fowls were nonoverlapping, 2 genes RYR2 and CAMK2D were revealed as strong selection candidates in all 4 sub-populations. These 2 previously identified high altitude associated genes indicated that the sub-populations responded to similar selection pressures in an independent but functionally similar fashion. Our results demonstrate robust population structure in Tibetan chickens that will help inform future genetic analyses on chickens and other domestic animals alike in Tibet, recommending thoughtful experimental design.

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

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

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.

Human adaptation to arsenic in Bolivians living in the Andes.

Humans living in the Andes Mountains have been historically exposed to arsenic from natural sources, including drinking water. Enzymatic methylation of arsenic allows it to be excreted more efficiently by the human body. Adaptation to high-arsenic environments via enhanced methylation and excretion of arsenic was first reported in indigenous women in the Argentinean Andes, but whether adaptation to arsenic is a general phenomenon across native populations from the Andes Mountains remains unclear. Therefore, we evaluated whether adaptation to arsenic has occurred in the Bolivian Andes by studying indigenous groups who belong to the Aymara-Quechua and Uru ethnicities and have lived in the Bolivian Andes for generations. Our population genetics methods, including genome-wide selection scans based on linkage disequilibrium patterns and allele frequency differences, in combination with targeted and whole-genome sequencing and genotype-phenotype association analyses, detected signatures of positive selection near the gene encoding arsenite methyltransferase (AS3MT), the main arsenic methylating enzyme. This was among the strongest selection signals (top 0.5% signals via locus-specific branch length and extended haplotype homozygosity tests) at a genome-wide level in the Bolivian study groups. We found a large haplotype block of 676 kb in the AS3MT region and identified candidate functional variants for further analysis. Moreover, our analyses revealed associations between AS3MT variants and the fraction of mono-methylated arsenic in urine and showed that the Bolivian study groups had the highest frequency of alleles associated with more efficient arsenic metabolism reported so far. Our data support the idea that arsenic exposure has been a driver for human adaptation to tolerate arsenic through more efficient arsenic detoxification in different Andean populations.

Genomics: Testing the limits of de-extinction.

Resurrecting extinct species through de-extinction by genome editing requires full and unbiased information from the extinct species' genome. A new study establishes a framework to assess how much of an extinct species genome can be recovered by ancient DNA sequencing and which factors influence recovery.

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.

Genetic substructure and complex demographic history of South African Bantu speakers.

South Eastern Bantu-speaking (SEB) groups constitute more than 80% of the population in South Africa. Despite clear linguistic and geographic diversity, the genetic differences between these groups have not been systematically investigated. Based on genome-wide data of over 5000 individuals, representing eight major SEB groups, we provide strong evidence for fine-scale population structure that broadly aligns with geographic distribution and is also congruent with linguistic phylogeny (separation of Nguni, Sotho-Tswana and Tsonga speakers). Although differential Khoe-San admixture plays a key role, the structure persists after Khoe-San ancestry-masking. The timing of admixture, levels of sex-biased gene flow and population size dynamics also highlight differences in the demographic histories of individual groups. The comparisons with five Iron Age farmer genomes further support genetic continuity over ~400 years in certain regions of the country. Simulated trait genome-wide association studies further show that the observed population structure could have major implications for biomedical genomics research in South Africa.

The Genetic Variation of Lactase Persistence Alleles in Sudan and South Sudan.

Lactase persistence (LP) is a well-studied example of a Mendelian trait under selection in some human groups due to gene-culture coevolution. We investigated the frequencies of genetic variants linked to LP in Sudanese and South Sudanese populations. These populations have diverse subsistence patterns, and some are dependent on milk to various extents, not only from cows but also from other livestock such as camels and goats. We sequenced a 316-bp region involved in regulating the expression of the LCT gene on chromosome 2, which encompasses five polymorphisms that have been associated with LP. Pastoralist populations showed a higher frequency of LP-associated alleles compared with nonpastoralist groups, hinting at positive selection also among northeast African pastoralists. Among the LP variants, the -14009:G variant occurs at the highest frequency among the investigated populations, followed by the -13915:G variant, which is likely of Middle Eastern origin, consistent with Middle Eastern gene flow to the Sudanese populations. There was no incidence of the "East African" LP allele (-14010:C) in the Sudanese and South Sudanese groups, and only one heterozygous individual for the "European" LP allele (-13910:T), suggesting limited recent admixture from these geographic regions. The Beja population of the Beni Amer show three different LP variants at substantial and similar levels, resulting in one of the greatest aggregation of LP variants among all populations across the world.

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.

Population collapse in Congo rainforest from 400 CE urges reassessment of the Bantu Expansion.

The present-day distribution of Bantu languages is commonly thought to reflect the early stages of the Bantu Expansion, the greatest migration event in African prehistory. Using 1149 radiocarbon dates linked to 115 pottery styles recovered from 726 sites throughout the Congo rainforest and adjacent areas, we show that this is not the case. Two periods of more intense human activity, each consisting of an expansion phase with widespread pottery styles and a regionalization phase with many more local pottery styles, are separated by a widespread population collapse between 400 and 600 CE followed by major resettlement centuries later. Coinciding with wetter climatic conditions, the collapse was possibly promoted by a prolonged epidemic. Comparison of our data with genetic and linguistic evidence further supports a spread-over-spread model for the dispersal of Bantu speakers and their languages.

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.

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.

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.

African population history: an ancient DNA perspective.

The history of human populations in Africa is complex and includes various demographic events that influenced patterns of genetic variation across the continent. Through genetic studies of modern-day, and most recently, ancient African genetic variation, it became evident that deep African history is captured by the relationships among hunter-gatherers. Furthermore, it was shown that agriculture had a large influence on the distribution of current-day Africans. These later population movements changed the demographic face of the continent and descendants of farming groups today form the majority populations across Africa. Ancient DNA methods are continually evolving, and we see evidence of this in how research has advanced in the last decade. With the increased availability of full genomic data from diverse sets of modern-day and prehistoric Africans we now have more power to infer human demography. Future ancient DNA research promises to reveal more detailed stories of human prehistory in 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.

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.

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.