THE EVOLVABILITY OF GROWTH FORM IN A CLONAL SEAWEED

Although modular construction is considered the key to adaptive growth or growth‐form plasticity in sessile taxa (e.g., plants, seaweeds and colonial invertebrates), the serial expression of genes in morphogenesis may compromise its evolutionary potential if growth forms emerge as integrated wholes from module iteration. To explore the evolvability of growth form in the red seaweed, Asparagopsis armata, we estimated genetic variances, covariances, and cross‐environment correlations for principal components of growth‐form variation in contrasting light environments. We compared variance–covariance matrices across environments to test environmental effects on heritable variation and examined the potential for evolutionary change in the direction of plastic responses to light. Our results suggest that growth form in Asparagopsis may constitute only a single genetic entity whose plasticity affords only limited evolutionary potential. We argue that morphological integration arising from modular construction may constrain the evolvability of growth form in Asparagopsis, emphasizing the critical distinction between genetic and morphological modularity in this and other modular taxa.     

A multivariate test of evolutionary constraints for thermal tolerance in Drosophila melanogaster

Exposure to extreme temperatures is increasingly likely to impose strong selection on many organisms in their natural environments. The ability of organisms to adapt to such selective pressures will be determined by patterns of genetic variation and covariation. Despite increasing interest in thermal adaptation, few studies have examined the extent to which the genetic covariance between traits might constrain thermal responses. Furthermore, it remains unknown whether sex‐specific genetic architectures will constrain responses to climatic selection. We used a paternal half‐sibling breeding design to examine whether sex‐specific genetic architectures and genetic covariances between traits might constrain evolutionary responses to warming climates in a population of Drosophila melanogaster. Our results suggest that the sexes share a common genetic underpinning for heat tolerance as indicated by a strong positive inter‐sexual genetic correlation. Further, we found no evidence in either of the sexes that genetic trade‐offs between heat tolerance and fitness will constrain responses to thermal selection. Our results suggest that neither trade‐offs, nor sex‐specific genetics, will significantly constrain an evolutionary response to climatic warming, at least in this population of D. melanogaster.     

Evolutionary potential of a widespread clonal grass under changing climate

Adaptive responses are probably the most effective long‐term responses of populations to climate change, but they require sufficient evolutionary potential upon which selection can act. This requires high genetic variance for the traits under selection and low antagonizing genetic covariances between the different traits. Evolutionary potential estimates are still scarce for long‐lived, clonal plants, although these species are predicted to dominate the landscape with climate change. We studied the evolutionary potential of a perennial grass, Festuca rubra, in western Norway, in two controlled environments corresponding to extreme environments in natural populations: cold–dry and warm–wet, the latter being consistent with the climatic predictions for the country. We estimated genetic variances, covariances, selection gradients and response to selection for a wide range of growth, resource acquisition and physiological traits, and compared their estimates between the environments. We showed that the evolutionary potential of F. rubra is high in both environments, and genetic covariances define one main direction along which selection can act with relatively few constraints to selection. The observed response to selection at present is not sufficient to produce genotypes adapted to the predicted climate change under a simple, space for time substitution model. However, the current populations contain genotypes which are pre‐adapted to the new climate, especially for growth and resource acquisition traits. Overall, these results suggest that the present populations of the long‐lived clonal plant may have sufficient evolutionary potential to withstand long‐term climate changes through adaptive responses.     

QUANTITATIVE GENETIC DIVERGENCE AND STANDING GENETIC (CO)VARIANCE IN THERMAL REACTION NORMS ALONG LATITUDE

Although the potential to adapt to warmer climate is constrained by genetic trade‐offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade‐offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region‐specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade‐offs in natural environments.     

THE EVOLUTIONARY STABILITY OF CROSS‐SEX,CROSS‐TRAIT GENETIC COVARIANCES

Although knowledge of the selective agents behind the evolution of sexual dimorphism has advanced considerably in recent years, we still lack a clear understanding of the evolutionary durability of cross‐sex genetic covariances that often constrain its evolution. We tested the relative stability of cross‐sex genetic covariances for a suite of homologous contact pheromones of the fruit fly Drosophila serrata, along a latitudinal gradient where these traits have diverged in mean. Using a Bayesian framework, which allowed us to account for uncertainty in all parameter estimates, we compared divergence in the total amount and orientation of genetic variance across populations, finding divergence in orientation but not total variance. We then statistically compared orientation divergence of within‐sex ( G ) to cross‐sex ( B ) covariance matrices. In line with a previous theoretical prediction, we find that the cross‐sex covariance matrix, B , is more variable than either within‐sex G matrix. Decomposition of B matrices into their symmetrical and nonsymmetrical components revealed that instability is linked to the degree of asymmetry. We also find that the degree of asymmetry correlates with latitude suggesting a role for spatially varying natural selection in shaping genetic constraints on the evolution of sexual dimorphism.     

Constraints on the evolution of function‐valued traits: a study of growth in Tribolium castaneum

Growth trajectories often impact individual fitness. They are continuous by nature and so are amenable to analysis using a function‐valued (FV) trait framework to reveal their underlying genetic architecture. Previous studies have found high levels of standing additive genetic (co)variance for growth trajectories despite the expectation that growth should be responding to frequent strong directional selection. In this study, the FV framework is used to estimate the additive genetic covariance function for growth trajectories in larval Tribolium castaneum to address questions about standing additive genetic (co)variance and possible evolutionary constraints on growth and to predict responses to four plausible selection regimes. Results show that additive genetic (co)variance is high at the early ages, but decreases towards later ages in the larval period. A selection gradient function of the same size and in the same direction of the first eigenfunction of the G‐function should give the maximal response. However, evolutionary constraints may be acting to keep this maximal response from being realized, through either conflicting effects on survivability and fecundity of larger body size, few evolutionary directions having sufficient additive variance for a response, genetic trade‐offs with other traits or physiological regulatory mechanisms. More light may be shed on these constraints through the development of more sophisticated statistical approaches and implementation of additional empirical studies to explicitly test for specific types of constraints.     

ANTAGONISTIC VERSUS NONANTAGONISTIC MODELS OF BALANCING SELECTION: CHARACTERIZING THE RELATIVE TIMESCALES AND HITCHHIKING EFFECTS OF PARTIAL SELECTIVE SWEEPS

Antagonistically selected alleles‐–those with opposing fitness effects between sexes, environments, or fitness components‐–represent an important component of additive genetic variance in fitness‐related traits, with stably balanced polymorphisms often hypothesized to contribute to observed quantitative genetic variation. Balancing selection hypotheses imply that intermediate‐frequency alleles disproportionately contribute to genetic variance of life‐history traits and fitness. Such alleles may also associate with population genetic footprints of recent selection, including reduced genetic diversity and inflated linkage disequilibrium at linked, neutral sites. Here, we compare the evolutionary dynamics of different balancing selection models, and characterize the evolutionary timescale and hitchhiking effects of partial selective sweeps generated under antagonistic versus nonantagonistic (e.g., overdominant and frequency‐dependent selection) processes. We show that the evolutionary timescales of partial sweeps tend to be much longer, and hitchhiking effects are drastically weaker, under scenarios of antagonistic selection. These results predict an interesting mismatch between molecular population genetic and quantitative genetic patterns of variation. Balanced, antagonistically selected alleles are expected to contribute more to additive genetic variance for fitness than alleles maintained by classic, nonantagonistic mechanisms. Nevertheless, classical mechanisms of balancing selection are much more likely to generate strong population genetic signatures of recent balancing selection.     

A MULTIVARIATE ANALYSIS OF GENETIC VARIATION IN THE ADVERTISEMENT CALL OF THE GRAY TREEFROG,HYLA VERSICOLOR

Genetic variation in sexual displays is crucial for an evolutionary response to sexual selection, but can be eroded by strong selection. Identifying the magnitude and sources of additive genetic variance underlying sexually selected traits is thus an important issue in evolutionary biology. We conducted a quantitative genetics experiment with gray treefrogs (Hyla versicolor) to investigate genetic variances and covariances among features of the male advertisement call. Two energetically expensive traits showed significant genetic variation: call duration, expressed as number of pulses per call, and call rate, represented by its inverse, call period. These two properties also showed significant genetic covariance, consistent with an energetic constraint to call production. Combining the genetic variance–covariance matrix with previous estimates of directional sexual selection imposed by female preferences predicts a limited increase in call duration but no change in call rate despite significant selection on both traits. In addition to constraints imposed by the genetic covariance structure, an evolutionary response to sexual selection may also be limited by high energetic costs of long‐duration calls and by preferences that act most strongly against very short‐duration calls. Meanwhile, the persistence of these preferences could be explained by costs of mating with males with especially unattractive calls.     

Genotype‐by‐environment interactions underlie the expression of pre‐ and post‐copulatory sexually selected traits in guppies

The role that genotype‐by‐environment interactions (GEIs) play in sexual selection has only recently attracted the attention of evolutionary biologists. Yet GEIs can have profound evolutionary implications by compromising the honesty of sexual signals, maintaining high levels of genetic variance underlying their expression and altering the patterns of genetic covariance among fitness traits. In this study, we test for GEIs in a highly sexually dimorphic freshwater fish, the guppy Poecilia reticulata. We conducted an experimental quantitative genetic study in which male offspring arising from a paternal half‐sibling breeding design were assigned to differing nutritional ‘environments’ (either high or low feed levels). We then determined whether the manipulation of diet quantity influenced levels of additive genetic variance and covariance for several highly variable and condition‐dependent pre‐ and post‐copulatory sexual traits. In accordance with previous work, we found that dietary limitation had strong phenotypic effects on numerous pre‐ and post‐copulatory sexual traits. We also report evidence for significant GEI for several of these traits, which in some cases (area of iridescence and sperm velocity) reflected a change in the rank order of genotypes across different nutritional environments (i.e. ecological crossover). Furthermore, we show that genetic correlations vary significantly between nutritional environments. Notably, a highly significant negative genetic correlation between iridescent coloration and sperm viability in the high food treatment broke down under dietary restriction. Taken together, these findings are likely to have important evolutionary implications for guppies; ecological crossover may influence sexual signal reliability in unstable (nutritional) environments and contribute towards the extreme levels of polymorphism in sexual traits typically reported for this species. Furthermore, the presence of environment‐specific genetic covariance suggests that trade‐offs measured in one environment may not be indicative of genetic constraints in others.     

TEMPORAL VARIATION FAVORS THE EVOLUTION OF GENERALISTS IN EXPERIMENTAL POPULATIONS OF DROSOPHILA MELANOGASTER

In variable environments, selection should favor generalists that maintain fitness across a range of conditions. However, costs of adaptation may generate fitness trade‐offs and lead to some compromise between specialization and generalization that maximizes fitness. Here, we evaluate the evolution of specialization and generalization in 20 populations of Drosophila melanogaster experimentally evolved in constant and variable thermal environments for 3 years. We developed genotypes from each population at two temperatures after which we measured fecundity across eight temperatures. We predicted that constant environments would select for thermal specialists and that variable environments would select for thermal generalists. Contrary to our predictions, specialists and generalists did not evolve in constant and spatially variable environments, respectively. However, temporal variation produced a type of generalist that has rarely been considered by theoretical models of developmental plasticity. Specifically, genotypes from the temporally variable selective environment were more fecund across all temperatures than were genotypes from other environments. These patterns suggest certain allelic effects and should inspire new directions for modeling adaptation to fluctuating environments.     

Very low levels of direct additive genetic variance in fitness and fitness components in a red squirrel population

A trait must genetically correlate with fitness in order to evolve in response to natural selection, but theory suggests that strong directional selection should erode additive genetic variance in fitness and limit future evolutionary potential. Balancing selection has been proposed as a mechanism that could maintain genetic variance if fitness components trade off with one another and has been invoked to account for empirical observations of higher levels of additive genetic variance in fitness components than would be expected from mutation–selection balance. Here, we used a long‐term study of an individually marked population of North American red squirrels (Tamiasciurus hudsonicus) to look for evidence of (1) additive genetic variance in lifetime reproductive success and (2) fitness trade‐offs between fitness components, such as male and female fitness or fitness in high‐ and low‐resource environments. “Animal model” analyses of a multigenerational pedigree revealed modest maternal effects on fitness, but very low levels of additive genetic variance in lifetime reproductive success overall as well as fitness measures within each sex and environment. It therefore appears that there are very low levels of direct genetic variance in fitness and fitness components in red squirrels to facilitate contemporary adaptation in this population.     

EVOLUTION OF TRADE-OFFS BETWEEN SEXUAL AND ASEXUAL PHASES AND THE ROLE OF REPRODUCTIVE PLASTICITY IN THE GENETIC ARCHITECTURE OF APHID LIFE HISTORIES

Life‐history theory postulates that evolution is constrained by trade‐offs (i.e., negative genetic correlations) among traits that contribute to fitness. However, in organisms with complex life cycles, trade‐offs may drastically differ between phases, putatively leading to different evolutionary trajectories. Here, we tested this possibility by examining changes in life‐history traits in an aphid species that alternates asexual and sexual reproduction in its life cycle. The quantitative genetics of reproductive and dispersal traits was studied in 23 lineages (genotypes) of the bird cherry‐oat aphid Rhopalosiphum padi, during both the sexual and asexual phases, which were induced experimentally under specific environmental conditions. We found large and significant heritabilities (broad‐sense) for all traits and several negative genetic correlations between traits (trade‐offs), which are related to reproduction (i.e., numbers of the various sexual or asexual morphs) or dispersal (i.e., numbers of winged or wingless morphs). These results suggest that R. padi exhibits lineage specialization both in reproductive and dispersal strategies. In addition, we found important differences in the structure of genetic variance–covariance matrices ( G ) between phases. These differences were due to two large, negative genetic correlations detected during the asexual phase only: (1) between fecundity and age at maturity and (2) between the production of wingless and winged parthenogenetic females. We propose that this differential expression in genetic architecture results from a reallocation scheme during the asexual phase, when sexual morphs are not produced. We also found significant G × E interaction and nonsignificant genetic correlations across phases, indicating that genotypes could respond independently to selection in each phase. Our results reveal a rather unique situation in which the same population and even the same genotypes express different genetic (co)variation under different environmental conditions, driven by optimal resource allocation criteria.     

QUANTIFYING EVOLUTIONARY POTENTIAL OF MARINE FISH LARVAE: HERITABILITY,SELECTION, AND EVOLUTIONARY CONSTRAINTS

For many marine fish, intense larval mortality may provide considerable opportunity for selection, yet much less is known about the evolutionary potential of larval traits. We combined field demographic studies and manipulative experiments to estimate quantitative genetic parameters for both larval size and swimming performance for a natural population of a common coral‐reef fish, the bicolor damselfish (Stegastes partitus). We also examined selection on larval size by synthesizing information from published estimates of selective mortality. We introduce a method that uses the Lande–Arnold framework for examining selection on quantitative traits to empirically reconstruct adaptive landscapes. This method allows the relationship between phenotypic value and fitness components to be described across a broad range of trait values. Our results suggested that despite strong viability selection for large larvae and moderate heritability (h²= 0.29), evolutionary responses of larvae would likely be balanced by reproductive selection favoring mothers that produce more, smaller offspring. Although long‐term evolutionary responses of larval traits may be constrained by size‐number trade‐offs, our results suggest that phenotypic variation in larval size may be an ecologically important source of variability in population dynamics through effects on larval survival and recruitment to benthic populations.     

STRONG SELECTION GENOME‐WIDE ENHANCES FITNESS TRADE‐OFFS ACROSS ENVIRONMENTS AND EPISODES OF SELECTION

Fitness trade‐offs across episodes of selection and environments influence life‐history evolution and adaptive population divergence. Documenting these trade‐offs remains challenging as selection can vary in magnitude and direction through time and space. Here, we evaluate fitness trade‐offs at the levels of the whole organism and the quantitative trait locus (QTL) in a multiyear field study of Boechera stricta (Brassicaceae), a genetically tractable mustard native to the Rocky Mountains. Reciprocal local adaptation was pronounced for viability, but not for reproductive components of fitness. Instead, local genomes had a fecundity advantage only in the high latitude garden. By estimating realized selection coefficients from individual‐level data on viability and reproductive success and permuting the data to infer significance, we examined the genetic basis of fitness trade‐offs. This analytical approach (Conditional Neutrality‐Antagonistic Pleiotropy, CNAP) identified genetic trade‐offs at a flowering phenology QTL (costs of adaptation) and revealed genetic trade‐offs across fitness components (costs of reproduction). These patterns would not have emerged from traditional ANOVA‐based QTL mapping. Our analytical framework can be applied to other systems to investigate fitness trade‐offs. This task is becoming increasingly important as climate change may alter fitness landscapes, potentially disrupting fitness trade‐offs that took many generations to evolve.     

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.     

THE GENOMIC ARCHITECTURE OF SEXUAL DIMORPHISM IN THE DIOECIOUS PLANT SILENE LATIFOLIA

Evaluating the genetic architecture of sexual dimorphism can aid our understanding of the extent to which shared genetic control of trait variation versus sex‐specific control impacts the evolutionary dynamics of phenotypic change within each sex. We performed a QTL analysis on Silene latifolia to evaluate the contribution of sex‐specific QTL to phenotypic variation in 46 traits, whether traits involved in trade‐offs had colocalized QTL, and whether the distribution of sex‐specific loci can explain differences between the sexes in their variance/covariance matrices. We used a backcross generation derived from two artificial‐selection lines. We found that sex‐specific QTL explained a significantly greater percent of the variation in sexually dimorphic traits than loci expressed in both sexes. Genetically correlated traits often had colocalized QTL, whose signs were in the expected direction. Lastly, traits with different genetic correlations within the sexes displayed a disproportionately high number of sex‐specific QTL, and more QTL co‐occurred in males than females, suggesting greater trait integration. These results show that sex differences in QTL patterns are congruent with theory on the resolution of sexual conflict and differences based on G ‐matrix results. They also suggest that trade‐offs and trait integration are likely to affect males more than females.     

Revealing hidden evolutionary capacity to cope with global change

The extent to which global change will impact the long‐term persistence of species depends on their evolutionary potential to adapt to future conditions. While the number of studies that estimate the standing levels of adaptive genetic variation in populations under predicted global change scenarios is growing all the time, few studies have considered multiple environments simultaneously and even fewer have considered evolutionary potential in multivariate context. Because conditions will not be constant, adaptation to climate change is fundamentally a multivariate process so viewing genetic variances and covariances over multivariate space will always be more informative than relying on bivariate genetic correlations between traits. A multivariate approach to understanding the evolutionary capacity to cope with global change is necessary to avoid misestimating adaptive genetic variation in the dimensions in which selection will act. We assessed the evolutionary capacity of the larval stage of the marine polychaete Galeolaria caespitosa to adapt to warmer water temperatures. Galeolaria is an important habitat‐forming species in Australia, and its earlier life‐history stages tend to be more susceptible to stress. We used a powerful quantitative genetics design that assessed the impacts of three temperatures on subsequent survival across over 30 000 embryos across 204 unique families. We found adaptive genetic variation in the two cooler temperatures in our study, but none in the warmest temperature. Based on these results, we would have concluded that this species has very little capacity to evolve to the warmest temperature. However, when we explored genetic variation in multivariate space, we found evidence that larval survival has the potential to evolve even in the warmest temperatures via correlated responses to selection across thermal environments. Future studies should take a multivariate approach to estimating evolutionary capacity to cope with global change lest they misestimate a species’ true adaptive potential.     

REACTION NORMS OF MORPHOLOGICAL AND LIFE-HISTORY TRAITS TO LIGHT AVAILABILITY IN IMPATIENS CAPENSIS

For plants, light availability is an important environmental factor that varies both within and between populations. Although the existence of sun and shade “ecotypes” is controversial, it is often assumed that trade-offs may exist between performance in sun and in shade. This study therefore investigated variation in reaction norms to light availability within and between two neighboring natural populations of the annual Impatiens capensis, one in full sun and the other in a forest understory. Seedlings were collected randomly from both populations and grown to maturity in a greenhouse under two light conditions: full light and 18% of full light. Selfed full-sib seed families were collected from plants from both populations grown in both parental light environments. To characterize family reaction norms, seedlings from each family were divided into the same two light treatments and individuals were scored for a variety of morphological and life-history traits. The maternal light environment had little impact on progeny reaction norms. However, the two study populations differed both qualitatively and quantitatively in plastic response to light availability (indicated by significant population x environment interactions in mixed-model ANCOVA). Much of this difference was attributable to population differences in light sensitivity of axillary meristem allocation patterns, which produced concurrent differences in reaction norms for a suite of developmentally linked traits. Within each population, different sets of traits displayed significant variation in plasticity (indicated by significant family x environment interactions). Thus, the genetic potential for evolutionary response to selection in heterogeneous light environments may differ dramatically between neighboring plant populations. Between-environment genetic correlations were largely positive in the woods population and positive or nonsignificant in the sun population; there was no evidence for performance trade-offs across environments or sun or shade “specialist” genotypes within either population. There was little evidence that population differences represented adaptive differentiation for sun or shade; rather, the results suggested the hypothesis of differential selection on patterns of meristem allocation caused by population differences in timing of mortality and intensity of competition.     

SENSITIVITY OF THE DISTRIBUTION OF MUTATIONAL FITNESS EFFECTS TO ENVIRONMENT,GENETIC BACKGROUND,AND ADAPTEDNESS: A CASE STUDY WITH DROSOPHILA

Heterogeneity in the fitness effects of individual mutations has been found across different environmental and genetic contexts. Going beyond effects on individual mutations, how is the distribution of selective effects, f(s), altered by changes in genetic and environmental context? In this study, we examined changes in the major features of f(s) by estimating viability selection on 36 individual mutations in Drosophila melanogaster across two different environments in two different genetic backgrounds that were either adapted or nonadapted to the two test environments. Both environment and genetic background affected selection on individual mutations. However, the overall distribution f(s) appeared robust to changes in genetic background but both the mean, E(s), and the variance, V(s) were dependent on the environment. Between these two properties, V(s) was more sensitive to environmental change. Contrary to predictions of fitness landscape theory, the match between genetic background and assay environment (i.e., adaptedness) had little effect on f(s).     

Negative frequency‐dependent selection maintains coexisting genotypes during fluctuating selection

Natural environments are rarely static; rather selection can fluctuate on timescales ranging from hours to centuries. However, it is unclear how adaptation to fluctuating environments differs from adaptation to constant environments at the genetic level. For bacteria, one key axis of environmental variation is selection for planktonic or biofilm modes of growth. We conducted an evolution experiment with Burkholderia cenocepacia, comparing the evolutionary dynamics of populations evolving under constant selection for either biofilm formation or planktonic growth with populations in which selection fluctuated between the two environments on a weekly basis. Populations evolved in the fluctuating environment shared many of the same genetic targets of selection as those evolved in constant biofilm selection, but were genetically distinct from the constant planktonic populations. In the fluctuating environment, mutations in the biofilm‐regulating genes wspA and rpfR rose to high frequency in all replicate populations. A mutation in wspA first rose rapidly and nearly fixed during the initial biofilm phase but was subsequently displaced by a collection of rpfR mutants upon the shift to the planktonic phase. The wspA and rpfR genotypes coexisted via negative frequency‐dependent selection around an equilibrium frequency that shifted between the environments. The maintenance of coexisting genotypes in the fluctuating environment was unexpected. Under temporally fluctuating environments, coexistence of two genotypes is only predicted under a narrow range of conditions, but the frequency‐dependent interactions we observed provide a mechanism that can increase the likelihood of coexistence in fluctuating environments.