SEED: A framework for integrating ecological stoichiometry and eco-evolutionary dynamics.

Characterising the extent and sources of intraspecific variation and their ecological consequences is a central challenge in the study of eco-evolutionary dynamics. Ecological stoichiometry, which uses elemental variation of organisms and their environment to understand ecosystem patterns and processes, can be a powerful framework for characterising eco-evolutionary dynamics. However, the current emphasis on the relative content of elements in the body (i.e. organismal stoichiometry) has constrained its application. Intraspecific variation in the rates at which elements are acquired, assimilated, allocated or lost is often greater than the variation in organismal stoichiometry. There is much to gain from studying these traits together as components of an 'elemental phenotype'. Furthermore, each of these traits can have distinct ecological effects that are underappreciated in the current literature. We propose a conceptual framework that explores how microevolutionary change in the elemental phenotype occurs, how its components interact with each other and with other traits, and how its changes can affect a wide range of ecological processes. We demonstrate how the framework can be used to generate novel hypotheses and outline pathways for future research that enhance our ability to explain, analyse and predict eco-evolutionary dynamics.

Life stage and taxonomy the most important factors determining vertebrate stoichiometry: A meta-analysis.

Whole-body elemental composition is a key trait for determining how organisms influence their ecosystems. Using mass-balance, ecological stoichiometry predicts that animals with higher concentrations of element X will selectively retain more X and will recycle less X in their waste than animals with lower X concentrations. These animals will also store high quantities of X during their lives and after their deaths (prior to full decomposition). Vertebrates may uniquely impact nutrient cycling because they store high quantities of phosphorus (P) in their bones. However, vertebrates have diverse body forms and invest variably in bone. Current analyses of vertebrate elemental content predominately evaluate fishes, typically neglecting other vertebrates and leaving much of the diversity unexplored. We performed a systematic review and identified 179 measurements of whole-body percent phosphorus (%P), percent nitrogen (%N), and N to P ratio (N:P) from 129 unique species of non-fish vertebrates (amphibians: 39 species; reptiles: 19 species; birds: 27 species; mammals: 46 species). We found that %P (mean: 1.94%; SD [standard deviation] = 0.77) and N:P (mean: 12.52) varied with taxonomy and life stage, while %N (mean: 10.51%; SD = 3.25) varied primarily with taxonomy. Habitat, diet, and size had small and inconsistent effects in different groups. Our study highlights two research gaps. Life stage, which is frequently neglected in stoichiometric studies, is an important factor determining vertebrate %P. Furthermore, amphibians dominate our dataset, while other vertebrate taxa are poorly represented in the current literature. Further research into these neglected vertebrate taxa is essential.

Phosphorus limitation does not drive loss of bony lateral plates in freshwater stickleback (Gasterosteus aculeatus).

Connecting the selective forces that drive the evolution of phenotypes to their underlying genotypes is key to understanding adaptation, but such connections are rarely tested experimentally. Threespine stickleback (Gasterosteus aculeatus) are a powerful model for such tests because genotypes that underlie putatively adaptive traits have been identified. For example, a regulatory mutation in the Ectodysplasin (Eda) gene causes a reduction in the number of bony armor plates, which occurs rapidly and repeatedly when marine sticklebacks invade freshwater. However, the source of selection on plate loss in freshwater is unknown. Here, we tested whether dietary reduction of phosphorus can account for selection on plate loss due to a growth advantage of low-plated fish in freshwater. We crossed marine fish heterozygous for the 16 kilobase freshwater Eda haplotype and compared the growth of offspring with different genotypes under contrasting levels of dietary phosphorus in both saltwater and freshwater. Eda genotype was not associated with growth differences in any treatment, or with mechanisms that could mitigate the impacts of phosphorus limitation, such as differential phosphorus deposition, phosphorus excretion, or intestine length. This study highlights the importance of experimentally testing the putative selective forces acting on phenotypes and their underlying genotypes in the wild.

Evaluating ecosystem effects of climate change on tropical island streams using high spatial and temporal resolution sampling regimes.

Climate change is expected to alter precipitation patterns worldwide, which will affect streamflow in riverine ecosystems. It is vital to understand the impacts of projected flow variations, especially in tropical regions where the effects of climate change are expected to be one of the earliest to emerge. Space-for-time substitutions have been successful at predicting effects of climate change in terrestrial systems by using a spatial gradient to mimic the projected temporal change. However, concerns have been raised that the spatial variability in these models might not reflect the temporal variability. We utilized a well-constrained rainfall gradient on Hawaii Island to determine (a) how predicted decreases in flow and increases in flow variability affect stream food resources and consumers and (b) if using a high temporal (monthly, four streams) or a high spatial (annual, eight streams) resolution sampling scheme would alter the results of a space-for-time substitution. Declines in benthic and suspended resource quantity (10- to 40-fold) and quality (shift from macrophyte to leaf litter dominated) contributed to 35-fold decreases in macroinvertebrate biomass with predicted changes in the magnitude and variability in the flow. Invertebrate composition switched from caddisflies and damselflies to taxa with faster turnover rates (mosquitoes, copepods). Changes in resource and consumer composition patterns were stronger with high temporal resolution sampling. However, trends and ranges of results did not differ between the two sampling regimes, indicating that a suitable, well-constrained spatial gradient is an appropriate tool for examining temporal change. Our study is the first to investigate resource to community wide effects of climate change on tropical streams on a spatial and temporal scale. We determined that predicted flow alterations would decrease stream resource and consumer quantity and quality, which can alter stream function, as well as biomass and habitat for freshwater, marine, and terrestrial consumers dependent on these resources.

Adaptation in temporally variable environments: stickleback armor in periodically breaching bar-built estuaries.

The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar-built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re-form within a few weeks or months, thereby re-isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within-population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure and then analysed how stickleback armor traits, the associated genes Eda and Pitx1 and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among-population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within-population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype and the presence of predators, thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within-population variation while reducing, but not eliminating, among-population variation driven by spatial environmental variation.

The influence of dietary and whole-body nutrient content on the excretion of a vertebrate consumer.

In many contexts, nutrient excretion by consumers can impact ecosystems by altering the availability of limiting nutrients. Variation in nutrient excretion can be predicted by mass balance models, most of which are premised on two key ideas: (1) consumers maintain fixed whole-body nutrient content (i.e., %N and %P), so-called fixed homeostasis; (2) if dietary nutrients are not matched to whole-body nutrients, excesses of any nutrient are released as excretion to maintain fixed homeostasis. Mass balance models thus predict that consumer excretion should be positively correlated with diet nutrients and negatively correlated with whole-body nutrients. Recent meta-analyses and field studies, however, have often failed to find these expected patterns, potentially because of a confounding influence-flexibility in whole-body nutrient content with diet quality (flexible homeostasis). Here, we explore the impact of flexible homeostasis on nutrient excretion by comparing the N and P excretion of four genetically diverged Trinidadian guppy (Poecilia reticulata) populations when reared on diets of variable P content. As predicted by mass balance, P excretion increased on the high-P diet, but, contrary to the notion of fixed homeostasis, guppy whole-body %P also increased on the high-P diet. While there was no overall correlation between excretion nutrients and whole-body nutrients, when the effect of diet on both whole-body and excretion nutrients was included, we detected the expected negative correlation between whole-body N:P and excretion N:P. This last result suggests that mass balance can predict excretion rates within species, but only if dietary effects on whole-body nutrient content are controlled. Flexible homeostasis can obscure patterns predicted by mass balance, creating an imperative to accurately capture an organism's diet quality in predicting its excretion rate.

Population variation in the trophic niche of the Trinidadian guppy from different predation regimes.

Population variation in trophic niche is widespread among organisms and is of increasing interest given its role in both speciation and adaptation to changing environments. Trinidadian guppies (Poecilia reticulata) inhabiting stream reaches with different predation regimes have rapidly evolved divergent life history traits. Here, we investigated the effects of both predation and resource availability on guppy trophic niches by evaluating their gut contents, resource standing stocks, and δ15N and δ13C stable isotopes across five streams during the wet season. We found that guppies from low predation (LP) sites had a consistently higher trophic position and proportion of invertebrates in their guts and assimilate less epilithon than guppies from high predation (HP) sites. Higher trophic position was also associated with lower benthic invertebrate availability. Our results suggest that LP guppies could be more efficient invertebrate consumers, possibly as an evolutionary response to greater intraspecific competition for higher quality food. This may be intensified by seasonality, as wet season conditions can alter resource availability, feeding rates, and the intensity of intraspecific competition. Understanding how guppy diets vary among communities is critical to elucidating the role of niche shifts in mediating the link between environmental change and the evolution of life histories.

Implications of guppy (Poecilia reticulata) life-history phenotype for mosquito control.

Guppies (Poecilia reticulata) are frequently introduced to both natural and artificial water bodies as a mosquito control. Laboratory studies have demonstrated that guppies can consume large numbers of larval mosquitoes. Our study investigates how intraspecific variability in guppy phenotype affects their importance as a mosquito biocontrol and how habitat conditions (natural ponds vs. water storage containers) may influence insect biomass and guppy feeding. Using a blocked experimental design, we established stream-side mesocosm ponds with half receiving gravel substrate to simulate pond-bottom habitat. To provide realistic diet choices and insect abundances, we allowed the mesocosms to colonize naturally with aquatic insect larvae for 1 month before introducing guppies. We tested two distinct guppy phenotypes (from high- and low-predation streams) alongside fish-free controls. After 1 month, we measured insect biomass in the mesocosms and examined guppy gut contents to document direct predation. While overall insect biomass was not significantly different across the three fish treatments, we observed a significant reduction in mosquito biomass in fish treatments compared to fish-free controls, as well as intraspecific differences in feeding. Overall insect biomass was significantly higher in mesocosms without gravel, while habitat condition had no effect on mosquito biomass. As guppy phenotype responds to changes in their environments, it is an important consideration for biocontrol policy to anticipate potential ecosystem effects. We close by relating our findings to other studies and by discussing the implications and potential risks of using guppies to control mosquitoes.

Stoichiometric traits of stickleback: Effects of genetic background, rearing environment, and ontogeny.

Phenotypes can both evolve in response to, and affect, ecosystem change, but few examples of diverging ecosystem-effect traits have been investigated. Bony armor traits of fish are good candidates for this because they evolve rapidly in some freshwater fish populations, and bone is phosphorus rich and likely to affect nutrient recycling in aquatic ecosystems. Here, we explore how ontogeny, rearing environment, and bone allocation among body parts affect the stoichiometric phenotype (i.e., stoichiometric composition of bodies and excretion) of threespine stickleback. We use two populations from distinct freshwater lineages with contrasting lateral plating phenotypes (full vs. low plating) and their hybrids, which are mostly fully plated. We found that ontogeny, rearing environment, and body condition were the most important predictors of organismal stoichiometry. Although elemental composition was similar between both populations and their hybrids, we found significant divergence in phosphorus allocation among body parts and in phosphorus excretion rates. Overall, body armor differences did not explain variation in whole body phosphorus, phosphorus allocation, or phosphorus excretion. Evolutionary divergence between these lineages in both allocation and excretion is likely to have important direct consequences for ecosystems, but may be mediated by evolution of multiple morphological or physiological traits beyond plating phenotype.

Fish introductions and light modulate food web fluxes in tropical streams: a whole-ecosystem experimental approach.

Decades of ecological study have demonstrated the importance of top-down and bottom-up controls on food webs, yet few studies within this context have quantified the magnitude of energy and material fluxes at the whole-ecosystem scale. We examined top-down and bottom-up effects on food web fluxes using a field experiment that manipulated the presence of a consumer, the Trinidadian guppy Poecilia reticulata, and the production of basal resources by thinning the riparian forest canopy to increase incident light. To gauge the effects of these reach-scale manipulations on food web fluxes, we used a nitrogen (15 N) stable isotope tracer to compare basal resource treatments (thinned canopy vs. control) and consumer treatments (guppy introduction vs. control). The thinned canopy stream had higher primary production than the natural canopy control, leading to increased N fluxes to invertebrates that feed on benthic biofilms (grazers), fine benthic organic matter (collector-gatherers), and organic particles suspended in the water column (filter feeders). Stream reaches with guppies also had higher primary productivity and higher N fluxes to grazers and filter feeders. In contrast, N fluxes to collector-gatherers were reduced in guppy introduction reaches relative to upstream controls. N fluxes to leaf-shredding invertebrates, predatory invertebrates, and the other fish species present (Hart's killifish, Anablepsoides hartii) did not differ across light or guppy treatments, suggesting that effects on detritus-based linkages and upper trophic levels were not as strong. Effect sizes of guppy and canopy treatments on N flux rates were similar for most taxa, though guppy effects were the strongest for filter feeding invertebrates while canopy effects were the strongest for collector-gatherer invertebrates. Combined, these results extend previous knowledge about top-down and bottom-up controls on ecosystems by providing experimental, reach-scale evidence that both pathways can act simultaneously and have equally strong influence on nutrient fluxes from inorganic pools through primary consumers.

Investment in boney defensive traits alters organismal stoichiometry and excretion in fish.

Understanding how trait diversification alters ecosystem processes is an important goal for ecological and evolutionary studies. Ecological stoichiometry provides a framework for predicting how traits affect ecosystem function. The growth rate hypothesis of ecological stoichiometry links growth and phosphorus (P) body composition in taxa where nucleic acids are a significant pool of body P. In vertebrates, however, most of the P is bound within bone, and organisms with boney structures can vary in terms of the relative contributions of bones to body composition. Threespine stickleback populations have substantial variation in boney armour plating. Shaped by natural selection, this variation provides a model system to study the links between evolution of bone content, elemental body composition, and P excretion. We measure carbon:nitrogen:P body composition from stickleback populations that vary in armour phenotype. We develop a mechanistic mass-balance model to explore factors affecting P excretion, and measure P excretion from two populations with contrasting armour phenotypes. Completely armoured morphs have higher body %P but excrete more P per unit body mass than other morphs. The model suggests that such differences are driven by phenotypic differences in P intake as well as body %P composition. Our results show that while investment in boney traits alters the elemental composition of vertebrate bodies, excretion rates depend on how acquisition and assimilation traits covary with boney trait investment. These results also provide a stoichiometric hypothesis to explain the repeated loss of boney armour in threespine sticklebacks upon colonizing freshwater ecosystems.

A test of the effects of timing of a pulsed resource subsidy on stream ecosystems.

Spatial resource subsidies can alter bottom-up and top-down forces of community regulation across ecosystem boundaries. Most subsidies are temporally variable, and recent theory has suggested that consumer-resource dynamics can be stabilized if the peak timing of a subsidy is desynchronized with that of prey productivity in the recipient ecosystem. However, magnitude of consumer responses per se could depend on the subsidy timing, which may be a critical component for community dynamics and ecosystem processes. The aim of this study was to test (i) whether a recipient consumer (cutthroat trout) responds differently to a resource subsidy occurring early in its growing season than to a subsidy occurring late in the season and, if this is the case, (ii) whether the timing-dependent consumer response has cascading effects on communities and ecosystem functions in streams. To test those hypotheses, we conducted a large-scale field experiment, in which we directly manipulated the timing of augmentation of the terrestrial invertebrates that enter stream (i.e. peak timing of June-August vs. August-October), keeping constant the total amounts of the invertebrates entered. We found large increases in the individual growth rate and population biomass of the cutthroat trout, in response to the early resource pulse, but not to the late pulse. This timing-dependent consumer response cascaded down to reduce benthic invertebrates and leaf breakdown rate, and increased water nutrient concentrations. Furthermore, the early resource pulse resulted in higher maturity rate of the cutthroat trout in the following spring, demonstrating the importance of the subsidy timing on long-term community dynamics via the consumer's numerical response. Our results emphasize the need to acknowledge timing-dependent consumer responses in understanding the effects of subsidies on communities and ecosystem processes. Elucidating the mechanisms by which consumers effectively exploit pulsed subsidies is an important avenue to better understand community dynamics in spatially coupled ecosystems.

Biodiversity and ecosystem risks arising from using guppies to control mosquitoes.

Deploying mosquito predators such as the guppy (Poecilia reticulata) into bodies of water where mosquitoes breed is a common strategy for limiting the spread of disease-carrying mosquitoes. Here, we draw on studies from epidemiology, conservation, ecology and evolution to show that the evidence for the effectiveness of guppies in controlling mosquitoes is weak, that the chances of accidental guppy introduction into local ecosystems are large, and that guppies can easily establish populations and damage these aquatic ecosystems. We highlight several knowledge and implementation gaps, and urge that this approach is either abandoned in favour of more effective strategies or that it is used much more rigorously. Controlling mosquitoes does not need to come at the expense of freshwater biodiversity.

Adaptive genetic variation mediates bottom-up and top-down control in an aquatic ecosystem.

Research in eco-evolutionary dynamics and community genetics has demonstrated that variation within a species can have strong impacts on associated communities and ecosystem processes. Yet, these studies have centred around individual focal species and at single trophic levels, ignoring the role of phenotypic variation in multiple taxa within an ecosystem. Given the ubiquitous nature of local adaptation, and thus intraspecific variation, we sought to understand how combinations of intraspecific variation in multiple species within an ecosystem impacts its ecology. Using two species that co-occur and demonstrate adaptation to their natal environments, black cottonwood (Populus trichocarpa) and three-spined stickleback (Gasterosteus aculeatus), we investigated the effects of intraspecific phenotypic variation on both top-down and bottom-up forces using a large-scale aquatic mesocosm experiment. Black cottonwood genotypes exhibit genetic variation in their productivity and consequently their leaf litter subsidies to the aquatic system, which mediates the strength of top-down effects from stickleback on prey abundances. Abundances of four common invertebrate prey species and available phosphorous, the most critically limiting nutrient in freshwater systems, are dictated by the interaction between genetic variation in cottonwood productivity and stickleback morphology. These interactive effects fit with ecological theory on the relationship between productivity and top-down control and are comparable in strength to the effects of predator addition. Our results illustrate that intraspecific variation, which can evolve rapidly, is an under-appreciated driver of community structure and ecosystem function, demonstrating that a multi-trophic perspective is essential to understanding the role of evolution in structuring ecological patterns.

Testing a 'genes-to-ecosystems' approach to understanding aquatic-terrestrial linkages.

A 'genes-to-ecosystems' approach has been proposed as a novel avenue for integrating the consequences of intraspecific genetic variation with the underlying genetic architecture of a species to shed light on the relationships among hierarchies of ecological organization (genes → individuals → communities → ecosystems). However, attempts to identify genes with major effect on the structure of communities and/or ecosystem processes have been limited and a comprehensive test of this approach has yet to emerge. Here, we present an interdisciplinary field study that integrated a common garden containing different genotypes of a dominant, riparian tree, Populus trichocarpa, and aquatic mesocosms to determine how intraspecific variation in leaf litter alters both terrestrial and aquatic communities and ecosystem functioning. Moreover, we incorporate data from extensive trait screening and genome-wide association studies estimating the heritability and genes associated with litter characteristics. We found that tree genotypes varied considerably in the quality and production of leaf litter, which contributed to variation in phytoplankton abundances, as well as nutrient dynamics and light availability in aquatic mesocosms. These 'after-life' effects of litter from different genotypes were comparable to the responses of terrestrial communities associated with the living foliage. We found that multiple litter traits corresponding with aquatic community and ecosystem responses differed in their heritability. Moreover, the underlying genetic architecture of these traits was complex, and many genes contributed only a small proportion to phenotypic variation. Our results provide further evidence that genetic variation is a key component of aquatic-terrestrial linkages, but challenge the ability to predict community or ecosystem responses based on the actions of one or a few genes.

Intraspecific variability modulates interspecific variability in animal organismal stoichiometry.

Interspecific differences in organismal stoichiometry (OS) have been documented in a wide range of animal taxa and are of significant interest for understanding evolutionary patterns in OS. In contrast, intraspecific variation in animal OS has generally been treated as analytical noise or random variation, even though available data suggest intraspecific variability in OS is widespread. Here, we assess how intraspecific variation in OS affects inferences about interspecific OS differences using two co-occurring Neotropical fishes: Poecilia reticulata and Rivulus hartii. A wide range of OS has been observed within both species and has been attributed to environmental differences among stream systems. We assess the contributions of species identity, stream system, and the interactions between stream and species to variability in N:P, C:P, and C:N. Because predation pressure can impact the foraging ecology and life-history traits of fishes, we compare predictors of OS between communities that include predators, and communities where predators are absent. We find that species identity is the strongest predictor of N:P, while stream or the interaction of stream and species contribute more to the overall variation in C:P and C:N. Interspecific differences in N:P, C:P, and C:N are therefore not consistent among streams. The relative contribution of stream or species to OS qualitatively changes between the two predation communities, but these differences do not have appreciable effects in interspecific patterns. We conclude that although species identity is a significant predictor of OS, intraspecific OS is sometimes sufficient to overwhelm or obfuscate interspecific differences in OS.

Effects of consumer interactions on benthic resources and ecosystem processes in a neotropical stream.

The effect of consumers on their resources has been demonstrated in many systems but is often confounded by trophic interactions with other consumers. Consumers may also have behavioral and life history adaptations to each other and to co-occurring predators that may additionally modulate their particular roles in ecosystems. We experimentally excluded large consumers from tile periphyton, leaves and natural benthic substrata using submerged electrified frames in three stream reaches with overlapping consumer assemblages in Trinidad, West Indies. Concurrently, we assessed visits to (non-electrified) control frames by the three most common large consumers-primarily insectivorous killifish (Rivulus hartii), omnivorous guppies (Poecilia reticulata) and omnivorous crabs (Eudaniela garmani). Consumers caused the greatest decrease in final chlorophyll a biomass and accrual rates the most in the downstream reach containing all three focal consumers in the presence of fish predators. Consumers also caused the greatest increase in leaf decay rates in the upstream reach containing only killifish and crabs. In the downstream reach where guppies co-occur with predators, we found significantly lower benthic invertebrate biomass in control relative to exclosure treatments than the midstream reach where guppies occur in the absence of predators. These data suggest that differences in guppy foraging, potentially driven by differences in their life history phenotype, may affect ecosystem structure and processes as much as their presence or absence and that interactions among consumers may further mediate their effects in these stream ecosystems.

Environmental and organismal predictors of intraspecific variation in the stoichiometry of a neotropical freshwater fish.

The elemental composition of animals, or their organismal stoichiometry, is thought to constrain their contribution to nutrient recycling, their interactions with other animals, and their demographic rates. Factors that affect organismal stoichiometry are generally poorly understood, but likely reflect elemental investments in morphological features and life history traits, acting in concert with the environmental availability of elements. We assessed the relative contribution of organismal traits and environmental variability to the stoichiometry of an insectivorous Neotropical stream fish, Rivulus hartii. We characterized the influence of body size, life history phenotype, stage of maturity, and environmental variability on organismal stoichiometry in 6 streams that differ in a broad suite of environmental variables. The elemental composition of R. hartii was variable, and overlapped with the wide range of elemental composition documented across freshwater fish taxa. Average %P composition was ∼3.2%(±0.6), average %N∼10.7%(±0.9), and average %C∼41.7%(±3.1). Streams were the strongest predictor of organismal stoichiometry, and explained up to 18% of the overall variance. This effect appeared to be largely explained by variability in quality of basal resources such as epilithon N:P and benthic organic matter C:N, along with variability in invertebrate standing stocks, an important food source for R. hartii. Organismal traits were weak predictors of organismal stoichiometry in this species, explaining when combined up to 7% of the overall variance in stoichiometry. Body size was significantly and positively correlated with %P, and negatively with N:P, and C:P, and life history phenotype was significantly correlated with %C, %P, C:P and C:N. Our study suggests that spatial variability in elemental availability is more strongly correlated with organismal stoichiometry than organismal traits, and suggests that the stoichiometry of carnivores may not be completely buffered from environmental variability. We discuss the relevance of these findings to ecological stoichiometry theory.