Evelyne Heyer

Geographic Patterns of (Genetic, Morphologic, Linguistic) Variation: How Barriers Can Be Detected by Using Monmonier’s Algorithm

When sampling locations are known, the association between genetic and geographic distances can be tested by spatial autocorrelation or regression methods. These tests give some clues to the possible shape of the genetic landscape. Nevertheless, correlation analyses fail when attempting to identify where genetic barriers exist, namely, the areas where a given variable shows an abrupt rate of change. To this end, a computational geometry approach is more suitable because it provides the locations and the directions of barriers and because it can show where geographic patterns of two or more variables are similar. In this frame we have implemented Monmoniers (1973) maximum difference algorithm in a new software package to identify genetic barriers. To provide a more realistic representation of the barriers in a genetic landscape, we implemented in the software a significance test by means of bootstrap matrices analysis. As a result, the noise associated with genetic markers can be visualized on a geographic map and the areas where genetic barriers are more robust can be identified. Moreover, this multiple matrices approach can visualize the patterns of variation associated with different markers in the same overall picture. This improved Monmoniers method is highly reliable and can be applied to nongenetic data whenever sampling locations and a distance matrix between corresponding data are available.

Genomic evidence for the Pleistocene and recent population history of Native Americans

Genetic history of Native Americans Several theories have been put forth as to the origin and timing of when Native American ancestors entered the Americas. To clarify this controversy, Raghavan et al. examined the genomic variation among ancient and modern individuals from Asia and the Americas. There is no evidence for multiple waves of entry or recurrent gene flow with Asians in northern populations. The earliest migrations occurred no earlier than 23,000 years ago from Siberian ancestors. Amerindians and Athabascans originated from a single population, splitting approximately 13,000 years ago. Science, this issue 10.1126/science.aab3884 Genetic variation within ancient and extant Native American populations informs on their migration into the Americas. INTRODUCTION The consensus view on the peopling of the Americas is that ancestors of modern Native Americans entered the Americas from Siberia via the Bering Land Bridge and that this occurred at least ~14.6 thousand years ago (ka). However, the number and timing of migrations into the Americas remain controversial, with conflicting interpretations based on anatomical and genetic evidence. RATIONALE In this study, we address four major unresolved issues regarding the Pleistocene and recent population history of Native Americans: (i) the timing of their divergence from their ancestral group, (ii) the number of migrations into the Americas, (iii) whether there was ~15,000 years of isolation of ancestral Native Americans in Beringia (Beringian Incubation Model), and (iv) whether there was post-Pleistocene survival of relict populations in the Americas related to Australo-Melanesians, as suggested by apparent differences in cranial morphologies between some early (“Paleoamerican”) remains and those of more recent Native Americans. We generated 31 high-coverage modern genomes from the Americas, Siberia, and Oceania; 23 ancient genomic sequences from the Americas dating between ~0.2 and 6 ka; and SNP chip genotype data from 79 present-day individuals belonging to 28 populations from the Americas and Siberia. The above data sets were analyzed together with published modern and ancient genomic data from worldwide populations, after masking some present-day Native Americans for recent European admixture. RESULTS Using three different methods, we determined the divergence time for all Native Americans (Athabascans and Amerindians) from their Siberian ancestors to be ~20 ka, and no earlier than ~23 ka. Furthermore, we dated the divergence between Athabascans (northern Native American branch, together with northern North American Amerindians) and southern North Americans and South and Central Americans (southern Native American branch) to be ~13 ka. Similar divergence times from East Asian populations and a divergence time between the two branches that is close in age to the earliest well-established archaeological sites in the Americas suggest that the split between the branches occurred within the Americas. We additionally found that several sequenced Holocene individuals from the Americas are related to present-day populations from the same geographical regions, implying genetic continuity of ancient and modern populations in some parts of the Americas over at least the past 8500 years. Moreover, our results suggest that there has been gene flow between some Native Americans from both North and South America and groups related to East Asians and Australo-Melanesians, the latter possibly through an East Asian route that might have included ancestors of modern Aleutian Islanders. Last, using both genomic and morphometric analyses, we found that historical Native American groups such as the Pericúes and Fuego-Patagonians were not “relicts” of Paleoamericans, and hence, our results do not support an early migration of populations directly related to Australo-Melanesians into the Americas. CONCLUSION Our results provide an upper bound of ~23 ka on the initial divergence of ancestral Native Americans from their East Asian ancestors, followed by a short isolation period of no more than ~8000 years, and subsequent entrance and spread across the Americas. The data presented are consistent with a single-migration model for all Native Americans, with later gene flow from sources related to East Asians and, indirectly, Australo-Melanesians. The single wave diversified ~13 ka, likely within the Americas, giving rise to the northern and southern branches of present-day Native Americans. Population history of present-day Native Americans. The ancestors of all Native Americans entered the Americas as a single migration wave from Siberia (purple) no earlier than ~23 ka, separate from the Inuit (green), and diversified into “northern” and “southern” Native American branches ~13 ka. There is evidence of post-divergence gene flow between some Native Americans and groups related to East Asians/Inuit and Australo-Melanesians (yellow). How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Maternal traces of deep common ancestry and asymmetric gene flow between Pygmy hunter–gatherers and Bantu-speaking farmers

Two groups of populations with completely different lifestyles—the Pygmy hunter–gatherers and the Bantu-speaking farmers—coexist in Central Africa. We investigated the origins of these two groups and the interactions between them, by analyzing mtDNA variation in 1,404 individuals from 20 farming populations and 9 Pygmy populations from Central Africa, with the aim of shedding light on one of the most fascinating cultural transitions in human evolution (the transition from hunting and gathering to agriculture). Our data indicate that this region was colonized gradually, with an initial L1c-rich ancestral population ultimately giving rise to current-day farmers, who display various L1c clades, and to Pygmies, in whom L1c1a is the only surviving clade. Detailed phylogenetic analysis of complete mtDNA sequences for L1c1a showed this clade to be autochthonous to Central Africa, with its most recent branches shared between farmers and Pygmies. Coalescence analyses revealed that these two groups arose through a complex evolutionary process characterized by (i) initial divergence of the ancestors of contemporary Pygmies from an ancestral Central African population no more than ≈70,000 years ago, (ii) a period of isolation between the two groups, accounting for their phenotypic differences, (iii) long-standing asymmetric maternal gene flow from Pygmies to the ancestors of the farming populations, beginning no more than ≈40,000 years ago and persisting until a few thousand years ago, and (iv) enrichment of the maternal gene pool of the ancestors of the farming populations by the arrival and/or subsequent demographic expansion of L0a, L2, and L3 carriers.

Genetic and demographic implications of the Bantu expansion: insights from human paternal lineages

The expansion of Bantu languages, which started around 5,000 years before present in west/central Africa and spread all throughout sub-Saharan Africa, may represent one of the major and most rapid demographic movements in the history of the human species. Although the genetic footprints of this expansion have been unmasked through the analyses of the maternally inherited mitochondrial DNA lineages, information on the genetic impact of this massive movement and on the genetic composition of pre-Bantu populations is still scarce. Here, we analyze an extensive collection of Y-chromosome markers--41 single nucleotide polymorphisms and 18 short tandem repeats--in 883 individuals from 22 Bantu-speaking agriculturalist populations and 3 Pygmy hunter-gatherer populations from Gabon and Cameroon. Our data reveal a recent origin for most paternal lineages in west Central African populations most likely resulting from the expansion of Bantu-speaking farmers that erased the more ancient Y-chromosome diversity found in this area. However, some traces of ancient paternal lineages are observed in these populations, mainly among hunter-gatherers. These results are at odds with those obtained from mtDNA analyses, where high frequencies of ancient maternal lineages are observed, and substantial maternal gene flow from hunter-gatherers to Bantu farmers has been suggested. These differences are most likely explained by sociocultural factors such as patrilocality. We also find the intriguing presence of paternal lineages belonging to Eurasian haplogroup R1b1*, which might represent footprints of demographic expansions in central Africa not directly related to the Bantu expansion.

Phylogenetic and familial estimates of mitochondrial substitution rates : Study of control region mutations in deep-rooting pedigrees

We studied mutations in the mtDNA control region (CR) using deep-rooting French-Canadian pedigrees. In 508 maternal transmissions, we observed four substitutions (0.0079 per generation per 673 bp, 95% CI 0.0023-0.186). Combined with other familial studies, our results add up to 18 substitutions in 1,729 transmissions (0.0104), confirming earlier findings of much greater mutation rates in families than those based on phylogenetic comparisons. Only 12 of these mutations occurred at independent sites, whereas three positions mutated twice each, suggesting that pedigree studies preferentially reveal a fraction of highly mutable sites. Fitting the data through use of a nonuniform rate model predicts the presence of 40 (95% CI 27-54) such fast sites in the whole CR, characterized by the mutation rate of 274 per site per million generations (95% CI 138-410). The corresponding values for hypervariable regions I (HVI; 1,729 transmissions) and II (HVII; 1,956 transmissions), are 19 and 22 fast sites, with rates of 224 and 274, respectively. Because of the high probability of recurrent mutations, such sites are expected to be of no or little informativity for the evaluation of mutational distances at the phylogenetic time scale. The analysis of substitution density in the alignment of 973 HVI and 650 HVII unrelated European sequences reveals that the bulk of the sites mutate at relatively moderate and slow rates. Assuming a star-like phylogeny and an average time depth of 250 generations, we estimate the rates for HVI and HVII at 23 and 24 for the moderate sites and 1.3 and 1.0 for the slow sites. The fast, moderate, and slow sites, at the ratio of 1:2:13, respectively, describe the mutation-rate heterogeneity in the CR. Our results reconcile the controversial rate estimates in the phylogenetic and familial studies; the fast sites prevail in the latter, whereas the slow and moderate sites dominate the phylogenetic-rate estimations.

From Social to Genetic Structures in Central Asia

Pastoral and farmer populations, who have coexisted in Central Asia since the fourth millennium B.C., present not only different lifestyles and means of subsistence but also various types of social organization. Pastoral populations are organized into so-called descent groups (tribes, clans, and lineages) and practice exogamous marriages (a man chooses a bride in a different lineage or clan). In Central Asia, these descent groups are patrilineal: The children are systematically affiliated with the descent groups of the father. By contrast, farmer populations are organized into families (extended or nuclear) and often establish endogamous marriages with cousins. This study aims at better understanding the impact of these differences in lifestyle and social organization on the shaping of genetic diversity. We show that pastoral populations exhibit a substantial loss of Y chromosome diversity in comparison to farmers but that no such a difference is observed at the mitochondrial-DNA level. Our analyses indicate that the dynamics of patrilineal descent groups, which implies different male and female sociodemographic histories, is responsible for these sexually-asymmetric genetic patterns. This molecular signature of the pastoral social organization disappears over a few centuries only after conversion to an agricultural way of life.

Sex-Specific Genetic Structure and Social Organization in Central Asia: Insights from a Multi-Locus Study

In the last two decades, mitochondrial DNA (mtDNA) and the non-recombining portion of the Y chromosome (NRY) have been extensively used in order to measure the maternally and paternally inherited genetic structure of human populations, and to infer sex-specific demography and history. Most studies converge towards the notion that among populations, women are genetically less structured than men. This has been mainly explained by a higher migration rate of women, due to patrilocality, a tendency for men to stay in their birthplace while women move to their husbands house. Yet, since population differentiation depends upon the product of the effective number of individuals within each deme and the migration rate among demes, differences in male and female effective numbers and sex-biased dispersal have confounding effects on the comparison of genetic structure as measured by uniparentally inherited markers. In this study, we develop a new multi-locus approach to analyze jointly autosomal and X-linked markers in order to aid the understanding of sex-specific contributions to population differentiation. We show that in patrilineal herder groups of Central Asia, in contrast to bilineal agriculturalists, the effective number of women is higher than that of men. We interpret this result, which could not be obtained by the analysis of mtDNA and NRY alone, as the consequence of the social organization of patrilineal populations, in which genetically related men (but not women) tend to cluster together. This study suggests that differences in sex-specific migration rates may not be the only cause of contrasting male and female differentiation in humans, and that differences in effective numbers do matter.

Sex-specific demographic behaviours that shape human genomic variation

In the human species, the two uniparental genetic systems (mitochondrial DNA and Y chromosome) exhibit contrasting diversity patterns. It has been proposed that sex‐specific behaviours, and in particular differences in migration rate between men and women, may explain these differences. The availability of high‐density genomic data and the comparison of genetic patterns on autosomal and sex chromosomes at global and local scales allow a reassessment of the extent to which sex‐specific behaviours shape our genome. In this article, we first review studies comparing the genetic patterns at uniparental and biparental genetic systems and assess the extent to which sex‐specific migration processes explain the differences between these genetic systems. We show that differences between male and female migration rates matter, but that they are certainly not the only contributing factor. In particular, differences in effective population size between men and women are also likely to account for these differences. Then, we present and discuss three anthropological processes that may explain sex‐specific differences in effective population size and thus human genomic variation: (i) variance in reproductive success arising from, for example, polygyny; (ii) descent rules; and (iii) transmission of reproductive success.

Variation in Rural African Gut Microbiota Is Strongly Correlated with Colonization by Entamoeba and Subsistence.

The human gut microbiota is impacted by host nutrition and health status and therefore represents a potentially adaptive phenotype influenced by metabolic and immune constraints. Previous studies contrasting rural populations in developing countries to urban industrialized ones have shown that industrialization is strongly correlated with patterns in human gut microbiota; however, we know little about the relative contribution of factors such as climate, diet, medicine, hygiene practices, host genetics, and parasitism. Here, we focus on fine-scale comparisons of African rural populations in order to (i) contrast the gut microbiota of populations inhabiting similar environments but having different traditional subsistence modes and either shared or distinct genetic ancestry, and (ii) examine the relationship between gut parasites and bacterial communities. Characterizing the fecal microbiota of Pygmy hunter-gatherers as well as Bantu individuals from both farming and fishing populations in Southwest Cameroon, we found that the gut parasite Entamoeba is significantly correlated with microbiome composition and diversity. We show that across populations, colonization by this protozoa can be predicted with 79% accuracy based on the composition of an individuals gut microbiota, and that several of the taxa most important for distinguishing Entamoeba absence or presence are signature taxa for autoimmune disorders. We also found gut communities to vary significantly with subsistence mode, notably with some taxa previously shown to be enriched in other hunter-gatherers groups (in Tanzania and Peru) also discriminating hunter-gatherers from neighboring farming or fishing populations in Cameroon.

Population genetic diversity of the NAT2 gene supports a role of acetylation in human adaptation to farming in Central Asia.

The arylamine N-acetyltransferase 2 (NAT2) enzyme detoxifies a wide spectrum of naturally occurring xenobiotics including carcinogens and drugs. Variation at the NAT2 gene has been linked to the human acetylation capacity, either ‘slow’ or ‘fast’, which modifies susceptibility to cancer and adverse drug reactions. We investigated the possible influence of natural selection in shaping the acetylation phenotype and the NAT2 gene variability in six Central Asian populations, who are either long-term sedentary agriculturalists (two Tajik populations), recent sedentary agriculturalists (Kazakhs, Uzbeks) or nomad pastoralists (two Kirghiz populations). To this end, we sequenced the entire NAT2 coding exon, as well as genotyping nine intergenic SNPs covering a 200-kb region. Our results revealed that the two Tajik populations exhibited significantly higher proportions of slow acetylators than the nomadic populations. In addition, sequence-based neutrality tests yielded significantly positive values in Central Asian populations following an agriculturalist lifestyle, due to an excess of haplotypes at intermediate frequencies. Taken together, our data suggest that balancing selection, and/or directional selection on standing low-frequency alleles, have shaped NAT2 genetic diversity and the human acetylation phenotype in Central Asian agriculturalists. These results further support the hypothesis that a major transition in human lifestyle, such as the emergence of farming has dramatically changed human chemical environments and the selective pressures they imposed.