Wright-Fisher revisited: the case of fertility correlation

Detecting Recent Positive Selection with a Single Locus Test Bipartitioning the Coalescent Tree

Many population genomic studies have been conducted in the past to search for traces of recent events of positive selection. These traces, however, can be obscured by temporal variation of population size or other demographic factors. To reduce the confounding impact of demography, the coalescent tree topology has been used as an additional source of information for detecting recent positive selection in a population or a species. Based on the branching pattern at the root, we partition the hypothetical coalescent tree, inferred from a sequence sample, into two subtrees. The reasoning is that positive selection could impose a strong impact on branch length in one of the two subtrees while demography has the same effect on average on both subtrees. Thus, positive selection should be detectable by comparing statistics calculated for the two subtrees. Simulations demonstrate that the proposed test based on these principles has high power to detect recent positive selection even when DNA polymorphism data from only one locus is available, and that it is robust to the confounding effect of demography. One feature is that all components in the summary statistics ([Formula: see text]) can be computed analytically. Moreover, misinference of derived and ancestral alleles is seen to have only a limited effect on the test, and it therefore avoids a notorious problem when searching for traces of recent positive selection.

Relaxed Selection During a Recent Human Expansion

Humans have colonized the planet through a series of range expansions, which deeply impacted genetic diversity in newly settled areas and potentially increased the frequency of deleterious mutations on expanding wave fronts. To test this prediction, we studied the genomic diversity of French Canadians who colonized Quebec in the 17th century. We used historical information and records from ∼4000 ascending genealogies to select individuals whose ancestors lived mostly on the colonizing wave front and individuals whose ancestors remained in the core of the settlement. Comparison of exomic diversity reveals that: (i) both new and low-frequency variants are significantly more deleterious in front than in core individuals, (ii) equally deleterious mutations are at higher frequencies in front individuals, and (iii) front individuals are two times more likely to be homozygous for rare very deleterious mutations present in Europeans. These differences have emerged in the past six to nine generations and cannot be explained by differential inbreeding, but are consistent with relaxed selection mainly due to higher rates of genetic drift on the wave front. Demographic inference and modeling of the evolution of rare variants suggest lower effective size on the front, and lead to an estimation of selection coefficients that increase with conservation scores. Even though range expansions have had a relatively limited impact on the overall fitness of French Canadians, they could explain the higher prevalence of recessive genetic diseases in recently settled regions of Quebec.

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.

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.

Impact of fertility transmission and other sociodemographic factors on reproductive success and coalescent trees

Summary Fertility transmission (FT) is a phenomenon with a cultural and/or genetic basis, whereby a positive correlation exists between the number of offspring of an individual and that of his/her parents. Theoretical studies using a haploid individual-based model have shown that FT increases the variance and intergenerational correlation in reproductive success and results in an imbalance in the coalescent tree of sampled genes. This phenomenon has been documented in several demographic studies conducted on the correlation in fertility between generations, or through the reconstruction of the genealogical trees of mitochondrial DNA sequences. However, as mtDNA is a single locus, potentially subject to other forces (e.g. natural selection), it is of interest to extend the theory of FT to nuclear loci. We show that because random mating between individuals leads to a mixing of their fertility profiles, FT in these cases will have less influence on the variance and intergenerational correlation of reproductive success. This, in turn, results in less impact on the shape of the coalescent trees. Nevertheless, in the presence of FT, high heterogeneity in reproductive success and homogamy for family size will increase the imbalance in the coalescent tree. Thus, FT should be easier to detect when occurring in conjunction with these other factors. We also show the utility of analysing different kinds of loci (X-linked, Y-linked, mitochondrial and autosomal) to assess whether FT is matrilineal, patrilineal or biparental. Finally, we demonstrate that the shape of the coalescent tree depends upon population size, in contrast to the classical Kingman's model.

Deep human genealogies reveal a selective advantage to be on an expanding wave front

Since their origin, human populations have colonized the whole planet, but the demographic processes governing range expansions are mostly unknown. We analyzed the genealogy of more than one million individuals resulting from a range expansion in Quebec between 1686 and 1960 and reconstructed the spatial dynamics of the expansion. We find that a majority of the present Saguenay Lac-Saint-Jean population can be traced back to ancestors having lived directly on or close to the wave front. Ancestors located on the front contributed significantly more to the current gene pool than those from the range core, likely due to a 20% larger effective fertility of women on the wave front. This fitness component is heritable on the wave front and not in the core, implying that this life-history trait evolves during range expansions.

Genetic traces of east-to-west human expansion waves in Eurasia

In this study, we describe the landscape of human demographic expansions in Eurasia using a large continental Y chromosome and mitochondrial DNA dataset. Variation at these two uniparentally-inherited genetic systems retraces expansions that occurred in the past 60 ky, and shows a clear decrease of expansion ages from east to west Eurasia. To investigate the demographic events at the origin of this westward decrease of expansion ages, the estimated divergence ages between Eurasian populations are compared with the estimated expansion ages within each population. Both markers suggest that the demographic expansion diffused from east to west in Eurasia in a demic way, i.e., through migrations of individuals (and not just through diffusion of new technologies), highlighting the prominent role of eastern regions within Eurasia during Palaeolithic times.

Male dominance rarely skews the frequency distribution of Y chromosome haplotypes in human populations

A central tenet of evolutionary social science holds that behaviors, such as those associated with social dominance, produce fitness effects that are subject to cultural selection. However, evidence for such selection is inconclusive because it is based on short-term statistical associations between behavior and fertility. Here, we show that the evolutionary effects of dominance at the population level can be detected using noncoding regions of DNA. Highly variable polymorphisms on the nonrecombining portion of the Y chromosome can be used to trace lines of descent from a common male ancestor. Thus, it is possible to test for the persistence of differential fertility among patrilines. We examine haplotype distributions defined by 12 short tandem repeats in a sample of 1269 men from 41 Indonesian communities and test for departures from neutral mutation-drift equilibrium based on the Ewens sampling formula. Our tests reject the neutral model in only 5 communities. Analysis and simulations show that we have sufficient power to detect such departures under varying demographic conditions, including founder effects, bottlenecks, and migration, and at varying levels of social dominance. We conclude that patrilines seldom are dominant for more than a few generations, and thus traits or behaviors that are strictly paternally inherited are unlikely to be under strong cultural selection.

Transmission of migration propensity increases genetic divergence between populations

The advent of molecular genetics has brought invaluable information, which is now routinely used by anthropologists in their attempt to reconstruct our demographic past. Since mitochondrial DNA loci are much more similar between populations than are Y-chromosome loci, it is suggested that women had a much higher migration rate than men throughout history. Based on an examination of intergenerational migration patterns in three large demographic databases, we bring this inference into question. In some early Canadian settlements (St. Lawrence Valley and Saguenay), and in the former Krummhörn region of northwest Germany, men whose fathers were migrants were more likely to migrate, while the migration probability of women was largely independent of that of their mothers. As a result, men's movements were less effective in preventing genetic differentiation between populations than women's movements. In order to account for the impact of transmission, we propose a slight modification of Wright's island model. We also address the relevance of this model at the regional scale, and we discuss the supporting historical and anthropological literature. We conclude that the widespread patrilocal rules of postmarital residence have generated both a higher female migration rate and a patrilineal dependency in the propensity to migrate.

Matrilineal fertility inheritance detected in hunter-gatherer populations using the imbalance of gene genealogies

Fertility inheritance, a phenomenon in which an individual's number of offspring is positively correlated with his or her number of siblings, is a cultural process that can have a strong impact on genetic diversity. Until now, fertility inheritance has been detected primarily using genealogical databases. In this study, we develop a new method to infer fertility inheritance from genetic data in human populations. The method is based on the reconstruction of the gene genealogy of a sample of sequences from a given population and on the computation of the degree of imbalance in this genealogy. We show indeed that this level of imbalance increases with the level of fertility inheritance, and that other phenomena such as hidden population structure are unlikely to generate a signal of imbalance in the genealogy that would be confounded with fertility inheritance. By applying our method to mtDNA samples from 37 human populations, we show that matrilineal fertility inheritance is more frequent in hunter-gatherer populations than in food-producer populations. One possible explanation for this result is that in hunter-gatherer populations, individuals belonging to large kin networks may benefit from stronger social support and may be more likely to have a large number of offspring.

Serial coalescent simulations suggest a weak genealogical relationship between Etruscans and modern Tuscans

The Etruscans, the only preclassical European population that has been genetically characterized so far, share only two haplotypes with their modern geographic counterparts, the Tuscans, who, nonetheless, appear to be their closest relatives. We modeled 10 demographic scenarios spanning the last 2,500 years and tested by serial coalescent simulation whether any are consistent with the patterns of genetic diversity observed within and between the Etruscan and the modern Tuscan populations. Models in which the Etruscans are the direct ancestors of modern Tuscans appear compatible with the observed data only when they also include a very high mutation rate and an ancient founder effect. A better fit was obtained when the ancient and the modern samples were extracted from two independently evolving populations, connected by little migration. Simulated and observed parameters were also similar for a scenario in which the ancient samples came from a subset, e.g., a social elite, genetically differentiated from the bulk of the Etruscan population. In principle, these results may be biased by factors such as gross and systematic errors in the ancient DNA sequences and failure to sample suitable modern individuals. If neither proves to be the case, this study strongly suggests that either the mitochondrial mutation rate is much higher than currently believed or the Etruscans left very few modern mitochondrial descendants.

Estimating sex-specific processes in human populations: Are XY-homologous markers an effective tool?

Homologous markers on the sex-specific regions of the X- and Y-chromosomes are differentially inherited through males and females, and have similar molecular characteristics. They may therefore be useful as a complement to the comparison of mtDNA and Y-chromosomal haplotypes for estimating sex-specific processes shaping human population structure. To test this idea, we analyzed XY-homologous microsatellite diversity in 33 human populations from Africa, Asia and Europe. Interpopulation comparisons suggest that the generally discordant pattern of genetic variation observed for X- and Y-linked markers could be an outcome of sex-specific migration processes (m(females)/m(males) approximately 3) or sex-specific demographic processes (N(females)/N(males) approximately 11) or a combination of both. However, intrapopulation diversity estimated by the X/Y ratio Watterson estimator (theta(H(Y))/theta(H(X))) suggests that the scenarios required to explain the global genetic variation of XY-homologous markers are many and complex, and that the sex-specific processes (effective population size and migration rate) shaping human population structures are likely to be specific to each population under study. XY-homologous markers provide an insight into the genuine complexity of sex-specific processes, and their further exploitation in human population studies seems worthwhile.

Detecting Population Growth, Selection and Inherited Fertility From Haplotypic Data in Humans

The frequency of a rare mutant allele and the level of allelic association between this allele and one or several closely linked markers are frequently measured in genetic epidemiology. Both quantities are related to the time elapsed since the appearance of the mutation in the population and the intrinsic growth rate of the mutation (which may be different from the average population growth rate). Here, we develop a method that uses these two kinds of genetic data to perform a joint estimation of the age of the mutation and the minimum growth rate that is compatible with its present frequency. In absence of demographic data, it provides a useful estimate of population growth rate. When such data are available, contrasts among estimates from several loci allow demographic processes, affecting all loci similarly, to be distinguished from selection, affecting loci differently. Testing these estimates on populations for which data are available for several disorders shows good congruence with demographic data in some cases whereas in others higher growth rates are obtained, which may be the result of selection or hidden demographic processes.