Genetic substructure and complex demographic history of South African Bantu speakers

9,000 years of genetic continuity in southernmost Africa demonstrated at Oakhurst rockshelter

Southern Africa has one of the longest records of fossil hominins and harbours the largest human genetic diversity in the world. Yet, despite its relevance for human origins and spread around the globe, the formation and processes of its gene pool in the past are still largely unknown. Here, we present a time transect of genome-wide sequences from nine individuals recovered from a single site in South Africa, Oakhurst Rockshelter. Spanning the whole Holocene, the ancient DNA of these individuals allows us to reconstruct the demographic trajectories of the indigenous San population and their ancestors during the last 10,000 years. We show that, in contrast to most regions around the world, the population history of southernmost Africa was not characterized by several waves of migration, replacement and admixture but by long-lasting genetic continuity from the early Holocene to the end of the Later Stone Age. Although the advent of pastoralism and farming substantially transformed the gene pool in most parts of southern Africa after 1,300 BP, we demonstrate using allele-frequency and identity-by-descent segment-based methods that the ‡Khomani San and Karretjiemense from South Africa still show direct signs of relatedness to the Oakhurst hunter-gatherers, a pattern obscured by recent, extensive non-Southern African admixture. Yet, some southern San in South Africa still preserve this ancient, Pleistocene-derived genetic signature, extending the period of genetic continuity until today.

BaTwa populations from Zambia retain ancestry of past hunter-gatherer groups

Sub-equatorial Africa is today inhabited predominantly by Bantu-speaking groups of Western African descent who brought agriculture to the Luangwa valley in eastern Zambia ~2000 years ago. Before their arrival the area was inhabited by hunter-gatherers, who in many cases were subsequently replaced, displaced or assimilated. In Zambia, we know little about the genetic affinities of these hunter-gatherers. We examine ancestry of two isolated communities in Zambia, known as BaTwa and possible descendants of recent hunter-gatherers. We genotype over two million genome-wide SNPs from two BaTwa populations (total of 80 individuals) and from three comparative farming populations to: (i) determine if the BaTwa carry genetic links to past hunter-gatherer-groups, and (ii) characterise the genetic affinities of past Zambian hunter-gatherer-groups. The BaTwa populations do harbour a hunter-gatherer-like genetic ancestry and Western African ancestry. The hunter-gatherer component is a unique local signature, intermediate between current-day Khoe-San ancestry from southern Africa and central African rainforest hunter-gatherer ancestry.

The Value Proposition of Coordinated Population Cohorts Across Africa

Building longitudinal population cohorts in Africa for coordinated research and surveillance can influence the setting of national health priorities, lead to the introduction of appropriate interventions, and provide evidence for targeted treatment, leading to better health across the continent. However, compared to cohorts from the global north, longitudinal continental African population cohorts remain scarce, are relatively small in size, and lack data complexity. As infections and noncommunicable diseases disproportionately affect Africa's approximately 1.4 billion inhabitants, African cohorts present a unique opportunity for research and surveillance. High genetic diversity in African populations and multiomic research studies, together with detailed phenotyping and clinical profiling, will be a treasure trove for discovery. The outcomes, including novel drug targets, biological pathways for disease, and gene-environment interactions, will boost precision medicine approaches, not only in Africa but across the globe.

Genetic and molecular architecture of complex traits

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

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.

Genome-wide association study meta-analysis of blood pressure traits and hypertension in sub-Saharan African populations: an AWI-Gen study

Most hypertension-related genome-wide association studies (GWASs) focus on non-African populations, despite hypertension (a major risk factor for cardiovascular disease) being highly prevalent in Africa. The AWI-Gen study GWAS meta-analysis for blood pressure (BP)-related traits (systolic and diastolic BP, pulse pressure, mean-arterial pressure and hypertension) from three sub-Saharan African geographic regions (N = 10,775), identifies two novel genome-wide significant signals (p < 5E-08): systolic BP near P2RY1 (rs77846204; intergenic variant, p = 4.95E-08) and pulse pressure near LINC01256 (rs80141533; intergenic variant, p = 1.76E-08). No genome-wide signals are detected for the AWI-Gen GWAS meta-analysis with previous African-ancestry GWASs (UK Biobank (African), Uganda Genome Resource). Suggestive signals (p < 5E-06) are observed for all traits, with 29 SNPs associating with more than one trait and several replicating known associations. Polygenic risk scores (PRSs) developed from studies on different ancestries have limited transferability, with multi-ancestry PRS providing better prediction. This study provides insights into the genetics of BP variation in African populations.

Whole genomes from Angola and Mozambique inform about the origins and dispersals of major African migrations

As the continent of origin for our species, Africa harbours the highest levels of diversity anywhere on Earth. However, many regions of Africa remain under-sampled genetically. Here we present 350 whole genomes from Angola and Mozambique belonging to ten Bantu ethnolinguistic groups, enabling the construction of a reference variation catalogue including 2.9 million novel SNPs. We investigate the emergence of Bantu speaker population structure, admixture involving migrations across sub-Saharan Africa and model the demographic histories of Angolan and Mozambican Bantu speakers. Our results bring together concordant views from genomics, archaeology, and linguistics to paint an updated view of the complexity of the Bantu Expansion. Moreover, we generate reference panels that better represents the diversity of African populations involved in the trans-Atlantic slave trade, improving imputation accuracy in African Americans and Brazilians. We anticipate that our collection of genomes will form the foundation for future African genomic healthcare initiatives.

Linking African ancestral substructure to prostate cancer health disparities

Prostate cancer (PCa) is a significant health burden in Sub-Saharan Africa, with mortality rates loosely linked to African ancestry. Yet studies aimed at identifying contributing risk factors are lacking within the continent and as such exclude for significant ancestral diversity. Here, we investigate a series of epidemiological demographic and lifestyle risk factors for 1387 men recruited as part of the multi-ethnic Southern African Prostate Cancer Study (SAPCS). We found poverty to be a decisive factor for disease grade and age at diagnosis, with other notably significant PCa associated risk factors including sexually transmitted diseases, erectile dysfunction, gynaecomastia, and vertex or complete pattern balding. Aligned with African American data, Black ethnicity showed significant risk for PCa diagnosis (OR = 1.44, 95% CI 1.05-2.00), and aggressive disease presentation (ISUP ≥ 4: OR = 2.25, 95% CI   1.49-3.40). New to this study, we demonstrate African ancestral population substructure associated PCa disparity, observing increased risk for advanced disease for the southern African Tsonga people (ISUP ≥ 4: OR = 3.43, 95% CI   1.62-7.27). Conversely, South African Coloured were less likely to be diagnosed with aggressive disease overall (ISUP ≥ 3: OR = 0.38, 95% 0.17-0.85). Understanding the basis for PCa health disparities calls for African inclusion, however, lack of available data has limited the power to begin discussions. Here, focusing on arguably the largest study of its kind for the African continent, we draw attention to the contribution of within African ancestral diversity as a contributing factor to PCa health disparities within the genetically diverse region of southern Africa.

Genome-wide association study of esophageal squamous cell cancer identifies shared and distinct risk variants in African and Chinese populations

Esophageal squamous cell carcinoma (ESCC) has a high disease burden in sub-Saharan Africa and has a very poor prognosis. Genome-wide association studies (GWASs) of ESCC in predominantly East Asian populations indicate a substantial genetic contribution to its etiology, but no genome-wide studies have been done in populations of African ancestry. Here, we report a GWAS in 1,686 African individuals with ESCC and 3,217 population-matched control individuals to investigate its genetic etiology. We identified a genome-wide-significant risk locus on chromosome 9 upstream of FAM120A (rs12379660, p = 4.58 × 10-8, odds ratio = 1.28, 95% confidence interval = 1.22-1.34), as well as a potential African-specific risk locus on chromosome 2 (rs142741123, p = 5.49 × 10-8) within MYO1B. FAM120A is a component of oxidative stress-induced survival signals, and the associated variants at the FAM120A locus co-localized with highly significant cis-eQTLs in FAM120AOS in both esophageal mucosa and esophageal muscularis tissue. A trans-ethnic meta-analysis was then performed with the African ESCC study and a Chinese ESCC study in a combined total of 3,699 ESCC-affected individuals and 5,918 control individuals, which identified three genome-wide-significant loci on chromosome 9 at FAM120A (rs12379660, pmeta = 9.36 × 10-10), chromosome 10 at PLCE1 (rs7099485, pmeta = 1.48 × 10-8), and chromosome 22 at CHEK2 (rs1033667, pmeta = 1.47 × 10-9). This indicates the existence of both shared and distinct genetic risk loci for ESCC in African and Asian populations. Our GWAS of ESCC conducted in a population of African ancestry indicates a substantial genetic contribution to ESCC risk in Africa.

The mutational profile in a South African cohort with inherited neuropathies and spastic paraplegia

Introduction

Limited diagnostics are available for inherited neuromuscular diseases (NMD) in South Africa and (excluding muscle disease) are mainly aimed at the most frequent genes underlying genetic neuropathy (GN) and spastic ataxias in Europeans. In this study, we used next-generation sequencing to screen 61 probands with GN, hereditary spastic paraplegia (HSP), and spastic ataxias for a genetic diagnosis.

Methods

After identifying four GN probands with PMP22 duplication and one spastic ataxia proband with SCA1, the remaining probands underwent whole exome (n = 26) or genome sequencing (n = 30). The curation of coding/splice region variants using gene panels was guided by allele frequencies from internal African-ancestry control genomes (n = 537) and the Clinical Genome Resource's Sequence Variant Interpretation guidelines.

Results

Of 32 GN probands, 50% had African-genetic ancestry, and 44% were solved: PMP22 (n = 4); MFN2 (n = 3); one each of MORC2, ATP1A1, ADPRHL2, GJB1, GAN, MPZ, and ATM. Of 29 HSP probands (six with predominant ataxia), 66% had African-genetic ancestry, and 48% were solved: SPG11 (n = 3); KIF1A (n = 2); and one each of SPAST, ATL1, SPG7, PCYT2, PSEN1, ATXN1, ALDH18A1, CYP7B1, and RFT1. Structural variants in SPAST, SPG11, SPG7, MFN2, MPZ, KIF5A, and GJB1 were excluded by computational prediction and manual visualisation.

Discussion

In this preliminary cohort screening panel of disease genes using WES/WGS data, we solved ~50% of cases, which is similar to diagnostic yields reported for global cohorts. However, the mutational profile among South Africans with GN and HSP differs substantially from that in the Global North.

Analysis of Structural Variants Previously Associated With ALS in Europeans Highlights Genomic Architectural Differences in Africans

Background and Objectives

Amyotrophic lateral sclerosis (ALS) is a degenerative condition of the brain and spinal cord in which protein-coding variants in known ALS disease genes explain a minority of sporadic cases. There is a growing interest in the role of noncoding structural variants (SVs) as ALS risk variants or genetic modifiers of ALS phenotype. In small European samples, specific short SV alleles in noncoding regulatory regions of SCAF4, SQSTM1, and STMN2 have been reported to be associated with ALS, and several groups have investigated the possible role of SMN1/SMN2 gene copy numbers in ALS susceptibility and clinical severity.

Methods

Using short-read whole genome sequencing (WGS) data, we investigated putative ALS-susceptibility SCAF4 (3'UTR poly-T repeat), SQSTM1 (intron 5 AAAC insertion), and STMN2 (intron 3 CA repeat) alleles in African ancestry patients with ALS and described the architecture of the SMN1/SMN2 gene region. South African cases with ALS (n = 114) were compared with ancestry-matched controls (n = 150), 1000 Genomes Project samples (n = 2,336), and H3Africa Genotyping Chip Project samples (n = 347).

Results

There was no association with previously reported SCAF4 poly-T repeat, SQSTM1 AAAC insertion, and long STMN2 CA alleles with ALS risk in South Africans (p > 0.2). Similarly, SMN1 and SMN2 gene copy numbers did not differ between South Africans with ALS and matched population controls (p > 0.9). Notably, 20% of the African samples in this study had no SMN2 gene copies, which is a higher frequency than that reported in Europeans (approximately 7%).

Discussion

We did not replicate the reported association of SCAF4, SQSTM1, and STMN2 short SVs with ALS in a small South African sample. In addition, we found no link between SMN1 and SMN2 copy numbers and susceptibility to ALS in this South African sample, which is similar to the conclusion of a recent meta-analysis of European studies. However, the SMN gene region findings in Africans replicate previous results from East and West Africa and highlight the importance of including diverse population groups in disease gene discovery efforts. The clinically relevant differences in the SMN gene architecture between African and non-African populations may affect the effectiveness of targeted SMN2 gene therapy for related diseases such as spinal muscular atrophy.

Fine-scale cultural variation reinforces genetic structure in England

OBJECTIVES

Genes and languages both contain signatures of human evolution, population movement, and demographic history. Cultural traits like language are transmitted by interactions between people, and these traits influence how people interact. In particular, if groups of people differentiate each other based on some qualities of their cultures, and if these qualities are passed to the next generation, then this differentiation can result in barriers to gene flow. Previous work finds such barriers to gene flow between groups that speak different languages, and we explore this phenomenon further: can more subtle cultural differences also produce genetic structure in a population? We focus on whether subtle, dialect-level linguistic differences in England have influenced genetic population structure, likely by affecting mating preferences.

MATERIALS AND METHODS

We analyze spatially dense linguistic and genetic data-both of which independently contain spatially structured variation in England-to examine whether the cultural differences represented by variation in English phonology colocalize with higher genetic rates of change.

RESULTS

We find that genetic variation and dialect markers have similar spatial distributions on a country-wide scale, and that throughout England, linguistic boundaries colocalize with the boundaries of genetic clusters found using fineSTRUCTURE.

DISCUSSION

This gene-language covariation, in the absence of geographic barriers that could coordinate cultural and genetic differentiation, suggests that similar social forces influenced both dialect boundaries and the genetic population structure of England.

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

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

Evolutionary Genetics and Admixture in African Populations

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

Dense sampling of ethnic groups within African countries reveals fine-scale genetic structure and extensive historical admixture

Previous studies have highlighted how African genomes have been shaped by a complex series of historical events. Despite this, genome-wide data have only been obtained from a small proportion of present-day ethnolinguistic groups. By analyzing new autosomal genetic variation data of 1333 individuals from over 150 ethnic groups from Cameroon, Republic of the Congo, Ghana, Nigeria, and Sudan, we demonstrate a previously underappreciated fine-scale level of genetic structure within these countries, for example, correlating with historical polities in western Cameroon. By comparing genetic variation patterns among populations, we infer that many northern Cameroonian and Sudanese groups share genetic links with multiple geographically disparate populations, likely resulting from long-distance migrations. In Ghana and Nigeria, we infer signatures of intermixing dated to over 2000 years ago, corresponding to reports of environmental transformations possibly related to climate change. We also infer recent intermixing signals in multiple African populations, including Congolese, that likely relate to the expansions of Bantu language-speaking peoples.

Characterization of CYP2D6 Pharmacogenetic Variation in Sub-Saharan African Populations

Cytochrome P450 2D6 (CYP2D6) is a key enzyme in drug response owing to its involvement in the metabolism of ~ 25% of clinically prescribed medications. The encoding CYP2D6 gene is highly polymorphic, and many pharmacogenetics studies have been performed worldwide to investigate the distribution of CYP2D6 star alleles (haplotypes); however, African populations have been relatively understudied to date. In this study, the distributions of CYP2D6 star alleles and predicted drug metabolizer phenotypes-derived from activity scores-were examined across multiple sub-Saharan African populations based on bioinformatics analysis of 961 high-depth whole genome sequences. This was followed by characterization of novel star alleles and suballeles in a subset of the participants via targeted high-fidelity Single-Molecule Real-Time resequencing (Pacific Biosciences). This study revealed varying frequencies of known CYP2D6 alleles and predicted phenotypes across different African ethnolinguistic groups. Twenty-seven novel CYP2D6 star alleles were predicted computationally and two of them were further validated. This study highlights the importance of studying variation in key pharmacogenes such as CYP2D6 in the African context to better understand population-specific allele frequencies. This will aid in the development of better genotyping panels and star allele detection approaches with a view toward supporting effective implementation of precision medicine strategies in Africa and across the African diaspora.

Interpreting mismatches between linguistic and genetic patterns among speakers of Tanimuka (Eastern Tukanoan) and Yukuna (Arawakan)

Northwestern Amazonia is home to a great degree of linguistic diversity, and the human societies in that region are part of complex networks of interaction that predate the arrival of Europeans. This study investigates the population and language contact dynamics between two languages found within this region, Yukuna and Tanimuka, which belong to the Arawakan and Tukanoan language families, respectively. We use evidence from linguistics, ethnohistory, ethnography and population genetics to provide new insights into the contact dynamics between these and other human groups in NWA. Our results show that the interaction between these groups intensified in the last 500 years, to the point that it is difficult to differentiate between them genetically. However, this close interaction has led to more substantial contact-induced language changes in Tanimuka than in Yukuna, consistent with a scenario of language shift and asymmetrical power relations.

Toad zoonyms mirror the linguistic and demographic history of Greece

The common toad (Bufo bufo) has been the subject of many folk tales and superstitions in Western Europe, and as a result, it is characterised by numerous common names (zoonyms). However, the zoonyms of the toad and its associated traditions have remained unexplored in the Balkans, one of Europe's linguistic hotspots. In the present study, it was attempted to fill this knowledge gap by focusing on Greece, where more than 7.700 individuals were interviewed both in the field and through online platforms, in order to document toad zoonyms from all varieties and dialects of Greek, as well as local non-Greek languages such as Arvanitika, South Slavic dialects, and Vlach. It was found that the academically unattested zoonyms of the toad provide an unmatched and previously unexplored linguistic and ethnographic tool, as they reflect the linguistic, demographic, and historical processes that shaped modern Greece. This is particularly pertinent in the 21st century, when a majority of the country's dialects and languages are in danger of imminent extinction-and some have already gone silent. Overall, the present study shows the significance of recording zoonyms of indigenous and threatened languages as excellent linguistic and ethnographic tools that safeguard our planet's ethnolinguistic diversity and enhance our understanding on how pre-industrial communities interacted with their local fauna. Furthermore, in contrast to all other European countries, which only possess one or only a few zoonyms for the toad, the Greek world boasts an unmatched 37 zoonyms, which attest to its role as a linguistic hotspot.

Coronavirus Host Genomics Study: South Africa (COVIGen-SA)

Host genetic factors are known to modify the susceptibility, severity, and outcomes of COVID-19 and vary across populations. However, continental Africans are yet to be adequately represented in such studies despite the importance of genetic factors in understanding Africa's response to the pandemic. We describe the development of a research resource for coronavirus host genomics studies in South Africa known as COVIGen-SA-a multicollaborator strategic partnership designed to provide harmonised demographic, clinical, and genetic information specific to Black South Africans with COVID-19. Over 2,000 participants have been recruited to date. Preliminary results on 1,354 SARS-CoV-2 positive participants from four participating studies showed that 64.7% were female, 333 had severe disease, and 329 were people living with HIV. Through this resource, we aim to provide insights into host genetic factors relevant to African-ancestry populations, using both genome-wide association testing and targeted sequencing of important genomic loci. This project will promote and enhance partnerships, build skills, and develop resources needed to address the COVID-19 burden and associated risk factors in South African communities.

A probable cis-acting genetic modifier of Huntington disease frequent in individuals with African ancestry

Huntington disease (HD)is a dominantly inherited neurodegenerative disorder caused by the expansion of a polyglutamine encoding CAG repeat in the huntingtin gene. Recently, it has been established that disease severity in HD is best predicted by the number of pure CAG repeats rather than total glutamines encoded. Along with uncovering DNA repair gene variants as trans-acting modifiers of HD severity, these data reveal somatic expansion of the CAG repeat as a key driver of HD onset. Using high-throughput DNA sequencing, we have determined the precise sequence and somatic expansion profiles of the HTT repeat tract of 68 HD-affected and 158 HD-unaffected African ancestry individuals. A high level of HTT repeat sequence diversity was observed, with three likely African-specific alleles identified. In the most common disease allele (30 out of 68), the typical proline-encoding CCGCCA sequence was absent. This CCGCCA-loss disease allele was associated with an earlier age of diagnosis of approximately 7.1 years and occurred exclusively on haplotype B2. Although somatic expansion was associated with an earlier age of diagnosis in the study overall, the CCGCCA-loss disease allele displayed reduced somatic expansion relative to the typical HTT expansions in blood DNA. We propose that the CCGCCA loss occurring on haplotype B2 is an African cis-acting modifier that appears to alter disease diagnosis of HD through a mechanism that is not driven by somatic expansion. The assessment of a group of individuals from an understudied population has highlighted population-specific differences that emphasize the importance of studying genetically diverse populations in the context of disease.

GWAS in the southern African context

Researchers would generally adjust for the possible confounding effect of population structure by considering global ancestry proportions or top principle components. Alternatively, researchers would conduct admixture mapping to increase the power to detect variants with an ancestry effect. This is sufficient in simple admixture scenarios, however, populations from southern Africa can be complex multi-way admixed populations. Duan et al. (2018) first described local ancestry adjusted allelic (LAAA) analysis as a robust method for discovering association signals, while producing minimal false positive hits. Their simulation study, however, was limited to a two-way admixed population. Realizing that their findings might not translate to other admixture scenarios, we simulated a three- and five-way admixed population to compare the LAAA model to other models commonly used in genome-wide association studies (GWAS). We found that, given our admixture scenarios, the LAAA model identifies the most causal variants in most of the phenotypes we tested across both the three-way and five-way admixed populations. The LAAA model also produced a high number of false positive hits which was potentially caused by the ancestry effect size that we assumed. Considering the extent to which the various models tested differed in their results and considering that the source of a given association is unknown, we recommend that researchers use multiple GWAS models when analysing populations with complex ancestry.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

A genome-wide association study for rheumatoid arthritis replicates previous HLA and non-HLA associations in a cohort from South Africa

The complex pathogenesis of rheumatoid arthritis (RA) is not fully understood, with few studies exploring the genomic contribution to RA in patients from Africa. We report a genome-wide association study (GWAS) of South-Eastern Bantu-Speaking South Africans (SEBSSAs) with seropositive RA (n = 531) and population controls (n = 2653). Association testing was performed using PLINK (logistic regression assuming an additive model) with sex, age, smoking and the first three principal components as covariates. The strong association with the Human Leukocyte Antigen (HLA) region, indexed by rs602457 (near HLA-DRB1), was replicated. An additional independent signal in the HLA region represented by the lead SNP rs2523593 (near the HLA-B gene; Conditional P-value = 6.4 × 10-10) was detected. Although none of the non-HLA signals reached genome-wide significance (P < 5 × 10-8), 17 genomic regions showed suggestive association (P < 5 × 10-6). The GWAS replicated two known non-HLA associations with MMEL1 (rs2843401) and ANKRD55 (rs7731626) at a threshold of P < 5 × 10-3 providing, for the first time, evidence for replication of non-HLA signals for RA in sub-Saharan African populations. Meta-analysis with summary statistics from an African-American cohort (CLEAR study) replicated three additional non-HLA signals (rs11571302, rs2558210 and rs2422345 around KRT18P39-NPM1P33, CTLA4-ICOS and AL645568.1, respectively). Analysis based on genomic regions (200 kb windows) further replicated previously reported non-HLA signals around PADI4, CD28 and LIMK1. Although allele frequencies were overall strongly correlated between the SEBSSA and the CLEAR cohort, we observed some differences in effect size estimates for associated loci. The study highlights the need for conducting larger association studies across diverse African populations to inform precision medicine-based approaches for RA in Africa.

Genomic analysis reveals geography rather than culture as the predominant factor shaping genetic variation in northern Kenyan human populations

OBJECTIVES

The aim of this study was to characterize the genetic relationships within and among four neighboring ethnolinguistic groups in northern Kenya in light of cultural relationships to understand the extent to which geography and culture shape patterns of genetic variation.

MATERIALS AND METHODS

We collected DNA and demographic information pertaining to aspects of social identity and heritage from 572 individuals across the Turkana, Samburu, Waso Borana, and Rendille of northern Kenya. We sampled individuals across a total of nine clans from these four groups and, additionally, three territorial sections within the Turkana and successfully genotyped 376 individuals.

RESULTS

Here we report that geography predominately shapes genetic variation within and among human groups in northern Kenya. We observed a clinal pattern of genetic variation that mirrors the overall geographic distribution of the individuals we sampled. We also found relatively higher rates of intermarriage between the Rendille and Samburu and evidence of gene flow between them that reflect these higher rates of intermarriage. Among the Turkana, we observed strong recent genetic substructuring based on territorial section affiliation. Within ethnolinguistic groups, we found that Y chromosome haplotypes do not consistently cluster by natal clan affiliation. Finally, we found that sampled populations that are geographically closer have lower genetic differentiation, and that cultural similarity does not predict genetic similarity as a whole across these northern Kenyan populations.

DISCUSSION

Overall, the results from this study highlight the importance of geography, even on a local geographic scale, in shaping observed patterns of genetic variation in human populations.

Meta-analysis of sub-Saharan African studies provides insights into genetic architecture of lipid traits

Genetic associations for lipid traits have identified hundreds of variants with clear differences across European, Asian and African studies. Based on a sub-Saharan-African GWAS for lipid traits in the population cross-sectional AWI-Gen cohort (N = 10,603) we report a novel LDL-C association in the GATB region (P-value=1.56 × 10⁻⁸). Meta-analysis with four other African cohorts (N = 23,718) provides supporting evidence for the LDL-C association with the GATB/FHIP1A region and identifies a novel triglyceride association signal close to the FHIT gene (P-value =2.66 × 10⁻⁸). Our data enable fine-mapping of several well-known lipid-trait loci including LDLR, PMFBP1 and LPA. The transferability of signals detected in two large global studies (GLGC and PAGE) consistently improves with an increase in the size of the African replication cohort. Polygenic risk score analysis shows increased predictive accuracy for LDL-C levels with the narrowing of genetic distance between the discovery dataset and our cohort. Novel discovery is enhanced with the inclusion of African data.

Genetic associations and polygenic scores for lipid traits have low transferability to African individuals. Here, the authors perform a large sub-Sarahan African lipid GWAS and find that larger datasets and better global representation in discovery GWAS help to bridge this gap.

A roadmap to increase diversity in genomic studies

Two decades ago, the sequence of the first human genome was published. Since then, advances in genome technologies have resulted in whole-genome sequencing and microarray-based genotyping of millions of human genomes. However, genetic and genomic studies are predominantly based on populations of European ancestry. As a result, the potential benefits of genomic research-including better understanding of disease etiology, early detection and diagnosis, rational drug design and improved clinical care-may elude the many underrepresented populations. Here, we describe factors that have contributed to the imbalance in representation of different populations and, leveraging our experiences in setting up genomic studies in diverse global populations, we propose a roadmap to enhancing inclusion and ensuring equal health benefits of genomics advances. Our Perspective highlights the importance of sincere, concerted global efforts toward genomic equity to ensure the benefits of genomic medicine are accessible to all.

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.

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

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

UGT1A1 regulatory variant with potential effect on efficacy of HIV and cancer drugs commonly prescribed in South Africa

Aim: Despite the high disease burden of human immunodeficiency virus (HIV) infection and colorectal cancer (CRC) in South Africa (SA), treatment-relevant pharmacogenetic variants are understudied. Materials & methods: Using publicly available genotype and gene expression data, a bioinformatic pipeline was developed to identify liver expression quantitative trait loci (eQTLs). Results: A novel cis-eQTL, rs28967009, was identified for UGT1A1, which is predicted to upregulate UGT1A1 expression thereby potentially affecting the metabolism of dolutegravir and irinotecan, which are extensively prescribed in SA for HIV and colorectal cancer treatment, respectively. Conclusion: As increased UGT1A1 expression could affect the clinical outcome of dolutegravir and irinotecan treatment by increasing drug clearance, patients with the rs28967009A variant may require increased drug doses to reach therapeutic levels or should be prescribed alternative drugs.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.