How habitat-modifying organisms structure the food web of two coastal ecosystems

The diversity and structure of ecosystems has been found to depend both on trophic interactions in food webs and on other species interactions such as habitat modification and mutualism that form non-trophic interaction networks. However, quantification of the dependencies between these two main interaction networks has remained elusive. In this study, we assessed how habitat-modifying organisms affect basic food web properties by conducting in-depth empirical investigations of two ecosystems: North American temperate fringing marshes and West African tropical seagrass meadows. Results reveal that habitat-modifying species, through non-trophic facilitation rather than their trophic role, enhance species richness across multiple trophic levels, increase the number of interactions per species (link density), but decrease the realized fraction of all possible links within the food web (connectance). Compared to the trophic role of the most highly connected species, we found this non-trophic effects to be more important for species richness and of more or similar importance for link density and connectance. Our findings demonstrate that food webs can be fundamentally shaped by interactions outside the trophic network, yet intrinsic to the species participating in it. Better integration of non-trophic interactions in food web analyses may therefore strongly contribute to their explanatory and predictive capacity.     

Systematic variation in food web body-size structure linked to external subsidies

The relationship between body mass (M) and size class abundance (N) depicts patterns of community structure and energy flow through food webs. While the general assumption is that M and N scale linearly (on log–log axes), nonlinearity is regularly observed in natural systems, and is theorized to be driven by nonlinear scaling of trophic level (TL) with M resulting in the rapid transfer of energy to consumers in certain size classes. We tested this hypothesis with data from 31 stream food webs. We predicted that allochthonous subsidies higher in the web results in nonlinear M–TL relationships and systematic abundance peaks in macroinvertebrate and fish size classes (latter containing salmonids), that exploit terrestrial plant material and terrestrial invertebrates, respectively. Indeed, both M–N and M–TL significantly deviated from linear relationships and the observed curvature in M–TL scaling was inversely related to that observed in M–N relationships. Systemic peaks in M–N, and troughs in M–TL occurred in size classes dominated by generalist invertebrates, and brown trout. Our study reveals how allochthonous resources entering high in the web systematically shape community size structure and demonstrates the relevance of a generalized metabolic scaling model for understanding patterns of energy transfer in energetically ‘open’ food webs.     

Exploring the drivers of variation in trophic mismatches: A systematic review of long‐term avian studies

1. Many organisms reproduce in seasonal environments, where selection on timing of reproduction is particularly strong as consumers need to synchronize reproduction with the peaked occurrence of their food. When a consumer species changes its phenology at a slower rate than its resources, this may induce a trophic mismatch, that is, offspring growing up after the peak in food availability, potentially leading to reductions in growth and survival. However, there is large variation in the degree of trophic mismatches as well as in its effects on reproductive output.
2. Here, we explore the potential causes for variation in the strength of trophic mismatches in published studies of birds. Specifically, we ask whether the changes in the degree of mismatch that have occurred over time can be explained by a bird''s (a) breeding latitude, (b) migration distance, and/or (c) life‐history traits.
3. We found that none of these three factors explain changes in the degree of mismatch over time. Nevertheless, food phenology did advance faster at more northerly latitudes, while shifts in bird phenology did not show a trend with latitude.
4. We argue that the lack of support in our results is attributable to the large variation in the metrics used to describe timing of food availability. We propose a pathway to improve the quantification of trophic mismatches, guided by a more rigorous understanding of links between consumers and their resources.

     

Combining food web and species distribution models for improved community projections

The ability to model biodiversity patterns is of prime importance in this era of severe environmental crisis. Species assemblage along environmental gradients is subject to the interplay of biotic interactions in complement to abiotic filtering and stochastic forces. Accounting for complex biotic interactions for a wide array of species remains so far challenging. Here, we propose using food web models that can infer the potential interaction links between species as a constraint in species distribution models. Using a plant–herbivore (butterfly) interaction dataset, we demonstrate that this combined approach is able to improve species distribution and community forecasts. The trophic interaction network between butterfly larvae and host plant was phylogenetically structured and driven by host plant nitrogen content allowing forecasting the food web model to unknown interactions links. This combined approach is very useful in rendering models of more generalist species that have multiple potential interaction links, where gap in the literature may occur. Our combined approach points toward a promising direction for modeling the spatial variation in entire species interaction networks.     

Rapid contemporary evolution and clonal food web dynamics

Character evolution that affects ecological community interactions often occurs contemporaneously with temporal changes in population size, potentially altering the very nature of those dynamics. Such eco-evolutionary processes may be most readily explored in systems with short generations and simple genetics. Asexual and cyclically parthenogenetic organisms such as microalgae, cladocerans and rotifers, which frequently dominate freshwater plankton communities, meet these requirements. Multiple clonal lines can coexist within each species over extended periods, until either fixation occurs or a sexual phase reshuffles the genetic material. When clones differ in traits affecting interspecific interactions, within-species clonal dynamics can have major effects on the population dynamics. We first consider a simple predator–prey system with two prey genotypes, parametrized with data from a well-studied experimental system, and explore how the extent of differences in defence against predation within the prey population determine dynamic stability versus instability of the system. We then explore how increased potential for evolution affects the community dynamics in a more general community model with multiple predator and multiple prey genotypes. These examples illustrate how microevolutionary ‘details’ that enhance or limit the potential for heritable phenotypic change can have significant effects on contemporaneous community-level dynamics and the persistence and coexistence of species.     

Interactions between meat intake and genetic variation in relation to colorectal cancer

Meat intake is associated with the risk of colorectal cancer. The objective of this systematic review was to evaluate interactions between meat intake and genetic variation in order to identify biological pathways involved in meat carcinogenesis. We performed a literature search of PubMed and Embase using “interaction”, “meat”, “polymorphisms”, and “colorectal cancer”, and data on meat–gene interactions were extracted. The studies were divided according to whether information on meat intake was collected prospectively or retrospectively. In prospective studies, interactions between meat intake and polymorphisms in PTGS2 (encoding COX-2), ABCB1, IL10, NFKB1, MSH3, XPC (P_int = 0.006, 0.01, 0.04, 0.03, 0.002, 0.01, respectively), but not IL1B, HMOX1, ABCC2, ABCG2, NR1I2 (encoding PXR), NR1H2 (encoding LXR), NAT1, NAT2, MSH6, or MLH1 in relation to CRC were found. Interaction between a polymorphism in XPC and meat was found in one prospective and one case–control study; however, the directions of the risk estimates were opposite. Thus, none of the findings were replicated. The results from this systematic review suggest that genetic variation in the inflammatory response and DNA repair pathway is involved in meat-related colorectal carcinogenesis, whereas no support for the involvement of heme and iron from meat or cooking mutagens was found. Further studies assessing interactions between meat intake and genetic variation in relation to CRC in large well-characterised prospective cohorts with relevant meat exposure are warranted.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-014-0448-9) contains supplementary material, which is available to authorized users.     

Gillespie eco‐evolutionary models (GEMs) reveal the role of heritable trait variation in eco‐evolutionary dynamics

Heritable trait variation is a central and necessary ingredient of evolution. Trait variation also directly affects ecological processes, generating a clear link between evolutionary and ecological dynamics. Despite the changes in variation that occur through selection, drift, mutation, and recombination, current eco‐evolutionary models usually fail to track how variation changes through time. Moreover, eco‐evolutionary models assume fitness functions for each trait and each ecological context, which often do not have empirical validation. We introduce a new type of model, Gillespie eco‐evolutionary models (GEMs), that resolves these concerns by tracking distributions of traits through time as eco‐evolutionary dynamics progress. This is done by allowing change to be driven by the direct fitness consequences of model parameters within the context of the underlying ecological model, without having to assume a particular fitness function. GEMs work by adding a trait distribution component to the standard Gillespie algorithm – an approach that models stochastic systems in nature that are typically approximated through ordinary differential equations. We illustrate GEMs with the Rosenzweig–MacArthur consumer–resource model. We show not only how heritable trait variation fuels trait evolution and influences eco‐evolutionary dynamics, but also how the erosion of variation through time may hinder eco‐evolutionary dynamics in the long run. GEMs can be developed for any parameter in any ordinary differential equation model and, furthermore, can enable modeling of multiple interacting traits at the same time. We expect GEMs will open the door to a new direction in eco‐evolutionary and evolutionary modeling by removing long‐standing modeling barriers, simplifying the link between traits, fitness, and dynamics, and expanding eco‐evolutionary treatment of a greater diversity of ecological interactions. These factors make GEMs much more than a modeling advance, but an important conceptual advance that bridges ecology and evolution through the central concept of heritable trait variation.     

Prey and predator density‐dependent interactions under different water volumes

Predation is a critical ecological process that directly and indirectly mediates population stabilities, as well as ecosystem structure and function. The strength of interactions between predators and prey may be mediated by multiple density dependences concerning numbers of predators and prey. In temporary wetland ecosystems in particular, fluctuating water volumes may alter predation rates through differing search space and prey encounter rates. Using a functional response approach, we examined the influence of predator and prey densities on interaction strengths of the temporary pond specialist copepod Lovenula raynerae preying on cladoceran prey, Daphnia pulex, under contrasting water volumes. Further, using a population dynamic modeling approach, we quantified multiple predator effects across differences in prey density and water volume. Predators exhibited type II functional responses under both water volumes, with significant antagonistic multiple predator effects (i.e., antagonisms) exhibited overall. The strengths of antagonistic interactions were, however, enhanced under reduced water volumes and at intermediate prey densities. These findings indicate important biotic and abiotic contexts that mediate predator–prey dynamics, whereby multiple predator effects are contingent on both prey density and search area characteristics. In particular, reduced search areas (i.e., water volumes) under intermediate prey densities could enhance antagonisms by heightening predator–predator interference effects.     

Evolution of parental care driven by mutual reinforcement of parental food provisioning and sibling competition

In mammals, altricial birds and some invertebrates, parents care for their offspring by providing them with food and protection until independence. Although parental food provisioning is often essential for offspring survival and growth, very little is known about the conditions favouring the evolutionary innovation of this key component of care. Here, we develop a mathematical model for the evolution of parental food provisioning. We find that this evolutionary innovation is favoured when the efficiency of parental food provisioning is high relative to the efficiency of offspring self-feeding and/or parental guarding. We also explore the coevolution between food provisioning and other components of parental care, as well as offspring behaviour. We find that the evolution of food provisioning prompts evolutionary changes in other components of care by allowing parents to choose safer nest sites, and that it promotes the evolution of sibling competition, which in turn further drives the evolution of parental food provisioning. This mutual reinforcement of parental care and sibling competition suggests that evolution of parental food provisioning should show a unidirectional trend from no parental food provisioning to full parental food provisioning.     

Impacts of consumer–resource interaction transitions on persistence and long‐term interaction outcomes of random ecological networks

Despite the prevalence of context‐dependent interaction transitions in ecological systems, their impacts on persistence and interaction diversity have scarcely been explored in complex ecological networks. By using multispecies bi‐directional and unidirectional consumer–resource models, representing a continuum of interaction transitions (sign change of interaction outcomes), we investigated the effects of structural interaction transitions on persistence (the fraction of remaining species) and long‐term interaction outcomes in random ecological networks. We found that high interaction strength of exploiting resources generally decreased persistence, and high strength of providing resources increased persistence when the strength of exploiting resources was low in more complex networks; also, the networks with high persistence had a high proportion of mutualistic interactions relative to antagonistic interactions present initially and over the long term. The shifting of interaction strengths shaped the long‐term interaction compositions. Meanwhile, population dynamics, especially species extinction, affected the difference between initial and long‐term interactions. Based on classical consumer–resource theory, these results establish a transitional continuum of interaction outcomes in ecological networks and imply a theoretical association among interaction transition, community persistence and interaction diversity.     

Warming shifts top-down and bottom-up control of pond food web structure and function

The effects of global and local environmental changes are transmitted through networks of interacting organisms to shape the structure of communities and the dynamics of ecosystems. We tested the impact of elevated temperature on the top-down and bottom-up forces structuring experimental freshwater pond food webs in western Canada over 16 months. Experimental warming was crossed with treatments manipulating the presence of planktivorous fish and eutrophication through enhanced nutrient supply. We found that higher temperatures produced top-heavy food webs with lower biomass of benthic and pelagic producers, equivalent biomass of zooplankton, zoobenthos and pelagic bacteria, and more pelagic viruses. Eutrophication increased the biomass of all organisms studied, while fish had cascading positive effects on periphyton, phytoplankton and bacteria, and reduced biomass of invertebrates. Surprisingly, virus biomass was reduced in the presence of fish, suggesting the possibility for complex mechanisms of top-down control of the lytic cycle. Warming reduced the effects of eutrophication on periphyton, and magnified the already strong effects of fish on phytoplankton and bacteria. Warming, fish and nutrients all increased whole-system rates of net production despite their distinct impacts on the distribution of biomass between producers and consumers, plankton and benthos, and microbes and macrobes. Our results indicate that warming exerts a host of indirect effects on aquatic food webs mediated through shifts in the magnitudes of top-down and bottom-up forcing.     

Drought has negative consequences on aphid fitness and plant vigor: Insights from a meta‐analysis

1. Aphids are abundant in natural and managed vegetation, supporting a diverse community of organisms and causing damage to agricultural crops. Due to a changing climate, periods of drought are anticipated to increase, and the potential consequences of this for aphid–plant interactions are unclear.
2. Using a meta‐analysis and synthesis approach, we aimed to advance understanding of how increased drought incidence will affect this ecologically and economically important insect group and to characterize any potential underlying mechanisms. We used qualitative and quantitative synthesis techniques to determine whether drought stress has a negative, positive, or null effect on aphid fitness and examined these effects in relation to (a) aphid biology, (b) geographical region, and (c) host plant biology.
3. Across all studies, aphid fitness is typically reduced under drought. Subgroup analysis detected no difference in relation to aphid biology, geographical region, or the aphid–plant combination, indicating the negative effect of drought on aphids is potentially universal. Furthermore, drought stress had a negative impact on plant vigor and increased plant concentrations of defensive chemicals, suggesting the observed response of aphids is associated with reduced plant vigor and increased chemical defense in drought‐stressed plants.
4. We propose a conceptual model to predict drought effects on aphid fitness in relation to plant vigor and defense to stimulate further research.

     

Ecological variation in South American geophagine cichlids arose during an early burst of adaptive morphological and functional evolution

Diversity and disparity are unequally distributed both phylogenetically and geographically. This uneven distribution may be owing to differences in diversification rates between clades resulting from processes such as adaptive radiation. We examined the rate and distribution of evolution in feeding biomechanics in the extremely diverse and continentally distributed South American geophagine cichlids. Evolutionary patterns in multivariate functional morphospace were examined using a phylomorphospace approach, disparity-through-time analyses and by comparing Brownian motion (BM) and adaptive peak evolutionary models using maximum likelihood. The most species-rich and functionally disparate clade (CAS) expanded more efficiently in morphospace and evolved more rapidly compared with both BM expectations and its sister clade (GGD). Members of the CAS clade also exhibited an early burst in functional evolution that corresponds to the development of modern ecological roles and may have been related to the colonization of a novel adaptive peak characterized by fast oral jaw mechanics. Furthermore, reduced ecological opportunity following this early burst may have restricted functional evolution in the GGD clade, which is less species-rich and more ecologically specialized. Patterns of evolution in ecologically important functional traits are consistent with a pattern of adaptive radiation within the most diverse clade of Geophagini.     

Revolution in food web analysis and trophic ecology: diet analysis by DNA and stable isotope analysis

Characterization of energy flow in ecosystems is one of the primary goals of ecology, and the analysis of trophic interactions and food web dynamics is key to quantifying energy flow. Predator‐prey interactions define the majority of trophic interactions and food web dynamics, and visual analysis of stomach, gut or fecal content composition is the technique traditionally used to quantify predator‐prey interactions. Unfortunately such techniques may be biased and inaccurate due to variation in digestion rates ( Sheppard & Hardwood 2005 ); however, those limitations can be largely overcome with new technology. In the last 20 years, the use of molecular genetic techniques in ecology has exploded ( King et al. 2008 ). The growing availability of molecular genetic methods and data has fostered the use of PCR‐based techniques to accurately distinguish and identify prey items in stomach, gut and fecal samples. In this month’s issue of Molecular Ecology Resources, Corse et al. (2010) describe and apply a new approach to quantifying predator‐prey relationships using an ecosystem‐level genetic characterization of available and consumed prey in European freshwater habitats ( Fig. 1a ). In this issue of Molecular Ecology, Hardy et al. (2010) marry the molecular genetic analysis of prey with a stable isotope (SI) analysis of trophic interactions in an Australian reservoir community ( Fig. 1b ). Both papers demonstrate novel and innovative approaches to an old problem – how do we effectively explore food webs and energy movement in ecosystems?

Figure 1 Open in figure viewer PowerPoint The aquatic habitats used for two studies of diet and trophic interactions that employed molecular genetic and stable isotope analyses. Panel a: Example of Rhone basin habitat (France) where fish diet was determined using PCR to classify prey to a series of ecological clades (photo by Emmanuel Corse). Panel b: A weir pool on the lower Murray River (Australia) where food web and prey use was evaluated using a combination of advanced molecular genetic and stable isotope analyses (photo credit: CSIRO).     

Long-term variation of link strength in a simple benthic food web

1. The predatory isopod Saduria entomon (L.) and its amphipod prey Monoporeia affinis (Lindstr?m) are key components of the food web in the northern Baltic Sea, together representing 80-90% of the macrobenthic biomass. We use 20 years of stomach content data for Saduria to investigate how diet dynamics affect the stability of the interaction between Saduria and Monoporeia. 2. Consumption of the main prey, Monoporeia, fitted a type III functional response. Consumption rates of the most important alternative prey, mysids, were found to be unrelated to mysid densities but negatively related to the density of Monoporeia. The fit of consumption data to a model that assumes passive prey selection was poor. Thus we conclude that some form of active choice is involved. 3. The effect of consumption of mysids, the alternative prey, on the stability of this system was investigated using a 'one predator-two prey' model with stochastic environmental variation. Analysis of the model suggests that feeding on mysids leads to a decreased extinction risk for the predator, Saduria, and reduced density oscillations for both Saduria and its main prey, Monoporeia.     

Offspring dynamics affect food provisioning,growth and mortality in a brood-caring spider

In brood-caring species, family members are faced with a conflict over resource distribution. While parents are selected to adapt the amount of care according to their offspring''s needs, offspring might be selected to demand more care than optimal for parents. Recent studies on birds have shown that the social network structure of offspring affects the amount of care and thus the fitness of families. Such a network structure of repeated interactions is probably influenced by within-brood relatedness. We experimentally manipulated the group composition in a brood-caring spider to test how the presence of unrelated spiderlings affects the dynamics between female and brood as well as within broods. Broods consisting of siblings grew better and had a lower mortality compared with mixed broods, no matter whether the caring female was a genetic or foster mother. Interestingly, we found that foster mothers lost weight when caring for sibling broods, whereas females caring for mixed broods gained weight. This indicates that females may be willing to share more prey when the brood contains exclusively siblings even if the entire brood is unrelated to the female. Resource distribution may thus be negotiated by offspring dynamics that could have a signalling function to females.     

Long-term variation in influenza A virus prevalence and subtype diversity in migratory mallards in northern Europe

Data on long-term circulation of pathogens in wildlife populations are seldom collected, and hence understanding of spatial–temporal variation in prevalence and genotypes is limited. Here, we analysed a long-term surveillance series on influenza A virus (IAV) in mallards collected at an important migratory stopover site from 2002 to 2010, and characterized seasonal dynamics in virus prevalence and subtype diversity. Prevalence dynamics were influenced by year, but retained a common pattern for all years whereby prevalence was low in spring and summer, but increased in early autumn with a first peak in August, and a second more pronounced peak during October–November. A total of 74 haemagglutinin (HA)/neuraminidase (NA) combinations were isolated, including all NA and most HA (H1–H12) subtypes. The most common subtype combinations were H4N6, H1N1, H2N3, H5N2, H6N2 and H11N9, and showed a clear linkage between specific HA and NA subtypes. Furthermore, there was a temporal structuring of subtypes within seasons based on HA phylogenetic relatedness. Dissimilar HA subtypes tended to have different temporal occurrence within seasons, where the subtypes that dominated in early autumn were rare in late autumn, and vice versa. This suggests that build-up of herd immunity affected IAV dynamics in this system.     

Archaea in metazoan diets: implications for food webs and biogeochemical cycling

Although the importance of trophic linkages, including ‘top-down forcing'', on energy flow and ecosystem productivity is recognized, the influence of metazoan grazing on Archaea and the biogeochemical processes that they mediate is unknown. Here, we test if: (1) Archaea provide a food source sufficient to allow metazoan fauna to complete their life cycle; (2) neutral lipid biomarkers (including crocetane) can be used to identify Archaea consumers; and (3) archaeal aggregates are a dietary source for methane seep metazoans. In the laboratory, we demonstrated that a dorvilleid polychaete, Ophryotrocha labronica, can complete its life cycle on two strains of Euryarchaeota with the same growth rate as when fed bacterial and eukaryotic food. Archaea were therefore confirmed as a digestible and nutritious food source sufficient to sustain metazoan populations. Both strains of Euryarchaeota used as food sources had unique lipids that were not incorporated into O. labronica tissues. At methane seeps, sulfate-reducing bacteria that form aggregations and live syntrophically with anaerobic-methane oxidizing Archaea contain isotopically and structurally unique fatty acids (FAs). These biomarkers were incorporated into tissues of an endolithofaunal dorvilleid polychaete species from Costa Rica (mean bulk δ¹³C=−92±4‰ polar lipids −116‰) documenting consumption of archaeal-bacterial aggregates. FA composition of additional soft-sediment methane seep species from Oregon and California provided evidence that consumption of archaeal-bacterial aggregates is widespread at methane seeps. This work is the first to show that Archaea are consumed by heterotrophic metazoans, a trophic process we coin as ‘archivory''.     

Novel feeding interactions amplify the impact of species redistribution on an Arctic food web

Species are redistributing globally in response to climate warming, impacting ecosystem functions and services. In the Barents Sea, poleward expansion of boreal species and a decreased abundance of Arctic species are causing a rapid borealization of the Arctic communities. This borealization might have profound consequences on the Arctic food web by creating novel feeding interactions between previously non co‐occurring species. An early identification of new feeding links is crucial to predict their ecological impact. However, detection by traditional approaches, including stomach content and isotope analyses, although fundamental, cannot cope with the speed of change observed in the region, nor with the urgency of understanding the consequences of species redistribution for the marine ecosystem. In this study, we used an extensive food web (metaweb) with nearly 2,500 documented feeding links between 239 taxa coupled with a trait data set to predict novel feeding interactions and to quantify their potential impact on Arctic food web structure. We found that feeding interactions are largely determined by the body size of interacting species, although species foraging habitat and metabolic type are also important predictors. Further, we found that all boreal species will have at least one potential resource in the Arctic region should they redistribute therein. During 2014–2017, 11 boreal species were observed in the Arctic region of the Barents Sea. These incoming species, which are all generalists, change the structural properties of the Arctic food web by increasing connectance and decreasing modularity. In addition, these boreal species are predicted to initiate novel feeding interactions with the Arctic residents, which might amplify their impact on Arctic food web structure affecting ecosystem functioning and vulnerability. Under the ongoing species redistribution caused by environmental change, we propose merging a trait‐based approach with ecological network analysis to efficiently predict the impacts of range‐shifting species on food webs.     

Intraspecific trait variation alters the outcome of competition in freshwater ciliates

Trait variation among heterospecific and conspecific organisms may substantially affect community and food web dynamics. While the relevance of competition and feeding traits have been widely studied for different consumer species, studies on intraspecific differences are more scarce, partly owing to difficulties in distinguishing different clones of the same species. Here, we investigate how intraspecific trait variation affects the competition between the freshwater ciliates Euplotes octocarinatus and Coleps hirtus in a nitrogen‐limited chemostat system. The ciliates competed for the microalgae Cryptomonas sp. (Cry) and Navicula pelliculosa (Nav), and the bacteria present in the cultures over a period of 33 days. We used monoclonal Euplotes and three different Coleps clones (Col 1, Col 2, and Col 3) in the experiment that could be distinguished by a newly developed rDNA‐based molecular assay based on the internal transcribed spacer (ITS) regions. While Euplotes feeds on Cry and on bacteria, the Coleps clones cannot survive on bacteria alone but feed on both Cry and Nav with clone‐specific rates. Experimental treatments comprised two‐species mixtures of Euplotes and one or all of the three different Coleps clones, respectively. We found intraspecific variation in the traits “selectivity” and “maximum ingestion rate” for the different algae to significantly affect the competitive outcome between the two ciliate species. As Nav quickly escaped top‐down control and likely reached a state of low food quality, ciliate competition was strongly determined by the preference of different Coleps clones for Cry as opposed to feeding on Nav. In addition, the ability of Euplotes to use bacteria as an alternative food source strengthened its persistence once Cry was depleted. Hence, trait variation at both trophic levels codetermined the population dynamics and the outcome of species competition.