Pyrosequencing the Manduca sexta larval midgut transcriptome: messages for digestion, detoxification and defence

Starvation-Induced Changes to the Midgut Proteome and Neuropeptides in Manduca sexta

Starvation is a complex physiological state that induces changes in protein expression to ensure survival. The insect midgut is sensitive to changes in dietary content as it is at the forefront of communicating information about incoming nutrients to the body via hormones. Therefore, a DIA proteomics approach was used to examine starvation physiology and, specifically, the role of midgut neuropeptide hormones in a representative lepidopteran, Manduca sexta. Proteomes were generated from midguts of M. sexta fourth-instar caterpillars, starved for 24 h and 48 h, and compared to fed controls. A total of 3047 proteins were identified, and 854 of these were significantly different in abundance. KEGG analysis revealed that metabolism pathways were less abundant in starved caterpillars, but oxidative phosphorylation proteins were more abundant. In addition, six neuropeptides or related signaling cascade proteins were detected. Particularly, neuropeptide F1 (NPF1) was significantly higher in abundance in starved larvae. A change in juvenile hormone-degrading enzymes was also detected during starvation. Overall, our results provide an exploration of the midgut response to starvation in M. sexta and validate DIA proteomics as a useful tool for quantifying insect midgut neuropeptide hormones.

Differential expression of immunity-related genes in larval Manduca sexta tissues in response to gut and systemic infection

Introduction

The midgut epithelium functions as tissue for nutrient uptake as well as physical barrier against pathogens. Additionally, it responds to pathogen contact by production and release of various factors including antimicrobial peptides, similar to the systemic innate immune response. However, if such a response is restricted to a local stimulus or if it appears in response to a systemic infection, too is a rather underexplored topic in insect immunity. We addressed the role of the midgut and the role of systemic immune tissues in the defense against gut-borne and systemic infections, respectively.

Methods

Manduca sexta larvae were challenged with DAP-type peptidoglycan bacteria - Bacillus thuringiensis for local gut infection and Escherichia coli for systemic stimulation. We compared the immune response to both infection models by measuring mRNA levels of four selected immunity-related genes in midgut, fat body, hematopoietic organs (HOs), and hemocytes, and determined hemolymph antimicrobial activity. Hemocytes and HOs were tested for presence and distribution of lysozyme mRNA and protein.

Results

The midgut and circulating hemocytes exhibited a significantly increased level of lysozyme mRNA in response to gut infection but did not significantly alter expression in response to a systemic infection. Conversely, fat body and HOs responded to both infection models by altered mRNA levels of at least one gene monitored. Most, but not all hemocytes and HO cells contain lysozyme mRNA and protein.

Discussion

These data suggest that the gut recruits immune-related tissues in response to gut infection whereas systemic infections do not induce a response in the midgut. The experimental approach implies a skewed cross-talk: An intestinal infection triggers immune activity in systemic immune organs, while a systemic infection does not elicit any or only a restricted immune response in the midgut. The HOs, which form and release hemocytes in larval M. sexta, i) synthesize lysozyme, and ii) respond to immune challenges by increased immune gene expression. These findings strongly suggest that they not only provide phagocytes for the cellular immune response but also synthesize humoral immune components.

BmTBP upregulates the transcription of BmSuc1 in silkworm (Bombyx mori) by binding to BmTfΙΙA-S

The BmSuc1 gene, which encodes a novel animal-type β-fructofuranosidase (EC 3.2.1.26), was first cloned and identified in silkworm (Bombyx mori). As an essential sucrase, the activity of BmSUC1 is unaffected by alkaloidal sugar mimics in mulberry leaves. This enzyme may also directly regulate the degree of sucrose hydrolysis in the silkworm midgut. In addition, BmSUC1 is involved in the synthesis of sericin 1 in the silk gland tissue. However, the mechanism underlying the regulation of BmSuc1 transcription remains unclear. In this study, we analyzed the BmSuc1 promoter activity using a dual-luciferase reporter assay and identified 4 regions that are critical for transcriptional activation. The gene encoding a predicted transcription factor (TATA-box-binding protein; BmTBP) capable of binding to the core promoter regions was cloned. A quantitative real-time polymerase chain reaction analysis indicated the gene was highly expressed in the midgut. Downregulating BmTBP expression via RNA interference decreased the expression of BmSuc1 at the transcript and protein levels. An electrophoretic mobility shift analysis and chromatin immunoprecipitation indicated that BmTBP can bind to the TATA-box cis-regulatory element in the BmSuc1 promoter. Furthermore, a bioinformatics-based analysis and a far-western blot revealed the interaction between BmTBP and another transcription factor (BmTfIIA-S). The luciferase reporter gene assay results confirmed that the BmTBP-BmTfIIA-S complex increases the BmSuc1 promoter activity. Considered together, these findings suggest that BmTBP regulates BmSuc1 expression through its interaction with BmTfIIA-S.

A salivary GMC oxidoreductase of Manduca sexta re-arranges the green leaf volatile profile of its host plant

Green leaf volatiles (GLVs) are short-chain oxylipins that are emitted from plants in response to stress. Previous studies have shown that oral secretions (OS) of the tobacco hornworm Manduca sexta, introduced into plant wounds during feeding, catalyze the re-arrangement of GLVs from Z-3- to E-2-isomers. This change in the volatile signal however is bittersweet for the insect as it can be used by their natural enemies, as a prey location cue. Here we show that (3Z):(2E)-hexenal isomerase (Hi-1) in M. sexta's OS catalyzes the conversion of the GLV Z-3-hexenal to E-2-hexenal. Hi-1 mutants that were raised on a GLV-free diet showed developmental disorders, indicating that Hi-1 also metabolizes other substrates important for the insect's development. Phylogenetic analysis placed Hi-1 within the GMCβ-subfamily and showed that Hi-1 homologs from other lepidopterans could catalyze similar reactions. Our results indicate that Hi-1 not only modulates the plant's GLV-bouquet but also functions in insect development.

BmSuc1 Affects Silk Properties by Acting on Sericin1 in Bombyx mori

BmSuc1, a novel animal-type β-fructofuranosidase (β-FFase, EC 3.2.1.26) encoding gene, was cloned and identified for the first time in the silkworm, Bombyx mori. BmSuc1 was specifically and highly expressed in the midgut and silk gland of Bombyx mori. Until now, the function of BmSuc1 in the silk gland was unclear. In this study, it was found that the expression changes of BmSuc1 in the fifth instar silk gland were consistent with the growth rate of the silk gland. Next, with the aid of the CRISPR/Cas9 system, the BmSuc1 locus was genetically mutated, and homozygous mutant silkworm strains with truncated β-FFase (BmSUC1) proteins were established. BmSuc1 mutant larvae exhibited stunted growth and decreased body weight. Interestingly, the molecular weight of part of Sericin1 (Ser1) in the silk gland of the mutant silkworms was reduced. The knockout of BmSuc1 reduced the sericin content in the silkworm cocoon shell, and the mechanical properties of the mutant line silk fibers were also negatively affected. These results reveal that BmSUC1 is involved in the synthesis of Ser1 protein in silk glands and helps to maintain the homeostasis of silk protein content in silk fibers and the mechanical properties of silk fibers, laying a foundation for the study of BmSUC1 regulation of silk protein synthesis in silk glands.

Quantity versus quality: Effects of diet protein-carbohydrate ratios and amounts on insect herbivore gene expression

Dietary protein and digestible carbohydrates are two key macronutrients for insect herbivores, but the amounts and ratios of these two macronutrients in plant vegetative tissues can be highly variable. Typically, insect herbivores regulate their protein-carbohydrate intake by feeding selectively on nutritionally complementary plant tissues, but this may not always be possible. Interestingly, lab experiments consistently demonstrate that performance - especially growth and survival - does not vary greatly when caterpillars and nymphal grasshoppers are reared on diets that differ in their protein-carbohydrate content. This suggests insect herbivores employ post-ingestive physiological mechanisms to compensate for variation in diet protein-carbohydrate profile. However, the molecular mechanisms that underlie this compensation are not well understood. Here we explore, for the first time in an insect herbivore, the transcriptional effects of two dietary factors: protein-to-carbohydrate ratio (p:c) and total macronutrient (p + c) content. Specifically, we reared Helicoverpa zea caterpillars on three diets that varied in diet p:c ratio and one diet that varied in total p + c concentration, all within an ecologically-relevant range. We observed two key findings. Caterpillars reared on diets with elevated total p + c content showed large differences in gene expression. In contrast, only small differences in gene expression were observed when caterpillars were reared on diets with different p:c ratios (spanning from protein-biased to carbohydrate-biased). The invariable expression of many metabolic genes across these variable diets suggests that H. zea caterpillars employ a strategy of constitutive expression to deal with protein-carbohydrate imbalances rather than diet-specific changes. This is further supported by two findings. First, few genes were uniquely associated with feeding on a protein- and carbohydrate-biased diet. Second, many differentially-expressed genes were shared across protein-biased, carbohydrate-biased, and concentrated diet treatments. Our study provides insights into the post-ingestive physiological mechanisms insect herbivores employ to regulate protein-carbohydrate intake. Most notably, it suggests that H. zea, and perhaps other generalist species, use similar post-ingestive mechanisms to deal with protein-carbohydrate imbalances - regardless of the direction of the imbalance.

Identification and Characterization of CYP6 Family Genes from the Oriental Fruit Moth (Grapholita molesta) and Their Responses to Insecticides

Simple Summary

Eight CYP6 family genes were identified in Grapholita molesta (Busck). The expression of individual gene members differed between the developmental stages and insect tissues. High expression was found in third/fourth-instar larvae and in the midgut and Malpighian tubules. The response patterns of the genes exhibited diverse response patterns to the three representative insecticides were diverse.

Abstract

Cytochrome P450 (CYP) monooxygenases comprise a superfamily of proteins that detoxify xenobiotics and plant secondary metabolites in insects. The CYP6 family is unique to the class Insecta, and its members participate in the metabolism of exogenous substances. In this study, we sequenced and characterized the full-length cDNAs of eight CYP6 family genes from Grapholita molesta (Busck), a global pest of pome fruits. P450 genes with the exception of CYP6AN35, which was most highly expressed in adults, consistently showed high expression in third- or fourth-instar larvae. The analysis of different tissues of adults showed that most of these genes were predominantly expressed in the midgut, Malpighian tubules, and/or fat body. The expression of these eight CYP6 genes was differentially affected by three representative insecticides: malathion (organophosphate), deltamethrin (pyrethroid), and chlorantraniliprole (carbamate). All eight CYP6 genes responded to malathion treatment. Only three CYP6 genes were highly expressed in deltamethrin-treated individuals. Chlorantraniliprole treatment exerted weak effects on gene expression. Interestingly, CYP6AN35 was a highly expression level in the adult head and its expression was induced by all three insecticides. CYP6AN35 may be a key gene in the metabolism of insecticides. This study provides a fundamental understanding of the functions of the CYP6 gene family in insecticide metabolism in G. molesta.

Transcriptomic analysis of salivary gland and proteomic analysis of oral secretion in Helicoverpa armigera under cotton plant leaves, gossypol, and tannin stresses

Gossypol and tannin are involved in important chemical defense processes in cotton plants. In this study, we used transcriptomics and proteomics to explore the changes in salivary gland functional genes and oral secretion (OS) proteins after feeding with artificial diet (containing gossypols and tannins) and cotton plant leaves. We found that dietary cotton plant leaves, gossypols and tannins exerted adverse impacts on the genes that regulated the functions of peptidase, GTPase, glycosyl hydrolases in the salivary glands of the Helicoverpa armigera (H. armigera). However, GST, UGT, hydrolases, and lipase genes were up-regulated to participate in the detoxification and digestive of H. armigera. The oral secretory proteins of H. armigera were significantly inhibited under the stress of gossypol and tannin, such as enzyme activity, but some proteins (such as PZC71358.1) were up-regulated and involved in immune and digestive functions. The combined analysis of transcriptomics and metabolomics showed a weak correlation, and the genes and proteins involved were mainly in digestive enzyme activities. Our work clarifies the deleterious physiological impacts of gossypols and tannins on H. armigera and reveals the mechanism by which H. armigera effectively mitigate the phytotoxic effects through detoxification and immune systems.

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.

Bacillus thuringiensis and Chlorantraniliprole Trigger the Expression of Detoxification-Related Genes in the Larval Midgut of Plutella xylostella

Background: Diamondback moth (DBM), Plutella xylostella (L.), has developed resistance to many insecticides. The molecular mechanism of DBM resistance to Bt-G033A combined with chlorantraniliprole (CL) remains undefined. Methods: In this study, field-resistant strains of Plutella xylostella to three pesticides, namely, Bacillus thuringiensis (Bt) toxin (Bt-G033A), CL, and a mixture of Bt + CL, were selected to evaluate the resistance level. Additionally, transcriptomic profiles of a susceptible (SS-DBM), field-resistant (FOH-DBM), Bt-resistant (Bt-DBM), CL-resistant (CL-DBM), and Bt + CL-resistant (BtC-DBM) strains were performed by comparative analysis to identify genes responsible for detoxification. Results: The Bt-G033A was the most toxic chemical to all the DBM strains among the three insecticides. The comparative analysis identified 25,518 differentially expressed genes (DEGs) between pairs/combinations of strains. DEGs were enriched in pathways related to metabolic and catalytic activity and ABC transporter in resistant strains. In total, 17 metabolic resistance-related candidate genes were identified in resistance to Bt-G033A, CL, and Bt + CL by co-expression network analysis. Within candidate genes, the majority was upregulated in key genes including cytochrome P450, glutathione S-transferase (GST), carboxylesterase, and acetylcholinesterase in CL- and BtC-resistant strains. Furthermore, aminopeptidase N (APN), alkaline phosphatase (ALP), cadherin, trypsin, and ABC transporter genes were eminent as Bt-resistance-related genes. Expression patterns of key genes by the quantitative real-time PCR (qRT-PCR) proved the credibility of transcriptome data and suggest their association in the detoxification process. Conclusion: To date, this study is the most comprehensive research presenting functional transcriptome analysis of DBM using Bt-G033A and CL combined insecticidal activity.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Intestinal proteases profiling from Anticarsia gemmatalis and their binding to inhibitors

Although the importance of intestinal hydrolases is recognized, there is little information on the intestinal proteome of lepidopterans such as Anticarsia gemmatalis. Thus, we carried out the proteomic analysis of the A. gemmatalis intestine to characterize the proteases by LC/MS. We examined the interactions of proteins identified with protease inhibitors (PI) using molecular docking. We found 54 expressed antigens for intestinal protease, suggesting multiple important isoforms. The hydrolytic arsenal featured allows for a more comprehensive understanding of insect feeding. The docking analysis showed that the soybean PI (SKTI) could bind efficiently with the trypsin sequences and, therefore, insect resistance does not seem to involve changing the sequences of the PI binding site. In addition, a SERPIN was identified and the interaction analysis showed the inhibitor binding site is in contact with the catalytic site of trypsin, possibly acting as a regulator. In addition, this SERPIN and the identified PI sequences can be targets for the control of proteolytic activity in the caterpillar intestine and serve as a support for the rational design of a molecule with greater stability, less prone to cleavage by proteases and viable for the control of insect pests such as A. gemmatalis.

Comparative proteomics reveals mechanisms that underlie insecticide resistance in Culex pipiens pallens Coquillett

Mosquito control based on chemical insecticides is considered as an important element of the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of insecticide resistance of important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. In contrast to target site resistance, other mechanisms are far from being fully understood. Global protein profiles among cypermethrin-resistant, propoxur-resistant, dimethyl-dichloro-vinyl-phosphate-resistant and susceptible strain of Culex pipiens pallens were obtained and proteomic differences were evaluated by using isobaric tags for relative and absolute quantification labeling coupled with liquid chromatography/tandem mass spectrometric analysis. A susceptible strain of Culex pipiens pallens showed elevated resistance levels after 25 generations of insecticide selection, through iTRAQ data analysis detected 2,502 proteins, of which 1,513 were differentially expressed in insecticide-selected strains compared to the susceptible strain. Finally, midgut differential protein expression profiles were analyzed, and 62 proteins were selected for verification of differential expression using iTRAQ and parallel reaction monitoring strategy, respectively. iTRAQ profiles of adaptation selection to three insecticide strains combined with midgut profiles revealed that multiple insecticide resistance mechanisms operate simultaneously in resistant insects of Culex pipiens pallens. Significant molecular resources were developed for Culex pipiens pallens, potential candidates were involved in metabolic resistance and reducing penetration or sequestering insecticide. Future research that is targeted towards RNA interference of the identified metabolic targets, such as cuticular proteins, cytochrome P450s, glutathione S-transferases and ribosomal proteins proteins and biological pathways (drug metabolism—cytochrome P450, metabolism of xenobiotics by cytochrome P450, oxidative phosphorylation, ribosome) could lay the foundation for a better understanding of the genetic basis of insecticide resistance in Culex pipiens pallens.

Global protein profiles were compared among a susceptible strain of Cx. pipiens pallens and strains that were cypermethrin-resistant, propoxur-resistant, and dimethyl-dichloro-vinyl-phosphate-resistant after 25 generations of selection by distinct chemical insecticide families, multiple mechanisms were found to operate simultaneously in resistant mosquitoes of Cx. pipiens pallens, including mechanisms to lower penetration of or sequester the insecticide or to increase biodegradation of the insecticide via subtle alterations in either the cuticular protein levels or the activities of detoxification enzymes (P450s and glutathione S-transferases).

Honey as a Functional Food for Apis mellifera

Although nectar is consumed, primarily as a supplemental food, by a broad range of insects spanning at least five orders, it is processed and stored by only a small number of species, most of which are bees and wasps in the superfamily Apoidea. Within this group, Apis mellifera has evolved remarkable adaptations facilitating nectar processing and storage; in doing so, this species utilizes the end product, honey, for diverse functions with few if any equivalents in other phytophagous insects. Honey and its phytochemical constituents, some of which likely derive from propolis, have functional significance in protecting honey bees against microbial pathogens, toxins, and cold stress, as well as in regulating development and adult longevity. The distinctive properties of A. mellifera honey appear to have arisen in multiple ways, including genome modification; partnerships with microbial symbionts; and evolution of specialized behaviors, including foraging for substances other than nectar. That honey making by A. mellifera involves incorporation of exogenous material other than nectar, as well as endogenous products such as antimicrobial peptides and royal jelly, suggests that regarding honey as little more than a source of carbohydrates for bees is a concept in need of revision.

Structural insight into the substrate specificity of Bombyx mori β-fructofuranosidase belonging to the glycoside hydrolase family 32

Sucrose-hydrolyzing enzymes are largely divided into β-fructofuranosidase and sucrose α-glucosidase. The domestic silkworm Bombyx mori possesses both enzymes, BmSUC1 and BmSUH, belonging to the glycoside hydrolase family 32 (GH32) and GH13, respectively. BmSUC1 was presumed to be acquired by horizontal gene transfer from bacteria based on phylogenetic analysis and related to tolerance to sugar-mimic alkaloids contained in mulberry latex. Here we investigated the substrate specificity of recombinant BmSUC1 that can hydrolyze not only sucrose but also fructooligosaccharides and fructans, and revealed that the enzyme was competitively inhibited by 1,4-dideoxy-1,4-imino-D-arabinitol, one of the alkaloids. Moreover, the crystal structures of BmSUC1 in apo form and complex with sucrose were determined, and the active site pocket was shallow and suitable for shorter substrates but was related to more relaxed substrate specificity than the strict sucrose α-glucosidase BmSUH. Considering together with the distribution of BmSUC1-orthologous genes in many lepidopterans, our results suggest that BmSUC1 contributes to the digestion of fructooligosaccharides and fructans derived from feed plants.

Multi-level analysis of reproduction in an Antarctic midge identifies female and male accessory gland products that are altered by larval stress and impact progeny viability

The Antarctic midge, Belgica antarctica, is a wingless, non-biting midge endemic to Antarctica. Larval development requires at least 2 years, but adults live only 2 weeks. The nonfeeding adults mate in swarms and females die shortly after oviposition. Eggs are suspended in a gel of unknown composition that is expressed from the female accessory gland. This project characterizes molecular mechanisms underlying reproduction in this midge by examining differential gene expression in whole males, females, and larvae, as well as in male and female accessory glands. Functional studies were used to assess the role of the gel encasing the eggs, as well as the impact of stress on reproductive biology. RNA-seq analyses revealed sex- and development-specific gene sets along with those associated with the accessory glands. Proteomic analyses were used to define the composition of the egg-containing gel, which is generated during multiple developmental stages and derived from both the accessory gland and other female organs. Functional studies indicate the gel provides a larval food source as well as a buffer for thermal and dehydration stress. All of these function are critical to juvenile survival. Larval dehydration stress directly reduces production of storage proteins and key accessory gland components, a feature that impacts adult reproductive success. Modeling reveals that bouts of dehydration may have a significant impact on population growth. This work lays a foundation for further examination of reproduction in midges and provides new information related to general reproduction in dipterans. A key aspect of this work is that reproduction and stress dynamics, currently understudied in polar organisms, are likely to prove critical in determining how climate change will alter their survivability.

Identification and Functional Characterization of a Sigma Glutathione S-Transferase CpGSTs2 Involved in λ-Cyhalothrin Resistance in the Codling Moth Cydia pomonella

The codling moth, Cydia pomonella (L.), is a quarantine pest of global significance impacting pome fruits and walnuts. It has evolved resistance to many commonly used insecticides including λ-cyhalothrin. Glutathione S-transferases (GSTs) are multifunctional enzymes playing a crucial role in the detoxification of insecticides in insects. However, the role of specific GST gene in λ-cyhalothrin resistance in C. pomonella is unclear. In this study, we identified three sigma-class genes (CpGSTs1, CpGSTs2, and CpGSTs3). These genes were ubiquitously expressed at all developmental stages, and of these, the expression level of CpGSTs2 in the larval stage was significantly higher than in the egg, pupal, and adult stages. Moreover, CpGSTs2 was predominantly expressed in the fat body while lower levels in the cuticle. In addition to exposure of larvae to LD₁₀ of λ-cyhalothrin elevating the expression level of CpGSTs2, mRNA levels of CpGSTs2 in a field population (ZW_R) from northeast China, which has developed moderate level resistance to λ-cyhalothrin, was significantly higher than that of susceptible strains. In vitro inhibition assays demonstrated that λ-cyhalothrin inhibited the conjugating activities of recombinant CpGSTs2, and metabolic assays indicated that λ-cyhalothrin could be depleted by recombinant CpGSTs2. These results bring evidence for the involvement of CpGSTs2 in C. pomonella in resistance to λ-cyhalothrin.

Digestion-related proteins in the tobacco hornworm, Manduca sexta

Food digestion is vital for the survival and prosperity of insects. Research on insect digestive enzymes yields knowledge of their structure and function, and potential targets of antifeedants to control agricultural pests. While such enzymes from pest species are more relevant for inhibitor screening, a systematic analysis of their counterparts in a model insect has broader impacts. In this context, we identified a set of 122 digestive enzyme genes from the genome of Manduca sexta, a lepidopteran model related to some major agricultural pests. These genes encode hydrolases of proteins (85), lipids (20), carbohydrates (16), and nucleic acids (1). Gut serine proteases (62) and their noncatalytic homologs (11) in the S1A subfamily are encoded by abundant transcripts whose levels correlate well with larval feeding stages. Aminopeptidases (10), carboxypeptidases (10), and other proteases (3) also participate in dietary protein digestion. A large group of 11 lipases as well as 9 esterases are probably responsible for digesting lipids in diets. The repertoire of carbohydrate hydrolases (16) is relatively small, including two amylases, three maltases, two sucrases, two α-glucosidases, and others. Lysozymes, peptidoglycan amidases, and β-1,3-glucanase may hydrolyze peptidoglycans and glucans to harvest energy and defend the host from microbes on plant leaves. One alkaline nuclease is associated with larval feeding, which is likely responsible for hydrolyzing denatured DNA and RNA undergoing autolysis at a high pH of midgut. Proteomic analysis of the ectoperitrophic fluid from feeding larvae validated at least 131 or 89% of the digestive enzymes and their homologs. In summary, this study provides for the first time a holistic view of the digestion-related proteins in a lepidopteran model insect and clues for comparative research in lepidopteran pests and beyond.

Serpin7 controls egg diapause of migratory locust (Locusta migratoria) by regulating polyphenol oxidase

Diapause is a state of arrested growth, which allows insects to adapt to diverse environments. Serine protease inhibitors (serpins) play an important role in various physiological processes, including blood coagulation, fibrinolysis, development, complement activation and extracellular matrix remodeling. We hypothesized that serpin may affect energy metabolism and thereby control diapause of migratory locust (Locusta migratoria) embryos by regulating protease cascades. A total of seven nonredundant serpin genes (named serpin1–serpin7) of L. migratoria were obtained through transcriptomic analysis. We further performed label‐free proteomic sequencing and analysis of diapause and nondiapause eggs of L. migratoria, revealing significant differences in serpin7 expression. A significant reduction in diapause rate under the short photoperiod was observed in insects treated with serpin7 double‐stranded RNA. Furthermore, knockdown of the serpin7 gene resulted in significant upregulation of the activity of polyphenol oxidase. We therefore propose that the observed serpin7 gene plays a crucial role in diapause, suggesting that control of energy metabolism may have potential as a future strategy for the reproductive control of insect pests.

Multifunctional cellulase enzymes are ancestral in Polyneoptera

Many hemimetabolous insects produce their own cellulase enzymes from the glycoside hydrolase family 9, first observed in termites and cockroaches. Phasmatodea have multiple cellulases, some of which are multifunctional and can degrade xylan or xyloglucan. To discover when these abilities evolved, we identified cellulases from the Polyneoptera sampled by the 1000 Insect Transcriptome and Evolution (1KITE) project, including all cockroach and termite transcriptomes. We hoped to identify what role enzyme substrate specificities had in the evolution of dietary specification, such as leaf-feeding or wood-feeding. Putative cellulases were identified from the transcriptomes and analysed phylogenetically. All cellulases were amplified from an exemplar set of Polyneoptera species using rapid amplification of cDNA ends PCR and heterologously expressed in an insect cell line, then tested against different polysaccharides for their digestive abilities. We identified several multifunctional xyloglucanolytic enzymes across Polyneoptera, plus a large group of cellulase-like enzymes found in nearly all insect orders with no discernible digestive ability. Multifunctional xylanolytic cellulases remain unique to Phasmatodea. The presence or absence of multifunctional enzymes does not impact dietary specification, but rather having multiple, multifunctional cellulase genes is an ancestral state for Polyneoptera and possibly Insecta. The prevalence of multifunctional cellulases in other animals demands further investigation.

Gene duplication and subsequent functional diversification of sucrose hydrolase in Papilio xuthus

Sucrose is the main product of photosynthesis in plants, providing a rich carbon and energy source for the physiological growth and development of insects. In a previous study, we identified a novel sucrose hydrolase (SUH) in the larval midgut of moths. Intriguingly, there are two copies of Suh, namely Suh1 and Suh2, in several species of butterflies. However, the biochemical characteristics of SUHs in butterflies remain unclear. In this study, we found that this duplication and subsequent diversification produced two Suh genes in Papilio xuthus. These two PxSuh genes were significantly divergent in terms of their expression pattern and enzyme properties. PxSuh messenger RNA expression was highest during the larval stage, reduced in the prepupal and pupal stages and, for PxSuh1, slightly increased again in the adult. The observed levels of PxSuh2 were overall below those of PxSuh1 amongst the development stages examined. Compared with PxSUH2, which has maintained the original gene function of maltose hydrolysis, PxSUH1 exhibits substrate specificity for sucrose with an optimum enzyme activity occurring at an alkaline pH. The data show that PxSuh1 is evolutionarily adapted for effective functioning in an alkaline digestive system. Furthermore, we find that functional diversification of Suh facilitates P. xuthus to digestive carbohydrate of host plants. Thus, our findings offer new insights into the ecological and evolutionary adaptation of digestive enzymes in butterflies.

Copper exposure enhances Spodoptera litura larval tolerance to β-cypermethrin

Environmental xenobiotics can influence the tolerance of insects to chemical insecticides. Heavy metals are widespread distributed, can be easily bio-accumulated in plants and subsequently within phytophagous insects via the food chains. However, less attention has been paid to the effect of heavy metal exposure on their insecticide tolerance. In this study, pre-exposure of copper (Cu, 25-100 mg kg^-1) significantly enhanced the subsequent tolerance of Spodoptera litura to β-cypermethrin, a widely used pyrethroid insecticide in crop field. Cytochrome P450 monooxygenases (CYPs) activities were cross-induced in larvae exposed to Cu and β-cypermethrin, while the activities of glutathione S-transferase (GST) and carboxylesterase (CarE) were not affected. Application of piperonyl butoxide (PBO), a P450 synergist, effectively impaired the tolerance to β-cypermethrin in Cu-exposed S. litura larvae with a synergistic ratio of 1.72, indicating that P450s contribute to larval tolerance to β-cypermethrin induced by Cu exposure. Among the four CYP6AB family genes examined, only larval midgut-specific CYP6AB12 was found to be cross-induced by Cu and β-cypermethrin. RNA interference (RNAi)-mediated silencing of CYP6AB12 effectively decreased the mRNA levels of the target gene, and significantly reduced the larval tolerance to β-cypermethrin following exposure to Cu. These results showed that pre-exposure of heavy metal Cu enhanced larval tolerance to β-cypermethrin in S. litura, possibly through the cross-induction of P450s. Our findings provide new insights on the relationship between heavy metals and chemical insecticides that may benefit both the risk evaluation of heavy metal contamination and development of pest management strategies.

Identification of a novel cytochrome P450 CYP3356A1 linked with insecticide detoxification in Bradysia odoriphaga

BACKGROUND

Cytochrome P450 monooxygenases play an important role in the metabolic detoxification of insecticides in insect pests. However, little is known about the role of a specific P450 gene and its responses to insecticide exposure in Bradysia odoriphaga, a major pest in Chinese chive production.

RESULTS

In this study, a novel P450 gene, CYP3356A1, was cloned from Bradysia odoriphaga. The full-length cDNA sequence of CYP3356A1 is 2153 bp and its open reading frame (ORF) encodes 508 amino acids. Quantitative real time PCR(qRT-PCR) analyses in different tissues showed that CYP3356A1 expression was the highest in the Malpighian tubule. Moreover, among the different developmental stages of the insect, the highest expression of CYP3356A1 was found in fourth-instar larvae. Expression of CYP3356A1 was upregulated by treatment with imidacloprid, thiamethoxam, and β-cypermethrin at median lethal concentrations (LC₅₀ ). RNA interference (RNAi)-mediated silencing of CYP3356A1 significantly increased mortality by 36.90%, 25.17%, and 36.73 when fourth-instar B. odoriphaga larvae were exposed to imidacloprid, thiamethoxam, and β-cypermethrin, respectively, at the LC₅₀ dose.

CONCLUSION

These results demonstrate that CYP3356A1 is related to the detoxification of imidacloprid, thiamethoxam, and β-cypermethrin in B. odoriphaga. Moreover, the study also increased our understanding of the molecular mechanisms of insecticide detoxification in this pest insect. © 2018 Society of Chemical Industry.

CpGSTd3 is a lambda-Cyhalothrin Metabolizing Glutathione S-Transferase from Cydia pomonella (L.)

Little is known about the role of specific delta GST genes in the detoxification of lambda-cyhalothrin in the global quarantine fruit pest codling moth, Cydia pomonella (L.). Real-time quantitative PCR shows that CpGSTd3 was ubiquitously expressed at all developmental stages and is most abundant in the larval stage and lowest in the egg stage; the mRNA level of CpGSTd3 is higher in the midgut and Malpighian tubules of fourth-instar larvae and abdomens of adults than in other tissues. Exposure of fourth-instar larvae to an LD₁₀ dosage of lambda-cyhalothrin significantly induced the transcript of CpGSTd3 at 3 h, but the mRNA level was down-regulated after 12 h of treatment. Recombinant CpGSTd3 expressed in Escherichia coli was able to catalyze the conjugation of 1-chloro-2,4-dinitrobenzene (CDNB) and with an IC₅₀ value of 0.65 mM for lambda-cyhalothrin. Metabolism assays indicate that recombinant CpGSTd3 could metabolize lambda-cyhalothrin. These results suggest that CpGSTd3 is probably a lambda-cyhalothrin metabolizing GST in C. pomonella.

RNA-Seq Analyses of Midgut and Fat Body Tissues Reveal the Molecular Mechanism Underlying Spodoptera litura Resistance to Tomatine

Plants produce secondary metabolites to provide chemical defense against herbivorous insects, whereas insects can induce the expression of detoxification metabolism-related unigenes in counter defense to plant xenobiotics. Tomatine is an important secondary metabolite in tomato (Lycopersicon esculentum L.) that can protect the plant from bacteria and insects. However, the mechanism underlying the adaptation of Spodoptera litura, a major tomato pest, to tomatine in tomato is largely unclear. In this study, we first found that the levels of tomatine in tomatoes subjected to S. litura treatment were significantly increased. Second, we confirmed the inhibitory effect of tomatine on S. litura by adding moderate amounts of commercial tomatine to an artificial diet. Then, we utilized RNA-Seq to compare the differentially expressed genes (DEGs) in the midgut and fat body tissues of S. litura exposed to an artificial diet supplemented with tomatine. In total, upon exposure to tomatine, 134 and 666 genes were upregulated in the S. litura midgut and fat body, respectively. These DEGs comprise a significant number of detoxification-related genes, including 7 P450 family genes, 8 glutathione S-transferases (GSTs) genes, 6 ABC transport enzyme genes, 9 UDP-glucosyltransferases genes and 3 carboxylesterases genes. Moreover, KEGG analysis demonstrated that the upregulated genes were enriched in xenobiotic metabolism by cytochrome P450s, ABC transporters and drug metabolism by other enzymes. Furthermore, as numerous GSTs were induced by tomatine in S. litura, we chose one gene, namely GSTS1, to confirm the detoxification function on tomatine. Expression profiling revealed that GSTS1 transcripts were mainly expressed in larvae, and the levels were the highest in the midgut. Finally, when larvae were injected with double-stranded RNA specific to GSTS1, the transcript levels in the midgut and fat body decreased, and the negative effect of the plant xenobiotic tomatine on larval growth was magnified. These results preliminarily clarified the molecular mechanism underlying the resistance of S. litura to tomatine, establishing a foundation for subsequent pest control.

New insight into foregut functions of xenobiotic detoxification in the cockroach Periplaneta americana

The physiological functions of insect foregut, especially in xenobiotic detoxification, are scarcely reported because of unimportance in appearance and insufficient molecular information. The cockroach Periplaneta americana, an entomological model organism, provides perfect material to study physiological functions of foregut tissue due to its architectural feature. Through Illumina sequencing of foregut tissue from P. americana individuals (control) or insects treated with cycloxaprid, as a novel neonicotinoid insecticide, 54 193 166 clean reads were obtained and further assembled into 53 853 unigenes with an average length of 366 bp. Furthermore, the number of unigenes involved in xenobiotic detoxification was analyzed, mainly including 70 cytochrome P450s, 12 glutathione S-transferases (GSTs), seven carboxylesterases (CarEs) and seven adenosine triphosphate-binding cassette (ABC) transporters. Compared to control, the expression of 22 xenobiotic detoxification unigenes was up-regulated after cycloxaprid application, mainly containing 18 P450s, one GST, two CarEs and one ABC adenosine triphosphate transporter, indicating that the oxidation-reduction was the major reactive process to cycloxaprid application. Through quantitative real-time polymerase chain reaction analysis, the expression of selected unigenes (six P450s, one GST and one CarE) was up-regulated at least two-fold following cycloxaprid treatment, and was generally in agreement with transcriptome data. Compared to the previous midgut transcriptome of P. americana, it looks like the expressive abundance of the xenobiotic detoxification unigenes might be important factors to the detoxifying functional differences between foregut and midgut. In conclusion, insect foregut would also play important roles in the physiological processes related to xenobiotic detoxification.

Targeting detoxification genes by phloem-mediated RNAi: A new approach for controlling phloem-feeding insect pests

Many phloem-feeding insects are considered severe pests of agriculture and are controlled mainly by chemical insecticides. Continued extensive use of these inputs is environmentally undesirable, and also leads to the development of insecticide resistance. Here, we used a plant-mediated RNA interference (RNAi) approach, to develop a new control strategy for phloem-feeding insects. The approach aims to silence "key" detoxification genes, involved in the insect's ability to neutralize defensive and toxic plant chemistry. We targeted a glutathione S-transferase (GST) gene, BtGSTs5, in the phloem-feeding whitefly Bemisia tabaci, a devastating global agricultural pest. We report three major findings. First, significant down regulation of the BtGSTs5 gene was obtained in the gut of B. tabaci when the insects were fed on Arabidopsis thaliana transgenic plants expressing dsRNA against BtGSTs5 under a phloem-specific promoter. This brings evidence that phloem-feeding insects can be efficiently targeted by plant-mediated RNAi. Second, in-silico and in-vitro analyses indicated that the BtGSTs5 enzyme can accept as substrates, hydrolyzed aliphatic- and indolic-glucosinolates, and produce their corresponding detoxified conjugates. Third, performance assays suggested that the BtGSTs5 gene silencing prolongs the developmental period of B. tabaci nymphs. Taken together, these findings suggest that BtGSTs5 is likely to play an important role in enabling B. tabaci to successfully feed on glucosinolate-producing plants. Targeting the gene by RNAi in Brassicaceae cropping systems, will likely not eliminate the pest populations from the fields but will significantly reduce their success over the growing season, support prominent activity of natural enemies, eventually allowing the establishment of stable and sustainable agroecosystem.

Alimentary Tract Transcriptome Analysis of the Tea Geometrid, Ectropis oblique (Lepidoptera: Geometridae)

Ectropis oblique Prout (Lepidoptera: Geometridae) is one of the main pests that damages the tea crop in Southeast Asia. To understand the molecular mechanisms of its feeding biology, transcriptomes of the alimentary tract (AT) and of the body minus the AT of E. oblique were successfully sequenced and analyzed in this study. A total of 36,950 unigenes from de novo sequences were assembled. After analysis using six annotation databases (e.g., Gene Ontology, Kyoto Encyclopedia of Genes and Genome, and NCBI nr), a series of putative genes were found for this insect species that were related to digestion, detoxification, the immune system, and Bacillus thuringiensis (Bt) receptors. From this series of genes, 21 were randomly selected to verify the relative expression levels of transcripts using quantitative real-time polymerase chain reaction. These results will provide an invaluable genomic resource for future studies on the molecular mechanisms of E. oblique, which will be useful in developing biological control strategies for this pest.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes

Background

The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation.

Method

A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome.

Results

In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified.

Conclusions

A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

BmSUC1 is essential for glycometabolism modulation in the silkworm, Bombyx mori

Sucrose is the most commonly transported sugar in plants and is easily assimilated by insects to fulfill the requirement of physiological metabolism. BmSuc1 is a novel animal β-fructofuranosidase (β-FFase, EC 3.2.1.26)-encoding gene that was firstly cloned and identified in silkworm, Bombyx mori. BmSUC1 was presumed to play an important role in the silkworm-mulberry enzymatic adaptation system by effectively hydrolyzing sucrose absorbed from mulberry leaves. However, this has not been proved with direct evidence thus far. In this study, we investigated sucrose hydrolysis activity in the larval midgut of B. mori by inhibition tests and found that sucrase activity mainly stemmed from β-FFase, not α-glucosidase. Next, we performed shRNA-mediated transgenic RNAi to analyze the growth characteristics of silkworm larvae and variations in glycometabolism in vivo in transgenic silkworms. The results showed that in the RNAi-BmSuc1 transgenic line, larval development was delayed, and their body size was markedly reduced. Finally, the activity of several disaccharidases alone in the midgut and the sugar distribution, total sugar and glycogen in the midgut, hemolymph and fat body were then determined and compared. Our results demonstrated that silencing BmSuc1 significantly reduced glucose and apparently activated maltase and trehalase in the midgut. Together with a clear decrease in both glycogen and trehalose in the fat body, we conclude that BmSUC1 acts as an essential sucrase by directly modulating the degree of sucrose hydrolysis in the silkworm larval midgut, and insufficient sugar storage in the fat body may be responsible for larval malnutrition and abnormal petite phenotypes.

D181A Site-Mutagenesis Enhances Both the Hydrolyzing and Transfructosylating Activities of BmSUC1, a Novel β-Fructofuranosidase in the Silkworm Bombyx mori

β-fructofuranosidase (β-FFase) belongs to the glycosyl-hydrolase family 32 (GH32), which can catalyze both the release of β-fructose from β-d-fructofuranoside substrates to hydrolyze sucrose and the synthesis of short-chain fructooligosaccharide (FOS). BmSuc1 has been cloned and identified from the silkworm Bombyx mori as a first animal type of β-FFase encoding gene. It was hypothesized that BmSUC1 plays an important role in the silkworm-mulberry adaptation system. However, there is little information about the enzymatic core sites of BmSUC1. In this study, we mutated three amino acid residues (D63, D181, and E234) that represent important conserved motifs for β-FFase activity in GH32 to alanine respectively by using site-directed mutagenesis. Recombinant proteins of three mutants and wild type BmSUC1 were obtained by using a Bac-to-Bac/BmNPV expression system and BmN cells. Enzymatic activity, kinetic properties, and substrate specificity of the four proteins were analyzed. High Performance Liquid Chromatography (HPLC) was used to compare the hydrolyzing and transfructosylating activities between D181A and wtBmSUC1. Our results revealed that the D63A and E234A mutations lost activity, suggesting that D63 and E234 are key amino acid residues for BmSUC1 to function as an enzyme. The D181A mutation significantly enhanced both hydrolyzing and transfructosylating activities of BmSUC1, indicating that D181 may not be directly involved in catalyzation. The results provide insight into the chemical catalyzation mechanism of BmSUC1 in B. mori. Up-regulated transfructosylating activity of BmSUC1 could provide new ideas for using B. mori β-FFase to produce functional FOS.

Nucleases as a barrier to gene silencing in the cotton boll weevil, Anthonomus grandis

RNA interference (RNAi) approaches have been applied as a biotechnological tool for controlling plant insect pests via selective gene down regulation. However, the inefficiency of RNAi mechanism in insects is associated with several barriers, including dsRNA delivery and uptake by the cell, dsRNA interaction with the cellular membrane receptor and dsRNA exposure to insect gut nucleases during feeding. The cotton boll weevil (Anthonomus grandis) is a coleopteran in which RNAi-mediated gene silencing does not function efficiently through dsRNA feeding, and the factors involved in the mechanism remain unknown. Herein, we identified three nucleases in the cotton boll weevil transcriptome denoted AgraNuc1, AgraNuc2, and AgraNuc3, and the influences of these nucleases on the gene silencing of A. grandis chitin synthase II (AgraChSII) were evaluated through oral dsRNA feeding trials. A phylogenetic analysis showed that all three nucleases share high similarity with the DNA/RNA non-specific endonuclease family of other insects. These nucleases were found to be mainly expressed in the posterior midgut region of the insect. Two days after nuclease RNAi-mediated gene silencing, dsRNA degradation by the gut juice was substantially reduced. Notably, after nucleases gene silencing, the orally delivered dsRNA against the AgraChSII gene resulted in improved gene silencing efficiency when compared to the control (non-silenced nucleases). The data presented here demonstrates that A. grandis midgut nucleases are effectively one of the main barriers to dsRNA delivery and emphasize the need to develop novel RNAi delivery strategies focusing on protecting the dsRNA from gut nucleases and enhancing its oral delivery and uptake to crop insect pests.

An analysis of 67 RNA-seq datasets from various tissues at different stages of a model insect, Manduca sexta

Background

Manduca sexta is a large lepidopteran insect widely used as a model to study biochemistry of insect physiological processes. As a part of its genome project, over 50 cDNA libraries have been analyzed to profile gene expression in different tissues and life stages. While the RNA-seq data were used to study genes related to cuticle structure, chitin metabolism and immunity, a vast amount of the information has not yet been mined for understanding the basic molecular biology of this model insect. In fact, the basic features of these data, such as composition of the RNA-seq reads and lists of library-correlated genes, are unclear. From an extended view of all insects, clear-cut tempospatial expression data are rarely seen in the largest group of animals including Drosophila and mosquitoes, mainly due to their small sizes.

Results

We obtained the transcriptome data, analyzed the raw reads in relation to the assembled genome, and generated heatmaps for clustered genes. Library characteristics (tissues, stages), number of mapped bases, and sequencing methods affected the observed percentages of genome transcription. While up to 40% of the reads were not mapped to the genome in the initial Cufflinks gene modeling, we identified the causes for the mapping failure and reduced the number of non-mappable reads to <8%. Similarities between libraries, measured based on library-correlated genes, clearly identified differences among tissues or life stages. We calculated gene expression levels, analyzed the most abundantly expressed genes in the libraries. Furthermore, we analyzed tissue-specific gene expression and identified 18 groups of genes with distinct expression patterns.

Conclusion

We performed a thorough analysis of the 67 RNA-seq datasets to characterize new genomic features of M. sexta. Integrated knowledge of gene functions and expression features will facilitate future functional studies in this biochemical model insect.

Electronic supplementary material

The online version of this article doi: (10.1186/s12864-017-4147-y) contains supplementary material, which is available to authorized users.

Host plant driven transcriptome plasticity in the salivary glands of the cabbage looper (Trichoplusia ni)

Generalist herbivores feed on a wide array of plants and need to adapt to varying host qualities and defenses. One of the first insect derived secretions to come in contact with the plant is the saliva. Insect saliva is potentially involved in both the pre-digestion of the host plant as well as induction/suppression of plant defenses, yet how the salivary glands respond to changes in host plant at the transcriptional level is largely unknown. The objective of this study was to determine how the labial salivary gland transcriptome varies according to the host plant on which the insect is feeding. In order to determine this, cabbage looper (Trichoplusia ni) larvae were reared on cabbage, tomato, and pinto bean artificial diet. Labial glands were dissected from fifth instar larvae and used to extract RNA for RNASeq analysis. Assembly of the resulting sequencing reads resulted in a transcriptome library for T. ni salivary glands consisting of 14,037 expressed genes. Feeding on different host plant diets resulted in substantial remodeling of the gland transcriptomes, with 4,501 transcripts significantly differentially expressed across the three treatment groups. Gene expression profiles were most similar between cabbage and artificial diet, which corresponded to the two diets on which larvae perform best. Expression of several transcripts involved in detoxification processes were differentially expressed, and transcripts involved in the spliceosome pathway were significantly downregulated in tomato-reared larvae. Overall, this study demonstrates that the transcriptomes of the salivary glands of the cabbage looper are strongly responsive to diet. It also provides a foundation for future functional studies that can help us understand the role of saliva of chewing insects in plant-herbivore interactions.

Differential expression of chemosensory-protein genes in midguts in response to diet of Spodoptera litura

While it has been well characterized that chemosensory receptors in guts of mammals have great influence on food preference, much remains elusive in insects. Insect chemosensory proteins (CSPs) are soluble proteins that could deliver chemicals to olfactory and gustatory receptors. Recent studies have identified a number of CSPs expressed in midgut in Lepidoptera insects, which started to reveal their roles in chemical recognition and stimulating appetite in midgut. In this study, we examined expression patterns in midgut of 21 Spodoptera litura CSPs (SlitCSPs) characterized from a previously reported transcriptome, and three CSPs were identified to be expressed highly in midgut. The orthologous relationships between midgut expressed CSPs in S. litura and those in Bombyx mori and Plutella xylostella also suggest a conserved pattern of CSP expression in midgut. We further demonstrated that the expression of midgut-CSPs may change in response to different host plants, and SlitCSPs could bind typical chemicals from host plant in vitro. Overall, our results suggested midgut expressed SlitCSPs may have functional roles, likely contributing to specialization and adaption to different ecosystems. Better knowledge of this critical component of the chemsensation signaling pathways in midguts may improve our understanding of food preference processes in a new perspective.

Phytocystatins: Defense Proteins against Phytophagous Insects and Acari

This review deals with phytocystatins, focussing on their potential role as defence proteins against phytophagous arthropods. Information about the evolutionary, molecular and biochemical features and inhibitory properties of phytocystatins are presented. Cystatin ability to inhibit heterologous cysteine protease activities is commented on as well as some approaches of tailoring cystatin specificity to enhance their defence function towards pests. A general landscape on the digestive proteases of phytophagous insects and acari and the remarkable plasticity of their digestive physiology after feeding on cystatins are highlighted. Biotechnological approaches to produce recombinant cystatins to be added to artificial diets or to be sprayed as insecticide-acaricide compounds and the of use cystatins as transgenes are discussed. Multiple examples and applications are included to end with some conclusions and future perspectives.

Transcriptome profiling of Diachasmimorpha longicaudata towards useful molecular tools for population management

Background

Diachasmimorpha longicaudata (Hymenoptera: Braconidae) is a solitary parasitoid of Tephritidae (Diptera) fruit flies of economic importance currently being mass-reared in bio-factories and successfully used worldwide. A peculiar biological aspect of Hymenoptera is its haplo-diploid life cycle, where females (diploid) develop from fertilized eggs and males (haploid) from unfertilized eggs. Diploid males were described in many species and recently evidenced in D. longicaudata by mean of inbreeding studies. Sex determination in this parasitoid is based on the Complementary Sex Determination (CSD) system, with alleles from at least one locus involved in early steps of this pathway. Since limited information is available about genetics of this parasitoid species, a deeper analysis on D. longicaudata’s genomics is required to provide molecular tools for achieving a more cost effective production under artificial rearing conditions.

Results

We report here the first transcriptome analysis of male-larvae, adult females and adult males of D. longicaudata using 454-pyrosequencing. A total of 469766 reads were analyzed and 8483 high-quality isotigs were assembled. After functional annotation, a total of 51686 unigenes were produced, from which, 7021 isotigs and 20227 singletons had at least one BLAST hit against the NCBI non-redundant protein database. A preliminary comparison of adult female and male evidenced that 98 transcripts showed differential expression profiles, with at least a 10-fold difference. Among the functionally annotated transcripts we detected four sequences potentially involved in sex determination and three homologues to two known genes involved in the sex determination cascade. Finally, a total of 4674SimpleSequence Repeats (SSRs) were in silico identified and characterized.

Conclusion

The information obtained here will significantly contribute to the development of D. longicaudata functional genomics, genetics and population-based genome studies. Thousands of new microsatellite markers were identified as toolkits for population genetics analysis. The transcriptome characterized here is the starting point to elucidate the molecular bases of the sex determination mechanism in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2759-2) contains supplementary material, which is available to authorized users.

Identification of the Mamestra configurata (Lepidoptera: Noctuidae) peritrophic matrix proteins and enzymes involved in peritrophic matrix chitin metabolism

The peritrophic matrix (PM) is essential for insect digestive system physiology as it protects the midgut epithelium from damage by food particles, pathogens, and toxins. The PM is also an attractive target for development of new pest control strategies due to its per os accessibility. To understand how the PM performs these functions, the molecular architecture of the PM was examined using genomic and proteomic approaches in Mamestra configurata (Lepidoptera: Noctuidae), a major pest of cruciferous oilseed crops in North America. Liquid chromatography-tandem mass spectrometry analyses of the PM identified 82 proteins classified as: (i) peritrophins, including a new class with a CBDIII domain; (ii) enzymes involved in chitin modification (chitin deacetylases), digestion (serine proteases, aminopeptidases, carboxypeptidases, lipases and α-amylase) or other reactions (β-1,3-glucanase, alkaline phosphatase, dsRNase, astacin, pantetheinase); (iii) a heterogenous group consisting of polycalin, REPATs, serpin, C-Type lectin and Lsti99/Lsti201 and 3 novel proteins without known orthologs. The genes encoding PM proteins were expressed predominantly in the midgut. cDNAs encoding chitin synthase-2 (McCHS-2), chitinase (McCHI), and β-N-acetylglucosaminidase (McNAG) enzymes, involved in PM chitin metabolism, were also identified. McCHS-2 expression was specific to the midgut indicating that it is responsible for chitin synthesis in the PM, the only chitinous material in the midgut. In contrast, the genes encoding the chitinolytic enzymes were expressed in multiple tissues. McCHS-2, McCHI, and McNAG were expressed in the midgut of feeding larvae, and NAG activity was present in the PM. This information was used to generate an updated model of the lepidopteran PM architecture.

Horizontal Gene Transfer Contributes to the Evolution of Arthropod Herbivory

Within animals, evolutionary transition toward herbivory is severely limited by the hostile characteristics of plants. Arthropods have nonetheless counteracted many nutritional and defensive barriers imposed by plants and are currently considered as the most successful animal herbivores in terrestrial ecosystems. We gather a body of evidence showing that genomes of various plant feeding insects and mites possess genes whose presence can only be explained by horizontal gene transfer (HGT). HGT is the asexual transmission of genetic information between reproductively isolated species. Although HGT is known to have great adaptive significance in prokaryotes, its impact on eukaryotic evolution remains obscure. Here, we show that laterally transferred genes into arthropods underpin many adaptations to phytophagy, including efficient assimilation and detoxification of plant produced metabolites. Horizontally acquired genes and the traits they encode often functionally diversify within arthropod recipients, enabling the colonization of more host plant species and organs. We demonstrate that HGT can drive metazoan evolution by uncovering its prominent role in the adaptations of arthropods to exploit plants.

Tissue-specific Proteogenomic Analysis of Plutella xylostella Larval Midgut Using a Multialgorithm Pipeline

The diamondback moth, Plutella xylostella (L.), is the major cosmopolitan pest of brassica and other cruciferous crops. Its larval midgut is a dynamic tissue that interfaces with a wide variety of toxicological and physiological processes. The draft sequence of the P. xylostella genome was recently released, but its annotation remains challenging because of the low sequence coverage of this branch of life and the poor description of exon/intron splicing rules for these insects. Peptide sequencing by computational assignment of tandem mass spectra to genome sequence information provides an experimental independent approach for confirming or refuting protein predictions, a concept that has been termed proteogenomics. In this study, we carried out an in-depth proteogenomic analysis to complement genome annotation of P. xylostella larval midgut based on shotgun HPLC-ESI-MS/MS data by means of a multialgorithm pipeline. A total of 876,341 tandem mass spectra were searched against the predicted P. xylostella protein sequences and a whole-genome six-frame translation database. Based on a data set comprising 2694 novel genome search specific peptides, we discovered 439 novel protein-coding genes and corrected 128 existing gene models. To get the most accurate data to seed further insect genome annotation, more than half of the novel protein-coding genes, i.e. 235 over 439, were further validated after RT-PCR amplification and sequencing of the corresponding transcripts. Furthermore, we validated 53 novel alternative splicings. Finally, a total of 6764 proteins were identified, resulting in one of the most comprehensive proteogenomic study of a nonmodel animal. As the first tissue-specific proteogenomics analysis of P. xylostella, this study provides the fundamental basis for high-throughput proteomics and functional genomics approaches aimed at deciphering the molecular mechanisms of resistance and controlling this pest.

Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp

Adaption to dietary changes is critical in the evolution of host plant ranges in polyphagous insects. We compared three taxa of lepidopteran herbivores from the predominantly generalist genus Spodoptera showing different degrees of polyphagy: Spodoptera littoralis, with a broad host range including both mono- and dicotyledonous plants, and two Spodoptera frugiperda strains [Corn (i.e. maize) (C) and Rice (R)] adapted primarily to different grass species. When feeding on maize we show a lower performance in the broad generalist taxon compared to the grass adapted taxa. Among these taxa, the maize adapted S. frugiperda C-strain generally performed better than the R-strain on maize leaves. On artificial pinto diet, all taxa performed well. Our RNA-Seq analysis of midgut transcriptomes from 3rd instar larvae feeding on maize showed broader transcriptional readjustments in the generalist S. littoralis compared to grass adapted S. frugiperda strains. Substantial alteration in the expression levels of midgut physiological function related transcripts, such as digestive and detoxifying enzymes, transporters, immunity, and peritrophic membrane associated transcripts, existed in all taxa. We found high background expression of UDP-glucosyl transferases, which are known to neutralize maize leaf toxins, in the maize adapted S. frugiperda C-strain, contributing to its fitness on maize compared to the R-strain. Our findings provide evidence for divergent diet specific response of digestive physiology within these Spodoptera taxa. Unexpectedly, the C- and R-strains of S. frugiperda fed on the same diet showed large differences in expression patterns between these two closely related taxa.

Deltamethrin-Mediated Toxicity and Cytomorphological Changes in the Midgut and Nervous System of the Mayfly Callibaetis radiatus

Immature instars of mayflies are important constituents of the food web in aquatic ecosystems (especially in Neotropical regions) and they are among the most susceptible arthropods to pyrethroid insecticides. These insecticides have been recognized as important stressors of freshwater ecosystems, but their cellular effects in aquatic insects have been neglected. Here, we assessed the susceptibility to deltamethrin (a typical type II pyrethroid) as well as the deltamethrin-mediated cytomorphological changes in the central nervous system and midgut of the mayfly Callibaetis radiatus. While the deltamethrin LC50 for 24 h of exposure was of 0.60 (0.46-0.78) μg of a.i/L, the survival of C. radiatus was significantly reduced in deltamethrin concentrations ≥ 0.25 μg a.i/L at 96 h of exposure. Sub-lethal deltamethrin exposure severely affected the cytomorphology of C. radiatus midgut (e.g., muscle layer retraction, cytoplasm vacuolation, nucleus and striated border disorganization) and also induced slight cytomorphological changes in the brain (e.g., presence of pyknotic nuclei) and in the thoracic ganglia (e.g., vacuolation of neurons and presence of pyknotic nuclei) of these insects. However, DNA damage was absent in all of these organs, suggesting that the sublethal cellular stress induced by deltamethrin might disrupt physiological processes (e.g., metabolism or electrical signal transmission) rather than cause cell death (e.g., apoptosis) in C. radiatus. Thus, our findings indicated that deltamethrin actions at cellular levels represent a clear indication of sublethal effects on the C. radiatus survival abilities.

A specialist herbivore pest adaptation to xenobiotics through up-regulation of multiple Cytochrome P450s

The adaptation of herbivorous insects to their host plants is hypothesized to be intimately associated with their ubiquitous development of resistance to synthetic pesticides. However, not much is known about the mechanisms underlying the relationship between detoxification of plant toxins and synthetic pesticides. To address this knowledge gap, we used specialist pest Colorado potato beetle (CPB) and its host plant, potato, as a model system. Next-generation sequencing (454 pyrosequencing) was performed to reveal the CPB transcriptome. Differential expression patterns of cytochrome P450 complement (CYPome) were analyzed between the susceptible (S) and imidacloprid resistant (R) beetles. We also evaluated the global transcriptome repertoire of CPB CYPome in response to the challenge by potato leaf allelochemicals and imidacloprid. The results showed that more than half (51.2%) of the CBP cytochrome P450 monooxygenases (P450s) that are up-regulated in the R strain are also induced by both host plant toxins and pesticide in a tissue-specific manner. These data suggest that xenobiotic adaptation in this specialist herbivore is through up-regulation of multiple P450s that are potentially involved in detoxifying both pesticide and plant allelochemicals.

Only a minority of broad-range detoxification genes respond to a variety of phytotoxins in generalist Bemisia tabaci species

Generalist insect can utilize two different modes for regulating their detoxification genes, the constitutive mode and the induced mode. Here, we used the Bemisia tabaci sibling species MEAM1 and MED, as a model system for studying constitutive and induced detoxification resistance and their associated tradeoffs. B. tabaci adults were allowed to feed through membranes for 24 h on diet containing only sucrose or sucrose with various phytotoxins. Quantitative real-time PCR analyses of 18 detoxification genes, indicated that relatively few transcripts were changed in both the MEAM1 and MED species, in response to the addition of phytotoxins to the diet. Induced transcription of detoxification genes only in the MED species, in response to the presence of indole-3-carbinol in the insect’s diet, was correlated with maintenance of reproductive performance in comparison to significant reduction in performance of the MEAM1 species. Three genes, COE2, CYP6-like 5 and BtGST2, responded to more than one compound and were highly transcribed in the insect gut. Furthermore, functional assays showed that the BtGST2 gene encodes a protein capable of interacting with both flavonoids and glucosinolates. In conclusion, several detoxification genes were identified that could potentially be involved in the adaptation of B. tabaci to its host plants.

GLUTATHIONE S-TRANSFERASE Genes IN THE RICE LEAFFOLDER, Cnaphalocrocis medinalis (LEPIDOPTERA: PYRALIDAE): IDENTIFICATION AND EXPRESSION PROFILES

In insects, glutathione S-transferases (GSTs) play critical roles in the detoxification of various insecticides, resulting in insecticide resistance. The rice leaffolder, Cnaphalocrocis medinalis, is an economically important pest of rice in Asia. GST genes have not been largely identified in this insect species. In the present study, by searching the transcriptome dataset, 25 candidate GST genes were identified in C. medinalis for the first time. Of these, 23 predicted GST proteins fell into five cytosolic classes (delta, epsilon, omega, sigma, and zeta), and two were assigned to the "unclassified" subgroup. Real-time quantitative PCR analysis showed that these GST genes were differentially expressed in various tissues, including the midgut, Malpighian tubules, and fat body of larvae, and the antenna, abdomen, and leg of adults, indicating diversified functions for these genes. Transcription levels of CmGSTd2, CmGSTe6, and CmGSTe7 increased significantly in larvae following exposure to chlorpyrifos, suggesting that these GST genes could be involved in the detoxification of this insecticide. The results of our study pave the way to a better understanding of the detoxification system of C. medinalis.

The plastic response of Manduca sexta to host and non-host plants

Specialist insect herbivores have evolved efficient ways to adapt to the major defenses of their host plants. Although Manduca sexta, specialized on Solanaceous plants, has become a model organism for insect molecular biology, little is known about its adaptive responses to the chemical defenses of its hosts. To study larval performance and transcriptomic responses to host and non-host plants, we conducted developmental assays and replicated RNAseq experiments with Manduca larvae fed on different Solanaceous plants as well as on a Brassicaceous non-host plant, Brassica napus. Manduca larvae developed fastest on Nicotiana attenuata, but no significant differences in performance were found on larvae fed on other Solanaceae or the non-host B. napus. The RNAseq experiments revealed that Manduca larvae display plastic responses at the gene expression level, and transcriptional signatures specific to the challenges of each host- and non-host plant. Our observations are not consistent with expectations that specialist herbivores would perform poorly on non-host plants. Instead, our findings demonstrate the ability of this specialized insect herbivore to efficiently use a larger repertoire of host plants than it utilizes in the field.