Widespread overharvesting of top consumers of the world’s ecosystems has “skewed” food webs, in terms of biomass and species richness, towards a generally greater domination at lower trophic levels. This skewing is exacerbated in locations where exotic species are predominantly low-trophic level consumers such as benthic macrophytes, detritivores, and filter feeders. However, in some systems where numerous exotic predators have been added, sometimes purposefully as in many freshwater systems, food webs are skewed in the opposite direction toward consumer dominance. Little is known about how such modifications to food web topology, e.g., changes in the ratio of predator to prey species richness, affect ecosystem functioning. We experimentally measured the effects of trophic skew on production in an estuarine food web by manipulating ratios of species richness across three trophic levels in experimental mesocosms. After 24 days, increasing macroalgal richness promoted both plant biomass and grazer abundance, although the positive effect on plant biomass disappeared in the presence of grazers. The strongest trophic cascade on the experimentally stocked macroalgae emerged in communities with a greater ratio of prey to predator richness (bottom-rich food webs), while stronger cascades on the accumulation of naturally colonizing algae (primarily microalgae with some early successional macroalgae that recruited and grew in the mesocosms) generally emerged in communities with greater predator to prey richness (the more top-rich food webs). These results suggest that trophic skewing of species richness and overall changes in food web topology can influence marine community structure and food web dynamics in complex ways, emphasizing the need for multitrophic approaches to understand the consequences of marine extinctions and invasions.
Citation: Reynolds PL, Bruno JF (2012) Effects of Trophic Skewing of Species Richness on Ecosystem Functioning in a Diverse Marine Community. PLoS ONE 7(5): e36196. doi:10.1371/journal.pone.0036196
Editor: Jon Moen, Umea University, Sweden
Received: January 23, 2012; Accepted: April 3, 2012; Published: May 31, 2012
Copyright: © 2012 Reynolds, Bruno. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: Funding support came from University of North Carolina’s Institute of Marine Sciences, National Science Foundation grant OCE-0623736 to Bruno, and a University of North Carolina Wilson Award to Reynolds. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Species losses from habitat destruction and overharvesting, and species gains from accidental and intentional introductions, are changing the topology of food webs with consequences for ecosystem functioning , . Despite ongoing losses of species richness at the global scale , the rate of species gain from introductions at local scales can outpace those lost to extinctions with potentially little effect on overall community diversity , . However, inherent species biases in extinction and invasion processes are altering the distribution of diversity across trophic levels ,  with potential effects on ecosystems . While natural food webs are thought to be slightly weighted toward greater species richness at lower trophic levels , biases in which species are more likely to be lost and gained can result in food webs skewed toward ratios of greater or lower predator to prey richness.
Local species extirpations are often biased toward species at higher trophic levels. Large consumers such as top predators are generally more likely to go extinct due to their characteristic small population sizes, small geographic ranges, slow population growth, and high susceptibility to over-harvesting and habitat loss , , , , . However, recent analyses on global fisheries stock assessments and landings indicate that small-bodied consumer populations are also highly susceptible to collapse, with substantial impacts on adjacent trophic levels in oceanic food webs .
Despite documented global declines in consumer guilds , there are many examples of local increases in predator richness due to species introductions and range shifts, such as those observed in fishes, snakes and birds in the Everglades , , lionfish in the Caribbean , , and seagrass-associated fishes in the Gulf of Mexico . Intentional or accidental releases of predatory game fish in streams and lakes modified for human recreation have also elevated local predator richness in North American and European waterways , . From brown tree snakes in Guam, Burmese pythons in the Everglades, Asian carp in the Great Lakes, cactus moths in Mexico, cane toads in Australia, and black rats and feral cats around the world, many striking examples of the negative impacts of invasive species come from the establishment of exotic consumers , . However, in general localized species gains may be inherently biased toward species at lower trophic levels due to mechanisms of transport and establishment. For example, in estuarine systems species at lower trophic levels, such as macroplanktivores and plants (particularly species found in ballast water), are more likely to be transported and become established . Highly evolved dispersal and hitchhiking abilities, as well as traits such as fast growth and reproduction of plants, fungi, and invertebrates may promote their transport as well as spread and establishment in new regions .
Regardless of the direction, trophic skew (a re-organization of trophic structure due to a change in the ratio of predator to prey richness , ) is changing the structural biodiversity of natural food webs with unknown consequences for ecosystem processes , , . However, understanding the consequences of shifts in species richness across multiple trophic levels is difficult because diversity effects are often contingent upon the presence and diversity of adjacent trophic levels , . Capture and consumption of prey from within a given trophic level can be influenced by richness at lower and higher trophic levels , and diversity can affect production through fundamentally different mechanisms across trophic levels . Thus explicit investigation of how species gains and losses, both consequences and drivers of global change, simultaneously affect ecosystem functioning is key to understanding current and future responses of natural ecosystems to trophic skew .
While a wealth of research over the past decade indicates that changes in biodiversity alter ecosystem functioning and services , , , most studies have focused on manipulating richness at only one or, more recently and rarely, two trophic levels . Manipulations of species richness across multiple trophic levels typically employ designs in which overall community richness varies across treatments , and/or feature limited species pools where results are largely ascribed to identity effects and changes in community composition . Because elevated system diversity is known to affect production and other ecosystem properties , it is often difficult to partition the effects of changes in predator to prey richness from overall changes in diversity at the community level. Expect perhaps work by Srivastava et al. , to date, no study has explicitly tested the effects of trophic skew on specific ecosystem properties including the strength of a trophic cascade.
We empirically tested how trophic skew can influence ecosystem functioning in a diverse tri-trophic estuarine food web in outdoor mesocosms by simultaneously manipulating plant and predator species richness while holding overall community richness constant. Using a diverse species pool, we created four food web structures (Fig. 1A) that reflected real or predicted degrees of trophic skewing of natural systems: 1) top-rich food webs (inverted triangle shape) with greater ratios of predator to prey richness, consistent with predictions based on accidental and intentional predator additions , 2) neutral (rectangular shape) with constant predator to prey richness ratios, and 3–4) two degrees of bottom-rich food webs (triangular shapes) with greater ratios of prey to predator richness, as may be typical for impacted estuarine food webs . We tested the direct and indirect effects of predator (top-rich) and plant (bottom-rich) diversity , ,  by comparing production and trophic cascade strength across these trophic skewing scenarios. With concurrent and opposite changes in prey and predator richness, it is possible that 1) the species richness at one trophic level may dominate and dictate final primary producer biomass, 2) effects of concurrent changes in richness may cancel out, resulting in constant primary production across different food web structures, or 3) the concurrent changes in richness may interact additively or nonadditively , .
A) Structure of experimental food webs with varying species richness per trophic level, i.e., trophic skew, and B) experimental species pool (predators: Hypleurochilus geminatus, Monacanthus hispidus, Fundulus heteroclitus, swimming crabs (Callinectes sapidus or C. similis), Lagodon rhomboids, Penaeus aztecus, Palaemonetes vulgaris, mud crabs (Panopeus herbstii, Eurypanopeus depressus or Neopanope sayi); grazers: Gammarus mucronatus, Elasmopus levis, Dulichiella appendiculata, Paracerceis caudata; macroalgae: Dictyota menstrualis, Codium fragile, Padina gymnospora, Sargassum filipendula, Ceramium sp., Gracilaria tikvahiae, G. verrucosa, Hypnea musciformis, Ulva lactuca. Species are not drawn to scale.
Based on known dynamics in this and similar experimental systems where predator richness promoted predation and increased top-down control of grazer populations , , we expected that top-rich food webs (greater predator to prey richness) would exhibit lower grazer abundances and higher algal biomass compared to bottom-rich food webs (lower predator to prey richness). However, algal richness is known to promote primary as well as secondary production , , and thus bottom-up food webs could have either higher or lower algal biomass depending on the intensity of grazing. If predator (top-down) richness effects drive this system, we expected top-rich food webs to have a stronger positive trophic cascade on algal biomass compared to bottom-rich food webs, and the reverse if algal (bottom-up) richness is more important.
The presence of consumers and the distribution of species richness across trophic levels (trophic skew) influenced final algal biomass and grazer abundance (Fig. 2, Table 1). Generally, grazers reduced the wet mass of experimentally stocked macroalgae (hereafter referred to as macroalgae) by 33% compared to grazer-free controls (LSM contrast, F1,83 = 42.64, P<0.001), but this effect disappeared when their predators were present, indicating a trophic cascade (Fig. 2A; Table 1; see Fig. S2 for images of final macroalgal biomass). Predators generally increased macroalgal biomass (LSM contrast, F1,83 = 10.44, P<0.002). The strongest positive trophic cascade on macroalgae appeared in the bottom-rich food webs (e.g., triangular shaped, Fig. 2D).
Final A) biomass of experimentally stocked macroalgae, B) colonizing algae (e.g., microalgae) chlorophyll a accumulation C) mesograzer abundance, and the trophic cascade strength on macroalgae (D) and microalgae (E) in experimental mesocosms across trophic skewing treatments after 24 days. Higher values in D and E indicate a stronger positive trophic cascade on the algae. Values are means±1SE.
Microalgae (primarily the chain-forming diatom Tabellaria sp.) and some early-successional Cladophora sp. and Ulva linza colonized and grew in all mesocosms. Chlorophyll a concentration, a proxy for growth of these naturally recruiting algae (hereafter referred to as “microalgae” due to their small size and dominance by diatoms), was affected by the presence of upper trophic levels (food chain length treatment) but not by trophic skew (Fig. 2B, Table 1). On average, grazers reduced chlorophyll a concentration by 83% in the absence of predators, and 66% in their presence (Fig. 2B). Predators generally promoted microalgal accumulation (LSM contrast, F1,83 = 4.19, P = 0.044). In contrast with the macroalgal results, the strongest trophic cascade on microalgal accumulation emerged in the top-rich food webs (inverted triangle) with a greater ratio of predator to prey species richness, although this response was highly variable (Fig. 2E).
Incidental grazer immigration was minimal across all food web structures (∼39 individuals per algae-only control mesocosm, Fig. 2C). Both the presence of predators and changes in predator to prey richness influenced final grazer abundance (Table 1) and community composition (Fig. 3). The inclusion of predators reduced grazer abundance by 91%. Grazer abundance was greatest in the two bottom-rich food web structures regardless of the presence of predators (LSM contrast F1,56 = 27.53, P<0.0001, Fig. 2C) suggesting that macroalgal richness promoted grazer population growth. Grazer communities in treatments in which grazers were initially added (e.g., excluding algae-only treatments) were dominated by the amphipod Elasmopus levis (Fig. 3). Experimental grazer communities were significantly affected by predator presence (PERMANOVA, F1,56 = 97.032, P<0.001) and by the trophic skew treatment (F3,56 = 4.621, P = 0.006), as well as by the interaction between the two (F3,56 = 3.656, P = 0.009; Fig. 4). However, data failed the test of multivariate homogeneity (e.g., permutation analysis of betadisper results for homogeneity across trophic skew treatments, F3,60 = 2.804, P = 0.031) and treatment effects on grazer community composition should be interpreted cautiously.
Final grazer abundance per species across trophic skew treatments in the presence and absence of predators. Values are means±1SE.
Plot of the nonmetric multidimensional scaling (NMDS) scores for the two most important axes for all replicates using final grazer species composition data. Circles and stars represent the absence and presence of predators, respectively. Trophic skew treatments are represented from light to dark as follows: top-rich, neutral, slightly bottom-rich, and bottom-rich.
Changing the structure of food webs, via reductions in length and diversity, can alter ecosystem functioning , . The results of our experiment suggest that trophic skewing of species richness (changes in the ratio of predator to prey richness) can also affect ecosystem function in a diverse estuarine food web. Specifically, this trophic skewing of species richness can affect primary and secondary production as well as the strength of a trophic cascade.
Similar to results from other aquatic and terrestrial studies, we found that increasing primary producer richness increased production (primary and secondary biomass) after 24 days in our experimental communities , , , , . The observed positive effects of increasing plant richness on plant biomass are thought to be largely due to complementarity and sampling effects, whereby the likelihood of incorporating facilitators and highly productive species, as well as heightening resource partitioning, increases with elevated richness .
Although increasing plant richness may increase plant biomass, these effects can be weak or undetectable in the presence of consumers due to grazing compensation , , , . Here we found that increasing macroalgal richness from two to six species elevated macroalgal biomass (Fig. 2A), but this effect was not evident in the presence of grazers. This adds to a growing list of experimental manipulations demonstrating reduced diversity effects on production with the inclusion of a higher trophic level , i.e., under ecologically realistic conditions. However, it is unclear whether the added production is transferred up food chains resulting in greater production at top trophic levels .
Although a recent meta-analysis supports the concept that prey richness is a strong predictor of consumer effects and that more diverse prey communities can maintain higher abundances by containing species that are highly tolerant or resistant to consumers , we did not see a suppression of grazer effects on macroalgal biomass with increasing macroalgal diversity (Fig. 2A). Increasing macroalgal richness did not decrease the strength of top-down control by consumers , , . It is possible, as with many short-term biodiversity studies in which plant biomass per species is initially low, that this experiment was not of sufficient duration to allow unpalatable algal species to replace the biomass of the more palatable species lost to herbivory . Given that macroalgal diversity effects in this system are known to vary in space and time , the lack of environmental variation within the mesocosms may also have reduced our ability to detect macroalgal richness effects on primary production in the presence of grazers , . Alternatively, as several of the macroalgal species in our species pool are known to be highly productive and readily consumed by these and other small crustacean grazers , , , , it is possible that increasing macroalgal richness increased overall macroalgal community palatability or the likelihood of incorporating a preferred macroalgal species (Table S1A), thereby promoting grazer populations which limited primary producer biomass. Elevating plant richness and the number of plant functional groups can promote herbivore diversity and overall herbivory , , potentially leading to the observed stronger top-down control in the most diverse macroalgal communities (Fig. 2A). This decoupling of plant richness and productivity in the presence of consumers may be a general trend in this and other systems with strong top down control , .
Generally, all grazer populations grew in response to the different experimental assemblages (Figs. 2C, 3). However, grazer abundance was greater in the more diverse macroalgal communities (Fig. 2C). Increased plant diversity can promote herbivore production directly through increased diet breadth (balanced diet hypothesis) and overall resource availability , , and indirectly by influencing herbivore-predator interactions through increased prey refuge availability and/or quality , . Although studies in similar and other systems found that predators, while capable of destabilizing herbivore populations, had little effect on plants , , , here we found that predator presence can affect primary producers (Fig. 2). Elevating predator richness can promote primary producers by increasing prey capture – due to diet complementarity , , –and prey antipredator behaviors such as reduced grazing, thereby releasing plants from herbivory and strengthening a trophic cascade , .
The effects of trophic skew on the strength of a trophic cascade in this model estuarine system, however, varied between the two primary producer groups (experimentally stocked macroalgae and naturally colonizing microalgae, Fig. 2D,E). Increasing prey richness relative to predator richness (e.g., from top-rich to bottom-rich food webs) strengthened the trophic cascade on macroalgae and generally promoted standing macroalgal biomass, but potentially weakened the trophic cascade on microalgae. A recent meta-analysis of work in detrital food webs, however, found that top-down effects of changing consumer (detritivore) richness had stronger effects on functioning (e.g., decomposition) than the bottom-up effects of changing resource (detritus) diversity . Contrasts between this work in a ‘brown’ versus our study in a marine ‘green’ web may reflect inherent differences in processes associated with the two systems including the role of dynamic responses of ‘live’ resources to their consumers .
Increasing predator richness in our experimental communities was likely correlated with increased likelihood and potential promotion of omnivory (both a sampling and nonadditive richness effect), thereby promoting overall consumption of macroalgae and appearing to weaken the trophic cascade on these producers in top-rich food webs , . Additionally, it is possible that the variable effects on producer functional type (experimentally stocked macroalgae versus the naturally colonizing microalgae) were due to changes in grazer foraging behavior. Grazer abundance increased in response to elevated plant richness (e.g., bottom-rich food webs), potentially increasing inter-and intraspecific competition and promoting grazer consumption of microalgae. Elevated predator richness in the top-rich food webs may have reduced grazer activity and increased use of macroalgal refuges such as Dictyota menstrualis and Ulva lactuca , , , reducing access to and overall consumption of microalgae in these treatments. Thus, trophic skew may affect different types and dynamics of primary producers, making it difficult to predict the overall effects of concurrent species gains and losses on primary production.
Community composition and identity effects may drive observed differences among experimentally skewed food webs. Elevated macroalgal richness, coupled with decreased predator richness, may have promoted grazer survivorship and population growth through increased refuge and food quality, and/or decreased predator efficiency of prey capture , . An increased likelihood of incorporating unpalatable algae due to higher algal richness and a decreased likelihood of omnivory due to lower predator richness (Table S1A, B) could further promote algae in bottom-up skewed communities. Reduced predator richness in these food webs could also reduce predator efficiency if intraspecific competition among predators is stronger than interspecific interactions. As food webs become skewed, the identity of the species being gained or lost (e.g., whether they are an omnivorous predator or a palatable plant) as well as the trophic level in which they reside may become increasingly important.
Biodiversity can significantly affect primary production, nutrient cycling and community composition. Control of algal blooms, the yield of important commercial and recreational fisheries, and other ecosystem services may depend not only on the maintenance of biodiversity, but also on its distribution throughout a given food web. As substantial losses of species at local and global scales as well as local gains of non-native species are expected to continue , , , , understanding the effects of trophic skew on ecosystem functioning may be important for predicting the potentially synergistic effects of species extinctions and invasions on ecosystem functioning and will be an important challenge for empirical and applied endeavors across systems.
Materials and Methods
Mesocosms and Experimental Design
The experiment was performed in outdoor, flow through mesocosms at the University of North Carolina at Chapel Hill’s Institute of Marine Sciences (IMS) in Morehead City, NC in July 2007. We independently manipulated secondary consumer (referred to as predator throughout) and macroalgal richness to create four different food web structures (trophic skew treatment) mimicking different trophic skewing scenarios, with constant total community richness (Fig. 1A). Experimental food webs were skewed to be top-rich (2 macroalgal: 4 grazer: 6 predator species), neutral (4∶4∶4), slightly bottom-rich (5∶4∶3), or bottom-rich (6∶4∶2). Trophic skew treatments were crossed by three manipulations of food chain length (algae only, algae+grazers, algae+grazers+predators) to compare changes in trophic cascade strength across experimental communities (n = 9; 108 mesocosms total).
We used a substitutive design, manipulating initial algal and predator richness and identity while holding biomass and density constant per mesocosm (35 g algae and 6 predator individuals, c. 18 g) at densities comparable with natural levels in North Carolina subtidal communities , , , . Algal and predator species composition per replicate were chosen randomly from a larger pool of nine macroalgae and eight predator species (Fig. 1B). Selection of species from these larger species pools allowed for a conservative test of richness effects in our system by varying community composition, but not richness, for replicates within a given treatment. This controlled for species identity and composition effects . Algal composition varied among replicates and levels of trophic skew, but was replicated across the different food chain length treatments to account for the effects of initial algal composition on the final response variables. Initial grazer richness, composition and abundance were replicated such that each treatment received a mixture of four grazer species. We chose a substitutive design because it does not confound diversity with density, as do additive designs , . Although replacement designs can potentially diminish intraspecific interactions by decreasing species-specific densities with increasing species richness , , they are useful in detecting complementarity effects .
Replicates were maintained in 30 L clear plastic aquaria provided with gravel-filtered seawater from Bogue Sound (see Fig. S1 for images of experimental mesocosms). Seawater flowed through 100 µm nylon mesh filter bags and was delivered through a dump bucket system that maintained aeration and approximated natural subtidal turbulence , , . This system limited outside grazer recruitment, but not the passage of microscopic algal propagules that colonized the mesocosms (referred to throughout as “microalgae” as they were primarily composed of diatoms). Mesocosms were covered with 5 cm opening Vexar mesh lids and were placed in water tables to maintain constant temperature. Light, temperature and salinity within the mesocosms closely approximated field conditions in the nearby Bogue Sound during the course of the experiments , . Mesocosms were rotated every 5 days to reduce positioning artifacts.
Study System and Organisms
The South Atlantic Bight hard-substratum communities are highly diverse, composed of tropical and temperate species of algae, invertebrates and fishes . Macro-and epiphytic algae, the main primary producers in this system, are intensely grazed by a diverse macroinvertebrate community composed largely of amphipods and isopods , , which are in turn consumed by an array of invertebrate and vertebrate predators including shrimp, crabs and fishes , . Experimental communities featured local algal, grazer, and predator species that commonly co-occur and typically dominate hard-substratum sites of North Carolina estuaries. Organisms were collected and cultured or maintained in outdoor water tables at IMS prior to experimentation.
Chosen macroalgal species (Fig. 1B) are common in NC estuaries, although their abundances fluctuate seasonally , , . We attached seven algal thalli haphazardly to 25×25 cm Vexar mesh screens (with 5 mm openings), which were secured to the bottom of each 30 L polypropylene mesocosm such that algae floated upright in natural orientation. Initial total algal biomass per mesocosm was held constant at 35 g with approximately 5 g per individual thalli attachment. Initial algal biomass was purposefully lower than field densities  in order to allow room for growth. Algal biomass was determined after first immersing the algae in seawater for at least 15 minutes and then spinning it 15 revolutions in a salad spinner to remove excess water. We dipped algae in a diluted pesticide, Sevin, and rinsed it with fresh seawater to remove existing mobile invertebrates and trace pesticides before placement in mesocosms , , . Any fouling organisms (e.g., tunicates) were removed by hand prior to weighing; algal pieces with encrusting invertebrates (e.g., bryozoans) were replaced.
Mesocosms received an initial equal volume of grazers from a mixture of three amphipods (Dulichiella appendiculata, Gammarus mucronatus, and Elasmopus levis) and one isopod (Paracerceis caudata) prior to predator additions. These mesograzers are common in NC estuaries, achieving densities of 10–145 individuals g−1 wet mass of the alga Sargassum filipendula . They also have short generation times, respond quickly to changes in habitat and predation, and consume various types of macro-and microalgae , , . Each subsequent week an additional equal volume of grazers was added to each mesocosm to mimic natural recruitment  and to remove the possibility of predator overexploitation (for a total of c. 120 herbivores per mesocosm overall). Volume additions were subsampled (n = 20) and composed mostly of E. levis for the initial additions, and D. appendiculata and P. caudata for the recruitment additions. The initial volume addition was supplemented with five individuals of each grazer species to ensure that all replicates received the same grazer richness. Grazers were stocked within the lower range of natural field densities to allow for natural reproduction and population growth throughout the experiment , , .
Predator assignments were chosen randomly from a pool of functionally distinct invertebrates and fishes including omnivorous and strictly carnivorous species (Table 1). Due to low field abundances it was impossible to collect enough of any one of the mud and swimming crab species. To resolve this issue without risking elevating richness, each replicate assigned to either of these crab groups received individuals of only one species for that group throughout the duration of the experiment. Each mesocosm received six individual predators, which is within natural field densities . Predators were collected within their respective average juvenile size classes. This was most important for L. rhomboides, which ontogenetically switches from a strict predator to an omnivore at 3.5 cm , , or around 2.5 g (feeding trial pilot study, n = 8). Total predator mass per mesocosm varied (0.29–3.71 g), but did not differ greatly among treatments. Mesocosms were checked daily and dead or stressed predators were replaced throughout the experiment. This predator press design maintained the potential for species interactions (e.g., intraguild predation), although it precluded direct, long-term effects of such encounters on lower trophic levels.
Responses and Analyses
After 24 days we quantified treatment effects on grazer abundance, the biomass of all experimentally stocked macroalgae, and the accumulation of naturally colonizing algae (“microalgae”). This endpoint was based on observable changes in algal growth among treatments and represented approximately two overlapping grazer generations . Grazers were preserved in 70% EtOH and later identified and counted. To assess microalgal production, we measured the chlorophyll a concentration from standardized 2 cm2 samples scraped from the side of each mesocosm. We extracted and quantified chlorophyll a concentration as in Bruno and O’Connor  to quantify microalgal accumulation. Trophic cascade strength was assessed as (Rp–Rg)/(Ra–Rg) where R is the macro-or microalgal response when in the presence of grazers (g) and predators (p), or the absence of both consumer groups(a).
We used a linear model with the fixed effects of trophic skew (4 levels) and food chain length (3 levels) and their interaction, with initial algal composition type as a blocking factor using the procedure PROC MIXED in SAS (v. 2.11, Cary, NC) to analyze treatment effects on all response variables. Algal composition type had no effect and was subsequently removed from the model and the analyses re-run without this blocking factor . Replicates were excluded from the analysis if the mesocosms cracked or received inadequate flow, or the predators exhibited chronic stress (13 mesocosms total). Algae-only replicates were excluded in the analysis of treatment effects on final grazer abundance as these mesocosms did not initially receive any grazers and total immigration was minimal (Fig. 3). Results for all analyses were transformed as necessary to meet the assumptions of normality and heteroscedacity .
Non-metric multidimensional scaling (NMDS)  conducted with the metaMDS function in the vegan package in R v. 2.14.0  was used to graphically examine changes in grazer community composition across trophic skew treatments in the absence and presence of their predators. Bray-Curtis distance was used as it performs well in simulations for ecological data . The results were plotted in 2-dimensions, and the envfit procedure in vegan  used to overlay species vectors. Replicate points which occur close together on the NMDS axes are similar in species composition. A permutational analysis of variance (PERMANOVA) using the adonis function in the vegan package was then used to test treatment effects on grazer community composition. Multivariate homogeneity of the treatments was assessed using the betadisper function in vegan, which is analogous to Levene’s test  for equality of variances.
Experimental mesocosm setup. Clockwise: mesocosms receiving filtered seawater, macroalgal community, and (courtesy of M. O’Connor) mesocosm side view.
Images of final experimental algal communities. Example algal communities after 24 days in experimental mesocosms exposed to different types and degrees of trophic skew. Images within the same column featured the same initial macroalgal community.
Experimental design. Species composition of A) macroalgal and B) predator community in each mesocosm (experimental unit). Initial macroalgal biomass and predator abundance per species changed with species richness in a substitutive design. Abundance of predators by species is denoted in parentheses when multiple individuals were present. Macroalgae species in shaded cells are known to be chemically defended and less preferred by most of our experimental mesograzers, while species in lighter cells are increasingly palatable. Predator species in shaded cells are omnivorous, while species in light cells are strict carnivores.
We thank J.E. Duffy and M.I. O’Connor for comments on this manuscript, J.E. Byrnes and J.J. Stachowicz for intellectual support, E. Moran and C. Nicolini for field assistance, A. Ciesielski and D. Crowley for laboratory assistance, and M. Whalen for images of Gammarus and Elasmopus. Research was conducted under IACUC protocol ID 07-142.0 at UNC-CH. The University of North Carolina at Chapel Hill’s Libraries provided support for open access publication.
Conceived and designed the experiments: PLR JFB. Performed the experiments: PLR. Analyzed the data: PLR JFB. Contributed reagents/materials/analysis tools: PLR JFB. Wrote the paper: PLR JFB.
- 1. Ray J, Redford KH, Steneck R, Berger J, editors. (2005) Large carnivores and the conservation of biodiversity. Washington, DC: Island Press.
- 2. Duffy JE (2003) Biodiversity loss, trophic skew and ecosystem functioning. Ecology Letters 6: 680–687.
- 3. (2005) Ecosystems and human well-being: synthesis. Washington, DC: Island Press.
- 4. Sax DF, Gaines SD (2003) Species diversity: From global decreases to local increases. Trends in Ecology & Evolution 18: 561–566.
- 5. Bardgett RD, Wardle DA (2010) Aboveground-belowground linkages: Biotic interactions, ecosystem processes, and global change. Oxford, UK: Oxford University Press.
- 6. Byrnes JE, Reynolds PL, Stachowicz JJ (2007) Invasions and extinctions reshape coastal marine food webs. PLoS ONE 2: e295.
- 7. Wardle DA, Bardgett RD, Callaway RM, Van der Putten WH (2011) Terrestrial ecosystem responses to species gains and losses. Science 332: 1273–1277.
- 8. Petchey OL, Downing AL, Mittelbach GG, Persson L, Steiner CF, et al. (2004) Species loss and the structure and functioning of multitrophic aquatic systems. Oikos 104: 467–478.
- 9. Tracy CR, George TL (1992) On the determinants of extinction. The American Naturalist 139: 102–122.
- 10. McKinney ML (1997) Extinction vulnerability and selectivity: Combining ecological and paleontological views. Annual Review of Ecology and Systematics 28: 495–516.
- 11. Terborgh J, Lopez L, V PNe, Rao M, Shahabuddin G, et al. (2001) Ecological meltdown in predator-free forest fragments. Science 294: 1923–1926.
- 12. Collette BB, Carpenter KE, Polidoro BA, Juan-Jordá MJ, Boustany A, et al. (2011) High value and long life–double jeopardy for tunas and billfishes. Science 333: 291–292.
- 13. Cardillo M, Mace GM, Jones KE, Bielby J, Bininda-Emonds ORP, et al. (2005) Multiple causes of high extinction risk in large mammal species. Science 309: 1239–1241.
- 14. Pinsky ML, Jensen OP, Ricard D, Palumbi SR (2011) Unexpected patterns of fisheries collapse in the world’s oceans. Proceedings of the National Academy of Sciences 108: 8317–8322.
- 15. Estes JA, Terborgh J, Brashares JS, Power ME, Berger J, et al. (2011) Trophic downgrading of planet earth. Science 333: 301–306.
- 16. Ferriter A, Serbesoff-King K, Bodle M, Goodyear C, Doren B, et al. (2004) Exotic species in the Everglades protection area. 2004 Everglades consolidated report. West Palm Beach, FL: South Florida Water Management District. pp. 11–15.
- 17. Krysko KL, Burgess JP, Rochford MR, Gillette CR, Cueva D, et al. (2011) Verified non-indigenous amphibians and reptiles in Florida from 1863 through 2010: Outlining the invasion process and identifying invasion pathways and stages. Zootaxa 3028: 1–64.
- 18. Albins MA, Hixon MA (2008) Invasive Indo-Pacific lionfish Pterois volitans reduce recruitment of Atlantic coral-reef fishes. Marine Ecology Progress Series 367: 233–238.
- 19. Lesser MP, Slattery M (2011) Phase shift to algal dominated communities at mesophotic depths associated with lionfish (Pterois volitans) invasion on a Bahamian coral reef. Biol Invasions 13: 1855–1868.
- 20. Fodrie FJ, Heck KL, Powers SP, Graham WM, Robinson KL (2010) Climate-related, decadal-scale assemblage changes of seagrass-associated fishes in the northern Gulf of Mexico. Global Change Biology 16: 48–59.
- 21. Eby LA, Roach WJR, Crowder LBC, Stanford JA (2006) Effects of stocking-up freshwater food webs. Trends in Ecology & Evolution 21: 576–584.
- 22. Gido KB, Brown JH (1999) Invasion of North American drainages by alien fish species. Freshwater Biology 42: 387–399.
- 23. Pitt WC, Witmer GW (2007) Invasive Predators: a synthesis of the past, present, and future. In: Elewa AMT, editor. Predation in organisms: A distinct phenomenon. Berlin: Springer. pp. 265–293.
- 24. Simberloff D, Rejmánek M, editors. (2011) Encyclopedia of biological invasions. Berkeley: University of California Press.
- 25. Sax DF, Stachowicz JJ, Gaines SD (2005) Species invasions: Insights into ecology, evolution, and biogeography. Sunderland, MA: Sinauer Associates Inc.
- 26. Dobson A, Lodge D, Alder J, Cumming GS, Keymer J, et al. (2006) Habitat loss, trophic collapse, and the decline of ecosystem services. Ecology 87: 1915–1924.
- 27. Lotze HK, Milewski I (2004) Two centuries of multiple human impacts and successive changes in a North Atlantic food web. Ecological Applications 14: 1428–1447.
- 28. Gamfeldt L, Hillebrand H, Jonsson PR (2005) Species richness changes across two trophic levels simultaneously affect prey and consumer biomass. Ecology Letters 8: 696–703.
- 29. Thebault E, Loreau M (2006) The relationship between biodiversity and ecosystem functioning in food webs. Ecological Research 21: 17–25.
- 30. Aquilino KM, Cardinale BJ, Ives AR (2005) Reciprocal effects of host plant and natural enemy diversity on herbivore suppression: an empirical study of a model tritrophic system. Oikos 108: 275–282.
- 31. Long ZT, Bruno JF, Duffy JE (2007) Biodiversity mediates productivity through different mechanisms at adjacent trophic levels. Ecology 88: 2821–2829.
- 32. Stachowicz JJ, Bruno JF, Duffy JE (2007) Understanding the effects of marine biodiversity on communities and ecosystems. Annual Review of Ecology, Evolution, and Systematics 38: 739–766.
- 33. Cardinale BJ, Matulich KL, Hooper DU, Byrnes JE, Duffy JE, et al. (2011) The functional role of producer diversity in ecosystems. American Journal of Botany 98: 572–592.
- 34. Schmid B, Balvanera P, Cardinale BJ, Godbold J, Pfisterer AB, et al. (2009) Consequences of species loss for ecosystem functioning: Meta-analyses of data from biodiversity experiments. In: Naeem S, Bunker DE, Hector A, Loreau M, Perrings C, editors. Biodiversity, ecosystem functioning, and human wellbeing: An ecological and economic perspective. Oxford, UK: Oxford University Press. pp. 14–29.
- 35. Cardinale B, Duffy E, Srivastava D, Loreau M, Thomas M, et al. (2009) Towards a food web perspective on biodiversity and ecosystem functioning. In: Naeem S, Bunker DE, Hector A, Loreau M, Perrings C, editors. Biodiversity, ecosystem functioning, and human wellbeing: An ecological and economic perspective. Oxford, UK: Oxford University Press. pp. 105–120.
- 36. Downing AL, Leibold MA (2002) Ecosystem consequences of species richness and composition in pond food webs. Nature 416: 837–841.
- 37. Bruno JF, Boyer KE, Duffy JE, Lee SC (2008) Relative and interactive effects of plant and grazer richness in a benthic marine community. Ecology 89: 2518–2528.
- 38. Naeem S, Thompson LJ, Lawler SP, Lawton JH, Woodfin RM (1994) Declining biodiversity can alter the performance of ecosystems. Nature 368: 734–737.
- 39. Srivastava DS, Cardinale BJ, Downing AL, Duffy JE, Jouseau C, et al. (2009) Diversity has stronger top-down than bottom-up effects on decomposition. Ecology 90: 1073–1083.
- 40. Ruzycki JR, Beauchamp DA, Yule DL (2003) Effects of introduced lake trout on native cutthroat trout in Yellowstone Lake. Ecological Applications 13: 23–37.
- 41. Dyer LA, Letourneau D (2003) Top-down and bottom-up diversity cascades in detrital vs. living food webs. Ecology Letters 6: 60–68.
- 42. Hunter MD, Price PW (1992) Playing chutes and ladders: heterogeneity and the relative roles of bottom-up and top-down forces in natural communities. Ecology 73:
- 43. Hairston NG, Smith FE, Slobodkin LB (1960) Community structure, population control, and competition. The American Naturalist 94: 425.
- 44. Douglass JG, Duffy JE, Bruno JF (2008) Herbivore and predator diversity interactively affect ecosystem properties in an experimental marine community. Ecology Letters 11: 598–608.
- 45. Bruno JF, O’Connor MI (2005) Cascading effects of predator diversity and omnivory in a marine food web. Ecology Letters 8: 1048–1056.
- 46. Bruno JF, Boyer KE, Duffy JE, Lee SC, Kertesz JS (2005) Effects of macroalgal species identity and richness on primary production in benthic marine communities. Ecology Letters 8: 1165–1174.
- 47. Duffy JE (2006) Biodiversity and the functioning of seagrass ecosystems. Marine Ecology Progress Series 311: 233–250.
- 48. Stachowicz JJ, Graham M, Bracken MES, Szoboszlai AI (2008) Diversity enhances cover and stability of seaweed assemblages: the role of heterogeneity and time. Ecology 89: 3008–3019.
- 49. Tilman D, Reich PB, Knops J, Wedin D, Mielke T, et al. (2001) Diversity and productivity in a long-term grassland experiment. Science 294: 843–845.
- 50. Hooper DU, Chapin ES, III , Ewel JJ, Hector A, Inchausti P, et al. (2005) Effects of biodiversity on ecosystem functioning: A consensus of current knowledge. Ecological Monographs 75: 3–35.
- 51. Cardinale BJ, Srivastava DS, Duffy JE, Wright JP, Downing AL, et al. (2006) Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature 443: 989–992.
- 52. Loreau M, Naeem S, Inchausti P, Bengtsson J, Grime JP, et al. (2001) Biodiversity and ecosystem functioning: current knowledge and future challenges. Science 294: 804–808.
- 53. Fox JW (2004) Effects of algal and herbivore diversity on the partitioning of biomass within and among trophic levels. Ecology 85: 549–559.
- 54. Naeem S, Knops JMH, Tilman D, Howe KM, Kennedy T, et al. (2000) Plant diversity increases resistance to invasion in the absence of covarying extrinsic factors. Oikos 91: 97–108.
- 55. Bruno JF, Cardinale BJ (2008) Cascading effects of predator richness. Frontiers in Ecology and the Environment 6: 539–546.
- 56. Edwards KE, Aquilino KM, Best RJ, Sellheim KL, Stachowicz JJ (2009) Prey diversity is associated with weaker consumer effects in a meta-analysis of benthic marine experiments. Ecology Letters 13: 194–201.
- 57. Leibold MA (1989) Resource edibility and the effects of predators and productivity on consumer-resource interactions. American Naturalist 134: 922–949.
- 58. Duffy JE (2002) Biodiversity and ecosystem function: The consumer connection. Oikos 99: 201–219.
- 59. Root RB (1973) Organization of a plant-arthropod association in simple and diverse habitats: the fauna of collards (Brassica oleracea). Ecological Monographs. pp. 95–124.
- 60. Jolliffe PA (2000) The replacement series. Journal of Ecology 88: 371–385.
- 61. Boyer KE, Kertesz JS, Bruno JF (2009) Environmental context influences the effects of species richness on productivity but not stability of marine macroalgal communities. Oikos 118: 1062–1072.
- 62. Stachowicz JJ, Best RJ, Bracken MES, Graham MH (2008) Complementarity in marine biodiversity manipulations: Reconciling divergent evidence from field and mesocosm experiments. PNAS 105: 18842–18847.
- 63. Griffin JN, Jenkins SR, Gamfeldt L, Jones D, Hawkins SJ, et al. (2009) Spatial heterogeneity increases the importance of species richness for an ecosystem process. Oikos 118: 1335–1342.
- 64. Cruz-Rivera E, Hay ME (2001) Macroalgal traits and the feeding and fitness of an herbivorous amphipod: The roles of selectivity, mixing, and compensation. Marine Ecology Progress Series 218: 249–266.
- 65. Cruz-Rivera E, Hay ME (2003) Prey nutritional quality interacts with chemical defenses to affect consumer feeding and fitness. Ecological Monographs 73: 483–506.
- 66. Haddad NM, Tilman D, Haarstad J, Ritchie M, Knops JMH (2001) Contrasting effects of plant richness and composition on insect communities: a field experiment. The American Naturalist 158: 17–35.
- 67. Scherber C, Mulcu A, Partsch S, Scheu S, Weisser WW (2006) The effects of plant diversity and insect herbivory on performance of individual plant species in experimental grassland. Journal of Ecology 94: 922–931.
- 68. Mulder CPH, Koricheva J, Huss-Danell K, Högberg P, Joshi J (1999) Insects affect relationships between plant species richness and ecosystem processes. Ecology Letters 2: 237–246.
- 69. DeMott WR (1998) Utilization of a cyanobacterium and a phosphorus-deficient green alga as complementary resources by daphnids. Ecology 79: 2463–2481.
- 70. Siemann E, Tilman D, Haartstad J, Ritchie M (1998) Experimental tests of the dependence of arthropod diversity on plant diversity. American Naturalist 152: 738–750.
- 71. Moran E, Reynolds PL, Ladwig LM, O’Connor MI, Long ZT, et al. (2010) Predation intensity is negatively related to plant species richness in a benthic marine community. Marine Ecology Progress Series 400: 277–282.
- 72. Long ZT, Bruno JF, Duffy EJ (2011) Food chain length and omnivory determine the stability of a marine subtidal food web. Journal of Animal Ecology 80: 586–594.
- 73. Halpern BS, Borer ET, Seabloom EW, Shurin JB (2005) Predator effects on herbivore and plant stability. Ecology Letters 8: 189–194.
- 74. Burkepile DE, Hay ME (2008) Herbivore species richness and feeding complementarity affect community structure and function on a coral reef. Proceedings of the National Academy of Sciences of the United States of America 105: 16201–16206.
- 75. Byrnes J, Stachowicz JJ, Hultgren KM, Hughes AR, Olyarnik SV, et al. (2006) Predator diversity strengthens trophic cascades in kelp forests by modifying herbivore behaviour. Ecology Letters 9: 61–71.
- 76. Finke DL, Denno RF (2004) Predator diversity dampens trophic cascades. Nature 429: 407–410.
- 77. Duffy JE, Hay ME (1994) Herbivore resistance to seaweed chemical defense: The roles of mobility and predation risk. Ecology 75: 1304–1319.
- 78. Holmlund MB, Peterson CH, Hay ME (1990) Does algal morphology affect amphipod susceptibility to fish predation? Journal of Experimental Marine Biology and Ecology 139: 65–84.
- 79. Hoffmann M, Hilton-Taylor C, Angulo A, Bohm M, Brooks TM, et al. (2010) The impact of conservation on the status of the world’s vertebrates. Science 330: 1503–1509.
- 80. Pereira HM, Leadley PW, Proenca V, Alkemade R, Scharlemann JPW, et al. (2010) Scenarios for global biodiversity in the 21st century. Science 330: 1496–1501.
- 81. Nelson WG (1979) An analysis of structural pattern in an eelgrass (Zostera marina L.) amphipod community. Journal of Experimental Marine Biology and Ecology 39: 231–264.
- 82. Powers SP, Peterson CH, Christian RR, Sullivan E, Powers MJ, et al. (2005) Effects of eutrophication on bottom habitat and prey resources of demersal fishes. Marine Ecology Progress Series 302: 233–243.
- 83. O’Connor MI, Bruno JF (2009) Predator richness has no effect in a diverse marine food web. Journal of Animal Ecology 78: 732–740.
- 84. Tilman D (1997) Distinguishing between the effects of species diversity and species composition. Oikos 80: 185.
- 85. Sih A, Englund G, Wooster D (1998) Emergent impacts of multiple predators on prey. Trends in Ecology and Evolution 13: 350–355.
- 86. Griffin JN, De la Haye KL, Hawkins SJ, Thompson RC, Jenkins SR (2008) Predator diversity and ecosystem functioning: Density modifies the effect of resource partitioning. Ecology 89: 298–305.
- 87. Byrnes JE, Stachowicz JJ (2009) The consequences of consumer diversity loss: Different answers from different experimental designs. Ecology 90: 2879–2888.
- 88. Duffy JE, Hay ME (2000) Strong impacts of grazing amphipods on the organization of a benthic community. Ecological Monographs 70: 237–263.
- 89. Reynolds PL, Sotka EE (2011) Non-consumptive predator effects indirectly influence marine plant biomass and palatability. Journal of Ecology 99: 1272–1281.
- 90. Hay ME, Sutherland JP (1988) The ecology of rubble structures of the South Atlantic Bight: A community profile. In: Washington , editor. D.C.: National Wetlands Research Center. 85(7.20) p.
- 91. Jernakoff P, Brearley A, Nielsen J (1996) Factors affecting grazer-epiphyte interactions in temperate seagrass meadows. Oceanography and Marine Biology 34: 109–162.
- 92. Nelson WG (1979) Experimental studies of selective predation on amphipods: Consequences for amphipod distribution and abundance. Journal of Experimental Marine Biology and Ecology 38: 225–245.
- 93. Powers SP, Kittinger JN (2002) Hydrodynamic mediation of predator-prey interactions: differential patterns of prey susceptibility and predator success explained by variation in water flow. Journal of Experimental Marine Biology and Ecology 273: 171–187.
- 94. Carpenter SR, Lodge DM (1986) Effects of submerged macrophytes on ecosystem processes. Aquatic Botany 26: 341–370.
- 95. Duffy J (1989) Ecology and evolution of herbivory by marine amphipods: University of North Carolina at Chapel Hill.
- 96. Edgar GJ (1992) Patterns of colonization of mobile epifauna in a Western Australian seagrass bed. Journal of Experimental Marine Biology and Ecology 157: 225–246.
- 97. Virnstein RW, Curran MC (1986) Colonization of artificial seagrass versus time and distance from source. Marine Ecology Progress Series 29: 279–288.
- 98. Lee SC, Bruno JF (2009) Propagule supply controls grazer community structure and primary production in a benthic marine ecosystem. Proceedings of the National Academy of Sciences 106: 7052–7057.
- 99. Duffy JE, Richardson JP, Canuel EA (2003) Grazer diversity effects on ecosystem functioning in seagrass beds. Ecology Letters 6: 637–645.
- 100. Duffy JE, Richardson JP, France KE (2005) Ecosystem consequences of diversity depend on food chain length in estuarine vegetation. Ecology Letters 8: 301–309.
- 101. Muncy RJ (1984) Species profiles: Life histories and environmental requirements of coastal fishes and invertebrates (Gulf of Mexico): Pinfish. U.S. Army Corps of Engineers publication TR ER-82–7. 18 p.
- 102. Stoner AW (1979) Species-specific predation on amphipod crustacea by the pinfish Lagodon rhomboides: Mediation by macrophyte standing crop. Marine Biology 55: 201–207.
- 103. Sotka EE, Reynolds PL (2011) Rapid experimental shifts in host use traits of a polyphagous marine herbivore reveals fitness costs on alternative hosts. Evolutionary Ecology 25: 1335–1355.
- 104. Searle SR, Casella G, McCulloch CE (1992) Variance components. New York: John Wiley & Sons.
- 105. Underwood AJ (1997) Experiments in ecology: Their logistical design and interpretation using analysis of variance. Cambridge, UK: Cambridge University Press.
- 106. Kruskal JB (1964) Nonmetric multidimensional scalng: a numerical method. Psychometrika 29: 115–129.
- 107. Oksanen J, Kindt R, Legendre P, O’Hara B, Simpson GL, et al. (2011) Vegan: community ecology package. R package version 1.17–6. Available: http://vegan.r-forge.r-project.org/. Accessed 4 April 2012.
- 108. Faith DP, Minchin PR, Belbin L (1987) Compositional dissimilarity as a robust measure of ecological distance. Plant Ecology 69: 57–68.
- 109. Levene H (1960) Contributions to probability and statistics: Essays in honor of Harold Hotelling. Stanford University Press.