Effects of Competition and Temporal Variation on the Evolutionary Potential of Two Native Bunchgrass Species

The capacity of restored plant populations to adapt to new environmental challenges depends on within‐population genetic variation. We examined how much genetic and environmentally based variation for fitness‐associated traits exists within populations of two native grasses commonly used for restoration in California. We were also interested in understanding how phenotypic expression of genetic variation for these traits varies with growth environment. Thirty maternal families of Elymus glaucus (Blue wild rye) and Nassella pulchra (Purple needlegrass) were sampled from both coastal and interior populations and reciprocally transplanted into three replicated common gardens with and without interspecific competition at each site. Reproductive output of families differed both among years and with competition treatments. Phenotypic expression of genetic variation in culm production differed among populations and was very low when families were grown with interspecific competition. Without interspecific competition, the degree of genetic determination peaked in year two in both species (8.4 and 15.1% in E. glaucus and N. pulchra, respectively). Significant genetic differences in reproduction and phenotypic plasticity were found among N. pulchra subpopulations sampled less than 3 km apart, further highlighting the importance of thoroughly sampling available genetic variation in populations used for restoration. The variable and generally low expression of genetic variation indicates that rates of adaptation in restored populations of these native grasses may vary temporally and may be especially slow within competitive environments.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

Variation in chronic radiation exposure does not drive life history divergence among Daphnia populations across the Chernobyl Exclusion Zone

Ionizing radiation is a mutagen with known negative impacts on individual fitness. However, much less is known about how these individual fitness effects translate into population‐level variation in natural environments that have experienced varying levels of radiation exposure. In this study, we sampled genotypes of the freshwater crustacean, Daphnia pulex, from the eight inhabited lakes across the Chernobyl Exclusion Zone (CEZ). Each lake has experienced very different levels of chronic radiation exposure since a nuclear power reactor exploded there over thirty years ago. The sampled Daphnia genotypes represent genetic snapshots of current populations and allowed us to examine fitness‐related traits under controlled laboratory conditions at UK background dose rates. We found that whilst there was variation in survival and schedules of reproduction among populations, there was no compelling evidence that this was driven by variation in exposure to radiation. Previous studies have shown that controlled exposure to radiation at dose rates included in the range measured in the current study reduce survival, or fecundity, or both. One limitation of this study is the lack of available sites at high dose rates, and future work could test life history variation in various organisms at other high radiation areas. Our results are nevertheless consistent with the idea that other ecological factors, for example competition, predation or parasitism, are likely to play a much bigger role in driving variation among populations than exposure to the high radiation dose rates found in the CEZ. These findings clearly demonstrate that it is important to examine the potential negative effects of radiation across wild populations that are subject to many and varied selection pressures as a result of complex ecological interactions.     

Natural quantitative genetic variance in plant growth differs in response to ecologically relevant temperature heterogeneity

Adaptation to large‐scale spatial heterogeneity in the environment accounts for a major proportion of genetic diversity within species. Theory predicts the erosion of adaptive genetic variation on a within‐population level, but considerable genetic diversity is often found locally. Genetic diversity could be expected to be maintained within populations in temporally or spatially variable conditions if genotypic rank orders vary across contrasting microenvironmental settings. Taking advantage of fine‐resolution environmental data, we tested the hypothesis that temperature heterogeneity among years could be one factor maintaining quantitative genetic diversity within a natural and genetically diverse plant population. We sampled maternal families of Boechera stricta, an Arabidopsis thaliana relative, at one location in the central Rocky Mountains and grew them in three treatments that, based on records from an adjacent weather station, simulated hourly temperature changes at the native site during three summers with differing mean temperatures. Treatment had a significant effect on all traits, with 2–3‐fold increase in above‐ and belowground biomass and the highest allocation to roots observed in the treatment simulating the warmest summer on record at the site. Treatment affected bivariate associations between traits, with the weakest correlation between above‐ and belowground biomass in the warmest treatment. The magnitude of quantitative genetic variation for all traits differed across treatments: Genetic variance of biomass was 0 in the warmest treatment, while highly significant diversity was found in average conditions, resulting in broad‐sense heritability of 0.31. Significant genotype × environment interactions across all treatments were found only in root‐to‐shoot ratio. Therefore, temperature variation among summers appears unlikely to account for the observed levels of local genetic variation in size in this perennial species, but may influence family rank order in growth allocation. Our results indicate that natural environmental fluctuations can have a large impact on the magnitude of within‐population quantitative genetic variance.     

Phenotypic selection on growth rhythm in whitebark pine under climatic conditions warmer than seed origins

Growth rhythm that is well synchronized with seasonal changes in local climatic conditions is understood to enhance fitness; however, rapid ongoing climate change threatens to disrupt this synchrony. To evaluate phenotypic selection on growth rhythm under expected warmer and drier future climate, seedlings from 49 populations of whitebark pine (Pinus albicaulis Engelm.) were grown and measured over more than 10 years in two common garden field experiments on sites that approximate the projected future climate of the seed origins. Selection on growth rhythm was assessed by relating individual plant fitness to timing and rate of shoot elongation. Differential survival clearly evidenced selection on growth rhythm. We detected directional and stabilizing selection that varied in magnitude between experimental sites and among years. The observed phenotypic selection supports the interpretation of clinal variation among populations within tree species as reflecting adaptive variation in response to past natural selection mediated by climate. To the extent that growth rhythm is heritable, results of the present study suggest evolution of whitebark pine toward a more distinct timing of shoot elongation and generally more rapid elongation in the immediate next generation under ongoing climate change in environments similar to the study sites.     

Genetic variation for sexual dimorphism in flower size within and between populations of gynodioecious Thymus vulgaris

There has been very little empirical study of quantitative genetic variation in flower size in sexually dimorphic plant species, despite the frequent occurrence of flower size differences between sexual phenotypes. In this study we quantify the nature of quantitative flower size variation in females and hermaphrodites of gynodioecious Thymus vulgaris. In a field study, females had significantly smaller flowers than hermaphrodites, and the degree of flower size dimorphism varied significantly among populations. To quantify the genetic basis of flower size variation we sampled maternal progeny from 10 F₀ females in three populations (across the range of variation in flower size in the field), performed controlled crosses on F₁ offspring in the glasshouse and grew F₂ progeny to flowering in uniform field conditions. A significant population * sex interaction was again observed, hence the degree of sexual dimorphism shows genetic variation among populations. A significant family * sex interaction was also observed, indicating that the degree of sexual dimorphism shows genetic variation among families. Females showed significantly greater variation among populations and among families than hermaphrodites. Female flower size varied significantly depending on the degree of stamen abortion, with morphologically intermediate females having flowers more similar to hermaphrodites than to other females. The frequency of female types that differ in the degree of stamen abortion varied among populations and families and mean family female flower size increased as the proportion of intermediate female types increased across families. Variation in the degree of flower size dimorphism thus appears to be a result of variation in the degree of stamen abortion in females, the potential causes of which are discussed.     

Patterns of genetic diversity vary among shoot and root functional traits in Norway spruce Picea abies along a latitudinal gradient

Roots constitute a major segment of plant biomass, and variation in belowground traits in situ correlates with environmental gradients at large spatial scales. Local adaptation of populations maintains intraspecific genetic variation in various shoot traits, but the contribution of genetic factors to adaptation to soil heterogeneity remains poorly known. I established a common-garden experiment with three Norway spruce Picea abies populations sampled between 60° and 67° N in Finland, each represented by 13 or 15 maternal families, to determine whether belowground traits are as genetically differentiated among populations as those in the shoot along a collective latitudinal gradient of temperature and soil heterogeneity. Two growing season simulations enabled testing for among-population differences in phenotypic plasticity. I phenotyped 777 first-year seedlings from shoot to root to capture functional traits that may influence survival in the wild: autumn phenology, shoot growth, root system size, root architecture, root morphology and growth allocation. All traits exhibited within-population genetic diversity, but among-population differentiation ranged from strong in shoot traits to nonexistent in root system architecture and morphology that are scaled to root system size. However, latitudinal trends characterised root-to-shoot ratio and root tip-to-shoot ratio that account for among-population differences in aboveground growth. Overall trait variability was multidimensional with variable among- versus within-population trends: for example, phenology and shoot growth covaried across populations, but their association within individual populations was variable. Shoot growth correlated positively with root system size, but not with root architecture or morphology. Finally, the two higher-latitude populations exhibited greater phenotypic plasticity in shoot traits and growth allocation. The results demonstrate varying patterns of genetic variation in functional traits of Norway spruce in the boreal zone, suggesting simultaneous adaptation to multiple environmental factors. Functional traits that exhibit phenotypic plasticity, genetic diversity and little covariation will promote long-term survival of populations in fluctuating environments.     

Quantifying natural seasonal variation in mutation parameters with mutation accumulation lines

Mutations create novel genetic variants, but their contribution to variation in fitness and other phenotypes may depend on environmental conditions. Furthermore, natural environments may be highly heterogeneous. We assessed phenotypes associated with survival and reproductive success in over 30,000 plants representing 100 mutation accumulation lines of Arabidopsis thaliana across four temporal environments at a single field site. In each of the four assays, environmental variance was substantially larger than mutational variance. For some traits, whether mutational variance was significantly varied between seasons. The founder genotype had mean trait values near the mean of the distribution of the mutation accumulation lines in all field experiments. New mutations also contributed more phenotypic variation than would be predicted, given phenotypic and sequence‐level divergence among natural populations of A. thaliana. The combination of large environmental variance with a mean effect of mutation near zero suggests that mutations could contribute substantially to standing genetic variation.     

Assortative mating by flowering time and its effect on correlated traits in variable environments

Reproductive timing is a key life‐history trait that impacts the pool of available mates, the environment experienced during flowering, and the expression of other traits through genetic covariation. Selection on phenology, and its consequences on other life‐history traits, has considerable implications in the context of ongoing climate change and shifting growing seasons. To test this, we grew field‐collected seed from the wildflower Mimulus guttatus in a greenhouse to assess the standing genetic variation for flowering time and covariation with other traits. We then created full‐sib families through phenological assortative mating and grew offspring in three photoperiod treatments representing seasonal variation in daylength. We find substantial quantitative genetic variation for the onset of flowering time, which covaried with vegetative traits. The assortatively‐mated offspring varied in their critical photoperiod by over two hours, so that families differed in their probability of flowering across treatments Allocation to flowering and vegetative growth changed across the daylength treatments, with consistent direction and magnitude of covariation among flowering time and other traits. Our results suggest that future studies of flowering time evolution should consider the joint evolution of correlated traits and shifting seasonal selection to understand how environmental variation influences life histories.     

Testing for local adaptation and evolutionary potential along altitudinal gradients in rainforest Drosophila: beyond laboratory estimates

Predicting how species will respond to the rapid climatic changes predicted this century is an urgent task. Species distribution models (SDMs) use the current relationship between environmental variation and species’ abundances to predict the effect of future environmental change on their distributions. However, two common assumptions of SDMs are likely to be violated in many cases: (i) that the relationship of environment with abundance or fitness is constant throughout a species’ range and will remain so in future and (ii) that abiotic factors (e.g. temperature, humidity) determine species’ distributions. We test these assumptions by relating field abundance of the rainforest fruit fly Drosophila birchii to ecological change across gradients that include its low and high altitudinal limits. We then test how such ecological variation affects the fitness of 35 D. birchii families transplanted in 591 cages to sites along two altitudinal gradients, to determine whether genetic variation in fitness responses could facilitate future adaptation to environmental change. Overall, field abundance was highest at cooler, high‐altitude sites, and declined towards warmer, low‐altitude sites. By contrast, cage fitness (productivity) increased towards warmer, lower‐altitude sites, suggesting that biotic interactions (absent from cages) drive ecological limits at warmer margins. In addition, the relationship between environmental variation and abundance varied significantly among gradients, indicating divergence in ecological niche across the species’ range. However, there was no evidence for local adaptation within gradients, despite greater productivity of high‐altitude than low‐altitude populations when families were reared under laboratory conditions. Families also responded similarly to transplantation along gradients, providing no evidence for fitness trade‐offs that would favour local adaptation. These findings highlight the importance of (i) measuring genetic variation in key traits under ecologically relevant conditions, and (ii) considering the effect of biotic interactions when predicting species’ responses to environmental change.     

Explaining the sawtooth: latitudinal periodicity in a circadian gene correlates with shifts in generation number

Many temperate insects take advantage of longer growing seasons at lower latitudes by increasing their generation number or voltinism. In some insects, development time abruptly decreases when additional generations are fit into the season. Consequently, latitudinal ‘sawtooth’ clines associated with shifts in voltinism are seen for phenotypes correlated with development time, like body size. However, latitudinal variation in voltinism has not been linked to genetic variation at specific loci. Here, we show a pattern in allele frequency among voltinism ecotypes of the European corn borer moth (Ostrinia nubilalis) that is reminiscent of a sawtooth cline. We characterized 145 autosomal and sex‐linked SNPs and found that period, a circadian gene that is genetically linked to a major QTL determining variation in post‐diapause development time, shows cyclical variation between voltinism ecotypes. Allele frequencies at an unlinked circadian clock gene cryptochrome1 were correlated with period. These results suggest that selection on development time to ‘fit’ complete life cycles into a latitudinally varying growing season produces oscillations in alleles associated with voltinism, primarily through changes at loci underlying the duration of transitions between diapause and other life history phases. Correlations among clock loci suggest possible coupling between the circadian clock and the circannual rhythms for synchronizing seasonal life history. We anticipate that latitudinal oscillations in allele frequency will represent signatures of adaptation to seasonal environments in other insects and may be critical to understanding the ecological and evolutionary consequences of variable environments, including response to global climate change.     

Persistence of transgenes in wild rice populations depends on the interaction between genetic background of recipients and environmental conditions

The persistence of transgenes in wild populations may cause unintended ecological consequences, and the possibility of transgenes' persistence and introgression is dependent on fitness performance of transgenic crop–wild hybrids. To investigate the effects of transgene and genotype × environment on the fitness of crop–wild rice hybrids, a total of 11 cross‐combination progenies between insect‐resistant transgene (CpTI and Bt/CpTI) rice lines and wild rice (Oryza rufipogon) were evaluated at different sites with contrasting insect treatments. The results showed that fitness performance varied between transgenic hybrids having different wild parents and under different environmental conditions, indicating that fitness effects of transgenes on hybrid progenies depend heavily on the genetic background of recipient plants and growing environment. Significant fitness advantages conferred by transgenes were found only in some hybrids under high insect pressure condition, demonstrating that the level of target insects in the field environment influences the persistence and spread of insect‐resistant transgenes in wild rice populations. These findings suggest that evolutionary fate of escaped transgenes is different in wild populations with diverse genetic backgrounds under various environmental conditions.     

Aphid genotypes vary in their response to the presence of fungal endosymbionts in host plants

Genetic variation for fitness‐relevant traits may be maintained in natural populations by fitness differences that depend on environmental conditions. For herbivores, plant quality and variation in chemical plant defences can maintain genetic variation in performance. Apart from plant secondary compounds, symbiosis between plants and endosymbiotic fungi (endophytes) can produce herbivore‐toxic compounds. We show that there is significant variation among aphid genotypes in response to endophytes by comparing life‐history traits of 37 clones of the bird cherry‐oat aphid Rhopalosiphum padi feeding on endophyte‐free and endophyte‐infected tall fescue Lolium arundinaceum. Clonal variation for life‐history traits was large, and most clones performed better on endophyte‐free plants. However, the clones differed in the relative performance across the two environments, resulting in significant genotype × environment interactions for all reproductive traits. These findings suggest that natural variation in prevalence of endophyte infection can contribute to the maintenance of genetic diversity in aphid populations.     

Quantitative genetic variance in experimental fly populations evolving with or without environmental heterogeneity

Heterogeneous environments are typically expected to maintain more genetic variation in fitness within populations than homogeneous environments. However, the accuracy of this claim depends on the form of heterogeneity as well as the genetic basis of fitness traits and how similar the assay environment is to the environment of past selection. Here, we measure quantitative genetic (QG) variance for three traits important for fitness using replicated experimental populations of Drosophila melanogaster evolving under four selective regimes: constant salt‐enriched medium (Salt), constant cadmium‐enriched medium (Cad), and two heterogeneous regimes that vary either temporally (Temp) or spatially (Spatial). As theory predicts, we found that Spatial populations tend to harbor more genetic variation than Temp populations or those maintained in a constant environment that is the same as the assay environment. Contrary to expectation, Salt populations tend to have more genetic variation than Cad populations in both assay environments. We discuss the patterns for QG variances across regimes in relation to previously reported data on genome‐wide sequence diversity. For some traits, the QG patterns are similar to the diversity patterns of ecological selected SNPs, whereas the QG patterns for some other traits resembled that of neutral SNPs.     

Adaptive phenotypic divergence in an annual grass differs across biotic contexts*

Climate is a powerful force shaping adaptation within species, yet adaptation to climate occurs against a biotic background: species interactions can filter fitness consequences of genetic variation by altering phenotypic expression of genotypes. We investigated this process using populations of teosinte, a wild annual grass related to maize (Zea mays ssp. mexicana), sampling plants from 10 sites along an elevational gradient as well as rhizosphere biota from three of those sites. We grew half‐sibling teosinte families in each biota to test whether trait divergence among teosinte populations reflects adaptation or drift, and whether rhizosphere biota affect expression of diverged traits. We further assayed the influence of rhizosphere biota on contemporary additive genetic variation. We found that adaptation across environment shaped divergence of some traits, particularly flowering time and root biomass. We also observed that different rhizosphere biota shifted expressed values of these traits within teosinte populations and families and altered within‐population genetic variance and covariance. In sum, our results imply that changes in trait expression and covariance elicited by rhizosphere communities could have played a historical role in teosinte adaptation to environments and that they are likely to play a role in the response to future selection.     

Circadian rhythms vary over the growing season and correlate with fitness components

Circadian clocks have evolved independently in all three domains of life, suggesting that internal mechanisms of time‐keeping are adaptive in contemporary populations. However, the performance consequences of either discrete or quantitative clock variation have rarely been tested in field settings. Clock sensitivity of diverse segregating lines to the environment remains uncharacterized as do the statistical genetic parameters that determine evolutionary potential. In field studies with Arabidopsis thaliana, we found that major perturbations to circadian cycle length (referred to as clock period) via mutation reduce both survival and fecundity. Subtler adjustments via genomic introgression of naturally occurring alleles indicated that clock periods slightly >24 hr were adaptive, consistent with prior models describing how well the timing of biological processes is adjusted within a diurnal cycle (referred to as phase). In segregating recombinant inbred lines (RILs), circadian phase varied up to 2 hr across months of the growing season, and both period and phase expressed significant genetic variances. Performance metrics including developmental rate, size and fruit set were described by principal components (PC) analyses and circadian parameters correlated with the first PC, such that period lengths slightly >24 hr were associated with improved performance in multiple RIL sets. These experiments translate functional analyses of clock behaviour performed in controlled settings to natural ones, demonstrating that quantitative variation in circadian phase is highly responsive to seasonally variable abiotic factors. The results expand upon prior studies in controlled settings, showing that discrete and quantitative variation in clock phenotypes correlates with performance in nature.     

Climate change shifts natural selection and the adaptive potential of the perennial forb Boechera stricta in the Rocky Mountains

Heritable genetic variation is necessary for populations to evolve in response to anthropogenic climate change. However, antagonistic genetic correlations among traits may constrain the rate of adaptation, even if substantial genetic variation exists. We examine potential genetic responses to selection by comparing multivariate genetic variance–covariances of traits and fitness (multivariate Robertson–Price identities) across different environments in a reciprocal transplant experiment of the forb Boechera stricta in the Rocky Mountains. By transplanting populations into four common gardens arrayed along an elevational gradient, and exposing populations to control and snow removal treatments, we simulated future and current climates and snowmelt regimes. Genetic variation in flowering and germination phenology declined in plants moved downslope to warmer, drier sites, suggesting that these traits may have a limited ability to evolve under future climates. Simulated climate change via snow removal altered the strength of selection on flowering traits, but we found little evidence that genetic correlations among traits are likely to affect the rate of adaptation to climate change. Overall, our results suggest that climate change may alter the evolutionary potential of B. stricta, but reduced expression of genetic variation may be a larger impediment to adaptation than constraints imposed by antagonistic genetic correlations.     

CHARACTERIZING SELECTION ON PHENOTYPIC PLASTICITY IN RESPONSE TO NATURAL ENVIRONMENTAL HETEROGENEITY

Adaptive genetic differentiation and adaptive phenotypic plasticity can increase the fitness of plant lineages in heterogeneous environments. We examine the relative importance of genetic differentiation and plasticity in determining the fitness of the annual plant, Erodium cicutarium, in a serpentine grassland in California. Previous work demonstrated that the serpentine sites within this mosaic display stronger dispersal‐scale heterogeneity than nonserpentine sites. We conducted a reciprocal transplant experiment among six sites to characterize selection on plasticity expressed by 180 full‐sibling families in response to natural environmental heterogeneity across these sites. Multivariate axes of environmental variation were constructed using a principal components analysis of soil chemistry data collected at every experimental block. Simple linear regressions were used to characterize the intercept, and slope (linear and curvilinear) of reaction norms for each full‐sibling family in response to each axis of environmental variation. Multiple linear regression analyses revealed significant selection on trait means and slopes of reaction norms. Multivariate analyses of variance demonstrated genetic differentiation between serpentine and nonserpentine lineages in the expression of plasticity in response to three of the five axes of environmental variation considered. In all but one case, serpentine genotypes expressed a stronger adaptive plastic response than nonserpentine genotypes.     

Postglacial colonization history reflects in the genetic structure of natural populations of Festuca rubra in Europe

We conducted a large‐scale population genetic survey of genetic diversity of the host grass Festuca rubra s.l., which fitness can be highly dependent on its symbiotic fungus Epichloë festucae, to evaluate genetic variation and population structure across the European range. The 27 studied populations have previously been found to differ in frequencies of occurrence of the symbiotic fungus E. festucae and ploidy levels. As predicted, we found decreased genetic diversity in previously glaciated areas in comparison with nonglaciated regions and discovered three major maternal genetic groups: southern, northeastern, and northwestern Europe. Interestingly, host populations from Greenland were genetically similar to those from the Faroe Islands and Iceland, suggesting gene flow also between those areas. The level of variation among populations within regions is evidently highly dependent on the postglacial colonization history, in particular on the number of independent long‐distance seed colonization events. Yet, also anthropogenic effects may have affected the population structure in F. rubra. We did not observe higher fungal infection rates in grass populations with lower levels of genetic variability. In fact, the fungal infection rates of E. festucae in relation to genetic variability of the host populations varied widely among geographical areas, which indicate differences in population histories due to colonization events and possible costs of systemic fungi in harsh environmental conditions. We found that the plants of different ploidy levels are genetically closely related within geographic areas indicating independent formation of polyploids in different maternal lineages.     

Heterosis is common and inbreeding depression absent in natural populations of Arabidopsis thaliana

The importance of genetic drift in shaping patterns of adaptive genetic variation in nature is poorly known. Genetic drift should drive partially recessive deleterious mutations to high frequency, and inter‐population crosses may therefore exhibit heterosis (increased fitness relative to intra‐population crosses). Low genetic diversity and greater genetic distance between populations should increase the magnitude of heterosis. Moreover, drift and selection should remove strongly deleterious recessive alleles from individual populations, resulting in reduced inbreeding depression. To estimate heterosis, we crossed 90 independent line pairs of Arabidopsis thaliana from 15 pairs of natural populations sampled across Fennoscandia and crossed an additional 41 line pairs from a subset of four of these populations to estimate inbreeding depression. We measured lifetime fitness of crosses relative to parents in a large outdoor common garden (8,448 plants in total) in central Sweden. To examine the effects of genetic diversity and genetic distance on heterosis, we genotyped parental lines for 869 SNPs. Overall, genetic variation within populations was low (median expected heterozygosity = 0.02), and genetic differentiation was high (median F_ST = 0.82). Crosses between 10 of 15 population pairs exhibited significant heterosis, with magnitudes of heterosis as high as 117%. We found no significant inbreeding depression, suggesting that the observed heterosis is due to fixation of mildly deleterious alleles within populations. Widespread and substantial heterosis indicates an important role for drift in shaping genetic variation, but there was no significant relationship between fitness of crosses relative to parents and genetic diversity or genetic distance between populations.