Plant physical defenses contribute to a latitudinal gradient in resistance to insect herbivory within a widespread perennial grass

PREMISE

Herbivore pressure can vary across the range of a species, resulting in different defensive strategies. If herbivory is greater at lower latitudes, plants may be better defended there, potentially driving a latitudinal gradient in defense. However, relationships that manifest across the entire range of a species may be confounded by differences within genetic subpopulations, which may obscure the drivers of these latitudinal gradients.

METHODS

We grew plants of the widespread perennial grass Panicum virgatum in a common garden that included genotypes from three genetic subpopulations spanning an 18.5° latitudinal gradient. We then assessed defensive strategies of these plants by measuring two physical resistance traits-leaf mass per area (LMA) and leaf ash, a proxy for silica-and multiple measures of herbivory by caterpillars of the generalist herbivore fall armyworm (Spodoptera frugiperda).

RESULTS

Across all genetic subpopulations, low-latitude plants experienced less herbivory than high-latitude plants. Within genetic subpopulations, however, this relationship was inconsistent-the most widely distributed and phenotypically variable subpopulation (Atlantic) exhibited more consistent latitudinal trends than either of the other two subpopulations. The two physical resistance traits, LMA and leaf ash, were both highly heritable and positively associated with resistance to different measures of herbivory across all subpopulations, indicating their importance in defense against herbivores. Again, however, these relationships were inconsistent within subpopulations.

CONCLUSIONS

Defensive gradients that occur across the entire species range may not arise within localized subpopulations. Thus, identifying the drivers of latitudinal gradients in herbivory defense may depend on adequately sampling the diversity within a species.

Demographic consequences of an extreme heat wave are mitigated by spatial heterogeneity in an annual monkeyflower

Heat waves are becoming more frequent and intense with climate change, but the demographic and evolutionary consequences of heat waves are rarely investigated in herbaceous plant species. We examine the consequences of a short but extreme heat wave in Oregon populations of the common yellow monkeyflower (Mimulus guttatus) by leveraging a common garden experiment planted with range-wide populations and observational studies of 11 local populations. In the common garden, 89% of seedlings died during the heat wave including >96% of seedlings from geographically local populations. Some populations from hotter and drier environments had higher fitness, however, others from comparable environments performed poorly. Observational studies of local natural populations drastically differed in the consequences of the heat wave-one population was completely extirpated and nearly half had a >50% decrease in fitness. However, a few populations had greater fitness during the heat wave year. Differences in mortality corresponded to the impact of the heat wave on soil moisture-retention of soil moisture throughout the heat wave led to greater survivorship. Our results suggest that not all populations experience the same intensity or degree of mortality during extreme events and such heterogeneity could be important for genetic rescue or to facilitate the distribution of adaptive variants throughout the region.

Unique drought resistance strategies occur among monkeyflower populations spanning an aridity gradient

PREMISE

Annual plants often exhibit drought-escape and avoidance strategies to cope with limited water availability. Determining the extent of variation and factors underlying the evolution of divergent strategies is necessary for determining population responses to more frequent and severe droughts.

METHODS

We leveraged five Mimulus guttatus populations collected across an aridity gradient within manipulative drought and quantitative genetics experiments to examine constitutive and terminal-drought induced responses in drought resistance traits.

RESULTS

Populations varied considerably in drought-escape- and drought-avoidance-associated traits. The most mesic population demonstrated a unique resource conservative strategy. Xeric populations exhibited extreme plasticity when exposed to terminal drought that included flowering earlier at shorter heights, increasing water-use efficiency, and shifting C:N ratios. However, plasticity responses also differed between populations, with two populations slowing growth rates and flowering at earlier nodes and another population increasing growth rate. While nearly all traits were heritable, phenotypic correlations differed substantially between treatments and often, populations.

CONCLUSIONS

Our results suggest drought resistance strategies of populations may be finely adapted to local patterns of water availability. Substantial plastic responses suggest that xeric populations can already acclimate to drought through plasticity, but populations not frequently exposed to drought may be more vulnerable.

Evolution in response to climate in the native and introduced ranges of a globally distributed plant

The extent to which species can adapt to spatiotemporal climatic variation in their native and introduced ranges remains unresolved. To address this, we examined how clines in cyanogenesis (hydrogen cyanide [HCN] production-an antiherbivore defense associated with decreased tolerance to freezing) have shifted in response to climatic variation in space and time over a 60-year period in both the native and introduced ranges of Trifolium repens. HCN production is a polymorphic trait controlled by variation at two Mendelian loci (Ac and Li). Using phenotypic assays, we estimated within-population frequencies of HCN production and dominant alleles at both loci (i.e., Ac and Li) from 10,575 plants sampled from 131 populations on five continents, and then compared these frequencies to those from historical data collected in the 1950s. There were no clear relationships between changes in the frequency of HCN production, Ac, or Li and changes in temperature between contemporary and historical samples. We did detect evidence of continued evolution to temperature gradients in the introduced range, whereby the slope of contemporary clines for HCN and Ac in relation to winter temperature became steeper than historical clines and more similar to native clines. These results suggest that cyanogenesis clines show no clear changes through time in response to global warming, but introduced populations continue to adapt to their contemporary environments.

Functional evidence supports adaptive plant chemical defense along a geographical cline

Environmental clines in organismal defensive traits are usually attributed to stronger selection by enemies at lower latitudes or near the host's range center. Nonetheless, little functional evidence has supported this hypothesis, especially for coevolving plants and herbivores. We quantified cardenolide toxins in seeds of 24 populations of common milkweed (Asclepias syriaca) across 13 degrees of latitude, revealing a pattern of increasing cardenolide concentrations toward the host's range center. The unusual nitrogen-containing cardenolide labriformin was an exception and peaked at higher latitudes. Milkweed seeds are eaten by specialist lygaeid bugs that are even more tolerant of cardenolides than the monarch butterfly, concentrating most cardenolides (but not labriformin) from seeds into their bodies. Accordingly, whether cardenolides defend seeds against these specialist bugs is unclear. We demonstrate that Oncopeltus fasciatus (Lygaeidae) metabolized two major compounds (glycosylated aspecioside and labriformin) into distinct products that were sequestered without impairing growth. We next tested several isolated cardenolides in vitro on the physiological target of cardenolides (Na⁺/K⁺-ATPase); there was little variation among compounds in inhibition of an unadapted Na⁺/K⁺-ATPase, but tremendous variation in impacts on that of monarchs and Oncopeltus. Labriformin was the most inhibitive compound tested for both insects, but Oncopeltus had the greater advantage over monarchs in tolerating labriformin compared to other compounds. Three metabolized (and stored) cardenolides were less toxic than their parent compounds found in seeds. Our results suggest that a potent plant defense is evolving by natural selection along a geographical cline and targets specialist herbivores, but is met by insect tolerance, detoxification, and sequestration.

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

Adaptive and nonadaptive causes of heterogeneity in genetic differentiation across the Mimulus guttatus genome

Genetic diversity becomes structured among populations over time due to genetic drift and divergent selection. Although population structure is often treated as a uniform underlying factor, recent resequencing studies of wild populations have demonstrated that diversity in many regions of the genome may be structured quite dissimilar to the genome-wide pattern. Here, we explored the adaptive and nonadaptive causes of such genomic heterogeneity using population-level, whole genome resequencing data obtained from annual Mimulus guttatus individuals collected across a rugged environment landscape. We found substantial variation in how genetic differentiation is structured both within and between chromosomes, although, in contrast to other studies, known inversion polymorphisms appear to serve only minor roles in this heterogeneity. In addition, much of the genome can be clustered into eight among-population genetic differentiation patterns, but only two of these clusters are particularly consistent with patterns of isolation by distance. By performing genotype-environment association analysis, we also identified genomic intervals where local adaptation to specific climate factors has accentuated genetic differentiation among populations, and candidate genes in these windows indicate climate adaptation may proceed through changes affecting specialized metabolism, drought resistance, and development. Finally, by integrating our findings with previous studies, we show that multiple aspects of plant reproductive biology may be common targets of balancing selection and that variants historically involved in climate adaptation among populations have probably also fuelled rapid adaptation to microgeographic environmental variation within sites.

Plantago spp. as Models for Studying the Ecology and Evolution of Species Interactions across Environmental Gradients

AbstractA central challenge in ecology and evolutionary biology is to understand how variation in abiotic and biotic factors combine to shape the distribution, abundance, and diversity of focal species. Environmental gradients, whether natural (e.g., latitude, elevation, ocean proximity) or anthropogenic (e.g., land-use intensity, urbanization), provide compelling settings for addressing this challenge. However, not all organisms are amenable to the observational and experimental approaches required for untangling the factors that structure species along gradients. Here we highlight herbaceous plants in the genus Plantago as models for studying the ecology and evolution of species interactions along abiotic gradients. Plantago lanceolata and P. major are native to Europe and Asia but distributed globally, and they are established models for studying population ecology and interactions with herbivores, pathogens, and soil microbes. Studying restricted range congeners in comparison with those cosmopolitan species can provide insight into abiotic and biotic determinants of range size and population structure. We highlight one such species, P. rugelii, which is endemic to eastern North America. We give an overview of the literature on these focal Plantago species and explain why they are logical candidates for studies of species interactions across environmental gradients. Finally, we emphasize collaborative and community science approaches that can facilitate such research and note the amenability of Plantago for authentic research projects in science education.

Regional differences in rapid evolution during severe drought

Climate change is increasing drought intensity, threatening biodiversity. Rapid evolution of drought adaptations might be required for population persistence, particularly in rear-edge populations that may already be closer to physiological limits. Resurrection studies are a useful tool to assess adaptation to climate change, yet these studies rarely encompass the geographic range of a species. Here, we sampled 11 populations of scarlet monkeyflower (Mimulus cardinalis), collecting seeds across the plants' northern, central, and southern range to track trait evolution from the lowest to the greatest moisture anomaly over a 7-year period. We grew families generated from these populations across well-watered and terminal drought treatments in a greenhouse and quantified five traits associated with dehydration escape and avoidance. When considering pre-drought to peak-drought phenotypes, we find that later date of flowering evolved across the range of M. cardinalis, suggesting a shift away from dehydration escape. Instead, traits consistent with dehydration avoidance evolved, with smaller and/or thicker leaves evolving in central and southern regions. The southern region also saw a loss of plasticity in these leaf traits by the peak of the drought, whereas flowering time remained plastic across all regions. This observed shift in traits from escape to avoidance occurred only in certain regions, revealing the importance of geographic context when examining adaptations to climate change.

Evolutionary history and ecology shape the diversity and abundance of phytochemical arsenals across monkeyflowers

We examine the extent to which phylogenetic effects and ecology are associated with macroevolutionary patterns of phytochemical defence production across the Mimulus phylogeny. We grew plants from 21 species representing the five major sections of the Mimulus phylogeny in a common garden to assess how the arsenals (NMDS groupings) and abundances (concentrations) of a phytochemical defence, phenylpropanoid glycosides (PPGs), vary across the phylogeny. Very few PPGs are widespread across the genus, but many are common to multiple sections of the genus. Phytochemical arsenals cluster among sections in an NMDS and are not associated with total concentration of PPGs. There is a strong phylogenetic signal for phytochemical arsenal composition across the Mimulus genus, whereas ecological variables such as growing season length, latitude, and elevation do not significantly influence arsenal. In contrast, there is little phylogenetic signal for total PPG concentration, and this trait is significantly influenced by several ecological factors. Phytochemical arsenals and abundances are influenced by plant life history form. Both phylogenetic effects and ecology are related to phytochemical patterns across species, albeit in different ways. The independence of phytochemical defence concentrations from arsenal compositions indicates that these aspects of defence may continue to evolve independently of one another.

Applying a Chemogeographic Strategy for Natural Product Discovery from the Marine Cyanobacterium Moorena bouillonii

The tropical marine cyanobacterium Moorena bouillonii occupies a large geographic range across the Indian and Western Tropical Pacific Oceans and is a prolific producer of structurally unique and biologically active natural products. An ensemble of computational approaches, including the creation of the ORCA (Objective Relational Comparative Analysis) pipeline for flexible MS¹ feature detection and multivariate analyses, were used to analyze various M. bouillonii samples. The observed chemogeographic patterns suggested the production of regionally specific natural products by M. bouillonii. Analyzing the drivers of these chemogeographic patterns allowed for the identification, targeted isolation, and structure elucidation of a regionally specific natural product, doscadenamide A (1). Analyses of MS² fragmentation patterns further revealed this natural product to be part of an extensive family of herein annotated, proposed natural structural analogs (doscadenamides B–J, 2–10); the ensemble of structures reflect a combinatorial biosynthesis using nonribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) components. Compound 1 displayed synergistic in vitro cancer cell cytotoxicity when administered with lipopolysaccharide (LPS). These discoveries illustrate the utility in leveraging chemogeographic patterns for prioritizing natural product discovery efforts.

The Genetic Architecture of Plant Defense Trade-offs in a Common Monkeyflower

Determining how adaptive combinations of traits arose requires understanding the prevalence and scope of genetic constraints. Frequently observed phenotypic correlations between plant growth, defenses, and/or reproductive timing have led researchers to suggest that pleiotropy or strong genetic linkage between variants affecting independent traits is pervasive. Alternatively, these correlations could arise via independent mutations in different genes for each trait and extensive correlational selection. Here we evaluate these alternatives by conducting a quantitative trait loci (QTL) mapping experiment involving a cross between 2 populations of common monkeyflower (Mimulus guttatus) that differ in growth rate as well as total concentration and arsenal composition of plant defense compounds, phenylpropanoid glycosides (PPGs). We find no evidence that pleiotropy underlies correlations between defense and growth rate. However, there is a strong genetic correlation between levels of total PPGs and flowering time that is largely attributable to a single shared QTL. While this result suggests a role for pleiotropy/close linkage, several other QTLs also contribute to variation in total PPGs. Additionally, divergent PPG arsenals are influenced by a number of smaller-effect QTLs that each underlie variation in 1 or 2 PPGs. This result indicates that chemical defense arsenals can be finely adapted to biotic environments despite sharing a common biochemical precursor. Together, our results show correlations between defense and life-history traits are influenced by pleiotropy or genetic linkage, but genetic constraints may have limited impact on future evolutionary responses, as a substantial proportion of variation in each trait is controlled by independent loci.

The improved resistance against gypsy moth in Larix olgensis seedlings exposed to Cd stress association with elemental and chemical defenses

BACKGROUND

Cadmium (Cd), as an environmental pollutant, can endanger various biological and chemical characteristics of plants in multiple aspects. In this study, the effects of Cd contamination or exposure for 30 and 60 days at 1, 2, 4 mg kg^-1 concentrations on the resistance of Larix olgensis seedlings to the gypsy moth (Lymantria dispar) larvae were investigated.

RESULTS

Our results showed that Cd stress did not significantly affect the growth and biomass parameters of the larch seedlings, which might be attributed to the scavenging mechanism of reactive oxygen species (e.g. superoxide dismutase and peroxidase). Regarding the phytochemical defense, we found that Cd stress significantly changed the contents or activities of protease inhibitors (such as trypsin and chymotrypsin inhibitors) and secondary metabolites (tannins and phenolic acids) in L. olgensis seedling needles; however, their response trends varied with Cd exposure concentrations with a significant increase at low concentrations and a significant decrease at high concentrations. Moreover, both chemical and elemental defenses contributed to the resistance of L. olgensis seedlings to the gypsy moth larvae, and their synergistic effects (between toxic elements and organic metabolites) could provide an overall improved defense of L. olgensis seedlings even at low concentrations of single components, resulting in a detrimental effect on the growth of gypsy moth larvae.

CONCLUSION

These findings call for an urgent need to adjust and optimize pest control strategies in heavy metal polluted areas based on the effects of heavy metal stress on woody plant resistance to pest insects. © 2019 Society of Chemical Industry.

Contrasting environmental factors drive local adaptation at opposite ends of an environmental gradient in the yellow monkeyflower (Mimulus guttatus)

PREMISE

Identifying the environmental factors responsible for natural selection across different habitats is crucial for understanding the process of local adaptation in plants. Despite its importance, few studies have successfully isolated the environmental factors driving local adaptation in nature. In this study, we evaluated the agents of selection responsible for local adaptation of the monkeyflower Mimulus guttatus to California's coastal and inland habitats.

METHODS

We implemented a manipulative reciprocal transplant experiment at coastal and inland sites, where we excluded aboveground stressors in an effort to elucidate their role in the evolution of local adaptation.

RESULTS

Excluding aboveground stressors, most likely a combination of salt spray and herbivory, completely rescued inland annual plant fitness when transplanted to coastal habitat. The exclosures in inland habitat provided a benefit to the performance of coastal perennial plants. However, the exclosures are unlikely to provide much fitness benefit to the coastal plants at the inland site because of their general inability to flower in time to escape from the summer drought.

CONCLUSIONS

Our study demonstrates that a distinct set of selective agents (aboveground vs. belowground) are responsible for local adaptation at opposite ends of an environmental gradient.

Environmental heterogeneity generates intrapopulation variation in life-history traits in an annual plant

Environmental variation affects a plant's life cycle by influencing the timing of germination and flowering, and the duration of the growing season. Yet we know little information about how environmental heterogeneity generates variation in germination schedules and the consequences for growth and fecundity through genetic and plastic responses. We use an annual population of Mimulus guttatus in which, in nature, seeds germinate in both fall and spring. We investigate whether there is a genetic basis to the timing of germination, the effect of germination timing on fecundity, and if growth and flowering respond plastically to compensate for different season lengths. Using sibling families grown in simulated seasonal conditions, we find that families do not differ in their propensity to germinate between seasons. However, the germination season affects subsequent growth and flowering time, with significant genotype-by-environment interactions (G × E). Most G × E is due to unequal variance between seasons, because the spring cohort harbours little genetic variance. Despite their different season lengths, the cohorts do not differ in flower number (fecundity). Heterogeneous environments with unpredictable risks may maintain promiscuous germination, which then affects flowering time. Therefore, if selection at particular life stages changes with climate change, there may be consequences for the entire life cycle.

A Molecular View of Plant Local Adaptation: Incorporating Stress-Response Networks

Ecological specialization in plants occurs primarily through local adaptation to different environments. Local adaptation is widely thought to result in costly fitness trade-offs that result in maladaptation to alternative environments. However, recent studies suggest that such trade-offs are not universal. Further, there is currently a limited understanding of the molecular mechanisms responsible for fitness trade-offs associated with adaptation. Here, we review the literature on stress responses in plants to identify potential mechanisms underlying local adaptation and ecological specialization. We focus on drought, high and low temperature, flooding, herbivore, and pathogen stresses. We then synthesize our findings with recent advances in the local adaptation and plant molecular biology literature. In the process, we identify mechanisms that could cause fitness trade-offs and outline scenarios where trade-offs are not a necessary consequence of adaptation. Future studies should aim to explicitly integrate molecular mechanisms into studies of local adaptation.

Comparative transcriptome analyses revealed differential strategies of roots and leaves from methyl jasmonate treatment Baphicacanthus cusia (Nees) Bremek and differentially expressed genes involved in tryptophan biosynthesis

Baphicacanthus cusia (Nees) Bremek (B. cusia) is an effective herb for the treatment of acute promyelocytic leukemia and psoriasis in traditional Chinese medicine. Methyl jasmonate (MeJA) is a well-known signaling phytohormone that triggers gene expression in secondary metabolism. Currently, MeJA-mediated biosynthesis of indigo and indirubin in B. cusia is not well understood. In this study, we analyzed the content of indigo and indirubin in leaf and root tissues of B. cusia with high-performance liquid chromatography and measured photosynthetic characteristics of leaves treated by MeJA using FluorCam6 Fluorometer and chlorophyll fluorescence using the portable photosynthesis system CIRAS-2. We performed de novo RNA-seq of B. cusia leaf and root transcriptional profiles to investigate differentially expressed genes (DEGs) in response to exogenous MeJA application. The amount of indigo in MeJA-treated leaves were higher than that in controled leaves (p = 0.004), and the amounts of indigo in treated roots was higher than that in controlled roots (p = 0.048); Chlorophyll fluorescence of leaves treated with MeJA were significantly decreased. Leaves treated with MeJA showed lower photosynthetic rate compared to the control in the absence of MeJA. Functional annotation of DEGs showed the DEGs related to growth and development processes were down-regulated in the treated leaves, while most of the unigenes involved in the defense response were up-regulated in treated roots. This coincided with the effects of MeJA on photosynthetic characteristics and chlorophyll fluorescence. The qRT-PCR results showed that MeJA appears to down-regulate the gene expression of tryptophan synthase β-subunits (trpA-β) in leaves but increased the gene expression of anthranilate synthase (trp 3) in roots responsible for increased indigo content. The results showed that MeJA suppressed leaf photosynthesis for B. cusia and this growth-defense trade-off may contribute to the improved adaptability of B. cusia in changing environments.

Population Variation, Environmental Gradients, and the Evolutionary Ecology of Plant Defense against Herbivory

A central tenet of plant defense theory is that adaptation to the abiotic environment sets the template for defense strategies, imposing a trade-off between plant growth and defense. Yet this trade-off, commonly found among species occupying divergent resource environments, may not occur across populations of single species. We hypothesized that more favorable climates and higher levels of herbivory would lead to increases in growth and defense across plant populations. We evaluated whether plant growth and defense traits covaried across 18 populations of showy milkweed (Asclepias speciosa) inhabiting an east-west climate gradient spanning 25° of longitude. A suite of traits impacting defense (e.g., latex, cardenolides), growth (e.g., size), or both (e.g., specific leaf area [SLA], trichomes) were measured in natural populations and in a common garden, allowing us to evaluate plastic and genetically based variation in these traits. In natural populations, herbivore pressure increased toward warmer sites with longer growing seasons. Growth and defense traits showed strong clinal patterns and were positively correlated. In a common garden, clines with climatic origin were recapitulated only for defense traits. Correlations between growth and defense traits were also weaker and more negative in the common garden than in the natural populations. Thus, our data suggest that climatically favorable sites likely facilitate the evolution of greater defense at minimal costs to growth, likely because of increased resource acquisition.

Assortment of Flowering Time and Immunity Alleles in Natural Arabidopsis thaliana Populations Suggests Immunity and Vegetative Lifespan Strategies Coevolve

The selective impact of pathogen epidemics on host defenses can be strong but remains transient. By contrast, life-history shifts can durably and continuously modify the balance between costs and benefits of immunity, which arbitrates the evolution of host defenses. Their impact on the evolutionary dynamics of host immunity, however, has seldom been documented. Optimal investment into immunity is expected to decrease with shortening lifespan, because a shorter life decreases the probability to encounter pathogens or enemies. Here, we document that in natural populations of Arabidopsis thaliana, the expression levels of immunity genes correlate positively with flowering time, which in annual species is a proxy for lifespan. Using a novel genetic strategy based on bulk-segregants, we partitioned flowering time-dependent from -independent immunity genes and could demonstrate that this positive covariation can be genetically separated. It is therefore not explained by the pleiotropic action of some major regulatory genes controlling both immunity and lifespan. Moreover, we find that immunity genes containing variants reported to impact fitness in natural field conditions are among the genes whose expression covaries most strongly with flowering time. Taken together, these analyses reveal that natural selection has likely assorted alleles promoting lower expression of immunity genes with alleles that decrease the duration of vegetative lifespan in A. thaliana and vice versa. This is the first study documenting a pattern of variation consistent with the impact that selection on flowering time is predicted to have on diversity in host immunity.

Competition drives trait evolution and character displacement between Mimulus species along an environmental gradient

Closely related species may evolve to coexist stably in sympatry through niche differentiation driven by in situ competition, a process termed character displacement. Alternatively, past evolution in allopatry may have already sufficiently reduced niche overlap to permit establishment in sympatry, a process called ecological sorting. The relative importance of each process to niche differentiation is contentious even though they are not mutually exclusive and are both mediated via multivariate trait evolution. We explore how competition has impacted niche differentiation in two monkeyflowers, Mimulus alsinoides and M. guttatus, which often co-occur. Through field observations, common gardens, and competition experiments, we demonstrate that M. alsinoides is restricted to marginal habitats in sympatry and that the impacts of character displacement on niche differentiation are complex. Competition with M. guttatus alters selection gradients and has favored taller M. alsinoides with earlier seasonal flowering at low elevation and floral shape divergence at high elevation. However, no trait exhibits the pattern typically associated with character displacement, higher divergence between species in sympatry than allopatry. Thus, although character displacement was unlikely the process driving initial divergence along niche axes necessary for coexistence, we conclude that competition in sympatry has likely driven trait evolution along additional niche axes.