Indirect effects drive coevolution in mutualistic networks

From insect-plant interactions to ecological networks

Abstract Since its inception, biodiversity has largely been understood as species diversity and assessed as such. Interactions among species or functional groups are gradually becoming part of an expanded concept of biodiversity. As a case study of the development of a research program in biodiversity, we summarize our multi-decade studies on interactions of Asteraceae and flowerhead-feeding insects in Brazil. Initially, host species were treated as independent replicates in order to assess the local and turnover components of their herbivore diversity. Research then expanded into sampling entire interactive communities of host plants and their associated herbivores in different localities and regions, enabling new research lines to be pursued. Interaction diversity could be assessed and factored into spatial and among-host components, suggesting a new field of interaction geography. Second, host specialization, a key component of interaction diversity, was reframed considering simultaneously relatedness and local availability of plant hosts. Third, with the influence of complex network theory, community-wide species interactions were probed for topological patterns. Having identified the modular structure of these plant-herbivore systems, later we demonstrated that they fit a compound hierarchical topology, in which interactions are nested within large-scale modules. In a brief survey of research funded by Fapesp, especially within the Biota-Fapesp program, we highlight several lines of internationally recognized research on interaction diversity, notably on plant-frugivore and plant-pollinator interactions, together with new theoretical models. The interplay of field studies with new theoretical and analytical approaches has established interaction diversity as an essential component for monitoring, conserving and restoring biodiversity in its broader sense.

Dietary abundance distributions: Dominance and diversity in vertebrate diets

Diet composition is among the most important yet least understood dimensions of animal ecology. Inspired by the study of species abundance distributions (SADs), we tested for generalities in the structure of vertebrate diets by characterising them as dietary abundance distributions (DADs). We compiled data on 1167 population-level diets, representing >500 species from six vertebrate classes, spanning all continents and oceans. DADs near-universally (92.5%) followed a hollow-curve shape, with scant support for other plausible rank-abundance-distribution shapes. This strong generality is inherently related to, yet incompletely explained by, the SADs of available food taxa. By quantifying dietary generalisation as the half-saturation point of the cumulative distribution of dietary abundance (sp50, minimum number of foods required to account for 50% of diet), we found that vertebrate populations are surprisingly specialised: in most populations, fewer than three foods accounted for at least half the diet. Variation in sp50 was strongly associated with consumer type, with carnivores being more specialised than herbivores or omnivores. Other methodological (sampling method and effort, taxonomic resolution), biological (body mass, frugivory) and biogeographic (latitude) factors influenced sp50 to varying degrees. Future challenges include identifying the mechanisms underpinning the hollow-curve DAD, its generality beyond vertebrates, and the biological determinants of dietary generalisation.

Phylogenetic congruence between Neotropical primates and plants is driven by frugivory

Seed dispersal benefits plants and frugivores, and potentially drives co-evolution, with consequences to diversification evidenced for, e.g., primates. Evidence for macro-coevolutionary patterns in multi-specific, plant-animal mutualisms is scarce, and the mechanisms driving them remain unexplored. We tested for phylogenetic congruences in primate-plant interactions and showed strong co-phylogenetic signals across Neotropical forests, suggesting that both primates and plants share evolutionary history. Phylogenetic congruence between Platyrrhini and Angiosperms was driven by the most generalist primates, modulated by their functional traits, interacting with a wide-range of Angiosperms. Consistently similar eco-evolutionary dynamics seem to be operating irrespective of local assemblages, since co-phylogenetic signal emerged independently across three Neotropical regions. Our analysis supports the idea that macroevolutionary, coevolved patterns among interacting mutualistic partners are driven by super-generalist taxa. Trait convergence among multiple partners within multi-specific assemblages appears as a mechanism favouring these likely coevolved outcomes.

Global trends in the trophic specialisation of flower-visitor networks are explained by current and historical climate

Trophic specialisation is known to vary across space, but the environmental factors explaining such variation remain elusive. Here we used a global dataset of flower-visitor networks to evaluate how trophic specialisation varies between latitudinal zones (tropical and temperate) and across elevation gradients, while considering the environmental variation inherent in these spatial gradients. Specifically, we assessed the role of current (i.e., net primary productivity, temperature, and precipitation) and historical (i.e., temperature and precipitation stability) environmental factors in structuring the trophic specialisation of floral visitors. Spatial variations in trophic specialisation were not explained by latitudinal zones or elevation. Moreover, regardless of network location on the spatial gradient, there was a tendency for higher trophic specialisation in sites with high productivity and precipitation, whereas historical temperature stability was related to lower trophic specialisation. We highlight that both energetic constraints in animal foraging imposed by climate and resource availability may drive the global variation in trophic specialisation.

Pollinator-Mediated Indirect Effects on Plant Fecundity Revealed by Network Indices

AbstractIndirect effects arise when one species influences how another species interacts with a third. Pollinator-mediated indirect effects are widespread in many plant communities and are often not restricted to plant species pairs. An analytical framework does not exist yet that allows for the evaluation of indirect effects through shared pollinators in a community context as well as their consequences for plant fitness. We used network indices describing pollinator sharing to assess the extent to which plant species affect and are affected by others in a pollination network from a species-rich dune community. For 23 plant species, we explore how these indices relate to plant fecundity (seeds/flower) over two years. We further linked plant traits and indices to uncover functional aspects of pollinator-mediated indirect interactions. Species frequently visited by shared pollinators showed higher fecundity and exhibited traits that increase pollinator attraction and generalization. Conversely, species whose shared pollinators frequently visited other plants had lower fecundity and more specialized traits. Thus, pollinator sharing benefited some species while others suffered reproductive disadvantages, consistent with competition. The framework developed here uses network tools to advance our understanding of how pollinator-mediated indirect interactions influence a species' relative reproductive success at the community level.

Consumers' active choice behaviour promotes coevolutionary units in antagonistic networks

Individual behaviour and local context can influence the evolution of ecological interactions and how they structure into networks. In trophic interactions, consumers can increase their fitness by actively choosing resources that they are more likely to explore successfully. Mathematical modelling is often employed in theoretical studies to understand the coevolutionary dynamics between consumers and resources. However, they often disregard the individual consumer behaviour since the complexity of these systems usually requires simplifying assumptions about interaction details. Using an individual-based model, we model a community of several species that interact antagonistically. Each individual has a trait (attack or defence) that is explicitly modelled and the probability of the interaction to occur successfully increases with increased trait-matching. In addition, consumers can actively choose resources that guarantee greater fitness. We show that active consumer choice can generate coevolutionary units over time. It means that the traits of both consumers and resources converge into multiple groups with similar traits and the species interactions stay restricted to these groups over time. We also observed that network structure is more dependent on the parameter that delimits active consumer choice than on the intensity of selective pressure. Thus, our results support the idea that consumer active choice behaviour plays an important role in the ecological and evolutionary processes that structure interacting communities.

The impact of individual variation on abrupt collapses in mutualistic networks

Individual variation is central to species involved in complex interactions with others in an ecological system. Such ecological systems could exhibit tipping points in response to changes in the environment, consequently leading to abrupt transitions to alternative, often less desirable states. However, little is known about how individual trait variation could influence the timing and occurrence of abrupt transitions. Using 101 empirical mutualistic networks, I model the eco‐evolutionary dynamics of such networks in response to gradual changes in strength of co‐evolutionary interactions. Results indicated that individual variation facilitates the timing of transition in such networks, albeit slightly. In addition, individual variation significantly increases the occurrence of large abrupt transitions. Furthermore, topological network features also positively influence the occurrence of such abrupt transitions. These findings argue for understanding tipping points using an eco‐evolutionary perspective to better forecast abrupt transitions in ecological systems.

The disruption of a keystone interaction erodes pollination and seed dispersal networks

Understanding the impacts of global change on ecological communities is a major challenge in modern ecology. The gain or loss of particular species and the disruption of key interactions are both consequences and drivers of global change that can lead to the disassembly of ecological networks. We examined whether the disruption of a hummingbird-mistletoe-marsupial mutualism by the invasion of non-native species can have cascading effects on both pollination and seed dispersal networks in the temperate forest of Patagonia, Argentina. We focused on network motifs, subnetworks composed of a small number of species exhibiting particular patterns of interaction, to examine the structure and diversity of mutualistic networks. We found that the hummingbird-mistletoe-marsupial mutualism plays a critical role in the community by increasing the complexity of pollination and seed dispersal networks through supporting a high diversity of interactions. Moreover, we found that the disruption of this tripartite mutualism by non-native ungulates resulted in diverse indirect effects that led to less complex pollination and seed dispersal networks. Our results demonstrate that the gains and losses of particular species and the alteration of key interactions can lead to cascading effects in the community through the disassembly of mutualistic networks.

The variable effects of global change on insect mutualisms

Insect mutualisms are essential for reproduction of many plants, protection of plants and other insects, and provisioning of nutrients for insects. Disruption of these mutualisms by global change can have important implications for ecosystem processes. Here, we assess the general effects of global change on insect mutualisms, including the possible impacts on mutualistic networks. We find that the effects of global change on mutualisms are extremely variable, making broad patterns difficult to detect. We require studies focusing on changes in cost-benefit ratios, effects of partner dependency, and degree of specialization to further understand how global change will influence insect mutualism dynamics. We propose that rapid coevolution is one avenue by which mutualists can ameliorate the effects of global change.

A keystone mutualism promotes resistance to invasion

It is not uncommon for one or a few species, and their interactions, to have disproportionate effects on other species in ecological communities. Such keystone interactions might affect how communities respond to the invasion of non-native species by preventing or inhibiting the establishment, spread or impact of non-native species. We explore whether a keystone mutualism among a hummingbird-mistletoe-marsupial promotes ecological resistance to an invasive pollinator, the bumblebee Bombus terrestris, by comparing data collected at sites prior to bumblebee invasion to data collected 11 years after the invasion in sites with and without the keystone mutualism. We built pollination networks and focused on network motifs, regarded as building blocks of networks, to identify the central pollinators and estimate the change in their interactions after invasion of B. terrestris. We also estimated the interaction rewiring across the season in post-invasion networks and tested it as a possible mechanism explaining how the keystone mutualism increased ecological resistance to invasion. We found two times more species in post-invasion sites with the keystone mutualism than in post-invasion sites without the keystone mutualism. Moreover, we found that invasive bumblebee reduced the strength and interaction niche of the five central pollinator species while increasing its own strength and interaction niche, suggesting a replacement of interactions. Also, we found that the keystone mutualism promoted resistance to B. terrestris invasion by reducing its negative impacts on central species. In the presence of the keystone mutualism, central species had three times more direct interactions than in sites without this keystone mutualism. The higher interaction rewiring, after invasion of B. terrestris, in sites with the keystone mutualism indicates greater chances of central pollinators to form new interactions and reduces their competence for resources with the non-native bumblebee. Our results demonstrate that a keystone mutualism can enhance community resistance against the impacts of a non-native invasive pollinator by increasing species diversity and promoting interaction rewiring in the community. This study suggests that the conservation of mutualisms, especially those considered keystone, could be essential for long-term preservation of natural communities under current and future impacts of global change.

The joint role of coevolutionary selection and network structure in shaping trait matching in mutualisms

Coevolution can sculpt remarkable trait similarity between mutualistic partners. Yet, it remains unclear which network topologies and selection regimes enhance trait matching. To address this, we simulate coevolution in topologically distinct networks under a gradient of mutualistic selection strength. We describe three main insights. First, trait matching is jointly influenced by the strength of mutualistic selection and the structural properties of the network where coevolution is unfolding. Second, the strength of mutualistic selection determines the network descriptors better correlated with higher trait matching. While network modularity enhances trait matching when coevolution is weak, network connectance does so when coevolution is strong. Third, the structural properties of networks outrank those of modules or species in determining the degree of trait matching. Our findings suggest networks can both enhance or constrain trait matching, depending on the strength of mutualistic selection.

Phenotypic and functional variation in venom and venom resistance of two sympatric rattlesnakes and their prey

Predator-prey interactions often lead to the coevolution of adaptations associated with avoiding predation and, for predators, overcoming those defences. Antagonistic coevolutionary relationships are often not simple interactions between a single predator and prey but rather a complex web of interactions between multiple coexisting species. Coevolution between venomous rattlesnakes and small mammals has led to physiological venom resistance in several mammalian taxa. In general, viperid venoms contain large quantities of snake venom metalloproteinase toxins (SVMPs), which are inactivated by SVMP inhibitors expressed in resistant mammals. We explored variation in venom chemistry, SVMP expression, and SVMP resistance across four co-distributed species (California Ground Squirrels, Bryant's Woodrats, Southern Pacific Rattlesnakes, and Red Diamond Rattlesnakes) collected from four different populations in Southern California. Our aim was to understand phenotypic and functional variation in venom and venom resistance in order to compare coevolutionary dynamics of a system involving two sympatric predator-prey pairs to past studies that have focused on single pairs. Proteomic analysis of venoms indicated that these rattlesnakes express different phenotypes when in sympatry, with Red Diamonds expressing more typical viperid venom (with a diversity of SVMPs) and Southern Pacifics expressing a more atypical venom with a broader range of non-enzymatic toxins. We also found that although blood sera from both mammals were generally able to inhibit SVMPs from both rattlesnake species, inhibition depended strongly on the snake population, with snakes from one geographic site expressing SVMPs to which few mammals were resistant. Additionally, we found that Red Diamond venom, rather than woodrat resistance, was locally adapted. Our findings highlight the complexity of coevolutionary relationships between multiple predators and prey that exhibit similar offensive and defensive strategies in sympatry.

Phylogenetic reconstruction of ancestral ecological networks through time for pierid butterflies and their host plants

The study of herbivorous insects underpins much of the theory that concerns the evolution of species interactions. In particular, Pieridae butterflies and their host plants have served as a model system for studying evolutionary arms races. To learn more about the coevolution of these two clades, we reconstructed ancestral ecological networks using stochastic mappings that were generated by a phylogenetic model of host-repertoire evolution. We then measured if, when, and how two ecologically important structural features of the ancestral networks (modularity and nestedness) evolved over time. Our study shows that as pierids gained new hosts and formed new modules, a subset of them retained or recolonised the ancestral host(s), preserving connectivity to the original modules. Together, host-range expansions and recolonisations promoted a phase transition in network structure. Our results demonstrate the power of combining network analysis with Bayesian inference of host-repertoire evolution to understand changes in complex species interactions over time.

Coevolutionary Arms Races and the Conditions for the Maintenance of Mutualism

AbstractEmpirical evidence suggests that coevolutionary arms races between flowering plants and their pollinators can occur in wild populations. In extreme cases, trait escalation may result in evolutionary switching from mutualism to parasitism. However, theoretical approaches to studying coevolution typically assume fixed types of ecological interactions and ignore the evolution of absolute fitness. Here, we introduce a novel approach to track the evolution of absolute fitness as a framework to determine when escalatory coevolution results in a switch from mutualism to parasitism. We apply our approach to two previously studied mechanisms mediating selection as a function of phenotype. Our results demonstrate that interactions mediated by a "bigger-is-better" mechanism evolve toward parasitism. In contrast, generalizing the classical trait-matching mechanism so that the fitness of each species is optimized when trait values mismatch by a particular amount, we find theoretical support for indefinite trait exaggeration that preserves mutualistic interactions. Building on our results, we discuss the consequences of coevolutionary arms races for the maintenance of cheating. Moving beyond pairwise interactions, we consider the ramifications of coevolution in a South African pollination network for the evolution of parasitism. Future work extending our approach beyond pairwise interactions can lead to a framework for understanding the evolution of parasitism in mutualistic networks and further insights into the structure and dynamic nature of ecological communities in general.

Effects of stochasticity on the length and behaviour of ecological transients

There is a growing recognition that ecological systems can spend extended periods of time far away from an asymptotic state, and that ecological understanding will therefore require a deeper appreciation for how long ecological transients arise. Recent work has defined classes of deterministic mechanisms that can lead to long transients. Given the ubiquity of stochasticity in ecological systems, a similar systematic treatment of transients that includes the influence of stochasticity is important. Stochasticity can of course promote the appearance of transient dynamics by preventing systems from settling permanently near their asymptotic state, but stochasticity also interacts with deterministic features to create qualitatively new dynamics. As such, stochasticity may shorten, extend or fundamentally change a system's transient dynamics. Here, we describe a general framework that is developing for understanding the range of possible outcomes when random processes impact the dynamics of ecological systems over realistic time scales. We emphasize that we can understand the ways in which stochasticity can either extend or reduce the lifetime of transients by studying the interactions between the stochastic and deterministic processes present, and we summarize both the current state of knowledge and avenues for future advances.

Phenological traits foster persistence of mutualistic networks by promoting facilitation

Morphological and phenological traits are key determinants of the structure of mutualistic networks. Both traits create forbidden links, but phenological traits can also decouple interaction in time. While such difference likely affects the indirect effects among species and consequently network persistence, it remains overlooked. Here, using a dynamic model, we show that networks structured by phenology favour facilitation over competition within guilds of pollinators and plants, thereby increasing network persistence, while the contrary holds for networks structured by morphology. We further show that such buffering of competition by phenological traits mostly beneficiate to specialists, the most vulnerable species otherwise, which propagate the most positive effects within guilds and promote nestedness. Our results indicate that beyond trophic mismatch, phenological shifts such as those induced by climate change are likely to affect indirect effects within mutualistic assemblages, with consequences for biodiversity.

Disentangling effects of multiple stressors on matter flow in a lake food web

Understanding the relative importance of multiple stressors is valuable to prioritize conservation and restoration measures. Yet, the effects of multiple stressors on ecosystem functioning remain largely unknown in many fresh waters. Here, we provided a methodology combining ecosystem modeling with linear regression to disentangle the effects of multiple stressors on matter flow, an important ecosystem function. Treating a shallow lake as the model ecosystem, we simulated matter flow dynamics during 1950s-2010s with different combinations of stressors by Ecopath with Ecosim (EwE) modeling and determined the relative importance of each stressor by generalized linear mixed models. We found that matter flow of the lake food web was highly dynamic, attributing to effects of multiple anthropogenic stressors. Biological invasion played the strongest role in driving the matter flow dynamics, followed by eutrophication, while biomanipulation (i.e., phytoplankton control by planktivorous fish stocking) was of little importance. Eutrophication had a stronger role on primary producers, pelagic food chain, and top predators, while biological invasion on consumers in the middle of food chains. The former was more important in driving the quantity of matter flow, while the latter on trophic transfer efficiencies. Scenario forecasting showed that reducing nutrients contents could largely shape the matter flow pattern, while biomanipulation had little effect. Our findings provided new insights into understanding the mechanistic links between anthropogenic stressors and ecosystem functioning by combining ecosystem modeling with linear regression.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Network analyses reveal the role of large snakes in connecting feeding guilds in a species-rich Amazonian snake community

In ecological communities, interactions between consumers and resources lead to the emergence of ecological networks and a fundamental problem to solve is to understand which factors shape network structure. Empirical and theoretical studies on ecological networks suggest predator body size is a key factor structuring patterns of interaction. Because larger predators consume a wider resource range, including the prey consumed by smaller predators, we hypothesized that variation in body size favors the rise of nestedness. In contrast, if resource consumption requires specific adaptations, predators are expected to consume distinct sets of resources, thus favoring modularity. We investigate these predictions by characterizing the trophic network of a species-rich Amazonian snake community (62 species). Our results revealed an intricate network pattern resulting from larger species feeding on higher diversity of prey and therefore promoting nestedness, whereas snakes with specific lifestyles and feeding on distinct resources, promoting modularity. Species removal simulations indicated that the nested structure is favored mainly by the presence of five species of the family Boidae, which because of their body size and generalist lifestyles connect modules in the network. Our study highlights the particular ways traits affect the structure of interactions among consumers and resources at the community level.

Genetic and plastic rewiring of food webs under climate change

Climate change is altering ecological and evolutionary processes across biological scales. These simultaneous effects of climate change pose a major challenge for predicting the future state of populations, communities and ecosystems. This challenge is further exacerbated by the current lack of integration of research focused on these different scales.

We propose that integrating the fields of quantitative genetics and food web ecology will reveal new insights on how climate change may reorganize biodiversity across levels of organization. This is because quantitative genetics links the genotypes of individuals to population‐level phenotypic variation due to genetic (G), environmental (E) and gene‐by‐environment (G × E) factors. Food web ecology, on the other hand, links population‐level phenotypes to the structure and dynamics of communities and ecosystems.

We synthesize data and theory across these fields and find evidence that genetic (G) and plastic (E and G × E) phenotypic variation within populations will change in magnitude under new climates in predictable ways. We then show how changes in these sources of phenotypic variation can rewire food webs by altering the number and strength of species interactions, with consequences for ecosystem resilience. We also find evidence suggesting there are predictable asymmetries in genetic and plastic trait variation across trophic levels, which set the pace for phenotypic change and food web responses to climate change. Advances in genomics now make it possible to partition G, E and G × E phenotypic variation in natural populations, allowing tests of the hypotheses we propose.

By synthesizing advances in quantitative genetics and food web ecology, we provide testable predictions for how the structure and dynamics of biodiversity will respond to climate change.

Differential effects of nematode infection on pollinating and non-pollinating fig wasps: Can shared antagonism provide net benefits to a mutualism?

Species pairs that form mutualistic associations are also components of broader organismal community networks. These interaction networks have shaped the evolution of individual mutualisms through interspecific interactions ranging from secondarily mutualistic to intensely antagonistic. Our understanding of this complex context remains limited because characterizing the impacts of species interacting with focal mutualists is often difficult. How is the fitness of mutualists impacted by the co-occurring interactive network of community associates? We investigated this context using a model interaction network comprised of a fig and fig wasp mutualist, eight non-pollinating fig wasp (NPFW) antagonists/commensals and a nematode previously believed to be associated only with the pollinator wasp mutualist. Through repeated sampling and field observations, we characterized the ecological roles of these mutualist-associated organisms to identify key antagonists. We then investigated how potential nematode infection of NPFWs could impact wasp survival across key life stages and, in turn, inferred how this influences the fitness of the fig-pollinator mutualists. Unexpectedly, we found all Ficus petiolaris-associated NPFWs to be the targets for nematode infection, with infection levels sometimes exceeding that of pollinators. Experimental data collected for the most abundant NPFW species suggest that nematode infection significantly reduces their longevity. Further, comparisons of nematode loads for emerging and successfully arriving NPFWs suggest that infection severely limits their dispersal ability. Through these observations, we conclude that this infection could impact NPFWs more severely than either mutualistic partner, suggesting a novel role of density-dependent facultative mutualism between figs, pollinator wasps and the nematode. This antagonist-mediated suppression of other network antagonists may present an ecologically common mechanism through which antagonists can present net benefits for mutualists' fitness.

Epizoochory in Parrots as an Overlooked Yet Widespread Plant-Animal Mutualism

Plant–animal interactions are key to sustaining whole communities and ecosystem function. However, their complexity may limit our understanding of the underlying mechanisms and the species involved. The ecological effects of epizoochory remain little known compared to other seed dispersal mechanisms given the few vectors identified. In addition, epizoochory is mostly considered non-mutualistic since dispersers do not obtain nutritional rewards. Here, we show a widespread but unknown mutualistic interaction between parrots and plants through epizoochory. Combining our observations with photos from web-sources, we recorded nearly 2000 epizoochory events in 48 countries across five continents, involving 116 parrot species and nearly 100 plant species from 35 families, including both native and non-native species. The viscid pulp of fleshy fruits and anemochorous structures facilitate the adherence of tiny seeds (mean 3.7 × 2.56 mm) on the surface of parrots while feeding, allowing the dispersion of these seeds over long distances (mean = 118.5 m). This parrot–plant mutualism could be important in ecosystem functioning across a wide diversity of environments, also facilitating the spread of exotic plants. Future studies should include parrots for a better understanding of plant dispersal processes and for developing effective conservation actions against habitat loss and biological invasions.

Quantification of the covariation of lake microbiomes and environmental variables using a machine learning-based framework

It is known that microorganisms are essential for the functioning of ecosystems, but the extent to which microorganisms respond to different environmental variables in their natural habitats is not clear. In the current study, we present a methodological framework to quantify the covariation of the microbial community of a habitat and environmental variables of this habitat. It is built on theoretical considerations of systems ecology, makes use of state-of-the-art machine learning techniques and can be used to identify bioindicators. We apply the framework to a data set containing operational taxonomic units (OTUs) as well as more than twenty physicochemical and geographic variables measured in a large-scale survey of European lakes. While a large part of variation (up to 61%) in many environmental variables can be explained by microbial community composition, some variables do not show significant covariation with the microbial lake community. Moreover, we have identified OTUs that act as "multitask" bioindicators, i.e., that are indicative for multiple environmental variables, and thus could be candidates for lake water monitoring schemes. Our results represent, for the first time, a quantification of the covariation of the lake microbiome and a wide array of environmental variables for lake ecosystems. Building on the results and methodology presented here, it will be possible to identify microbial taxa and processes that are essential for functioning and stability of lake ecosystems.

An ecological approach to structural flexibility in online communication systems

Human cognitive abilities are limited resources. Today, in the age of cheap information-cheap to produce, to manipulate, to disseminate-this cognitive bottleneck translates into hypercompetition for rewarding outcomes among actors. These incentives push actors to mutualistically interact with specific memes, seeking the virality of their messages. In turn, memes' chances to persist and spread are subject to changes in the communication environment. In spite of all this complexity, here we show that the underlying architecture of empirical actor-meme information ecosystems evolves into recurring emergent patterns. We then propose an ecology-inspired modelling framework, bringing to light the precise mechanisms causing the observed flexible structural reorganisation. The model predicts-and the data confirm-that users' struggle for visibility induces a re-equilibration of the network's mesoscale towards self-similar nested arrangements. Our final microscale insights suggest that flexibility at the structural level is not mirrored at the dynamical one.

A white noise approach to evolutionary ecology

Although the evolutionary response to random genetic drift is classically modelled as a sampling process for populations with fixed abundance, the abundances of populations in the wild fluctuate over time. Furthermore, since wild populations exhibit demographic stochasticity and since random genetic drift is in part due to demographic stochasticity, theoretical approaches are needed to understand the role of demographic stochasticity in eco-evolutionary dynamics. Here we close this gap for quantitative characters evolving in continuously reproducing populations by providing a framework to track the stochastic dynamics of abundance density across phenotypic space using stochastic partial differential equations. In the process we develop a set of heuristics to operationalize the powerful, but abstract theory of white noise and diffusion-limits of individual-based models. Applying these heuristics, we obtain stochastic ordinary differential equations that generalize classical expressions of ecological quantitative genetics. In particular, by supplying growth rate and reproductive variance as functions of abundance densities and trait values, these equations track population size, mean trait and additive genetic variance responding to mutation, demographic stochasticity, random genetic drift, deterministic selection and noise-induced selection. We demonstrate the utility of our approach by formulating a model of diffuse coevolution mediated by exploitative competition for a continuum of resources. In addition to trait and abundance distributions, this model predicts interaction networks defined by niche-overlap, competition coefficients, or selection gradients. Using a high-richness approximation, we find linear selection gradients and competition coefficients are uncorrelated, but magnitudes of linear selection gradients and quadratic selection gradients are both positively correlated with competition coefficients. Hence, competing species that strongly affect each other's abundance tend to also impose selection on one another, but the directionality is not predicted. This approach contributes to the development of a synthetic theory of evolutionary ecology by formalizing first principle derivations of stochastic models tracking feedbacks of biological processes and the patterns of diversity they produce.

Habitat generalist species constrain the diversity of mimicry rings in heterogeneous habitats

How evolution creates and maintains trait patterns in species-rich communities is still an unsolved topic in evolutionary ecology. One classical example of community-level pattern is the unexpected coexistence of different mimicry rings, each of which is a group of mimetic species with the same warning signal. The coexistence of different mimicry rings in a community seems paradoxical because selection among unpalatable species should favor convergence to a single warning pattern. We combined mathematical modeling based on network theory and numerical simulations to explore how different types of selection, such as mimetic and environmental selections, and habitat use by mimetic species influence the formation of coexisting rings. We show that when habitat and mimicry are strong sources of selection, the formation of multiple rings takes longer due to conflicting selective pressures. Moreover, habitat generalist species decrease the distinctiveness of different mimicry rings’ patterns and a few habitat generalist species can generate a “small-world effect”, preventing the formation of multiple mimicry rings. These results may explain why the coexistence of mimicry rings is more common in groups of animals that tend towards habitat specialism, such as butterflies.

A community perspective on the concept of marine holobionts: current status, challenges, and future directions

Host-microbe interactions play crucial roles in marine ecosystems. However, we still have very little understanding of the mechanisms that govern these relationships, the evolutionary processes that shape them, and their ecological consequences. The holobiont concept is a renewed paradigm in biology that can help to describe and understand these complex systems. It posits that a host and its associated microbiota with which it interacts, form a holobiont, and have to be studied together as a coherent biological and functional unit to understand its biology, ecology, and evolution. Here we discuss critical concepts and opportunities in marine holobiont research and identify key challenges in the field. We highlight the potential economic, sociological, and environmental impacts of the holobiont concept in marine biological, evolutionary, and environmental sciences. Given the connectivity and the unexplored biodiversity specific to marine ecosystems, a deeper understanding of such complex systems requires further technological and conceptual advances, e.g., the development of controlled experimental model systems for holobionts from all major lineages and the modeling of (info)chemical-mediated interactions between organisms. Here we propose that one significant challenge is to bridge cross-disciplinary research on tractable model systems in order to address key ecological and evolutionary questions. This first step is crucial to decipher the main drivers of the dynamics and evolution of holobionts and to account for the holobiont concept in applied areas, such as the conservation, management, and exploitation of marine ecosystems and resources, where practical solutions to predict and mitigate the impact of human activities are more important than ever.

Trait positions for elevated invasiveness in adaptive ecological networks

Our ability to predict the outcome of invasion declines rapidly as non-native species progress through intertwined ecological barriers to establish and spread in recipient ecosystems. This is largely due to the lack of systemic knowledge on key processes at play as species establish self-sustaining populations within the invaded range. To address this knowledge gap, we present a mathematical model that captures the eco-evolutionary dynamics of native and non-native species interacting within an ecological network. The model is derived from continuous-trait evolutionary game theory (i.e., Adaptive Dynamics) and its associated concept of invasion fitness which depicts dynamic demographic performance that is both trait mediated and density dependent. Our approach allows us to explore how multiple resident and non-native species coevolve to reshape invasion performance, or more precisely invasiveness, over trait space. The model clarifies the role of specific traits in enabling non-native species to occupy realised opportunistic niches. It also elucidates the direction and speed of both ecological and evolutionary dynamics of residing species (natives or non-natives) in the recipient network under different levels of propagule pressure. The versatility of the model is demonstrated using four examples that correspond to the invasion of (i) a horizontal competitive community; (ii) a bipartite mutualistic network; (iii) a bipartite antagonistic network; and (iv) a multi-trophic food web. We identified a cohesive trait strategy that enables the success and establishment of non-native species to possess high invasiveness. Specifically, we find that a non-native species can achieve high levels of invasiveness by possessing traits that overlap with those of its facilitators (and mutualists), which enhances the benefits accrued from positive interactions, and by possessing traits outside the range of those of antagonists, which mitigates the costs accrued from negative interactions. This ‘central-to-reap, edge-to-elude’ trait strategy therefore describes the strategic trait positions of non-native species to invade an ecological network. This model provides a theoretical platform for exploring invasion strategies in complex adaptive ecological networks.

Frugivory Specialization in Birds and Fruit Chemistry Structure Mutualistic Networks across the Neotropics

AbstractThe interaction between fruit chemistry and the physiological traits of frugivores is expected to shape the structure of mutualistic seed dispersal networks, but it has been understudied compared with the role of morphological trait matching in structuring interaction patterns. For instance, highly frugivorous birds (i.e., birds that have fruits as the main component of their diets), which characteristically have fast gut passage times, are expected to avoid feeding on lipid-rich fruits because of the long gut retention times associated with lipid digestion. Here, we compiled data from 84 studies conducted in the Neotropics that used focal plant methods to record 35,815 feeding visits made by 317 bird species (155 genera in 28 families) to 165 plant species (82 genera in 48 families). We investigated the relationship between the degree of frugivory of birds (i.e., how much of their diet is composed of fruit) at the genus level and their visits to plant genera that vary in fruit lipid content. We used a hierarchical modeling of species communities approach that accounted for the effects of differences in body size, bird and plant phylogeny, and spatial location of study sites. We found that birds with a low degree of frugivory (e.g., predominantly insectivores) tend to have the highest increase in visitation rates as fruits become more lipid rich, while birds that are more frugivorous tend to increase visits at a lower rate or even decrease visitation rates as lipids increase in fruits. This balance between degree of frugivory and visitation rates to lipid-poor and lipid-rich fruits provides a mechanism to explain specialized dispersal systems and the occurrence of certain physiological nutritional filters, ultimately helping us to understand community-wide interaction patterns between birds and plants.

Tricky partners: native plants show stronger interaction preferences than their exotic counterparts

In ecological networks, neutral predictions suggest that species' interaction frequencies are proportional to their relative abundances. Deviations from neutral predictions thus correspond to interaction preferences (when positive) or avoidances (when negative), driven by nonneutral (e.g., niche-based) processes. Exotic species interact with many partners with which they have not coevolved, and it remains unclear whether this systematically influences the strength of neutral processes on interactions, and how these interaction-level differences scale up to entire networks. To fill this gap, we compared interactions between plants and frugivorous birds at nine forest sites in New Zealand varying in the relative abundance and composition of native and exotic species, with independently sampled data on bird and plant abundances from the same sites. We tested if the strength and direction of interaction preferences differed between native and exotic species. We further evaluated whether the performance of neutral predictions at the site level was predicted by the proportion of exotic interactions in each network from both bird and plant perspectives, and the species composition in each site. We found that interactions involving native plants deviated more strongly from neutral predictions than did interactions involving exotics. This "pickiness" of native plants could be detrimental in a context of global biotic homogenization where they could be increasingly exposed to novel interactions with neutrally interacting mutualists. However, the realization of only a subset of interactions in different sites compensated for the neutrality of interactions involving exotics, so that neutral predictions for whole networks did not change systematically with the proportion of exotic species or species composition. Therefore, the neutral and niche processes that underpin individual interactions may not scale up to entire networks. This shows that seemingly simplistic neutral assumptions entail complex processes and can provide valuable understanding of community assembly or invasion dynamics.

The Structure of Ecological Networks Across Levels of Organization

Interactions connect the units of ecological systems, forming networks. Individual-based networks characterize variation in niches among individuals within populations. These individual-based networks merge with each other, forming species-based networks and food webs that describe the architecture of ecological communities. Networks at broader spatiotemporal scales portray the structure of ecological interactions across landscapes and over macroevolutionary time. Here, I review the patterns observed in ecological networks across multiple levels of biological organization. A fundamental challenge is to understand the amount of interdependence as we move from individual-based networks to species-based networks and beyond. Despite the uneven distribution of studies, regularities in network structure emerge across scales due to the fundamental architectural patterns shared by complex networks and the interplay between traits and numerical effects. I illustrate the integration of these organizational scales by exploring the consequences of the emergence of highly connected species for network structures across scales.

The Evolution of Mutualistic Dependence

While the importance of mutualisms across the tree of life is recognized, it is not understood why some organisms evolve high levels of dependence on mutualistic partnerships, while other species remain autonomous or retain or regain minimal dependence on partners. We identify four main pathways leading to the evolution of mutualistic dependence. Then, we evaluate current evidence for three predictions: ( a) Mutualisms with different levels of dependence have distinct stabilizing mechanisms against exploitation and cheating, ( b) less dependent mutualists will return to autonomy more often than those that are highly dependent, and ( c) obligate mutualisms should be less context dependent than facultative ones. Although we find evidence supporting all three predictions, we stress that mutualistic partners follow diverse paths toward—and away from—dependence. We also highlight the need to better examine asymmetry in partner dependence. Recognizing how variation in dependence influences the stability, breakdown, and context dependence of mutualisms generates new hypotheses regarding how and why the benefits of mutualistic partnerships differ over time and space.

A need to consider the evolutionary genetics of host-symbiont mutualisms

Despite the ubiquity and importance of mutualistic interactions, we know little about the evolutionary genetics underlying their long-term persistence. As in antagonistic interactions, mutualistic symbioses are characterized by substantial levels of phenotypic and genetic diversity. In contrast to antagonistic interactions, however, we, by and large, do not understand how this variation arises, how it is maintained, nor its implications for future evolutionary change. Currently, we rely on phenotypic models to address the persistence of mutualistic symbioses, but the success of an interaction almost certainly depends heavily on genetic interactions. In this review, we argue that evolutionary genetic models could provide a framework for understanding the causes and consequences of diversity and why selection may favour processes that maintain variation in mutualistic interactions.

Tipping point and noise-induced transients in ecological networks

A challenging and outstanding problem in interdisciplinary research is to understand the interplay between transients and stochasticity in high-dimensional dynamical systems. Focusing on the tipping-point dynamics in complex mutualistic networks in ecology constructed from empirical data, we investigate the phenomena of noise-induced collapse and noise-induced recovery. Two types of noise are studied: environmental (Gaussian white) noise and state-dependent demographic noise. The dynamical mechanism responsible for both phenomena is a transition from one stable steady state to another driven by stochastic forcing, mediated by an unstable steady state. Exploiting a generic and effective two-dimensional reduced model for real-world mutualistic networks, we find that the average transient lifetime scales algebraically with the noise amplitude, for both environmental and demographic noise. We develop a physical understanding of the scaling laws through an analysis of the mean first passage time from one steady state to another. The phenomena of noise-induced collapse and recovery and the associated scaling laws have implications for managing high-dimensional ecological systems.

Accelerating homogenization of the global plant-frugivore meta-network

Introductions of species by humans are causing the homogenization of species composition across biogeographic barriers^1-3. The ecological and evolutionary consequences of introduced species derive from their effects on networks of species interactions^4,5, but we lack a quantitative understanding of the impacts of introduced species on ecological networks and their biogeographic patterns globally. Here we address this data gap by analysing mutualistic seed-dispersal interactions from 410 local networks, encompassing 24,455 unique pairwise interactions between 1,631 animal and 3,208 plant species. We show that species introductions reduce biogeographic compartmentalization of the global meta-network, in which nodes are species and links are interactions observed within any local network. This homogenizing effect extends across spatial scales, decreasing beta diversity among local networks and modularity within networks. The prevalence of introduced interactions is directly related to human environmental modifications and is accelerating, having increased sevenfold over the past 75 years. These dynamics alter the coevolutionary environments that mutualists experience⁶, and we find that introduced species disproportionately interact with other introduced species. These processes are likely to amplify biotic homogenization in future ecosystems⁷ and may reduce the resilience of ecosystems by allowing perturbations to propagate more quickly and exposing disparate ecosystems to similar drivers. Our results highlight the importance of managing the increasing homogenization of ecological complexity.

Estimating the risk of species interaction loss in mutualistic communities

Interactions between species generate the functions on which ecosystems and humans depend. However, we lack an understanding of the risk that interaction loss poses to ecological communities. Here, we quantify the risk of interaction loss for 4,330 species interactions from 41 empirical pollination and seed dispersal networks across 6 continents. We estimate risk as a function of interaction vulnerability to extinction (likelihood of loss) and contribution to network feasibility, a measure of how much an interaction helps a community tolerate environmental perturbations. Remarkably, we find that more vulnerable interactions have higher contributions to network feasibility. Furthermore, interactions tend to have more similar vulnerability and contribution to feasibility across networks than expected by chance, suggesting that vulnerability and feasibility contribution may be intrinsic properties of interactions, rather than only a function of ecological context. These results may provide a starting point for prioritising interactions for conservation in species interaction networks in the future.

A study of 4,330 species interactions from 41 empirical pollination and seed dispersal networks across six continents reveals that species interactions which are most vulnerable to extinction are also the most important for ecological community stability; moreover, vulnerable interactions that are important for stability tend to be important and vulnerable wherever they occur.

Co-adaptation enhances the resilience of mutualistic networks

Mutualistic networks, which describe the ecological interactions between multiple types of species such as plants and pollinators, play a paramount role in the generation of Earth's biodiversity. The resilience of a mutualistic network denotes its ability to retain basic functionality when errors and failures threaten the persistence of the community. Under the disturbances of mass extinctions and human-induced disasters, it is crucial to understand how mutualistic networks respond to changes, which enables the system to increase resilience and tolerate further damages. Despite recent advances in the modelling of the structure-based adaptation, we lack mathematical and computational models to describe and capture the co-adaptation between the structure and dynamics of mutualistic networks. In this paper, we incorporate dynamic features into the adaptation of structure and propose a co-adaptation model that drastically enhances the resilience of non-adaptive and structure-based adaptation models. Surprisingly, the reason for the enhancement is that the co-adaptation mechanism simultaneously increases the heterogeneity of the mutualistic network significantly without changing its connectance. Owing to the broad applications of mutualistic networks, our findings offer new ways to design mechanisms that enhance the resilience of many other systems, such as smart infrastructures and social-economical systems.

Scoring Species for Synthetic Community Design: Network Analyses of Functional Core Microbiomes

Constructing biological communities is a major challenge in both basic and applied sciences. Although model synthetic communities with a few species have been constructed, designing systems consisting of tens or hundreds of species remains one of the most difficult goals in ecology and microbiology. By utilizing high-throughput sequencing data of interspecific association networks, we here propose a framework for exploring “functional core” species that have great impacts on whole community processes and functions. The framework allows us to score each species within a large community based on three criteria: namely, topological positions, functional portfolios, and functional balance within a target network. The criteria are measures of each species’ roles in maximizing functional benefits at the community or ecosystem level. When species with potentially large contributions to ecosystem-level functions are screened, the framework also helps us design “functional core microbiomes” by focusing on properties of species groups (modules) within a network. When embedded into agroecosystems or human gut, such functional core microbiomes are expected to organize whole microbiome processes and functions. An application to a plant-associated microbiome dataset actually highlighted potential functional core microbes that were known to control rhizosphere microbiomes by suppressing pathogens. Meanwhile, an example of application in mouse gut microbiomes called attention to poorly investigated bacterial species, whose potential roles within gut microbiomes deserve future experimental studies. The framework for gaining “bird’s-eye” views of functional cores within networks is applicable not only to agricultural and medical data but also to datasets produced in food processing, brewing, waste water purification, and biofuel production.

Robustness to extinction and plasticity derived from mutualistic bipartite ecological networks

Understanding the response of ecological networks to perturbations and disruptive events is needed to anticipate the biodiversity loss and extinction cascades. Here, we study how network plasticity reshapes the topology of mutualistic networks in response to species loss. We analyze more than one hundred empirical mutualistic networks and considered random and targeted removal as mechanisms of species extinction. Network plasticity is modeled as either random rewiring, as the most parsimonious approach, or resource affinity-driven rewiring, as a proxy for encoding the phylogenetic similarity and functional redundancy among species. This redundancy should be positively correlated with the robustness of an ecosystem, as functions can be taken by other species once one of them is extinct. We show that effective modularity, i.e. the ability of an ecosystem to adapt or restructure, increases with increasing numbers of extinctions, and with decreasing the replacement probability. Importantly, modularity is mostly affected by the extinction rather than by rewiring mechanisms. These changes in community structure are reflected in the robustness and stability due to their positive correlation with modularity. Resource affinity-driven rewiring offers an increase of modularity, robustness, and stability which could be an evolutionary favored mechanism to prevent a cascade of co-extinctions.

Dynamical network analysis reveals key microRNAs in progressive stages of lung cancer

Non-coding RNAs are fundamental to the competing endogenous RNA (CeRNA) hypothesis in oncology. Previous work focused on static CeRNA networks. We construct and analyze CeRNA networks for four sequential stages of lung adenocarcinoma (LUAD) based on multi-omics data of long non-coding RNAs (lncRNAs), microRNAs and mRNAs. We find that the networks possess a two-level bipartite structure: common competing endogenous network (CCEN) composed of an invariant set of microRNAs over all the stages and stage-dependent, unique competing endogenous networks (UCENs). A systematic enrichment analysis of the pathways of the mRNAs in CCEN reveals that they are strongly associated with cancer development. We also find that the microRNA-linked mRNAs from UCENs have a higher enrichment efficiency. A key finding is six microRNAs from CCEN that impact patient survival at all stages, and four microRNAs that affect the survival from a specific stage. The ten microRNAs can then serve as potential biomarkers and prognostic tools for LUAD.

Persistence and Oscillations of Plant-Pollinator-Herbivore Systems

This paper considers plant-pollinator-herbivore systems where the plant produces food for the pollinator, the pollinator provides pollination service for the plant in return, while the herbivore consumes both the food and the plant itself without providing pollination service. Based on these resource-consumer interactions, we form a plant-pollinator-herbivore model which includes the intermediary food. Using qualitative method and Kuznetsov theorem, we show global dynamics of the subsystems, uniform persistence of the whole system and periodic oscillation by Hopf bifurcation. Rigorous analysis on the system demonstrates mechanisms by which varying parameters could make the system transition between extinction of herbivore, coexistence of the three species at steady states, coexistence in periodic oscillations and extinction of pollinator. It is shown that (i) in plant-pollinator interactions, the plant would produce food; (ii) in plant-herbivore interactions, the plant would produce toxin; (iii) in the presence of both pollinator and herbivore, the plant would produce both food and toxin, and intermediate productions are analytically given by which the plant can reach its maximal density; and (iv) an appropriate toxin production could drive the herbivore into extinction, an unappropriate one would drive the pollinator into extinction, while too much toxin production will drive the plant itself into extinction. The analysis leads to explanations for experimental observations and provides new insights.

Plant ontogeny determines strength and associated plant fitness consequences of plant-mediated interactions between herbivores and flower visitors

Plants show ontogenetic variation in growth-defence strategies to maximize reproductive output within a community context. Most work on plant ontogenetic variation in growth-defence trade-offs has focussed on interactions with antagonistic insect herbivores. Plants respond to herbivore attack with phenotypic changes. Despite the knowledge that plant responses to herbivory affect plant mutualistic interactions with pollinators required for reproduction, indirect interactions between herbivores and pollinators have not been included in the evaluation of how ontogenetic growth-defence trajectories affect plant fitness.In a common garden experiment with the annual Brassica nigra, we investigated whether exposure to various herbivore species on different plant ontogenetic stages (vegetative, bud or flowering stage) affects plant flowering traits, interactions with flower visitors and results in fitness consequences for the plant.Effects of herbivory on flowering plant traits and interactions with flower visitors depended on plant ontogeny. Plant exposure in the vegetative stage to the caterpillar Pieris brassicae and aphid Brevicoryne brassicae led to reduced flowering time and flower production, and resulted in reduced pollinator attraction, pollen beetle colonization, total seed production and seed weight. When plants had buds, infestation by most herbivore species tested reduced flower production and pollen beetle colonization. Pollinator attraction was either increased or reduced. Plants infested in the flowering stage with P. brassicae or Lipaphis erysimi flowered longer, while infestation by any of the herbivore species tested increased the number of flower visits by pollinators.Our results show that the outcome of herbivore-flower visitor interactions in B. nigra is specific for the combination of herbivore species and plant ontogenetic stage. Consequences of herbivory for flowering traits and reproductive output were strongest when plants were attacked early in life. Such differences in selection pressures imposed by herbivores to specific plant ontogenetic stages may drive the evolution of distinct ontogenetic trajectories in growth-defence-reproduction strategies and include indirect interactions between herbivores and flower visitors. Synthesis. Plant ontogeny can define the direct and indirect consequences of herbivory. Our study shows that the ontogenetic stage of plant individuals determined the effects of herbivory on plant flowering traits, interactions with flower visitors and plant fitness.