Chiara Batini

Insights into the Demographic History of African Pygmies from Complete Mitochondrial Genomes

Pygmy populations are among the few hunter-gatherers currently living in sub-Saharan Africa and are mainly represented by two groups, Eastern and Western, according to their current geographical distribution. They are scattered across the Central African belt and surrounded by Bantu-speaking farmers, with whom they have complex social and economic interactions. To investigate the demographic history of Pygmy groups, a population approach was applied to the analysis of 205 complete mitochondrial DNA (mtDNA) sequences from ten central African populations. No sharing of maternal lineages was observed between the two Pygmy groups, with haplogroup L1c being characteristic of the Western group but most of Eastern Pygmy lineages falling into subclades of L0a, L2a, and L5. Demographic inferences based on Bayesian coalescent simulations point to an early split among the maternal ancestors of Pygmies and those of Bantu-speaking farmers (∼ 70,000 years ago [ya]). Evidence for population growth in the ancestors of Bantu-speaking farmers has been observed, starting ∼ 65,000 ya, well before the diffusion of Bantu languages. Subsequently, the effective population size of the ancestors of Pygmies remained constant over time and ∼ 27,000 ya, coincident with the Last Glacial Maximum, Eastern and Western Pygmies diverged, with evidence of subsequent migration only among the Western group and the Bantu-speaking farmers. Western Pygmies show signs of a recent bottleneck 4,000-650 ya, coincident with the diffusion of Bantu languages, whereas Eastern Pygmies seem to have experienced a more ancient decrease in population size (20,000-4,000 ya). In conclusion, the results of this first attempt at analyzing complete mtDNA sequences at the population level in sub-Saharan Africa not only support previous findings but also offer new insights into the demographic history of Pygmy populations, shedding new light on the ancient peopling of the African continent.

The Y-Chromosome Tree Bursts into Leaf: 13,000 High-Confidence SNPs Covering the Majority of Known Clades

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Large-scale recent expansion of European patrilineages shown by population resequencing

The proportion of Europeans descending from Neolithic farmers ∼10 thousand years ago (KYA) or Palaeolithic hunter-gatherers has been much debated. The male-specific region of the Y chromosome (MSY) has been widely applied to this question, but unbiased estimates of diversity and time depth have been lacking. Here we show that European patrilineages underwent a recent continent-wide expansion. Resequencing of 3.7 Mb of MSY DNA in 334 males, comprising 17 European and Middle Eastern populations, defines a phylogeny containing 5,996 single-nucleotide polymorphisms. Dating indicates that three major lineages (I1, R1a and R1b), accounting for 64% of our sample, have very recent coalescent times, ranging between 3.5 and 7.3 KYA. A continuous swathe of 13/17 populations share similar histories featuring a demographic expansion starting ∼2.1–4.2 KYA. Our results are compatible with ancient MSY DNA data, and contrast with data on mitochondrial DNA, indicating a widespread male-specific phenomenon that focuses interest on the social structure of Bronze Age Europe.

Signatures of the Preagricultural Peopling Processes in Sub-Saharan Africa as Revealed by the Phylogeography of Early Y Chromosome Lineages

The study of Y chromosome variation has helped reconstruct demographic events associated with the spread of languages, agriculture, and pastoralism in sub-Saharan Africa, but little attention has been given to the early history of the continent. In order to overcome this lack of knowledge, we carried out a phylogeographic analysis of haplogroups A and B in a broad data set of sub-Saharan populations. These two lineages are particularly suitable for this objective because they are the two most deeply rooted branches of the Y chromosome genealogy. Their distribution is almost exclusively restricted to sub-Saharan Africa where their frequency peaks at 65% in groups of foragers. The combined high-resolution single nucleotide polymorphism analysis with short tandem repeats variation of their subclades reveals strong geographic and population structure for both haplogroups. This has allowed us to identify specific lineages related to regional preagricultural dynamics in different areas of sub-Saharan Africa. In addition, we observed signatures of relatively recent contact, both among Pygmies and between them and Khoisan speaker groups from southern Africa, thus contributing to the understanding of the complex evolutionary relationships among African hunter-gatherers. Finally, by revising the phylogeography of the very early human Y chromosome lineages, we have obtained support for the role of southern Africa as a sink, rather than a source, of the first migrations of modern humans from eastern and central parts of the continent. These results open new perspectives on the early history of Homo sapiens in Africa, with particular attention to areas of the continent where human fossil remains and archaeological data are scant.

The Bantu expansion revisited: a new analysis of Y chromosome variation in Central Western Africa

The current distribution of Bantu languages is commonly considered to be a consequence of a relatively recent population expansion (3–5 kya) in Central Western Africa. While there is a substantial consensus regarding the centre of origin of Bantu languages (the Benue River Valley, between South East Nigeria and Western Cameroon), the identification of the area from where the population expansion actually started, the relation between the processes leading to the spread of languages and peoples and the relevance of local migratory events remain controversial. In order to shed new light on these aspects, we studied Y chromosome variation in a broad dataset of populations encompassing Nigeria, Cameroon, Gabon and Congo. Our results evidence an evolutionary scenario which is more complex than had been previously thought, pointing to a marked differentiation of Cameroonian populations from the rest of the dataset. In fact, in contrast with the current view of Bantu speakers as a homogeneous group of populations, we observed an unexpectedly high level of interpopulation genetic heterogeneity and highlighted previously undetected diversity for lineages associated with the diffusion of Bantu languages (E1b1a (M2) sub‐branches). We also detected substantial differences in local demographic histories, which concord with the hypotheses regarding an early diffusion of Bantu languages into the forest area and a subsequent demographic expansion and migration towards eastern and western Africa.

HLA-DRB1 and -DQB1 loci in three west African ethnic groups: Genetic relationship with sub-Saharan African and European populations

The Fulani of west Africa have been shown to be less susceptible to malaria and to mount a stronger immune response to malaria than sympatric ethnic groups. The analysis of HLA diversity is useful for the assessment of the genetic distance between the Fulani and sympatric populations, which represents the necessary theoretical background for the investigation of genetic determinants of susceptibility to malaria. We assessed the polymorphism of HLA-DRB1 and -DQB1 loci and analyzed the distribution of alleles/haplotypes in Fulani, Mossi, and Rimaibé from Burkina Faso. We then investigated the genetic relationship of these three ethnic groups with other sub-Saharan African populations as well as with Europeans. We confirmed that the Fulani from Burkina Faso are genetically distinct from sympatric Mossi and Rimaibé. Furthermore the Fulani from Burkina Faso are close to those from The Gambia and, intriguingly, share the distribution of specific alleles with east African populations (Amhara and Oromo). It is noteworthy that the HLA-DRB1*04 and -DQB1*02 alleles, which are implicated in the development of several autoimmune diseases, are present at high frequency in the Fulani, suggesting their potential involvement in the enhanced immune reactivity observed in this population.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale-multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8-15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

The influence of habitats on female mobility in Central and Western Africa inferred from human mitochondrial variation

BackgroundWhen studying the genetic structure of human populations, the role of cultural factors may be difficult to ascertain due to a lack of formal models. Linguistic diversity is a typical example of such a situation. Patrilocality, on the other hand, can be integrated into a biological framework, allowing the formulation of explicit working hypotheses. The present study is based on the assumption that patrilocal traditions make the hypervariable region I of the mtDNA a valuable tool for the exploration of migratory dynamics, offering the opportunity to explore the relationships between genetic and linguistic diversity. We studied 85 Niger-Congo-speaking patrilocal populations that cover regions from Senegal to Central African Republic. A total of 4175 individuals were included in the study.ResultsBy combining a multivariate analysis aimed at investigating the population genetic structure, with a Bayesian approach used to test models and extent of migration, we were able to detect a stepping-stone migration model as the best descriptor of gene flow across the region, with the main discontinuities corresponding to forested areas.ConclusionsOur analyses highlight an aspect of the influence of habitat variation on human genetic diversity that has yet to be understood. Rather than depending simply on geographic linear distances, patterns of female genetic variation vary substantially between savannah and rainforest environments. Our findings may be explained by the effects of recent gene flow constrained by environmental factors, which superimposes on a background shaped by pre-agricultural peopling.

A multi-perspective view of genetic variation in Cameroon.

In this study, we report the genetic variation of autosomal and Y-chromosomal microsatellites in a large Cameroon population dataset (a total of 11 populations) and jointly analyze novel and previous genetic data (mitochondrial DNA and protein coding loci) taking geographic and cultural factors into consideration. The complex pattern of genetic variation of Cameroon can in part be described by contrasting two geographic areas (corresponding to the northern and southern part of the country), which differ substantially in environmental, biological, and cultural aspects. Northern Cameroon populations show a greater within- and among-group diversity, a finding that reflects the complex migratory patterns and the linguistic heterogeneity of this area. A striking reduction of Y-chromosomal genetic diversity was observed in some populations of the northern part of the country (Podokwo and Uldeme), a result that seems to be related to their demographic history rather than to sampling issues. By exploring patterns of genetic, geographic, and linguistic variation, we detect a preferential correlation between genetics and geography for mtDNA. This finding could reflect a female matrimonial mobility that is less constrained by linguistic factors than in males. Finally, we apply the island model to mitochondrial and Y-chromosomal data and obtain a female-to-male migration Nnu ratio that was more than double in the northern part of the country. The combined effect of the propensity to inter-populational admixture of females, favored by cultural contacts, and of genetic drift acting on Y-chromosomal diversity could account for the peculiar genetic pattern observed in northern Cameroon.

Inherited chromosomally integrated human herpesvirus 6 genomes are ancient, intact and potentially able to reactivate from telomeres

ABSTRACT The genomes of human herpesvirus 6A (HHV-6A) and HHV-6B have the capacity to integrate into telomeres, the essential capping structures of chromosomes that play roles in cancer and ageing. About 1% of people worldwide are carriers of chromosomally integrated HHV-6 (ciHHV-6), which is inherited as a genetic trait. Understanding the consequences of integration for the evolution of the viral genome, for the telomere, and for the risk of disease associated with carrier status is hampered by a lack of knowledge about ciHHV-6 genomes. Here, we report an analysis of 28 ciHHV-6 genomes and show that they are significantly divergent from the few modern nonintegrated HHV-6 strains for which complete sequences are currently available. In addition, ciHHV-6B genomes in Europeans are more closely related to each other than to ciHHV-6B genomes from China and Pakistan, suggesting regional variation of the trait. Remarkably, at least one group of European ciHHV-6B carriers has inherited the same ciHHV-6B genome, integrated in the same telomere allele, from a common ancestor estimated to have existed 24,500 ± 10,600 years ago. Despite the antiquity of some, and possibly most, germ line HHV-6 integrations, the majority of ciHHV-6B (95%) and ciHHV-6A (72%) genomes contain a full set of intact viral genes and therefore appear to have the capacity for viral gene expression and full reactivation. IMPORTANCE Inheritance of HHV-6A or HHV-6B integrated into a telomere occurs at a low frequency in most populations studied to date, but its characteristics are poorly understood. However, stratification of ciHHV-6 carriers in modern populations due to common ancestry is an important consideration for genome-wide association studies that aim to identify disease risks for these people. Here, we present full sequence analysis of 28 ciHHV-6 genomes and show that ciHHV-6B in many carriers with European ancestry most likely originated from ancient integration events in a small number of ancestors. We propose that ancient ancestral origins for ciHHV-6A and ciHHV-6B are also likely in other populations. Moreover, despite their antiquity, all of the ciHHV-6 genomes appear to retain the capacity to express viral genes, and most are predicted to be capable of full viral reactivation. These discoveries represent potentially important considerations in immunocompromised patients, in particular in organ transplantation and in stem cell therapy.