Cytochrome P450 diversification and hostplant utilization patterns in specialist and generalist moths: Birth, death and adaptation

MicroRNA bmo-miR-31-5p inhibits apoptosis and promotes BmNPV proliferation by targeting the CYP9e2 gene of Bombyx mori

BACKGROUND

Bombyx mori nuclear polyhedrosis virus (BmNPV), as a typical baculovirus, is the primary pathogen that infects the silkworm B. mori, a lepidopteran species. Owing to the high biological safety of BmNPV in infecting insects, it is commonly utilized as a biological insecticide for pest control. Apoptosis is important in the interaction between the host and pathogenic microorganisms. MicroRNAs (miRNAs) influence immune responses and promote stability of the immune system via apoptosis. Therefore, the study of apoptosis-related miRNA in silkworms during virus infection can not only provide support for standardizing the prevention and control of diseases and insect pests, but also reduce the economic losses to sericulture caused by the misuse of biological pesticides.

RESULTS

Through transcriptome sequencing, we identified a miRNA, miR-31-5p, and demonstrated that it can inhibit apoptosis in silkworm cells and promote the proliferation of BmNPV in BmE-SWU1 cells. We identified a target gene of miR-31-5p, B. mori cytochrome P450 9e2 (BmCYP9e2), and demonstrated that it can promote apoptosis in silkworm cells and inhibit the proliferation of BmNPV. Moreover, we constructed transgenic silkworm strains with miR-31-5p knockout and confirmed that they can inhibit the proliferation of BmNPV.

CONCLUSION

These data indicate that miR-31-5p may exert functions of inhibiting apoptosis and promoting virus proliferation by regulating BmCYP9e2. The findings demonstrate how miRNAs influence host cell apoptosis and how they are involved in the host immune system response to viruses, providing important insights into the applications of biological insecticides for pest control. © 2024 Society of Chemical Industry.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti

BACKGROUND

Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking.

RESULTS

In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds.

CONCLUSIONS

We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Girdling behavior of the longhorn beetle modulates the host plant to enhance larval performance

BACKGROUND

Preingestive behavioral modulations of herbivorous insects on the host plant are abundant over insect taxa. Those behaviors are suspected to have functions such as deactivation of host plant defenses, nutrient accumulation, or modulating plant-mediated herbivore interactions. To understand the functional consequence of behavioral modulation of insect herbivore, we studied the girdling behavior of Phytoecia rufiventris Gautier (Lamiinae; Cerambycidae) on its host plant Erigeron annuus L. (Asteraceae) that is performed before endophytic oviposition in the stem.

RESULTS

The girdling behavior significantly increased the larval performance in both field monitoring and lab experiment. The upper part of the girdled stem exhibited lack of jasmonic acid induction upon larval attack, lowered protease inhibitor activity, and accumulated sugars and amino acids in compared to non-girdled stem. The girdling behavior had no effect on the larval performance of a non-girdling longhorn beetle Agapanthia amurensis, which also feeds on the stem of E. annuus during larval phase. However, the girdling behavior decreased the preference of A. amurensis females for oviposition, which enabled P. rufiventris larvae to avoid competition with A. amurensis larvae.

CONCLUSIONS

In conclusion, the girdling behavior modulates plant physiology and morphology to provide a modulated food source for larva and hide it from the competitor. Our study implies that the insect behavior modulations can have multiple functions, providing insights into adaptation of insect behavior in context of plant-herbivore interaction.

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

Biodegradation of Polystyrene by Galleria mellonella: Identification of Potential Enzymes Involved in the Degradative Pathway

Galleria mellonella is a lepidopteran whose larval stage has shown the ability to degrade polystyrene (PS), one of the most recalcitrant plastics to biodegradation. In the present study, we fed G. mellonella larvae with PS for 54 days and determined candidate enzymes for its degradation. We first confirmed the biodegradation of PS by Fourier transform infrared spectroscopy- Attenuated total reflectance (FTIR-ATR) and then identified candidate enzymes in the larval gut by proteomic analysis using liquid chromatography with tandem mass spectrometry (LC-MS/MS). Two of these proteins have structural similarities to the styrene-degrading enzymes described so far. In addition, potential hydrolases, isomerases, dehydrogenases, and oxidases were identified that show little similarity to the bacterial enzymes that degrade styrene. However, their response to a diet based solely on polystyrene makes them interesting candidates as a potential new group of polystyrene-metabolizing enzymes in eukaryotes.

Evolutionary genomics of three agricultural pest moths reveals rapid evolution of host adaptation and immune-related genes

BACKGROUND

Understanding the genotype of pest species provides an important baseline for designing integrated pest management (IPM) strategies. Recently developed long-read sequence technologies make it possible to compare genomic features of nonmodel pest species to disclose the evolutionary path underlying the pest species profiles. Here we sequenced and assembled genomes for 3 agricultural pest gelechiid moths: Phthorimaea absoluta (tomato leafminer), Keiferia lycopersicella (tomato pinworm), and Scrobipalpa atriplicella (goosefoot groundling moth). We also compared genomes of tomato leafminer and tomato pinworm with published genomes of Phthorimaea operculella and Pectinophora gossypiella to investigate the gene family evolution related to the pest species profiles.

RESULTS

We found that the 3 solanaceous feeding species, P. absoluta, K. lycopersicella, and P. operculella, are clustered together. Gene family evolution analyses with the 4 species show clear gene family expansions on host plant-associated genes for the 3 solanaceous feeding species. These genes are involved in host compound sensing (e.g., gustatory receptors), detoxification (e.g., ABC transporter C family, cytochrome P450, glucose-methanol-choline oxidoreductase, insect cuticle proteins, and UDP-glucuronosyl), and digestion (e.g., serine proteases and peptidase family S1). A gene ontology enrichment analysis of rapid evolving genes also suggests enriched functions in host sensing and immunity.

CONCLUSIONS

Our results of family evolution analyses indicate that host plant adaptation and pathogen defense could be important drivers in species diversification among gelechiid moths.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Tetraose steroidal glycoalkaloids from potato provide resistance against Alternaria solani and Colorado potato beetle

Plants with innate disease and pest resistance can contribute to more sustainable agriculture. Natural defence compounds produced by plants have the potential to provide a general protective effect against pathogens and pests, but they are not a primary target in resistance breeding. Here, we identified a wild relative of potato, Solanum commersonii, that provides us with unique insight in the role of glycoalkaloids in plant immunity. We cloned two atypical resistance genes that provide resistance to Alternaria solani and Colorado potato beetle through the production of tetraose steroidal glycoalkaloids (SGA). Moreover, we provide in vitro evidence to show that these compounds have potential against a range of different (potato pathogenic) fungi. This research links structural variation in SGAs to resistance against potato diseases and pests. Further research on the biosynthesis of plant defence compounds in different tissues, their toxicity, and the mechanisms for detoxification, can aid the effective use of such compounds to improve sustainability of our food production.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Metabolic consequences of various fruit-based diets in a generalist insect species

Most phytophagous insect species exhibit a limited diet breadth and specialize on a few or a single host plant. In contrast, some species display a remarkably large diet breadth, with host plants spanning several families and many species. It is unclear, however, whether this phylogenetic generalism is supported by a generic metabolic use of common host chemical compounds ('metabolic generalism') or alternatively by distinct uses of diet-specific compounds ('multi-host metabolic specialism')? Here, we simultaneously investigated the metabolomes of fruit diets and of individuals of a generalist phytophagous species, Drosophila suzukii, that developed on them. The direct comparison of metabolomes of diets and consumers enabled us to disentangle the metabolic fate of common and rarer dietary compounds. We showed that the consumption of biochemically dissimilar diets resulted in a canalized, generic response from generalist individuals, consistent with the metabolic generalism hypothesis. We also showed that many diet-specific metabolites, such as those related to the particular color, odor, or taste of diets, were not metabolized, and rather accumulated in consumer individuals, even when probably detrimental to fitness. As a result, while individuals were mostly similar across diets, the detection of their particular diet was straightforward. Our study thus supports the view that dietary generalism may emerge from a passive, opportunistic use of various resources, contrary to more widespread views of an active role of adaptation in this process. Such a passive stance towards dietary chemicals, probably costly in the short term, might favor the later evolution of new diet specializations.

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Metabolization and sequestration of plant specialized metabolites in insect herbivores: Current and emerging approaches

Herbivorous insects encounter diverse plant specialized metabolites (PSMs) in their diet, that have deterrent, anti-nutritional, or toxic properties. Understanding how they cope with PSMs is crucial to understand their biology, population dynamics, and evolution. This review summarizes current and emerging cutting-edge methods that can be used to characterize the metabolic fate of PSMs, from ingestion to excretion or sequestration. It further emphasizes a workflow that enables not only to study PSM metabolism at different scales, but also to tackle and validate the genetic and biochemical mechanisms involved in PSM resistance by herbivores. This review thus aims at facilitating research on PSM-mediated plant-herbivore interactions.

Chromosomal-level reference genome of the moth Heortia vitessoides (Lepidoptera: Crambidae), a major pest of agarwood-producing trees

The moth Heortia vitessoides Moore (Lepidoptera: Crambidae) is a major pest of ecologically, commercially and culturally important agarwood-producing trees in the genus Aquilaria. In particular, H. vitessoides is one of the most destructive defoliating pests of the incense tree Aquilaria sinesis, which produces a valuable fragrant wood used as incense and in traditional Chinese medicine [33]. Nevertheless, a genomic resource for H. vitessoides is lacking. Here, we present a chromosomal-level assembly for H. vitessoides, consisting of a 517 megabase (Mb) genome assembly with high physical contiguity (scaffold N50 of 18.2 Mb) and high completeness (97.9% complete BUSCO score). To aid gene annotation, 8 messenger RNA transcriptomes from different developmental stages were generated, and a total of 16,421 gene models were predicted. Expansion of gene families involved in xenobiotic metabolism and development were detected, including duplications of cytosolic sulfotransferase (SULT) genes shared among lepidopterans. In addition, small RNA sequencing of 5 developmental stages of H. vitessoides facilitated the identification of 85 lepidopteran conserved microRNAs, 94 lineage-specific microRNAs, as well as several microRNA clusters. A large proportion of the H. vitessoides genome consists of repeats, with a 29.12% total genomic contribution from transposable elements, of which long interspersed nuclear elements (LINEs) are the dominant component (17.41%). A sharp decrease in the genome-wide percentage of LINEs with lower levels of genetic distance to family consensus sequences suggests that LINE activity has peaked in H. vitessoides. In contrast, opposing patterns suggest a substantial recent increase in DNA and LTR element activity. Together with annotations of essential sesquiterpenoid hormonal pathways, neuropeptides, microRNAs and transposable elements, the high-quality genomic and transcriptomic resources we provide for the economically important moth H. vitessoides provide a platform for the development of genomic approaches to pest management, and contribute to addressing fundamental research questions in Lepidoptera.

Phylogenomic and functional characterization of an evolutionary conserved cytochrome P450-based insecticide detoxification mechanism in bees

Bee pollinator pesticide risk assessment is a regulatory requirement for pesticide registration and is largely based on experimental data collected for surrogate species such as the western honeybee. Recently, CYP9Q3, a honeybee cytochrome P450 enzyme, has been shown to efficiently detoxify certain insecticides such as the butenolide flupyradifurone and the neonicotinoid thiacloprid. Here we analyzed genomic data for 75 bee species and demonstrated by the recombinant expression of 26 CYP9Q3 putative functional orthologs that this detoxification principle is an evolutionary conserved mechanism across bee families. Our toxicogenomics approach has the potential to inform pesticide risk assessment for nonmanaged bee species that are not accessible for acute toxicity testing.

The regulatory process for assessing the risks of pesticides to bees relies heavily on the use of the honeybee, Apis mellifera, as a model for other bee species. However, the validity of using A. mellifera as a surrogate for other Apis and non-Apis bees in pesticide risk assessment has been questioned. Related to this line of research, recent work on A. mellifera has shown that specific P450 enzymes belonging to the CYP9Q subfamily act as critically important determinants of insecticide sensitivity in this species by efficiently detoxifying certain insecticide chemotypes. However, the extent to which the presence of functional orthologs of these enzymes is conserved across the diversity of bees is unclear. Here we used a phylogenomic approach to identify > 100 putative CYP9Q functional orthologs across 75 bee species encompassing all major bee families. Functional analysis of 26 P450s from 20 representative bee species revealed that P450-mediated detoxification of certain systemic insecticides, including the neonicotinoid thiacloprid and the butenolide flupyradifurone, is conserved across all major bee pollinator families. However, our analyses also reveal that CYP9Q-related genes are not universal to all bee species, with some Megachilidae species lacking such genes. Thus, our results reveal an evolutionary conserved capacity to metabolize certain insecticides across all major bee families while identifying a small number of bee species where this function may have been lost. Furthermore, they illustrate the potential of a toxicogenomic approach to inform pesticide risk assessment for nonmanaged bee species by predicting the capability of bee pollinator species to break down synthetic insecticides.

Caterpillar-Induced Rice Volatile (E)-β-Farnesene Impairs the Development and Survival of Chilo suppressalis Larvae by Disrupting Insect Hormone Balance

Significant research progress has recently been made on establishing the roles of tps46 in rice defense. (E)-β-farnesene (Eβf) is a major product of tps46 activity but its physiological functions and potential mechanisms against Chilo suppressalis have not yet been clarified. In the present study, C. suppressalis larvae were artificially fed a diet containing 0.8 g/kg Eβf and the physiological performance of the larvae was evaluated. In response to Eβf treatment, the average 2nd instar duration significantly increased from 4.78 d to 6.31 d while that of the 3rd instar significantly increased from 5.70 d to 8.00 d compared with the control. There were no significant differences between the control and Eβf-fed 4th and 5th instars in terms of their durations. The mortalities of the 2nd and 3rd Eβf-fed instars were 21.00-fold and 6.39-fold higher, respectively, than that of the control. A comparative transcriptome analysis revealed that multiple differentially expressed genes are involved in insect hormone biosynthesis. An insect hormone assay on the 3rd instars disclosed that Eβf disrupted the balance between the juvenile hormone and ecdysteroid levels. Eβf treatment increased the juvenile hormones titers but not those of the ecdysteroids. The qPCR results were consistent with those of the RNA-Seq. The foregoing findings suggested that Eβf impairs development and survival in C. suppressalis larvae by disrupting their hormone balance. Moreover, Eβf altered the pathways associated with carbohydrate and xenobiotic metabolism as well as those related to cofactors and vitamins in C. suppressalis larvae. The discoveries of this study may contribute to the development and implementation of an integrated control system for C. suppressalis infestations in rice.

Expanding the Menu: Are Polyphagy and Gene Family Expansions Linked across Lepidoptera?

Evolutionary expansions and contractions of gene families are often correlated with key innovations and/or ecological characteristics. In butterflies and moths (Lepidoptera), expansions of gene families involved in detoxification of plant specialized metabolites are hypothesized to facilitate a polyphagous feeding style. However, analyses supporting this hypothesis are mostly based on a limited number of lepidopteran species. We applied a phylogenomics approach, using 37 lepidopteran genomes, to analyze if gene family evolution (gene gain and loss) is associated with the evolution of polyphagy. Specifically, we compared gene counts and evolutionary gene gain and loss rates of gene families involved in adaptations with plant feeding. We correlated gene evolution to host plant family range (phylogenetic diversity) and specialized metabolite content of plant families (functional metabolite diversity). We found a higher rate for gene loss than gene gain in Lepidoptera, a potential consequence of genomic rearrangements and deletions after (potentially small-scale) duplication events. Gene family expansions and contractions varied across lepidopteran families, and were associated to host plant use and specialization levels. Within the family Noctuidae, a higher expansion rate for gene families involved in detoxification can be related to the large number of polyphagous species. However, gene family expansions are observed in both polyphagous and monophagous lepidopteran species and thus seem to be species-specific in the taxa sampled. Nevertheless, a significant positive correlation of gene counts of the carboxyl- and choline esterase and glutathione-S-transferase detoxification gene families with the level of polyphagy was identified across Lepidoptera.

Silencing of cytochrome P450 gene CYP321A1 effects tannin detoxification and metabolism in Spodoptera litura

Cytochrome P450 monooxygenase (P450 or CYP) plays an important role in the metabolism of insecticides and plant allelochemicals by insects. CYP321B1, a novel Spodoptera litura P450 gene, was identified and characterized. CYP321B1 contains a 1488 bp open reading frame (ORF) that encodes a 495 amino acid protein. In fourth instar larvae, the highest CYP321B1 expression levels were found in the midgut and fat body. In the tannin feeding test, tannin can significantly induce the expression of CYP321B1 in the midgut and fat body of 4th instar larvae. To verify the function of CYP321B1, RNA interference and metabolome analysis were performed. The results showed that silencing CYP321B1 significantly reduced the rate of weight gain under tannin induction. Metabolome analysis showed silencing affected 47 different metabolites, mainly involved in secondary metabolite biosynthesis and amino acid metabolism, including amino acids, lipid fatty acids, organic acids and their derivatives. Henoxyacetic acid and cysteamine are the most highly regulated metabolites, respectively. These findings demonstrate that CYP321B1 plays an important role in tannin detoxification and metabolism. Functional knowledge about metabolite detoxification genes in this major herbivorous insect pest can provide new insights into this biological process and provide new targets for agricultural pest control.

The Birth-and-Death Evolution of Cytochrome P450 Genes in Bees

The birth-and-death model of multigene family evolution describes how gene families evolve and diversify through duplication and deletion. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Extensive studies of insect P450 genes have demonstrated their role in insecticide resistance. Bees are thought to experience toxin exposure through their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some species. Here, we describe the repertoire of P450 genes in the orchid bee Euglossa dilemma. Male orchid bees form perfume bouquets used in courtship displays by collecting volatile compounds, resulting in exposure to compounds known to be toxic. In addition, we conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We find that social behavior and resin collection are not correlated with the repertoire of P450 present in a bee species. However, our analyses revealed that P450 clades can be classified as stable and unstable, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures and exhibit signals of adaptive evolution. This work highlights the complexity of multigene family evolution, revealing that multiple factors contribute to the diversification, stability, and dynamics of this gene family. Furthermore, we provide a resource for future detailed studies investigating the function of different P450s in economically important bee species.

Integrated miRNA and transcriptome profiling to explore the molecular determinism of convergent adaptation to corn in two lepidopteran pests of agriculture

Background

The degree to which adaptation to same environment is determined by similar molecular mechanisms, is a topic of broad interest in evolutionary biology, as an indicator of evolutionary predictability. We wished to address if adaptation to the same host plant in phytophagous insects involved related gene expression patterns. We compared sRNA-Seq and RNA-Seq data between two pairs of taxa of Ostrinia and Spodoptera frugiperda sharing maize as host-plant. For the latter, we had previously carried out a reciprocal transplant experiment by feeding of the larvae of the Corn strain (Sf-C) and the Rice strain (Sf-R) on corn versus rice and characterized the mRNA and miRNA responses.

Results

First, we predicted the genes encoding miRNA in Ostrinia nubilalis (On) and O. scapulalis (Os). Respectively 67 and 65 known miRNA genes, as well as 196 and 190 novel ones were predicted with Os genome using sncRNAs extracted from whole larvae feeding on corn or mugwort. In On, a read counts analysis showed that 37 (55.22%) known miRNAs and 19 (9.84%) novel miRNAs were differentially expressed (DE) on mugwort compared to corn (in Os, 25 known miRs (38.46%) and 8 novel ones (4.34%)). Between species on corn, 8 (12.5%) known miRNAs and 8 (6.83%) novel ones were DE while only one novel miRNA showed expression variation between species on mugwort. Gene target prediction led to the identification of 2953 unique target genes in On and 2719 in Os, among which 11.6% (344) were DE when comparing species on corn. 1.8% (54) of On miR targets showed expression variation upon a change of host-plant.

We found molecular changes matching convergent phenotype, i.e., a set of nine miRNAs that are regulated either according to the host-plant both in On and Sf-C or between them on the same plant, corn. Among DE miR target genes between taxa, 13.7% shared exactly the same annotation between the two pairs of taxa and had function related to insect host-plant interaction.

Conclusion

There is some similarity in underlying genetic mechanisms of convergent evolution of two distant Lepidopteran species having adopted corn in their host range, highlighting possible adaptation genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07905-7.

Why are there not more herbivorous insect species?

Insect species richness is estimated to exceed three million species, of which roughly half is herbivorous. Despite the vast number of species and varied life histories, the proportion of herbivorous species among plant-consuming organisms is lower than it could be due to constraints that impose limits to their diversification. These include ecological factors, such as vague interspecific competition; anatomical and physiological limits, such as neural limits and inability of handling a wide range of plant allelochemicals; phylogenetic constraints, like niche conservatism; and most importantly, a low level of concerted genetic variation necessary to a phyletic conversion. It is suggested that diversification ultimately depends on what we call the intrinsic trend of diversification of the insect genome. In support of the above, we survey the major types of host-specificity, the mechanisms and constraints of host specialization, possible pathways of speciation, and hypotheses concerning insect diversification.

CRISPR/Cas9 mutagenesis reveals a role for ABCB1 in gut immune responses to Vibrio diazotrophicus in sea urchin larvae

The ABC transporter ABCB1 plays an important role in the disposition of xenobiotics. Embryos of most species express high levels of this transporter in early development as a protective mechanism, but its native substrates are not known. Here, we used larvae of the sea urchin Strongylocentrotus purpuratus to characterize the early life expression and role of Sp-ABCB1a, a homolog of ABCB1. The results indicate that while Sp-ABCB1a is initially expressed ubiquitously, it becomes enriched in the developing gut. Using optimized CRISPR/Cas9 gene editing methods to achieve high editing efficiency in the F0 generation, we generated ABCB1a crispant embryos with significantly reduced transporter efflux activity. When infected with the opportunistic pathogen Vibrio diazotrophicus, Sp-ABCB1a crispant larvae demonstrated significantly stronger gut inflammation, immunocyte migration and cytokine Sp-IL-17 induction, as compared with infected control larvae. The results suggest an ancestral function of ABCB1 in host-microbial interactions, with implications for the survival of invertebrate larvae in the marine microbial environment.

Signatures of selection and evolutionary relevance of cytochrome P450s in plant-insect interactions

Enzymes in the cytochrome P450 (P450) superfamily have important functions ranging from those that are essential for the physiology and development of the individual to those that mediate interactions between individuals and their biotic environment. Until recently the study of P450s had focused on single functions, substrates, or pathways. Recent advances in sequencing, genome assembly, and phylogenetic methods have returned emphasis to the adaptive value of these enzymes in the context of herbivory. Comparisons of whole repertoires of P450s across related species reveal that P450s capable of metabolizing xenobiotics have an increased rate of gains compared to losses after gene duplications. In plants, studies have focused on enzymes and end-functions that have converged to provide increased resistance to herbivory. This review summarizes the latest findings related to the ecological value of P450s in the interactions between phytophagous insects and their host plants.

Cytochrome P450-mediated mycotoxin metabolism by plant-feeding insects

Mycotoxins are secondary metabolites produced primarily by filamentous fungi that when consumed cause pathological responses in animal hosts or consumers. Defined functionally rather than structurally, mycotoxins derive from numerous primary metabolic pathways. Through opportunistic or mutualistic associations, insect herbivores inflict damage that can predispose plants to infection by mycotoxin-producing phytopathogens, resulting in economically significant contamination. The few cytochrome P450 subfamilies implicated in mycotoxin detoxification by insects, including CYP6 and CYP9, are also known to detoxify phytochemicals. Some insect P450s bioactivate, rather than detoxify, mycotoxins, suggestive of an 'escalation' in arms-race interactions between these herbivores and fungi. Characterizing insect P450s that detoxify mycotoxins can be useful for developing biological remediation technologies and for ensuring the safety of insects reared for human or livestock consumption.

Differential regulation of cytochrome P450 genes associated with biosynthesis and detoxification in bifenthrin-resistant populations of navel orangewom (Amyelois transitella)

Pyrethroid resistance was first reported in 2013 for the navel orangeworm, Amyelois transitella, but the genetic underpinnings of pyrethroid resistance are unknown. We investigated the role of cytochrome P450 monooxygenases (P450s) belonging to the CYP3 and CYP4 clans using colonies derived from individuals collected in 2016 from almond orchards in two counties. One colony (ALM) originated from an almond orchard in Madera County with no reported pyrethroid resistance and the second colony (R347) originated from the same Kern County orchard where pyrethroid resistance was first reported. We used high-throughput quantitative real-time PCR (qRT-PCR) analyses of 65 P450s in the CYP3 and CYP4 clans of A. transitella to identify P450s induced by bifenthrin and associated with pyrethroid resistance. Nine P450s were constitutively overexpressed in R347 compared to ALM, including CYP6AE54 (11.7-fold), belonging to a subfamily associated with metabolic pesticide detoxification in Lepidoptera and CYP4G89 (33-fold) belonging to a subfamily associated with cuticular hydrocarbon (CHC) synthesis and resistance via reduced pesticide penetrance. Cuticular hydrocarbons analysis revealed that R347 produced twice as many total CHCs in the egg and adult stages as ALM. Topical toxicity bioassays for R347 determined that egg mortality was reduced at low bifenthrin concentrations and larval mortality was reduced at high concentrations of bifenthrin compared to ALM. Our discovery of both changes in metabolism and production of CHCs for R347 have implications for the possible decreased efficacy of other classes of insecticide used to control this insect. The threat of widespread pyrethroid resistance combined with the potential for cross-resistance to develop through the mechanism of reduced penetrance warrants developing management strategies that facilitate insecticide passage across the cuticle.

Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

The mega-diversity of herbivorous insects is attributed to their co-evolutionary associations with plants. Despite abundant studies on insect-plant interactions, we do not know whether host-plant shifts have impacted both genomic adaptation and species diversification over geological times. We show that the antagonistic insect-plant interaction between swallowtail butterflies and the highly toxic birthworts began 55 million years ago in Beringia, followed by several major ancient host-plant shifts. This evolutionary framework provides a valuable opportunity for repeated tests of genomic signatures of macroevolutionary changes and estimation of diversification rates across their phylogeny. We find that host-plant shifts in butterflies are associated with both genome-wide adaptive molecular evolution (more genes under positive selection) and repeated bursts of speciation rates, contributing to an increase in global diversification through time. Our study links ecological changes, genome-wide adaptations and macroevolutionary consequences, lending support to the importance of ecological interactions as evolutionary drivers over long time periods.

Roles of the variable P450 substrate recognition sites SRS1 and SRS6 in esfenvalerate metabolism by CYP6AE subfamily enzymes in Helicoverpa armigera

The cotton bollworm P450s of the clustered CYP6AE subfamily share high sequence identities but differ dramatically in their capacity to metabolize xenobiotics, especially esfenvalerate. Among them, CYP6AE17 has the highest sequence identity with CYP6AE18 but shows ~7-fold higher metabolic efficiency. CYP6AE11 is most active towards esfenvalerate but CYP6AE20 is inactive even though the enzymes share 54.8% sequence identity. Sequence analysis revealed the SRS1 (Substrate Recognition Site) and SRS6 between CYP6AE17 and CYP6AE18, and SRS1 between CYP6AE11 and CYP6AE20 are the most variable among all six SRSs. In order to identify the key factors that underlie the observed catalytic difference, we exchanged these SRS sequences between two pairs of P450s and studied the activity of the resulting hybrid mutants or chimeras. In vitro metabolism showed that the CYP6AE17/18 chimeras had 2- and 14-fold decreased activities and the CYP6AE18/17 chimeras had 6- and 10-fold increased activities to esfenvalerate. Meanwhile, after exchanging SRS1 with each other, the CYP6AE11/20 chimera folded incorrectly but the CYP6AE20/11 chimera gained moderate activity to esfenvalerate. Molecular modelling showed that amino acids variants within SRS1 or SRS6 change the shape and chemical environment of the active sites, which may affect the ligand-binding interactions. These results indicate that the protein structure variation resulting from the sequence diversity of SRSs promotes the evolution of insect chemical defense and contributes to the development of insect resistance to pesticides.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

Selective Sweeps in a Nutshell: The Genomic Footprint of Rapid Insecticide Resistance Evolution in the Almond Agroecosystem

Among the most familiar forms of human-driven evolution on ecological time scales is the rapid acquisition of resistance to pesticides by insects. Since the widespread adoption of synthetic organic insecticides in the mid-twentieth century, over 500 arthropod species have evolved resistance to at least one insecticide. Efforts to determine the genetic bases of insecticide resistance have historically focused on individual loci, but the availability of genomic tools has facilitated the screening of genome-wide characteristics. We resequenced three contemporary populations of the navel orangeworm (Amyelois transitella), the principal pest of almond orchards in California, differing in bifenthrin resistance status to examine insecticide-induced changes in the population genomic landscape of this species. We detected an exceptionally large region with virtually no polymorphisms, extending to up to 1.3 Mb in the resistant population. This selective sweep includes genes associated with pyrethroid and DDT resistance, including a cytochrome P450 gene cluster and the gene encoding the voltage-gated sodium channel para. Moreover, the sequence along the sweep is nearly identical in the genome assembled from a population founded in 1966, suggesting that the foundation for insecticide resistance may date back a half-century, when California’s Central Valley experienced massive area-wide applications of DDT for pest control.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Substrate-specificity of cytochrome P450-mediated detoxification as an evolutionary strategy for specialization on furanocoumarin-containing hostplants: CYP6AE89 in parsnip webworms

The parsnip webworm, Depressaria pastinacella, is restricted to two hostplant genera containing six structurally diverse furanocoumarins. Of these, imperatorin is detoxified by a specialized cytochrome P450, CYP6AB3. A previous whole-larva transcriptome analysis confirmed the presence of nine transcripts that belong to the CYP6AE subfamily. Here, by examining midgut-specific gene expression patterns we determined that CYP6AE89 transcripts were highly expressed and furanocoumarin-inducible. Computer docking and energy-minimization of a CYP6AE89 model with all six furanocoumarins showed that 5-methoxylated bergapten and 8-methoxylated xanthotoxin had the smallest distances from the heme to the proton-donor residue in the catalytic I-helix, and that the 5,8-dimethoxylated isopimpinellin and bergapten had the smallest energy-minimized distance from the heme oxygen to the furan ring double bond. To evaluate this prediction, we expressed the CYP6AE89 protein in an Escherichia coli system, and used it to detect high catalytic activity against the two mono-methoxylated linear furanocoumarins - bergapten and xanthotoxin - and weak activity against isopimpinellin. Thus, CYP6AE89, like CYP6AB3, is probably specialized for detoxifying only a subset of hostplant furanocoumarins. A maximum-likelihood tree built with six representative lepidopterans with manually annotated cytochrome P450s shows that CYP6AE89 may have evolved much faster than the other CYP6AE proteins, possibly indicative of host selection pressure.

Genomic and transcriptomic analyses of glutathione S-transferases in an endoparasitoid wasp, Pteromalus puparum

Pteromalus puparum is a gregarious pupal endoparasitoid with a wide host range. It deposits eggs into pierid and papilionid butterfly pupae. Glutathione S-transferases (GSTs) are a family of multifunctional detoxification enzymes that act in xenobiotic metabolism in insects. Insect genome projects have facilitated identification and characterization of GST family members. We identified 20 putative GSTs in the P. puparum genome, including 19 cytosolic and one microsomal. Phylogenetic analysis showed that P. puparum GSTs are clustered into Hymenoptera-specific branches. Transcriptomic data of embryos, larvae, female pupae, male pupae, female adults, male adults, venom glands, carcass, salivary glands, and ovaries revealed stage-, sex-, and tissue-specific expression patterns of GSTs in P. puparum. This is the most comprehensive study of genome-wide identification, characterization, and expression profiling of GST family in hymenopterans. Our results provide valuable information for understanding the metabolic adaptation of this wasp.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Molecular Insights into the Insensitivity of Lepidopteran Pests to Cycloxaprid

Cycloxaprid (CYC) is effective in the control of hemipteran pests, but its bioactivity against lepidopteran pests is still unclear. Here, the bioactivity of CYC against lepidopteran pests was found to be much worse than that against hemipteran insects. To reveal the mechanism, the transcriptomes of CYC-treated and untreated Ostrinia furnacalis larvae were compared. Among the top 20 differentially expressed genes, 11 encode proteins involved in cuticle formation, while only one encodes a detoxifying enzyme. Thus, the cuticle appears to be important for the insensitivity of O. furnacalis to CYC. A pretreatment of O. furnacalis larvae with methoprene enhanced the bioactivity of CYC by 1.12-fold. Moreover, mixtures of CYC with graphene oxide increased the bioactivity of CYC by 1.88-fold. Because lepidopteran and hemipteran insects often harm crops at the same time, the work can help make full use of CYC and reduce the environmental impacts of using multiple pesticides.

How interactions with plant chemicals shape insect genomes

The transition to herbivory by insects is associated with distinct genomic signatures. Sequenced genomes of extant herbivore species reveal the result of these transitions, but in lieu of comparisons between herbivorous and non-herbivorous lineages that diverged recently, such datasets have shed less light on the evolutionary genomic processes involved in diet shifts to or from herbivory. Here, we propose that the comparative genomics of diet shifts between closely related insect herbivores and non-herbivores, and within densely-sampled clades of herbivores, will help reveal the extent to which herbivory evolves through the co-option and subtle remodeling of widely-conserved gene families with functions ancestrally distinct from phytophagy.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Transcriptional Responses of the Trichoplusia ni Midgut to Oral Infection by the Baculovirus Autographa californica Multiple Nucleopolyhedrovirus

Baculoviruses are large double-stranded DNA viruses that are virulent pathogens of certain insect species. In a natural host, Trichoplusia ni, infection by the model baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) begins when the occluded form of the virus disassembles in the midgut and virions infect midgut epithelial cells to establish the primary phase of the infection. To better understand the primary phase of the AcMNPV infection cycle, newly molted 5th-instar T. ni larvae were orally infected with AcMNPV occlusion bodies and the transcriptional responses of the T. ni midgut were analyzed at various times from 0 to 72 h postinfection, using transcriptome sequencing analysis and a T. ni reference genome. The numbers of differentially expressed host genes increased as the infection progressed, and we identified a total of 3,372 differentially expressed T. ni transcripts in the AcMNPV-infected midgut. Genes encoding orthologs of HMG176, atlastin, and CPH43 were among the most dramatically upregulated in response to AcMNPV infection. A number of cytochrome P450 genes were downregulated in response to infection. We also identified the effects of AcMNPV infection on a large variety of genes associated with innate immunity. This analysis provides an abundance of new and detailed information on host responses to baculovirus infection during the primary phase of the infection in the midgut and will be important for understanding how baculoviruses establish productive infections in the organism.IMPORTANCE Baculoviruses are virulent pathogens of a number of important insect pest species. In the host Trichoplusia ni, infection begins in the midgut when infectious virions of the occlusion-derived virus (ODV) phenotype enter and subsequently replicate in cells of the midgut epithelium. A second virion phenotype (budded virus [BV]) is produced there, and BV mediates systemic infection of the animal. Most prior detailed studies of baculovirus infections have focused on BV infections of cultured cells. In this study, we examined the transcriptional responses of the T. ni midgut to infection by ODV of the baculovirus AcMNPV and identified a variety of host genes that respond dramatically to viral infection. Understanding the transcriptional responses of the host midgut to viral infection is critically important for understanding the biphasic infection in the animal as a whole.

CYPs in different families are involved in the divergent regio-specific epoxidation of alkenyl sex pheromone precursors in moths

In moth species that utilize alkenyl sex pheromones, the epoxidation of alkenes confers further diversity on the chemical structures of pheromone components. Hc_epo1 (CYP341B14), the first pheromone gland (PG)-specific epoxidase identified from the fall webworm Hyphantria cunea (Erebidae), specifically epoxidizes the Z9 double bond in the triene precursor, (3Z,6Z,9Z)-3,6,9-henicosatriene (Z3,Z6,Z9-21:H). In the present study, we identified a novel PG-specific epoxidase, As_epo1, from the Japanese giant looper Ascotis selenaria (Geometridae), which secretes cis-3,4-epoxy-(6Z,9Z)-6,9-nonadecadiene (epo3,Z6,Z9-19:H) as the main sex pheromone component. A functional assay using the Sf9 insect cell line-baculovirus expression system showed that As_epo1 specifically epoxidizes the Z3 double bond in the pheromone precursor triene, (3Z,6Z,9Z)-3,6,9-nonadecatriene (Z3,Z6,Z9-19:H). As_epo1 also Z3-specifically epoxidized a triene with a longer carbon chain, Z3,Z6,Z9-21:H, which does not occur in this species. A phylogenetic analysis indicated that As_epo1 belonged to the CYP340 family, not the CYP341 family to which Hc_epo1 belongs. These results suggest that moth PG-specific epoxidases with divergent regio-specificities have evolved independently.

Functional characterization of the epoxidase gene, Li_epo1 (CYP341B14), involved in generation of epoxyalkene pheromones in the mulberry tiger moth Lemyra imparilis

Epoxidation of alkenes derived from essential fatty acids is a key step in the biosynthesis of sex pheromones in moth species that utilize alkenyl sex pheromones. The position of the epoxy ring in the pheromone molecule differs depending on the species, thereby conferring diversities on sex pheromones. To date, only one pheromone gland (PG)-specific epoxidase, Hc_epo1 (CYP341B14), has been reported. Hc_epo1, which was identified from an arctiid moth Hyphantria cunea, catalyzes the epoxidation of a double bond at position 9 of the triene, Z3,Z6,Z9-21:H. In the present study, we investigated the PG-specific epoxidase from another arctiid, the mulberry tiger moth Lemyra imparilis, in order to verify whether cytochrome P450 in the CYP341B subfamily, to which Hc_epo1 belongs to, is responsible for the epoxidation of pheromone precursors at position 9 in moths other than H. cunea. A fragment of the Hc_epo1 homolog was amplified from cDNA prepared from the PG of L. imparilis by PCR with degenerate primers. The deduced amino acid sequence of the subsequently cloned homolog, Li_epo1, showed 88.5% identity to Hc_epo1. A functional assay using the Sf9 insect cell line-baculovirus expression system showed that Li_epo1 exhibited epoxidase activity with high selectivity to the double bond at position 9 of two trienes, Z3,Z6,Z9-21:H and Z3,Z6,Z9-23:H, precursors of epoxy diene sex pheromone components in L. imparilis.

Detoxifying enzyme complements and host use phenotypes in 160 insect species

We use the genomes of 160 insect species to test the hypothesis that the size of detoxifying enzyme families is greater in species using more chemically diverse food resources. Phylogenetically appropriate contrasts in subsamples of the data generally support the hypothesis. We find relatively high numbers of cytochrome P450, glutathione S-transferase and carboxyl/choline esterase genes in omnivores and herbivores feeding on chemically complex tissues and relatively low numbers of these genes in specialists on relatively simple diets, including plant sap, nectar and pollen, and blood. Among Lepidoptera feeding on green plant tissue and Condylognatha feeding on sap we also find more of these genes in highly polyphagous species, many of which are major agricultural pests. These genomic signatures of food resource use are consistent with the hypothesis that some taxa are preadapted for insecticide resistance evolution.

Genome Size Reversely Correlates With Host Plant Range in Helicoverpa Species

In organisms with very low percentages of transposable elements (TEs), genome size may positively or negatively correlate with host range, depending on whether host adaptation or host modification is the main route to host generalism. To test if this holds true for insect herbivores with greater percentages of TEs, we conducted flow cytometry to measure the endopolyploidy levels and C-values of the host modification (salivary gland and mandibular gland in head), host adaptation (midgut), and host use-independent tissues (male gonad, hemolymph, and Malpighian tubules) of the generalist Helicoverpa armigera and the head of its older specialist sister H. assulta. Larval salivary gland displayed a consecutive chain of endopolyploidy particles from 8Cx to higher than 32Cx and larval head and midgut had endopolyploidy nuclei clusters of 16Cx and 32Cx, whereas larval male gonad, hemolymph, and Malpighian tubules possessed no endopolyploidy nuclei of higher than 8Cx. The estimated genome size of the Solanaceae plant specialist H. assulta is 430 Mb, significantly larger than that of its older generalist sister Heliothis virescens (408 Mb) and those of its two generalist descendants H. armigera (394 Mb) and H. zea (363 Mb). These data not only reveal a negative correlation between host plant range and genome size in this terminal lineage, but also imply that Helicoverpa species appear to depend more on host modification than on host adaptation to achieve polyphagy.

Molecular Evolution of the Major Arthropod Chemoreceptor Gene Families

The evolutionary origins of the three major families of chemoreceptors in arthropods-the odorant receptor (OR), gustatory receptor (GR), and ionotropic receptor (IR) families-occurred at the base of the Insecta, Animalia, and Protostomia, respectively. Comparison of receptor family sizes across arthropods reveals a generally positive correlation with their widely disparate complexity of chemical ecology. Closely related species reveal the ongoing processes of gene family evolution, including gene duplication, divergence, pseudogenization, and loss, that mediate these larger patterns. Sets of paralogous receptors within species reveal positive selection on amino acids in regions likely to contribute to ligand binding and specificity. Ligands of many ORs and some GRs and IRs have been identified; however, ligand identification for many more chemoreceptors is needed, as are structures for the OR/GR superfamily, to improve our understanding of the molecular evolution of these ecologically important receptors in arthropods.

Preferences of Specialist and Generalist Mammalian Herbivores for Mixtures Versus Individual Plant Secondary Metabolites

Herbivores that forage on chemically defended plants consume complex mixtures of plant secondary metabolites (PSMs). However, the mechanisms by which herbivores tolerate mixtures of PSMs are relatively poorly understood. As such, it remains difficult to predict how PSMs, singly or as complex mixtures, influence diet selection by herbivores. Although relative rates of detoxification of PSMs have been used to explain tolerance of PSMs by dietary specialist herbivores, few studies have used the rate of detoxification of individual PSMs to understand dietary preferences of individual herbivores for individual versus mixtures of PSMs. We coupled in vivo experiments using captive feeding trials with in vitro experiments using enzymatic detoxification assays to evaluate the dietary preferences and detoxification capacities of pygmy rabbits (Brachylagus idahoensis), dietary specialists on sagebrush (Artemisia spp.), and mountain cottontails (Sylvilagus nuttallii), dietary generalists. We compared preference for five single PSMs in sagebrush compared to a mixture containing those same five PSMs. We hypothesized that relative preference for individual PSMs would coincide with faster detoxification capacity for those PSMs by specialists and generalists. Pygmy rabbits generally showed little preference among individual PSMs compared to mixed PSMs, whereas mountain cottontails exhibited stronger preferences. Pygmy rabbits had faster detoxification capacities for all PSMs and consumed higher concentrations of individual PSMs versus a mixture than cottontails. However, detoxification capacity for an individual PSM did not generally coincide with preferences or avoidance of individual PSMs by either species. Cottontails avoided, but pygmy rabbits preferred, camphor, the PSM with the slowest detoxification rate by both species. Both species avoided β-pinene despite it having one of the fastest detoxification rate. Taken together our in vivo and in vitro results add to existing evidence that detoxification capacity is higher in dietary specialist than generalist herbivores. However, results also suggest that alternative mechanisms such as absorption and the pharmacological action of individual or mixtures of PSMs may play a role in determining preference of PSMs within herbivore species.

miRNA-Mediated Interactions in and between Plants and Insects

Our understanding of microRNA (miRNA) regulation of gene expression and protein translation, as a critical area of cellular regulation, has blossomed in the last two decades. Recently, it has become apparent that in plant-insect interactions, both plants and insects use miRNAs to regulate their biological processes, as well as co-opting each others' miRNA systems. In this review article, we discuss the current paradigms of miRNA-mediated cellular regulation and provide examples of plant-insect interactions that utilize this regulation. Lastly, we discuss the potential biotechnological applications of utilizing miRNAs in agriculture.

Does host plant adaptation lead to pesticide resistance in generalist herbivores?

Most herbivorous arthropods feed on one or a few closely related plant species; however, certain insect and mite species have a greatly expanded host range. Several of these generalists also show a remarkable propensity to evolve resistance to chemical pesticides. In this review, we ask if the evolution of mechanisms to tolerate the diversity of plant secondary metabolites that generalist herbivores encounter, has pre-adapted them to resist synthetic pesticides. Critical examination of the evidence suggests that a generalist life-style per se is not a predictor of rapid resistance evolution to pesticides. Rather the prevalence of pesticide resistance in generalist herbivores probably reflects their economic importance as pests and thus the strong selection imposed by intensive pesticide use.