Larval transcriptomes reflect the evolutionary history of plant-insect associations.

In this study, we investigated whether patterns of gene expression in larvae feeding on different plants can explain important aspects of the evolution of insect-plant associations, such as phylogenetic conservatism of host use and re-colonization of ancestral hosts that have been lost from the host repertoire. To this end, we performed a phylogenetically informed study comparing the transcriptomes of 4 nymphalid butterfly species in Polygonia and the closely related genus Nymphalis. Larvae were reared on Urtica dioica, Salix spp., and Ribes spp. Plant-specific gene expression was found to be similar across butterfly species, even in the case of host plants that are no longer used by two of the butterfly species. These results suggest that plant-specific transcriptomes can be robust over evolutionary time. We propose that adaptations to particular larval food plants can profitably be understood as an evolved set of modules of co-expressed genes, promoting conservatism in host use and facilitating re-colonization. Moreover, we speculate that the degree of overlap between plant-specific transcriptomes may correlate with the strength of trade-offs between plants as resources and hence to the probability of colonizing hosts and complete host shifts.

Phylogenetic reconstruction of ancestral ecological networks through time for pierid butterflies and their host plants.

The study of herbivorous insects underpins much of the theory that concerns the evolution of species interactions. In particular, Pieridae butterflies and their host plants have served as a model system for studying evolutionary arms races. To learn more about the coevolution of these two clades, we reconstructed ancestral ecological networks using stochastic mappings that were generated by a phylogenetic model of host-repertoire evolution. We then measured if, when, and how two ecologically important structural features of the ancestral networks (modularity and nestedness) evolved over time. Our study shows that as pierids gained new hosts and formed new modules, a subset of them retained or recolonised the ancestral host(s), preserving connectivity to the original modules. Together, host-range expansions and recolonisations promoted a phase transition in network structure. Our results demonstrate the power of combining network analysis with Bayesian inference of host-repertoire evolution to understand changes in complex species interactions over time.

Bayesian Inference of Ancestral Host-Parasite Interactions under a Phylogenetic Model of Host Repertoire Evolution.

Intimate ecological interactions, such as those between parasites and their hosts, may persist over long time spans, coupling the evolutionary histories of the lineages involved. Most methods that reconstruct the coevolutionary history of such interactions make the simplifying assumption that parasites have a single host. Many methods also focus on congruence between host and parasite phylogenies, using cospeciation as the null model. However, there is an increasing body of evidence suggesting that the host ranges of parasites are more complex: that host ranges often include more than one host and evolve via gains and losses of hosts rather than through cospeciation alone. Here, we develop a Bayesian approach for inferring coevolutionary history based on a model accommodating these complexities. Specifically, a parasite is assumed to have a host repertoire, which includes both potential hosts and one or more actual hosts. Over time, potential hosts can be added or lost, and potential hosts can develop into actual hosts or vice versa. Thus, host colonization is modeled as a two-step process that may potentially be influenced by host relatedness. We first explore the statistical behavior of our model by simulating evolution of host-parasite interactions under a range of parameter values. We then use our approach, implemented in the program RevBayes, to infer the coevolutionary history between 34 Nymphalini butterfly species and 25 angiosperm families. Our analysis suggests that host relatedness among angiosperm families influences how easily Nymphalini lineages gain new hosts. [Ancestral hosts; coevolution; herbivorous insects; probabilistic modeling.].

Structural plasticity of olfactory neuropils in relation to insect diapause.

Many insects that live in temperate zones spend the cold season in a state of dormancy, referred to as diapause. As the insect must rely on resources that were gathered before entering diapause, keeping a low metabolic rate is of utmost importance. Organs that are metabolically expensive to maintain, such as the brain, can therefore become a liability to survival if they are too large.Insects that go through diapause as adults generally do so before entering the season of reproduction. This order of events introduces a conflict between maintaining low metabolism during dormancy and emerging afterward with highly developed sensory systems that improve fitness during the mating season.We investigated the timing of when investments into the olfactory system are made by measuring the volumes of primary and secondary olfactory neuropils in the brain as they fluctuate in size throughout the extended diapause life-period of adult Polygonia c-album butterflies.Relative volumes of both olfactory neuropils increase significantly during early adult development, indicating the importance of olfaction to this species, but still remain considerably smaller than those of nondiapausing conspecifics. However, despite butterflies being kept under the same conditions as before the dormancy, their olfactory neuropil volumes decreased significantly during the postdormancy period.The opposing directions of change in relative neuropil volumes before and after diapause dormancy indicate that the investment strategies governing structural plasticity during the two life stages could be functionally distinct. As butterflies were kept in stimulus-poor conditions, we find it likely that investments into these brain regions rely on experience-expectant processes before diapause and experience-dependent processes after diapause conditions are broken.As the shift in investment strategies coincides with a hard shift from premating season to mating season, we argue that these developmental characteristics could be adaptations that mitigate the trade-off between dormancy survival and reproductive fitness.

Unifying host-associated diversification processes using butterfly-plant networks.

Explaining the exceptional diversity of herbivorous insects is an old problem in evolutionary ecology. Here we focus on the two prominent hypothesised drivers of their diversification, radiations after major host switch or variability in host use due to continuous probing of new hosts. Unfortunately, current methods cannot distinguish between these hypotheses, causing controversy in the literature. Here we present an approach combining network and phylogenetic analyses, which directly quantifies support for these opposing hypotheses. After demonstrating that each hypothesis produces divergent network structures, we then investigate the contribution of each to diversification in two butterfly families: Pieridae and Nymphalidae. Overall, we find that variability in host use is essential for butterfly diversification, while radiations following colonisation of a new host are rare but can produce high diversity. Beyond providing an important reconciliation of alternative hypotheses for butterfly diversification, our approach has potential to test many other hypotheses in evolutionary biology.

Host use dynamics in a heterogeneous fitness landscape generates oscillations in host range and diversification.

Colonization of novel hosts is thought to play an important role in parasite diversification, yet little consensus has been achieved about the macroevolutionary consequences of changes in host use. Here, we offer a mechanistic basis for the origins of parasite diversity by simulating lineages evolved in silico. We describe an individual-based model in which (i) parasites undergo sexual reproduction limited by genetic proximity, (ii) hosts are uniformly distributed along a one-dimensional resource gradient, and (iii) host use is determined by the interaction between the phenotype of the parasite and a heterogeneous fitness landscape. We found two main effects of host use on the evolution of a parasite lineage. First, the colonization of a novel host allowed parasites to explore new areas of the resource space, increasing phenotypic and genotypic variation. Second, hosts produced heterogeneity in the parasite fitness landscape, which led to reproductive isolation and therefore, speciation. As a validation of the model, we analyzed empirical data from Nymphalidae butterflies and their host plants. We then assessed the number of hosts used by parasite lineages and the diversity of resources they encompass. In both simulated and empirical systems, host diversity emerged as the main predictor of parasite species richness.

Embracing Colonizations: A New Paradigm for Species Association Dynamics.

Parasite-host and insect-plant research have divergent traditions despite the fact that most phytophagous insects live parasitically on their host plants. In parasitology it is a traditional assumption that parasites are typically highly specialized; cospeciation between parasites and hosts is a frequently expressed default expectation. Insect-plant theory has been more concerned with host shifts than with cospeciation, and more with hierarchies among hosts than with extreme specialization. We suggest that the divergent assumptions in the respective fields have hidden a fundamental similarity with an important role for potential as well as actual hosts, and hence for host colonizations via ecological fitting. A common research program is proposed which better prepares us for the challenges from introduced species and global change.

Experience-dependent mushroom body plasticity in butterflies: consequences of search complexity and host range.

An ovipositing insect experiences many sensory challenges during her search for a suitable host plant. These sensory challenges become exceedingly pronounced when host range increases, as larger varieties of sensory inputs have to be perceived and processed in the brain. Neural capacities can be exceeded upon information overload, inflicting costs on oviposition accuracy. One presumed generalist strategy to diminish information overload is the acquisition of a focused search during its lifetime based on experiences within the current environment, a strategy opposed to a more genetically determined focus expected to be seen in relative specialists. We hypothesized that a broader host range is positively correlated with mushroom body (MB) plasticity, a brain structure related to learning and memory. To test this hypothesis, butterflies with diverging host ranges (Polygonia c-album, Aglais io and Aglais urticae) were subjected to differential environmental complexities for oviposition, after which ontogenetic MB calyx volume differences were compared among species. We found that the relative generalist species exhibited remarkable plasticity in ontogenetic MB volumes; MB growth was differentially stimulated based on the complexity of the experienced environment. For relative specialists, MB volume was more canalized. All in all, this study strongly suggests an impact of host range on brain plasticity in Nymphalid butterflies.

Polyphagy and diversification in tussock moths: Support for the oscillation hypothesis from extreme generalists.

Theory on plasticity driving speciation, as applied to insect-plant interactions (the oscillation hypothesis), predicts more species in clades with higher diversity of host use, all else being equal. Previous support comes mainly from specialized herbivores such as butterflies, and plasticity theory suggests that there may be an upper host range limit where host diversity no longer promotes diversification. The tussock moths (Erebidae: Lymantriinae) are known for extreme levels of polyphagy. We demonstrate that this system is also very different from butterflies in terms of phylogenetic signal for polyphagy and for use of specific host orders. Yet we found support for the generality of the oscillation hypothesis, in that clades with higher diversity of host use were found to contain more species. These clades also consistently contained the most polyphagous single species. Comparing host use in Lymantriinae with related taxa shows that the taxon indeed stands out in terms of the frequency of polyphagous species. Comparative evidence suggests that this is most probably due to its nonfeeding adults, with polyphagy being part of a resulting life history syndrome. Our results indicate that even high levels of plasticity can drive diversification, at least when the levels oscillate over time.

Why stay in a bad relationship? The effect of local host phenology on a generalist butterfly feeding on a low-ranked host.

BACKGROUND: In plant-feeding insects, the evolutionary retention of polyphagy remains puzzling. A better understanding of the relationship between these organisms and changes in the metabolome of their host plants is likely to suggest functional links between them, and may provide insights into how polyphagy is maintained. RESULTS: We investigated the phenological change of Cynoglossum officinale, and how a generalist butterfly species, Vanessa cardui, responded to this change. We used untargeted metabolite profiling to map plant seasonal changes in both primary and secondary metabolites. We compared these data to differences in larval performance on vegetative plants early and late in the season. We also performed two oviposition preference experiments to test females' ability to choose between plant developmental stages (vegetative and reproductive) early and late in the season. We found clear seasonal changes in plant primary and secondary metabolites that correlated with larval performance. The seasonal change in plant metabolome reflected changes in both nutrition and toxicity and resulted in zero survival in the late period. However, large differences among families in larval ability to feed on C. officinale suggest that there is genetic variation for performance on this host. Moreover, females accepted all plants for oviposition, and were not able to discriminate between plant developmental stages, in spite of the observed overall differences in metabolite profile potentially associated with differences in suitability as larval food. CONCLUSIONS: In V. cardui, migratory behavior, and thus larval feeding times, are not synchronized with plant phenology at the reproductive site. This lack of synchronization, coupled with the observed lack of discriminatory oviposition, obviously has potential fitness costs. However, this "opportunistic" behavior may as well function as a source of potential host plant evolution, promoting for example the acceptance of new plants.

On oscillations and flutterings-A reply to Hamm and Fordyce.

The diversification of plant-feeding insects is seen as a spectacular example of evolutionary radiation. Hence, developing hypotheses to explain this diversification, and methods to test them, is an important undertaking. Some years ago, we presented the oscillation hypothesis as a general process that could drive diversification of this and similar interactions, through repeated expansions and contractions of host ranges. Hamm and Fordyce recently presented a study with the outspoken intention of testing this hypothesis where they concluded that the oscillation hypothesis was not supported. We point out several problems with their study, owing both to a misrepresentation of our hypothesis and to the methods. We provide a clarifying description of the oscillation hypothesis, and detail some predictions that follow from it. A reanalysis of the data demonstrated a troubling sensitivity of the "SSE" class of models to small changes in model specification, and we caution against using them for tests of trait-based diversification. Future tests of the hypothesis also need to better acknowledge the processes behind the host range oscillations. We suspect that doing so will resolve some of the apparent conflicts between our hypothesis and the view presented by Hamm and Fordyce.

It's All in the Mix: Blend-Specific Behavioral Response to a Sexual Pheromone in a Butterfly.

Among insects, sexual pheromones are typically mixtures of two to several components, all of which are generally required to elicit a behavioral response. Here we show for the first time that a complete blend of sexual pheromone components is needed to elicit a response also in a butterfly. Males of the Green-veined White, Pieris napi, emit an aphrodisiac pheromone, citral, from wing glands. This pheromone is requisite for females to accept mating with a courting male. Citral is a mixture of the two geometric isomers geranial (E-isomer) and neral (Z-isomer) in an approximate 1:1 ratio. We found that both these compounds are required to elicit acceptance behavior, which indicates synergistic interaction between processing of the isomers. Using functional Ca(2+) imaging we found that geranial and neral evoke significantly different but overlapping glomerular activity patterns in the antennal lobe, which suggests receptors with different affinity for the two isomers. However, these glomeruli were intermingled with glomeruli responding to, for example, plant-related compounds, i.e., no distinct subpopulation of pheromone-responding glomeruli as in moths and other insects. In addition, these glomeruli showed lower specificity than pheromone-activated glomeruli in moths. We could, however, not detect any mixture interactions among four identified glomeruli, indicating that the synergistic effect may be generated at a higher processing level. Furthermore, correlations between glomerular activity patterns evoked by the single isomers and the blend did not change over time.

Evolutionary history of host use, rather than plant phylogeny, determines gene expression in a generalist butterfly.

BACKGROUND: Although most insect species are specialized on one or few groups of plants, there are phytophagous insects that seem to use virtually any kind of plant as food. Understanding the nature of this ability to feed on a wide repertoire of plants is crucial for the control of pest species and for the elucidation of the macroevolutionary mechanisms of speciation and diversification of insect herbivores. Here we studied Vanessa cardui, the species with the widest diet breadth among butterflies and a potential insect pest, by comparing tissue-specific transcriptomes from caterpillars that were reared on different host plants. We tested whether the similarities of gene-expression response reflect the evolutionary history of adaptation to these plants in the Vanessa and related genera, against the null hypothesis of transcriptional profiles reflecting plant phylogenetic relatedness. RESULT: Using both unsupervised and supervised methods of data analysis, we found that the tissue-specific patterns of caterpillar gene expression are better explained by the evolutionary history of adaptation of the insects to the plants than by plant phylogeny. CONCLUSION: Our findings suggest that V. cardui may use two sets of expressed genes to achieve polyphagy, one associated with the ancestral capability to consume Rosids and Asterids, and another allowing the caterpillar to incorporate a wide range of novel host-plants.

Specialist and generalist oviposition strategies in butterflies: maternal care or precocious young?

Herbivorous insects specialized on a narrow set of plants are believed to be better adapted to their specific hosts. This hypothesis is supported by observations of herbivorous insect species with a broader diet breadth which seemingly pay a cost through decreased oviposition accuracy. Despite many studies investigating female oviposition behavior, there is a lack of knowledge on how larvae cope behaviorally with their mothers' egg-laying strategies. We have examined a unique system of five nymphalid butterfly species with different host plant ranges that all feed on the same host plant. The study of this system allowed us to compare at the species level how oviposition preference is related to neonate larval responses in several disadvantageous situations. We found a general co-adaptation between female and larval abilities, where species with more discriminating females had larvae that were less able to deal with a suboptimal initial feeding site. Conversely, relatively indiscriminate females had more precocious larvae with better abilities to cope with suboptimal sites. Despite similarities between the tested species with similar host ranges, there were also striking differences. Generalist and specialist species can be found side by side in many clades, with each clade having a specific evolutionary history. Such clade-specific, phylogenetically determined preconditions apparently have affected how precisely a broad or narrow diet breadth can be realized.

Plant fertilization interacts with life history: variation in stoichiometry and performance in nettle-feeding butterflies.

Variation in food stoichiometry affects individual performance and population dynamics, but it is also likely that species with different life histories should differ in their sensitivity to food stoichiometry. To address this question, we investigated the ability of the three nettle-feeding butterflies (Aglais urticae, Polygonia c-album, and Aglais io) to respond adaptively to induced variation in plant stoichiometry in terms of larval performance. We hypothesized that variation in larval performance between plant fertilization treatments should be functionally linked to species differences in host plant specificity. We found species-specific differences in larval performance between plant fertilization treatments that could not be explained by nutrient limitation. We showed a clear evidence of a positive correlation between food stoichiometry and development time to pupal stage and pupal mass in A. urticae. The other two species showed a more complex response. Our results partly supported our prediction that host plant specificity affects larval sensitivity to food stoichiometry. However, we suggest that most of the differences observed may instead be explained by differences in voltinism (number of generations per year). We believe that the potential of some species to respond adaptively to variation in plant nutrient content needs further attention in the face of increased eutrophication due to nutrient leakage from human activities.

Effect of host plant and immune challenge on the levels of chemosensory and odorant-binding proteins in caterpillar salivary glands.

More than half of the proteome from mandibular glands in caterpillars is represented by chemosensory proteins. Based on sequence similarity, these proteins are putative transporters of ligands to gustatory receptors in sensory organs of insects. We sought to determine whether these proteins are inducible by comparing, both qualitatively and quantitatively, the salivary (mandibular and labial) proteomes from caterpillars (Vanessa cardui) reared on different plants and artificial diet containing either bacteria or bacterial cell-walls. We included a treatment where the caterpillars were switched from feeding on artificial diet to plant material at some point in their development. Additionally, we evaluated the degree of overlap between the proteomes in the hemolymph-filled coelom and salivary glands of caterpillars reared on plant material. We found that the quality and quantity of the identified proteins differed clearly between hemolymph-filled coelome, labial and mandibular glands. Our results indicated that even after molting and two-day feeding on a new diet, protein production is affected by the previous food source used by the caterpillar. Candidate proteins involved in chemosensory perception by insects were detected: three chemosensory (CSPs) and two odorant-binding proteins (OBPs). Using the relative amounts of these proteins across tissues and treatments as criteria for their classification, we detected hemolymph- and mandibular gland-specific CSPs and observed that their levels were affected by caterpillar diet. Moreover, we could compare the protein and transcript levels across tissues and treatment for at least one CSP and one OBP. Therefore, we have identified specific isoforms for testing the role of CSPs and OBPs in plant and pathogen recognition. We detected catalase, immune-related protein and serine proteases and their inhibitors in high relative levels in the mandibular glands in comparison to the labial glands. These findings suggest that the mandibular glands of caterpillars may play an important role protecting the caterpillar from oxidative stress, pathogens and aiding in digestion. Contamination with hemolymph proteins during dissection of salivary glands from caterpillars may occur but it is not substantial since the proteomes from hemolymph, mandibular and labial glands were easily discriminated from each other by principal component analysis of proteomic data.

Host plant choice in the comma butterfly-larval choosiness may ameliorate effects of indiscriminate oviposition.

In most phytophagous insects, the larval diet strongly affects future fitness and in species that do not feed on plant parts as adults, larval diet is the main source of nitrogen. In many of these insect-host plant systems, the immature larvae are considered to be fully dependent on the choice of the mothers, who, in turn, possess a highly developed host recognition system. This circumstance allows for a potential mother-offspring conflict, resulting in the female maximizing her fecundity at the expense of larval performance on suboptimal hosts. In two experiments, we aimed to investigate this relationship in the polyphagous comma butterfly, Polygonia c-album, by comparing the relative acceptance of low- and medium-ranked hosts between females and neonate larvae both within individuals between life stages, and between mothers and their offspring. The study shows a variation between females in oviposition acceptance of low-ranked hosts, and that the degree of acceptance in the mothers correlates with the probability of acceptance of the same host in the larvae. We also found a negative relationship between stages within individuals as there was a higher acceptance of lower ranked hosts in females who had abandoned said host as a larva. Notably, however, neonate larvae of the comma butterfly did not unconditionally accept to feed from the least favorable host species even when it was the only food source. Our results suggest the possibility that the disadvantages associated with a generalist oviposition strategy can be decreased by larval participation in host plant choice.

Host plant utilization, host range oscillations and diversification in nymphalid butterflies: a phylogenetic investigation.

It has been suggested that phenotypic plasticity is a major factor in the diversification of life, and that variation in host range in phytophagous insects is a good model for investigating this claim. We explore the use of angiosperm plants as hosts for nymphalid butterflies, and in particular the evidence for past oscillations in host range and how they are linked to host shifts and to diversification. At the level of orders of plants, a relatively simple pattern of host use and host shifts emerges, despite the 100 million years of history of the family Nymphalidae. We review the evidence that these host shifts and the accompanying diversifications were associated with transient polyphagous stages, as suggested by the "oscillation hypothesis." In addition, we investigate all currently polyphagous nymphalid species and demonstrate that the state of polyphagy is rare, has a weak phylogenetic signal, and a very apical distribution in the phylogeny; we argue that these are signs of its transient nature. We contrast our results with data from the bark beetles Dendroctonus, in which a more specialized host use is instead the apical state. We conclude that plasticity in host use is likely to have contributed to diversification in nymphalid butterflies.

Implications of a temperature increase for host plant range: predictions for a butterfly.

Although changes in phenology and species associations are relatively well-documented responses to global warming, the potential interactions between these phenomena are less well understood. In this study, we investigate the interactions between temperature, phenology (in terms of seasonal timing of larval growth) and host plant use in the polyphagous butterfly Polygonia c-album. We found that the hierarchy of larval performance on three natural host plants was not modified by a temperature increase as such. However, larval performance on each host plant and temperature treatment was affected by rearing season. Even though larvae performed better at the higher temperature regardless of the time of the rearing, relative differences between host plants changed with the season. For larvae reared late in the season, performance was always better on the herbaceous plant than on the woody plants. In this species, it is likely that a prolonged warming will lead to a shift from univoltinism to bivoltinism. The demonstrated interaction between host plant suitability and season means that such a shift is likely to lead to a shift in selective regime, favoring specialization on the herbaceous host. Based on our result, we suggest that host range evolution in response to temperature increase would in this species be highly contingent on whether the population undergoes a predicted shift from one to two generations. We discuss the effect of global warming on species associations and the outcome of asynchrony in rates of phenological change.

Mechanisms of macroevolution: polyphagous plasticity in butterfly larvae revealed by RNA-Seq.

Transcriptome studies of insect herbivory are still rare, yet studies in model systems have uncovered patterns of transcript regulation that appear to provide insights into how insect herbivores attain polyphagy, such as a general increase in expression breadth and regulation of ribosomal, digestion- and detoxification-related genes. We investigated the potential generality of these emerging patterns, in the Swedish comma, Polygonia c-album, which is a polyphagous, widely-distributed butterfly. Urtica dioica and Ribes uva-crispa are hosts of P. c-album, but Ribes represents a recent evolutionary shift onto a very divergent host. Utilizing the assembled transcriptome for read mapping, we assessed gene expression finding that caterpillar life-history (i.e. 2nd vs. 4th-instar regulation) had a limited influence on gene expression plasticity. In contrast, differential expression in response to host-plant identified genes encoding serine-type endopeptidases, membrane-associated proteins and transporters. Differential regulation of genes involved in nucleic acid binding was also observed suggesting that polyphagy involves large scale transcriptional changes. Additionally, transcripts coding for structural constituents of the cuticle were differentially expressed in caterpillars in response to their diet indicating that the insect cuticle may be a target for plant defence. Our results state that emerging patterns of transcript regulation from model species appear relevant in species when placed in an evolutionary context.