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.     

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.     

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.     

Brucellosis, botflies, and brainworms: the impact of edge habitats on pathogen transmission and species extinction

Ecological interactions between species that prefer different habitat types but come into contact in edge regions at the interfaces between habitat types are modeled via reaction-diffusion systems. The primary sort of interaction described by the models is competition mediated by pathogen transmission. The models are somewhat novel because the spatial domains for the variables describing the population densities of the interacting species overlap but do not coincide. Conditions implying coexistence of the two species or the extinction of one species are derived. The conditions involve the principal eigenvalues of elliptic operators arising from linearizations of the model system around equilibria with only one species present. The conditions for persistence or extinction are made explicit in terms of the parameters of the system and the geometry of the underlying spatial domains via estimates of the principal eigenvalues. The implications of the models with respect to conservation and refuge design are discussed. Received: 10 June 1999 / Revised version: 7 July 2000 / Published online: 20 December 2000     

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.     

Comparison of quantitative genetic parameters between two natural populations of a selfing plant species, Medicago truncatula Gaertn.

 In this paper we compare mean values, heritability estimates, coefficient of genetic variation, and genetic correlations among several fitness components of two natural populations of a selfing plant species, Medicago truncatula L. It is shown that the population that had been found most polymorphic for molecular markers in a previous study was also the most variable for quantitative characters. Depending on the traits, the larger heritabilities in this population were due to either larger coefficients of genetic variances or smaller coefficients of environmental variances. Whereas genetic and phenotypic correlation matrices were very similar within each population, they were quite different between populations. In particular, although a positive correlation between age and size at maturity was found in both populations, the correlation between age at maturity and reproductive success was negative in the more variable population (late flowering plant, with a larger size at flowering, produced fewer pods), whereas no correlation was observed in the less variable population. We suggest that while in the less variable population all individuals have a high reproductive effort, several strategies coexist in the more variable population, with some early-flowering genotypes showing a high reproductive effort and other late-flowering genotypes showing a larger competitive ability through increased vegetative growth. Received: 25 June 1996 / Accepted: 11 October 1996     

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within‐group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within‐population selection‐induced reductions in variation. An important implication of this ‘sweep‐before‐differentiation’ model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.     

Geographic and population variation in pollinating seed-consuming interactions between senita cacti (Lophocereus schottii) and senita moths (Upiga virescens)

Interspecific interactions can vary within and among populations and geographic locations. This variation can subsequently influence the evolution and coevolution of species interactions. We investigated population and geographic variation in traits important to pollinating seed-consuming interactions between the senita cactus (Lophocereus schottii) and its obligate pollinating moth (Upiga virescens), both of which are geographically restricted to the Sonoran Desert. Female moths actively pollinate senita flowers and oviposit onto flowers. Their larvae consume developing seeds and fruit of flowers pollinated by females. Traits important to this interaction include fruit set from moth pollination, fruit survivorship, and costs of fruit consumption by larvae. We studied these traits for five populations at two widely separated geographic locations. On average, 37% of flowers set fruit, 22% of flowers produced mature fruit, and larvae consumed 25% of immature fruit pollinated by female senita moths. Senita cactus and senita moth interactions were strongly mutualistic in all populations that we studied. Although one population had statistically lower fruit set and fruit production than the other four, all five populations were qualitatively similar in fruit production, costs, and patterns of fruit survivorship. Hand-pollination experiments suggested that fruit set was resource-limited in all but this one population. Apparent pollen limitation in the one population explains the quantitative differences in fruit set and fruit survivorship among the populations. As predicted by theory and exemplified by the senita mutualism, specialized and/or obligate interactions vary little among populations and geographic locations. Received: 5 March 1999 / Accepted: 25 June 1999     

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.     

An explicit representation of the Luria–Delbrück distribution

The probability distribution of the number of mutant cells in a growing single-cell population is presented in explicit form. We use a discrete model for mutation and population growth which in the limit of large cell numbers and small mutation rates reduces to certain classical models of the Luria-Delbrück distribution. Our results hold for arbitrarily large values of the mutation rate and for cell populations of arbitrary size. We discuss the influence of cell death on fluctuation experiments and investigate a version of our model that accounts for the possibility that both daughter cells of a non-mutant cell might be mutants. An algorithm is presented for the quick calculation of the distribution. Then, we focus on the derivation of two essentially different limit laws, the first of which applies if the population size tends to infinity while the mutation rate tends to zero such that the product of mutation rate times population size converges. The second limit law emerges after a suitable rescaling of the distribution of non-mutant cells in the population and applies if the product of mutation rate times population size tends to infinity. We discuss the distribution of mutation events for arbitrary values of the mutation rate and cell populations of arbitrary size, and, finally, consider limit laws for this distribution with respect to the behavior of the product of mutation rate times population size. Thus, the present paper substantially extends results due to Lea and Coulson (1949), Bartlett (1955), Stewart et al. (1990), and others.     

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.     

Coevolution of dispersal in a parasitoid–host system

Interspecific interactions and the evolution of dispersal are both of interest when considering the potential impact of habitat fragmentation on community ecology, but the interaction between these processes is not well studied. We address this by considering the coevolution of dispersal strategies in a host–parasitoid system. An individual-based host–parasitoid metapopulation model was constructed for a patchy environment, allowing for evolution in dispersal rates of both species. Highly rarefied environments with few suitable patches selected against dispersal in both species, as did relatively static environments. Provided that parasitoids persist, all the variables studied led to stable equilibria in dispersal rates for both species. There was a tendency toward higher dispersal rates in parasitoids because of the asymmetric relationships of the two species to the patches: vacant patches are most valuable for hosts, but unsuitable for parasitoids, which require an established host population to reproduce. High host dispersal rate was favoured by high host population growth rate, and in the parasitoid by high growth rates in both species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Interspecific competition among macroparasites in a density-dependent host population

 We analyze the dynamics of a community of macroparasite species that share the same host. Our work extends an earlier framework for a host species that would grow exponentially in the absence of parasitism, to one where an uninfected host population is regulated by factors other than parasites. The model consists of one differential equation for each parasite species and a single density-dependent nonlinear equation for the host. We assume that each parasite species has a negative binomial distribution within the host and there is zero covariance between the species (exploitation competition). New threshold conditions on model parameters for the coexistence and competitive exclusion of parasite species are derived via invadibility and stability analysis of corresponding equilibria. The main finding is that the community of parasite species coexisting at the stable equilibrium is obtained by ranking the species according t! o th e minimum host density H ^* above which a parasite species can grow when rare: the lower H ^* , the higher the competitive ability. We also show that ranking according to the basic reproduction number Q ₀ does not in general coincide with ranking according to H ^* . The second result is that the type of interaction between host and parasites is crucial in determining the competitive success of a parasite species, because frequency-dependent transmission of free-living stages enhances the invading ability of a parasite species while density-dependent transmission makes a parasite very sensitive to other competing species. Finally, we show that density dependence in the host population entails a simplification of the portrait of possible outcomes with respect to previous studies, because all the cases resulting in the exponential growth of host and parasite populations are eliminated.. Received: 24 June 1996 / Revised version: 28 April 1998     

Density dependent selection incorporating intraspecific competition 1. A haploid model

A haploid model is introduced and analyzed in which intraspecific competition is incorporated within a density dependent framework. It is assumed that each genotype has a unique carrying capacity corresponding to the equilibrium population size when fixed for that type. Each genotypic fitness at a single multi-allelic locus is a function of a distinctive effective population size formed by adding the numbers of each genotype present, weighted by an intraspecific competition coefficient. As a result, the fitnesses depend upon the relative frequencies of the various genotypes as well as the total population size. Intergenotypic interactions can have a profound effect upon the outcome of the population. In particular, when the density effect of one individual upon another depends upon their respective genotypes, a unique stable interior equilibrium is possible in which all alleles are present. This stands in contrast to the purely density dependent haploid system in which the only possible stable state corresponds to fixation for the type with the highest carrying capacity. In the present model selective advantage is determined by a balance between carrying capacity and sensitivity to density pressures from other genotypes. Fixation for the genotype with the highest carrying capacity, for instance, will not be stable if it exerts a sufficiently weak competitive effect upon the other genotypes. In the diallelic case, maintenance of both alleles at a stable equilibrium requires that the net intragenotypic competition between individuals of like genotype be stronger than that between unlike types. As for purely density regulated systems, there may be no stable equilibria and/or regular and chaotic cycling may occur. The results may also be interpreted in terms of a discrete time model of interspecific competition with each haplotype representing a different species.     

Temperature effects on gametophyte life-history traits and geographic distribution of two cryptic kelp species

A major determinant of the geographic distribution of a species is expected to be its physiological response to changing abiotic variables over its range. The range of a species often corresponds to the geographic extent of temperature regimes the organism can physiologically tolerate. Many species have very distinct life history stages that may exhibit different responses to environmental factors. In this study we emphasized the critical role of the haploid microscopic stage (gametophyte) of the life cycle to explain the difference of edge distribution of two related kelp species. Lessonia nigrescens was recently identified as two cryptic species occurring in parapatry along the Chilean coast: one located north and the other south of a biogeographic boundary at latitude 29-30°S. Six life history traits from microscopic stages were identified and estimated under five treatments of temperature in eight locations distributed along the Chilean coast in order to (1) estimate the role of temperature in the present distribution of the two cryptic L. nigrescens species, (2) compare marginal populations to central populations of the two cryptic species. In addition, we created a periodic matrix model to estimate the population growth rate (λ) at the five temperature treatments. Differential tolerance to temperature was demonstrated between the two species, with the gametophytes of the Northern species being more tolerant to higher temperatures than gametophytes from the south. Second, the two species exhibited different life history strategies with a shorter haploid phase in the Northern species contrasted with considerable vegetative growth in the Southern species haploid stage. These results provide strong ecological evidence for the differentiation process of the two cryptic species and show local adaptation of the life cycle at the range limits of the distribution. Ecological and evolutionary implications of these findings are discussed.     

Arabidopsis thaliana as a model for the study of plant–virus co-evolution

Understanding plant–virus coevolution requires wild systems in which there is no human manipulation of either host or virus. To develop such a system, we analysed virus infection in six wild populations of Arabidopsis thaliana in Central Spain. The incidence of five virus species with different life-styles was monitored during four years, and this was analysed in relation to the demography of the host populations. Total virus incidence reached 70 per cent, which suggests a role of virus infection in the population structure and dynamics of the host, under the assumption of a host fitness cost caused by the infection. Maximum incidence occurred at early growth stages, and co-infection with different viruses was frequent, two factors often resulting in increased virulence. Experimental infections under controlled conditions with two isolates of the most prevalent viruses, cauliflower mosaic virus and cucumber mosaic virus, showed that there is genetic variation for virus accumulation, although this depended on the interaction between host and virus genotypes. Comparison of Q_ST-based genetic differentiations between both host populations with F_ST genetic differentiation based on putatively neutral markers suggests different selection dynamics for resistance against different virus species or genotypes. Together, these results are compatible with a hypothesis of plant–virus coevolution.     

Intersimple sequence repeat (ISSR) variation in Lactoris fernandeziana (Lactoridaceae), a rare endemic of the Juan Fernández Archipelago, Chile

Sixteen populations and 89 individuals of Lactoris fernandeziana were examined for variation in intersimple sequence repeat (ISSR) banding patterns. The species is a rare endemic of Masatierra Island in the Juan Fernández Archipelago, and is the only member of the endemic family Lactoridaceae. Five populations showed a single genotype whereas the other 11 populations had from two to 16 multilocus genotypes. Over 73% of the ISSR diversity occurred across populations, with only about 27% within populations. Diversity among populations results from the presence of different subsets of loci within each population rather than unique loci within populations; only two populations displayed novel loci, with one and three in each. Levels of differentiation at ISSR loci among populations are not correlated with geographic distance on Masatierra; rather, the pattern of variation is mosaic. The presence of differentiated local populations is concordant with the geitonogamous breeding system of the species and suggests low levels of long distance pollen or seed dispersal. The mosaic pattern of ISSR variation on Masatierra may result, in part, from drift and inbreeding in small populations following fragmentation of a once more continuous distribution of Lactoris with the formation of canyons by erosion. Also, the generation of new ISSR loci by mutation could occur with rare, sporadic gene flow among populations accounting for the mosaic pattern of variation and the paucity of unique alleles within populations. The ISSR results for Lactoris suggest that studies of morphological, ecological and physiological features may elucidate differentiation among populations of L. fernandeziana . Field studies have demonstrated that plants occur both in the dense forest understory and in the full sunlight in forest openings.     

PLOIDY AND EVOLUTION BY SEXUAL SELECTION: A COMPARISON OF HAPLOID AND DIPLOID FEMALE CHOICE MODELS NEAR FIXATION EQUILIBRIA

We compare the stability properties of haploid and diploid models of Fisherian sexual selection (with male contribution limited to sperm) by examining both models at equilibria for which a male trait is fixed or absent. Haploid and diploid two locus diallelic models share the property that the stability of such fixation equilibria is determined by the relationship between the harmonic mean of relative preference values for the common male trait, weighted by the frequency of the preferences, and the relative viability associated with the common male trait. When diploid females with heterozygotic-based preferences express preference strengths intermediate between homozygote-based preferences, then boundary equilibria of haploid and diploid models share many stability properties. However, even with intermediate heterozygote preferences, haploid and diploid models do differ: (1) for a particular frequency of the preference allele, both fixation boundaries can be stable for the diploid model, and (2) with over- or underdominance at the preference locus (a possibility precluded in the haploid model), a fixation boundary in the diploid model may show two switches in its stability state for increasing frequencies of one of the preference alleles. These differences are due not just to the impossibility of dominance in haploid models, but also to the larger number of diploid genotypes.     

Variation in reproductive isolation across a species range

Reproductive isolation is often variable within species, a phenomenon that while largely ignored by speciation studies, can be leveraged to gain insight into the potential mechanisms driving the evolution of genetic incompatibilities. We used experimental greenhouse crosses to characterize patterns of reproductive isolation among three divergent genetic lineages of Campanulastrum americanum that occur in close geographic proximity in the Appalachian Mountains. Substantial, asymmetrical reproductive isolation for survival due to cytonuclear incompatibility was found among the lineages (up to 94% reduction). Moderate reductions in pollen viability, as well as cytoplasmic male sterility, were also found between some Mountain populations. We then compared these results to previously established patterns of reproductive isolation between these Mountain lineages and a fourth, widespread Western lineage to fully characterize reproductive isolation across the complete geographic and genetic range of C. americanum. Reproductive isolation for survival and pollen viability was consistent across studies, indicating the evolution of the underlying genetic incompatibilities is primarily determined by intrinsic factors. In contrast, reproductive isolation for germination was only found when crossing Mountain populations with the Western lineage, suggesting the underlying genetic incompatibility is likely influenced by environmental or demographic differences between the two lineages. Cytoplasmic male sterility was also limited in occurrence, being restricted to a handful of Mountain populations in a narrow geographic range. These findings illustrate the complexity of speciation by demonstrating multiple, independent genetic incompatibilities that lead to a mosaic of genetic divergence and reproductive isolation across a species range.