Evolution,plasticity and evolving plasticity of phenology in the tree species Alnus glutinosa

Both traits and the plasticity of these traits are subject to evolutionary change and therefore affect the long‐term persistence of populations and their role in local communities. We subjected clones from 12 different populations of Alnus glutinosa, located along a latitudinal gradient, to two different temperature treatments, to disentangle the distribution of genetic variation in timing of bud burst and bud burst plasticity within and among genotypes, populations, and regions. We calculated heritability and evolvability estimates for bud burst and bud burst plasticity and assessed the influence of divergent selection relative to neutral drift. We observed higher levels of heritability and evolvability for bud burst than for its plasticity, whereas the total phenological heritability and evolvability (i.e. combining timing of bud burst and bud burst plasticity) suggest substantial evolutionary potential with respect to phenology. Earlier bud burst was observed for the low‐latitudinal populations than for the populations from higher latitudes, whereas the high‐latitudinal populations did not show the expected delayed bud burst. This countergradient variation can be due to evolution towards increased phenological plasticity at higher latitudes. However, because we found little evidence for adaptive differences in phenological plasticity across the latitudinal gradient, we suggest differential frost tolerance as the most likely explanation for the observed phenological patterns in A. glutinosa.     

Relaxed predation results in reduced phenotypic integration in a suite of dragonflies

Although changes in magnitude of single traits responding to selective agents have been studied intensively, little is known about selection shaping networks of traits and their patterns of covariation. However, this is central for our understanding of phenotypic evolution as traits are embedded in a multivariate environment with selection affecting a multitude of traits simultaneously rather than individually. Here, we investigate inter‐ and intraspecific patterns of trait integration (trait correlations) in the larval abdomen of dragonflies as a response to a change in predator selection. Species of the dragonfly genus Leucorrhinia underwent a larval habitat shift from predatory fish to predatory dragonfly‐dominated lakes with an associated relaxation in selection pressure from fish predation. Our results indicate that the habitat‐shift‐induced relaxed selection pressure caused phenotypic integration of abdominal traits to be reduced. Intraspecific findings matched patterns comparing species from both habitats with higher abdominal integration in response to predatory fish. This higher integration is probably a result of faster burst swimming speed. The abdomen holds the necessary morphological machinery to successfully evade predatory fish via burst swimming. Hence, abdominal traits have to function in a tight coordinated manner, as maladaptive variation and consequently nonoptimal burst swimming would cause increased mortality. In predatory dragonfly‐dominated lakes, no such strong link between burst swimming and mortality is present. Our findings highlight the importance of studying multivariate trait relationships as a response to selection for understanding patterns of phenotypic diversification.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Evolutionary dynamics of the leaf phenological cycle in an oak metapopulation along an elevation gradient

It has been predicted that environmental changes will radically alter the selective pressures on phenological traits. Long‐lived species, such as trees, will be particularly affected, as they may need to undergo major adaptive change over only one or a few generations. The traits describing the annual life cycle of trees are generally highly evolvable, but nothing is known about the strength of their genetic correlations. Tight correlations can impose strong evolutionary constraints, potentially hampering the adaptation of multivariate phenological phenotypes. In this study, we investigated the evolutionary, genetic and environmental components of the timing of leaf unfolding and senescence within an oak metapopulation along an elevation gradient. Population divergence, estimated from in situ and common‐garden data, was compared to expectations under neutral evolution, based on microsatellite markers. This approach made it possible (1) to evaluate the influence of genetic correlation on multivariate local adaptation to elevation and (2) to identify traits probably exposed to past selective pressures due to the colder climate at high elevation. The genetic correlation was positive but very weak, indicating that genetic constraints did not shape the local adaptation pattern for leaf phenology. Both spring and fall (leaf unfolding and senescence, respectively) phenology timings were involved in local adaptation, but leaf unfolding was probably the trait most exposed to climate change‐induced selection. Our data indicated that genetic variation makes a much smaller contribution to adaptation than the considerable plastic variation displayed by a tree during its lifetime. The evolutionary potential of leaf phenology is, therefore, probably not the most critical aspect for short‐term population survival in a changing climate.     

Reaction norms of Arabidopsis (Brassicaceae). III. Response to nutrients in 26 populations from a worldwide collection

The study of phenotypic plasticity, the ability of a given genotype to express different phenotypes as environments change, is becoming a central focus of ecological genetics and evolutionary theory. To help address the most pressing questions about plasticity (its genetic control, ecological relevance, and macroevolutionary consequences) we advocate the use of Arabidopsis thaliana (and eventually other related species of the same genus) as a model system. In this study we present experimental data concerning: (a) the extent of reaction norm variation to two levels of nutrients in a worldwide collection of 26 A. thaliana populations; and (b) the existence of multivariate associations among key phenotypic characters, and their reaction to changes in the environment. We found significant among-population genetic variation for eight of the nine traits measured, as well as plasticity in four traits. Five traits showed significant differences in genetic variation between the two environments. The multivariate association of the nine traits defines four major groups of covarying characters, each of which may be plastic or not, depending on the particular population. The use of populations that can be easily obtained by any researcher, because they are part of a worldwide collection, implies that it will be easy to build on our results during future investigations of phenotypic plasticity in this species.     

Phenological matching rather than genetic variation in host preference underlies geographical variation in host plants used by orange tip butterflies

An insect species that shows variation in host species association across its geographical range may do so either because of local adaptation in host plant preference of the insect or through environmentally or genetically induced differences in the plants, causing variation in host plant suitability between regions. In the present study, we experimentally investigate the host plant preference of Anthocharis cardamines (orange tip butterfly) in two populations from the UK and two from Sweden. Previous reports indicate that A. cardamines larvae are found on different host plant species in different regions of the UK, and some variation has been reported in Sweden. Host plant choice trials showed that females prefer to oviposit on plants in an earlier phenological stage, as well as on larger plants. When controlling for plant phenological stage and size, the host species had no statistically significant effect on the choice of the females. Moreover, there were no differences in host plant species preference among the four butterfly populations. Based on our experiment, the oviposition choice by A. cardamines mainly depends on the phenological stage and the size of the host plant. This finding supports the idea that the geographical patterns of host–plant association of A. cardamines in the UK and Sweden are consequences of the phenology and availability of the local hosts, rather than regional genetic differences in the host species preference of the butterfly.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.     

Host‐associated genomic differentiation in congeneric butterflies: now you see it,now you do not

Ecotypic variation among populations may become associated with widespread genomic differentiation, but theory predicts that this should happen only under particular conditions of gene flow, selection and population size. In closely related species, we might expect the strength of host‐associated genomic differentiation (HAD) to be correlated with the degree of phenotypic differentiation in host‐adaptive traits. Using microsatellite and Amplified Fragment Length Polymorphism (AFLP) markers, and controlling for isolation by distance between populations, we sought HAD in two congeneric species of butterflies with different degrees of host plant specialization. Prior work on Euphydryas editha had shown strong interpopulation differentiation in host‐adapted traits, resulting in incipient reproductive isolation among host‐associated ecotypes. We show here that Euphydryas aurinia had much weaker host‐associated phenotypic differentiation. Contrary to our expectations, we detected HAD in Euphydryas aurinia, but not in E. editha. Even within an E. aurinia population that fed on both hosts, we found weak but significant sympatric HAD that persisted in samples taken 9 years apart. The finding of significantly stronger HAD in the system with less phenotypic differentiation may seem paradoxical. Our findings can be explained by multiple factors, ranging from differences in dispersal or effective population size, to spatial variation in genomic or phenotypic traits and to structure induced by past histories of host‐adapted populations. Other infrequently measured factors, such as differences in recombination rates, may also play a role. Our result adds to recent work as a further caution against assumptions of simple relationships between genomic and adaptive phenotypic differentiation.     

Adcyap1 polymorphism covaries with breeding latitude in a Nearctic migratory songbird,the Wilson's warbler (Cardellina pusilla)

Understanding the genetic background of complex behavioral traits, showing multigenic control and extensive environmental effects, is a challenging task. Among such traits, migration is known to show a large additive genetic component. Yet, the identification of specific genes or gene regions explaining phenotypic variance in migratory behavior has received less attention. Migration ultimately depends on seasonal cycles, and polymorphism at phenological candidate genes may underlie variation in timing of migration or other aspects of migratory behavior. In this study of a Nearctic–Neotropical migratory songbird, the Wilson's warbler (Cardellina pusilla), we investigated the association between polymorphism at two phenological candidate genes, Clock and Adcyap1, and two aspects of the migratory phenotype, timing of spring migration through a stopover site and inferred latitude of the breeding destination. The breeding destination of migrating individuals was identified using feather deuterium ratio (δ²H), which reliably reflects breeding latitude throughout the species' western breeding range. Ninety‐eight percent of the individuals were homozygous at Clock, and the rare heterozygotes did not deviate from homozygous migration phenology. Adcyap1 was highly polymorphic, and allele size was not significantly associated with migration date. However, Adcyap1 allele size significantly positively predicted the inferred breeding latitude of males but not of females. Moreover, we found a strong positive association between inferred breeding latitude and Adcyap1 allele size in long‐distance migrating birds from the northern sector of the breeding range (western Canada), while this was not the case in short‐distance migrating birds from the southern sector of the breeding range (coastal California). Our findings support previous evidence for a role of Adcyap1 in shaping the avian migratory phenotype, while highlighting that patterns of phenological candidate gene–phenotype associations may be complex, significantly varying between geographically distinct populations and even between the sexes.     

Polymorphism at the Clock gene predicts phenology of long‐distance migration in birds

Dissecting phenotypic variance in life history traits into its genetic and environmental components is at the focus of evolutionary studies and of pivotal importance to identify the mechanisms and predict the consequences of human‐driven environmental change. The timing of recurrent life history events (phenology) is under strong selection, but the study of the genes that control potential environmental canalization in phenological traits is at its infancy. Candidate genes for circadian behaviour entrained by photoperiod have been screened as potential controllers of phenological variation of breeding and moult in birds, with inconsistent results. Despite photoperiodic control of migration is well established, no study has reported on migration phenology in relation to polymorphism at candidate genes in birds. We analysed variation in spring migration dates within four trans‐Saharan migratory species (Luscinia megarhynchos; Ficedula hypoleuca; Anthus trivialis; Saxicola rubetra) at a Mediterranean island in relation to Clock and Adcyap1 polymorphism. Individuals with larger number of glutamine residues in the poly‐Q region of Clock gene migrated significantly later in one or, respectively, two species depending on sex and whether the within‐individual mean length or the length of the longer Clock allele was considered. The results hinted at dominance of the longer Clock allele. No significant evidence for migration date to covary with Adcyap1 polymorphism emerged. This is the first evidence that migration phenology is associated with Clock in birds. This finding is important for evolutionary studies of migration and sheds light on the mechanisms that drive bird phenological changes and population trends in response to climate change.     

Selection and plasticity both account for interannual variation in life‐history phenology in an annual prairie legume

As the environment changes, so too must plant communities and populations if they are to persist. Life‐history transitions and their timing are often the traits that are most responsive to changing environmental conditions. To compare the contributions of plasticity and natural selective response to variation in germination and flowering phenology, we performed a quantitative genetic study of phenotypic selection on Chamaecrista fasciculata (Fabaceae) across two consecutive years in a restored tallgrass prairie. The earliest dates of germination and flowering were recorded for two parental cohorts and one progeny cohort in an experimental garden. Environmental differences between years were the largest contributors to phenological variation in this population. In addition, there was substantial heritability for flowering time and statistically significant selection for advancement of flowering. Comparison between a progeny cohort and its preselection parental cohort indicated a change in mean flowering time consistent with the direction of selection. Selection on germination time was weaker than that on flowering time, while environmental effects on germination time were stronger. The response to selection on flowering time was detectable when accounting for the effect of the environment on phenotypic differences, highlighting the importance of controlling for year‐to‐year environmental variation in quantitative genetic studies.     

Parallel adaptation to climate above the 35th parallel

Independent or parallel evolution of similar traits is key to understanding the genetics and limitations of adaptation. Adaptation from the same genetic changes in different populations defines parallel evolution. Such genetic changes can derive from standing ancestral variation or de novo mutations and excludes instances of adaptive introgression. In this issue of Molecular Ecology, Walden et al.(2020) investigate the scale of parallel climate adaptation from standing genetic variation between two North American Arabidopsis lyrata lineages, each formed by a distinct evolutionary history during the last glacial cycle. By identifying adaptive variants correlated with three ecologically significant climatic gradients, they show that instead of the same genetic variants or even genes, parallel evolution is only observed at the level of biological processes. The evolution of independent adaptive variants to climate in two genetically close lineages is explained by their different post‐glacial demographic histories. Separate glacial refugia and strong population bottlenecks were probably sufficient to change the landscape of shared allele frequencies, hindering the possibility of parallel evolution.     

Patterns of univariate and multivariate plasticity to elevated carbon dioxide in six European populations of Arabidopsis thaliana

The impact of elevated carbon dioxide on plants is a growing concern in evolutionary ecology and global change biology. Characterizing patterns of phenotypic integration and multivariate plasticity to elevated carbon dioxide can provide insights into ecological and evolutionary dynamics in future human‐altered environments. Here, we examined univariate and multivariate responses to carbon enrichment in six functional traits among six European accessions of Arabidopsis thaliana. We detected phenotypic plasticity in both univariate and multivariate phenotypes, but did not find significant variation in plasticity (genotype by environment interactions) within or among accessions. Eigenvector, eigenvalue variance, and common principal components analyses showed that elevated carbon dioxide altered patterns of trait covariance, reduced the strength of phenotypic integration, and decreased population‐level differentiation in the multivariate phenotype. Our data suggest that future carbon dioxide conditions may influence evolutionary dynamics in natural populations of A. thaliana.     

Geographical differentiation in floral traits across the distribution range of the Patagonian oil-secreting Calceolaria polyrhiza: do pollinators matter?

Background and Aims

The underlying evolutionary processes of pollinator-driven floral diversification are still poorly understood. According to the Grant–Stebbins model speciation begins with adaptive local differentiation in the response to spatial heterogeneity in pollinators. Although this crucial process links the micro- and macroevolution of floral adaptation, it has received little attention. In this study geographical phenotypic variation was investigated in Patagonian Calceolaria polyrhiza and its pollinators, two oil-collecting bee species that differ in body size and geographical distribution.

Methods

Patterns of phenotypic variation were examined together with their relationships with pollinators and abiotic factors. Six floral and seven vegetative traits were measured in 45 populations distributed across the entire species range. Climatic and edaphic parameters were determined for 25 selected sites, 2–16 bees per site of the most frequent pollinator species were captured, and a critical flower–bee mechanical fitting trait involved in effective pollination was measured. Geographical patterns of phenotypic and environmental variation were examined using uni- and multivariate analyses. Decoupled geographical variation between corolla area and floral traits related to the mechanical fit of pollinators was explored using a Mantel test.

Key Results

The body length of pollinators and the floral traits related to mechanical fit were strongly correlated with each other. Geographical variation of the mechanical-fit-related traits was decoupled from variation in corolla size; the latter had a geographical pattern consistent with that of the vegetative traits and was mainly affected by climatic gradients.

Conclusions

The results are consistent with pollinators playing a key role in shaping floral phenotype at a geographical scale and promoting the differentiation of two floral ecotypes. The relationship between the critical floral-fit-related trait and bee length remained significant even in models that included various environmental variables and an allometric predictor (corolla area). The abiotic environment also has an important role, mainly affecting floral size. Decoupled geographical variation between floral mechanical-fit-related traits and floral size would represent a strategy to maintain plant–pollinator phenotypic matching in this environmentally heterogeneous area.     

Behavioral constraints on local adaptation and counter‐gradient variation: Implications for climate change

Resource allocation to growth, reproduction, and body maintenance varies within species along latitudinal gradients. Two hypotheses explaining this variation are local adaptation and counter‐gradient variation. The local adaptation hypothesis proposes that populations are adapted to local environmental conditions and are therefore less adapted to environmental conditions at other locations. The counter‐gradient variation hypothesis proposes that one population out performs others across an environmental gradient because its source location has greater selective pressure than other locations. Our study had two goals. First, we tested the local adaptation and counter‐gradient variation hypotheses by measuring effects of environmental temperature on phenotypic expression of reproductive traits in the burying beetle, Nicrophorus orbicollis Say, from three populations along a latitudinal gradient in a common garden experimental design. Second, we compared patterns of variation to evaluate whether traits covary or whether local adaptation of traits precludes adaptive responses by others. Across a latitudinal range, N. orbicollis exhibits variation in initiating reproduction and brood sizes. Consistent with local adaptation: (a) beetles were less likely to initiate breeding at extreme temperatures, especially when that temperature represents their source range; (b) once beetles initiate reproduction, source populations produce relatively larger broods at temperatures consistent with their local environment. Consistent with counter‐gradient variation, lower latitude populations were more successful at producing offspring at lower temperatures. We found no evidence for adaptive variation in other adult or offspring performance traits. This suite of traits does not appear to coevolve along the latitudinal gradient. Rather, response to selection to breed within a narrow temperature range may preclude selection on other traits. Our study highlights that N. orbicollis uses temperature as an environmental cue to determine whether to initiate reproduction, providing insight into how behavior is modified to avoid costly reproductive attempts. Furthermore, our results suggest a temperature constraint that shapes reproductive behavior.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Root architecture,plant size and soil nutrient variation in natural populations of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Phenotypic variation in ecologically important traits may vary at large and small geographic scales, and may be shaped by natural selection. Here our explicit aim is to evaluate phenotypic differentiation among local populations and examine its relationship with ecological edaphic and climatic features that could lead to local adaptation. We characterized six populations of the model plant Arabidopsis thaliana over 3 years in the field in its native range. At each site, we measured edaphic conditions and aboveground and belowground phenotypes. In addition, we grew plants from the six characterized populations in a common greenhouse along with an additional fifteen populations from the Iberian Peninsula to examine evolutionary and ecological differentiation among populations, and relationships between geographic and ecological distance to phenotypic differences among populations. Significant differences in aboveground and belowground traits, population density, and micro- and macro-nutrient soil concentrations were found among the field populations. In particular, root architectural traits differed significantly among field populations. Complex patterns of ecological differences among population and plant phenotypes emerged when examining edaphic conditions in the Extremadura region, and geographic and climate variables at a broader scale of the Iberian Peninsula. We report levels of phenotypic variation at the local scale comparable to those found at broad geographic scales and report that local edaphic conditions contribute to population-level phenotypic variation in root and shoot traits. To our knowledge, these are the first reports of among population root architectural variation from natural field populations for this model organism. We demonstrate how ecological features, such as soil nutrients, can be associated with the phenotypic variation of A. thaliana measured in natural populations and may contribute to adaptive differentiation at a local scale.     

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

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.