Population bottlenecks increase additive genetic variance but do not break a selection limit in rain forest Drosophila

According to neutral quantitative genetic theory, population bottlenecks are expected to decrease standing levels of additive genetic variance of quantitative traits. However, some empirical and theoretical results suggest that, if nonadditive genetic effects influence the trait, bottlenecks may actually increase additive genetic variance. This has been an important issue in conservation genetics where it has been suggested that small population size might actually experience an increase rather than a decrease in the rate of adaptation. Here we test if bottlenecks can break a selection limit for desiccation resistance in the rain forest-restricted fly Drosophila bunnanda. After one generation of single-pair mating, additive genetic variance for desiccation resistance increased to a significant level, on average higher than for the control lines. Line crosses revealed that both dominance and epistatic effects were responsible for the divergence in desiccation resistance between the original control and a bottlenecked line exhibiting increased additive genetic variance for desiccation resistance. However, when bottlenecked lines were selected for increased desiccation resistance, there was only a small shift in resistance, much less than predicted by the released additive genetic variance. The small selection response in the bottlenecked lines was no greater than that observed in the control lines. Thus bottlenecks might produce a statistically detectable change in additive genetic variance but this change has no impact on the response to selection.     

The change in quantitative genetic variation with inbreeding

Inbreeding is known to reduce heterozygosity of neutral genetic markers, but its impact on quantitative genetic variation is debated. Theory predicts a linear decline in additive genetic variance (V(A)) with increasing inbreeding coefficient (F) when loci underlying the trait act additively, but a nonlinear hump-shaped relationship when dominance and epistasis are important. Predictions for heritability (h2) are similar, although the exact shape depends on the value of h2 in the absence of inbreeding. We located 22 published studies in which the level of genetic variation in experimentally inbred populations (measured by V(A) or h2) was compared with that in outbred control populations. For life-history traits, the data strongly supported a nonlinear change in genetic variation with increasing F. V(A) and h2 were, respectively, 244% and 50% higher at F = 0.4 than in outbred populations, and dominance plus epistatic variance together exceeded additive variance by a factor of four. For nonfitness traits the decline was linear and estimates of nonadditive variance were small. These results confirm that population bottlenecks frequently increase V(A) in some traits, and imply that life-history traits are underlain by substantial dominance or epistasis. However, the importance of drift-induced genetic variation in conservation or evolutionary biology is questionable, in part because inbreeding depression usually accompanies inbreeding.     

EPISTASIS AND THE EVOLUTION OF ADDITIVE GENETIC VARIANCE IN POPULATIONS THAT PASS THROUGH A BOTTLENECK

Traditional models of genetic drift predict a linear decrease in additive genetic variance for populations passing through a bottleneck. This perceived lack of heritable variance limits the scope of founder-effect models of speciation. We produced 55 replicate bottleneck populations maintained at two male-female pairs through four generations of inbreeding (average F = 0.39). These populations were formed from an F₂ intercross of the LG/J and SM/J inbred mouse strains. Two contemporaneous control strains maintained with more than 60 mating pairs per generation were formed from this same source population. The average level of within-strain additive genetic variance for adult body weight was compared between the control and experimental lines. Additive genetic variance for adult body weight within experimental bottleneck strains was significantly higher than expected under an additive genetic model This enhancement of additive genetic variance under inbreeding is likely to be due to epistasis, which retards or reverses the loss of additive genetic variance under inbreeding for adult body weight in this population. Therefore, founder-effect speciation processes may not be constrained by a loss of heritable variance due to population bottlenecks.     

Comparing Mutational Variabilities

We have reviewed the available data on V(M), the amount of genetic variation in phenotypic traits produced each generation by mutation. We use these data to make several qualitative tests of the mutation-selection balance hypothesis for the maintenance of genetic variance (MSB). To compare V(M) values, we use three dimensionless quantities: mutational heritability, V(M)/V(E); the mutational coefficient of variation, CV(M); and the ratio of the standing genetic variance to V(M), V(G)/V(M). Since genetic coefficients of variation for life history traits are larger than those for morphological traits, we predict that under MSB, life history traits should also have larger CV(M). This is confirmed; life history traits have a median CV(M) value more than six times higher than that for morphological traits. V(G)/V(M) approximates the persistence time of mutations under MSB in an infinite population. In order for MSB to hold, V(G)/V(M) must be small, substantially less than 1000, and life history traits should have smaller values than morphological traits. V(G)/V(M) averages about 50 generations for life history traits and 100 generations for morphological traits. These observations are all consistent with the predictions of a mutation-selection balance model.     

The Effect of an Experimental Bottleneck upon Quantitative Genetic Variation in the Housefly

Effects of a population bottleneck (founder-flush cycle) upon quantitative genetic variation of morphometric traits were examined in replicated experimental lines of the housefly founded with one, four or 16 pairs of flies. Heritability and additive genetic variances for eight morphometric traits generally increased as a result of the bottleneck, but the pattern of increase among bottleneck sizes differed among traits. Principal axes of the additive genetic correlation matrix for the control line yielded two suites of traits, one associated with general body size and another set largely independent of body size. In the former set containing five of the traits, additive genetic variance was greatest in the bottleneck size of four pairs, whereas in the latter set of two traits the largest additive genetic variance occurred in the smallest bottleneck size of one pair. One trait exhibited changes in additive genetic variance intermediate between these two major responses. These results were inconsistent with models of additive effects of alleles within loci or of additive effects among loci. An observed decline in viability measures and body size in the bottleneck lines also indicated that there was nonadditivity of allelic effects for these traits. Several possible nonadditive models were explored that increased additive genetic variance as a result of a bottleneck. These included a model with complete dominance, a model with overdominance and a model incorporating multiplicative epistasis.     

Effects of bottlenecks on quantitative genetic variation in the butterfly Bicyclus anynana

The effects of a single population bottleneck of differing severity on heritability and additive genetic variance was investigated experimentally using a butterfly. An outbred laboratory stock was used to found replicate lines with one pair, three pairs and 10 pairs of adults, as well as control lines with approximately 75 effective pairs. Heritability and additive genetic variance of eight wing pattern characters and wing size were estimated using parent-offspring covariances in the base population and in all daughter lines. Individual morphological characters and principal components of the nine characters showed a consistent pattern of treatment effects in which average heritability and additive genetic variance was lower in one pair and three pair lines than in 10 pair and control lines. Observed losses in heritability and additive genetic variance were significantly greater than predicted by the neutral additive model when calculated with coefficients of inbreeding estimated from demographic parameters alone. However, use of molecular markers revealed substantially more inbreeding, generated by increased variance in family size and background selection. Conservative interpretation of a statistical analysis incorporating this previously undetected inbreeding led to the conclusion that the response to inbreeding of the morphological traits studied showed no significant departure from the neutral additive model. This result is consistent with the evidence for minimal directional dominance for these traits. In contrast, egg hatching rate in the same experimental lines showed strong inbreeding depression, increased phenotypic variance and rapid response to selection, highly indicative of an increase in additive genetic variance due to dominance variance conversion.     

Epistasis Is a Major Determinant of the Additive Genetic Variance in Mimulus guttatus

The influence of genetic interactions (epistasis) on the genetic variance of quantitative traits is a major unresolved problem relevant to medical, agricultural, and evolutionary genetics. The additive genetic component is typically a high proportion of the total genetic variance in quantitative traits, despite that underlying genes must interact to determine phenotype. This study estimates direct and interaction effects for 11 pairs of Quantitative Trait Loci (QTLs) affecting floral traits within a single population of Mimulus guttatus. With estimates of all 9 genotypes for each QTL pair, we are able to map from QTL effects to variance components as a function of population allele frequencies, and thus predict changes in variance components as allele frequencies change. This mapping requires an analytical framework that properly accounts for bias introduced by estimation errors. We find that even with abundant interactions between QTLs, most of the genetic variance is likely to be additive. However, the strong dependency of allelic average effects on genetic background implies that epistasis is a major determinant of the additive genetic variance, and thus, the population’s ability to respond to selection.     

AN ANALYSIS OF SELECTIONAL RESPONSE IN RELATION TO A POPULATION BOTTLENECK

Selection for increased morphometric shape (ratio of wing length to thorax width) was compared between control (nonbottlenecked) populations and bottlenecked populations founded with two male–female pairs of flies. Contrary to neutral expectation, selectional response was not reduced in bottlenecked populations, and the mean realized heritabilities and additive genetic variances were higher for the bottlenecked lines than for the nonbottlenecked lines. Additive genetic variances based on these realized heritabilities were consistent with independent estimates of genetic variances based on parent–offspring covariances. Joint scaling tests applied to the crosses between selected lines and their controls revealed significant nonadditive components of genetic variance in the ancestor, which were not detected in the crosses involving bottlenecked lines. The nonbottlenecked lines responded principally by changes in one trait or the other (wing length or thorax width) but not in both, and regardless of which trait responded, larger trait size was dominant and epistatic to smaller size. Stabilizing selection for morphometric shape in the ancestor likely molded the genetic architecture to include nonadditive genetic effects.     

Heritability and genetic correlations of personality traits in a wild population of yellow‐bellied marmots (Marmota flaviventris)

Describing and quantifying animal personality is now an integral part of behavioural studies because individually distinctive behaviours have ecological and evolutionary consequences. Yet, to fully understand how personality traits may respond to selection, one must understand the underlying heritability and genetic correlations between traits. Previous studies have reported a moderate degree of heritability of personality traits, but few of these studies have either been conducted in the wild or estimated the genetic correlations between personality traits. Estimating the additive genetic variance and covariance in the wild is crucial to understand the evolutionary potential of behavioural traits. Enhanced environmental variation could reduce heritability and genetic correlations, thus leading to different evolutionary predictions. We estimated the additive genetic variance and covariance of docility in the trap, sociability (mirror image stimulation), and exploration and activity in two different contexts (open‐field and mirror image simulation experiments) in a wild population of yellow‐bellied marmots (Marmota flaviventris). We estimated both heritability of behaviours and of personality traits and found nonzero additive genetic variance in these traits. We also found nonzero maternal, permanent environment and year effects. Finally, we found four phenotypic correlations between traits, and one positive genetic correlation between activity in the open‐field test and sociability. We also found permanent environment correlations between activity in both tests and docility and exploration in the MIS test. This is one of a handful of studies to adopt a quantitative genetic approach to explain variation in personality traits in the wild and, thus, provides important insights into the potential variance available for selection.     

Comparing evolvabilities: common errors surrounding the calculation and use of coefficients of additive genetic variation

In 1992, David Houle showed that measures of additive genetic variation standardized by the trait mean, CV(A) (the coefficient of additive genetic variation) and its square (I(A) ), are suitable measures of evolvability. CV(A) has been used widely to compare patterns of genetic variation. However, the use of CV(A) s for comparative purposes relies critically on the correct calculation of this parameter. We reviewed a sample of quantitative genetic studies, focusing on sire models, and found that 45% of studies use incorrect methods for calculating CV(A) and that practices that render these coefficients meaningless are frequent. This may have important consequences for conclusions drawn from comparative studies. Our results are suggestive of a broader problem because miscalculation of the additive genetic variance from a sire model is prevalent among the studies sampled, implying that other important quantitative genetic parameters might also often be estimated incorrectly. We discuss the most prominent issues affecting the use of CV(A) and I(A) , including scale effects, data transformation, and the comparison of traits with different dimensions. Our aim is to increase awareness of the potential mistakes surrounding the calculation and use of evolvabilities, and to compile general guidelines for calculating, reporting, and interpreting these useful measures in future studies.     

Effects of genetic drift on variance components under a general model of epistasis

We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(delta(z)) = sigma(n)k=1 FkV(A(k)), where n is the number of loci contributing to the trait variance, V(A(1)) = V(A) is the additive genetic variance, and V(A(k)) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted (V*(A)), is closely related to var(delta(z)); (V*(A)) = (1 - F) sigma(n)k=1 kFk-1V(A(k)). Thus, epistasis inflates the expected additive variance above V(A)(1 - F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates (V*(AA. This is not true in general with diploidy, because dominance alone can reduce (V*(A)) below V(A)(1 - F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(delta(z)) = sigma(n)k=1 (2F)kV(A(k)) and (V(A)) = (1 - F) sigma(n)k=1 k(2F)k-1V(A(k)). With dominance (and even without epistasis), var(delta(z)) and (V*(A)) no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to (V*(A)) approximately equal to (1 - F)V(A) + 4FV(AA) + 2FV(D) + 2FC(AD), where C(AD) is a sum of two terms describing covariances between additive effects and dominance and additive X dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase V(A), we question the evolutionary importance of such increases.     

Selection and genetic (co)variance in bighorn sheep

Genetic theory predicts that directional selection should deplete additive genetic variance for traits closely related to fitness, and may favor the maintenance of alleles with antagonistically pleiotropic effects on fitness-related traits. Trait heritability is therefore expected to decline with the degree of association with fitness, and some genetic correlations between selected traits are expected to be negative. Here we demonstrate a negative relationship between trait heritability and association with lifetime reproductive success in a wild population of bighorn sheep (Ovis canadensis) at Ram Mountain, Alberta, Canada. Lower heritability for fitness-related traits, however, was not wholly a consequence of declining genetic variance, because those traits showed high levels of residual variance. Genetic correlations estimated between pairs of traits with significant heritability were positive. Principal component analyses suggest that positive relationships between morphometric traits constitute the main axis of genetic variation. Trade-offs in the form of negative genetic or phenotypic correlations among the traits we have measured do not appear to constrain the potential for evolution in this population.     

Pleiotropic Models of Polygenic Variation, Stabilizing Selection, and Epistasis

We show that in polymorphic populations many polygenic traits pleiotropically related to fitness are expected to be under apparent ``stabilizing selection' independently of the real selection acting on the population. This occurs, for example, if the genetic system is at a stable polymorphic equilibrium determined by selection and the nonadditive contributions of the loci to the trait value either are absent, or are random and independent of those to fitness. Stabilizing selection is also observed if the polygenic system is at an equilibrium determined by a balance between selection and mutation (or migration) when both additive and nonadditive contributions of the loci to the trait value are random and independent of those to fitness. We also compare different viability models that can maintain genetic variability at many loci with respect to their ability to account for the strong stabilizing selection on an additive trait. Let V(m) be the genetic variance supplied by mutation (or migration) each generation, V(g) be the genotypic variance maintained in the population, and n be the number of the loci influencing fitness. We demonstrate that in mutation (migration)-selection balance models the strength of apparent stabilizing selection is order V(m)/V(g). In the overdominant model and in the symmetric viability model the strength of apparent stabilizing selection is approximately 1/(2n) that of total selection on the whole phenotype. We show that a selection system that involves pairwise additive by additive epistasis in maintaining variability can lead to a lower genetic load and genetic variance in fitness (approximately 1/(2n) times) than an equivalent selection system that involves overdominance. We show that, in the epistatic model, the apparent stabilizing selection on an additive trait can be as strong as the total selection on the whole phenotype.     

The changes in genetic and environmental variance with inbreeding in Drosophila melanogaster

We performed a large-scale experiment on the effects of inbreeding and population bottlenecks on the additive genetic and environmental variance for morphological traits in Drosophila melanogaster. Fifty-two inbred lines were created from the progeny of single pairs, and 90 parent-offspring families on average were measured in each of these lines for six wing size and shape traits, as well as 1945 families from the outbred population from which the lines were derived. The amount of additive genetic variance has been observed to increase after such population bottlenecks in other studies; in contrast here the mean change in additive genetic variance was in very good agreement with classical additive theory, decreasing proportionally to the inbreeding coefficient of the lines. The residual, probably environmental, variance increased on average after inbreeding. Both components of variance were highly variable among inbred lines, with increases and decreases recorded for both. The variance among lines in the residual variance provides some evidence for a genetic basis of developmental stability. Changes in the phenotypic variance of these traits are largely due to changes in the genetic variance.     

Quantitative genetics of seed size variation inPhlox

Summary Because seed size is often associated with survival and reproduction in plant populations, genetic variation for seed size may be reduced or eliminated by natural selection. To test this hypothesis we assessed genetic sources of variation in seed size in a population ofPhlox drummondii to determine whether genetic differences among seeds influence the size they attain. A diallel cross among 12 plants from a population at Bastrop, Texas, USA allowed us to partition variance in the mass of seeds among several genetic and parental effects. We found no evidence of additive genetic variance or dominance genetic variance for seed mass in the contribution of plants to their offspring. Extranuclear maternal effects accounted for 56% of the variance in seed mass. A small interaction was observed between seed genotype and maternal plant. Results of this study support theory that predicts little genetic variation for traits associated with fitness.     

The evolution of male genitalia: patterns of genetic variation and covariation in the genital sclerites of the dung beetle Onthophagus taurus

Three main hypotheses, have been invoked to explain divergent genital evolution, the lock and key, pleiotropy, and sexual selection hypotheses, each of which make different predictions about how genital traits are inherited. Here we used a half-sib breeding design to examine the patterns of genetic variation and covariation between male genital sclerites, and their covariance with general body morphology in the dung beetle Onthophagus taurus. We found CV(A)'s and CV(P)'s were similar for both genital and general morphological traits and that CV(R)'s were large for both trait types. We found that male genital sclerites were negatively genetically correlated with general morphological traits. Variation in male genital morphology has direct implications for a male's fertilization success and the resulting sexual selection acting on male genitalia is predicted to maintain high levels of additive genetic variance. Contrary to this prediction, we found that individual genital sclerites all had low levels of additive genetic variance and large maternal and environmental sources of variation. Our data suggest that the genital sclerites in O. taurus are not inherited independently but as a genetically integrated unit. More importantly, the way the different sclerites function to influence male fertilization success reflects this genetic integration. Even though levels of V(A) in individual genital sclerites may be low, there may still be sufficient V(A) in multivariate trait space for selection to generate evolutionary change in the overall morphology of male genitalia.     

The evolution of compensation to herbivory in scarlet gilia, Ipomopsis aggregata: herbivore-imposed natural selection and the quantitative genetics of tolerance

Tolerance is the ability of plants to maintain fitness after experiencing herbivore damage. We investigated scarlet gilia tolerance to browsing in the framework of phenotypic plasticity using both an operational and candidate trait approach. Individuals from full-sib families were split into an artificial clipping treatment, a natural-damage treatment, or left as controls. We tested for genetic variation in tolerance by evaluating family x herbivory treatment interactions on fitness in a mixed model analysis of variance. In addition, we used selection analyses to assess the function of flowering phenology and compensatory regrowth (via branch production) as candidate tolerance traits. We found a strong detrimental fitness effect of browsing and considerable variation among sire half-sib families in levels of tolerance (25% to 63% of the fitness of controls). There was no evidence of overcompensation at either the population or family level and no additive genetic variation in operationally defined tolerance. Phenotypic selection analyses provide evidence that early flowering and compensatory regrowth function as tolerance characters. We found strong linear and correlational selection for early flowering and increased branch production for damaged plants and linear selection for apical dominance (reduced branchiness) and early flowering in control plants. Moreover, reduced phenological delay and increased plasticity in branch production were correlated with tolerance. We detected significant additive genetic variation in flowering phenology in both treatments and a positive genetic correlation between the phenology of control and damaged plants. We found significant additive genetic variation in branch production in undamaged and naturally damaged plants, but not in clipped plants. Damaged plants exhibited marginally significant additive genetic variance in fitness, although its heritability was very low (approximately 3.6%). We failed to find additive genetic variation in the fitness of control plants. Our results suggest that tolerance traits are under herbivore-imposed natural selection in this population, but that responses to selection are limited by available genetic variation and selective constraints.     

Population mixing and the adaptive divergence of quantitative traits in discrete populations: a theoretical framework for empirical tests

Empirical tests for the importance of population mixing in constraining adaptive divergence have not been well grounded in theory for quantitative traits in spatially discrete populations. We develop quantitative-genetic models to examine the equilibrium difference between two populations that are experiencing different selective regimes and exchanging individuals. These models demonstrate that adaptive divergence is negatively correlated with the rate of population mixing (m, most strongly so when m is low), positively correlated with the difference in phenotypic optima between populations, and positively correlated with the amount of additive genetic variance (G, most strongly so when G is low). The approach to equilibrium is quite rapid (fewer than 50 generations for two populations to evolve 90% of the distance to equilibrium) when either heritability or mixing are not too low (h2 > 0.2 or m > 0.05). The theory can be used to aid empirical tests that: (1) compare observed divergence to that predicted using estimates of population mixing, additive genetic variance/covariance, and selection; (2) test for a negative correlation between population mixing and adaptive divergence across multiple independent population pairs; and (3) experimentally manipulate the rate of mixing. Application of the first two of these approaches to data from two well-studied natural systems suggests that population mixing has constrained adaptive divergence for color patterns in Lake Erie water snakes (Nerodia sipedon), but not for trophic traits in sympatric pairs of benthic and limnetic stickleback (Gasterosteus aculeatus). The theoretical framework we outline should provide an improved basis for future empirical tests of the role of population mixing in adaptive divergence.     

Natural selection on floral traits of female Silene dioica by a sexually transmitted disease

Floral traits endowing high reproductive fitness can also affect the probability of plants contracting sexually transmitted diseases. We explore how variations in floral traits influence the fitness of Silene dioica females in their interactions with pollinators carrying pollen or spores of the sterilizing anther-smut fungus Microbotryum violaceum. We collected healthy and infected plants in a highly diseased population and grew them under conditions that 'cure' infected individuals, and used standard regression methods to detect natural selection on floral traits. Narrow-sense heritabilities, coefficients of additive genetic variation (CV(A)) and genetic correlations among traits were estimated from paternal half-sib groups. Pollinator preferences imposed strong direct and directional selection on traits affecting female attractiveness and pollen-/spore-capturing abilities. Levels of additive genetic variance were high in these traits, suggesting that rapid responses to selection are possible. By considering our results in the light of spatial and temporal heterogeneity resulting from the colonization dynamics typical for this species, we suggest that the conflicting selective effects of pollen/spore loads lead to the maintenance of genetic variation in these traits.     

Population structure and genetic diversity in Swainson’s Hawks (Buteo swainsoni): implications for conservation

Swainson’s Hawks (Buteo swainsoni) are large raptors with a breeding distribution extending across much of western North America where they were historically considered one of the most abundant raptors. Swainson’s Hawks have declined precipitously in many parts of their range during the 20th century, and the historical range in California has been much reduced. In the Central Valley of California (CV), Swainson’s Hawks exhibit behavioral and morphological characteristics apparently different from other regions. To describe the genetic diversity and population structure of Swainson’s Hawks throughout their range, 19 microsatellite loci and 416 base pairs of the mitochondrial control region were analyzed. Microsatellite diversity appears high throughout the contemporary range. A Bayesian model-based analysis of microsatellite genotypes revealed clusters associated with the CV and the Great Basin/Great Plains region of North America (GBGP) with overlap between regions. F _ST estimates suggest limited differentiation among Swainson’s Hawks with isolation by distance. A heterozygote excess indicated a recent reduction in effective population size of Swainson’s Hawks across all regions. Control region data revealed no population structure and provided evidence of historic population expansion in the GBGP. In the CV a weaker signature of population expansion was detected, possibly altered by recent declines. While genetic data suggests recent gene-flow across regions, apparent differences between the CV and GBGP in traits with potential fitness consequences (migratory behavior and morphology) along with marked decline in numbers in California call for careful conservation, management, and monitoring of Swainson’s Hawks in the CV.