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.     

On the coevolution of pathogen and host: II. Selfing hosts and haploid pathogens

The theory of discrete time coevolution is applied to the problem of maintenance of genetic polymorphism with selfing hosts and haploid pathogens. It is shown that the usual simplifying assumption, discrete synchroized generations with no intraspecific frequency-dependent selection, precludes stability. This situation is not corrected by the incorporation of special features such as mutation, alternate hosts, partial outcrossing of the hosts, or genetic recombination in the pathogen population.     

Effects of environmental heterogeneity on victim-exploiter coevolution

We study victim-exploiter coevolution in a spatially heterogeneous island model. In each species, fitness consequences of between-species interactions are controlled by a single haploid diallelic locus. Our emphasis is on the conditions for the maintenance of genetic variation, the dynamic patterns observed, the extent of local adaptation and genetic differentiation between different demes, and on how different parameters (such as the strength and heterogeneity in selection, migration rates, and the number of sites) affect the dynamic and static behavior of the system. We show that under spatially homogeneous selection the maintenance of genetic variation is possible through asynchronous nonlinear dynamics where the allele frequencies in a majority of demes quickly synchronize but the rest do not. Spatially heterogeneous selection can maintain genetic variation even if migration rates are maximal. This happens in an oscillatory way. Genetic variation is most likely to be maintained at high levels if the heterogeneity in selection is large. If there are some restrictions on migration, genetic variation can be maintained at a stable equilibrium. This behavior is most likely at intermediate migration rates. In this case, the system can exhibit high spatial subdivision as measured by F(ST) values but relatively low local adaptation.     

Natural variation in the Pto disease resistance gene within species of wild tomato (Lycopersicon). II. Population genetics of Pto

Disease resistance to the bacterial pathogen Pseudomonas syringae pv. tomato (Pst) in the host species Lycopersicon esculentum, the cultivated tomato, and the closely related L. pimpinellifolium is triggered by the physical interaction between the protein products of the host resistance (R) gene Pto and the pathogen avirulence genes AvrPto and AvrPtoB. Sequence variation at the Pto locus was surveyed in natural populations of seven species of Lycopersicon to test hypotheses of host-parasite coevolution and functional adaptation of the Pto gene. Pto shows significantly higher nonsynonymous polymorphism than 14 other non-R-gene loci in the same samples of Lycopersicon species, while showing no difference in synonymous polymorphism, suggesting that the maintenance of amino acid polymorphism at this locus is mediated by pathogen selection. Also, a larger proportion of ancestral variation is maintained at Pto as compared to these non-R-gene loci. The frequency spectrum of amino acid polymorphisms known to negatively affect Pto function is skewed toward low frequency compared to amino acid polymorphisms that do not affect function or silent polymorphisms. Therefore, the evolution of Pto appears to be influenced by a mixture of both purifying and balancing selection.     

Balancing selection at the tomato RCR3 Guardee gene family maintains variation in strength of pathogen defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant-pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the "Guard-Hypothesis," R proteins (the "guards") can sense modification of target molecules in the host (the "guardees") by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the "guardee-effector" interface for pathogen recognition, natural selection acts on the "guard-guardee" interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen.     

Hot Spots, Cold Spots, and the Geographic Mosaic Theory of Coevolution

Species interactions commonly coevolve as complex geographic mosaics of populations shaped by differences in local selection and gene flow. We use a haploid matching-alleles model for coevolution to evaluate how a pair of species coevolves when fitness interactions are reciprocal in some locations ("hot spots") but not in others ("cold spots"). Our analyses consider mutualistic and antagonistic interspecific interactions and a variety of gene flow patterns between hot and cold spots. We found that hot and cold spots together with gene flow influence coevolutionary dynamics in four important ways. First, hot spots need not be ubiquitous to have a global influence on evolution, although rare hot spots will not have a disproportionate impact unless selection is relatively strong there. Second, asymmetries in gene flow can influence local adaptation, sometimes creating stable equilibria at which species experience minimal fitness in hot spots and maximal fitness in cold spots, or vice versa. Third, asymmetries in gene flow are no more important than asymmetries in population regulation for determining the maintenance of local polymorphisms through coevolution. Fourth, intraspecific allele frequency differences among hot and cold spot populations evolve under some, but not all, conditions. That is, selection mosaics are indeed capable of producing spatially variable coevolutionary outcomes across the landscapes over which species interact. Altogether, our analyses indicate that coevolutionary trajectories can be strongly shaped by the geographic distribution of coevolutionary hot and cold spots, and by the pattern of gene flow among populations.     

Coevolution in temporally variable environments

Many potentially mutualistic interactions are conditional, with selection that varies between mutualism and antagonism over space and time. We develop a genetic model of temporally variable coevolution that incorporates stochastic fluctuations between mutualism and antagonism. We use this model to determine conditions necessary for the coevolution of matching traits between a host and a conditional mutualist. Using an analytical approximation, we show that matching traits will coevolve when the geometric mean interaction is mutualistic. When this condition does not hold, polymorphism and trait mismatching are maintained, and coevolutionary cycles may result. Numerical simulations verify this prediction and suggest that it remains robust in the presence of temporal autocorrelation. These results are compared with those from spatial models with unrestricted movement. The comparisons demonstrate that gene flow is unnecessary for generating empirical patterns predicted by the geographic mosaic theory of coevolution.     

Compensatory cis-trans evolution and the dysregulation of gene expression in interspecific hybrids of Drosophila

Hybrids between species are often characterized by novel gene-expression patterns. A recent study on allele-specific gene expression in hybrids between species of Drosophila revealed cases in which cis- and trans-regulatory elements within species had coevolved in such a way that changes in cis-regulatory elements are compensated by changes in trans-regulatory elements. We hypothesized that such coevolution should often lead to gene misexpression in the hybrid. To test this hypothesis, we estimated allele-specific expression and overall expression levels for 31 genes in D. melanogaster, D. simulans, and their F1 hybrid. We found that 13 genes with cis-trans compensatory evolution are in fact misexpressed in the hybrid. These represent candidate genes whose dysregulation might be the consequence of coevolution of cis- and trans-regulatory elements within species. Using a mathematical model for the regulation of gene expression, we explored the conditions under which cis-trans compensatory evolution can lead to misexpression in interspecific hybrids.     

The maintenance of nucleocytoplasmic polymorphism in a metapopulation: the case of gynodioecy

In gynodioecious species, gender is generally determined by epistatic interactions between cytoplasmic and nuclear loci. However, theoretical studies suggest that, for a joint polymorphism at both cytoplasmic and nuclear loci to be maintained in a panmictic population, selection must act differently on the various genotypes that determine the same gender. Here we show that, in a metapopulation with local extinction and restricted gene flow, nucleocytoplasmic polymorphism can be maintained without these differences. We use deterministic simulations. We assume that gene flow occurred only at recolonization. Founder effects create genetic variance between populations in the metapopulation, and local population growth is faster when the local frequency of females is high. Group selection phenomena are involved in the maintenance of the joint polymorphism in the metapopulation. The frequency of females in the metapopulation at equilibrium is higher than in a panmictic population with the same genetic system. However, these conclusions hold only if nuclear alleles restoring male fertility are dominant.     

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.     

Comparing environmental and genetic variance as adaptive response to fluctuating selection

Phenotypic variation within populations has two sources: genetic variation and environmental variation. Here, we investigate the coevolution of these two components under fluctuating selection. Our analysis is based on the lottery model in which genetic polymorphism can be maintained by negative frequency-dependent selection, whereas environmental variation can be favored due to bet-hedging. In our model, phenotypes are characterized by a quantitative trait under stabilizing selection with the optimal phenotype fluctuating in time. Genotypes are characterized by their phenotypic offspring distribution, which is assumed to be Gaussian with heritable variation for its mean and variance. Polymorphism in the mean corresponds to genetic variance while the width of the offspring distribution corresponds to environmental variance. We show that increased environmental variance is favored whenever fluctuations in the selective optima are sufficiently strong. Given the environmental variance has evolved to its optimum, genetic polymorphism can still emerge if the distribution of selective optima is sufficiently asymmetric or leptokurtic. Polymorphism evolves in a diagonal direction in trait space: one type becomes a canalized specialist for the more common ecological conditions and the other type a de-canalized bet-hedger thriving on the less-common conditions. All results are based on analytical approximations, complemented by individual-based simulations.     

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.     

Genetic covariation of the marine fungal symbiont Haloguignardia irritans (Ascomycota, Pezizomycotina) with its algal hosts Cystoseira and Halidrys (Phaeophyceae, Fucales) along the west coast of North America

The fungal endophyte Haloguignardia irritans induces gall formation on the brown algal genera Cystoseira and Halidrys occurring from Oregon to Baja California, Mexico. Here we examine genetic covariation and compare rDNA phylogenies to investigate the coevolutionary histories of H. irritans and its algal hosts. Despite recognition of H. irritans as a single morphological species, internal transcribed spacer rDNA sequences representative of its geographic range are characterized by sequence variation at the intraspecific to intrageneric levels. An assessment of parallel cladogenesis between endophyte and host phylogenies provides evidence for a combination of independent fungal divergence and host jumping, similar to that observed in terrestrial lichens. Our results suggest that reduced gene flow due to geographic isolation is a major contributing factor to more concerted covariation observed at one island site, rather than to differences among algal host species alone. Because geography and its effects on gene flow can create heterogeneous mosaics of coevolution for symbioses in terrestrial environments, our results support the notion that conservation efforts toward the maintenance of genetic diversity in marine environments should likewise consider geographic complexity and its effects on coevolving marine species.     

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.     

Host–parasite genotypic interactions in the honey bee: the dynamics of diversity

Parasites are thought to be a major driving force shaping genetic variation in their host, and are suggested to be a significant reason for the maintenance of sexual reproduction. A leading hypothesis for the occurrence of multiple mating (polyandry) in social insects is that the genetic diversity generated within‐colonies through this behavior promotes disease resistance. This benefit is likely to be particularly significant when colonies are exposed to multiple species and strains of parasites, but host–parasite genotypic interactions in social insects are little known. We investigated this using honey bees, which are naturally polyandrous and consequently produce genetically diverse colonies containing multiple genotypes (patrilines), and which are also known to host multiple strains of various parasite species. We found that host genotypes differed significantly in their resistance to different strains of the obligate fungal parasite that causes chalkbrood disease, while genotypic variation in resistance to the facultative fungal parasite that causes stonebrood disease was less pronounced. Our results show that genetic variation in disease resistance depends in part on the parasite genotype, as well as species, with the latter most likely relating to differences in parasite life history and host–parasite coevolution. Our results suggest that the selection pressure from genetically diverse parasites might be an important driving force in the evolution of polyandry, a mechanism that generates significant genetic diversity in social insects.     

Population structure affected by excess gene flow in self-pollinating Elymus nutans and E. burchan-buddae (Triticeae: Poaceae)

Seed-mediated gene flow can considerably affect population genetic structure of strictly self-pollinating species, but little is known on the extent and nature of such gene flow among pastoral plant populations. Molecular fingerprints provide a powerful tool to address the relevant issues. Genetic structure of 22 populations of two self-pollinating pasture species, Elymus nutans and E. burchan-buddae, collected from various altitudes of the Qinghai–Tibetan Plateau was studied using fluorescence-based amplified fragment length polymorphism technique. Analysis of molecular variance revealed 42.97% and 37.63% among-population variation for the two Elymus species, respectively, indicating that the majority of the total variation presented within populations. This result contradicts the common genetic variation pattern for a selfing plant species: lower genetic variation within populations. Further analysis suggested higher level of gene flow among populations within the same region than among different regions across the sampled area for the two Elymus species. STRUCTURE analyses of the Elymus populations indicated an evident admixture genetic structure, particularly among neighboring populations from the same region, supporting the hypothesis of considerable seed dispersal among populations. The excess within-region gene flow of E. nutans and E. burchan-buddae might be caused by grazing animals that promote seed dispersal when moved around the pastoral lands during foraging. The among-population gene flow promulgated by grazing animals may promote the maintenance of genetic diversity in the pasture species, particularly in small and fragmented populations within a given region.     

Asymmetry in host and parasitoid diffuse coevolution: when the red queen has to keep a finger in more than one pie

BACKGROUND: Coevolution between pairs of antagonistic species is generally considered an endless "arms race" between attack and defense traits to counteract the adaptive responses of the other species. PRESENTATION OF THE HYPOTHESIS: When more than two species are involved, diffuse coevolution of hosts and parasitoids could be asymmetric because consumers can choose their prey whereas preys do not choose their predator. This asymmetry may lead to differences in the rate of evolution of the antagonistic species in response to selection. The more long-standing the coevolution of a given pair of antagonistic populations, the higher should be the fitness advantage for the consumer. Therefore, the main prediction of the hypothesis is that the consumer trophic level is more likely to win the coevolution race. TESTING THE HYPOTHESIS: We propose testing the asymmetry hypothesis by focusing on the tritrophic system plant/aphid/aphid parasitoid. The analysis of the genetic variability in the virulence of several parasitoid populations and in the defenses of several aphid species or several clones of the same aphid species could be compared. Moreover, the analysis of the neutral population genetic structure of the parasitoid as a function of the aphid host, the plant host and geographic isolation may complement the detection of differences between host and parasitoid trophic specialization. IMPLICATIONS OF THE HYPOTHESIS: Genetic structures induced by the arms race between antagonistic species may be disturbed by asymmetry in coevolution, producing neither rare genotype advantages nor coevolutionary hotspots. Thus this hypothesis profoundly changes our understanding of coevolution and may have important implications in terms of pest management.     

Antagonistic coevolution mediated by phenotypic differences between quantitative traits

Many well-studied coevolutionary interactions between predators and prey or hosts and parasites are mediated by quantitative traits. In some interactions, such as those between cuckoos and their hosts, interactions are mediated by the degree of phenotype matching among species, and a significant body of theory has been developed to predict the coevolutionary dynamics and outcomes of such interactions. In a large number of other cases, however, interactions are mediated by the extent to which the phenotype of one species exceeds that of the other. For these cases-which are arguably more numerous-few theoretical predictions exist for coevolutionary dynamics and outcomes. Here we develop and analyze mathematical models of interspecific interactions mediated by the extent to which the quantitative trait of one species exceeds that of the other. Our results identify important differences from previously studied models based on trait matching. First, our results show that cyclical dynamics are possible only if the strength of coevolutionary selection exceeds a threshold and stabilizing selection acts on the interacting traits. Second, our results demonstrate that significant levels of genetic polymorphism can be maintained only when cyclical dynamics occur. This result leads to the unexpected prediction that maintenance of genetic polymorphism is enhanced by strong selection. Finally, our results demonstrate that there is no a priori reason to expect the traits of interacting species should match in any literal sense, even in the absence of gene flow among populations.     

Genetic Variation in a Heterogeneous Environment. I. Temporal Heterogeneity and the Absolute Dominance Model

The conditions for a stable polymorphism and the equilibrium gene frequency in an infinite population are compared when there is spatial or temporal environmental heterogeneity for the absolute dominance model. For temporal variation the conditions for stability are more restrictive and the equilibrium gene frequency is often at a low gene frequency. In a finite population, temporal environmental heterogeneity for the absolute dominance model was found to be quite ineffective in maintaining genetic variation and is often less effective than no selection at all. For comparison, the maximum maintenance for temporal variation is related to the overdominant model. In general, cyclic environmental variation was found to be more effective at maintaining genetic variation than where the environment varies stochastically. The importance of temporal environmental variation and the maintenance of genetic variation is discussed.