Metapopulation structure of a seed-predator weevil and its host plant in arms race coevolution

Fine-Scale Genetic Structure of Curculio chinensis (Coleoptera: Curculionidae) Based on Mitochondrial COI: The Role of Host Specificity and Spatial Distance

The Camellia weevil, Curculio chinensis (Chevrolat, 1978), is a dominant oligophagous pest that bores into the fruit of oil-tea Camellia. Genetic differentiation among populations in various hosts can easily occur, which hinders research on pest management. In this study, the genetic structure, genetic diversity, and phylogenetic structure of local C. chinensis populations were examined using 147 individuals (from 6 localities in Jiangxi), based on 2 mitochondrial COI markers. Results indicated that the C. chinensis population in Jiangxi exhibits a high haplotype diversity, especially for the populations from Cam. meiocarpa plantations. Structural differentiation was observed between Haplogroup 1 (73 individuals from Ganzhou, Jian, and Pingxiang) in the monoculture plantations of Cam. meiocarpa and Haplogroup 2 (75 individuals from Pingxiang and Jiujiang) in Cam. oleifera. Two haplogroups have recently undergone a demographic expansion, and Haplogroup 1 has shown a higher number of effective migrants than Haplogroup 2. This suggests that C. chinensis has been spreading from Cam. meiocarpa plantations to other oil-tea Camellia, such as Cam. oleifera. The increased cultivation of oil-tea Camellia in Jiangxi has contributed to a unique genetic structure within the C. chinensis population. This has, in turn, expanded the distribution of C. chinensis and increased migration between populations.

Conspicuous coloration of toxin-resistant predators implicates additional trophic interactions in a predator-prey arms race

Antagonistic coevolution between natural enemies can produce highly exaggerated traits, such as prey toxins and predator resistance. This reciprocal process of adaptation and counter-adaptation may also open doors to other evolutionary novelties not directly involved in the phenotypic interface of coevolution. We tested the hypothesis that predator-prey coevolution coincided with the evolution of conspicuous coloration on resistant predators that retain prey toxins. In western North America, common garter snakes (Thamnophis sirtalis) have evolved extreme resistance to tetrodotoxin (TTX) in the coevolutionary arms race with their deadly prey, Pacific newts (Taricha spp.). TTX-resistant snakes can retain large amounts of ingested TTX, which could serve as a deterrent against the snakes' own predators if TTX toxicity and resistance are coupled with a conspicuous warning signal. We evaluated whether arms race escalation covaries with bright red coloration in snake populations across the geographic mosaic of coevolution. Snake colour variation departs from the neutral expectations of population genetic structure and covaries with escalating clines of newt TTX and snake resistance at two coevolutionary hotspots. In the Pacific Northwest, bright red coloration fits an expected pattern of an aposematic warning to avian predators: TTX-resistant snakes that consume highly toxic newts also have relatively large, reddish-orange dorsal blotches. Snake coloration also seems to have evolved with the arms race in California, but overall patterns are less intuitively consistent with aposematism. These results suggest that interactions with additional trophic levels can generate novel traits as a cascading consequence of arms race coevolution across the geographic mosaic.

Plant‐associate interactions and diversification across trophic levels

Interactions between species are widely understood to have promoted the diversification of life on Earth, but how interactions spur the formation of new species remains unclear. Interacting species often become locally adapted to each other, but they may also be subject to shared dispersal limitations and environmental conditions. Moreover, theory predicts that different kinds of interactions have different effects on diversification. To better understand how species interactions promote diversification, we compiled population genetic studies of host plants and intimately associated herbivores, parasites, and mutualists. We used Bayesian multiple regressions and the BEDASSLE modeling framework to test whether host and associate population structures were correlated over and above the potentially confounding effects of geography and shared environmental variation. We found that associates' population structure often paralleled their hosts' population structure, and that this effect is robust to accounting for geographic distance and climate. Associate genetic structure was significantly explained by plant genetic structure somewhat more often in antagonistic interactions than in mutualistic ones. This aligns with a key prediction of coevolutionary theory that antagonistic interactions promote diversity through local adaptation of antagonists to hosts, while mutualistic interactions more often promote diversity via the effect of hosts' geographic distribution on mutualists' dispersal.

Effects of Fruit Sizes of Two Camellia Trees on the Larval Sizes of Curculio styracis (Roelofs, 1875): Testing the Endoparasitoid Body Size Hypothesis

Simple Summary

In endoparasitoids that feed within small discrete resource patches, such as seeds or fruits, body size could be subject to a trade-off: larger size could lead to increased overall fitness but could simultaneously increase the risk of resource depletion and starvation, resulting in a body size just below the host holding capacity. We analyzed the relationship of the larval size of the within-fruits-developing curculionid beetle Curculio styracis (Roelofs, 1875) and the size of the fruits of its two congeneric host species of Camellia to test this hypothesis. A logistic model can most accurately describe larval size in association with host-fruit size after a series of models were tested. Based on the characteristics of the optimal model, the hypothesis seemed to be confirmed because larvae that developed in host plant with larger fruits had a larger size, and larval size in both host species remained only a little below the host-fruit capacity. The novelty of the study is that this hypothesis is being tested in a more formal way using appropriate mathematical models.

Abstract

The endoparasitoid body size hypothesis suggests that the size of larvae that develop in a single host should be subject to a trade-off: larger size could lead to increase overall fitness but could simultaneously increase the risk of resource depletion and starvation, resulting in a body size just below the host holding capacity. However, this hypothesis has not been rigorously tested using mathematical models thus far. The camellia weevil, C.styracis (Coleoptera: Curculionidae), is a notorious pest attacking fruits of Camellia oleifera Abel. and C. meiocarpa Hu., in which the larvae develop within a single fruit and larval development is limited by the available food resources. We developed a feasible method to test this hypothesis. First, five models were used to describe the relationship between larval mass and host size. Then, the minimum fruit threshold that had to be met for ad libitum larval development and the corresponding larval size (W_a) of this threshold were calculated based on the characteristics of the optimal model. Finally, the difference between the measured larval size and the predicted larval size (W_a) was determined. The results showed that (1) the data were better described by a logistic function than any other equation; (2) larval size in both host plants increased with increasing fruit size until leveling off when the fruits were large enough to allow unconstrained larval development; (3) larval size remained just below the host-fruit holding capacity, as there was no difference between the measured and predicted larval sizes (W_a); and (4) larvae developed in host plant with larger fruits had a larger size. These results confirmed the endoparasitoid body size hypothesis.

The geographic mosaic of arms race coevolution is closely matched to prey population structure

Reciprocal adaptation is the hallmark of arms race coevolution. Local coadaptation between natural enemies should generate a geographic mosaic pattern where both species have roughly matched abilities across their shared range. However, mosaic variation in ecologically relevant traits can also arise from processes unrelated to reciprocal selection, such as population structure or local environmental conditions. We tested whether these alternative processes can account for trait variation in the geographic mosaic of arms race coevolution between resistant garter snakes (Thamnophis sirtalis) and toxic newts (Taricha granulosa). We found that predator resistance and prey toxin levels are functionally matched in co-occurring populations, suggesting that mosaic variation in the armaments of both species results from the local pressures of reciprocal selection. By the same token, phenotypic and genetic variation in snake resistance deviates from neutral expectations of population genetic differentiation, showing a clear signature of adaptation to local toxin levels in newts. Contrastingly, newt toxin levels are best predicted by genetic differentiation among newt populations, and to a lesser extent, by the local environment and snake resistance. Exaggerated armaments suggest that coevolution occurs in certain hotspots, but prey population structure seems to be of particular influence on local phenotypic variation in both species throughout the geographic mosaic. Our results imply that processes other than reciprocal selection, like historical biogeography and environmental pressures, represent an important source of variation in the geographic mosaic of coevolution. Such a pattern supports the role of "trait remixing" in the geographic mosaic theory, the process by which non-adaptive forces dictate spatial variation in the interactions among species.

Approximate Bayesian estimation of coevolutionary arms races

Exaggerated traits involved in species interactions have long captivated the imagination of evolutionary biologists and inspired the durable metaphor of the coevolutionary arms race. Despite decades of research, however, we have only a handful of examples where reciprocal coevolutionary change has been rigorously established as the cause of trait exaggeration. Support for a coevolutionary mechanism remains elusive because we lack generally applicable tools for quantifying the intensity of coevolutionary selection. Here we develop an approximate Bayesian computation (ABC) approach for estimating the intensity of coevolutionary selection using population mean phenotypes of traits mediating interspecific interactions. Our approach relaxes important assumptions of a previous maximum likelihood approach by allowing gene flow among populations, variable abiotic environments, and strong coevolutionary selection. Using simulated data, we show that our ABC method accurately infers the strength of coevolutionary selection if reliable estimates are available for key background parameters and ten or more populations are sampled. Applying our approach to the putative arms race between the plant Camellia japonica and its seed predatory weevil, Curculio camelliae, provides support for a coevolutionary hypothesis but fails to preclude the possibility of unilateral evolution. Comparing independently estimated selection gradients acting on Camellia pericarp thickness with values simulated by our model reveals a correlation between predicted and observed selection gradients of 0.941. The strong agreement between predicted and observed selection gradients validates our method.

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.

The geographic mosaic of coevolution in mutualistic networks

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

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

Seed Predators, not Herbivores, Exert Natural Selection on Solidago spp. in an Urban Archipelago

The effects of urbanization on biodiversity are well established, as a growing city will reduce the size and diversity of patches of native plants. Recolonization of old patches and discovery of new ones by arthropod herbivores should occur as predicted by island biogeography theory. Although colonization represents an increase in biodiversity, such arrivals may exert new forms of natural selection on plants through herbivory and seed predation. Using a single species of old-field aster (Solidago altissima L.), we found that the level of natural selection by seed predators and herbivores follows patterns of island biogeography, with lower amounts of damage on smaller islands, where there are fewer species, and hypothetically smaller populations of arthropods. We also found that in an urban system, levels of herbivory are far below the tolerance levels of Solidago, and that seed predators are likely to be the only arthropod to cause reduced fitness. The pattern seen also implies that as a patch of Solidago grows through clonal expansion, it will come under higher selective pressure.

Phylogeography of specialist weevil Trichobaris soror: a seed predator of Datura stramonium

Can the genetic structure of a specialist weevil be explained by the geological history of their distribution zone? We analyze the genetic variation of the weevil Trichobaris soror, a specialist seed predator of Datura stramonium, in order to address this question. For the phylogeographic analysis we used the COI gene, and assessed species identity in weevil populations through geometric morphometric approach. In total, we found 53 haplotypes in 413 samples, whose genetic variation supports the formation of three groups: (1) the Transmexican Volcanic Belt (TVB group), (2) the Sierra Madre Sur (SMS group) and (3) the Balsas Basin (BB group). The morphometric analysis suggests that BB group is probably not T. soror. Our results have two implications: first, the phylogeographic pattern of T. soror is explained by both the formation of the geological provinces where it is currently distributed and the coevolution with its host plant, because the TVB and SMS groups could be separated due to the discontinuity of altitude between the geological provinces, but the recent population expansion of TVB group and the high frequency of only one haplotype can be due to specialization to the host plant. Second, we report a new record of a different species of weevil in BB group parasitizing D. stramonium fruits.

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

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

Local adaptation: simultaneously considering herbivores and their host plants

• Although a major expectation of coevolutionary theory between plants and herbivores is the occurrence of reciprocal local adaptation, this has remained almost untested. Thus, we evaluated the presence and variation in the patterns of reciprocal local adaptation between an herbivorous insect and its host plant. • Two four-by-four cross-infestation experiments were performed under similar abiotic conditions. The first one was done under laboratory conditions to estimate herbivore individual performance while the second one was performed in a common garden to simultaneously estimate herbivore population growth rate as well as seed production and plant defenses (resistance and tolerance to herbivory). • The patterns of population differentiation for the herbivore and the plant were not independent of each other, showing all the possible outcomes from locally adapted to maladapted populations. These results indicate differences in the magnitude of local adaptation. While an association between resistance and herbivore performance was observed, there was no clear pattern between tolerance and herbivore local adaptation. • Our results demonstrated the occurrence of reciprocal local adaptation following the pattern expected by theory: when the herbivores or the plants were adapted, the other species was non-adapted or even maladapted.

Arms race between weevil rostrum length and camellia pericarp thickness: Geographical cline and theory

The geographical cline of the coevolving traits of weevil rostrum (mouthpart) length and camellia pericarp (fruit coat) thickness provides an opportunity to test the arms race theory of defense (pericarp thickness) and countermeasure (rostrum length) between antagonistically interacting species. By extending the previous model for the coevolution of quantitative traits to introduce nonlinear costs for exaggerated traits, the generation overlap, and density-dependent regulation in the host, we studied the evolutionarily stable (ES) pericarp thickness in the Japanese camellia (Camellia japonica) and the ES rostrum length in the camellia-weevil (Curculio camelliae). The joint monomorphic ES system has a robust outcome with nonlinear costs, and we analyzed how the traits of both species at evolutionary equilibrium depend on demographic parameters. If camellia demographic parameters vary latitudinally, data collected over the geographical scale of rostrum length and pericarp thickness should lie on an approximately linear curve with the slope less than that of the equiprobability line A/B of boring success, where A and B are coefficients for the logistic regression of boring success to pericarp thickness and rostrum length, respectively. This is a robust prediction as long as the cost of rostrum length is nonlinear (accelerating). As a result, boring success should be lower in populations with longer rostrum length, as reported in the weevil-camellia system (Toju, H., and Sota, T., 2006a. Imbalance of predator and prey armament: Geographic clines in phenotypic interface and natural selection. American Naturalist 167, 105-117). The nonlinearity (exponent) for the cost of rostrum length estimated from the geographical cline data for the weevil-camellia system was 2.2, suggesting nonlinearity between quadratic and cubic forms.