Molecular signatures of selection associated with host plant differences in Pieris butterflies.

Adaptive traits that enable organisms to conquer novel niches and experience subsequent diversification are ecologically and evolutionarily important. The larvae of Pieris butterflies express nitrile-specifier proteins (NSPs), a key innovation for overcoming the glucosinolate (GLS)-myrosinase-based defence system of their Brassicales host plants. Nitrile-specifier proteins are a member of the NSP-like gene family, which includes the major allergen (MA) protein, a paralog of NSP with a GLS-disarming function, and a single domain major allergen (SDMA) protein, whose function is unknown. The arms-race between GLS-based defences and the NSP-like gene family is suggested to mediate diversification in both Pierid butterflies and Brassicales plants. Here, we tested whether the expected strong selection on NSP-like gene family correlates with shifts in host plant spectra among Pierid butterflies. We combined feeding experiments using 25 Brassicaceae plants and five Pieris species with larval transcriptome data to investigate the patterns of selection acting on NSP-like gene family members. Although we observed significantly elevated nonsynonymous to synonymous substitution rate ratios in NSPs on branches associated with changes in patterns of host plant usage, no such pattern was observed in MAs or SDMAs. Furthermore, we found evidence for positive selection of NSP at a phylogenetic branch which reflects different host plant spectra. Our data indicate that the NSP-related gene members have evolved differently: NSPs have accumulated more amino acid changes in response to shifting preferences for host plants, whereas MAs and SDMAs appear to be more conserved. Further detailed functional assays of these genes would provide important insights to understand their role in the chemical arms-race between Pieris butterflies and their Brassicales host plants.

Differential regulation of host plant adaptive genes in Pieris butterflies exposed to a range of glucosinolate profiles in their host plants.

Specialist herbivores have often evolved highly sophisticated mechanisms to counteract defenses mediated by major plant secondary-metabolites. Plant species of the herbivore host range often display high chemical diversity and it is not well understood how specialist herbivores respond to this chemical diversity. Pieris larvae overcome toxic products from glucosinolate hydrolysis, the major chemical defense of their Brassicaceae hosts, by expressing nitrile-specifier proteins (NSP) in their gut. Furthermore, Pieris butterflies possess so-called major allergen (MA) proteins, which are multi-domain variants of a single domain major allergen (SDMA) protein expressed in the guts of Lepidopteran larvae. Here we show that Pieris larvae fine-tune NSP and MA gene expression depending on the glucosinolate profiles of their Brassicaceae hosts. Although the role of MA is not yet fully understood, the expression levels of NSP and MA in larvae that fed on plants whose glucosinolate composition varied was dramatically changed, whereas levels of SDMA expression remained unchanged. In addition, we found a similar regulation pattern among these genes in larvae feeding on Arabidopsis mutants with different glucosinolate profiles. Our results demonstrate that Pieris larvae appear to use different host plant adaptive genes to overcome a wide range of glucosinolate profiles in their host plants.

An Insect Counteradaptation against Host Plant Defenses Evolved through Concerted Neofunctionalization.

Antagonistic chemical interactions between herbivorous insects and their host plants are often thought to coevolve in a stepwise process, with an evolutionary innovation on one side being countered by a corresponding advance on the other. Glucosinolate sulfatase (GSS) enzyme activity is essential for the Diamondback moth, Plutella xylostella, to overcome a highly diversified secondary metabolite-based host defense system in the Brassicales. GSS genes are located in an ancient cluster of arylsulfataselike genes, but the exact roles of gene copies and their evolutionary trajectories are unknown. Here, we combine a functional investigation of duplicated insect arylsulfatases with an analysis of associated nucleotide substitution patterns. We show that the Diamondback moth genome encodes three GSSs with distinct substrate spectra and distinct expression patterns in response to glucosinolates. Contrary to our expectations, early functional diversification of gene copies was not indicative of a coevolutionary arms race between host and herbivore. Instead, both copies of a duplicated arylsulfatase gene evolved concertedly in the context of an insect host shift to acquire novel detoxifying functions under positive selection, a pattern of duplicate gene retention that we call "concerted neofunctionalization."

Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore.

Phytophagous insects have to deal with toxic defense compounds from their host plants. Although it is known that insects have evolved genes and mechanisms to detoxify plant allochemicals, how specialist and generalist precisely respond to specific secondary metabolites at the molecular level is less understood. Here we studied the larval performance and transcriptome of the generalist moth Heliothis virescens and the specialist butterfly Pieris brassicae feeding on Arabidopsis thaliana genotypes with different glucosinolate (GS) levels. H. virescens larvae gained significantly more weight on the GS-deficient mutant quadGS compared to wild-type (Col-0) plants. On the contrary, P. brassicae was unaffected by the presence of GS and performed equally well on both genotypes. Strikingly, there was a considerable differential gene expression in H. virescens larvae feeding on Col-0 compared to quadGS. In contrast, compared to H. virescens, P. brassicae displayed a much-reduced transcriptional activation when fed on both plant genotypes. Transcripts coding for putative detoxification enzymes were significantly upregulated in H. virescens, along with digestive enzymes and transposable elements. These data provide an unprecedented view on transcriptional changes that are specifically activated by GS and illustrate differential molecular responses that are linked to adaptation to diet in lepidopteran herbivores.

Gossypol toxicity and detoxification in Helicoverpa armigera and Heliothis virescens.

Gossypol is a polyphenolic secondary metabolite produced by cotton plants, which is toxic to many organisms. Gossypol's aldehyde groups are especially reactive, forming Schiff bases with amino acids of proteins and cross-linking them, inhibiting enzyme activities and contributing to toxicity. Very little is known about gossypol's mode of action and its detoxification in cotton-feeding insects that can tolerate certain concentrations of this compound. Here, we tested the toxicity of gossypol and a gossypol derivative lacking free aldehyde groups (SB-gossypol) toward Helicoverpa armigera and Heliothis virescens, two important pests on cotton plants. Larval feeding studies with these two species on artificial diet supplemented with gossypol or SB-gossypol revealed no detectable toxicity of gossypol, when the aldehyde groups were absent. A cytochrome P450 enzyme, CYP6AE14, is upregulated in H. armigera feeding on gossypol, and has been claimed to directly detoxify gossypol. However, using in vitro assays with heterologously expressed CYP6AE14, no metabolites of gossypol were detected, and further studies suggest that gossypol is not a direct substrate of CYP6AE14. Furthermore, larvae feeding on many other plant toxins also upregulate CYP6AE14. Our data demonstrate that the aldehyde groups are critical for the toxicity of gossypol when ingested by H. armigera and H. virescens larvae, and suggest that CYP6AE14 is not directly involved in gossypol metabolism, but may play a role in the general stress response of H. armigera larvae toward plant toxins.

Immune modulation enables a specialist insect to benefit from antibacterial withanolides in its host plant.

The development of novel plant chemical defenses and counter adaptations by herbivorous insect could continually drive speciation, producing more insect specialists than generalists. One approach to test this hypothesis is to compare closely related generalist and specialist species to reveal the associated costs and benefits of these different adaptive strategies. We use the specialized moth Heliothis subflexa, which feeds exclusively on plants in the genus Physalis, and its close generalist relative H. virescens. Specialization on Physalis plants necessitates the ability to tolerate withanolides, the secondary metabolites of Physalis species that are known to have feeding deterrent and immune inhibiting properties for other insects. Here we find that only H. subflexa benefits from the antibacterial properties of withanolides, and thereby gains a higher tolerance of the pathogen Bacillus thuringiensis. We argue that the specialization in H. subflexa has been guided to a large extent by a unique role of plant chemistry on ecological immunology.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta.

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Potential detoxification of gossypol by UDP-glycosyltransferases in the two Heliothine moth species Helicoverpa armigera and Heliothis virescens.

The cotton bollworm Helicoverpa armigera and the tobacco budworm Heliothis virescens are closely related generalist insect herbivores and serious pest species on a number of economically important crop plants including cotton. Even though cotton is well defended by its major defensive compound gossypol, a toxic sesquiterpene dimer, larvae of both species are capable of developing on cotton plants. In spite of severe damage larvae cause on cotton plants, little is known about gossypol detoxification mechanisms in cotton-feeding insects. Here, we detected three monoglycosylated and up to five diglycosylated gossypol isomers in the feces of H. armigera and H. virescens larvae fed on gossypol-supplemented diet. Candidate UDP-glycosyltransferase (UGT) genes of H. armigera were selected by microarray studies and in silico analyses and were functionally expressed in insect cells. In enzymatic assays, we show that UGT41B3 and UGT40D1 are capable of glycosylating gossypol mainly to the diglycosylated gossypol isomer 5 that is characteristic for H. armigera and is absent in H. virescens feces. In conclusion, our results demonstrate that gossypol is partially metabolized by UGTs via glycosylation, which might be a crucial step in gossypol detoxification in generalist herbivores utilizing cotton as host plant.

The butterfly plant arms-race escalated by gene and genome duplications.

Coevolutionary interactions are thought to have spurred the evolution of key innovations and driven the diversification of much of life on Earth. However, the genetic and evolutionary basis of the innovations that facilitate such interactions remains poorly understood. We examined the coevolutionary interactions between plants (Brassicales) and butterflies (Pieridae), and uncovered evidence for an escalating evolutionary arms-race. Although gradual changes in trait complexity appear to have been facilitated by allelic turnover, key innovations are associated with gene and genome duplications. Furthermore, we show that the origins of both chemical defenses and of molecular counter adaptations were associated with shifts in diversification rates during the arms-race. These findings provide an important connection between the origins of biodiversity, coevolution, and the role of gene and genome duplications as a substrate for novel traits.

Molecular mechanisms of insect adaptation to plant secondary compounds.

During feeding, herbivorous insects are exposed to an array of plant defensive compounds. In this review, we examine molecular mechanisms of insect adaptation to these toxic metabolites. We discuss both the importance of evolutionary variation of existing detoxification gene families, as well as the evolution of novel mechanisms through gene recruitment, neofunctionalization and horizontal gene transfer. The ability of insects to cope with the chemical diversity of their host plants and the different mechanisms that insects use to resist these toxins open new avenues for understanding fundamental aspects of insect-plant coevolutionary adaptation.

Microevolutionary dynamics of a macroevolutionary key innovation in a Lepidopteran herbivore.

BACKGROUND: A molecular population genetics understanding is central to the study of ecological and evolutionary functional genomics. Population genetics identifies genetic variation and its distribution within and among populations, it reveals the demographic history of the populations studied, and can provide indirect insights into historical selection dynamics. Here we use this approach to examine the demographic and selective dynamics acting of a candidate gene involved in plant-insect interactions. Previous work documents the macroevolutionary and historical ecological importance of the nitrile-specifier protein (Nsp), which facilitated the host shift of Pieridae butterflies onto Brassicales host plants approximately 80 Myr ago. RESULTS: Here we assess the microevolutionary dynamics of the Nsp gene by studying the within and among-population variation at Nsp and reference genes in the butterfly Pieris rapae (Small Cabbage White). Nsp exhibits unexpectedly high amounts of amino acid polymorphism, unequally distributed across the gene. The vast majority of genetic variation exists within populations, with little to no genetic differentiation among four populations on two continents. A comparison of synonymous and nonsynonymous substitutions in 70 randomly chosen genes among P. rapae and its close relative Pieris brassicae (Large Cabbage White) finds Nsp to have a significantly relaxed functional constraint compared to housekeeping genes. We find strong evidence for a recent population expansion and no role for strong purifying or directional selection upon the Nsp gene. CONCLUSIONS: The microevolutionary dynamics of the Nsp gene in P. rapae are dominated by recent population expansion and variation in functional constraint across the repeated domains of the Nsp gene. While the high amounts of amino acid diversity suggest there may be significant functional differences among allelic variants segregating within populations, indirect tests of selection could not conclusively identify a signature of historical selection. The importance of using this information for planning future studies of potential performance and fitness consequences of the observed variation is discussed.

Phylogenetic relatedness and host plant growth form influence gene expression of the polyphagous comma butterfly (Polygonia c-album).

BACKGROUND: The mechanisms that shape the host plant range of herbivorous insect are to date not well understood but knowledge of these mechanisms and the selective forces that influence them can expand our understanding of the larger ecological interaction. Nevertheless, it is well established that chemical defenses of plants influence the host range of herbivorous insects. While host plant chemistry is influenced by phylogeny, also the growth forms of plants appear to influence the plant defense strategies as first postulated by Feeny (the "plant apparency" hypothesis). In the present study we aim to investigate the molecular basis of the diverse host plant range of the comma butterfly (Polygonia c-album) by testing differential gene expression in the caterpillars on three host plants that are either closely related or share the same growth form. RESULTS: In total 120 genes were identified to be differentially expressed in P. c-album after feeding on different host plants, 55 of them in the midgut and 65 in the restbody of the caterpillars. Expression patterns could be confirmed with an independent method for 14 of 27 tested genes. Pairwise similarities in upregulation in the midgut of the caterpillars were higher between plants that shared either growth form or were phylogenetically related. No known detoxifying enzymes were found to be differently regulated in the midgut after feeding on different host plants. CONCLUSION: Our data suggest a complex picture of gene expression in response to host plant feeding. While each plant requires a unique gene regulation in the caterpillar, both phylogenetic relatedness and host plant growth form appear to influence the expression profile of the polyphagous comma butterfly, in agreement with phylogenetic studies of host plant utilization in butterflies.

Immunity or digestion: glucanase activity in a glucan-binding protein family from Lepidoptera.

The cell surfaces of microorganisms display distinct molecular patterns formed from lipopolysaccharides, peptidoglycans, or beta1,3-glucans. Binding of these surfaces by pattern recognition proteins such as beta1,3-glucan recognition proteins (betaGRPs) activates the immune response in arthropods. We identified a 40-kDa beta1,3-glucan-binding protein with sequence similarity to previously characterized lepidopteran betaGRPs from hemolymph, but unlike these it is secreted into the larval gut lumen and is an active beta1,3-glucanase. This glucanase was not detected in hemolymph. Its mRNA is constitutively and predominantly expressed in the midgut and is induced there when larvae feed on a diet containing bacteria. Homologs of this predominantly midgut-expressed gene from many Lepidoptera possess key residues shown to be part of the active site of other glucanases, and form a cluster that is distinct from previously described betaGRPs. In addition, this group includes proteins from insects such as the Anopheles gambiae GNBP subgroup B for which a catalytic role has not been previously suspected. The current domain classification does not distinguish between the catalytic and noncatalytic clades, and should be revised. The noncatalytic betaGRPs may be evolutionarily derived from this newly described enzyme family that continues to function catalytically in digestion and/or pathogen defense.