Polygenic Adaptation: Integrating Population Genetics and Gene Regulatory Networks.

The adaptation of populations to local environments often relies on the selection of optimal values for polygenic traits. Here, we first summarize the results obtained from different quantitative genetics and population genetics models, about the genetic architecture of polygenic traits and their response to directional selection. We then highlight the contribution of systems biology to the understanding of the molecular bases of polygenic traits and the evolution of gene regulatory networks involved in these traits. Finally, we discuss the need for a unifying framework merging the fields of population genetics, quantitative genetics and systems biology to better understand the molecular bases of polygenic traits adaptation.

An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions.

Species generally undergo a complex demographic history consisting, in particular, of multiple changes in population size. Genome-wide sequencing data are potentially highly informative for reconstructing this demographic history. A crucial point is to extract the relevant information from these very large data sets. Here, we design an approach for inferring past demographic events from a moderate number of fully sequenced genomes. Our new approach uses Approximate Bayesian Computation, a simulation-based statistical framework that allows 1) identifying the best demographic scenario among several competing scenarios and 2) estimating the best-fitting parameters under the chosen scenario. Approximate Bayesian Computation relies on the computation of summary statistics. Using a cross-validation approach, we show that statistics such as the lengths of haplotypes shared between individuals, or the decay of linkage disequilibrium with distance, can be combined with classical statistics (e.g., heterozygosity and Tajima's D) to accurately infer complex demographic scenarios including bottlenecks and expansion periods. We also demonstrate the importance of simultaneously estimating the genotyping error rate. Applying our method on genome-wide human-sequence databases, we finally show that a model consisting in a bottleneck followed by a Paleolithic and a Neolithic expansion is the most relevant for Eurasian populations.

Sahelian pastoralism from the perspective of variants associated with lactase persistence.

OBJECTIVES: Archeological evidence shows that first nomadic pastoralists came to the African Sahel from northeastern Sahara, where milking is reported by ~7.5 ka. A second wave of pastoralists arrived with the expansion of Arabic tribes in 7th-14th century CE. All Sahelian pastoralists depend on milk production but genetic diversity underlying their lactase persistence (LP) is poorly understood. MATERIALS AND METHODS: We investigated SNP variants associated with LP in 1,241 individuals from 29 mostly pastoralist populations in the Sahel. Then, we analyzed six SNPs in the neighboring fragment (419 kb) in the Fulani and Tuareg with the -13910*T mutation, reconstructed haplotypes, and calculated expansion age and growth rate of this variant. RESULTS: Our results reveal a geographic localization of two different LP variants in the Sahel: -13910*T west of Lake Chad (Fulani and Tuareg pastoralists) and -13915*G east of there (mostly Arabic-speaking pastoralists). We show that -13910*T has a more diversified haplotype background among the Fulani than among the Tuareg and that the age estimate for expansion of this variant among the Fulani (~8.5 ka) corresponds to introduction of cattle to the area. CONCLUSIONS: This is the first study showing that the "Eurasian" LP allele -13910*T is widespread both in northern Europe and in the Sahel; however, it is limited to pastoralists in the Sahel. Since the Fulani haplotype with -13910*T is shared with contemporary Eurasians, its origin could be in a region encompassing the Near East and northeastern Africa in a population ancestral to both Saharan pastoralists and European farmers.

Neutral Theory: From Complex Population History to Natural Selection and Sociocultural Phenomena in Human Populations.

Here, we present a synthetic view on how Kimura's Neutral theory has helped us gaining insight on the different evolutionary forces that shape human evolution. We put this perspective in the frame of recent emerging challenges: the use of whole genome data for reconstructing population histories, natural selection on complex polygenic traits, and integrating cultural processes in human evolution.

An evaluation of the methods to estimate effective population size from measures of linkage disequilibrium.

In 1971, John Sved derived an approximate relationship between linkage disequilibrium (LD) and effective population size for an ideal finite population. This seminal work was extended by Sved and Feldman (Theor Pop Biol 4, 129, 1973) and Weir and Hill (Genetics 95, 477, 1980) who derived additional equations with the same purpose. These equations yield useful estimates of effective population size, as they require a single sample in time. As these estimates of effective population size are now commonly used on a variety of genomic data, from arrays of single nucleotide polymorphisms to whole genome data, some authors have investigated their bias through simulation studies and proposed corrections for different mating systems. However, the cause of the bias remains elusive. Here, we show the problems of using LD as a statistical measure and, analogously, the problems in estimating effective population size from such measure. For that purpose, we compare three commonly used approaches with a transition probability-based method that we develop here. It provides an exact computation of LD. We show here that the bias in the estimates of LD and effective population size are partly due to low-frequency markers, tightly linked markers or to a small total number of crossovers per generation. These biases, however, do not decrease when increasing sample size or using unlinked markers. Our results show the issues of such measures of effective population based on LD and suggest which of the method here studied should be used in empirical studies as well as the optimal distance between markers for such estimates.

Inferring Population Size History from Large Samples of Genome-Wide Molecular Data - An Approximate Bayesian Computation Approach.

Inferring the ancestral dynamics of effective population size is a long-standing question in population genetics, which can now be tackled much more accurately thanks to the massive genomic data available in many species. Several promising methods that take advantage of whole-genome sequences have been recently developed in this context. However, they can only be applied to rather small samples, which limits their ability to estimate recent population size history. Besides, they can be very sensitive to sequencing or phasing errors. Here we introduce a new approximate Bayesian computation approach named PopSizeABC that allows estimating the evolution of the effective population size through time, using a large sample of complete genomes. This sample is summarized using the folded allele frequency spectrum and the average zygotic linkage disequilibrium at different bins of physical distance, two classes of statistics that are widely used in population genetics and can be easily computed from unphased and unpolarized SNP data. Our approach provides accurate estimations of past population sizes, from the very first generations before present back to the expected time to the most recent common ancestor of the sample, as shown by simulations under a wide range of demographic scenarios. When applied to samples of 15 or 25 complete genomes in four cattle breeds (Angus, Fleckvieh, Holstein and Jersey), PopSizeABC revealed a series of population declines, related to historical events such as domestication or modern breed creation. We further highlight that our approach is robust to sequencing errors, provided summary statistics are computed from SNPs with common alleles.

Using genome-wide measures of coancestry to maintain diversity and fitness in endangered and domestic pig populations.

Conservation and breeding programs aim at maintaining the most diversity, thereby avoiding deleterious effects of inbreeding while maintaining enough variation from which traits of interest can be selected. Theoretically, the most diversity is maintained using optimal contributions based on many markers to calculate coancestries, but this can decrease fitness by maintaining linked deleterious variants. The heterogeneous patterns of coancestry displayed in pigs make them an excellent model to test these predictions. We propose methods to measure coancestry and fitness from resequencing data and use them in population management. We analyzed the resequencing data of Sus cebifrons, a highly endangered porcine species from the Philippines, and genotype data from the Pietrain domestic breed. By analyzing the demographic history of Sus cebifrons, we inferred two past bottlenecks that resulted in some inbreeding load. In Pietrain, we analyzed signatures of selection possibly associated with commercial traits. We also simulated the management of each population to assess the performance of different optimal contribution methods to maintain diversity, fitness, and selection signatures. Maximum genetic diversity was maintained using marker-by-marker coancestry, and least using genealogical coancestry. Using a measure of coancestry based on shared segments of the genome achieved the best results in terms of diversity and fitness. However, this segment-based management eliminated signatures of selection. We demonstrate that maintaining both diversity and fitness depends on the genomic distribution of deleterious variants, which is shaped by demographic and selection histories. Our findings show the importance of genomic and next-generation sequencing information in the optimal design of breeding or conservation programs.

Genetic and linguistic histories in Central Asia inferred using approximate Bayesian computations.

Linguistic and genetic data have been widely compared, but the histories underlying these descriptions are rarely jointly inferred. We developed a unique methodological framework for analysing jointly language diversity and genetic polymorphism data, to infer the past history of separation, exchange and admixture events among human populations. This method relies on approximate Bayesian computations that enable the identification of the most probable historical scenario underlying each type of data, and to infer the parameters of these scenarios. For this purpose, we developed a new computer program PopLingSim that simulates the evolution of linguistic diversity, which we coupled with an existing coalescent-based genetic simulation program, to simulate both linguistic and genetic data within a set of populations. Applying this new program to a wide linguistic and genetic dataset of Central Asia, we found several differences between linguistic and genetic histories. In particular, we showed how genetic and linguistic exchanges differed in the past in this area: some cultural exchanges were maintained without genetic exchanges. The methodological framework and the linguistic simulation tool developed here can be used in future work for disentangling complex linguistic and genetic evolutions underlying human biological and cultural histories.

Trait level analysis of multitrait population projection matrices.

In most matrix population projection models, individuals are characterized according to, usually, one or two traits such as age, stage, size or location. A broad theory of multitrait population projection matrices (MPPMs) incorporating larger number of traits was long held back by time and space computational complexity issues. As a consequence, no study has yet focused on the influence of the structure of traits describing a life-cycle on population dynamics and life-history evolution. We present here a novel vector-based MPPM building methodology that allows to computationally-efficiently model populations characterized by numerous traits with large distributions, and extend sensitivity analyses for these models. We then present a new method, the trait level analysis consisting in folding an MPPM on any of its traits to create a matrix with alternative trait structure (the number of traits and their characteristics) but similar asymptotic properties. Adding or removing one or several traits to/from the MPPM and analyzing the resulting changes in spectral properties, allows investigating the influence of the trait structure on the evolution of traits. We illustrate this by modeling a 3-trait (age, parity and fecundity) population designed to investigate the implications of parity-fertilitytrade-offs in a context of fecundity heterogeneity in humans. The trait level analysis, comparing models of the same population differing in trait structures, demonstrates that fertility selection gradients differ between cases with or without parity-fertility trade-offs. Moreover it shows that age-specific fertility has seemingly very different evolutionary significance depending on whether heterogeneity is accounted for. This is because trade-offs can vary strongly in strength and even direction depending on the trait structure used to model the population.

The historical spread of Arabian Pastoralists to the eastern African Sahel evidenced by the lactase persistence -13,915*G allele and mitochondrial DNA.

OBJECTIVES: Thanks to the ability to digest lactose, Arabian nomads had become less dependent upon their sedentary neighbors and some of these populations spread to Africa. When and by which route they migrated to their current locations have previously been addressed only by historical and archaeological data. METHODS: To address the question of Arab expansion into Africa, we collected samples from several Arabic populations, especially the Baggara in Chad and Sudan. We analyzed mutations associated with lactase persistence and reconstructed the surrounding haplotypes defined by SNP polymorphisms. We also sequenced their mitochondrial DNA to investigate relative proportions of sub-Saharan and Eurasian origins. RESULTS: We estimated the expansion age of the -13,915*G mutation in four different Arabian datasets. The oldest age was identified in Yemen (1,356-1,799 ya) and the youngest in a Sudanese group of Rashaayda Arabs (219-312 ya). We also found a negative correlation between the frequency of the -13,915*G allele and the frequency of sub-Saharan mtDNA haplotypes. CONCLUSIONS: Even if the age of the most recent common ancestor of -13,915*G is ∼4 ka as shown in a previous study, our results suggest that its spread to Africa was more recent, which is consistent with the migrations of Arabic tribes. Because the incidence of sub-Saharan mtDNA haplotypes is negatively correlated with the occurrence of -13,915*G, we suggest that the decrease of its frequency in Africa has been caused by progressive admixture of the Arabian nomads with sub-Saharan populations.

Statistical inference on genetic data reveals the complex demographic history of human populations in central Asia.

The demographic history of modern humans constitutes a combination of expansions, colonizations, contractions, and remigrations. The advent of large scale genetic data combined with statistically refined methods facilitates inference of this complex history. Here we study the demographic history of two genetically admixed ethnic groups in Central Asia, an area characterized by high levels of genetic diversity and a history of recurrent immigration. Using Approximate Bayesian Computation, we infer that the timing of admixture markedly differs between the two groups. Admixture in the traditionally agricultural Tajiks could be dated back to the onset of the Neolithic transition in the region, whereas admixture in Kyrgyz is more recent, and may have involved the westward movement of Turkic peoples. These results are confirmed by a coalescent method that fits an isolation-with-migration model to the genetic data, with both Central Asian groups having received gene flow from the extremities of Eurasia. Interestingly, our analyses also uncover signatures of gene flow from Eastern to Western Eurasia during Paleolithic times. In conclusion, the high genetic diversity currently observed in these two Central Asian peoples most likely reflects the effects of recurrent immigration that likely started before historical times. Conversely, conquests during historical times may have had a relatively limited genetic impact. These results emphasize the need for a better understanding of the genetic consequences of transmission of culture and technological innovations, as well as those of invasions and conquests.

Inference of sex-specific expansion patterns in human populations from Y-chromosome polymorphism.

Studying the current distribution of genetic diversity in humans has important implications for our understanding of the history of our species. We analyzed a set of linked STR and SNP loci from the paternally inherited Y chromosome to infer the past demography of 55 African and Eurasian populations, using both the parametric and nonparametric coalescent-based methods implemented in the BEAST application. We inferred expansion events in most sedentary farmer populations, while we found constant effective population sizes for both nomadic hunter-gatherers and seminomadic herders. Our results differed, on several aspects, from previous results on mtDNA and autosomal markers. First, we found more recent expansion patterns in Eurasia than in Africa. This discrepancy, substantially stronger than the ones found with the other kind of markers, may result from a lower effective population size for men, which might have made male-transmitted markers more sensitive to the out-of-Africa bottleneck. Second, we found expansion signals only for sedentary farmers but not for nomadic herders in Central Asia, while these signals were found for both kind of populations in this area when using mtDNA or autosomal markers. Expansion signals in this area may result from spatial expansion processes and may have been erased for the Y chromosome among the herders because of restricted male gene flow.

Patrilineal populations show more male transmission of reproductive success than cognatic populations in Central Asia, which reduces their genetic diversity.

OBJECTIVE: The extent to which social organization of human societies impacts the patterns of genetic diversity remains an open question. Here, we investigate the transmission of reproductive success in patrilineal and cognatic populations from Central Asia using a coalescent approach. METHODS: We performed a study on the mitochondrial DNA (mtDNA) and Y chromosome polymorphism of patrilineal and cognatic populations from Central Asia. We reconstructed the gene genealogies in each population for both kind of markers and inferred the imbalance level of these genealogies, a parameter directly related to the level of transmission of reproductive success. RESULTS: This imbalance level appeared much stronger for the Y chromosome in patrilineal populations than in cognatic populations, while no difference was found for mtDNA. Furthermore, we showed that this imbalance level correlates negatively with Y-chromosomal, mtDNA, and autosomal genetic diversity. CONCLUSIONS: This shows that patrilineality might be one of the factors explaining the male transmission of reproductive success, which, in turn, lead to a reduction of genetic diversity. Thus, notwithstanding the fact that our population genetic approach clearly shows that there is a strong male-biased transmission of reproductive success in patrilineal societies, it also highlights the fact that a social process such as cultural transmission of reproductive success could play an important role in shaping human genetic diversity, although we cannot formally exclude that this transmission has also a genetic component.

History of expansion and anthropogenic collapse in a top marine predator of the Black Sea estimated from genetic data.

Two major ecological transitions marked the history of the Black Sea after the last Ice Age. The first was the postglacial transition from a brackish-water to a marine ecosystem dominated by porpoises and dolphins once this basin was reconnected back to the Mediterranean Sea (ca. 8,000 y B.P.). The second occurred during the past decades, when overfishing and hunting activities brought these predators close to extinction, having a deep impact on the structure and dynamics of the ecosystem. Estimating the extent of this decimation is essential for characterizing this ecosystem's dynamics and for formulating restoration plans. However, this extent is poorly documented in historical records. We addressed this issue for one of the main Black Sea predators, the harbor porpoise, using a population genetics approach. Analyzing its genetic diversity using an approximate Bayesian computation approach, we show that only a demographic expansion (at most 5,000 y ago) followed by a contemporaneous population collapse can explain the observed genetic data. We demonstrate that both the postglacial settlement of harbor porpoises in the Black Sea and the recent anthropogenic activities have left a clear footprint on their genetic diversity. Specifically, we infer a strong population reduction (~90%) that occurred within the past 5 decades, which can therefore clearly be related to the recent massive killing of small cetaceans and to the continuing incidental catches in commercial fisheries. Our study thus provides a quantitative assessment of these demographically catastrophic events, also showing that two separate historical events can be inferred from contemporary genetic data.

Human genetic data reveal contrasting demographic patterns between sedentary and nomadic populations that predate the emergence of farming.

Demographic changes are known to leave footprints on genetic polymorphism. Together with the increased availability of large polymorphism data sets, coalescent-based methods allow inferring the past demography of populations from their present-day patterns of genetic diversity. Here, we analyzed both nuclear (20 noncoding regions) and mitochondrial (HVS-I) resequencing data to infer the demographic history of 66 African and Eurasian human populations presenting contrasting lifestyles (nomadic hunter-gatherers, nomadic herders, and sedentary farmers). This allowed us to investigate the relationship between lifestyle and demography and to address the long-standing debate about the chronology of demographic expansions and the Neolithic transition. In Africa, we inferred expansion events for farmers, but constant population sizes or contraction events for hunter-gatherers. In Eurasia, we inferred higher expansion rates for farmers than herders with HVS-I data, except in Central Asia and Korea. Although isolation and admixture processes could have impacted our demographic inferences, these processes alone seem unlikely to explain the contrasted demographic histories inferred in populations with different lifestyles. The small expansion rates or constant population sizes inferred for herders and hunter-gatherers may thus result from constraints linked to nomadism. However, autosomal data revealed contraction events for two sedentary populations in Eurasia, which may be caused by founder effects. Finally, the inferred expansions likely predated the emergence of agriculture and herding. This suggests that human populations could have started to expand in Paleolithic times, and that strong Paleolithic expansions in some populations may have ultimately favored their shift toward agriculture during the Neolithic.

Sociocultural behavior, sex-biased admixture, and effective population sizes in Central African Pygmies and non-Pygmies.

Sociocultural phenomena, such as exogamy or phylopatry, can largely determine human sex-specific demography. In Central Africa, diverging patterns of sex-specific genetic variation have been observed between mobile hunter-gatherer Pygmies and sedentary agricultural non-Pygmies. However, their sex-specific demography remains largely unknown. Using population genetics and approximate Bayesian computation approaches, we inferred male and female effective population sizes, sex-specific migration, and admixture rates in 23 Central African Pygmy and non-Pygmy populations, genotyped for autosomal, X-linked, Y-linked, and mitochondrial markers. We found much larger effective population sizes and migration rates among non-Pygmy populations than among Pygmies, in agreement with the recent expansions and migrations of non-Pygmies and, conversely, the isolation and stationary demography of Pygmy groups. We found larger effective sizes and migration rates for males than for females for Pygmies, and vice versa for non-Pygmies. Thus, although most Pygmy populations have patrilocal customs, their sex-specific genetic patterns resemble those of matrilocal populations. In fact, our results are consistent with a lower prevalence of polygyny and patrilocality in Pygmies compared with non-Pygmies and a potential female transmission of reproductive success in Pygmies. Finally, Pygmy populations showed variable admixture levels with the non-Pygmies, with often much larger introgression from male than from female lineages. Social discrimination against Pygmies triggering complex movements of spouses in intermarriages can explain these male-biased admixture patterns in a patrilocal context. We show how gender-related sociocultural phenomena can determine highly variable sex-specific demography among populations, and how population genetic approaches contrasting chromosomal types allow inferring detailed human sex-specific demographic history.

Postglacial climate changes and rise of three ecotypes of harbour porpoises, Phocoena phocoena, in western Palearctic waters.

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base-pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23-19 kilo-years ago, kyrbp). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp) than either does with the European continental shelf ecotype (c. 28 kyrbp), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post-glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.

Genetic signature of a range expansion and leap-frog event after the recent invasion of Europe by the grapevine downy mildew pathogen Plasmopara viticola.

Biologic invasions can have important ecological, economic and social consequences, particularly when they involve the introduction and spread of plant invasive pathogens, as they can threaten natural ecosystems and jeopardize the production of human food. Examples include the grapevine downy mildew, caused by the oomycete Plasmopara viticola, an invasive species native to North America, introduced into Europe in the 1870s. We investigated the introduction and spread of this invasive pathogen, by analysing its genetic structure and diversity in a large sample from European vineyards. Populations of P. viticola across Europe displayed little genetic diversity, consistent with the occurrence of a bottleneck at the time of introduction. Bayesian coalescent analyses revealed a clear population expansion signal in the genetic data. We detected a weak, but significant, continental-wide population structure, with two geographically and genetically distinct clusters in Western and Eastern European vineyards. Approximate Bayesian computation, analyses of clines of genetic diversity and of isolation-by-distance patterns provided evidence for a wave of colonization moving in an easterly direction across Europe. This is consistent with historical reports, first mentioning the introduction of the disease in Bordeaux vineyards (France) and sub-sequently documenting its rapid spread across Europe. This initial introduction in the west was probably followed by a 'leap-frog' event into Eastern Europe, leading to the formation of the two genetic clusters we detected. This study shows that recent population genetics methods within the Bayesian and coalescence frameworks are extremely powerful for increasing our understanding of pathogen population dynamics and invasion histories.

Consecutive five-year analysis of paternal and maternal gene flow and contributions of gametic heterogeneities to overall genetic composition of dispersed seeds of Pinus densiflora (Pinaceae).

PREMISE OF THE STUDY: Genetic variability in monoecious woody plant populations results from the assemblage of individuals issued from asymmetrical male and female reproductive functions, produced during spatially and temporarily heterogeneous reproductive and dispersal events. Here we investigated the dispersal patterns and levels of genetic diversity and differentiation of both paternal and maternal gametes in a natural population of Pinus densiflora at the multiple-year scale as long as five consecutive years. • METHODS: We analyzed the paternity and maternity for 1576 seeds and 454 candidate adult trees using nuclear DNA polymorphisms of diploid biparental embryos and haploid maternal megagametophytes at eight microsatellite loci. • KEY RESULTS: Despite the low levels of genetic differentiation among gamete groups, a two-way AMOVA analysis showed that the parental origin (paternal vs. maternal gametes), the year of gamete production and their interaction had significant effects on the genetic composition of the seeds. While maternal gamete groups showed a significant FST value across the 5 years, this was not true for their paternal counterparts. Within the population, we found that the relative reproductive contributions of the paternal vs. the maternal parent differed among adult trees, the maternal contributions showing a larger year-to-year fluctuation. • CONCLUSIONS: The overall genetic variability of dispersed seeds appeared to result from two sources of heterogeneity: the difference between paternal and maternal patterns of reproduction and gamete dispersal and year-to-year heterogeneity of reproduction of adult trees, especially in their maternal reproduction.