Price Equations for Understanding the Response of Ecosystem Function to Community Change.

AbstractThe relationship between biodiversity and ecosystem function (BEF) remains unclear in many natural ecosystems, partially for lack of theoretical and analytical tools that match common characteristics of observational community data. The ecological Price equation promises to meet this need by organizing many different species-level changes into a few ecologically meaningful categories that sum to total ecosystem function change. Current versions of the ecological Price equation focus on species richness and presence-absence. However, abundance and relative abundance are better estimated in samples and are likely showing a stronger response to global change. Here, we present a novel, abundance-based version of the ecological Price equation in both discrete and continuous forms and explain the similarities and differences between this method and a related, previously developed richness-based method. We also present new empirical techniques for applying the Price equation to ecological data. Our two demonstration analyses reveal how additive effects of increasing abundance on total function are modified by concurrent selection effects due to shifts in species' composition as well as intraspecific change in species' per capita function. The ecological Price equations derived here complement existing approaches and together offer BEF researchers analytical tools and a unifying framework for studying BEF in observational community data.

Many bee species, including rare species, are important for function of entire plant-pollinator networks.

It is important to understand how biodiversity, including that of rare species, affects ecosystem function. Here, we consider this question with regard to pollination. Studies of pollination function have typically focused on pollination of single plant species, or average pollination across plants, and typically find that pollination depends on a few common species. Here, we used data from 11 plant-bee visitation networks in New Jersey, USA, to ask whether the number of functionally important bee species changes as we consider function separately for each plant species in increasingly diverse plant communities. Using rarefaction analysis, we found the number of important bee species increased with the number of plant species. Overall, 2.5 to 7.6 times more bee species were important at the community scale, relative to the average plant species in the same community. This effect did not asymptote in any of our datasets, suggesting that even greater bee biodiversity is needed in real-world systems. Lastly, on average across plant communities, 25% of bee species that were important at the community scale were also numerically rare within their network, making this study one of the strongest empirical demonstrations to date of the functional importance of rare species.

The relative importance of pollinator abundance and species richness for the temporal variance of pollination services.

The relationship between biodiversity and the stability of ecosystem function is a fundamental question in community ecology, and hundreds of experiments have shown a positive relationship between species richness and the stability of ecosystem function. However, these experiments have rarely accounted for common ecological patterns, most notably skewed species abundance distributions and non-random extinction risks, making it difficult to know whether experimental results can be scaled up to larger, less manipulated systems. In contrast with the prolific body of experimental research, few studies have examined how species richness affects the stability of ecosystem services at more realistic, landscape scales. The paucity of these studies is due in part to a lack of analytical methods that are suitable for the correlative structure of ecological data. A recently developed method, based on the Price equation from evolutionary biology, helps resolve this knowledge gap by partitioning the effect of biodiversity into three components: richness, composition, and abundance. Here, we build on previous work and present the first derivation of the Price equation suitable for analyzing temporal variance of ecosystem services. We applied our new derivation to understand the temporal variance of crop pollination services in two study systems (watermelon and blueberry) in the mid-Atlantic United States. In both systems, but especially in the watermelon system, the stronger driver of temporal variance of ecosystem services was fluctuations in the abundance of common bee species, which were present at nearly all sites regardless of species richness. In contrast, temporal variance of ecosystem services was less affected by differences in species richness, because lost and gained species were rare. Thus, the findings from our more realistic landscapes differ qualitatively from the findings of biodiversity-stability experiments.

Measuring partner choice in plant-pollinator networks: using null models to separate rewiring and fidelity from chance.

Recent studies of mutualistic networks show that interactions between partners change across years. Both biological mechanisms and chance could drive these patterns, but the relative importance of these factors has not been separated. We established a field experiment consisting of 102 monospecific plots of 17 native plant species, from which we collected 6713 specimens of 52 bee species over four years. We used these data and a null model to determine whether bee species' foraging choices varied more or less over time beyond the variation expected by chance. Thus we provide the first quantitative definition of rewiring and fidelity as these terms are used in the literature on interaction networks. All 52 bee species varied in plant partner choice across years, but for 27 species this variation was indistinguishable from random partner choice. Another 11 species showed rewiring, varying more across years than expected by chance, while 14 species showed fidelity, indicating that they both prefer certain plant species and are consistent in those preferences across years. Our study shows that rewiring and fidelity both exist in mutualist networks, but that once sampling effects have been accounted for, they are less common than has been reported in the ecological literature.

Biodiversity-ecosystem function relationships change in sign and magnitude across the Hill diversity spectrum.

Motivated by accelerating anthropogenic extinctions, decades of biodiversity-ecosystem function (BEF) experiments show that ecosystem function declines with species loss from local communities. Yet, at the local scale, changes in species' total and relative abundances are more common than species loss. The consensus best biodiversity measures are Hill numbers, which use a scaling parameter, ℓ, to emphasize rarer versus more common species. Shifting that emphasis captures distinct, function-relevant biodiversity gradients beyond species richness. Here, we hypothesized that Hill numbers that emphasize rare species more than richness does may distinguish large, complex and presumably higher-functioning assemblages from smaller and simpler ones. In this study, we tested which values of ℓ produce the strongest BEF relationships in community datasets of ecosystem functions provided by wild, free-living organisms. We found that ℓ values that emphasized rare species more than richness does most often correlated most strongly with ecosystem functions. As emphasis shifted to more common species, BEF correlations were often weak and/or negative. We argue that unconventional Hill diversities that shift emphasis towards rarer species may be useful for describing biodiversity change, and that employing a wide spectrum of Hill numbers can clarify mechanisms underlying BEF relationships. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Life history and chemical defense interact to drive patterns of adaptation in an annual monkeyflower.

Although chemical defenses and herbivore pressure are widely established as key targets and agents of selection, their roles in local adaptation and determining potential evolutionary responses to changing climates are often neglected. Here, we explore fitness differences between 11 rangewide M. guttatus populations in a field common garden experiment and assess the agents and targets of selection driving relative fitness patterns. We use piecewise structural equation models to disentangle associations between chemical defenses, (phenylpropanoid glycosides; PPGs), and life history traits with herbivory and fitness. While the historical environment of populations is not predictive of fitness differences between populations, >90% of variation in fitness can be predicted by the flowering time and foliar PPG defense arsenal of a population. Piecewise structural equation models indicate that life history traits, particularly earlier flowering time, are strongly and directly linked to fitness. However, herbivory, particularly fruit predation, is also an important agent of selection that creates indirect links between fitness and both chemical defenses and life history traits. Our results emphasize the multivariate nature of the agents and targets of selections in producing adaptation and suggest that future responses to selection must navigate a complex fitness landscape.

Rare and declining bee species are key to consistent pollination of wildflowers and crops across large spatial scales.

Biodiversity promotes ecosystem function (EF) in experiments, but it remains uncertain how biodiversity loss affects function in larger-scale natural ecosystems. In these natural ecosystems, rare and declining species are more likely to be lost, and function needs to be maintained across space and time. Here, we explore the importance of rare and declining bee species to the pollination of three wildflowers and three crops using large-scale (72 sites across 5000 km2 ), multi-year datasets. Half of the sampled bee species (82/164) were rare or declining, but these species provided only ~15% of overall pollination. To determine the number of species important to EF, we used two methods of "scaling up," both of which have previously been used for biodiversity-function analysis. First, we summed bee species' contributions to pollination across space and time and then found the minimum set of species needed to provide a threshold level of function across all sites; according to this method, effectively no rare and declining bee species were important to pollination. Second, we account for the "insurance value" of biodiversity by finding the minimum set of bee species needed to simultaneously provide a threshold level of function at each site in each year. The second method leads to the conclusion that 25 rare and eight declining bee species (36% and 53% of all rare and declining bee species, respectively) are included in the minimum set. Our findings provide some of the strongest evidence yet that rare and declining species are key to meeting threshold levels of EF, thereby providing a more direct link between real-world biodiversity loss and EF.

Welcome to the neighbourhood: interspecific genotype by genotype interactions in Solidago influence above- and belowground biomass and associated communities.

Intra- and interspecific plant-plant interactions are fundamental to patterns of community assembly and to the mixture effects observed in biodiversity studies. Although much research has been conducted at the species level, very little is understood about how genetic variation within and among interacting species may drive these processes. Using clones of both Solidago altissima and Solidago gigantea, we found that genotypic variation in a plant's neighbours affected both above- and belowground plant traits, and that genotype by genotype interactions between neighbouring plants impacted associated pollinator communities. The traits for which focal plant genotypic variation explained the most variation varied by plant species, whereas neighbour genotypic variation explained the most variation in coarse root biomass. Our results provide new insight into genotypic and species diversity effects in plant-neighbour interactions, the extended consequences of diversity effects, and the potential for evolution in response to competitive or to facilitative plant-neighbour interactions.

Aphid and ladybird beetle abundance depend on the interaction of spatial effects and genotypic diversity.

Intraspecific variation and genotypic diversity of host-plants can affect the structure of associated arthropod communities and the dynamics of populations. Similarly, neighboring plants can also affect interactions between host-plants and their associated arthropods. However, most studies on the effects of host-plant genotypes have largely ignored the potential effects of neighboring host-plants on arthropod communities. In this study, we used a common garden experiment to ask how spatial effects of neighboring patches, along with genotype identity and genotypic diversity in tall goldenrod (Solidago altissima), affect the abundances of a common goldenrod herbivore (Uroleucon nigrotuberculatum) and their dominant predator (Harmonia axyridis, a ladybird beetle). Aphid abundance varied 80-fold among genotypes, while ladybird beetle abundance was not affected by genotype identity. Additionally, there were strong effects of neighboring plots: aphid abundance in a focal plot was positively correlated to aphid abundance in nearby plots, suggesting strong spatial patterning in the abundance of aphids. Neither aphid nor ladybird beetle abundance was affected by genotypic diversity. However, focal plot genotypic diversity mediated the strength of the neighborhood effect (i.e., strong effects for genotype polyculture focal plots and weak effects for genotype monoculture focal plots). Our results show that aphids were directly influenced by host-plant genotype identity while ladybird beetles responded mainly to prey abundance, and suggest that genotypic diversity can influence the effects of spatial processes on the plant-herbivore interactions.

Evolutionary history determines how plant productivity responds to phylogenetic diversity and species richness.

The relationship between biodiversity and ecosystem function has received a great deal of attention in ecological research and recent results, from re-analyses, suggest that ecosystem function improves with increases in phylogenetic diversity. However, many of these results have been generalized across a range of different species and clades, and plants with different evolutionary histories could display different relationships between biodiversity and ecosystem function. To experimentally test this hypothesis, we manipulated species richness and phylogenetic diversity using 26 species from two subgenera of the genus Eucalyptus (subgenus Eucalyptus and subgenus Symphyomyrtus). We found that plant biomass (a measurement of ecosystem function) sometimes, but not always, responded to increases in species richness and phylogenetic diversity. Specifically, Symphyomyrtus plants showed a positive response while no comparable effect was observed for Eucalyptus plants, showing that responses to biodiversity can vary across different phylogenetic groups. Our results show that the impacts of evolutionary history may complicate the relationship between the diversity of plant communities and plant biomass.

Belowground interactions shift the relative importance of direct and indirect genetic effects.

Intraspecific genetic variation can affect decomposition, nutrient cycling, and interactions between plants and their associated belowground communities. However, the effects of genetic variation on ecosystems can also be indirect, meaning that genes in a focal plant may affect ecosystems by altering the phenotype of interacting (i.e., neighboring) individuals. We manipulated genotype identity, species identity, and the possibility of belowground interactions between neighboring Solidago plants. We hypothesized that, because our plants were nitrogen (N) limited, the most important interactions between focal and neighbor plants would occur belowground. More specifically, we hypothesized that the genotypic identity of a plant's neighbor would have a larger effect on belowground biomass than on aboveground biomass, but only when neighboring plants were allowed to interact belowground. We detected species- and genotype-level variation for aboveground biomass and ramet production. We also found that belowground biomass and ramet production depended on the interaction of neighbor genotype identity and the presence or absence of belowground interactions. Additionally, we found that interspecific indirect genetic effects (IIGEs; changes in focal plant traits due to the genotype identity of a heterospecific neighbor) had a greater effect size on belowground biomass than did focal genotype; however, this effect only held in pots that allowed belowground interactions. These results expand the types of natural processes that can be attributed to genotypes by showing that, under certain conditions, a plant's phenotype can be strongly determined by the expression of genes in its neighbor. By showing that IIGEs are dependent upon plants being able to interact belowground, our results also provide a first step for thinking about how genotype-based, belowground interactions influence the evolutionary outcomes of plant-neighbor interactions.

The afterlife of interspecific indirect genetic effects: genotype interactions alter litter quality with consequences for decomposition and nutrient dynamics.

Aboveground-belowground linkages are recognized as divers of community dynamics and ecosystem processes, but the impacts of plant-neighbor interactions on these linkages are virtually unknown. Plant-neighbor interactions are a type of interspecific indirect genetic effect (IIGE) if the focal plant's phenotype is altered by the expression of genes in a neighboring heterospecific plant, and IIGEs could persist after plant senescence to affect ecosystem processes. This perspective can provide insight into how plant-neighbor interactions affect evolution, as IIGEs are capable of altering species interactions and community composition over time. Utilizing genotypes of Solidago altissima and Solidago gigantea, we experimentally tested whether IIGEs that had affected living focal plants would affect litter decomposition rate, as well as nitrogen (N) and phosphorous (P) dynamics after the focal plant senesced. We found that species interactions affected N release and genotype interactions affected P immobilization. From a previous study we knew that neighbor genotype influenced patterns of biomass allocation for focal plants. Here we extend those previous results to show that these changes in biomass allocation altered litter quality, that then altered rates of decomposition and nutrient cycling. Our results provide insights into above- and belowground linkages by showing that, through their effects on plant litter quality (e.g., litter lignin:N), IIGEs can have afterlife effects, tying plant-neighbor interactions to ecosystem processes. This holistic approach advances our understanding of decomposition and nutrient cycling by showing that evolutionary processes (i.e., IIGEs) can influence ecosystem functioning after plant senescence. Because plant traits are determined by the combined effects of genetic and environmental influences, and because these traits are known to affect decomposition and nutrient cycling, we suggest that ecosystem processes can be described as gene-less products of genetic interactions among the species comprising ecological communities.

Non-additive effects of genotypic diversity increase floral abundance and abundance of floral visitors.

BACKGROUND: In the emerging field of community and ecosystem genetics, genetic variation and diversity in dominant plant species have been shown to play fundamental roles in maintaining biodiversity and ecosystem function. However, the importance of intraspecific genetic variation and diversity to floral abundance and pollinator visitation has received little attention. METHODOLOGY/PRINCIPAL FINDINGS: Using an experimental common garden that manipulated genotypic diversity (the number of distinct genotypes per plot) of Solidago altissima, we document that genotypic diversity of a dominant plant can indirectly influence flower visitor abundance. Across two years, we found that 1) plant genotype explained 45% and 92% of the variation in flower visitor abundance in 2007 and 2008, respectively; and 2) plant genotypic diversity had a positive and non-additive effect on floral abundance and the abundance of flower visitors, as plots established with multiple genotypes produced 25% more flowers and received 45% more flower visits than would be expected under an additive model. CONCLUSIONS/SIGNIFICANCE: These results provide evidence that declines in genotypic diversity may be an important but little considered factor for understanding plant-pollinator dynamics, with implications for the global decline in pollinators due to reduced plant diversity in both agricultural and natural ecosystems.