Adaptive generalization in pollination systems: Hawkmoths increase fitness to long-tubed flowers, but secondary pollinators remain important.

Selection on floral traits by animal pollinators is important in the evolution of flowering plants, yet whether floral divergence requires specialized pollination remains uncertain. Longer floral tubes, a trait associated with long-tongued pollinators, can also exclude other pollinators from accessing rewards, a potential mechanism for specialization. Across most of its range, Castilleja sessiliflora displays much longer corollas than most Castilleja species, though tube length varies geographically and correlates partially with hawkmoth visitation. To assess whether long corolla tubes reflect adaptation to hawkmoth pollinators, we performed a day/night pollinator exclusion experiment in nine natural populations that varied in corolla length across the range of C. sessiliflora and short-tubed members of the parapatric C. purpurea complex. We compared the fitness contributions of nocturnal and diurnal visitors, revealing that long-tubed populations visited predominantly by hawkmoths experienced greater fruit set at night, in contrast with short-tubed populations or those visited mainly by diurnal pollinators. Next, leveraging a range-wide multiyear dataset of pollinator visitation to these species, we identify that hawkmoth visitation is associated with increased fitness in long-tubed populations overall, and that long tubes are associated with less diverse visitor assemblages. Thus, long corollas represent an adaptation to hawkmoth pollination at the exclusion of diverse pollinators. Nonetheless, while hawkmoths were scarce in the northern range, secondary diurnal pollinators contributed to fruit set across the range, providing reproductive assurance despite possible trait mismatch. This study illustrates adaptive generalization in pollination systems and that floral divergence may proceed along a continuum of generalized and specialized pollinator interactions.

A phylogeny of the evening primrose family (Onagraceae) using a target enrichment approach with 303 nuclear loci.

BACKGROUND: The evening primrose family (Onagraceae) includes 664 species (803 taxa) with a center of diversity in the Americas, especially western North America. Ongoing research in Onagraceae includes exploring striking variation in floral morphology, scent composition, and breeding system, as well as the role of these traits in driving diversity among plants and their interacting pollinators and herbivores. However, these efforts are limited by the lack of a comprehensive, well-resolved phylogeny. Previous phylogenetic studies based on a few loci strongly support the monophyly of the family and the sister relationship of the two largest tribes but fail to resolve several key relationships. RESULTS: We used a target enrichment approach to reconstruct the phylogeny of Onagraceae using 303 highly conserved, low-copy nuclear loci. We present a phylogeny for Onagraceae with 169 individuals representing 152 taxa sampled across the family, including extensive sampling within the largest tribe, Onagreae. Deep splits within the family are strongly supported, whereas relationships among closely related genera and species are characterized by extensive conflict among individual gene trees. CONCLUSIONS: This phylogenetic resource will augment current research projects focused throughout the family in genomics, ecology, coevolutionary dynamics, biogeography, and the evolution of characters driving diversification in the family.

Hawkmoth and bee pollinators impact pollen dispersal at the landscape but not local scales in two species of Oenothera.

PREMISE: Animal pollinators play an important role in pollen dispersal. Here, we assessed differences in pollen and seed dispersal and the role of pollinator functional groups with different foraging behaviors in generating patterns of genetic diversity over similar geographic ranges for two closely related taxa. We focused on two members of Oenothera section Calylophus (Onagraceae) that co-occur on gypsum outcrops throughout the northern Chihuahuan Desert but differ in floral phenotype and primary pollinator: Oenothera gayleana (bee) and O. hartwegii subsp. filifolia (hawkmoth). METHODS: We measured breeding system and floral traits and studied gene flow and population differentiation at the local (<13 km; four populations) and landscape (60-440 km; five populations) scales using 10-11 nuclear (pollen dispersal) and three plastid (seed dispersal) microsatellite markers. RESULTS: Both taxa were self-incompatible and floral traits were consistent with expectations for different pollinators. Seed and pollen dispersal patterns were distinctly different for both species. We found no evidence of genetic structure at the local scale but did at the landscape scale; O. gayleana showed greater differentiation and significant isolation by distance than in O. hartwegii subsp. filifolia. The plastid data were consistent with gravity dispersal of seeds and suggest that pollen dispersal is the principal driver of genetic structure in both species. CONCLUSIONS: We demonstrated that pollinator functional groups can impact genetic differentiation in different and predictable ways. Hawkmoths, with larger foraging distances, can maintain gene flow across greater spatial scales than bees.

Target Enrichment and Extensive Population Sampling Help Untangle the Recent, Rapid Radiation of Oenothera Sect. Calylophus.

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking noncoding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic, and morphological diversity in the group. We found that the combination of exons and flanking noncoding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation. [Gypsum endemism; Oenothera sect. Calylophus; Onagraceae; phylogenomics; pollinator shift; recent radiation; target enrichment.].

Mutualism has its limits: consequences of asymmetric interactions between a well-defended plant and its herbivorous pollinator.

Concern for pollinator health often focuses on social bees and their agricultural importance at the expense of other pollinators and their ecosystem services. When pollinating herbivores use the same plants as nectar sources and larval hosts, ecological conflicts emerge for both parties, as the pollinator's services are mitigated by herbivory and its larvae are harmed by plant defences. We tracked individual-level metrics of pollinator health-growth, survivorship, fecundity-across the life cycle of a pollinating herbivore, the common hawkmoth, Hyles lineata, interacting with a rare plant, Oenothera harringtonii, that is polymorphic for the common floral volatile (R)-(-)-linalool. Linalool had no impact on floral attraction, but its experimental addition suppressed oviposition on plants lacking linalool. Plants showed robust resistance against herbivory from leaf-disc to whole-plant scales, through poor larval growth and survivorship. Higher larval performance on other Oenothera species indicates that constitutive herbivore resistance by O. harringtonii is not a genus-wide trait. Leaf volatiles differed among populations of O. harringtonii but were not induced by larval herbivory. Similarly, elagitannins and other phenolics varied among plant tissues but were not herbivore-induced. Our findings highlight asymmetric plant-pollinator interactions and the importance of third parties, including alternative larval host plants, in maintaining pollinator health. This article is part of the theme issue 'Natural processes influencing pollinator health: from chemistry to landscapes'.

Evolution of selfing syndrome and its influence on genetic diversity and inbreeding: A range-wide study in Oenothera primiveris.

PREMISE: To avoid inbreeding depression, plants have evolved diverse breeding systems to favor outcrossing, such as self-incompatibility. However, changes in biotic and abiotic conditions can result in selective pressures that lead to a breakdown in self-incompatibility. The shift to increased selfing is commonly associated with reduced floral features, lower attractiveness to pollinators, and increased inbreeding. We tested the hypothesis that the loss of self-incompatibility, a shift to self-fertilization (autogamy), and concomitant evolution of the selfing syndrome (reduction in floral traits associated with cross-fertilization) will lead to increased inbreeding and population differentiation in Oenothera primiveris. Across its range, this species exhibits a shift in its breeding system and floral traits from a self-incompatible population with large flowers to self-compatible populations with smaller flowers. METHODS: We conducted a breeding system assessment, evaluated floral traits in the field and under controlled conditions, and measured population genetic parameters using RADseq data. RESULTS: Our results reveal a bimodal transition to the selfing syndrome from the west to the east of the range of O. primiveris. This shift includes variation in the breeding system and the mating system, a reduction in floral traits (flower diameter, herkogamy, and scent production), a shift to greater autogamy, reduced genetic diversity, and increased inbreeding. CONCLUSIONS: The observed variation highlights the importance of range-wide studies to understand breeding system variation and the evolution of the selfing syndrome within populations and species.

Extensive population-level sampling reveals clinal variation in (R)-(-)-linalool produced by the flowers of an endemic evening primrose, Oenothera harringtonii.

The study of floral trait diversity has a long history due to its role in angiosperm diversification. While many studies have focused on visual traits including morphology and color, few have included floral scent despite its importance in pollination. Of the studies that have included floral scent, sampling has been limited and rarely explores variation at the population level. We studied geographic variation in the flowers of Oenothera harringtonii, a rare plant endemic to a vulnerable shortgrass prairie habitat, whose population structure and conservation status are well studied. The self-incompatible flowers of O. harringtonii open at dusk, produce nectar and a strong fragrance, and are pollinated by hawkmoths. We collected floral trait (morphology, scent chemistry and emission rates) data from 650 individuals from 19 wild populations to survey floral variation across the entire range of this species. Similarly, we collected floral data from 49 individuals grown in a greenhouse common garden, to assess whether variation observed in the field is consistent when environment factors (temperature, watering regime, soil) are standardized. We identified 35 floral volatiles representing 5 biosynthetic classes. Population differentiation was stronger for floral scent chemistry than floral morphology. (R)-(-)-linalool was the most important floral trait differentiating populations, exhibiting clinal variation across the distribution of O. harringtonii without any correlated shifts in floral morphology. Populations in the north and west produced (R)-(-)-linalool consistently, those in the east and south largely lacked it, and populations at the center of the distribution were polymorphic. Floral scent emissions in wild populations varied across four years but chemical composition was largely consistent over time. Similarly, volatile emission rates and chemical composition in greenhouse-grown plants were consistent with those of wild populations of origin. Our data set, which represents the most extensive population-level survey of floral scent to date, indicates that such sampling may be needed to capture potentially adaptive geographic variation in wild populations.

Differential gene expression associated with a floral scent polymorphism in the evening primrose Oenothera harringtonii (Onagraceae).

BACKGROUND: Plant volatiles play an important role in both plant-pollinator and plant-herbivore interactions. Intraspecific polymorphisms in volatile production are ubiquitous, but studies that explore underlying differential gene expression are rare. Oenothera harringtonii populations are polymorphic in floral emission of the monoterpene (R)-(-)-linalool; some plants emit (R)-(-)-linalool (linalool+ plants) while others do not (linalool- plants). However, the genes associated with differential production of this floral volatile in Oenothera are unknown. We used RNA-Seq to broadly characterize differential gene expression involved in (R)-(-)-linalool biosynthesis. To identify genes that may be associated with the polymorphism for this trait, we used RNA-Seq to compare gene expression in six different Oenothera harringtonii tissues from each of three linalool+ and linalool- plants. RESULTS: Three clusters of differentially expressed genes were enriched for terpene synthase activity: two were characterized by tissue-specific upregulation and one by upregulation only in plants with flowers that produce (R)-(-)-linalool. A molecular phylogeny of all terpene synthases identified two putative (R)-(-)-linalool synthase transcripts in Oenothera harringtonii, a single allele of which is found exclusively in linalool+ plants. CONCLUSIONS: By using a naturally occurring polymorphism and comparing different tissues, we were able to identify candidate genes putatively involved in the biosynthesis of (R)-(-)-linalool. Expression of these genes in linalool- plants, while low, suggests a regulatory polymorphism, rather than a population-specific loss-of-function allele. Additional terpene biosynthesis-related genes that are up-regulated in plants that emit (R)-(-)-linalool may be associated with herbivore defense, suggesting a potential economy of scale between plant reproduction and defense.

Incomplete reproductive isolation and low genetic differentiation despite floral divergence across varying geographic scales in Castilleja.

PREMISE: Divergence depends on the strength of selection and frequency of gene flow between taxa, while reproductive isolation relies on mating barriers and geographic distance. Less is known about how these processes interact at early stages of speciation. Here, we compared population-level differentiation in floral phenotype and genetic sequence variation among recently diverged Castilleja to explore patterns of diversification under different scenarios of reproductive isolation. METHODS: Using target enrichment enabled by the Angiosperms353 probe set, we assessed genetic distance among 50 populations of four Castilleja species. We investigated whether patterns of genetic divergence are explained by floral trait variation or geographic distance in two focal groups: the widespread C. sessiliflora and the more restricted C. purpurea species complex. RESULTS: We document that C. sessiliflora and the C. purpurea complex are characterized by high diversity in floral color across varying geographic scales. Despite phenotypic divergence, groups were not well supported in phylogenetic analyses, and little genetic differentiation was found across targeted Angiosperms353 loci. Nonetheless, a principal coordinate analysis of single nucleotide polymorphisms revealed differentiation within C. sessiliflora across floral morphs and geography and less differentiation among species of the C. purpurea complex. CONCLUSIONS: Patterns of genetic distance in C. sessiliflora suggest species cohesion maintained over long distances despite variation in floral traits. In the C. purpurea complex, divergence in floral color across narrow geographic clines may be driven by recent selection on floral color. These contrasting patterns of floral and genetic differentiation reveal that divergence can arise via multiple eco-evolutionary paths.

Phylogeny, host use, and diversification in the moth family Momphidae (Lepidoptera: Gelechioidea).

Insect herbivores and their hostplants constitute much of Earth's described biological diversity, but how these often-specialized associations diversify is not fully understood. We combined detailed hostplant data and comparative phylogenetic analyses of the lepidopteran family Momphidae to explore how shifts in the use of hostplant resources, not just hostplant taxon, contribute to the diversification of a phytophagous insect lineage. We inferred two phylogenetic hypotheses emphasizing relationships among species in the nominate genus, Mompha Hübner. A six-gene phylogeny was constructed with reared exemplars and collections from hostplants in the family Onagraceae from western and southwestern USA, and a cytochrome c oxidase subunit 1 (COI) phylogeny was inferred from collections and publicly available accessions in the Barcode of Life Data System. Species delimitation analyses combined with morphological data revealed ca. 56 undescribed species-level taxa, many of which are hostplant specialists on Onagraceae in the southwestern USA. Our phylogenetic reconstructions divided Momphidae into six major clades: 1) an Onagraceae flower- and fruit-boring clade, 2) a Melastomataceae-galling clade, 3) a leafmining clade A, 4) a leafmining clade B, 5) a Zapyrastra Meyrick clade, and 6) a monobasic lineage represented by Mompha eloisella (Clemens). Ancestral trait reconstructions using the COI phylogeny identified leafmining on Onagraceae as the ancestral state for Momphidae. Our study finds that shifts along three hostplant resource axes (plant taxon, plant tissue type, and larval feeding mode) have contributed to the evolutionary success and diversification of momphids.

Pollinator identity and spatial isolation influence multiple paternity in an annual plant.

The occurrence and extent of multiple paternity is an important component of variation in plant mating dynamics. However, links between pollinator activity and multiple paternity are generally lacking, especially for plant species that attract functionally diverse floral visitors. In this study, we separated the influence of two functionally distinct floral visitors (hawkmoths and solitary bees) and characterized their impacts on multiple paternity in a self-incompatible, annual forb, Oenothera harringtonii (Onagraceae). We also situated pollinator-mediated effects in a spatial context by linking variation in multiple paternity to variation in plant spatial isolation. We documented pronounced differences in the number of paternal sires as function of pollinator identity: on average, the primary pollinator (hawkmoths) facilitated mating with nearly twice as many pollen donors relative to the secondary pollinator (solitary bees). This effect was consistent for both isolated and nonisolated individuals, but spatial isolation imposed pronounced reductions on multiple paternity regardless of pollinator identity. Considering that pollinator abundance and pollen dispersal distance did not vary significantly with pollinator identity, we attribute variation in realized mating dynamics primarily to differences in pollinator morphology and behaviour as opposed to pollinator abundance or mating incompatibility arising from underlying spatial genetic structure. Our findings demonstrate that functionally distinct pollinators can have strongly divergent effects on polyandry in plants and further suggest that both pollinator identity and spatial heterogeneity have important roles in plant mating dynamics.

Herbivory as an important selective force in the evolution of floral traits and pollinator shifts.

Floral trait evolution is frequently attributed to pollinator-mediated selection but herbivores can play a key role in shaping plant reproductive biology. Here we examine the role of florivores in driving floral trait evolution and pollinator shifts in a recently radiated clade of flowering plants, Oenothera sect. Calylophus We compare florivory by a specialist, internal feeder, Mompha, on closely related hawkmoth- and bee-pollinated species and document variation in damage based on floral traits within sites, species and among species. Our results show that flowers with longer floral tubes and decreased floral flare have increased Mompha damage. Bee-pollinated flowers, which have substantially smaller floral tubes, experience on average 13% less Mompha florivory than do hawkmoth-pollinated flowers. The positive association between tube length and Mompha damage is evident even within sites of some species, suggesting that Mompha can drive trait differentiation at microevolutionary scales. Given that there are at least two independent shifts from hawkmoth to bee pollination in this clade, florivore-mediated selection on floral traits may have played an important role in facilitating morphological changes associated with transitions from hawkmoth to bee pollination.

Land-use change has no detectable effect on reproduction of a disturbance-adapted, hawkmoth-pollinated plant species.

PREMISE OF THE STUDY: Land-use change is cited as a primary driver of global biodiversity loss, with myriad consequences for species, populations, and ecosystems. However, few studies have examined its impact on species interactions, particularly pollination. Furthermore, when the effects of land-use change on pollination have been studied, the focus has largely been on species pollinated by diurnal pollinators, namely, bees and butterflies. Here, we focus on Oenothera harringtonii, a night-flowering, disturbance-adapted species that has experienced a range-wide gradient of land-use change. We tested the hypothesis that the negative impacts of land-use change are mitigated by long-distance pollination. METHODS: Our study included both temporal (4 yr) and spatial (19 populations range-wide, and 1, 2, and 5 km from the population center) data, providing a comprehensive understanding of the role of land-use change on pollination biology and reproduction. KEY RESULTS: We first confirmed that O. harringtonii is self-incompatible and reliant on pollinators for reproduction. We then showed that hawkmoths (primarily Hyles lineata) are highly reliable and effective pollinators in both space and time. Unlike other studies, we did not detect an effect of population size, increased isolation, or a reduction in suitable habitat in areas with evidence of land-use change on pollination (visitation, pollen removal and deposition). Furthermore, the proportion of suitable habitat and other fragmentation metrics examined were not associated with population size or density in this plant species. CONCLUSIONS: We conclude that nocturnal pollination of Oenothera harringtonii via hawkmoths is robust to the negative impacts of land-use change.

Microsatellites for Oenothera gayleana and O. hartwegii subsp. filifolia (Onagraceae), and their utility in section Calylophus.

PREMISE OF THE STUDY: Eleven nuclear and four plastid microsatellite markers were screened for two gypsum endemic species, Oenothera gayleana and O. hartwegii subsp. filifolia, and tested for cross-amplification in the remaining 11 taxa within Oenothera sect. Calylophus (Onagraceae). METHODS AND RESULTS: Microsatellite markers were tested in two to three populations spanning the ranges of both O. gayleana and O. hartwegii subsp. filifolia. The nuclear microsatellite loci consisted of both di- and trinucleotide repeats with one to 17 alleles per population. Several loci showed significant deviation from Hardy-Weinberg equilibrium, which may be evidence of chromosomal rings. The plastid microsatellite markers identified one to seven haplotypes per population. The transferability of these markers was confirmed in all 11 taxa within Oenothera sect. Calylophus. CONCLUSIONS: The microsatellite loci characterized here are the first developed and tested in Oenothera sect. Calylophus. These markers will be used to assess whether pollinator foraging distance influences population genetic parameters in predictable ways.

Local topography shapes fine-scale spatial genetic structure in the Arkansas Valley evening primrose, Oenothera harringtonii (Onagraceae).

Identifying factors that shape the spatial distribution of genetic variation is crucial to understanding many population- and landscape-level processes. In this study, we explore fine-scale spatial genetic structure in Oenothera harringtonii (Onagraceae), an insect-pollinated, gravity-dispersed herb endemic to the grasslands of south-central and southeastern Colorado, USA. We genotyped 315 individuals with 11 microsatellite markers and utilized a combination of spatial autocorrelation analyses and landscape genetic models to relate life history traits and landscape features to dispersal processes. Spatial genetic structure was consistent with theoretical expectations of isolation by distance, but this pattern was weak (Sp = 0.00374). Anisotropic analyses indicated that spatial genetic structure was markedly directional, in this case consistent with increased dispersal along prominent slopes. Landscape genetic models subsequently confirmed that spatial genetic variation was significantly influenced by local topographic heterogeneity, specifically that geographic distance, elevation and aspect were important predictors of spatial genetic structure. Among these variables, geographic distance was ~68% more important than elevation in describing spatial genetic variation, and elevation was ~42% more important than aspect after removing the effect of geographic distance. From these results, we infer a mechanism of hydrochorous seed dispersal along major drainages aided by seasonal monsoon rains. Our findings suggest that landscape features may shape microevolutionary processes at much finer spatial scales than typically considered, and stress the importance of considering how particular dispersal vectors are influenced by their environmental context.

Characterization of microsatellite loci in Castilleja sessiliflora and transferability to 24 Castilleja species (Orobanchaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed in the hemiparasitic perennial forb Castilleja sessiliflora to investigate patterns of gene flow and genetic diversity within and among populations. • METHODS AND RESULTS: Twelve polymorphic loci were identified in C. sessiliflora and tested on three populations (32 individuals each) sampled across the range of the species. The loci amplified di- and trinucleotide repeats with 3-14 alleles per locus. To assess cross-amplification, primer pairs were also tested on 24 additional Castilleja species that represent the morphological and geographic diversity of the genus. We provide reports of their effectiveness in all 25 taxa. • CONCLUSIONS: These results indicate the utility of these primers in C. sessiliflora for future studies of genetic structure and gene flow, as well as their widespread applicability in other members of the diverse and complex genus Castilleja.

Microsatellite primers in Oenothera harringtonii (Onagraceae), an annual endemic to the shortgrass prairie of Colorado.

PREMISE OF THE STUDY: Microsatellite markers were developed in the annual herb, Oenothera harringtonii, to investigate patterns of genetic diversity, gene flow, and parentage within and among populations of this Colorado endemic. METHODS AND RESULTS: Ten polymorphic loci were identified in O. harringtonii and tested in four populations sampled across the range of the species. These loci contained trinucleotide repeats with 7-29 alleles per locus. Nine of the 10 loci also amplified in O. caespitosa subsp. macroglottis, O. caespitosa subsp. marginata, and O. caespitosa subsp. navajoensis. In addition, we optimized three markers developed for O. biennis and provide reports of their effectiveness in all four taxa. CONCLUSIONS: These results indicate the utility of these markers in O. harringtonii for future studies of genetic structure, gene flow, and parentage as well as their applicability in other members of the O. caespitosa species complex.

Nitrogen deposition, competition and the decline of a regionally threatened legume, Desmodium cuspidatum.

Increased nitrogen (N) deposition, resulting from the combustion of fossil fuels, production of synthetic fertilizers, growth of N(2)-fixing crops and high-intensity agriculture, is one of the anthropogenic factors most likely to cause global biodiversity changes over the next century. This influence may be especially large in temperate zone forests, which are highly N limited and occur in regions with the highest levels of N deposition. Within these ecosystems, N(2)-fixing plants, including legumes, may be more sensitive to N deposition than other plant species. Though it has long been recognized that the competitive edge conferred by N(2)-fixation diminishes with increasing soil N availability, the conservation implications of increased N deposition on native N(2)-fixers have received less attention. We focus on Desmodium cuspidatum, which has experienced dramatic population losses in the last 30-40 years in the northeastern United States. We explore competition between this regionally threatened legume and a common non-N(2)-fixing neighbor, Solidago canadensis, across a gradient of N deposition. Our data show that increased N deposition may be detrimental to N(2)-fixers such as D. cuspidatum in two ways: (1) biomass accumulation in the non-N(2)-fixer, S. canadensis, responds more strongly to increasing N deposition, and (2) S. canadensis competes strongly for available mineral nitrogen and can assimilate N previously fixed by D. cuspidatum, resulting in D. cuspidatum relying more heavily on energetically expensive N(2)-fixation when grown with S. canadensis. N deposition may thus reduce or eliminate the competitive advantage of N(2)-fixing species growing in N-limited ecosystems.