Population genetics of autopolyploids under a mixed mating model and the estimation of selfing rate

Nowadays, the population genetics analysis of autopolyploid species faces many difficulties due to (i) limited development of population genetics tools under polysomic inheritance, (ii) difficulties to assess allelic dosage when genotyping individuals and (iii) a form of inbreeding resulting from the mechanism of ‘double reduction’. Consequently, few data analysis computer programs are applicable to autopolyploids. To contribute bridging this gap, this article first derives theoretical expectations for the inbreeding and identity disequilibrium coefficients under polysomic inheritance in a mixed mating model. Moment estimators of these coefficients are proposed when exact genotypes or just markers phenotypes (i.e. allelic dosage unknown) are available. This led to the development of estimators of the selfing rate based on adult genotypes or phenotypes and applicable to any even‐ploidy level. Their statistical performances and robustness were assessed by numerical simulations. Contrary to inbreeding‐based estimators, the identity disequilibrium‐based estimator using phenotypes is robust (absolute bias generally < 0.05), even in the presence of double reduction, null alleles or biparental inbreeding due to isolation by distance. A fairly good precision of the selfing rate estimates (root mean squared error < 0.1) is already achievable using a sample of 30–50 individuals phenotyped at 10 loci bearing 5–10 alleles each, conditions reachable using microsatellite markers. Diallelic markers (e.g. SNP) can also perform satisfactorily in diploids and tetraploids but more polymorphic markers are necessary for higher ploidy levels. The method is implemented in the software SPAGeDi and should contribute to reduce the lack of population genetics tools applicable to autopolyploids.     

Effect of partial selfing and polygenic selection on establishment in a new habitat

This article analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the inbreeding history model to characterize mutation‐selection balance in a large, partially selfing source population under selection involving multiple nonidentical loci. I then use individual‐based simulations to study the eco‐evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long‐term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection is discussed.     

HIGH INBREEDING DEPRESSION,SELECTIVE INTERFERENCE AMONG LOCI,AND THE THRESHOLD SELFING RATE FOR PURGING RECESSIVE LETHAL MUTATIONS

The evolutionary dynamics of recessive or slightly dominant lethal mutations in partially self-fertilizing plants are analyzed using two models. In the identity-equilibrium model, lethals occur at a finite number of unlinked loci among which genotype frequencies are independent in mature plants. In the Kondrashov model, lethals occur at an infinite number of unlinked loci with identity disequilibrium produced by partial selfing. If the genomic mutation rate to (nearly) recessive lethal alleles is sufficiently high, such that the mean number of lethals (or lethal equivalents) per mature plant maintained at equilibrium under complete outcrossing exceeds 10, selective interference among loci creates a sharp discontinuity in the mean number of lethals maintained as a function of the selfing rate. Virtually no purging of the lethals occurs unless the selfing rate closely approaches or exceeds a threshold selfing rate, at which there is a precipitous drop in the mean number of lethals maintained. Identity disequilibrium lowers the threshold selfing rate by increasing the ratio of variance to mean number of lethals per plant, increasing the opportunity for selection. This theory helps to explain observations on plant species that display very high inbreeding depression despite intermediate selfing rates.     

INFERENCES ABOUT INBREEDING DEPRESSION BASED ON CHANGES OF THE INBREEDING COEFFICIENT

In a mixed-mating population, the fitness of selfed individuals, relative to outcrossed individuals, can be estimated from observed changes of the inbreeding coefficient F. In general, three measurements of F, or two measurements of F and one of the selfing rate, spanning two generations are needed. If however, adult F is assumed constant among generations, the adult F and selfing rate of one generation are sufficient to estimate relative fitness. Estimates of the relative fitness of selfed individuals for 14 Mimulus populations, assuming equilibrium of adult F, averaged 0.38, which is significantly lower than the 0.50 threshold needed to favor selfing genotypes. However, estimates for some populations showed large variance and over the 14 populations, relative fitness did not correlate with selfing rate. Violations of the equilibrium assumption causes a positive bias of the estimate and a spurious negative correlation of estimates with selfing rate. Estimates of relative fitness based on two generations of F do not suffer these problems. The need for large sample sizes, multiallelic loci, and moderate levels of natural selfing limits the usefulness of using changes of F to infer inbreeding depression.     

Quantifying inbreeding in natural populations of hermaphroditic organisms

We review molecular methods for estimating selfing rates and inbreeding in populations. Two main approaches are available: the population structure approach (PSA) and progeny-array approach (PAA). The PSA approach relies on single-generation samples and produces estimates that integrate the inbreeding history over several generations, but is based on strong assumptions (for example, inbreeding equilibrium). The PSA has classically relied on single-locus inbreeding coefficients averaged over loci. Unfortunately PSA estimates are very sensitive to technical problems such as the occurrence of null alleles at one or more of the loci. Consequently inbreeding might be substantially overestimated, especially in outbred populations. However, the robustness of the PSA has recently been greatly improved by the development of multilocus methods free of such bias. The PAA, on the other hand, is based on the comparison between offspring and mother genotypes. As a consequence, PAA estimates do not reflect long-term inbreeding history but only recent mating events of the maternal individuals studied ('here and now' selfing). In addition to selfing rates, the PAA allows estimating other mating system parameters, including biparental inbreeding and the correlation of selfing among sibs. Although PAA estimates could also be biased by technical problems, incompatibilities between the mother's genotype and her offspring allow the identification and correction of such bias. For all methods, we provide guidelines on the required number of loci and sample sizes. We conclude that the PSA and PAA are equally robust, provided multilocus information is used. Although experimental constraints may make the PAA more demanding, especially in animals, the two methods provide complementary information, and can fruitfully be conducted together.     

Estimation of quantitative genetic parameters using marker-inferred relatedness in Japanese flounder: a case study of upward bias

Marker-based methods for estimating heritability have been proposed as an effective means to study quantitative traits in long-lived organisms and natural populations. However, practical examinations to evaluate the usefulness and robustness of a regression method are limited. Using several quantitative traits of Japanese flounder Paralichthys olivaceus, the present study examined the influence of relatedness estimator and population structure on the estimation of heritability and genetic correlation under a regression method with 7 microsatellite loci. Significant heritability and genetic correlation were detected for several quantitative traits in 2 laboratory populations but not in a natural population. In the laboratory populations, upward bias in heritability appeared depending on the relatedness estimators and the populations. Upward bias in heritability increased with decreasing the actual variance of relatedness, suggesting that the estimates of heritability under the regression method tend to be overestimated due to the underestimation of the actual variance of relatedness. Therefore, relationship structure and precise estimation of relatedness are critical for applying this method.     

Estimating selfing rates from reconstructed pedigrees using multilocus genotype data

Several methods have been developed to estimate the selfing rate of a population from a sample of individuals genotyped for several marker loci. These methods can be based on homozygosity excess (or inbreeding), identity disequilibrium, progeny array (PA) segregation or population assignment incorporating partial selfing. Progeny array-based method is generally the best because it is not subject to some assumptions made by other methods (such as lack of misgenotyping, absence of biparental inbreeding and presence of inbreeding equilibrium), and it can reveal other facets of a mixed-mating system such as patterns of shared paternity. However, in practice, it is often difficult to obtain PAs, especially for animal species. In this study, we propose a method to reconstruct the pedigree of a sample of individuals taken from a monoecious diploid population practicing mixed mating, using multilocus genotypic data. Selfing and outcrossing events are then detected when an individual derives from identical parents and from two distinct parents, respectively. Selfing rate is estimated by the proportion of selfed offspring in the reconstructed pedigree of a sample of individuals. The method enjoys many advantages of the PA method, but without the need of a priori family structure, although such information, if available, can be utilized to improve the inference. Furthermore, the new method accommodates genotyping errors, estimates allele frequencies jointly and is robust to the presence of biparental inbreeding and inbreeding disequilibrium. Both simulated and empirical data were analysed by the new and previous methods to compare their statistical properties and accuracies.     

Evolutionary implications of a high selfing rate in the freshwater snail Lymnaea truncatula

Abstract Self-compatible hermaphroditic organisms that mix self-fertilization and outcrossing are of great interest for investigating the evolution of mating systems. We investigate the evolution of selfing in Lymnaea truncatula , a self-compatible hermaphroditic freshwater snail. We first analyze the consequences of selfing in terms of genetic variability within and among populations and then investigate how these consequences along with the species ecology (harshness of the habitat and parasitism) might govern the evolution of selfing. Snails from 13 localities (classified as temporary or permanent depending on their water availability) were sampled in western Switzerland and genotyped for seven microsatellite loci. F_IS (estimated on adults) and progeny array analyses (on hatchlings) provided similar selfing rate estimates of 80%. Populations presented a low polymorphism and were highly differentiated (F_ST= 0.58). Although the reproductive assurance hypothesis would predict higher selfing rate in temporary populations, no difference in selfing level was observed between temporary and permanent populations. However, allelic richness and gene diversity declined in temporary habitats, presumably reflecting drift. Infection levels varied but were not simply related to either estimated population selfing rate or to differences in heterozygosity. These findings and the similar selfing rates estimated for hatchlings and adults suggest that within-population inbreeding depression is low in L. truncatula.     

A Markov chain Monte Carlo approach for joint inference of population structure and inbreeding rates from multilocus genotype data

Nonrandom mating induces correlations in allelic states within and among loci that can be exploited to understand the genetic structure of natural populations (Wright 1965). For many species, it is of considerable interest to quantify the contribution of two forms of nonrandom mating to patterns of standing genetic variation: inbreeding (mating among relatives) and population substructure (limited dispersal of gametes). Here, we extend the popular Bayesian clustering approach STRUCTURE (Pritchard et al. 2000) for simultaneous inference of inbreeding or selfing rates and population-of-origin classification using multilocus genetic markers. This is accomplished by eliminating the assumption of Hardy-Weinberg equilibrium within clusters and, instead, calculating expected genotype frequencies on the basis of inbreeding or selfing rates. We demonstrate the need for such an extension by showing that selfing leads to spurious signals of population substructure using the standard STRUCTURE algorithm with a bias toward spurious signals of admixture. We gauge the performance of our method using extensive coalescent simulations and demonstrate that our approach can correct for this bias. We also apply our approach to understanding the population structure of the wild relative of domesticated rice, Oryza rufipogon, an important partially selfing grass species. Using a sample of n = 16 individuals sequenced at 111 random loci, we find strong evidence for existence of two subpopulations, which correlates well with geographic location of sampling, and estimate selfing rates for both groups that are consistent with estimates from experimental data (s approximately 0.48-0.70).     

RAD‐sequencing for estimating genomic relatedness matrix‐based heritability in the wild: A case study in roe deer

Estimating the evolutionary potential of quantitative traits and reliably predicting responses to selection in wild populations are important challenges in evolutionary biology. The genomic revolution has opened up opportunities for measuring relatedness among individuals with precision, enabling pedigree‐free estimation of trait heritabilities in wild populations. However, until now, most quantitative genetic studies based on a genomic relatedness matrix (GRM) have focused on long‐term monitored populations for which traditional pedigrees were also available, and have often had access to knowledge of genome sequence and variability. Here, we investigated the potential of RAD‐sequencing for estimating heritability in a free‐ranging roe deer (Capreolous capreolus) population for which no prior genomic resources were available. We propose a step‐by‐step analytical framework to optimize the quality and quantity of the genomic data and explore the impact of the single nucleotide polymorphism (SNP) calling and filtering processes on the GRM structure and GRM‐based heritability estimates. As expected, our results show that sequence coverage strongly affects the number of recovered loci, the genotyping error rate and the amount of missing data. Ultimately, this had little effect on heritability estimates and their standard errors, provided that the GRM was built from a minimum number of loci (above 7,000). Genomic relatedness matrix‐based heritability estimates thus appear robust to a moderate level of genotyping errors in the SNP data set. We also showed that quality filters, such as the removal of low‐frequency variants, affect the relatedness structure of the GRM, generating lower h² estimates. Our work illustrates the huge potential of RAD‐sequencing for estimating GRM‐based heritability in virtually any natural population.     

Locus-specific heritability estimation via the bootstrap in linkage scans for quantitative trait loci

BACKGROUND/AIMS: In genome-wide linkage analysis of quantitative trait loci (QTL), locus-specific heritability estimates are biased when the original data are used to both localize linkage and estimate effects, due to maximization of the LOD score over the genome. Positive bias is increased by adoption of stringent significance levels to control genome-wide type I error. We propose multi-locus bootstrap resampling estimators for bias reduction in the situation in which linkage peaks at more than one QTL are of interest. METHODS: Bootstrap estimates were based on repeated sample splitting in the original dataset. We conducted simulation studies in nuclear families with 0 to 5 QTLs and applied the methods in a genome-wide analysis of a blood pressure phenotype in extended pedigrees from the Framingham Heart Study (FHS). RESULTS: Compared to na?ve estimates in the original simulation samples, bootstrap estimates had reduced bias and smaller mean squared error. In the FHS pedigrees, the bootstrap yielded heritability estimates as much as 70% smaller than in the original sample. CONCLUSIONS: Because effect estimates obtained in an initial study are typically inflated relative to those expected in an independent replication study, successful replication will be more likely when sample size requirements are based on bias-reduced estimates.     

INBREEDING DEPRESSION IN TWO MIMULUS TAXA MEASURED BY MULTIGENERATIONAL CHANGES IN THE INBREEDING COEFFICIENT

In mixed-mating plant populations, one can estimate the relative fitness of selfed progeny w by measuring the inbreeding coefficient F and selfing rate s of adults of one generation, together with F of adults in the following generation (after selection). In the first application of this multigenerational method, we estimated F and s for adults over three consecutive generations in adjacent populations of two annual Mimulus taxa: the outbreeding M. guttatus and the inbreeding M. platycalyx. This gave estimates of w for the last two generations. Although average multilocus selfing rates were high in both taxa (0.63 in M. guttatus; 0.84 in M. platycalyx), the relative fitness of selfed progeny averaged only 0.19 in M. guttatus and 0.32 in M. platycalyx. An alternative estimator for w that incorporates biparental inbreeding gave even lower estimates of w. These values are significantly below the 0.5 threshold thought to favor selfing, and show that partially selfing populations can harbor substantial genetic load. In accordance with the purging hypothesis, the more highly selfing M. platycalyx showed marginally lower inbreeding depression than M. guttatus in both years (P = 0.08). Inbreeding depression and selfing rates also varied among years in concert among taxa. Several sources of bias are discussed, but computer simulations indicate it is unlikely that w is biased downwards by linkage of marker loci to load loci.     

Estimation of genealogical coancestry in plant species using a pedigree reconstruction algorithm and application to an oil palm breeding population

Key message

Explicit pedigree reconstruction by simulated annealing gave reliable estimates of genealogical coancestry in plant species, especially when selfing rate was lower than 0.6, using a realistic number of markers. Genealogical coancestry information is crucial in plant breeding to estimate genetic parameters and breeding values. The approach of Fernández and Toro (Mol Ecol 15:1657–1667, 2006) to estimate genealogical coancestries from molecular data through pedigree reconstruction was limited to species with separate sexes. In this study it was extended to plants, allowing hermaphroditism and monoecy, with possible selfing. Moreover, some improvements were made to take previous knowledge on the population demographic history into account. The new method was validated using simulated and real datasets. Simulations showed that accuracy of estimates was high with 30 microsatellites, with the best results obtained for selfing rates below 0.6. In these conditions, the root mean square error (RMSE) between the true and estimated genealogical coancestry was small (<0.07), although the number of ancestors was overestimated and the selfing rate could be biased. Simulations also showed that linkage disequilibrium between markers and departure from the Hardy–Weinberg equilibrium in the founder population did not affect the efficiency of the method. Real oil palm data confirmed the simulation results, with a high correlation between the true and estimated genealogical coancestry (>0.9) and a low RMSE (<0.08) using 38 markers. The method was applied to the Deli oil palm population for which pedigree data were scarce. The estimated genealogical coancestries were highly correlated (>0.9) with the molecular coancestries using 100 markers. Reconstructed pedigrees were used to estimate effective population sizes. In conclusion, this method gave reliable genealogical coancestry estimates. The strategy was implemented in the software MOLCOANC 3.0.     

Genetic variability in natural populations of Arabidopsis thaliana in northern Europe

Ten populations of the model plant Arabidopsis thaliana were collected along a north-south gradient in Norway and screened for microsatellite polymorphisms in 25 loci and variability in quantitative traits. Overall, the average levels of genetic diversity were found to be relatively high in these populations, compared to previously published surveys of within population variability. Six of the populations were polymorphic at microsatellite loci, resulting in an overall proportion of polymorphic loci of 18%, and a relatively high gene diversity for a selfing species (HE = 0.06). Of the overall variability, 12% was found within populations. Two of six polymorphic populations contained heterozygous individuals. Both FST and phylogenetic analyses showed no correlation between geographical and genetic distances. Haplotypic diversity patterns suggested postglacial colonization of Scandinavia from a number of different sources. Heritable variation was observed for many of the studied quantitative traits, with all populations showing variability in at least some traits, even populations with no microsatellite variability. There was a positive association between variability in quantitative traits and microsatellites within populations. Several quantitative traits exhibited QST values significantly less than FST, suggesting that selection may be acting to retard differentiation for these traits.     

Maximum likelihood estimation of population growth rates based on the coalescent.

We describe a method for co-estimating 4Nemu (four times the product of effective population size and neutral mutation rate) and population growth rate from sequence samples using Metropolis-Hastings sampling. Population growth (or decline) is assumed to be exponential. The estimates of growth rate are biased upwards, especially when 4Nemu is low; there is also a slight upwards bias in the estimate of 4Nemu itself due to correlation between the parameters. This bias cannot be attributed solely to Metropolis-Hastings sampling but appears to be an inherent property of the estimator and is expected to appear in any approach which estimates growth rate from genealogy structure. Sampling additional unlinked loci is much more effective in reducing the bias than increasing the number or length of sequences from the same locus.     

COVARIATION OF SELFING RATES WITH PARENTAL GENE FIXATION INDICES WITHIN POPULATIONS OF MIMULUS GUTTATUS

Wright's gene fixation index F and two single-locus effective selfing rates—the selfing rate at loci with fixed alleles, and the selfing rate at loci without fixed alleles—were estimated in five populations of Mimulus guttatus. These two effective selfing rates describe the inbreeding observed at a single locus when both uniparental and biparental inbreeding are practiced. Estimates were made using progeny arrays assayed for six allozyme loci and two morphological loci exhibiting dominance. The average of the two selfing rates computed for subpopulations (ca. 10 m diameter) ranged from 24% to 59%, with a mean of 37%. When computed for populations (ca. 1 km diameter), average selfing rates were about 10% higher. In four populations, the selfing rate at loci with fixed alleles was higher than the selfing rate at loci without fixed alleles. Thus, the covariance of selfing with parental gene fixation was positive. In one of the populations, estimates for individual plants sampled along a transect gave positive correlations for selfing rates and for gene-fixation indices between adjacent plants. A highly positive correlation between selfing rate and gene fixation of individual plants was also observed. In another population, the covariance of selfing with gene fixation was higher for a locus causing leaf spots than for allozyme loci. This covariance is partially caused by 1) variation in homozygosity among neighborhoods and 2) biparental inbreeding within neighborhoods. The consequences of this covariance are discussed.     

On estimating the proportion of variance in a phenotypic trait attributable to a measured locus

The measured genotype approach can be used to estimate the variance contributions of specific candidate loci to quantitative traits of interest. We show here that both the naive estimate of measured-locus heritability, obtained by invoking infinite-sample theory, and an estimate obtained from a bias-corrected variance estimate based on finite-sample theory, produce biased estimates of heritability. We identify the sources of bias, and quantify their effects. The two sources of bias are: (1) the estimation of heritability from population samples as the ratio of two variances, and (2) the existence of sampling error. We show that neither heritability estimator is less biased (in absolute value) than the other in all situations, and the choice of an ideal estimator is therefore a function of the sample size and magnitude of the locus-specific contribution to the overall phenotypic variance. In most cases the bias is small, so that the practical implications of using either estimator are expected to be minimal.     

Deleterious Mutations as an Evolutionary Factor. II. Facultative Apomixis and Selfing

A population with u deleterious mutations per genome per generation is considered in which only those individuals that carry less than a critical number k of mutations are viable. Besides a large number of loci subject to mutation and selection, the genome contains one or two special loci responsible for the mode of reproduction. Amphimixis vs. apomixis and amphimixis vs. selfing are considered separately. In the first case, the genome degradation rate v (= u/square root k) is found to play the decisive role, as in the case of recombination. When v greater than 1.25, obligate amphimixis is established. If v decreases below this value, the alleles with first low and then larger penetrance are fixed, until alleles conferring obligate asexual reproduction become advantageous. The proportion of resources allocated to produce seeds also increases with decrease of v. These results are unlikely to depend on the genetic basis of the mode of reproduction. The result of competition between outcrossing and selfing depends on both u and k, as well as on whether the mutations are recessive. The alleles for selfing with low penetrance are selected against if the mutations are at all recessive. The fitness of alleles with high penetrance depends primarily on u, decreasing when u increases. There may exist conditions when only the alleles providing intermediate selfing rates can be fixed in a population. In other cases a population may exist with either obligate outcrossing or selfing at a high rate. Thus, truncation selection against deleterious mutations may be a factor supporting obligate or facultative sex despite the twofold advantage of apomixis or selfing.     

Modes and rates of selfing and associated inbreeding depression in the self-incompatible plant Senecio squalidus (Asteraceae): a successful colonizing species in the British Isles

The strength of the self-incompatibility (SI) response in Senecio squalidus was measured across its British range. Geographic variation in SI was investigated and the extent and inheritance of pseudo-self-compatibility (PSC) and inbreeding depression were determined. Mean self-fruit-set per capitulum was calculated for individuals and sample populations. The heritability of PSC and the magnitude of inbreeding depression were assessed by comparing selfing rates and fitness trait values between SI and PSC parent-progeny lines. SI was found to be strongly expressed in S. squalidus throughout its British range, with only 3.1% of the individuals sampled showing PSC. This PSC had relatively low heritability with stronger expression of SI in selfed progeny relative to PSC parents. Inbreeding depression was shown to be great in S. squalidus, with mean life history stage values ranging from 0.18 to 0.25. The strength of SI in S. squalidus appears not to have weakened in response to its rapid colonization of Britain. The avoidance of inbreeding depression is likely to be the primary factor maintaining strong SI in this successful colonizing species.     

How geitonogamous selfing affects sex allocation in hermaphrodite plants

Does the mode of self-pollination affect the evolutionarily stable allocation to male vs. female function? We distinguish the following scenarios. (1) An ‘autogamous’ species, in which selfing occurs within the flower prior to opening. The pollen used in selfing is a constant fraction of all pollen grains produced. (2) A species with ‘abiotic pollination’, in which selfing occurs when pollen dispersed in one flower lands on the stigma of a nearby flower on the same plant (geitonogamy). The selfing rate increases with male allocation but a higher selfing rate does not mean a reduced export of pollen. (3) An ‘animal-pollinated’ species with geitonogamous selfing. Here the selfing rate also increases with male allocation, but pollen export to other plants in the population is a decelerating function of the number of simultaneously open flowers. In all three models selfing selects for increased female allocation. For model 3 this contradicts the general opinion that geitonogamous selfing does not affect evolutionarily stable allocations. In all models, the parent benefits more from a female-biased allocation than any other individual in the population. In addition, in models 2 and 3, greater male allocation results in more local mate competition. In model 3 and in model 2 with low levels of inbreeding depression, hermaphroditism is evolutionarily stable. In model 2 with high inbreeding depression, the population converges to a fitness minimum for the relative allocation to male function. In this case the fitness set is bowed inwards, corresponding with accelerating fitness gain curves. If the selfing rate increases with plant size, this is a sufficient condition for size-dependent sex allocation (more allocation towards seeds in large plants) to evolve. We discuss our results in relation to size-dependent sex allocation in plants and in relation to the evolution of dioecy.