Phenotypic mismatches reveal escape from arms-race coevolution

Because coevolution takes place across a broad scale of time and space, it is virtually impossible to understand its dynamics and trajectories by studying a single pair of interacting populations at one time. Comparing populations across a range of an interaction, especially for long-lived species, can provide insight into these features of coevolution by sampling across a diverse set of conditions and histories. We used measures of prey traits (tetrodotoxin toxicity in newts) and predator traits (tetrodotoxin resistance of snakes) to assess the degree of phenotypic mismatch across the range of their coevolutionary interaction. Geographic patterns of phenotypic exaggeration were similar in prey and predators, with most phenotypically elevated localities occurring along the central Oregon coast and central California. Contrary to expectations, however, these areas of elevated traits did not coincide with the most intense coevolutionary selection. Measures of functional trait mismatch revealed that over one-third of sampled localities were so mismatched that reciprocal selection could not occur given current trait distributions. Estimates of current locality-specific interaction selection gradients confirmed this interpretation. In every case of mismatch, predators were “ahead” of prey in the arms race; the converse escape of prey was never observed. The emergent pattern suggests a dynamic in which interacting species experience reciprocal selection that drives arms-race escalation of both prey and predator phenotypes at a subset of localities across the interaction. This coadaptation proceeds until the evolution of extreme phenotypes by predators, through genes of large effect, allows snakes to, at least temporarily, escape the arms race.     

The coevolutionary dynamics of antagonistic interactions mediated by quantitative traits with evolving variances

Quantitative traits frequently mediate coevolutionary interactions between predator and prey or parasite and host. Previous efforts to understand and predict the coevolutionary dynamics of these interactions have generally assumed that standing genetic variation is fixed or absent altogether. We develop a genetically explicit model of coevolution that bridges the gap between these approaches by allowing genetic variation itself to evolve. Analysis of this model shows that the evolution of genetic variance has important consequences for the dynamics and outcome of coevolution. Of particular importance is our demonstration that coevolutionary cycles can emerge in the absence of stabilizing selection, an outcome not possible in previous models of coevolution mediated by quantitative traits. Whether coevolutionary cycles evolve depends upon the strength of selection, the number of loci, and the rate of mutation in each of the interacting species. Our results also generate novel predictions for the expected sign and magnitude of linkage disequilibria in each species.     

Red Queens and ESS: the coevolution of evolutionary rates

Summary The Red Queen principle states that a set of interacting species reaches an evolutionary equilibrium at which all their rates of coevolution exactly balance each other. The lag-load model, which is one way of searching for Red Queens, has, by itself, previously predicted that they do not exist. But this model has assumed that infinite maladaptedness is possible. The lag-load model is improved by assuming that once the lag load of all but one species is determined, so is that of the final species. This assumption eliminates the possibility of infinite maladaptedness. Its result is to allow the lag-load model to yield Red Queen coevolution. It does this whether or not speciation and extinction rates are included. Thus the lag-load model is harmonized with the earlier Red Queen model derived from studies of predation.Because of the intercorrelation of phenotypic traits, the predatory model concluded that the eventual stable rate of coevolution must be zero (except for intermittent bursts after some correlation or compromise is successfully broken). Another model that predicts stable coevolutionary rates of zero is that of evolutionarily stable strategies (ESS).Red Queen assumes that the more extreme a phenotypic trait is, the better it is, and that there are no constraints on the growth of such a phenotypic trait value. Such traits are the key to the Red Queen prediction of progressive coevolution. ESS models make no such assumptions. Eliminating unbounded traits from the model of predator-victim evolution changed its prediction from progressive coevolution to stasis. Before this paper, no model had dealt simultaneously with both unbounded and constrained traits.To handle both sorts of phenotypic traits at the same time in the same model, we abandoned lag load as a measure of evolutionary rate (lag loads do not uniquely determine phenotype). Instead, we used the traditional assumption that rate is proportional to the slope of the adaptive landscape. A model, relying on continuous evolutionary game theory, was developed and simulated under various conditions in two or three species sets, with up to five independent traits coevolving simultaneously. The results were: (1) there was always a set of equilibrium densities eventually achieved by coevolution; if the population interaction represented by this stable coevolutionary state is also stable, then the system should persist whether it evolves further or not; (2) whenever traits were present which were unbounded and best at their most extreme values, then a Red Queen emerged; (3) whenever traits were present which were correlated with each other or constrained below infinity, then an ESS emerged; (4) if both types were present, both results occurred: Red Queen in the unbounded traits and ESS in the constrained ones.Because unbounded traits may not exist, the Red Queen may have no domain. But the domain of ESS is real. ESS should lead to the evolutionary pattern called punctuated equilibrium. The changes in design rules which punctuate stasis should lead to an ever-expanding independence of traits from each other, i.e. to more and more refined differentiation. A single set of design rules which governs a set of species is called a fitness-generating function. Such functions may help to define the concepts of adaptive zone and ecological guild.     

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.     

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 ecology and evolution of host-plant resistance to insects

Genetic techniques have yielded new insights into plant-herbivore coevolution. Quantitative genetic tests of herbivory theory reveal that in some cases insect herbivores impose selection on resistance traits. Also, some resistance traits are costly while others appear not to be, and genetic models can explain these results. Genetic variation in plant resistance influences insect community structure by modifying interactions of herbivores with competitors and natural enemies. Therefore, models of multispecies coevolution are more realistic than pairwise coevolutionary models. Ecological genetics will facilitate further theoretical and empirical exploration of multispecies coevolution of plants and herbivores.     

Sensory exploitation and sexual conflict

Much of the literature on male-female coevolution concerns the processes by which male traits and female preferences for these can coevolve and be maintained by selection. There has been less explicit focus on the origin of male traits and female preferences. Here, I argue that it is important to distinguish origin from subsequent coevolution and that insights into the origin can help us appreciate the relative roles of various coevolutionary processes for the evolution of diversity in sexual dimorphism. I delineate four distinct scenarios for the origin of male traits and female preferences that build on past contributions, two of which are based on pre-existing variation in quality indicators among males and two on exploitation of pre-existing sensory biases among females. Recent empirical research, and theoretical models, suggest that origin by sensory exploitation has been widespread. I argue that this points to a key, but perhaps transient, role for sexually antagonistic coevolution (SAC) in the subsequent evolutionary elaboration of sexual traits, because (i) sensory exploitation is often likely to be initially costly for individuals of the exploited sex and (ii) the subsequent evolution of resistance to sensory exploitation should often be associated with costs due to selective constraints. A review of a few case studies is used to illustrate these points. Empirical data directly relevant to the costs of being sensory exploited and the costs of evolving resistance is largely lacking, and I stress that such data would help determining the general importance of sexual conflict and SAC for the evolution of sexual dimorphism.     

The evolution of optimal female mating rate changes the coevolutionary dynamics of female resistance and male persistence

Mating decisions usually involve conflict of interests between sexes. Accordingly, males benefit from increased number of matings, whereas costs of mating favour a lower mating rate for females. The resulting sexual conflict underlies the coevolution of male traits that affect male mating success ('persistence') and female traits that affect female mating patterns ('resistance'). Theoretical studies on the coevolutionary dynamics of male persistence and female resistance assumed that costs of mating and, consequently, the optimal female mating rate are evolutionarily constant. Costs of mating, however, are often caused by male 'persistence' traits that determine mating success. Here, we present a model where the magnitude of costs of mating depend on, and evolve with, male persistence. We find that allowing costs of mating to depend on male persistence results in qualitatively different coevolutionary dynamics. Specifically, we find that male traits such as penis spikes that harm females are not predicted to exhibit runaway selection with female resistance, in contrast to previous theory that predicts indefinite escalation. We argue that it is essential to determine when and to what extent costs of mating are caused by male persistence in order to understand and accurately predict coevolutionary dynamics of traits involved in mating decisions.     

Long‐term coevolution between avian brood parasites and their hosts

Coevolutionary theory predicts that the most common long‐term outcome of the relationships between brood parasites and their hosts should be coevolutionary cycles based on a dynamic change selecting the currently least‐defended host species, given that when well‐defended hosts are abandoned, hosts will be selected to decrease their defences as these are usually assumed to be costly. This is assumed to be the case also in brood parasite‐host systems. Here I examine the frequency of the three potential long‐term outcomes of brood parasite–host coevolution (coevolutionary cycles, lack of rejection, and successful resistance) in 182 host species. The results of simple exploratory comparisons show that coevolutionary cycles are very scarce while the lack of rejection and successful resistance, which are considered evolutionary enigmas, are much more frequent. I discuss these results considering (i) the importance of different host defences at all stages of the breeding cycle, (ii) the role of phenotypic plasticity in long‐term coevolution, and (iii) the evolutionary history of host selection. I suggest that in purely antagonistic coevolutionary interactions, such as those involving brood parasites and their hosts, that although cycles will exist during an intermediate phase of the interactions, the arms race will end with the extinction of the host or with the host acquiring successful resistance. As evolutionary time passes, this resistance will force brood parasites to use previously less suitable host species. Furthermore, I present a model that represents the long‐term trajectories and outcomes of coevolutionary interactions between brood parasites and their hosts with respect to the evolution of egg‐rejection defence. This model suggests that as an increasing number of species acquire successful resistance, other unparasitized host species become more profitable and their parasitism rate and the costs imposed by brood parasitism at the population level will increase, selecting for the evolution of host defences. This means that although acceptance is adaptive when the parasitism rate and the costs of parasitism are very low, this cannot be considered to represent an evolutionary equilibrium, as conventional theory has done to date, because it is not stable.     

Tolerance traits and the stability of mutualism

Identifying factors which allow the evolution and persistence of cooperative interactions between species is a fundamental issue in evolutionary ecology. Various hypotheses have been suggested which generally focus on mechanisms that allow cooperative genotypes in different species to maintain interactions over space and time. Here, we emphasise the fact that even within mutualisms (interactions with net positive fitness effects for both partners), there may still be inherent costs, such as the occasional predation by ants upon aphids. Individuals engaged in mutualisms benefit from minimising these costs as long as it is not at the expense of breaking the interspecific interaction, which offers a net positive benefit. The most common and obvious defence traits to minimise interspecific interaction costs are resistance traits, which act to reduce encounter rate between two organisms. Tolerance traits, in contrast, minimise fitness costs to the actor, but without reducing encounter rate. Given that, by definition, it is beneficial to remain in mutualistic interactions, the only viable traits to minimise costs are tolerance-based 'defence' strategies. Thus, we propose that tolerance traits are an important factor promoting stability in mutualisms. Furthermore, because resistance traits tend to propagate coevolutionary arms races between antagonists, whilst tolerance traits do not, we also suggest that tolerance-based defence strategies may be important in facilitating the transition from antagonistic interactions into mutualisms. For example, the mutualism between ants and aphids has been suggested to have evolved from parasitism. We describe how phenotypic plasticity in honeydew production may be a tolerance trait that has prevented escalation into an antagonistic arms race and instead led to mutualistic coevolution.     

The influence of a competitor on the geographic mosaic of coevolution between crossbills and lodgepole pine

The geographic mosaic theory of coevolution posits that the form of selection between interacting species varies across a landscape with coevolution important and active in some locations (i.e., coevolutionary hotspots) but not in others (i.e., coevolutionary coldspots). We tested the hypothesis that the presence of red squirrels (Tamiasciurus hudsonicus) affects the occurrence of coevolution between red crossbills (Loxia curvirostra complex) and Rocky Mountain lodgepole pine (Pinus contorta ssp. latifolia) and thereby provides a mechanism giving rise to a geographic mosaic of selection. Red squirrels are the predominant predispersal seed predator and selective agent on lodgepole pine cones. However, in four isolated mountain ranges east and west of the Rocky Mountains, red squirrels are absent and red crossbills are the main predispersal seed predator. These isolated populations of pine have apparently evolved without Tamiasciurus for about 10,000 to 12,000 years. Based on published morphological, genetic, and paleobotanical studies, we infer that cone traits in these isolated populations that show parallel differences from cones in the Rocky Mountains have changed in parallel. We used data on crossbill and conifer cone morphology and feeding preferences and efficiency to detect whether red crossbills and lodgepole pine exhibit reciprocal adaptations, which would imply coevolution. Cone traits that act to deter Tamiasciurus and result in high ratios of cone mass to seed mass were less developed in the isolated populations. Cone traits that act to deter crossbills include larger and thicker scales and perhaps increased overlap between successive scales and were enhanced in the isolated populations. In the larger, isolated mountain ranges crossbills have evolved deeper, shorter, and therefore more decurved bills to exploit these cones. This provides crossbills with higher feeding rates, and the change in bill shape has improved efficiency by reducing the concomitant increases in body mass and daily energy expenditures that would have resulted if only bill size had increased. These parallel adaptations and counter adaptations in red crossbills and lodgepole pine are interpreted as reciprocal adaptations and imply that these crossbills and pine are in coevolutionary arms races where red squirrels are absent (i.e., coevolutionary hotspots) but not where red squirrels are present (i.e., coevolutionary cold-spots).     

Genetic correlations and the coevolutionary dynamics of three-species systems

The majority of species interact with at least several others. We develop simple genetic models of coevolution between three species where interactions are mediated by quantitative traits. We assume that one of the species has two quantitative traits, each of which governs its interaction with one of the other two species. We use this model to explore how genetic correlations between the two traits in the multivariate species shape the evolutionary dynamics and outcomes of three species interactions. Our results suggest that genetic correlations are most important when at least one of the interactions is between a predator and prey or parasite and host. In these cases, genetic correlations between traits lead to a wide variety of novel coevolutionary outcomes and dynamics. In particular, genetic correlations can affect the existence and stability of coevolutionary equilibrium points, and they can lead to recurrent or permanent maladaptation. When the three species interact only as competitors or mutualists, however, genetic correlations have no effect on the outcome of coevolution. In all cases, our results reveal the surprising conclusion that both positive and negative genetic correlations between traits have qualitatively identical effects on coevolutionary dynamics.     

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.     

The measurement of coevolution in the wild

Coevolution has long been thought to drive the exaggeration of traits, promote major evolutionary transitions such as the evolution of sexual reproduction and influence epidemiological dynamics. Despite coevolution’s long suspected importance, we have yet to develop a quantitative understanding of its strength and prevalence because we lack generally applicable statistical methods that yield numerical estimates for coevolution’s strength and significance in the wild. Here, we develop a novel method that derives maximum likelihood estimates for the strength of direct pairwise coevolution by coupling a well‐established coevolutionary model to spatially structured phenotypic data. Applying our method to two well‐studied interactions reveals evidence for coevolution in both systems. Broad application of this approach has the potential to further resolve long‐standing evolutionary debates such as the role species interactions play in the evolution of sexual reproduction and the organisation of ecological communities.     

DISENTANGLING THE CONTRIBUTION OF SEXUAL SELECTION AND ECOLOGY TO THE EVOLUTION OF SIZE DIMORPHISM IN PINNIPEDS

The positive relationship between sexual size dimorphism (SSD) and harem size across pinnipeds is often cited as a textbook example of sexual selection. It assumes that female aggregation selected for large male size via male–male competition. Yet, it is also conceivable that SSD evolved prior to polygyny due to ecological forces. We analyzed 11 life‐history traits in 35 pinniped species to determine their coevolutionary dynamics and infer their most likely evolutionary trajectories contrasting these two hypotheses. We find support for SSD having evolved prior to changes in the mating system, either as a consequence of niche partitioning during aquatic foraging or in combination with sexual selection on males to enforce copulations on females. Only subsequently did polygyny evolve, leading to further coevolution as the strength of sexual selection intensified. Evolutionary sequence analyses suggest a polar origin of pinnipeds and indicate that SSD and polygyny are intrinsically linked to a suite of ecological and life‐history traits. Overall, this study calls for the inclusion of ecological variables when studying sexual selection and argues for caution when assuming causality between coevolving traits. It provides novel insights into the role of sexual selection for the coevolutionary dynamics of SSD and mating system.     

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.     

Population mixing promotes arms race host–parasite coevolution

The consequences of host–parasite coevolution are highly contingent on the qualitative coevolutionary dynamics: whether selection fluctuates (fluctuating selection dynamic; FSD), or is directional towards increasing infectivity/resistance (arms race dynamic; ARD). Both genetics and ecology can play an important role in determining whether coevolution follows FSD or ARD, but the ecological conditions under which FSD shifts to ARD, and vice versa, are not well understood. The degree of population mixing is thought to increase host exposure to parasites, hence selecting for greater resistance and infectivity ranges, and we hypothesize this promotes ARD. We tested this by coevolving bacteria and viruses in soil microcosms and found that population mixing shifted bacteria–virus coevolution from FSD to ARD. A simple theoretical model produced qualitatively similar results, showing that mechanisms that increase host exposure to parasites tend to push dynamics towards ARD. The shift from FSD to ARD with increased population mixing may help to explain variation in coevolutionary dynamics between different host–parasite systems, and more specifically the observed discrepancies between laboratory and field bacteria–virus coevolutionary studies.     

Detecting sexual conflict and sexually antagonistic coevolution

We begin by providing an operational definition of sexual conflict that applies to both inter- and intralocus conflict. Using this definition, we examine a series of simple coevolutionary models to elucidate fruitful approaches for detecting interlocus sexual conflict and resultant sexually antagonistic coevolution. We then use published empirical examples to illustrate the utility of these approaches. Three relevant attributes emerge. First, the dynamics of sexually antagonistic coevolution may obscure the conflict itself. Second, competing models of inter-sexual coevolution may yield similar population patterns near equilibria. Third, a variety of evolutionary forces underlying competing models may be acting simultaneously near equilibria. One main conclusion is that studies of emergent patterns in extant populations (e.g. studies of population and/or female fitness) are unlikely to allow us to distinguish among competing coevolutionary models. Instead, we need more research aimed at identifying the forces of selection acting on shared traits and sexually antagonistic traits. More specifically, we need a greater number of functional studies of female traits as well as studies of the consequences of both male and female traits for female fitness. A mix of selection and manipulative studies on these is likely the most promising route.     

Inferring coevolution in a plant–pollinator network

Mutualistic interactions are at the core of community dynamics, determining dispersal, colonization and differential survival and reproduction among individuals and species. Mutualistic interactions therefore affect the fitness of interaction partners, hence modifying their respective evolutionary trajectories reciprocally, potentially leading to coevolution. Although mathematical models predict coevolution in mutualistic interaction networks, no empirical data are available. By taking into account the patterns of interactions and reconstructing evolutionary change in plant and pollinator traits, we tested the hypothesis that coevolution occurs between plants and insects that interact more frequently, or more symmetrically. To test this hypothesis, we built an interaction network with data from five flowering seasons, measured plant and insect morphology, mapped morphology on the plant and insect phylogenies, and reconstructed ancestral character changes based on maximum parsimony. We calculated an index, called the coevolutionary ratio, which represents the amount of correlated change in traits that mediate the interaction between plants and pollinators (i.e. proboscis versus corolla length, and body width and corolla aperture). Our results suggest that high frequency of interaction, i.e. the number of times two species interact, does not lead to coevolution. Instead, symmetry of interaction strength, i.e. the reciprocal similarity in the mutual effect of interaction partners, may lead to coevolution, in spite of a pervasive lack of reciprocal specialization and high interaction frequency. Although the statistical signal is quite weak, our results hold for three statistical tests of very different nature. The most specialized species, expected to be under directional selection, do not show more evolutionary change than do generalist species, expected to be under different, perhaps opposing, selective pressures. By dissecting the complexity of an interaction network we show that coevolution may partially shape functional morphology of interaction partners, thus providing the closest evidence to date of mutualistic adaptation of organisms within a community.     

The potential of a population genomics approach to analyse geographic mosaics of plant--insect coevolution

A central issue in the evolutionary ecology of species interactions is coevolution, which involves the reciprocal selection between individuals of interacting species. Understanding the importance of coevolution in shaping species interactions requires the consideration of spatial variation in their strength. This is exactly what the, recently developed, geographic mosaic theory of coevolution addresses. Another major development in the study of population ecology is the introduction of the population genomics approach in this field of research. This approach addresses spatial processes through molecular methods. It is of particular interest that population genomics is especially applicable to natural populations of non-model species. We describe how population genomics can be used in the context of the geographic mosaic of coevolution, specifically to identify coevolutionary hot-spots, and to attribute genetic variation found at specific loci to processes of selection versus trait remixing. The proposed integration of the population genomics approach with the conceptual framework of the geographic mosaic of coevolution is illustrated with a few selected, particularly demonstrative, examples from the realm of insect--plant interactions.