De novo characterization of transcriptome and gene expression dynamics in epidermis during the larval-pupal metamorphosis of common cutworm

Transcriptome-Wide Evaluation Characterization of microRNAs and Assessment of Their Functional Roles as Regulators of Diapause in Ostrinia furnacalis Larvae (Lepidoptera: Crambidae)

microRNAs (miRNAs) function as vital regulators of diapause in insects through their ability to post-transcriptionally suppress target gene expression. In this study, the miRNA of Ostrinia furnacalis, an economically important global crop pest species, was characterized. For the included analyses, 9 small RNA libraries were constructed using O. furnacalis larvae in different diapause states (non-diapause, ND; diapause, D; diapause-termination, DT). The results identified 583 total miRNAs, of which 256 had previously been identified, whereas 327 were novel. Furthermore, comparison analysis revealed that 119 and 27 miRNAs were differentially expressed in the D vs. ND and DT vs. D, respectively. Moreover, the expression patterns of their miRNAs were also analyzed. GO and KEGG analysis of the target genes of differentially expressed miRNAs highlighted the importance of these miRNAs as diapause regulators in O. furnacalis, especially through metabolic processes, endocrine processes, 20-hydroxyecdysone, and circadian clock signaling pathways. In summary, this study highlighted the involvement of specific miRNAs in the control of diapause in O. furnacalis. To the best of our knowledge, this is the first study to identify miRNA expression patterns in O. furnacalis, thereby providing reference and novel evidence enhancing our current understanding of how small RNAs influence insect diapause.

Cytoplasmic FKBPs are involved in molting and metamorphosis through regulating the nuclear localization of EcR

Molting and metamorphosis are important physiological processes in insects that are tightly controlled by ecdysone receptor (EcR) through the 20-hydroxyecdysone (20E) signaling pathway. EcR is a steroid nuclear receptor (SR). Several FK506-binding proteins (FKBPs) have been identified from the mammal SR complex, and are thought to be involved in the subcellular trafficking of SR. However, their roles in insects are poorly understood. To explore whether FKBPs are involved in insect molting or metamorphosis, we injected an FKBP inhibitor (FK506) into a lepidopteran insect, Spodoptera litura, and found that molting was inhibited in 61.11% of the larvae, and that the time for larvae to pupate was significantly extended. A total of 10 FKBP genes were identified from the genome of S. litura and were clustered into 2 distinct groups, according to their subcellular localization, with FKBP13 and FKBP14 belonging to the endoplasmic reticulum (ER) group and with the other members belonging to the cytoplasmic (Cy) group. All the CyFKBPs were significantly upregulated in the prepupal or pupal stages, with the opposite being observed for the ER group members. FK506 completely blocked the transfer of EcR to the nucleus under 20E induction, and significantly downregulated the transcriptional expression of many 20E signaling genes. A similar phenomenon was observed after RNA interference of 2 CyFKBPs (FKBP45 and FKBP12b), but not for FKBP13. Taken together, our data indicate that the cytoplasmic FKBPs, especially FKBP45 and FKBP12b, mediate the nuclear localization of EcR, thereby regulating the 20E signaling and ultimately affecting molting and metamorphosis in insects.

Tyrosine hydroxylase plays crucial roles in larval cuticle formation and larval-pupal tanning in the rice stem borer, Chilo suppressalis

The striped stem borer, Chilo suppressalis (Walker), a notorious pest infesting rice, has evolved a high level of resistance to many commonly used insecticides. In this study, we investigate whether tyrosine hydroxylase (TH), which is required for larval development and cuticle tanning in many insects, could be a potential target for the control of C. suppressalis. We identified and characterized the full-length cDNA (CsTH) of C. suppressalis. The complete open reading frame of CsTH (MW690914) was 1683 bp in length, encoding a protein of 560 amino acids. Within the first to the sixth larval instars, CsTH was high in the first day just after molting, and lower in the ensuing days. From the wandering stage to the adult stage, levels of CSTH began to rise and reached a peak at the pupal stage. These patterns suggested a role for the gene in larval development and larval-pupal cuticle tanning. When we injected dsCsTH or 3-iodotyrosine (3-IT) as a TH inhibitor or fed a larva diet supplemented with 3-IT, there were significant impairments in larval development and larval-pupal cuticle tanning. Adult emergence was severely impaired, and most adults died. These results suggest that CsTH might play a critical role in larval development as well as larval-pupal tanning and immunity in C. suppressalis, and this gene could form a potential novel target for pest control.

Genome-wide analysis of long noncoding RNAs and their association in regulating the metamorphosis of the Sarcophaga peregrina (Diptera: Sarcophagidae)

BACKGROUND

The flesh fly, Sarcophaga peregrina (Diptera: Sarcophagidae), is an important hygiene pest, that causes myiasis in humans and other mammals, typically livestock, and as a vector for various parasitic agents, including bacteria, viruses, and parasites. The role of long non-coding RNAs (lncRNAs) in regulating gene expression during metamorphosis of the flesh fly has not been well established.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed genome-wide identification and characterization of lncRNAs from the early pupal stage (1-days pupae), mid-term pupal stage (5-days pupae), and late pupal stage (9-days pupae) of S. peregrina by RNA-seq, and a total of 6921 lncRNAs transcripts were identified. RT-qPCR and enrichment analyses revealed the differentially expressed lncRNAs (DE lncRNAs) that might be associated with insect metamorphosis development. Furthermore, functional analysis revealed that the DE lncRNA (SP_lnc5000) could potentially be involved in regulating the metamorphosis of S. peregrina. RNA interference of SP_lnc5000 caused reduced expression of metamorphosis-related genes in 20-hydroxyecdysone (20E) signaling (Br-c, Ftz-F1), cuticle tanning pathway (TH, DOPA), and chitin related pathway (Cht5). Injection of dsSP_lnc5000 in 3rd instar larvae of S. peregrina resulted in deformed pupae, stagnation of pupal-adult metamorphosis, and a decrease in development time of pupal, pupariation rates and eclosion rates. Hematoxylin-eosin staining (H&E), scanning electron microscope (SEM) observation and cuticle hydrocarbons (CHCs) analysis indicated that SP_lnc5000 had crucial roles in the metamorphosis developmental by modulating pupal cuticular development.

CONCLUSIONS/SIGNIFICANCE

We established that the lncRNA SP_lnc5000 potentially regulates the metamorphosis of S. peregrina by putatively affecting the structure and composition of the pupal cuticle. This study enhances our understanding of lncRNAs as regulators of metamorphosis in S. peregrina, and provide valuable insights into the identification of potential targets for vector control and the development of effective strategies for controlling the spread of myiasis and parasitic diseases.

Cytoplasmic Hsp70s promote EcR transport into the nucleus by responding to various stimuli

Metamorphosis is one of the most important physiological processes in insects, which is coordinated by juvenile hormone (JH) and 20-hydroxyecdysone (20E). Ecdysone receptor (EcR) is a steroid receptor (SR), which usually presents in cytoplasm and transfers into nucleus after binding to 20E. Heat shock proteins (Hsps) are suggested to be important members of the SR complex. However, their role in nucleocytoplasmic shuttle of the EcR remains unclear. In the present study, we found that apoptozole (Hsp70 inhibitor) suppressed the larval molting by decreasing the expression of ecdysone signaling genes. Two cytoplasmic (Cy) Hsp70s (Hsp72 and Hsp73) interacted with both EcR and ultraspiracle (USP, the heterodimer partner of EcR). By immunohistochemistry experiments, we revealed that CyHsp70 co-localized with EcR in the cytoplasm, and that both apoptozole and interfering of CyHsp70 significantly inhibited the process of EcR entering the nucleus under 20E induction, while reducing the expression of ecdysone signaling genes. Interestingly, the nuclear localization of EcR was also promoted by two other stimuli, including JH and heat stress, and this promotion was inhibited by apoptozole. This implies that various stimuli can induce EcR entry into the nucleus, and that this process is mediated by CyHsp70. Curiously, neither JH nor heat stress activated the ecdysone signaling genes; instead, they have a significant inhibitory effect on them. Taken together, it seems that Cytoplasmic Hsp70s promote EcR transport into the nucleus by responding to various stimuli, and that the biological effects of various stimuli passing through the EcR are different. Thus, our data provide a new viewpoint to understand the mechanism of nucleocytoplasmic shuttle of EcR.

Transcriptome profiling for developmental stages Protaetia brevitarsis seulensis with focus on wing development and metamorphosis

A white-spotted flower chafer Protaetia brevitarsis seulensis widely distributed in Asian countries is traditionally used in oriental medicine. This study explored gene expression abundance with respect to wing development and metamorphosis in P. b. seulensis based on the large-scale RNA-seq data. The transcriptome assembly consists of 23,551 high-quality transcripts which are approximately 96.7% covered. We found 265 wing development genes, 19 metamorphosis genes, and 1,314 candidates. Of the 1,598 genes, 1,594 are included exclusively in cluster 4 with similar gene co-expression patterns. The network centrality analyses showed that wing development- and metamorphosis-related genes have a high degree of betweenness centrality and are expressed most highly in eggs, moderately in pupa and adults, and lowest in larva. This study provides some meaningful clues for elucidating the genetic modulation mechanism of wing development and metamorphosis in P. b. seulensis.

Pupal Age Estimation of Sarcophaga peregrina (Diptera: Sarcophagidae) at Different Constant Temperatures Utilizing ATR-FTIR Spectroscopy and Cuticular Hydrocarbons

Sarcophaga peregrina (Robineau-Desvoidy, 1830) (Diptera: Sarcophagidae) is a forensically important flesh fly that has potential value in estimating the PMImin. The precise pupal age estimation has great implications for PMImin estimation. During larval development, the age determination is straightforward by the morphological changes and variation of length and weight, however, the pupal age estimation is more difficult due to anatomical and morphological changes not being visible. Thus, it is necessary to find new techniques and methods that can be implemented by standard experiments for accurate pupal age estimation. In this study, we first investigated the potential of attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy and cuticular hydrocarbons (CHCs) for the age estimations of S. peregrina pupae at different constant temperatures (20 °C, 25 °C, and 30 °C). The orthogonal projections latent structure discrimination analysis (OPLS-DA) classification model was used to distinguish the pupae samples of different developmental ages. Then, a multivariate statistical regression model, partial least squares (PLS), was established with the spectroscopic and hydrocarbon data for pupal age estimations. We identified 37 CHCs with a carbon chain length between 11 and 35 in the pupae of S. peregrina. The results of the OPLS-DA model show a significant separation between different developmental ages of pupae (R2X > 0.928, R2Y > 0.899, Q2 > 0.863). The PLS model had a satisfactory prediction with a good fit between the actual and predicted ages of the pupae (R2 > 0.927, RMSECV < 1.268). The results demonstrate that the variation tendencies of spectroscopy and hydrocarbons were time-dependent, and ATR-FTIR and CHCs may be optimal for the age estimations of pupae of forensically important flies with implications for PMImin estimation in forensic practice.

RNA Sequencing Reveals the Potential Adaptation Mechanism to Different Hosts of Grapholita molesta

Grapholita molesta is an important fruit tree worldwide pest which feeds on hosts extensively and does serious harm. In this paper, the growth and development parameters and protease activities of G. molesta fed on different hosts were compared. Using Illumina RNA sequencing technology, 18 midgut samples from five different hosts (apple, pear, plum, peach and peach shoots) and artificial diet were sequenced and compared with the reference genome, resulting in 15269 genes and 2785 predicted new genes. From 15 comparative combinations, DEGs were found from 286 to 4187 in each group, with up-regulated genes from 107 to 2395 and down-regulated genes from 83 to 2665. KEGG pathway analysis showed that DEGs were associated with amino acid metabolism, starch and sucrose metabolism, carbohydrate metabolism, and hydrolase activity. A total of 31 co-expression gene modules of different hosts were identified by WGCNA. qRT-PCR showed that the expression pattern of the trypsin gene was consistent with RNA sequencing. In this study, growth and development parameters, protease activity, DEGs, enrichment analysis and qRT-PCR were combined to reveal the adaptation process to different hosts of G. molesta in many aspects. The results of this study provide a basis for further exploration of the molecular mechanism of host adaptation of G. molesta.

Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development

BACKGROUND

Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study.

RESULTS

A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura.

CONCLUSIONS

Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

Temporal Expression Profiles Reveal Potential Targets during Postembryonic Development of Forensically Important Sarcophaga peregrina (Diptera: Sarcophagidae)

Sarcophaga peregrina (Robineau-Desvoidy, 1830) is a species of medical and forensic importance. In order to investigate the molecular mechanism during postembryonic development and identify specific genes that may serve as potential targets, transcriptome analysis was used to investigate its gene expression dynamics from the larval to pupal stages, based on our previous de novo-assembled genome of S. peregrina. Totals of 2457, 3656, 3764, and 2554 differentially expressed genes were identified. The specific genes encoding the structural constituent of cuticle were significantly differentially expressed, suggesting that degradation and synthesis of cuticle-related proteins might actively occur during metamorphosis. Molting (20-hydroxyecdysone, 20E) and juvenile (JH) hormone pathways were significantly enriched, and gene expression levels changed in a dynamic pattern during the developmental stages. In addition, the genes in the oxidative phosphorylation pathway were significantly expressed at a high level during the larval stage, and down-regulated from the wandering to pupal stages. Weighted gene co-expression correlation network analysis (WGCNA) further demonstrated the potential regulation mechanism of tyrosine metabolism in the process of puparium tanning. Moreover, 10 consistently up-regulated genes were further validated by qRT-PCR. The utility of the models was then examined in a blind study, indicating the ability to predict larval development. The developmental, stage-specific gene profiles suggest novel molecular markers for age prediction of forensically important flies.

Expression profiles of tyrosine metabolic pathway genes and functional analysis of DOPA decarboxylase in puparium tanning of Bactrocera dorsalis (Hendel)

BACKGROUND

Tanning is an important physiological process with critical roles in cuticle pigmentation and sclerotization. Previous studies have shown that insect cuticle tanning is closely associated with the tyrosine metabolism pathway, which consists of a series of enzymes.

RESULTS

In this study, 24 tyrosine metabolism pathway genes were identified in the oriental fruit fly Bactrocera dorsalis (Hendel) genome. Gene expression profiles throughout 15 developmental stages of B. dorsalis were established based on our previous RNA sequencing data, and we found that 13 enzyme genes could be involved in the process of pupariation. Accordingly, a tyrosine-mediated tanning pathway during the pupariation of B. dorsalis was predicted and a critical enzyme, 3,4-dihydroxyphenylalanine (DOPA) decarboxylase (DDC), was used to explore its possible roles in formation of the puparium. First, a real-time quantitative polymerase chain reaction confirmed that BdDDC had an epidermis-specific expression pattern, and was highly expressed during larval metamorphosis in B. dorsalis. Subsequent disruption of BdDDC by feeding 5-day-old larvae with DDC inhibitor (l-α-methyl-DOPA) could lead to: (i) a significant decrease in BdDDC enzyme activity and dopamine concentration; (ii) defects in puparium pigmentation; (iii) impairment of the morphology and less thickness of the puparium; and (iv) lower pupal weight and obstacles to eclosion.

CONCLUSION

This study provided a potential tyrosine metabolic pathway that was responsible for insect tanning during pupariation, and the BdDDC enzyme has been shown to have crucial roles in larval-pupal tanning of B. dorsalis. © 2021 Society of Chemical Industry.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

Knockdown of specific cuticular proteins analogous to peritrophin 3 genes disrupt larval and ovarian development in Bactrocera dorsalis (Diptera: Tephritidae)

Cuticular proteins (CPs) are critical components of the insect cuticle and play important roles in maintaining normal insect development and defense against various environmental stresses. The oriental fruit fly (Bactrocera dorsalis) is one of the most destructive pests worldwide, and its eight CPs analogous to peritrophin 3 (BdCPAP3) family genes have been identified in our previous study. In the present study, we further explored the possible roles of CPAP3 genes in B. dorsalis development. Each sequence of BdCPAP3 genes contained three conserved ChtBD2 (chitin-binding) domains. Spatial and temporal expression patterns revealed that the four BdCPAP3 genes (BdCPAP3-A1, B, E, and E2) might play important roles in larval pupariation of B. dorsalis. Moreover, treatment with a juvenile hormone analog (methoprene) significantly restricted expression of these four CPAP3 genes, whereas treatment with 20-hydroxy-ecdysone induced expression. The RNA interference (RNAi) results revealed that down-regulated CPAP3 genes led to significant delay of pupariation, and injection of dsBdCPAP3-E into 5-d-old B. dorsalis larvae caused approximately 40% mortality. Interestingly, we also confirmed that BdCPAP3-D2 was involved in B. dorsalis ovarian development. This study showed that some specific CPAP3 genes had crucial roles in B. dorsalis development, and these CP genes could be used as potential targets to control this pest via RNAi.

Identification and expression analysis of cuticular protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

Structural cuticular proteins (CPs) are major components of the insect cuticle, and they play critical roles in insect development and insecticide resistance. Here, a total of 196 CP genes were successfully annotated in the Plutella xylostella genome. On the basis of motif analysis, these CPs were classified into 10 different families, including 122 CPR, 12 CPAP1, 8 CPAP3, 9 CPLCP, 2 Tweedle, 1 CPF, 1 CPFL, 1 CPCFC, 17 CPG and 2 18 aa proteins, and the remaining 21 unclassified CPs were classed as cuticular proteins hypothetical (CPH). A phylogenetic analysis of CPs from different insects revealed species-specific clades of RR-1 and RR-2 genes, suggesting that CP gene duplication might occur independently among insect taxa, while we also found that some other CPs (such as CPAP1 and CPAP3) had a closer relationship based on their conserved domain architecture. Using available RNAseq libraries, the expression profiles of the CPs were analyzed over the four developmental stages of the insect (i.e., egg, larva, pupa, and adult), revealing stage-specific expression patterns for the CPs. In a chlorpyrifos resistant strain, 18 CP genes were found to be more than two-fold upregulated compared to the susceptible control strain, and qRT-PCR analysis showed that these CP genes were overexpressed after exposure to chlorpyrifos, suggesting a potential role in the molecular mechanism of insecticide resistance in P. xylostella. This study provides the tools and molecular basis to study the role of CPs in the post-embryonal development and the mechanisms of insecticide resistance of P. xylostella.

Transcriptional identification of differentially expressed genes during the prepupal-pupal transition in the oriental armyworm, Mythimna separata (Walker) (Lepidoptera: Noctuidae)

The oriental armyworm, Mythimna separata (Walker) is a serious pest of agriculture that does particular damage to Gramineae crops in Asia, Europe, and Oceania. Metamorphosis is a key developmental stage in insects, although the genes underlying the metamorphic transition in M. separata remain largely unknown. Here, we sequenced the transcriptomes of five stages; mature larvae (ML), wandering (W), and pupation (1, 5, and 10 days after pupation, designated P1, P5, and P10) to identify transition-associated genes. Four libraries were generated, with 22,884, 23,534, 26,643, and 33,238 differentially expressed genes (DEGs) for the ML-vs-W, W-vs-P1, P1-vs-P5, and P5-vs-P10, respectively. Gene ontology enrichment analysis of DEGs showed that genes regulating the biosynthesis of the membrane and integral components of the membrane, which includes the cuticular protein (CP), 20-hydroxyecdysone (20E), and juvenile hormone (JH) biosynthesis, were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that DEGs were enriched in the metabolic pathways. Of these DEGs, thirty CP, seventeen 20E, and seven JH genes were differentially expressed across the developmental stages. For transcriptome validation, ten CP, 20E, and JH-related genes were selected and verified by real-time PCR quantitative. Collectively, our results provided a basis for further studies of the molecular mechanism of metamorphosis in M. separata.

RNA-seq analysis of gene expression changes in cuticles during the larval-pupal metamorphosis of Plutella xylostella

The diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae) is a holometabolous insect that its cuticles must undergo the significant changes during the larval-pupal metamorphosis development. To elucidate these changes at molecular levels, RNA-seq analysis of cuticles from LLS (later fourth instar larval stage), PPS (prepupal stage) and PS (pupal stage) were performed in P. xylostella. In this paper, a total of 17,710 transcripts were obtained in the larval-pupal transition of P. xylostella, and out of which 2293 (881 up-regulated and 1412 down-regulated) and 2989 transcripts (2062 up-regulated and 927 down-regulated) were identified to be differentially expressed between LLS and PPS, as well as PPS and PS, respectively. The further GO and KEGG analysis of differentially expressed genes (DEGs) revealed that the 'structural constituent of cuticle', 'chitin metabolic process', 'chitin binding', 'tyrosine metabolism' and 'insect hormone biosynthesis' pathways were significantly enriched, indicating these pathways might be involved in the process of larval pupation in P. xylostella. Then, we found some genes that encoded cuticular proteins, chitinolytic enzymes, chitin synthesis enzymes, and cuticle tanning proteins changed their expression levels remarkably, indicating these genes might play important roles in the restruction (degradation and biosynthesis) of insect cuticles during the larval metamorphosis. Additionally, the significant changes in the mRNA levels of 20-hydroxyecdysone (20E) and juvenile hormone (JH) related genes suggested their crucial roles in regulating cuticle remodeling during the larval metamorphosis of P. xylostella. In conclusion, the present study provide us the comprehensive gene expression profiles to explore the molecular mechanisms of cuticle metamorphosis in P. xylostella, which laid a molecular basis to study roles of specific pathways and genes in insect development.

Gene Expression Differences Between Developmental Stages of the Fall Armyworm (Spodoptera frugiperda)

The fall armyworm (Spodoptera frugiperda) is one of the most significant agricultural pests in the world and invaded China in early 2019. We sampled and sequenced RNA-seq data from 15 individuals across different developmental stages. Developmental stages were the larval stage (5th instar larvae and 6th instar larvae), chrysalis stage, and adult stage (female adult and male adult). Individual samples were mainly clustered by developmental stages and we then identified variation between developmental stages of differentially expressed transcripts (DETs). There were 2136 upregulated DETs and 1391 downregulated DETs in the larval stage when comparing larval and chrysalis stages. In the comparison between the chrysalis and adult stages, there were 2033 upregulated DETs and 1391 downregulated DETs in the chrysalis stage. In total, 19,195 abundantly expressed transcripts were obtained and 10% of them were DETs. We then obtained stage-specific DETs to investigate the potential function of the fall armyworm during different developmental stages. We also constructed our annotation background set for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. This indicated that the fall armyworm may undergo active metabolism during its lifespan, even in the chrysalis stage. And it also may experience detoxifying and xenobiotic metabolism throughout its life, especially in the larval stage, which partially explains the difficulty to eradicate using chemical control. Our study is the first insight into the developmental patterns of the fall armyworm and we also provide the fundamental information about enhanced drug resistance at the level of transcriptome. These results are beneficial for a future investigation related to the eradication and/or control stage.

RNA-seq Analysis Reveals Mitochondrial and Cuticular Protein Genes Are Associated with Phosphine Resistance in the Rusty Grain Beetle (Coleoptera:Laemophloeidae)

The rusty grain beetle, Cryptolestes ferrugineus (Stephens), is a serious pest of stored grain, which has developed high levels of resistance to phosphine. In this study, five geographically distant populations of C. ferrugineus had been collected in China, specifically in granaries where phosphine fumigant is used for pest control, and they showed a high resistance ratio up to 1,907 (LC50 = 21.0 mg/liter). Then, a reference transcriptome was constructed to use as a basis for investigating the molecular mechanisms of phosphine resistance in this species, which consisted of 47,006 unigenes with a mean length of 1,090. Subsequently, the RNA-Seq analysis of individuals from the most susceptible and resistant populations led to the identification of 54 genes that are differentially expressed. GO and KEGG analysis demonstrated that genes associated with mitochondrial and respiration functions were significantly enriched. Also, the 'structural constituent of cuticle' term was annotated in the GO enrichment analysis and further qRT-PCR confirmed that the expression levels of nine cuticular protein genes were significantly increased in the resistant population. In conclusion, we present here a transcriptome-wide overview of gene expression changes between resistant and susceptible populations of C. ferrugineus, and this in turn documents that mitochondria and cuticular protein genes may play together a crucial role in phosphine resistance. Further gene function analysis should enable the provision of advice to expedite resistance management decisions.

Multiple microRNAs control ecdysone signaling in the midgut of Spodoptera litura

Metamorphosis is one of the most important physiological processes in insects. It is regulated by a serial of ecdysone cascade genes. Recently, lots of microRNAs (miRNAs) were investigated in insects; however, their function in metamorphosis is largely unknown. In the present study, the dynamics of a small RNA population was investigated by RNA sequencing from the midgut of a lepidopteran pest Spodoptera litura during larval-pupal metamorphosis. A total of 101 miRNAs were identified, and 75 miRNAs were differentially expressed during the metamorphic process. The relationship between these differentially expressed miRNAs and 12 ecdysone cascade genes was analyzed by four classical software programs, and a multiple-to-multiple regulatory network was found to exist between these miRNAs and their targets. Among them, miR-14-3p and its two targets (EcR and E75) were chosen for further validation. MiR-14-3p had higher expression level in the 6th instar larvae as compared with either the prepupae or pupae, which was opposite to that of both EcR and E75, two ecdysone cascade genes. Luciferase reporter assay confirmed that both EcR and E75 were regulated by miR-14-3p. Interestingly, the 3' untranslated regions are nearly identical to each other among different transcript variants of the ecdysone cascade genes, including EcR, USP, E75, E74, E78, E93, Hr3, Hr4, Hr39, Krh1 and Ftzf1. Thus, different transcript variants of one ecdysone cascade gene could be regulated by the same miRNA. The above data suggest that the ecdysone signaling pathway is under the tight control of miRNA. These findings expand our understanding of the mechanism of insect metamorphosis and may also provide a novel possibility for the control of pest insects in the future.

Identification and temporal expression profiles of cuticular proteins in the endoparasitoid wasp, Microplitis mediator

Recently, parasitoid wasp species Microplitis mediator has evoked increasing research attention due to its possible use in the control of Lepidoptera insects. Because insect development involves changes in cuticle composition, identification and expression analysis of M. mediator cuticular proteins may clarify the mechanisms involved in parasite development processes. We found 70 cuticular proteins from the M. mediator transcriptome and divided them into seven distinct families. Expression profiling indicated that most of these cuticular protein genes have expression peaks specific for one particular developmental stage of M. mediator. Eggs and pupae have the highest number of transcriptionally active cuticular protein genes (47 and 52 respectively). Only 12 of these genes maintained high expression activity during late larval development. Functional analysis of two larval proteins, MmCPR3 and MmCPR14, suggested their important role in the proper organization of the cuticle layers of larvae. During M. mediator larval development, normal cuticle formation can be supported by a limited number of cuticular proteins.

Transcriptome-Wide Identification of MicroRNAs and Analysis of Their Potential Roles in Development of Indian Meal Moth (Lepidoptera:Pyralidae)

MicroRNAs (miRNAs) have been reported to play indispensable roles in regulating various developmental processes via the posttranscriptional repression of target genes in insect species. In the present paper, we studied the miRNAs in Indian meal moth (Plodia interpunctella (Hübener)), one of the most economically important stored grains pests around the world. In total, 12 small RNA libraries from four developmental stages of P. interpunctella were constructed, and 178 known and 23 novel miRNAs were identified. In addition, the expression profiles of these miRNAs were assessed across different developmental stages and miRNAs that were highly expressed in eggs, larvae, pupae, and adults were identified. Specifically, 100, 61, and 52 miRNAs were differentially expressed between eggs and larvae, larvae and pupae, and pupae and adults, respectively. The KEGG and GO analysis of the predicted target genes suggested the essential roles of miRNAs in the regulation of complex development of P. interpunctella. Importantly, we also found a set of miRNAs might be involved in the larval metamorphic molting process, with their expressions increasing and then decreasing during the larva-pupa-adult stages of P. interpunctella. In conclusion, the current paper has discovered numerous miRNAs, and some key miRNAs that might be responsible for regulating development in P. interpunctella. To our knowledge, this is the first study to document miRNAs and their expression patterns in interpunctella, and those findings would lay an important molecular foundation for future functional analysis of these miRNAs in P. interpunctella.

Genome-Wide Analysis of MicroRNAs in Relation to Pupariation in Oriental Fruit Fly

Insect metamorphosis is a complex process involving drastic morphological and physiological changes. microRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play key roles in regulating various biological processes, including metamorphosis, by post-transcriptional repression of mRNAs. The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive insect pests in many Asian countries and the Pacific Islands. The regulatory role of miRNAs in B. dorsalis metamorphosis is unclear. To better understand the molecular regulatory mechanisms of miRNAs in pupariation, Illumina sequencing of the wandering stage (WS), the late WS and the white puparium stage of B. dorsalis were performed. Two hundred forty-nine miRNAs, including 184 known miRNAs and 65 novel miRNAs, were obtained. Among these miRNAs, 19 miRNAs were differentially expressed in pupariation, and eight miRNAs showed relative high expression levels (>50 TPM), of which five differentially expressed miRNAs (DEMs) had target differentially expressed genes (DEGs) predicted by the expected miRNA-mRNA negative regulation pattern using the Illumina HiSeq data. Four sets of DEMs and their predicted target DEGs were confirmed by qPCR. Of the four miRNAs, two miRNAs were down-regulated: miR-981, which may target pdpc, and Bdo-novel-mir-55, which potentially regulates spsX1, psB/C, and chit3. The other two miRNAs were up-regulated: let-7a-3p, which possibly controls lap, and Bdo-novel-mir-24, which may regulate ipc and sp1/2. This study provides a useful resource to elucidate the regulatory role of miRNAs and understand the molecular mechanisms of metamorphosis.

Exploring integument transcriptomes, cuticle ultrastructure, and cuticular hydrocarbons profiles in eusocial and solitary bee species displaying heterochronic adult cuticle maturation

Differences in the timing of exoskeleton melanization and sclerotization are evident when comparing eusocial and solitary bees. This cuticular maturation heterochrony may be associated with life style, considering that eusocial bees remain protected inside the nest for many days after emergence, while the solitary bees immediately start outside activities. To address this issue, we characterized gene expression using large-scale RNA sequencing (RNA-seq), and quantified cuticular hydrocarbon (CHC) through gas chromatography-mass spectrometry in comparative studies of the integument (cuticle plus its underlying epidermis) of two eusocial and a solitary bee species. In addition, we used transmission electron microscopy (TEM) for studying the developing cuticle of these and other three bee species also differing in life style. We found 13,200, 55,209 and 30,161 transcript types in the integument of the eusocial Apis mellifera and Frieseomelitta varia, and the solitary Centris analis, respectively. In general, structural cuticle proteins and chitin-related genes were upregulated in pharate-adults and newly-emerged bees whereas transcripts for odorant binding proteins, cytochrome P450 and antioxidant proteins were overrepresented in foragers. Consistent with our hypothesis, a distance correlation analysis based on the differentially expressed genes suggested delayed cuticle maturation in A. mellifera in comparison to the solitary bee. However, this was not confirmed in the comparison with F. varia. The expression profiles of 27 of 119 genes displaying functional attributes related to cuticle formation/differentiation were positively correlated between A. mellifera and F. varia, and negatively or non-correlated with C. analis, suggesting roles in cuticular maturation heterochrony. However, we also found transcript profiles positively correlated between each one of the eusocial species and C. analis. Gene co-expression networks greatly differed between the bee species, but we identified common gene interactions exclusively between the eusocial species. Except for F. varia, the TEM analysis is consistent with cuticle development timing adapted to the social or solitary life style. In support to our hypothesis, the absolute quantities of n-alkanes and unsaturated CHCs were significantly higher in foragers than in the earlier developmental phases of the eusocial bees, but did not discriminate newly-emerged from foragers in C. analis. By highlighting differences in integument gene expression, cuticle ultrastructure, and CHC profiles between eusocial and solitary bees, our data provided insights into the process of heterochronic cuticle maturation associated to the way of life.

RNA-seq analysis of gene expression changes during pupariation in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

Background

The oriental fruit fly, Bactrocera dorsalis (Hendel) has been considered to be one of the most important agricultural pest around the world. As a holometabolous insect, larvae must go through a metamorphosis process with dramatic morphological and structural changes to complete their development. To better understand the molecular mechanisms of these changes, RNA-seq of B. dorsalis from wandering stage (WS), late wandering stage (LWS) and white puparium stage (WPS) were performed.

Results

In total, 11,721 transcripts were obtained, out of which 1914 genes (578 up-regulated and 1336 down-regulated) and 2047 genes (655 up-regulated and 1392 down-regulated) were found to be differentially expressed between WS and LWS, as well as between WS and WPS, respectively. Of these DEGs, 1862 and 1996 genes were successfully annotated in various databases. The analysis of RNA-seq data together with qRT-PCR validation indicated that during this transition, the genes in the oxidative phosphorylation pathway, and genes encoding P450s, serine protease inhibitor, and cuticular proteins were down-regulated, while the serine protease genes were up-regulated. Moreover, we found some 20-hydroxyecdysone (20E) biosynthesis and signaling pathway genes had a higher expression in the WS, while the genes responsible for juvenile hormone (JH) synthesis, degradation, signaling and transporter pathways were down-regulated, suggesting these genes might be involved in the process of larval pupariation in B. dorsalis. For the chitinolytic enzymes, the genes encoding chitinases (chitinase 2, chitinase 5, chitinase 8, and chitinase 10) and chitin deacetylase might play the crucial role in the degradation of insect chitin with their expressions significantly increased during the transition. Here, we also found that chitin synthase 1A might be involved in the chitin synthesis of cuticles during the metamorphosis in B. dorsalis.

Conclusions

Significant changes at transcriptional level were identified during the larval pupariation of B. dorsalis. Importantly, we also obtained a vast quantity of RNA-seq data and identified metamorphosis associated genes, which would all help us to better understand the molecular mechanism of metamorphosis process in B. dorsalis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5077-z) contains supplementary material, which is available to authorized users.

Genome-wide annotation of cuticular proteins in the oriental fruit fly (Bactrocera dorsalis), changes during pupariation and expression analysis of CPAP3 protein genes in response to environmental stresses

Cuticular proteins (CPs) are essential components of the insect cuticle as they create a structural and protective shield and may have a role in insect development. In this paper, we studied the CPs in the oriental fruit fly (Bactrocera dorsalis), one of the most economically important pests in the Tephritidae family around the world. The availability of a complete genome sequence (NCBI Assembly: ASM78921v2) allowed the identification of 164 CP genes in B. dorsalis. Comparative analysis of the CPs in B. dorsalis with those in the model insect Drosophila melanogaster and the closely related Ceratitis capitata, and CPs from mosquitoes, Lepidoptera, Hymenoptera and Coleoptera identified Diptera-specific genes and cuticle development patterns. Analysis of their evolutionary relationship revealed that some CP families had evolved according to the phylogeny of the different insect species, while others shared a closer relationship based on domain architecture. Subsequently, transcriptome analysis showed that while most of the CPs (60-100% of the family members) are expressed in the epidermis, some were also present in internal organs such as the fat body and the reproductive organs. Furthermore, the study of the expression profiles throughout development revealed a profound change in the expression of CPs during the formation of the puparium (pupariation). Further analysis of the expression profiles of the CPAP3 genes under various environmental stresses revealed them to be involved in the response to pesticides and arid and extreme temperatures conditions. In conclusion, the data provide a particular overview of CPs and their evolutionary and transcriptional dynamics, and in turn they lay a molecular foundation to explore their roles in the unique developmental process of insect metamorphosis and stress responses.

Tyrosine hydroxylase coordinates larval-pupal tanning and immunity in oriental fruit fly (Bactrocera dorsalis)

BACKGROUND

The oriental fruit fly Bactrocera dorsalis (Hendel), a notorious world pest infesting fruits and vegetables, has evolved a high level of resistance to many commonly used insecticides. In this study, we investigate whether tyrosine hydroxylase (TH) that is required for cuticle tanning (sclerotization and pigmentation) in many insects, could be a potential target in controlling B. dorsalis.

RESULTS

We cloned TH cDNA (BdTH) of B. dorsalis. The complete open reading frame of BdTH (KY911196) was 1737 bp in length, encoding a protein of 578 amino acids. Quantitative real-time PCR confirmed that BdTH was highly expressed in the epidermis of 3rd instar larvae, and its expression increased prior to pupation, suggesting a role in larval-pupal cuticle tanning. When we injected dsBdTH or 3-iodo-tyrosine (3-IT) as a TH inhibitor or fed insect diet supplemented with 3-IT, there was significant impairment of larval-pupal cuticle tanning and a severe obstacle to eclosion in adults followed by death in most. Furthermore, injection of Escherichia coli into larvae fed 3-IT resulted in 92% mortality and the expressions of four antimicrobial peptide genes were significantly downregulated.

CONCLUSION

These results suggest that BdTH might play a critical role in larval-pupal tanning and immunity of B. dorsalis, and could be used as a potential novel target for pest control. © 2017 Society of Chemical Industry.

De novo transcriptome of the mayfly Cloeon viridulum and transcriptional signatures of Prometabola

Mayflies (Ephemeroptera) display many primitive characters and a unique type of metamorphosis (Prometabola). However, information on the genomes and transcriptomes of this insect group is limited. The RNA sequencing study presented here generated the first de novo transcriptome assembly of Cloeon viridulum (Ephemeroptera: Baetidae), and compared gene expression signatures among the young larva (YL), mature larva (ML), subimago (SI), and imago (IM) stages of this mayfly. The transcriptome, based on 88 Gb of sequence data, comprised a set of 81,185 high quality transcripts. The number of differentially expressed genes (DEGs) in YL vs. ML, ML vs. SI, and SI vs. IM, was 4,825, 1,584, and 1,278, respectively, according to the reads per kilobase of transcript per million mapped reads analysis, assuming a false discovery rate <0.05 and a fold change >2. Gene enrichment analysis revealed that these DEGs were enriched in the "chitin metabolic process", "germ cell development", "steroid hormone biosynthesis", and "cutin, suberine, and wax biosynthesis" pathways. Finally, the expression pattern of a selected group of candidate signature genes for Prometabola, including vestigial, methoprene-tolerant, wingless, and broad-complex were confirmed by quantitative real time-PCR analysis. The Q-PCR analysis of larval, subimaginal, and imaginal stages of C. viridulum suggests that the development of mayflies more closely resembles hemimetamorphosis than holometamorphosis.

Cloning and structural characterization of juvenile hormone diol kinase in Spodoptera litura

Juvenile hormone (JH) is one of the key insect hormones that regulate metamorphosis. Juvenile hormone diol kinase (JHDK) is an enzyme involved in JH metabolism and catalyzes JH diol to form a polar end product, JH diol phosphate that has no JH activity. In this study, a JHDK complementary DNA (cDNA) was cloned from Spodoptera litura and the structure and expression of the gene was characterized. The cDNA was 714 base pairs in length and encoded a protein of 183 amino acids with a molecular mass of 21 kDa and an isoelectric point of 4.55. Based on the structure, three putative calcium binding motifs and guanosine triphosphate-binding motifs were predicted in the protein. Modeling of the 3-D structure showed that the protein consisted of eight α-helixes linked with loops, with no β-sheets. The gene was expressed in the epidermis, fat body and midgut of fifth and sixth instar larvae. The expression level in the epidermis was lower than in the fat body and midgut. The gene was expressed at higher levels at the early stages than in the later stages of fifth and sixth instar midgut and fat body. The results suggest that this gene may be involved in the regulation of the JH titer in larvae of S. litura.

Sequencing and de novo assembly of the Asian gypsy moth transcriptome using the Illumina platform

The Asian gypsy moth (Lymantria dispar) is a serious pest of forest and shade trees in many Asian and some European countries. However, there have been few studies of L. dispar genetic information and comprehensive genetic analyses of this species are needed in order to understand its genetic and metabolic sensitivities, such as the molting mechanism during larval development. In this study, high-throughput sequencing technology was used to sequence the transcriptome of the Asian subspecies of the gyspy moth, after which a comprehensive analysis of chitin metabolism was undertaken. We generated 37,750,380 high-quality reads and assembled them into contigs. A total of 37,098 unigenes were identified, of which 15,901 were annotated in the NCBI non-redundant protein database and 9,613 were annotated in the Swiss-Prot database. We mapped 4,329 unigenes onto 317 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. Chitin metabolism unigenes were found in the transcriptome and the data indicated that a variety of enzymes was involved in chitin catabolic and biosynthetic pathways.

De Novo Assembly and Developmental Transcriptome Analysis of the Small White Butterfly Pieris rapae

The small white butterfly Pieris rapae is one of the most destructive pests of Brassicaceae. Yet little is understood about its genes involved in development. To facilitate research on P. rapae, we sequenced the transcriptome of P. rapae during six developmental stages, including the egg, three larval stages, the pupa, and the adult. In total, 240 million high-quality reads were obtained. De novo assembly generated 96,069 unigenes with an average length of 1353 nt. Of these, 31,629 unigenes had homologs as determined by a blastx search against the NR database with a cut-off e-value of 10⁻⁵. Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted to functionally annotate those genes. Then, 849 genes involved in seven canonical development signaling pathway were identified, including dozens of key genes such as Hippo, Notch, and JAK2. A total of 21,883 differentially expressed (cut-off of 2-fold) unigenes were detected across the developmental stages, most of which were found between the egg and first larval stages. Interestingly, only 34 differentially expressed unigenes, most of which are cuticle protein related genes, were detected with a cut-off of 2¹⁰-fold. Furthermore, we identified 32 heat shock protein (Hsp) genes that were expressed with complete open reading frames. Based on phylogenetic trees of the Hsp genes, we found that Hsp genes with close evolutionary relationships had similar expression pattern. Additionally, partial pattern recognition receptors genes were found to be developmental regulated. This study provides comprehensive sequence resources for P. rapae and numerous differential expressed genes, and these findings will lay the foundation for future functional genomics studies on this species.

Dynamics and regulation of glycolysis-tricarboxylic acid metabolism in the midgut of Spodoptera litura during metamorphosis

Significant changes usually take place in the internal metabolism of insects during metamorphosis. The glycolysis-tricarboxylic acid (glycolysis-TCA) pathway is important for energy metabolism. To elucidate its dynamics, the mRNA levels of genes involved in this pathway were examined in the midgut of Spodoptera litura during metamorphosis, and the pyruvate content was quantified. The expression patterns of these genes in response to starvation were examined, and the interaction between protein phosphatase 1 (PP1) and phosphofructokinase (PFK) was studied. The results revealed that the expression or activities of most glycolytic enzymes was down-regulated in prepupae and then recovered in some degree in pupae, and all TCA-related genes were remarkably suppressed in both the prepupae and pupae. Pyruvate was enriched in the pupal midgut. Taken together, these results suggest that insects decrease both glycolysis and TCA in prepupae to save energy and then up-regulate glycolysis but down-regulate TCA in pupae to increase the supply of intermediates for construction of new organs. The expression of all these genes were down-regulated by starvation, indicating that non-feeding during metamorphosis may be a regulator of glycolysis-TCA pathway in the midgut. Importantly, interaction between PP1 and PFK was identified and is suggested to be involved in the regulation of glycolysis.

Proteomic analysis reveals that COP9 signalosome complex subunit 7A (CSN7A) is essential for the phase transition of migratory locust

The migratory locust displays a reversible, density-dependent transition between the two phases of gregaria and solitaria. This phenomenon is a typical kind of behavior plasticity. Here, we report that COP9 signalosome complex subunit 7A (CSN7A) is involved in the regulation of locust phase transition. Firstly, 90 proteins were identified to express differentially between the two phases by quantitative proteomic analysis. Gregaria revealed higher levels in proteins related to structure formation, melanism and energy metabolism, whereas solitaria had more abundant proteins related to digestion, absorption and chemical sensing. Subsequently, ten proteins including CSN7A were found to reveal differential mRNA expression profiles between the two phases. The CSN7A had higher mRNA level in the gregaria as compared with the solitaria, and the mRNA amount in the gregaria decreased remarkably during the 32 h-isolation. However, the mRNA level in the solitaria kept constant during the crowding rearing. Finally and importantly, RNA interference of CSN7A in gregaria resulted in obvious phase transition towards solitaria within 24 h. It suggests that CSN7A plays an essential role in the transition of gregaria towards solitaria in the migratory locust. To our knowledge, it's the first time to report the role of CSN in behavior plasticity of animals.

Transcriptome analysis of Bombyx mori larval midgut during persistent and pathogenic cytoplasmic polyhedrosis virus infection

Many insects can be persistently infected with viruses but do not show any obvious adverse effects with respect to physiology, development or reproduction. Here, Bombyx mori strain Daizo, persistently infected with cytoplasmic polyhedrosis virus (BmCPV), was used to study the host's transcriptional response after pathogenic infection with the same virus in midgut tissue of larvae persistently and pathogenically infected as 2nd and 4th instars. Next generation sequencing revealed that from 13,769 expressed genes, 167 were upregulated and 141 downregulated in both larval instars following pathogenic infection. Several genes that could possibly be involved in B. mori immune response against BmCPV or that may be induced by the virus in order to increase infectivity were identified, whereas classification of differentially expressed transcripts (confirmed by qRT-PCR) resulted in gene categories related to physical barriers, immune responses, proteolytic/metabolic enzymes, heat-shock proteins, hormonal signaling and uncharacterized proteins. Comparison of our data with the available literature (pathogenic infection of persistently vs. non-persistently infected larvae) unveiled various similarities of response in both cases, which suggests that pre-existing persistent infection does not affect in a major way the transcriptome response against pathogenic infection. To investigate the possible host's RNAi response against BmCPV challenge, the differential expression of RNAi-related genes and the accumulation of viral small RNAs (vsRNAs) were studied. During pathogenic infection, siRNA-like traces like the 2-fold up-regulation of the core RNAi genes Ago-2 and Dcr-2 as well as a peak of 20 nt small RNAs were observed. Interestingly, vsRNAs of the same size were detected at lower rates in persistently infected larvae. Collectively, our data provide an initial assessment of the relative significance of persistent infection of silkworm larvae on the host response following pathogenic infection with CPV, while they also highlight the relative importance of RNAi as an antiviral mechanism.

Expression, subcellular localization and protein-protein interaction of four isoforms of EcR/USP in the common cutworm

Ecdysone receptor (EcR) and ultraspiracle (USP) form heterodimers to mediate ecdysteroid signaling during molting and metamorphosis. Various EcR/USP heterodimers have been reported. However, it is unclear what kind of EcR/USP combination is adopted by lepidopteran insects during the larval-pupal metamorphosis and whether the EcR/USP heterodimer varies among different tissues. To address these questions, two isoforms of each EcR and USP were cloned from the common cutworm, their messenger RNA expression patterns were examined by real-time quantitative polymerase chain reaction in different tissues during the larval-pupal metamorphosis and in the midgut in response to hormonal induction. Furthermore, their subcellular localization and protein-protein interaction were explored by transient expression and far-western blotting, respectively. All the four genes were significantly up-regulated in prepuae and/or pupae. The expression profiles of EcRB1 and USP1 were nearly identical to each other in the epidermis, fat body and midgut, and a similar situation also applied to EcRA and USP2. The three genes responded to 20-hydroxyecdysone (20E) induction except for USP2, and USP1 could be up-regulated by both 20E and juvenile hormone. The four proteins mainly localized in the nucleus and the nuclear localization was promoted by 20E. The protein-protein interaction between each EcR and USP was found in vitro. These results suggest that two types of EcR/USP heterodimer (EcRA/USP2 and EcRB1/USP1) may exist simultaneously in the common cutworm, and the latter should play more important roles during the larval-pupal metamorphosis. In addition, the types of EcR/USP heterodimer do not vary in the tissues which undergo histolysis and regeneration during metamorphosis.

Protein profiles of Chinese white wax scale, Ericerus pela, at the male pupal stage by high-throughput proteomics

The Chinese white wax scale insect (Ericerus pela) is sexually dimorphic with holometabolous males and hemimetabolous females. Holometabolous insects were assumed to originate from hemimetabolous ancestors. Therefore, the male pupal stage is a major innovation compared with hemimetabolous female insects. Here, the protein profiles of the male pupae were obtained by high-throughput proteomics and analyzed using bioinformatics methods. A total of 1,437 peptides were identified and assigned to 677 protein groups. Most of the proteins had molecular weights below 40 kDa and isoelectric points from 4 to 7. Gene Ontology terms were assigned to 331 proteins, including metabolic process, developmental process, and cellular process. Kyoto Encyclopedia of Genes and Genomes annotations identified 142 pathways and most proteins were assigned to metabolism events. Pathways involved in cell growth and death, signal transduction, folding, and sorting and degradation were also identified. Six proteins that had undergone positive selection were classified into four groups, protein biosynthesis, protein degeneration, signal transduction, and detoxification. Many of the high-abundance proteins were enzymes involved in carbohydrate, lipid, and amino acid metabolism; signal transduction; degradation; and immunization, which indicated that metabolism, disruption, and development occurred intensely at the pupal stage. These processes are closely related to the physiological status of pupae. The results also suggested that these related proteins may be fundamental factors in the formation of pupae. This study describes pupal characterization at the molecular level and provides a basis for further physiological studies.

Function of the nuclear receptor FTZ-F1 during the pupal stage in Drosophila melanogaster

The nuclear receptor βFTZ-F1 is expressed in most cells in a temporally specific manner, and its expression is induced immediately after decline in ecdysteroid levels. This factor plays important roles during embryogenesis, larval ecdysis, and early metamorphic stages. However, little is known about the expression pattern, regulation and function of this receptor during the pupal stage. We analyzed the expression pattern and regulation of ftz-f1 during the pupal period, as well as the phenotypes of RNAi knockdown or mutant animals, to elucidate its function during this stage. Western blotting revealed that βFTZ-F1 is expressed at a high level during the late pupal stage, and this expression is dependent on decreasing ecdysteroid levels. By immunohistological analysis of the late pupal stage, FTZ-F1 was detected in the nuclei of most cells, but cytoplasmic localization was observed only in the oogonia and follicle cells of the ovary. Both the ftz-f1 genetic mutant and temporally specific ftz-f1 knockdown using RNAi during the pupal stage showed defects in eclosion and in the eye, the antennal segment, the wing and the leg, including bristle color and sclerosis. These results suggest that βFTZ-F1 is expressed in most cells at the late pupal stage, under the control of ecdysteroids and plays important roles during pupal development.

Physiological effect of mild thermal stress and its induction of gene expression in the common cutworm, Spodoptera litura

Heat shock protein (Hsp) and its cognate protein (Hsc) play important roles in helping insects survive extreme temperatures. However, high level of Hsp expression usually brings negative physiological effects on organisms. The mechanism of this trade-off is unclear. In this study, a lepidopteran insect, the common cutworm Spodoptera litura, was stressed at different temperatures, and the impact on both thermotolerance and fecundity was examined. The mRNA levels of four Hsp/Hscs (Hsp90, Hsc90, Hsp70 and Hsc70) and two ecdysone receptors (EcRs, EcRA and EcRB1) in different stresses and during the larval-pupal metamorphosis were determined. The results revealed that the pre-acclamation at mild stress increased the thermotolerance but decreased the egg production in adults. During the stress process, the mRNA levels of all the Hsp/Hsc and ecdysone receptor genes were significantly up-regulated. The two Hsp/Hsc70s and EcRs revealed consistent expression profiles with each other during the larval-pupal metamorphosis. Co-immunoprecipitation and Western blotting analysis indicated that Hsp/Hsc70 interacted with EcRs. RNAi of Hsc70 decreased the mRNA levels of two 20E-induced genes such as E74B and E75. Hsp70 transferred from the cytoplasm to nucleus in response to cold stress. These data together suggest that Hsp/Hsc70 might be involved in the regulation of 20E signaling, and the protein-protein interaction between Hsp/Hsc70 and EcRs probably act as a bridge mediating the trade-off between high thermotolerance and physiological defects.