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.

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.

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.

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.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics.

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation.

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500-6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east-west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

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.

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.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

Ancient substructure in early mtDNA lineages of southern Africa.

Among the deepest-rooting clades in the human mitochondrial DNA (mtDNA) phylogeny are the haplogroups defined as L0d and L0k, which are found primarily in southern Africa. These lineages are typically present at high frequency in the so-called Khoisan populations of hunter-gatherers and herders who speak non-Bantu languages, and the early divergence of these lineages led to the hypothesis of ancient genetic substructure in Africa. Here we update the phylogeny of the basal haplogroups L0d and L0k with 500 full mtDNA genome sequences from 45 southern African Khoisan and Bantu-speaking populations. We find previously unreported subhaplogroups and greatly extend the amount of variation and time-depth of most of the known subhaplogroups. Our major finding is the definition of two ancient sublineages of L0k (L0k1b and L0k2) that are present almost exclusively in Bantu-speaking populations from Zambia; the presence of such relic haplogroups in Bantu speakers is most probably due to contact with ancestral pre-Bantu populations that harbored different lineages than those found in extant Khoisan. We suggest that although these populations went extinct after the immigration of the Bantu-speaking populations, some traces of their haplogroup composition survived through incorporation into the gene pool of the immigrants. Our findings thus provide evidence for deep genetic substructure in southern Africa prior to the Bantu expansion that is not represented in extant Khoisan populations.

Cervical necrotizing fasciitis of nonodontogenic origin: case report and review of literature.

Cervical necrotizing fasciitis (CNF) is a potentially fatal infection characterized by generalized necrosis of the cervical fascia that progresses rapidly. The incidence of this entity corresponds to 2.6% of all infections of the head and neck. The most frequent primary origin is dental infection, although other causes exist that should be evaluated.Delay in the diagnosis of this entity may lead to rapid progression and fatal outcome. Patients often present immunosuppression or systemic diseases that predispose them to this pathology. Cervical necrotizing fasciitis is associated with mortality rates of 7% to 20% depending on the extension of the cervical lesion. The highest rates correspond to cases that progress to mediastinitis or septic shock, which are the main and most frequent complications. Early detection and adequate emergency treatment are critical in the management of these patients and may reduce morbimortality and improve survival. The emergency services should be prepared to manage such cases efficiently, through a multidisciplinary treatment by coordinating emergency surgery with critical support and clinical stabilization of patients.We present a case of CNF of non odontogenic origin managed in our hospital.

Flowering Coriander (Coriandrum sativum) Strips Do Not Enhance Ecosystem Services in Azorean Orchards.

The effect of flower strips on ecosystem services (ESs) and disservices (EDs) is routinely assessed following changes in service provider densities without measuring the associated levels of ES/EDs. By using the sentinel approach (i.e., exposing a plant, seeds, and prey models in a standardized way), we tested how coriander (Coriandrum sativum) strips planted in mixed orchards on Terceira Island (Azores, Portugal) affected herbivory on lettuce plants, seed predation on wheat and weed seeds, and predation on artificial caterpillars. Vertebrates had more influence than invertebrates on ESs/EDs. Herbivory (ED) after 2 weeks was similar in the coriander and the control plots (mean ± SD; 2.3% ± 3.3% vs. 2.2% ± 2.9%, n = 32 for both). Seed predation was higher in the control than in the coriander plots for both grain (ED; 30.8% ± 38.9% vs. 15.3% ± 10.8%, n = 18 for both) and weed seeds (ES; 2.5% ± 4.1% vs. 0.4% ± 0.5%, n = 18 for both). Vertebrate predation (ES) rates after 48 h were significantly higher in the control (estimate 9%, 95% CI: 4-20%) than in the coriander plots (3%, 1-8%), while no difference was observed for invertebrate predation. Coriander strips did not support increased ES/reduced ED levels in this setting. The tools used can be effective to quantitatively compare multiple ESs/EDs under different farming management strategies.

aMeta: an accurate and memory-efficient ancient metagenomic profiling workflow.

Analysis of microbial data from archaeological samples is a growing field with great potential for understanding ancient environments, lifestyles, and diseases. However, high error rates have been a challenge in ancient metagenomics, and the availability of computational frameworks that meet the demands of the field is limited. Here, we propose aMeta, an accurate metagenomic profiling workflow for ancient DNA designed to minimize the amount of false discoveries and computer memory requirements. Using simulated data, we benchmark aMeta against a current state-of-the-art workflow and demonstrate its superiority in microbial detection and authentication, as well as substantially lower usage of computer memory.

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.

CA 15-3 prognostic biomarker in SARS-CoV-2 pneumonia.

The severity of lung involvement is the main prognostic factor in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Carbohydrate antigen 15-3 (CA 15-3), a marker of lung damage and fibrosis, could help predict the prognosis of SARS-CoV-2 pneumonia. This was a retrospective and observational study. CA 15-3 was analyzed in the blood samples of patients consecutively admitted for SARS-CoV-2 pneumonia and whose blood samples were available in the biobank. Other prognostic markers were also measured (interleukin 6 [IL6], C-reactive protein [CRP], D-dimer, troponin T, and NT-ProBNP). The occurrence of in-hospital complications was registered, including death, the need for medical intensive care, and oxygen therapy at discharge. In this study, 539 patients were recruited (54.9% men, mean age: 59.6 ± 16.4 years). At admission, the mean concentrations of CA 15-3 was 20.5 ± 15.8 U/mL, and the concentration was correlated with male sex, older age, and other severity markers of coronavirus disease of 2019 (COVID-19) (IL6, CRP, D-dimer, troponine T, and NT-ProBNP). CA 15-3 levels were higher in patients who died (n = 56, 10.4%) (35.33 ± 30.45 vs. 18.8 ± 12.11, p < 0.001), who required intensive medical support (n = 78, 14.4%; 31.17 ± 27.83 vs. 18.68 ± 11.83; p < 0.001), and who were discharged with supplemental oxygen (n = 64, 13.3%; 22.65 ± 14.41 vs. 18.2 ± 11.7; p = 0.011). Elevated CA 15-3 levels (above 34.5 U/mL) were a strong predictor of a complicated in-hospital course, in terms of a higher risk of death (adjusted odds ratio [OR] 3.74, 95% confidence interval [CI]: 1.22-11.9, p = 0.022) and need for intensive care (adjusted OR 4.56, 95% CI: 1.37-15.8) after adjusting for all other risk factors. The degree of lung damage and fibrosis evaluated in terms of CA 15-3 concentrations may allow early identification of the increased risk of complications in patients with SARS-CoV-2 pneumonia.

The performance of common SNP arrays in assigning African mitochondrial haplogroups.

BACKGROUND: Mitochondrial haplogroup assignment is an important tool for forensics and evolutionary genetics. African populations are known to display a high diversity of mitochondrial haplogroups. In this research we explored mitochondrial haplogroup assignment in African populations using commonly used genome-wide SNP arrays. RESULTS: We show that, from eight commonly used SNP arrays, two SNP arrays outperform the other arrays when it comes to the correct assignment of African mitochondrial haplogroups. One array enables the recognition of 81% of the African mitochondrial haplogroups from our compiled dataset of full mitochondrial sequences. Other SNP arrays were able to assign 4-62% of the African mitochondrial haplogroups present in our dataset. We also assessed the performance of available software for assigning mitochondrial haplogroups from SNP array data. CONCLUSIONS: These results provide the first cross-checked quantification of mitochondrial haplogroup assignment performance from SNP array data. Mitochondrial haplogroup frequencies inferred from most common SNP arrays used for human population analysis should be considered with caution.

Philippine Ayta possess the highest level of Denisovan ancestry in the world.

Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ∼2.3 million genotypes from 118 ethnic groups of the Philippines, including 25 diverse self-identified Negrito populations, along with high-coverage genomes of Australopapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world-∼30%-40% greater than that of Australians and Papuans-consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in 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.

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.