THE CHANGE IN QUANTITATIVE GENETIC VARIATION WITH INBREEDING

The potential for floral evolution in response to competing selection pressures following the loss of hawkmoth pollination in Ruellia humilis

PREMISE

In the absence of hawkmoth pollinators, chasmogamous (CH) flowers of Ruellia humilis self-pollinate by two secondary mechanisms. Other floral visitors might exert selection on CH floral traits to restore outcrossing, but at the same time preferential predation of CH seeds generates selection to increase the allocation of resources to cleistogamous (CL) flowers.

METHODS

To assess the potential for an evolutionary response to these competing selection pressures, we estimated additive genetic variances ( σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ ) and covariances for 14 reproductive traits and three fitness components in a Missouri population lacking hawkmoth pollinators.

RESULTS

We found significant σ A 2 ${\sigma }_{{\rm{A}}}^{2}$ for all 11 floral traits and two measures of resource allocation to CL flowers, indicating the potential for a short-term response to selection on most reproductive traits. Selection generated by seed predators is predicted to increase the percentage of CL flowers by 0.24% per generation, and mean stigma-anther separation is predicted to decrease as a correlated response, increasing the fraction of plants that engage in prior selfing. However, the initial response to this selection is opposed by strong directional dominance.

CONCLUSIONS

The predicted evolutionary decrease in the number of CH flowers available for potential outcrossing, combined with the apparent preclusion of potential diurnal pollinators by the pollen-harvesting activities of sweat bees, suggest that 100% cleistogamy is the likely outcome of evolution in the absence of hawkmoths. However, rare mutations with large effects, such as delaying budbreak until after sunrise, could provide pathways for the restoration of outcrossing that are not reachable by gradual quantitative-genetic evolution.

Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species

About 50 y ago, Crow and Kimura [An Introduction to Population Genetics Theory (1970)] and Ohta and Kimura [Genet. Res. 22, 201-204 (1973)] laid the foundations of conservation genetics by predicting the relationship between population size and genetic marker diversity. This work sparked an enormous research effort investigating the importance of population dynamics, in particular small population size, for population mean performance, population viability, and evolutionary potential. In light of a recent perspective [J. C. Teixeira, C. D. Huber, Proc. Natl. Acad. Sci. U.S.A. 118, 10 (2021)] that challenges some fundamental assumptions in conservation genetics, it is timely to summarize what the field has achieved, what robust patterns have emerged, and worthwhile future research directions. We consider theory and methodological breakthroughs that have helped management, and we outline some fundamental and applied challenges for conservation genetics.

Gene action, genetic variation, and GWAS: A user-friendly web tool

Fisher's partitioning of genotypic values and genetic variance is highly relevant in the current era of genome-wide association studies (GWASs). However, despite being more than a century old, a number of persistent misconceptions related to nonadditive genetic effects remain. We developed a user-friendly web tool, the Falconer ShinyApp, to show how the combination of gene action and allele frequencies at causal loci translate to genetic variance and genetic variance components for a complex trait. The app can be used to demonstrate the relationship between a SNP effect size estimated from GWAS and the variation the SNP generates in the population, i.e., how locus-specific effects lead to individual differences in traits. In addition, it can also be used to demonstrate how within and between locus interactions (dominance and epistasis, respectively) usually do not lead to a large amount of nonadditive variance relative to additive variance, and therefore, that these interactions usually do not explain individual differences in a population.

The Contemporary Evolution of Fitness

The rate of evolution of population mean fitness informs how selection acting in contemporary populations can counteract environmental change and genetic degradation (mutation, gene flow, drift, recombination). This rate influences population increases (e.g., range expansion), population stability (e.g., cryptic eco-evolutionary dynamics), and population recovery (i.e., evolutionary rescue). We review approaches for estimating such rates, especially in wild populations. We then review empirical estimates derived from two approaches: mutation accumulation (MA) and additive genetic variance in fitness (IAw). MA studies inform how selection counters genetic degradation arising from deleterious mutations, typically generating estimates of <1% per generation. IAw studies provide an integrated prediction of proportional change per generation, nearly always generating estimates of <20% and, more typically, <10%. Overall, considerable, but not unlimited, evolutionary potential exists in populations facing detrimental environmental or genetic change. However, further studies with diverse methods and species are required for more robust and general insights.

The influence of genetic drift and selection on quantitative traits in a plant pathogenic fungus

Genetic drift and selection are ubiquitous evolutionary forces acting to shape genetic variation in populations. While their relative importance has been well studied in plants and animals, less is known about their relative importance in fungal pathogens. Because agro-ecosystems are more homogeneous environments than natural ecosystems, stabilizing selection may play a stronger role than genetic drift or diversifying selection in shaping genetic variation among populations of fungal pathogens in agro-ecosystems. We tested this hypothesis by conducting a Q_ST/F_ST analysis using agricultural populations of the barley pathogen Rhynchosporium commune. Population divergence for eight quantitative traits (Q_ST) was compared with divergence at eight neutral microsatellite loci (F_ST) for 126 pathogen strains originating from nine globally distributed field populations to infer the effects of genetic drift and types of selection acting on each trait. Our analyses indicated that five of the eight traits had Q_ST values significantly lower than F_ST, consistent with stabilizing selection, whereas one trait, growth under heat stress (22°C), showed evidence of diversifying selection and local adaptation (Q_ST>F_ST). Estimates of heritability were high for all traits (means ranging between 0.55–0.84), and average heritability across traits was negatively correlated with microsatellite gene diversity. Some trait pairs were genetically correlated and there was significant evidence for a trade-off between spore size and spore number, and between melanization and growth under benign temperature. Our findings indicate that many ecologically and agriculturally important traits are under stabilizing selection in R. commune and that high within-population genetic variation is maintained for these traits.

Pedigree-free estimates of heritability in the wild: promising prospects for selfing populations

Estimating the genetic variance available for traits informs us about a population's ability to evolve in response to novel selective challenges. In selfing species, theory predicts a loss of genetic diversity that could lead to an evolutionary dead-end, but empirical support remains scarce. Genetic variability in a trait is estimated by correlating the phenotypic resemblance with the proportion of the genome that two relatives share identical by descent ('realized relatedness'). The latter is traditionally predicted from pedigrees (Φ A : expected value) but can also be estimated using molecular markers (average number of alleles shared). Nevertheless, evolutionary biologists, unlike animal breeders, remain cautious about using marker-based relatedness coefficients to study complex phenotypic traits in populations. In this paper, we review published results comparing five different pedigree-free methods and use simulations to test individual-based models (hereafter called animal models) using marker-based relatedness coefficients, with a special focus on the influence of mating systems. Our literature review confirms that Ritland's regression method is unreliable, but suggests that animal models with marker-based estimates of relatedness and genomic selection are promising and that more testing is required. Our simulations show that using molecular markers instead of pedigrees in animal models seriously worsens the estimation of heritability in outcrossing populations, unless a very large number of loci is available. In selfing populations the results are less biased. More generally, populations with high identity disequilibrium (consanguineous or bottlenecked populations) could be propitious for using marker-based animal models, but are also more likely to deviate from the standard assumptions of quantitative genetics models (non-additive variance).

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Anthropogenically induced adaptation to invade (AIAI): contemporary adaptation to human-altered habitats within the native range can promote invasions

Adaptive evolution is currently accepted as playing a significant role in biological invasions. Adaptations relevant to invasions are typically thought to occur either recently within the introduced range, as an evolutionary response to novel selection regimes, or within the native range, because of long-term adaptation to the local environment. We propose that recent adaptation within the native range, in particular adaptations to human-altered habitat, could also contribute to the evolution of invasive populations. Populations adapted to human-altered habitats in the native range are likely to increase in abundance within areas frequented by humans and associated with human transport mechanisms, thus enhancing the likelihood of transport to a novel range. Given that habitats are altered by humans in similar ways worldwide, as evidenced by global environmental homogenization, propagules from populations adapted to human-altered habitats in the native range should perform well within similarly human-altered habitats in the novel range. We label this scenario 'Anthropogenically Induced Adaptation to Invade'. We illustrate how it differs from other evolutionary processes that may occur during invasions, and how it can help explain accelerating rates of invasions.

Genetic erosion impedes adaptive responses to stressful environments

Biodiversity is increasingly subjected to human-induced changes of the environment. To persist, populations continually have to adapt to these often stressful changes including pollution and climate change. Genetic erosion in small populations, owing to fragmentation of natural habitats, is expected to obstruct such adaptive responses: (i) genetic drift will cause a decrease in the level of adaptive genetic variation, thereby limiting evolutionary responses; (ii) inbreeding and the concomitant inbreeding depression will reduce individual fitness and, consequently, the tolerance of populations to environmental stress. Importantly, inbreeding generally increases the sensitivity of a population to stress, thereby increasing the amount of inbreeding depression. As adaptation to stress is most often accompanied by increased mortality (cost of selection), the increase in the ‘cost of inbreeding’ under stress is expected to severely hamper evolutionary adaptive processes. Inbreeding thus plays a pivotal role in this process and is expected to limit the probability of genetically eroded populations to successfully adapt to stressful environmental conditions. Consequently, the dynamics of small fragmented populations may differ considerably from large nonfragmented populations. The resilience of fragmented populations to changing and deteriorating environments is expected to be greatly decreased. Alleviating inbreeding depression, therefore, is crucial to ensure population persistence.

The consequences of rare sexual reproduction by means of selfing in an otherwise clonally reproducing species

Clonal reproduction of diploids leads to an increase in heterozygosity over time. A single round of selfing will then create new homozygotic genotypes. Given the same allele frequencies, heritable genetic variation is larger when there are more extreme, i.e. homozygotic genotypes. So after a long clonal expansion, one round of selfing increases heritable genetic variation, but any fully or partially recessive deleterious alleles simultaneously impose a fitness cost. Here we calculate that the cost of selfing in the yeast Saccharomyces is experienced only by a minority of zygotes. This allows a round of selfing to act as an evolutionary capacitor to unlock genetic variation previously found in a cryptic heterozygous form. We calculate the evolutionary consequences rather than the evolutionary causes of sex. We explore a range of parameter values describing sexual frequencies, focusing especially on the parameter values known for wild Saccharomyces. Our results are largely robust to many other parameter value choices, so long as meiosis is rare relative to the strength of selection on heterozygotes. Results may also be limited to organisms with a small number of genes. We therefore expect the same phenomenon in some other species with similar reproductive strategies.

The effect of a population bottleneck on the evolution of genetic variance/covariance structure

It is well known that standard population genetic theory predicts decreased additive genetic variance (V(a) ) following a population bottleneck and that theoretical models including interallelic and intergenic interactions indicate such loss may be avoided. However, few empirical data from multicellular model systems are available, especially regarding variance/covariance (V/CV) relationships. Here, we compare the V/CV structure of seventeen traits related to body size and composition between control (60 mating pairs/generation) and bottlenecked (2 mating pairs/generation; average F = 0.39) strains of mice. Although results for individual traits vary considerably, multivariate analysis indicates that V(a) in the bottlenecked populations is greater than expected. Traits with patterns and amounts of epistasis predictive of enhanced V(a) also show the largest deviations from additive expectations. Finally, the correlation structure of weekly weights is not significantly different between control and experimental lines but correlations between necropsy traits do differ, especially those involving the heart, kidney and tail length.

Predicting the probability of outbreeding depression

Fragmentation of animal and plant populations typically leads to genetic erosion and increased probability of extirpation. Although these effects can usually be reversed by re-establishing gene flow between population fragments, managers sometimes fail to do so due to fears of outbreeding depression (OD). Rapid development of OD is due primarily to adaptive differentiation from selection or fixation of chromosomal variants. Fixed chromosomal variants can be detected empirically. We used an extended form of the breeders' equation to predict the probability of OD due to adaptive differentiation between recently isolated population fragments as a function of intensity of selection, genetic diversity, effective population sizes, and generations of isolation. Empirical data indicated that populations in similar environments had not developed OD even after thousands of generations of isolation. To predict the probability of OD, we developed a decision tree that was based on the four variables from the breeders' equation, taxonomic status, and gene flow within the last 500 years. The predicted probability of OD in crosses between two populations is elevated when the populations have at least one of the following characteristics: are distinct species, have fixed chromosomal differences, exchanged no genes in the last 500 years, or inhabit different environments. Conversely, the predicted probability of OD in crosses between two populations of the same species is low for populations with the same karyotype, isolated for <500 years, and that occupy similar environments. In the former case, we recommend crossing be avoided or tried on a limited, experimental basis. In the latter case, crossing can be carried out with low probability of OD. We used crosses with known results to test the decision tree and found that it correctly identified cases where OD occurred. Current concerns about OD in recently fragmented populations are almost certainly excessive.

The relative importance of factors determining genetic drift: mating system, spatial genetic structure, habitat and census size in Arabidopsis lyrata

• The mating system, dispersal and census size are predicted to determine the magnitude of genetic drift, but little is known about their relative importance in nature. • We estimated the contributions of several population-level features to genetic drift in 18 populations of Arabidopsis lyrata. The factors were outcrossing rate, within-population spatial genetic structure, census size and substrate type. The expected heterozygosity (H(E)) at 10 microsatellite loci was taken to reflect the effective population size (N(e)) and the strength of genetic drift. • The mating system explained most of the variation in H(E) (60%), followed by substrate (10%), genetic structure (9%) and census size (6%). The most outcrossing population had a +0.32 higher predicted H(E) than the most selfing population; the estimated N(e) of selfing populations was less than half that of outcrossing populations. Rocky outcrops supported populations with a +0.14 higher H(E) than did sandy substrates. The most structured population had a +0.24 higher H(E) than the least structured population, and the largest population had a +0.18 higher H(E) than the smallest population. • This study illustrates the importance of outcrossing, genetic structure and the physical environment--together with census size--in maintaining H(E), and suggests that multiple population-level characteristics influence N(e) and the action of genetic drift.

Patterns of genetic variation and covariation in ejaculate traits reveal potential evolutionary constraints in guppies

Ejaculates comprise multiple and potentially interacting traits that determine male fertility and sperm competitiveness. Consequently, selection on these traits is likely to be intense, but the efficacy of selection will depend critically on patterns of genetic variation and covariation underlying their expression. In this study, I provide a prospective quantitative genetic analysis of ejaculate traits in the guppy Poecilia reticulata, a highly promiscuous live-bearing fish. I used a standard paternal half-sibling breeding design to characterize patterns of genetic (co)variation in components of sperm length and in vitro sperm performance. All traits exhibited high levels of phenotypic and additive genetic variation, and in several cases, patterns of genetic variation was consistent with Y-linkage. There were also highly significant negative genetic correlations between the various measures of sperm length and sperm performance. In particular, the length of the sperm's midpiece was strongly, negatively and genetically correlated with sperm's swimming velocity-an important determinant of sperm competitiveness in this and other species. Other components of sperm length, including the flagellum and head, were independently and negatively genetically correlated with the proportion of live sperm in the ejaculate (sperm viability). Whether these relationships represent evolutionary trade-offs depends on the precise relationships between these traits and competitive fertilization rates, which have yet to be fully resolved in this (and indeed most) species. Nevertheless, these prospective analyses point to potential constraints on ejaculate evolution and may explain the high level of phenotypic variability in ejaculate traits in this species.

Demographic factors and genetic variation influence population persistence under environmental change

Population persistence has been studied in a conservation context to predict the fate of small or declining populations. Persistence models have explored effects on extinction of random demographic and environmental fluctuations, but in the face of directional environmental change they should also integrate factors affecting whether a population can adapt. Here, we examine the population-size dependence of demographic and genetic factors and their likely contributions to extinction time under scenarios of environmental change. Parameter estimates were derived from experimental populations of the rainforest species, Drosophila birchii, held in the lab for 10 generations at census sizes of 20, 100 and 1000, and later exposed to five generations of heat-knockdown selection. Under a model of directional change in the thermal environment, rapid extinction of populations of size 20 was caused by a combination of low growth rate (r) and high stochasticity in r. Populations of 100 had significantly higher reproductive output, lower stochasticity in r and more additive genetic variance (V(A)) than populations of 20, but they were predicted to persist less well than the largest size class. Even populations of 1000 persisted only a few hundred generations under realistic estimates of environmental change because of low V(A) for heat-knockdown resistance. The experimental results document population-size dependence of demographic and adaptability factors. The simulations illustrate a threshold influence of demographic factors on population persistence, while genetic variance has a more elastic impact on persistence under environmental change.

Individual phenotype, kinship, and the occurrence of inbreeding in song sparrows

Inbreeding load, a key parameter in evolutionary ecology, is frequently estimated by regressing fitness (or related traits) on inbreeding coefficient across population members. This approach assumes that inbreeding occurs randomly with respect to an individual's intrinsic ability to produce fit offspring; estimated loads might otherwise be biased by covariation between inbreeding and individual quality. This assumption, however, is rarely validated. We tested whether, in free-living song sparrows Melospiza melodia, an individual's observed kinship with its social mate (and hence the degree of inbreeding in which an individual participated) was correlated with specific phenotypic traits that are likely to indicate individual quality. Males (and to some extent females) that hatched earlier within their cohort, had shorter tarsi, bred earlier during their first year, or survived fewer years paired with more closely related mates and therefore produced relatively inbred offspring. These correlations arose because males with specific phenotypes were more closely related to the female population (and therefore more likely to pair with closer relatives under random pairing), and because males with specific phenotypes paired with closer relatives than expected. Such correlations could bias estimated inbreeding loads, and should be considered in quantitative genetic analyses of phenotypic variance in populations in which inbreeding occurs.

Mate choice for non-additive genetic benefits: a resolution to the lek paradox

In promiscuous mating systems, females often show a consistent preference to mate with one or a few males, presumably to acquire heritable genetic benefits for their offspring. However, strong directional selection should deplete additive genetic variation in fitness and consequently any benefit to expressing the preference by females (referred to as the lek paradox). Here, we provide a novel resolution that examines non-additive genetic benefits, such as overdominance or inbreeding, as a source of genetic variation. Focusing on the inbreeding coefficient f and overdominance effects, we use dynamic models to show that (1) f can be inherited from sire to offspring, (2) populations with females that express a mating preferences for outbred males (low f) maintain higher genetic variation than populations with females that mate randomly, and (3) preference alleles for outbred males can invade populations even when the alleles are associated with a fecundity cost. We show that non-additive genetic variation due to overdominance can be converted to additive genetic variation and becomes "heritable" when the frequencies of alternative homozygous genotypes at fitness loci deviate from equality. Unlike previous models that assume an infinite population size, we now show that genetic drift in finite populations can lead to the necessary deviations in the frequencies of homozygous genotypes. We also show that the "heritability of f," and hence the benefit to a mating preference for non-additive genetic benefits, is highest in small populations and populations in which a smaller number of loci contribute to fitness via overdominance. Our model contributes to the solution of the lek paradox.

The role of gene flow in shaping genetic structures of the subtropical conifer species Araucaria angustifolia

The morphological features of pollen and seed of Araucaria angustifolia have led to the proposal of limited gene dispersal for this species. We used nuclear microsatellite and AFLP markers to assess patterns of genetic variation in six natural populations at the intra- and inter-population level, and related our findings to gene dispersal in this species. Estimates of both fine-scale spatial genetic structure (SGS) and migration rate suggest relatively short-distance gene dispersal. However, gene dispersal differed among populations, and effects of more efficient dispersal within population were observed in at least one stand. In addition, even though some seed dispersal may be aggregated in this principally barochorous species, reasonable secondary seed dispersal, presumably facilitated by animals, and overlap of seed shadows within populations is suggested. Overall, no correlation was observed between levels of SGS and inbreeding, density or age structure, except that a higher level of SGS was revealed for the population with a higher number of juvenile individuals. A low estimate for the number of migrants per generation between two neighbouring populations implies limited gene flow. We expect that stepping-stone pollen flow may have contributed to low genetic differentiation among populations observed in a previous survey. Thus, strategies for maintenance of gene flow among remnant populations should be considered in order to avoid degrading effects of population fragmentation on the evolution of A. angustifolia.

Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions

Invasive species are predicted to suffer from reductions in genetic diversity during founding events, reducing adaptive potential. Integrating evidence from two literature reviews and two case studies, we address the following questions: How much genetic diversity is lost in invasions? Do multiple introductions ameliorate this loss? Is there evidence for loss of diversity in quantitative traits? Do invaders that have experienced strong bottlenecks show adaptive evolution? How do multiple introductions influence adaptation on a landscape scale? We reviewed studies of 80 species of animals, plants, and fungi that quantified nuclear molecular diversity within introduced and source populations. Overall, there were significant losses of both allelic richness and heterozygosity in introduced populations, and large gains in diversity were rare. Evidence for multiple introductions was associated with increased diversity, and allelic variation appeared to increase over long timescales (~100 years), suggesting a role for gene flow in augmenting diversity over the long-term. We then reviewed the literature on quantitative trait diversity and found that broad-sense variation rarely declines in introductions, but direct comparisons of additive variance were lacking. Our studies of Hypericum canariense invasions illustrate how populations with diminished diversity may still evolve rapidly. Given the prevalence of genetic bottlenecks in successful invading populations and the potential for adaptive evolution in quantitative traits, we suggest that the disadvantages associated with founding events may have been overstated. However, our work on the successful invader Verbascum thapsus illustrates how multiple introductions may take time to commingle, instead persisting as a 'mosaic of maladaptation' where traits are not distributed in a pattern consistent with adaptation. We conclude that management limiting gene flow among introduced populations may reduce adaptive potential but is unlikely to prevent expansion or the evolution of novel invasive behaviour.