Dos and don'ts of testing the geographic mosaic theory of coevolution

The geographic mosaic theory of coevolution is stimulating much new research on interspecific interactions. We provide a guide to the fundamental components of the theory, its processes and main predictions. Our primary objectives are to clarify misconceptions regarding the geographic mosaic theory of coevolution and to describe how empiricists can test the theory rigorously. In particular, we explain why confirming the three main predicted empirical patterns (spatial variation in traits mediating interactions among species, trait mismatching among interacting species and few species-level coevolved traits) does not provide unequivocal support for the theory. We suggest that strong empirical tests of the geographic mosaic theory of coevolution should focus on its underlying processes: coevolutionary hot and cold spots, selection mosaics and trait remixing. We describe these processes and discuss potential ways each can be tested.     

The evolutionary response of predators to dangerous prey: hotspots and coldspots in the geographic mosaic of coevolution between garter snakes and newts

The "geographic mosaic" approach to understanding coevolution is predicated on the existence of variable selection across the landscape of an interaction between species. A range of ecological factors, from differences in resource availability to differences in community composition, can generate such a mosaic of selection among populations, and thereby differences in the strength of coevolution. The result is a mixture of hotspots, where reciprocal selection is strong, and coldspots, where reciprocal selection is weak or absent, throughout the ranges of species. Population subdivision further provides the opportunity for nonadaptive forces, including gene flow, drift, and metapopulation dynamics, to influence the coevolutionary interaction between species. Some predicted results of this geographic mosaic of coevolution include maladapted or mismatched phenotypes, maintenance of high levels of polymorphism, and prevention of stable equilibrium trait combinations. To evaluate the potential for the geographic mosaic to influence predator-prey coevolution, we investigated the geographic pattern of genetically determined TTX resistance in the garter snake Thamnophis sirtalis over much of the range of its ecological interaction with toxic newts of genus Taricha. We assayed TTX resistance in over 2900 garter snakes representing 333 families from 40 populations throughout western North America. Our results provide dramatic evidence that geographic structure is an important component in coevolutionary interactions between predators and prey. Resistance levels vary substantially (over three orders of magnitude) among populations and over short distances. The spatial array of variation is consistent with two areas of intense evolutionary response by predators ("hotspots") surrounded by clines of decreasing resistance. Some general predictions of the geographic mosaic process are supported, including clinal variation in phenotypes, polymorphism in some populations, and divergent outcomes of the interaction between predator and prey. Conversely, our data provide little support for one of the major predictions, mismatched values of interacting traits. Two lines of evidence suggest selection is paramount in determining population variation in resistance. First, phylogenetic information indicates that two hotspots of TTX resistance have evolved independently. Second, in the one region that TTX levels in prey have been quantified, resistance and toxicity levels match almost perfectly over a wide phenotypic and geographic range. However, these results do not preclude the role the nonadaptive forces in generating the overall geographic mosaic of TTX resistance. Much work remains to fill in the geographic pattern of variation among prey populations and, just as importantly, to explore the variation in the ecology of the interaction that occurs within populations.     

Genetic dimension of the coevolution of virulence-resistance in Drosophila -- parasitoid wasp relationships

Variations observed in parasite virulence and host resistance may be the outcome of coevolutionary processes. Recent theoretical developments have led to a 'geographic mosaic theory' of coevolution according to which there are some localities where reciprocal selection occurs (hot spots) and others where it is strongly reduced (cold spots). Studies of host-parasitoid systems back this up, revealing a geographical variation of traits subjected to antagonistic selection governed by variations in the strength of the ecological interactions. A more detailed analysis of the genetic basis of these geographic variations in a model system -- the interaction between Drosophila melanogaster and its specific parasitoid Leptopilina boulardi -- suggests that cold spots and hot spots are also driven by the amount of genetic variation available for the trait considered. Our approach, based on isolating reference strains, has been found to predict the result of sympatric interactions and it will be helpful in identifying the selective forces responsible for the coevolution. In this model, host resistance to a standardised reference strain is a weak predictor of the outcome of interactions in the field, and the main parameter accounting for the geographic variations is the number of host species available, with less parasitoid virulence towards D. melanogaster being found in areas displaying a more diversified host community.     

Increasing productivity accelerates host-parasite coevolution

Host-parasite coevolution is believed to influence a range of evolutionary and ecological processes, including population dynamics, evolution of diversity, sexual reproduction and parasite virulence. The impact of coevolution on these processes will depend on its rate, which is likely to be affected by the energy flowing through an ecosystem, or productivity. We addressed how productivity affected rates of coevolution during a coevolutionary arms race between experimental populations of bacteria and their parasitic viruses (phages). As hypothesized, the rate of coevolution between bacterial resistance and phage infectivity increased with increased productivity. This relationship can in part be explained by reduced competitiveness of resistant bacteria in low compared with high productivity environments, leading to weaker selection for resistance in the former. The data further suggest that variation in productivity can generate variation in selection for resistance across landscapes, a result that is crucial to the geographic mosaic theory of coevolution.     

Coevolution and maladaptation

Many of the most commonly cited examples of exquisite adaptationare of coevolved symbioses. As we learn more about the coevolutionaryprocess, however, it is becoming increasingly evident that coevolutionmay also keep populations moderately maladapted much of thetime. As a result, coevolving populations may only rarely occupyadaptive peaks, because the selective landscape is under continualchange through reciprocal selection on the species themselves.These shifting patterns of coadaptation are further shaped bythe geographic structure of most species. Selection mosaicsacross landscapes and coevolutionary hotspots can favor differentevolutionary trajectories in different populations. The combinedaction of gene flow, random genetic drift, and local extinctionof populations may then continually remold these local patterns,creating a geographic mosaic in the degrees of maladaptationfound within local interactions. Recent mathematical modelsof the geographic mosaic of coevolution suggest that complexmosaics of maladaptation are a likely consequence of spatiallystructured species interactions. These models indicate thatthe spatial structure of maladaptation may depend upon the typeof coevolutionary interaction, the underlying selection mosaic,and patterns of gene flow across landscapes. By maintaininglocal polymorphisms and driving the divergence of populations,coevolution may produce spatial patterns of maladaptation thatare a source of ongoing innovation and diversification in speciesinteractions.     

The maintenance of genetic diversity under host–parasite coevolution in finite,structured populations

As a corollary to the Red Queen hypothesis, host–parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host–parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.     

Natural selection drives the fine-scale divergence of a coevolutionary arms race involving a long-mouthed weevil and its obligate host plant

Background

One of the major recent advances in evolutionary biology is the recognition that evolutionary interactions between species are substantially differentiated among geographic populations. To date, several authors have revealed natural selection pressures mediating the geographically-divergent processes of coevolution. How local, then, is the geographic structuring of natural selection in coevolutionary systems?

Results

I examined the spatial scale of a "geographic selection mosaic," focusing on a system involving a seed-predatory insect, the camellia weevil (Curculio camelliae), and its host plant, the Japanese camellia (Camellia japonica). In this system, female weevils excavate camellia fruits with their extremely-long mouthparts to lay eggs into seeds, while camellia seeds are protected by thick pericarps. Quantitative evaluation of natural selection demonstrated that thicker camellia pericarps are significantly favored in some, but not all, populations within a small island (Yakushima Island, Japan; diameter ca. 30 km). At the extreme, camellia populations separated by only several kilometers were subject to different selection pressures. Interestingly, in a population with the thickest pericarps, camellia individuals with intermediate pericarp thickness had relatively high fitness when the potential costs of producing thick pericarps were considered. Also importantly, some parameters of the weevil - camellia interaction such as the severity of seed infestation showed clines along temperature, suggesting the effects of climate on the fine-scale geographic differentiation of the coevolutionary processes.

Conclusion

These results show that natural selection can drive the geographic differentiation of interspecific interactions at surprisingly small spatial scales. Future studies should reveal the evolutionary/ecological outcomes of the "fine scale geographic mosaics" in biological communities.     

Coevolutionary interactions between a haploid species and a diploid species

We investigate a general model describing coevolutionary interaction between a haploid population and a diploid population, each with two alleles at a single locus. Both species are allowed to evolve, with the fitness of the genotypes of each species assumed to depend linearly on the frequencies of the genotypes of the other species. We explore the resulting outcomes of these interactions, in particular determining the location of equilibria under various conditions. The coevolution here is much more complex than that between two haploid populations and allows for the possibility of two polymorphic equilibria. To allow for further analysis, we construct a semi-symmetric model. The variety of outcomes possible even in this second model provides support for the geographic mosaic theory of coevolution by suggesting the possibility of small local populations coevolving to very different outcomes, leading to a shifting geographic mosaic as neighboring populations interact with each other through migration.     

Coevolutionary clines across selection mosaics

Abstract. Much of the dynamics of coevolution may be driven by the interplay between geographic variation in reciprocal selection (selection mosaics) and the homogenizing action of gene flow. We develop a genetic model of geographically structured coevolution in which gene flow links coevolving communities that may differ in both the direction and magnitude of reciprocal selection. The results show that geographically structured coevolution may lead to allele-frequency clines within both interacting species when fitnesses are spatially uniform or spatially heterogeneous. Furthermore, the results show that the behavior and shape of clines differ dramatically among different types of coevolutionary interaction. Antagonistic interactions produce dynamic clines that change shape rapidly through time, producing shifting patterns of local adaptation and maladaptation. Unlike antagonistic interactions, mutualisms generate stable equilibrium patterns that lead to fixed spatial patterns of adaptation. Interactions that vary between mutualism and antagonism produce both equilibrium and dynamic clines. Furthermore, the results demonstrate that these interactions may allow mutualisms to persist throughout the geographic range of an interaction, despite pockets of locally antagonistic selection. In all cases, the coevolved spatial patterns of allele frequencies are sensitive to the relative contributions of gene flow, selection, and overall habitat size, indicating that the appropriate scale for studies of geographically structured coevolution depends on the relative contributions of each of these factors.     

Gene flow and geographically structured coevolution

Many of the dynamic properties of coevolution may occur at the level of interacting populations, with local adaptation acting as a force of diversification, as migration between populations homogenizes these isolated interactions. This interplay between local adaptation and migration may be particularly important in structuring interactions that vary from mutualism to antagonism across the range of an interacting set of species, such as those between some plants and their insect herbivores, mammals and trypanosome parasites, and bacteria and plasmids that confer antibiotic resistance. Here we present a simple geographically structured genetic model of a coevolutionary interaction that varies between mutualism and antagonism among communities linked by migration. Inclusion of geographic structure with gene flow alters the outcomes of local interactions and allows the maintenance of allelic polymorphism across all communities under a range of selection intensities and rates of migration. Furthermore, inclusion of geographic structure with gene flow allows fixed mutualisms to be evolutionarily stable within both communities, even when selection on the interaction is antagonistic within one community. Moreover, the model demonstrates that the inclusion of geographic structure with gene flow may lead to considerable local maladaptation and trait mismatching as predicted by the geographic mosaic theory of coevolution.     

Genetic diversity,population structure and sex‐biased dispersal in three co‐evolving species

Genetic diversity and spatial structure of populations are important for antagonistic coevolution. We investigated genetic variation and population structure of three closely related European ant species: the social parasite Harpagoxenus sublaevis and its two host species Leptothorax acervorum and Leptothorax muscorum. We sampled populations in 12 countries and analysed eight microsatellite loci and an mtDNA sequence. We found high levels of genetic variation in all three species, only slightly less variation in the host L. muscorum. Using a newly introduced measure of differentiation (Jost’s D_est ), we detected strong population structuring in all species and less male‐biased dispersal than previously thought. We found no phylogeographic patterns that could give information on post‐glacial colonization routes – northern populations are as variable as more southern populations. We conclude that conditions for Thompson’s geographic mosaic of coevolution are ideal in this system: all three species show ample genetic variation and strong population structure.     

Parasite local adaptation in a geographic mosaic

A central prediction of the geographic mosaic theory of coevolution is that coevolving interspecific interactions will show varying degrees of local maladaptation. According to the theory, much of this local maladaptation is driven by selection mosaics and spatially intermingled coevolutionary hot and cold spots, rather than a simple balance between gene flow and selection. Here I develop a genetic model of host-parasite coevolution that is sufficiently general to incorporate selection mosaics, coevolutionary hot and cold spots, and a diverse array of genetic systems of infection/resistance. Results from this model show that the selection mosaics experienced by the interacting species are an important determinant of the sign and magnitude of local maladaptation. In some cases, this effect may be stronger than a previously described effect of relative rates of parasite and host gene flow. These results provide the first theoretical evidence that selection mosaics and coevolutionary hot and cold spots per se determine the magnitude and sign of local maladaptation. At the same time, however, these results demonstrate that coevolution in a geographic mosaic can lead to virtually any pattern of local adaptation or local maladaptation. Consequently, empirical studies that describe only patterns of local adaptation or maladaptation do not provide evidence either for or against the theory.     

Emergent Dynamics of Fairness in the Spatial Coevolution of Proposer and Responder Species in the Ultimatum Game

While spatially local interactions are ubiquitous between coevolving species sharing recourses (e.g., plant-insect interactions), their effects on such coevolution processes of strategies involving the share of a resource are still not clearly understood. We construct a two-dimensional spatial model of the coevolution of the proposer and responder species in the ultimatum game (UG), in which a pair of proposer and responder individuals at each site plays the UG. We investigate the effects of the locality of interactions and the intensity of selection on the emergence of fairness between these species. We show that the lower intensity of selection favors fair strategies in general, and there are no significant differences in the evolution of fairness between the cases with local and global interactions when the intensity of selection is low. However, as the intensity of selection becomes higher, the spatially local interactions contribute to the evolution of fairer strategies more than the global interactions, even though fair strategies become more difficult to evolve. This positive effect of spatial interactions is expected to be due to the mutual benefit of fairness for both proposer and responder species in future generations, which brings about a dynamic evolution process of fairness.     

Toward an ecological synthesis: a case for habitat selection

Habitat selection, and its associated density and frequency-dependent evolution, has a profound influence on such vital phenomena as population regulation, species interactions, the assembly of ecological communities, and the origin and maintenance of biodiversity. Different strategies of habitat selection, and their importance in ecology and evolution, can often be revealed simply by plots of density in adjacent habitats. For individual species, the strategies are closely intertwined with mechanisms of population regulation, and with the persistence of populations through time. For interacting species, strategies of habitat selection are not only responsible for species coexistence, but provide one of the most convenient mechanisms for measuring competition, and the various community structures caused by competitive interactions. Other kinds of interactions, such as those between predators and prey, demonstrate that an understanding of the coevolution of habitat-selection strategies among strongly interacting species is essential to properly interpret their spatial and temporal dynamics. At the evolutionary scale, the frequency dependence associated with habitat selection may often allow populations to diverge and diversify into separate species. Habitat selection thereby demonstrates how we can map microevolutionary strategies in behavior onto their population and community consequences, and from there, onto macroevolutionary patterns of speciation and adaptive radiation. We can anticipate that future studies of habitat selection will not only help us complete those maps, but that they will also continue to enrich the panoply of ideas that shape evolutionary ecology.     

Spatial geographic mosaic in an aquatic predator-prey network

The geographic mosaic theory of coevolution predicts 1) spatial variation in predatory structures as well as prey defensive traits, and 2) trait matching in some areas and trait mismatching in others mediated by gene flow. We examined gene flow and documented spatial variation in crushing resistance in the freshwater snails Mexipyrgus churinceanus, Mexithauma quadripaludium, Nymphophilus minckleyi, and its relationship to the relative frequency of the crushing morphotype in the trophically polymorphic fish Herichthys minckleyi. Crushing resistance and the frequency of the crushing morphotype did show spatial variation among 11 naturally replicated communities in the Cuatro Ciénegas valley in Mexico where these species are all endemic. The variation in crushing resistance among populations was not explained by geographic proximity or by genetic similarity in any species. We detected clear phylogeographic patterns and limited gene flow for the snails but not for the fish. Gene flow among snail populations in Cuatro Ciénegas could explain the mosaic of local divergence in shell strength and be preventing the fixation of the crushing morphotype in Herichthys minckleyi. Finally, consistent with trait matching across the mosaic, the frequency of the fish morphotype was negatively correlated with shell crushing resistance likely reflecting the relative disadvantage of the crushing morphotype in communities where the snails exhibit relatively high crushing resistance.     

Geographical ecology and variation of plant-seed disperser interactions: southern Spanish junipers and frugivorous thrushes

Spatial and temporal predictability in the mutual selective pressures of plants and frugivorous birds is a prerequisite for coevolution to occur. I examine the interaction patterns of strongly frugivorous thrushes (Turdus spp.) and their major winter food plants (Juniperus spp., Cupressaceae) and how they vary in space and time. Spatial congruency, rarely considered in seed dispersal studies, is studied at three spatial scales: 1) the total species range; 2) regional distribution; and 3) local abundance and its variation between seasons. Southern Spanish frugivorous thrushes and junipers show very low congruence in distribution patterns at each of these scales. Most juniper species show geographic distributions that are nested within the geographic ranges of thrush species. Bird species showed greater habitat breadth values than plants and were found in a greater percentage of localities. The local bird abundance was strongly correlated across years and sites with the local availability of juniper cones. Cone production varied markedly between years, but the rankings for different species in different years were statistically concordant at mid-elevation and lowland sites. Both bird abundance and cone production showed greater temporal than spatial variability. Variation of cone productions at both temporal and spatial scales was greater than variability in bird abundance. Species with strong interactions of mutual dependence showed very low values of biogeographic congruence, caused by differences in geographic range and habitat specificity. This obviously limits the possibilities for pairwise, specific coevolution to occur. However, mutual effects of species groups are possible to the extent that the component species are ecologically interchangeable in their selective effects and other constraints on coevolution are not operating. The approach used here to examine the patterns of species interactions at different biogeographic scales might prove useful in comparative studies of plant-animal interactions.     

生态基因组学研究进展

生态基因组学是一个整合生态学、分子遗传学和进化基因组学的新兴交叉学科。生态基因组学将基因组学的研究手段和方法引入生态学领域,通过将群体基因组学、转录组学、蛋白质组学等手段与方法将个体、种群及群落、生态系统不同层次的生态学相互作用整合起来,确定在生态学响应及相互作用中具有重要意义的关键的基因和遗传途径,阐明这些基因及遗传途径变异的程度及其生态和进化后果的特征,从基因水平探索有机体响应天然环境(包括生物与非生物的环境因子)的遗传学机制。生态基因组学的研究对象可以分为模式生物与非模式生物两大类。拟南芥、酿酒酵母等模式生物在生态基因组学领域发挥了重要作用。随着越来越多基因组学技术的开发与完善,越来越多的非模式生物生态基因组学的研究将为生态学的发展提供重要的理论与实践依据。生态基因组学最核心的方法包括寻找序列变异、研究基因差异表达和分析基因功能等方法。生态基因组学已广泛渗透到生态学的相关领域中,将会在生物对环境的响应、物种间的相互作用、进化生态学、全球变化生态学、入侵生态学、群落生态学等研究领域发挥更大的作用。     

INTEGRATING LANDSCAPE GENOMICS AND SPATIALLY EXPLICIT APPROACHES TO DETECT LOCI UNDER SELECTION IN CLINAL POPULATIONS

Uncovering the genetic basis of adaptation hinges on the ability to detect loci under selection. However, population genomics outlier approaches to detect selected loci may be inappropriate for clinal populations or those with unclear population structure because they require that individuals be clustered into populations. An alternate approach, landscape genomics, uses individual‐based approaches to detect loci under selection and reveal potential environmental drivers of selection. We tested four landscape genomics methods on a simulated clinal population to determine their effectiveness at identifying a locus under varying selection strengths along an environmental gradient. We found all methods produced very low type I error rates across all selection strengths, but elevated type II error rates under “weak” selection. We then applied these methods to an AFLP genome scan of an alpine plant, Campanula barbata, and identified five highly supported candidate loci associated with precipitation variables. These loci also showed spatial autocorrelation and cline patterns indicative of selection along a precipitation gradient. Our results suggest that landscape genomics in combination with other spatial analyses provides a powerful approach for identifying loci potentially under selection and explaining spatially complex interactions between species and their environment.     

Spatial and temporal dynamics in a sexual selection mosaic

Selective regimes and phenotypic optima could either change smoothly and in a clinal fashion or be spatially organized in a more unpredictable mosaic pattern over the geographic landscape. When natural or sexual selection is driven by intra- or interspecific biotic interactions, fine-grained spatial variation in selective regimes could result in selection mosaics rather than clinal variation in selection. We investigated temporal variation and spatial organization in sexual selection on male body size along an ecological coastal-inland gradient of a polymorphic damselfly Ischnura elegans. Body size increased in a clinal fashion along this gradient: animals were smaller in size at the coast, but became larger in the inland areas. In contrast, the sexual selection regimes on male body size showed evidence of more fine-grained spatial organization with no evidence for a clinal pattern and low spatial autocorrelations between populations. These spatially fine-grained sexual selection regimes varied in sign and magnitude and were driven by a combination of the densities of heritable female color morphs and local female body sizes. We suggest that the spatial organization of the selective regimes can be interpreted as a sexual selection mosaic that is influenced by highly localized density- and frequency-dependent social interactions.     

Coevolution across landscapes: a spatially explicit model of parasitoid-host coevolution

The geographic mosaic theory of coevolution suggests that population spatial structure may have a strong impact on coevolutionary dynamics. Therefore, coevolution must be studied across geographic scales, not just in single populations. To examine the impact of movement rate on coevolutionary dynamics, we developed a spatially explicit model of host–parasitoid coevolution. We described space as a coupled-map lattice and assumed that resistance (defined as the ability of a host to encapsulate a parasitoid egg) and virulence (defined as the successful parasitization of a host) traits were graded and costly. The model explicitly detailed population and evolutionary dynamics. When holding all parameters constant and varying only the movement rate of the host and parasitoid, profoundly different dynamics were observed. We found that fluctuations in the mean levels of resistance and virulence in the global population were greatest when the movement rate of the host and parasitoid was high. In addition, we found that the variation in resistance and virulence levels among neighboring patches was greatest when the movement rates of the host and parasitoid was low. However, as the distance among patches increased, so did the variation in resistance and virulence levels regardless of movement rate. These generalizations did not hold when spatial patterns in the distribution of resistance and virulence traits, such as spirals, were observed. Finally, we found that the evolution of resistance and virulence caused the abundance of hosts to increase and the abundance of parasitoids to decrease. As a result, the spatial distribution of hosts and parasitoids was influenced.