Indirect effects shape species fitness in coevolved mutualistic networks.

Ecological interactions are one of the main forces that sustain Earth's biodiversity. A major challenge for studies of ecology and evolution is to determine how these interactions affect the fitness of species when we expand from studying isolated, pairwise interactions to include networks of interacting species1-4. In networks, chains of effects caused by a range of species have an indirect effect on other species they do not interact with directly, potentially affecting the fitness outcomes of a variety of ecological interactions (such as mutualism)5-7. Here we apply analytical techniques and numerical simulations to 186 empirical mutualistic networks and show how both direct and indirect effects alter the fitness of species coevolving in these networks. Although the fitness of species usually increased with the number of mutualistic partners, most of the fitness variation across species was driven by indirect effects. We found that these indirect effects prevent coevolving species from adapting to their mutualistic partners and to other sources of selection pressure in the environment, thereby decreasing their fitness. Such decreases are distributed in a predictable way within networks: peripheral species receive more indirect effects and experience higher reductions in fitness than central species. This topological effect was also evident when we analysed an empirical study of an invasion of pollination networks by honeybees. As honeybees became integrated as a central species within networks, they increased the contribution of indirect effects on several other species, reducing their fitness. Our study shows how and why indirect effects can govern the adaptive landscape of species-rich mutualistic assemblages.

The ripple effects of clines from coevolutionary hotspots to coldspots.

Coevolution has the potential to alter not only the ecological interactions of coevolving partners, but also their interactions with yet other species. The effects of coevolution may ripple throughout networks of interacting species, cascading across trophic levels, swamping competitors, or facilitating survival or reproduction of yet other species linked only indirectly to the coevolving partners. These ripple effects of coevolution may differ among communities, amplifying how the coevolutionary process produces geographic mosaics of traits and outcomes in interactions among species. In a From the Cover article in this issue of Molecular Ecology, Hague et al. (2022) provide a clear example, using the well-studied interactions between Pacific newts (Taricha spp.) and their common garter snake (Thamnophis sirtalis) predators in western North America. Pacific newts harbour tetrodotoxin (TTX), which is highly toxic to vertebrate predators. In coevolutionary hotspots, extreme escalation of toxicity in the newts and resistance to toxicity in the snakes have resulted in snake populations that retain high levels of TTX. In two geographic regions, snakes in these hotspot populations have evolved bright, aposematic colours that may act as warning signals to their own vertebrate predators. The warning signals and toxin-resistance alleles in the snake populations decrease clinally away from the coevolutionary hotpots, shaped by a geographically variable mix of selection imposed by the snakes' prey and by their own predators.

Components of local adaptation and divergence in pollination efficacy in a coevolving species interaction.

Selection leading to adaptation to interactions may generate rapid evolutionary feedbacks and drive diversification of species interactions. The challenge is to understand how the many traits of interacting species combine to shape local adaptation in ways directly or indirectly resulting in diversification. We used the well-studied interactions between Lithophragma plants (Saxifragaceae) and Greya moths (Prodoxidae) to evaluate how plants and moths together contributed to local divergence in pollination efficacy. Specifically, we studied L. bolanderi and its two specialized Greya moth pollinators in two contrasting environments in the Sierra Nevada in California. Both moths pollinate L. bolanderi during nectaring, one of them-G. politella-also while ovipositing through the floral corolla into the ovary. First, field surveys of floral visitors and the presence of G. politella eggs and larvae in developing capsules showed that one population was visited only by G. politella and few other pollinators, whereas the other was visited by both Greya species and other pollinators. Second, L. bolanderi in these two natural populations differed in several floral traits putatively important for pollination efficacy. Third, laboratory experiments with greenhouse-grown plants and field-collected moths showed that L. bolanderi was more efficiently pollinated by local compared to nonlocal nectaring moths of both species. Pollination efficacy of ovipositing G. politella was also higher for local moths for the L. bolanderi population, which relies more heavily on this species in nature. Finally, time-lapse photography in the laboratory showed that G. politella from different populations differed in oviposition behavior, suggesting the potential for local adaptation also among Greya populations. Collectively, our results are a rare example of components of local adaptation contributing to divergence in pollination efficacy in a coevolving interaction and, thus, provide insights into how geographic mosaics of coevolution may lead to coevolutionary diversification in species interactions.

Genetic correlations and ecological networks shape coevolving mutualisms.

Ecological interactions shape the evolution of multiple species traits in populations. These traits are often linked to each other through genetic correlations, affecting how each trait evolves through selection imposed by interacting partners. Here, we integrate quantitative genetics, coevolutionary theory and network science to explore how trait correlations affect the coevolution of mutualistic species not only in pairs of species but also in species-rich networks across space. We show that genetic correlations may determine the pace of coevolutionary change, affect species abundances and fuel divergence among populations of the same species. However, this trait divergence promoted by genetic correlations is partially buffered by the nested structure of species-rich mutualisms. Our study, therefore, highlights how coevolution and its ecological consequences may result from conflicting processes at different levels of organisation, ranging from genes to communities.

Coevolution Creates Complex Mosaics across Large Landscapes.

The spatial distribution of populations can influence the evolutionary outcome of species interactions. The variation in direction and strength of selection across local communities creates geographic selection mosaics that, when combined with gene flow and genomic processes such as genome duplication or hybridization, can fuel ongoing coevolution. A fundamental problem to solve is how coevolution proceeds when many populations that vary in their ecological outcomes are connected across large landscapes. Here we use a lattice model to explore this problem. Our results show that the complex interrelationships among the elements of the geographic mosaic of coevolution can lead to the formation of clusters of populations with similar phenotypes that are larger than expected by local selection. Our results indicate that neither the spatial distribution of phenotypes nor the spatial differences in magnitude and direction of selection alone dictate coevolutionary dynamics: the geographic mosaic of coevolution affects formation of phenotypic clusters, which in turn affect the spatial and temporal dynamics of coevolution. Because the formation of large phenotypic clusters depends on gene flow, we predict that current habitat fragmentation will change the outcomes of geographic mosaics, coupling spatial patterns in selection and phenotypes.

Extreme diversification of floral volatiles within and among species of Lithophragma (Saxifragaceae).

A major challenge in evolutionary biology is to understand how complex traits of multiple functions have diversified and codiversified across interacting lineages and geographic ranges. We evaluate intra- and interspecific variation in floral scent, which is a complex trait of documented importance for mutualistic and antagonistic interactions between plants, pollinators, and herbivores. We performed a large-scale, phylogenetically structured study of an entire plant genus (Lithophragma, Saxifragaceae), of which several species are coevolving with specialized pollinating floral parasites of the moth genus Greya (Prodoxidae). We sampled 94 Lithophragma populations distributed across all 12 recognized Lithophragma species and subspecies, and four populations of related saxifragaceous species. Our results reveal an unusually high diversity of floral volatiles among populations, species, and clades within the genus. Moreover, we found unexpectedly major changes at each of these levels in the biosynthetic pathways used by local populations in their floral scents. Finally, we detected significant, but variable, genus- and species-level patterns of ecological convergence in the floral scent signal, including an impact of the presence and absence of two pollinating Greya moth species. We propose that one potential key to understanding floral scent variation in this hypervariable genus is its geographically diverse interactions with the obligate specialized Greya moths and, in some species and sites, more generalized copollinators.

The geographic mosaic of coevolution in mutualistic networks.

Ecological interactions shape adaptations through coevolution not only between pairs of species but also through entire multispecies assemblages. Local coevolution can then be further altered through spatial processes that have been formally partitioned in the geographic mosaic theory of coevolution. A major current challenge is to understand the spatial patterns of coadaptation that emerge across ecosystems through the interplay between gene flow and selection in networks of interacting species. Here, we combine a coevolutionary model, network theory, and empirical information on species interactions to investigate how gene flow and geographical variation in selection affect trait patterns in mutualistic networks. We show that gene flow has the surprising effect of favoring trait matching, especially among generalist species in species-rich networks typical of pollination and seed dispersal interactions. Using an analytical approximation of our model, we demonstrate that gene flow promotes trait matching by making the adaptive landscapes of different species more similar to each other. We use this result to show that the progressive loss of gene flow associated with habitat fragmentation may undermine coadaptation in mutualisms. Our results therefore provide predictions of how spatial processes shape the evolution of species-rich interactions and how the widespread fragmentation of natural landscapes may modify the coevolutionary process.

Indirect effects drive coevolution in mutualistic networks.

Ecological interactions have been acknowledged to play a key role in shaping biodiversity. Yet a major challenge for evolutionary biology is to understand the role of ecological interactions in shaping trait evolution when progressing from pairs of interacting species to multispecies interaction networks. Here we introduce an approach that integrates coevolutionary dynamics and network structure. Our results show that non-interacting species can be as important as directly interacting species in shaping coevolution within mutualistic assemblages. The contribution of indirect effects differs among types of mutualism. Indirect effects are more likely to predominate in nested, species-rich networks formed by multiple-partner mutualisms, such as pollination or seed dispersal by animals, than in small and modular networks formed by intimate mutualisms, such as those between host plants and their protective ants. Coevolutionary pathways of indirect effects favour ongoing trait evolution by promoting slow but continuous reorganization of the adaptive landscape of mutualistic partners under changing environments. Our results show that coevolution can be a major process shaping species traits throughout ecological networks. These findings expand our understanding of how evolution driven by interactions occurs through the interplay of selection pressures moving along multiple direct and indirect pathways.

Species-rich networks and eco-evolutionary synthesis at the metacommunity level.

Understanding how ecological and evolutionary processes interdependently structure biosphere dynamics is a major challenge in the era of worldwide ecosystem degradation. However, our knowledge of 'eco-evolutionary feedbacks' depends largely on findings from simple systems representing limited spatial scales and involving few species. Here we review recent conceptual developments for the understanding of multispecies coevolutionary processes and then discuss how new lines of concepts and methods will accelerate the integration of ecology and evolutionary biology. To build a research workflow for integrating insights into spatiotemporal dynamics of species-rich systems, we focus on the roles of 'metacommunity hub' species, whose population size and/or genetic dynamics potentially control landscape- or regional-scale phenomena. As large amounts of network data are becoming available with high-throughput sequencing of various host-symbiont, prey-predator, and symbiont-symbiont interactions, we suggest it is now possible to develop bases for the integrated understanding and management of species-rich ecosystems.

Diversification of Trait Combinations in Coevolving Plant and Insect Lineages.

Closely related species often have similar traits and sometimes interact with the same species. A crucial problem in evolutionary ecology is therefore to understand how coevolving species diverge when they interact with a set of closely related species from another lineage rather than with a single species. We evaluated geographic differences in the floral morphology of all woodland star plant species (Lithophragma, Saxifragaceae) that are pollinated by Greya (Prodoxidae) moths. Flowers of each woodland star species differed depending on whether plants interact locally with one, two, or no pollinating moth species. Plants of one species grown in six different environments showed few differences in floral traits, suggesting that the geographic differences are not due significantly to trait plasticity. Greya moth populations also showed significant geographic divergence in morphology, depending on the local host and on whether the moth species co-occurred locally. Divergence in the plants and the moths involved shifts in combinations of partially correlated traits, rather than any one trait. The results indicate that the geographic mosaic of coevolution can be amplified as coevolving lineages diversify into separate species and come together in different combinations in different ecosystems.

Network Structure and Selection Asymmetry Drive Coevolution in Species-Rich Antagonistic Interactions.

Ecological interactions shape and are shaped by the evolution of interacting species. Mathematical models and empirical work have explored the multiple ways coevolution could occur in small sets of species, revealing that the addition of even one species can change the coevolutionary dynamics of a pairwise interaction. As a consequence, one of the current challenges in evolutionary biology is to understand how species-rich assemblages evolve and coevolve as networks of interacting species. We combined an adaptive network framework, a trait evolutionary model, and data on network structure to study how network organization affects and is affected by selection in antagonistic interactions such as parasitism, predation, and herbivory. We explored how selection imposed by interactions shapes the evolution of attack and defense traits, parameterizing our models with structural information from 31 empirical assemblages of antagonistic species. In the simulations, the form of coevolution in antagonistic interactions is affected by the intensity and asymmetry of the selection imposed by the interacting partners. Transient escalation in attack and defensive traits was the most prevalent form of coevolutionary dynamics, especially in networks formed by modules of highly interacting species. Fluctuating evolution of traits was observed when the intensity of selection was higher in exploiters than in victims and was especially favored in nested networks. At the species level, highly connected species experienced higher temporal variation in selection regardless of the network structure, resulting in high trait mismatching with their partners. The mismatched patterns of highly connected species, in turn, may explain the emergence of modularity in antagonistic interactions in which selection is stronger on exploiters than on their victims. Our results highlight the roles of different aspects of network structure on antagonistic coevolution: nestedness shapes coevolutionary dynamics, whereas modularity emerges as one outcome of coevolutionary dynamics.

Nutrient availability affects floral scent much less than other floral and vegetative traits in Lithophragma bolanderi.

Background and Aims: Many plant-pollinator interactions are mediated by floral scents that can vary among species, among populations within species and even among individuals within populations. This variation could be innate and unaffected by the environment, but, because many floral volatiles have amino-acid precursors, scent variation also could be affected by differences in nutrient availability among environments. In plants that have coevolved with specific pollinators, natural selection is likely to favour low phenotypic plasticity in floral scent even under different conditions of nutrient availability if particular scents or scent combinations are important for attracting local pollinators. Methods: Clonal pairs of multiple seed-families of two Lithophragma bolanderi (Saxifragaceae) populations were subjected to a high and a low nutrient treatment. These plants are pollinated primarily by host-specific Greya moths. It was evaluated how nutrient treatment affected variation in floral scent relative to other vegetative and reproductive traits. Key Results: Floral scent strength (the per-flower emission rate) and composition were unaffected by nutrient treatment, but low-nutrient plants produced fewer and lighter leaves, fewer scapes and fewer flowers than high-nutrient plants. The results held in both populations, which differed greatly in the number and composition of floral scents produced. Conclusions: The results reveal a strong genetic component both to scent composition and emission level, and partly contrasts with the only previous study that has assessed the susceptibility of floral volatile signals to variation in the abundance of nutrients. These results, and the tight coevolutionary relationship between Lithophragma plants and their specialized Greya moth pollinators, indicate that reproductive traits important to coevolving interactions, such as the floral scent of L. bolanderi, may be locally specialized and more canalized than other traits important for plant fitness.

Gene flow and metacommunity arrangement affects coevolutionary dynamics at the mutualism-antagonism interface.

Interspecific interactions are affected by community context and, as a consequence, show spatial variation in magnitude and sign. The selective forces imposed by interactions at the mutualism-antagonism interface are a consequence of the traits involved and their matching between species. If mutualistic and antagonistic communities are linked by gene flow, coevolution between a pair of interacting species is influenced by how selection varies in space. Here we investigate the effects of metacommunity arrangement, i.e. patterns of connection between communities and the number of communities, on the coevolutionary dynamics between two species for which the sign and magnitude of the interaction varies across the landscape. We quantify coevolutionary outcome as an index that can be decomposed into the contribution of intraspecific genetic diversity and interspecific interaction. We show that polymorphisms and mismatches are an expected outcome, which is influenced by spatial structure, interaction strength and the degree of gene flow. The index describes how variation is distributed within and between species, and provides information on the directionality of the mismatches and polymorphisms. Finally, we argue that depending on metacommunity arrangement, some communities have disproportionate roles in maintaining genetic diversity, with implications for the coevolution of interacting species in a fragmented landscape.

Unravelling Darwin's entangled bank: architecture and robustness of mutualistic networks with multiple interaction types.

Trying to unravel Darwin's entangled bank further, we describe the architecture of a network involving multiple forms of mutualism (pollination by animals, seed dispersal by birds and plant protection by ants) and evaluate whether this multi-network shows evidence of a structure that promotes robustness. We found that species differed strongly in their contributions to the organization of the multi-interaction network, and that only a few species contributed to the structuring of these patterns. Moreover, we observed that the multi-interaction networks did not enhance community robustness compared with each of the three independent mutualistic networks when analysed across a range of simulated scenarios of species extinction. By simulating the removal of highly interacting species, we observed that, overall, these species enhance network nestedness and robustness, but decrease modularity. We discuss how the organization of interlinked mutualistic networks may be essential for the maintenance of ecological communities, and therefore the long-term ecological and evolutionary dynamics of interactive, species-rich communities. We suggest that conserving these keystone mutualists and their interactions is crucial to the persistence of species-rich mutualistic assemblages, mainly because they support other species and shape the network organization.

Coevolution, local adaptation and ecological speciation.

Coevolution is one of the major processes organizing the earth's biodiversity, but it remains unclear when and how it may generate species diversity. The study by Parchman et al. () in this issue of Molecular Ecology provides the clearest evidence to date that divergent local adaptation in a coevolving interaction may lead to speciation on one side of an interaction but not necessarily on the other side. Red crossbills in North America have diversified into ecotypes that specialize on different conifer species, use different calls and vary in the extent to which they are nomadic or sedentary. This new study evaluated genomic divergence among nine crossbill ecotypes. The authors found low overall genomic divergence among many of the ecotypes, but the sedentary South Hills crossbills, which are specialized to eat the seeds of a unique population of lodgepole pines, showed substantial divergence from other crossbills at a small number of genomic regions. These results corroborate past studies showing local coadaptation of the morphological traits of South Hills crossbills and lodgepole pines, and premating isolation of the South Hills crossbills from other populations. Together, the past and new results suggest that local coevolution with lodgepole pines has led to reduced gene flow between South Hills crossbills and other crossbills.

Below-ground plant-fungus network topology is not congruent with above-ground plant-animal network topology.

In nature, plants and their pollinating and/or seed-dispersing animals form complex interaction networks. The commonly observed pattern of links between specialists and generalists in these networks has been predicted to promote species coexistence. Plants also build highly species-rich mutualistic networks below ground with root-associated fungi, and the structure of these plant-fungus networks may also affect terrestrial community processes. By compiling high-throughput DNA sequencing data sets of the symbiosis of plants and their root-associated fungi from three localities along a latitudinal gradient, we uncovered the entire network architecture of these interactions under contrasting environmental conditions. Each network included more than 30 plant species and hundreds of mycorrhizal and endophytic fungi belonging to diverse phylogenetic groups. The results were consistent with the notion that processes shaping host-plant specialization of fungal species generate a unique linkage pattern that strongly contrasts with the pattern of above-ground plant-partner networks. Specifically, plant-fungus networks lacked a "nested" architecture, which has been considered to promote species coexistence in plant-partner networks. Rather, the below-ground networks had a conspicuous "antinested" topology. Our findings lead to the working hypothesis that terrestrial plant community dynamics are likely determined by the balance between above-ground and below-ground webs of interspecific interactions.

Assembly of complex plant-fungus networks.

Species in ecological communities build complex webs of interaction. Although revealing the architecture of these networks is fundamental to understanding ecological and evolutionary dynamics in nature, it has been difficult to characterize the structure of most species-rich ecological systems. By overcoming this limitation through next-generation sequencing technology, we herein uncover the network architecture of below-ground plant-fungus symbioses, which are ubiquitous to terrestrial ecosystems. The examined symbiotic network of a temperate forest in Japan includes 33 plant species and 387 functionally and phylogenetically diverse fungal taxa, and the overall network architecture differs fundamentally from that of other ecological networks. In contrast to results for other ecological networks and theoretical predictions for symbiotic networks, the plant-fungus network shows moderate or relatively low levels of interaction specialization and modularity and an unusual pattern of 'nested' network architecture. These results suggest that species-rich ecological networks are more architecturally diverse than previously recognized.

Floral scent contributes to interaction specificity in coevolving plants and their insect pollinators.

Chemical defenses, repellents, and attractants are important shapers of species interactions. Chemical attractants could contribute to the divergence of coevolving plant-insect interactions, if pollinators are especially responsive to signals from the local plant species. We experimentally investigated patterns of daily floral scent production in three Lithophragma species (Saxifragaceae) that are geographically isolated and tested how scent divergence affects attraction of their major pollinator-the floral parasitic moth Greya politella (Prodoxidae). These moths oviposit through the corolla while simultaneously pollinating the flower with pollen adhering to the abdomen. The complex and species-specific floral scent profiles were emitted in higher amounts during the day, when these day-flying moths are active. There was minimal divergence found in petal color, which is another potential floral attractant. Female moths responded most strongly to scent from their local host species in olfactometer bioassays, and were more likely to oviposit in, and thereby pollinate, their local host species in no-choice trials. The results suggest that floral scent is an important attractant in this interaction. Local specialization in the pollinator response to a highly specific plant chemistry, thus, has the potential to contribute importantly to patterns of interaction specificity among coevolving plants and highly specialized pollinators.