The extent of introgression between incipient Clarkia species is determined by temporal environmental variation and mating system.

Introgression is pervasive across the tree of life but varies across taxa, geography, and genomic regions. However, the factors modulating this variation and how they may be affected by global change are not well understood. Here, we used 200 genomes and a 15-y site-specific environmental dataset to investigate the effects of environmental variation and mating system divergence on the magnitude of introgression between a recently diverged outcrosser-selfer pair of annual plants in the genus Clarkia. These sister taxa diverged very recently and subsequently came into secondary sympatry where they form replicated contact zones. Consistent with observations of other outcrosser-selfer pairs, we found that introgression was asymmetric between taxa, with substantially more introgression from the selfer to the outcrosser. This asymmetry was caused by a bias in the direction of initial F1 hybrid formation and subsequent backcrossing. We also found extensive variation in the outcrosser's admixture proportion among contact zones, which was predicted nearly entirely by interannual variance in spring precipitation. Greater fluctuations in spring precipitation resulted in higher admixture proportions, likely mediated by the effects of spring precipitation on the expression of traits that determine premating reproductive isolation. Climate-driven hybridization dynamics may be particularly affected by global change, potentially reshaping species boundaries and adaptation to novel environments.

Bet Hedging Is Not Sufficient to Explain Germination Patterns of a Winter Annual Plant.

AbstractBet hedging consists of life history strategies that buffer against environmental variability by trading off immediate and long-term fitness. Delayed germination in annual plants is a classic example of bet hedging and is often invoked to explain low germination fractions. We examined whether bet hedging explains low and variable germination fractions among 20 populations of the winter annual plant Clarkia xantiana ssp. xantiana that experience substantial variation in reproductive success among years. Leveraging 15 years of demographic monitoring and 3 years of field germination experiments, we assessed the fitness consequences of seed banks and compared optimal germination fractions from a density-independent bet-hedging model to observed germination fractions. We did not find consistent evidence of bet hedging or the expected trade-off between arithmetic and geometric mean fitness, although delayed germination increased long-term fitness in 7 of 20 populations. Optimal germination fractions were two to five times higher than observed germination fractions, and among-population variation in germination fractions was not correlated with risks across the life cycle. Our comprehensive test suggests that bet hedging is not sufficient to explain the observed germination patterns. Understanding variation in germination strategies will likely require integrating bet hedging with complementary forces shaping the evolution of delayed germination.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge.

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

Direct tracking of pollen with quantum dots reveals surprising uniformity in dispersal distance across 11 populations of an annual plant.

PREMISE: Pollen movement is a crucial component of dispersal in seed plants. Although pollen dispersal is well studied, methodological constraints have made it challenging to directly track pollen flow within multiple populations across landscapes. We labeled pollen with quantum dots, a new technique that overcomes past limitations, to evaluate the spatial scale of pollen dispersal and its relationship with conspecific density within 11 populations of Clarkia xantiana subsp. xantiana, a bee-pollinated annual plant. METHODS: We used experimental arrays in two years to track pollen movement across distances of 5-35 m within nine populations and across distances of 10-70 m within two additional populations. We tested for distance decay of pollen dispersal, whether conspecific density modulated dispersal distance, and whether dispersal kernels varied among populations across an environmentally complex landscape. RESULTS: Labeled pollen receipt did not decline with distance over 35 m within eight of nine populations or over 70 m within either of two populations. Pollen receipt increased with conspecific density. Overall, dispersal kernels were consistent across populations. CONCLUSIONS: The surprising uniformity in dispersal distance within different populations was likely influenced by low precipitation and plant density in our study years. This suggests that spatiotemporal variation in the abiotic environment substantially influences the extent of gene flow within and among populations.

Rapid evolution during climate change: demographic and genetic constraints on adaptation to severe drought.

Populations often vary in their evolutionary responses to a shared environmental perturbation. A key hurdle in building more predictive models of rapid evolution is understanding this variation-why do some populations and traits evolve while others do not? We combined long-term demographic and environmental data, estimates of quantitative genetic variance components, a resurrection experiment and individual-based evolutionary simulations to gain mechanistic insights into contrasting evolutionary responses to a severe multi-year drought. We examined five traits in two populations of a native California plant, Clarkia xantiana, at three time points over 7 years. Earlier flowering phenology evolved in only one of the two populations, though both populations experienced similar drought severity and demographic declines and were estimated to have considerable additive genetic variance for flowering phenology. Pairing demographic and experimental data with evolutionary simulations suggested that while seed banks in both populations probably constrained evolutionary responses, a stronger seed bank in the non-evolving population resulted in evolutionary stasis. Gene flow through time via germ banks may be an important, underappreciated control on rapid evolution in response to extreme environmental perturbations.

The efficacy of selection may increase or decrease with selfing depending upon the recombination environment.

Much theory has focused on how a population's selfing rate affects the ability of natural selection to remove deleterious mutations from a population. However, most such theory has focused on mutations of a given dominance and fitness effect in isolation. It remains unclear how selfing affects the purging of deleterious mutations in a genome-wide context where mutations with different selection and dominance coefficients co-segregate. Here, we use individual-based forward simulations and analytical models to investigate how mutation, selection and recombination interact with selfing rate to shape genome-wide patterns of mutation accumulation and fitness. In addition to recovering previously described results for how selfing affects the efficacy of selection against mutations of a given dominance class, we find that the interaction of purifying selection against mutations of different dominance classes changes with selfing and recombination rates. In particular, when recombination is low and recessive deleterious mutations are common, outcrossing populations transition from purifying selection to pseudo-overdominance, dramatically reducing the efficacy of selection. At these parameter combinations, the efficacy of selection remains low until populations hit a threshold selfing rate, above which it increases. In contrast, selection is more effective in outcrossing than (partial) selfing populations when recombination rates are moderate to high and recessive deleterious mutations are rare.

The spatial scale of adaptation in a native annual plant and its implications for responses to climate change.

Spatial patterns of adaptation provide important insights into agents of selection and expected responses of populations to climate change. Robust inference into the spatial scale of adaptation can be gained through reciprocal transplant experiments that combine multiple source populations and common gardens. Here, we examine the spatial scale of local adaptation of the North American annual plant common ragweed, Ambrosia artemisiifolia, using data from four common gardens with 22 source populations sampled from across a ∼1200 km latitudinal gradient within the native range. We found evidence of local adaptation at the northernmost common garden, but maladaptation at the two southern gardens, where more southern source populations outperformed local populations. Overall, the spatial scale of adaptation was large-at the three gardens where distance between source populations and gardens explained variation in fitness, it took an average of 820 km for fitness to decline to 50% of its predicted maximum. Taken together, these results suggest that climate change has already caused maladaptation, especially across the southern portion of the range, and may result in northward range contraction over time.

Of mutualism and migration: will interactions with novel ericoid mycorrhizal communities help or hinder northward Rhododendron range shifts?

Rapid climate change imperils many small-ranged endemic species as the climate envelopes of their native ranges shift poleward. In addition to abiotic changes, biotic interactions are expected to play a critical role in plant species' responses. Below-ground interactions are of particular interest given increasing evidence of microbial effects on plant performance and the prevalence of mycorrhizal mutualisms. We used greenhouse mesocosm experiments to investigate how natural northward migration/assisted colonization of Rhododendron catawbiense, a small-ranged endemic eastern U.S. shrub, might be influenced by novel below-ground biotic interactions from soils north of its native range, particularly with ericoid mycorrhizal fungi (ERM). We compared germination, leaf size, survival, and ERM colonization rates of endemic R. catawbiense and widespread R. maximum when sown on different soil inoculum treatments: a sterilized control; a non-ERM biotic control; ERM communities from northern R. maximum populations; and ERM communities collected from the native range of R. catawbiense. Germination rates for both species when inoculated with congeners' novel soils were significantly higher than when inoculated with conspecific soils, or non-mycorrhizal controls. Mortality rates were unaffected by treatment, suggesting that the unexpected reciprocal effect of each species' increased establishment in association with heterospecific ERM could have lasting demographic effects. Our results suggest that seedling establishment of R. catawbiense in northern regions outside its native range could be facilitated by the presence of extant congeners like R. maximum and their associated soil microbiota. These findings have direct relevance to the potential for successful poleward migration or future assisted colonization efforts.

Plant-soil interactions limit lifetime fitness outside a native plant's geographic range margin.

Plant species' distributions are often thought to overwhelmingly reflect their climatic niches. However, climate represents only a fraction of the n-dimensional environment to which plant populations adapt, and studies are increasingly uncovering strong effects of nonclimatic factors on species' distributions. We used a manipulative, factorial field experiment to quantify the effects of soil environment and precipitation (the putatively overriding climatic factor) on plant lifetime fitness outside the geographic range boundary of a native California annual plant. We grew plants outside the range edge in large mesocosms filled with soil from either within or outside the range, and plants were subjected to either a low (ambient) or high (supplemental) spring precipitation treatment. Soil environment had large effects on plant lifetime fitness that were similar in magnitude to the effects of precipitation. Moreover, mean fitness of plants grown with within-range soil in the low precipitation treatment approximated that of plants grown with beyond-range soil in the high precipitation treatment. The positive effects of within-range soil persisted in the second, wetter year of the experiment, though the magnitude of the soil effect was smaller than in the first, drier year. These results are the first we know of to quantify the effects of edaphic variation on plant lifetime fitness outside a geographic range limit and highlight the need to include factors other than climate in models of species' distributions.

Microbes, mutualism, and range margins: testing the fitness consequences of soil microbial communities across and beyond a native plant's range.

Interactions between plants and soil fungi and bacteria are ubiquitous and have large effects on individual plant fitness. However, the degree to which spatial variation in soil microbial communities modulates plant species' distributions remains largely untested. Using the California native plant Clarkia xantiana ssp. xantiana we paired glasshouse and field reciprocal transplants of plant populations and soils to test whether plant-microbe interactions affect the plant's geographic range limit and whether there is local adaptation between plants and soil microbe communities. In the field and glasshouse, one of the two range interior inocula had a positive effect on plant fitness. In the field, this benefit was especially pronounced at the range edge and beyond, suggesting possible mutualist limitation. In the glasshouse, soil inocula from beyond-range tended to increase plant growth, suggesting microbial enemy release beyond the range margin. Amplicon sequencing revealed stark variation in microbial communities across the range boundary. Plants dispersing beyond their range limit are likely to encounter novel microbial communities. In C. x. xantiana, our results suggest that range expansion may be facilitated by fewer pathogens, but could also be hindered by a lack of mutualists. Both negative and positive plant-microbe interactions will likely affect contemporary range shifts.

Consequences of ignoring dispersal variation in network models for landscape connectivity.

Habitat loss and fragmentation can negatively influence population persistence and biodiversity, but the effects can be mitigated if species successfully disperse between isolated habitat patches. Network models are the primary tool for quantifying landscape connectivity, yet in practice, an overly simplistic view of species dispersal is applied. These models often ignore individual variation in dispersal ability under the assumption that all individuals move the same fixed distance with equal probability. We developed a modeling approach to address this problem. We incorporated dispersal kernels into network models to determine how individual variation in dispersal alters understanding of landscape-level connectivity and implemented our approach on a fragmented grassland landscape in Minnesota. Ignoring dispersal variation consistently overestimated a population's robustness to local extinctions and underestimated its robustness to local habitat loss. Furthermore, a simplified view of dispersal underestimated the amount of habitat substructure for small populations but overestimated habitat substructure for large populations. Our results demonstrate that considering biologically realistic dispersal alters understanding of landscape connectivity in ecological theory and conservation practice.

Consecuencias de la Omisión de la Variación en la Dispersión en los Modelos de Redes para la Conectividad de Paisajes Resumen La pérdida y la fragmentación del hábitat pueden influir negativamente la persistencia de poblaciones y biodiversidad. Sin embargo, estos efectos pueden ser mitigados si las especies tienen una dispersión exitosa entre los fragmentos aislados de hábitat. Los modelos de redes son la herramienta principal para la cuantificación de la conectividad del paisaje, no obstante en la práctica, se tiende a usar una visión excesivamente simplista de la dispersión de especies. Es común que estos modelos ignoren la variación que existe entre individuos en sus habilidades de dispersión y que asuman que todos los individuos se pueden mover la misma distancia y con la misma probabilidad. En este estudio, desarrollamos una estrategia de modelaje para (minimizar o aminorar) estas limitaciones incorporando kernels de dispersión dentro de los modelos de redes para determinar cómo la variación individual de la dispersión altera el entendimiento de la conectividad a nivel de paisaje. Como un ejemplo, implementamos esta estrategia en un paisaje de pastizal fragmentado en Minnesota. Omitir la variación en la dispersión generó una sobreestimación sistemática de la robustez de la población ante las extinciones locales y una subestimación de la robustez ante la pérdida local del hábitat. Además, una visión simplificada de la dispersión subestimó la complejidad de hábitat para las poblaciones pequeñas, sin emgargo sobreestimó la complejidad para las poblaciones grandes. Nuestros resultados demuestran que incorporar parámetros que describan una dispersión biológica realista tiene implicaciones importantes en la teoría de conectividad de paisajes e implementación de practicas de conservación.

The opportunity for outcrossing varies across the geographic range of the primarily selfing Clarkia xantiana ssp. parviflora.

PREMISE: The timing of self-fertilization has potentially important consequences for the trajectory of mating system evolution, the opportunity for outcrossing, and the maintenance of genetic variation in populations. For primarily selfing taxa, it remains poorly understood as to how floral variation influences the opportunity for outcrossing and whether those attributes vary among populations across geographic ranges. METHODS: We examined variation in floral traits (herkogamy, protandry, flower size, stigma stage at anthesis, timing of stigma receptivity) in seven populations of Clarkia xantiana ssp. parviflora, a primarily selfing taxon, spanning from the western to eastern margins of its distribution. We also performed experimental emasculations and pollinations (followed by stigma severing) to quantify the extent of opportunities for outcrossing across flower development. RESULTS: There was marked among-population variation in all floral traits, particularly between far eastern and western populations. Emasculation experiments showed that the eastern populations had minimal autonomous selfing, but western populations had high rates of selfing within 24 h after anthesis. Population variation in autofertility was significantly predicted by floral trait variation, especially protandry and petal size. CONCLUSIONS: Greater protandry both extended the time over which outcrossing could potentially occur and reduced the probability of autonomous selfing, suggesting that there may be a tradeoff that results in fitness loss when pollinator visits are not common. The east-west pattern of differentiation in some floral traits parallels that of postglacial range expansion, suggesting that selection on the mating system may have been strong in the process of range expansion.

Pleiotropy facilitates local adaptation to distant optima in common ragweed (Ambrosia artemisiifolia).

Pleiotropy, the control of multiple phenotypes by a single locus, is expected to slow the rate of adaptation by increasing the chance that beneficial alleles also have deleterious effects. However, a prediction arising from classical theory of quantitative trait evolution states that pleiotropic alleles may have a selective advantage when phenotypes are distant from their selective optima. We examine the role of pleiotropy in regulating adaptive differentiation among populations of common ragweed (Ambrosia artemisiifolia); a species that has recently expanded its North American range due to human-mediated habitat change. We employ a phenotype-free approach by using connectivity in gene networks as a proxy for pleiotropy. First, we identify loci bearing footprints of local adaptation, and then use genotype-expression mapping and co-expression networks to infer the connectivity of the genes. Our results indicate that the putatively adaptive loci are highly pleiotropic, as they are more likely than expected to affect the expression of other genes, and they reside in central positions within the gene networks. We propose that the conditionally advantageous alleles at these loci avoid the cost of pleiotropy by having large phenotypic effects that are beneficial when populations are far from their selective optima. We further use evolutionary simulations to show that these patterns are in agreement with a model where populations face novel selective pressures, as expected during a range expansion. Overall, our results suggest that highly connected genes may be targets of positive selection during environmental change, even though they likely experience strong purifying selection in stable selective environments.

Context Dependence of Local Adaptation to Abiotic and Biotic Environments: A Quantitative and Qualitative Synthesis.

Understanding how spatially variable selection shapes adaptation is an area of long-standing interest in evolutionary ecology. Recent meta-analyses have quantified the extent of local adaptation, but the relative importance of abiotic and biotic factors in driving population divergence remains poorly understood. To address this gap, we combined a quantitative meta-analysis and a qualitative metasynthesis to (1) quantify the magnitude of local adaptation to abiotic and biotic factors and (2) characterize major themes that influence the motivation and design of experiments that seek to test for local adaptation. Using local-foreign contrasts as a metric of local adaptation (or maladaptation), we found that local adaptation was greater in the presence than in the absence of a biotic interactor, especially for plants. We also found that biotic environments had stronger effects on fitness than abiotic environments when ignoring whether those environments were local versus foreign. Finally, biotic effects were stronger at low latitudes, and abiotic effects were stronger at high latitudes. Our qualitative analysis revealed that the lens through which local adaptation has been examined differs for abiotic and biotic factors. It also revealed biases in the design and implementation of experiments that make quantitative results challenging to interpret and provided directions for future research.

Maladaptation beyond a geographic range limit driven by antagonistic and mutualistic biotic interactions across an abiotic gradient.

Species' geographic range limits often result from maladaptation to the novel environments beyond the range margin. However, we rarely know which aspects of the n-dimensional environment are driving this maladaptation. Especially of interest is the influence of abiotic versus biotic factors in delimiting species' distributions. We conducted a 2-year reciprocal transplant experiment involving manipulations of the biotic environment to explore how spatiotemporal gradients in precipitation, fatal mammalian herbivory, and pollination affected lifetime fitness within and beyond the range of the California annual plant, Clarkia xantiana ssp. xantiana. In the first, drier year of the experiment, fitness outside the range edge was limited mainly by low precipitation, and there was some evidence for local adaptation within the range. In the second, wetter year, we did not observe abiotic limitations to plant fitness outside the range; instead biotic interactions, especially herbivory, limited fitness outside the range. Together, protection from herbivory and supplementation of pollen resulted in three- to sevenfold increases in lifetime fitness outside the range margin in the abiotically benign year. Overall, our work demonstrates the importance of biotic interactions, particularly as they interact with the abiotic environment, in determining fitness beyond geographic range boundaries.

Does adaptation to historical climate shape plant responses to future rainfall patterns? A rainfall manipulation experiment with common ragweed.

Climate change is affecting both the volume and distribution of precipitation, which in turn is expected to affect the growth and reproduction of plant populations. The near ubiquity of local adaptation suggests that adaptive differentiation may have important consequences for how populations are affected by and respond to changing precipitation. Here, we manipulated rainfall in a common garden to examine how differentiation among populations of common ragweed, Ambrosia artemisiifolia (Asteraceae) affects responses to water availability expected under climate change. We collected seeds from 26 populations along gradients of historical rainfall and used event-based rainout shelters and watering additions to simulate drier summer conditions and more extreme rainfall events, respectively. Ambrosia artemisiifolia had higher fitness on average under reduced rainfall, suggesting it may spread and become more abundant in areas projected to become hotter and drier during the summer months. We also found strong evidence for phenotypic and fitness clines across both latitude and longitude, and that phenological responses and fitness effects of altered rainfall depended on seed source or historical climate. The effect of rainfall treatment on female fitness was highest in western and mid longitudes, but there was little effect on eastern populations. Across latitude, the effect of rainfall treatment on male fitness was highest in southern populations. These phenology and fitness clines suggest that adaptive differentiation across the species' range has the potential to shape future responses of A. artemisiifolia populations to climate change, particularly altered patterns of rainfall.

Biotic Interactions Contribute to the Geographic Range Limit of an Annual Plant: Herbivory and Phenology Mediate Fitness beyond a Range Margin.

Species' geographic distributions have already shifted during the Anthropocene. However, we often do not know what aspects of the environment drive range dynamics, much less which traits mediate organisms' responses to these environmental gradients. Most studies focus on possible climatic limits to species' distributions and have ignored the role of biotic interactions, despite theoretical support for their importance in setting distributional limits. We used field experiments and simulations to estimate contributions of mammalian herbivory to a range boundary in the Californian annual plant Clarkia xantiana ssp. xantiana. A steep gradient of increasing probability of herbivory occurred across the boundary, and a reanalysis of prior transplant experiments revealed that herbivory drove severalfold declines in lifetime fitness at and beyond the boundary. Simulations showed that populations could potentially persist beyond the range margin in the absence of herbivory. Using data from a narrowly sympatric subspecies, Clarkia xantiana parviflora, we also showed that delayed phenology is strongly associated with C. xantiana ssp. xantiana's susceptibility to herbivory and low fitness beyond its border. Overall, our results provide some of the most comprehensive evidence to date of how the interplay of demography, traits, and spatial gradients in species interactions can produce a geographic range limit, and they lend empirical support to recent developments in range limits theory.

Seed predation increases from the Arctic to the Equator and from high to low elevations.

Species interactions have long been predicted to increase in intensity toward the tropics and low elevations because of gradients in climate, productivity, or biodiversity. Despite their importance for understanding global ecological and evolutionary processes, plant-animal interaction gradients are particularly difficult to test systematically across large geographic gradients, and evidence from smaller, disparate studies is inconclusive. By systematically measuring postdispersal seed predation using 6995 standardized seed depots along 18 mountains in the Pacific cordillera, we found that seed predation increases by 17% from the Arctic to the Equator and by 17% from 4000 meters above sea level to sea level. Clines in total predation, likely driven by invertebrates, were consistent across treeline ecotones and within continuous forest and were better explained by climate seasonality than by productivity, biodiversity, or latitude. These results suggest that species interactions play predictably greater ecological and evolutionary roles in tropical, lowland, and other less seasonal ecosystems.

Species distribution models throughout the invasion history of Palmer amaranth predict regions at risk of future invasion and reveal challenges with modeling rapidly shifting geographic ranges.

Palmer amaranth (Amaranthus palmeri) is an annual plant native to the desert Southwest of the United States and Mexico and has become invasive and caused large economic losses across much of the United States. In order to examine the temporal and spatial dynamics of past invasion, and to predict future invasion, we developed a broad array of species distribution models (SDMs). In particular, we constructed sequential SDMs throughout the invasion history and asked how well those predicted future invasion (1970 to present). We showed that invasion occurred from a restricted set of environments in the native range to a diverse set in the invaded range. Spatial autocorrelation analyses indicated that rapid range expansion was facilitated by stochastic, long-distance dispersal events. Regardless of SDM approach, all SDMs built using datasets from early in the invasion (1970-2010) performed poorly and failed to predict most of the current invaded range. Together, these results suggest that climate is unlikely to have influenced early stages of range expansion. SDMs that incorporated data from the most recent sampling (2011-2017) performed considerably better, predicted high suitability in regions that have recently become invaded, and identified mean annual temperature as a key factor limiting northward range expansion. Under future climates, models predicted both further northward range expansion and significantly increased suitability across large portions of the U.S. Overall, our results indicate significant challenges for SDMs of invasive species far from climate equilibrium. However, our models based on recent data make more robust predictions for northward range expansion of A. palmeri with climate change.

Mycorrhizal interactions do not influence plant-herbivore interactions in populations of Clarkia xantiana ssp. xantiana spanning from center to margin of the geographic range.

Multispecies interactions can be important to the expression of phenotypes and in determining patterns of individual fitness in nature. Many plants engage in symbiosis with arbuscular mycorrhizal fungi (AMF), but the extent to which AMF modulate other species interactions remains poorly understood. We examined multispecies interactions among plants, AMF, and insect herbivores under drought stress using a greenhouse experiment and herbivore choice assays. The experiment included six populations of Clarkia xantiana (Onagraceae), which span a complex environmental gradient in the Southern Sierra Nevada of California. Clarkia xantiana's developing fruits are commonly attacked by grasshoppers at the end of the growing season, and the frequency of attack is more common in populations from the range center than range margin. We found that AMF negatively influenced all metrics of plant growth and reproduction across all populations, presumably because plants supplied carbon to AMF but did not benefit substantially from resources potentially supplied by the AMF. The fruits of plants infected with AMF did not differ from those without AMF in their resistance to grasshoppers. There was significant variation among populations in damage from herbivores but did not reflect the center-to-margin pattern of herbivory observed in the field. In sum, our results do not support the view that AMF interactions modulate plant-herbivore interactions in this system.