Insight into the molecular mechanism of phosphine toxicity provided by functional analysis of cytochrome b5 fatty acid desaturase and dihydrolipoamide dehydrogenase in the red flour beetle, Tribolium castaneum.

Phosphine is the dominant chemical used in postharvest pest control. Widespread and highly frequent use of phosphine has been selected for pest insects, including Tribolium castaneum, which is highly resistant. Lipid peroxidation and reactive oxygen species (ROS) are two major factors determining phosphine toxicity; however, the mechanisms of production of these two factors in phosphine toxicity are still unknown. Here, we first determined the time course of phosphine-induced lipid peroxidation and ROS production in T. castaneum. Our results showed that lipid peroxidation occurs before ROS in the process of phosphine toxicity, and fumigated beetles with higher resistance levels were associated with weaker activity on lipid peroxidation and ROS. A significant decline in lipid peroxidation was observed in fumigated individuals after knockdown of cytochrome b5 fatty acid desaturase (Cyt-b5-r) via RNA interference (RNAi), indicating that Cyt-b5-r is critical for triggering phosphine-induced lipid peroxidation. Moreover, significant decreases in both ROS and mortality were detected in fumigated T. castaneum adults fed melatonin for 7 days, an inhibitor of lipid peroxidation. Cyt-b5-r RNAi also inhibited ROS production and mortality in phosphine-treated beetles. Meanwhile, a significant decrease in ROS production (68.4%) was detected in dihydrolipoamide dehydrogenase (DLD) knockdown individuals with phenotypes susceptible to phosphine, suggesting that lipid peroxidation initiates ROS with the expression of DLD. However, a significant increase in ROS (122.1%) was detected in the DLD knockdown beetles with strongly resistant phenotypes, indicating that the DLD-involved pathway may not be the only mechanism of ROS generation in phosphine toxicity and the existence of a moonlighting role in downregulating ROS in strongly resistant T. castaneum.

Identification of cuticular protein genes and analysis of their roles in phosphine resistance of the rusty grain beetle Cryptolestes ferrugineus.

The rusty grain beetle, Cryptolestes ferrugineus (Stephens) is one of the most economically important stored grain pests, and it has evolved the high resistance to phosphine. Cuticular proteins (CPs) are the major structural components of insect cuticle, and previous studies have confirmed that CPs were involved in insecticide resistance. However, the CPs of C. ferrugineus are still poorly characterized, and thus we conducted transcriptome-wide identification of CP genes and analyze their possible relationships with phosphine resistance in this pest. In this study, a total of 122 putative CPs were annotated in the C. ferrugineus transcriptome data by blasting with the known CPs of Tribolium castaneum. The analysis of conserved motifs revealed these CPs of C. ferrugineus belonging to 9 different families, including 87 CPR, 13 CPAP1, 7 CPAP3, 3 Tweedle, 1 CPLCA, 1 CPLCG, 5 CPLCP, 2 CPCFC, and 3 CPFL proteins. The further phylogenetic analysis showed the different evolutionary patterns of CPs. Namely, we found some CPs (CPR family) formed species-specific protein clusters, indicating these CPs might occur independently among insect taxa, and while some other CPs (CPAP1 and CPAP3 family) shared a closer correlation based on the architecture of protein domains. Subsequently, the previous RNA-seq data were applied to establish the expression profiles of CPs in a phosphine susceptible and resistant populations of C. ferrugineus, and a large amount of CP genes were found to be over-expressed in resistant insects. Lastly, an up-regulated CP gene (CPR family) was selected for the further functional analysis, and after this gene was silenced via RNA interference (RNAi), the sensitivity to phosphine was significantly enhanced in C. ferrugineus. In conclusion, the present results provided us an overview of C. ferrugineus CPs, and which suggested that the CPs might play the critical roles in phosphine resistance.

Amplification refractory mutation system based real-time PCR (ARMS-qPCR) for rapid resistance characterization of Tribolium castaneum to phosphine.

Resistance of Tribolium castaneum to phosphine is related to point mutations in DNA code corresponding to amino acid changes associated with a core metabolic enzyme dihydrolipoamide dehydrogenase (DLD), but the mutation patterns vary among different resistant populations. Thus, there is a great need to develop a cost-effective method to detect core mutations in T. castaneum, which would be the key factor to understand the molecular basis of phosphine resistance. Amplification refractory mutation system-based quantitative Real-Time PCR (ARMS-qPCR) is an ideal method that can rapidly detect point mutations. Here, the P45S and G131D mutations existed in the DLD of T. castaneum selected from strong Chinese resistance phenotypes, and the DLD P45S mutation, which represents a strong phosphine resistance allele, was confirmed as the most abundant mutation to determine strong resistance genotypes. Our study found that 85 out of 120 beetles carried the P45S resistance allele, including 51 homozygous and 34 heterozygous individuals. Moreover, there was a strong linear relationship (R2 = 0.917) between the resistance ratio and the resistance allele frequency among the strongly resistant populations. Our data showed that the ARMS-qPCR method that we developed could rapidly determine strong resistance phenotypes of T. castaneum to phosphine by detecting the DLD P45S mutation. These results not only provide a detailed example for developing an ARMS-qPCR-based method to characterize pesticide resistance, but also support further elucidation of the molecular basis of phosphine resistance.

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.

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.

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.

RNA interference of Argonaute-1 delays ovarian development in the oriental fruit fly, Bactrocera dorsalis (Hendel).

BACKGROUND: With the development of rapid resistance, new modes of action for pesticides are needed for insect control, such as RNAi-based biopesticides targeting essential genes. To explore the function of Argonaute-1 (Ago-1) and potential miRNAs in ovarian development of Bactrocera dorsalis, an important agricultural pest, and to develop a novel control strategy for the pest, BdAgo-1 was first identified in B. dorsalis. RESULTS: Spatiotemporal expression analysis indicated that BdAgo-1 had a relatively high transcriptional level in the ovarian tissues of adult female B. dorsalis during the sexual maturation period. RNA interference (RNAi) experiment showed that BdAgo-1 knockdown significantly decreased the expression levels of ovarian development-related genes and delayed ovarian development. Although RNAi-mediated silencing of Ago-1 led to a reduced ovary surface area, a subsequent oviposition assay revealed that the influence was minimal over a longer time period. Small RNA libraries were constructed and sequenced from different ovarian developmental stages of B. dorsalis adults. Among 161 identified miRNAs, 84 miRNAs were differentially expressed during the three developmental stages of the B. dorsalis ovary. BdAgo-1 silencing caused significant down-regulation of seven differentially expressed miRNAs (DEMs) showing relatively high expression levels (>1000 TPM (Transcripts per kilobase of exon model per million mapped reads)). The expression patterns of these seven core DEMs and their putative target genes were analyzed in the ovaries of B. dorsalis. CONCLUSION: The results indicate that Ago-1 and Ago-1-dependent miRNAs are indispensable for normal ovarian development in B. dorsalis and help identify miRNA targets useful for control of this pest.

Ovary-Specific Transcriptome and Essential Role of Nanos in Ovary Development in the Oriental Fruit Fly (Diptera: Tephritidae).

We used transcriptome analysis to research ovary development in Bactrocera dorsalis (Hendel). The ovary transcriptome of B. dorsalis yielded 66,463,710 clean reads that were assembled into 23,822 unigenes. After aligning to the Nr database in NCBI, 15,473 (64.95%) of the unigenes were matched to identified proteins. As determined by BLAST search, 11,043 (46.36%), 6,102 (25.61%), and 12,603 (52.90%) unigenes were each allocated to clusters via gene ontology, orthologous groups, and SwissProt, respectively. The Kyoto encyclopedia database of genes and genomes (KEGG) was further used to annotate these sequences, and 11,068 unigenes were mapped to 255 known pathways. Afterward, the genes that were possibly involved in oogenesis and ovary development were obtained from the transcriptome data and analyzed. Interestingly, seven ovary-specific genes were identified, including a Nanos gene that is involved in maintaining the primordial germ cells in many insects. Therefore, we further focused on the function of the BdNanos gene, and the gene was injected into B. dorsalis. As expected, the knocking down of Nanos gene expression led to significant inhibition of ovary development, suggesting an important role of this gene in the reproductive process of B. dorsalis. In summary, the present study provides an important reference for identifying the molecular mechanisms of oogenesis and ovary development in B. dorsalis. The BdNanos gene is crucial for ovary development in B. dorsalis and is therefore a potential new pest control target.

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.

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.

Molecular Characterization and Expression Profiling of Nuclear Receptor Gene Families in Oriental Fruit Fly, Bactrocera Dorsalis (Hendel).

The oriental fruit fly (Bactrocera dorsalis) is a pest that causes large economic losses in the fruit and vegetable industry, so its control is a major challenge. Nuclear receptors (NRs) are a superfamily of ligand-dependent transcription factors that directly combine with DNA to regulate the expression of downstream target genes. NRs are closely associated with multiple physiological processes such as metabolism, reproduction, and development. Through sequence searches and analysis, we identified 21 B. dorsalis NR genes, all of which contained at least one of the two characteristic binding domains. On the basis of the conserved sequences and phylogenetic relationships, we divided the 21 NR genes into seven subfamilies. All members of the NR0 subfamily and BdHR83, which belonged to the NR2E group, lacked ligand-binding domains. The BdDSF and BdHR51, which also belonged to the NR2Egroup, and BdE78 (which belonged to the NR1E group) all lacked DNA-binding domains. The BdDSF and BdHR83 sequences were incomplete, and were not successfully amplified. Development- and tissue-specific expression profiling demonstrated that the transcript levels of the 19 NR genes varied considerably among eggs, larva, pupae, and adults, as well as among larval and adult male and female tissues. Our results will contribute to a better understanding of NR evolution and expand our knowledge of B. dorsalis physiology.

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.

Potential targets for controlling Bactrocera dorsalis using cuticle- and hormone-related genes revealed by a developmental transcriptome analysis.

BACKGROUND: The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important agricultural pest and has developed resistance to many insecticides. To investigate vital genes participating in metamorphosis for development of additional control strategies, a comprehensive transcriptome analysis covering ten developmental stages of B. dorsalis was performed. RESULTS: There were 2132, 952, 1062, 2301 and 1333 differentially expressed genes identified during hatching, 1st-instar larval molting, 2nd-instar larval molting, pupariation and emergence, respectively. Further expression analyses indicated that genes in hormone- (20-hydroxyecdysone and juvenile hormone) and cuticle- (chitin and cuticle protein) related pathways were essential for metamorphosis in B. dorsalis. Among chitinase (Cht) genes, BdCht-5, -8 and -10 were differentially expressed during larval-larval, larval-pupal and pupal-adult moltings. However, BdCht7 was differentially expressed during egg-larval and larval-larval moltings. Knockdown of BdCht7 at the 1st-instar larval stage disrupted normal development of larvae and was lethal to B. dorsalis. Among cuticle protein (CP) genes, 15 genes (BdCPLCG-1, BdCPLCP-2, BdCPAP1-B2, BdRR1-21, BdRR1-31, BdRR2-15, BdRR2-26, BdRR2-30, BdRR2-32, BdTweedle-9, BdTweedle-24, BdRR2-10, BdCPAP3-C1, BdRR1-34 and BdRR1-41) were differentially expressed during four of five types of moltings. Among hormone-relative genes, BdJHBP-4, -9 and -13 were differentially expressed during all five types of moltings, whereas BdJHBP-5, -12 and BdHR4 were differentially expressed during four of five types of moltings. CONCLUSION: This study reveals critical genes involved in development and metamorphosis of B. dorsaslis, and BdCht7 is dispensable for larval survival. It also provides comprehensive transcriptome information for finding more molecular targets to control this pest. © 2020 Society of Chemical Industry.

Assessment of Bactrocera dorsalis (Diptera: Tephritidae) Diets on Adult Fecundity and Larval Development: Insights Into Employing the Sterile Insect Technique.

Bactrocera dorsalis (Hendel) is a notorious insect pest that attacks diverse vegetables and fruits worldwide. The sterile insect technique has been developed as an environmentally friendly and effective control method that depends on the mass production of target flies. Because dietary yeast (protein) and sucrose (carbohydrate) are important in adult diets, yeast:sucrose (Y:S) mixtures are crucial for the mass-rearing of B. dorsalis. In this study, we found adult diets with different ratios of yeast to sucrose-influenced fecundity, and an extremely high or low Y:S ratios significantly decreased egg production of B. dorsalis. Additionally, the maximum oviposition efficiency was realized at dietary yeast to sucrose ratios of 1:1 and 1:3, suggesting their potential use to produce more eggs for the mass production of B. dorsalis. Here, new gel diets having different yeast concentrations (g/L water) were also assessed for rearing B. dorsalis larvae. Gel diets containing 20 g/L yeast led to a higher pupation, pupal weight and adult eclosion rate, and a shorter developmental time than other yeast concentrations. Moreover, the present gel diet also resulted in greater pupal production and adult emergence rates than previously used liquid and solid artificial diets, revealing that it is suitable for rearing B. dorsalis larvae. This research provides a useful reference on artificial diets mixtures for mass rearing B. dorsalis, which is critical for employing the sterile insect technique.

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.

Transcriptome Profiling of the Whitefly Bemisia tabaci MED in Response to Single Infection of Tomato yellow leaf curl virus, Tomato chlorosis virus, and Their Co-infection.

Tomato yellow leaf curl virus (TYLCV) and Tomato chlorosis virus (ToCV) are two of the most devastating cultivated tomato viruses, causing significant crop losses worldwide. As the vector of both TYLCV and ToCV, the whitefly Bemisia tabaci Mediterranean (MED) is mainly responsible for the rapid spread and mixed infection of TYLCV and ToCV in China. However, little is known concerning B. tabaci MED's molecular response to TYLCV and ToCV infection or their co-infection. We determined the transcriptional responses of the whitefly MED to TYLCV infection, ToCV infection, and TYLCV&ToCV co-infection using Illumina sequencing. In all, 78, 221, and 60 differentially expressed genes (DEGs) were identified in TYLCV-infected, ToCV-infected, and TYLCV&ToCV co-infected whiteflies, respectively, compared with non-viruliferous whiteflies. Differentially regulated genes were sorted according to their roles in detoxification, stress response, immune response, transport, primary metabolism, cell function, and total fitness in whiteflies after feeding on virus-infected tomato plants. Alterations in the transcription profiles of genes involved in transport and energy metabolism occurred between TYLCV&ToCV co-infection and single infection with TYLCV or ToCV; this may be associated with the adaptation of the insect vector upon co-infection of the two viruses. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses demonstrated that the single infection with TYLCV or ToCV and the TYLCV&ToCV co-infection could perturb metabolic processes and metabolic pathways. Taken together, our results provide basis for further exploration of the molecular mechanisms of the response to TYLCV, ToCV single infection, and TYLCV&ToCV co-infection in B. tabaci MED, which will add to our knowledge of the interactions between plant viruses and insect vectors.

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.

Corazonin Signaling Is Required in the Male for Sperm Transfer in the Oriental Fruit Fly Bactrocera dorsalis.

Corazonin (Crz) is a widely distributed neuropeptide (or neurohormone) in insects with diverse physiological functions. The present study aimed to reveal the functions of Crz and its receptor (CrzR) in the regulation of sexual behavior and fertility in male Bactrocera dorsalis. Tissue-specific expression analyses showed that the BdCrz transcript was most abundant in the central nervous system (CNS), and the BdCrzR transcript was most abundant in both the fat body and CNS. Immunochemical localization confirmed that three pairs of Crz-immunoreactive neurons are located in the dorsolateral protocerebrum region of male adult brain. Importantly, RNAi-mediated Crz knockdown lengthened mating duration in males, and knockdown of Crz or CrzR strongly decreased male fertility in the following 3 days, while the courtship behavior and mating efficiency were not affected. The reduced number of sperm in the reproductive organs of mated females indicated that Crz knockdown in males reduced sperm transfer. The findings of this study indicate that Crz contributes to the reproductive physiology of the oriental fruit fly B. dorsalis by regulating sperm transfer in male adults.

Involvement of Met and Kr-h1 in JH-Mediated Reproduction of Female Bactrocera dorsalis (Hendel).

Juvenile hormone (JH) prevents metamorphosis during insect larval stages and promotes adult reproductive processes. Krüppel-homolog 1 (Kr-h1), a zinc finger transcription factor assumed to be induced by JH via the JH receptor methoprene-tolerant (Met), mediates the antimetamorphic effect of JH in insects, but its function in JH-mediated reproductive processes has not been fully explored. In this study, Met and Kr-h1 involved in the JH signaling pathway were first cloned and identified from the oriental fruit fly, Bactrocera dorsalis, an important pest infesting fruit and vegetables worldwide. Subsequent spatiotemporal expression analysis revealed that Met and Kr-h1 were both highly expressed in 7-day-old adults and fat body of female adults, respectively. Treatment with a JH analog (methoprene) significantly induced the expression of JH signaling and vitellogenin (Vg) genes and accelerated ovary development. RNA interference (RNAi) further revealed that either Met or Kr-h1 depletion at the adult stage of B. dorsalis impeded ovary development, with significantly lower egg production noted as well. In addition, rescue through methoprene application after RNAi stimulated the expression of JH signaling and Vg genes. Although there were still differences in ovary phenotype between rescued insects and the pre-RNAi control, ovary redevelopment with a larger surface area was observed, consistent with the spatiotemporal expression and phenotypes recorded in the original methoprene experiment. Our data reveal the involvement of Met and Kr-h1 in insect vitellogenesis and egg production, thus indicating the crucial role of the JH signaling pathway in insect reproduction.

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.