In silico cloning and annotation of genes involved in the digestion, detoxification and RNA interference mechanism in the midgut of Bactrocera dorsalis [Hendel (Diptera: Tephritidae)]

CYP4G100 contributes to desiccation resistance by mediating cuticular hydrocarbon synthesis in Bactrocera dorsalis

The oriental fruit fly Bactrocera dorsalis (Hendel) is expanding its distribution to higher latitudes. Our goal in this study was to understand how B. dorsalis adapts to higher latitude environments that are more arid than tropical regions. Cuticular hydrocarbons (CHCs) on the surface of the epicuticle in insects act as a hydrophobic barrier against water loss. The essential decarbonylation reaction in CHC synthesis is catalysed by CYP4G, a cytochrome P450 subfamily protein. Hence, in B. dorsalis it is necessary to clarify the function of the CYP4G gene and its role in desiccation resistance. CYP4G100 was identified in the B. dorsalis genome. The complete open reading frame (ORF) encodes a CYP4 family protein (552 amino acid residues) that has the CYP4G-specific insertion. This CYP4G gene was highly expressed in adults, especially in the oenocyte-rich peripheral fat body. The gene can be induced by desiccation treatment, suggesting its role in CHC synthesis and waterproofing. Silencing of CYP4G100 resulted in a decrease of CHC levels and the accumulation of triglycerides. It also increased water loss and resulted in higher desiccation susceptibility. CYP4G100 is involved in hydrocarbon synthesis and contributes to cuticle waterproofing to help B. dorsalis resist desiccation in arid environments.

Sublethal Dose of β-Cypermethrin Impairs the Olfaction of Bactrocera dorsalis by Suppressing the Expression of Chemosensory Genes

The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive fruit insect pests. β-cypermethrin has been widely used in the orchard to control this major insect. Based on the resistance monitoring in 2011, B. dorsalis developed significant resistance against β-cypermethrin in fields. This indicated that the B. dorsalis has been exposed to sublethal concentrations of β-cypermethrin in the field for a long time. Thus, it is urgent to understand the sublethal effects of β-cypermethrin on this fly to guide the rational use of an insecticide. According to the olfactory preference assays and electroantennogram (EAG) recording, the B. dorsalis after β-cypermethrin exposure (LD₃₀ = 10 ng/fly) severely decreased the ability to perceive the tested odorants. Moreover, we then performed quantitative real-time PCR and found the chemosensory genes including odorant receptor co-receptor (BdorORco) and ionotropic receptor co-receptors (BdorIRcos) were obviously suppressed. Our results demonstrated that the sublethal dose of β-cypermethrin impairs the olfaction of the pest insects by suppressing the expression of chemosensory genes (BdorORco and BdorIRcos), which expanded our knowledge of the sublethal effects of the pesticide on insects.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

Overexpression of PxαE14 Contributing to Detoxification of Multiple Insecticides in Plutella xylostella (L.)

The diamondback moth, Plutella xylostella (L.), has evolved with varying degrees of resistance to almost all major classes of insecticides and has become the most resistant pest worldwide. The multiresistance to different types of insecticides has been frequently reported in P. xylostella, but little is known about the mechanism. In this study, a carboxylesterase (CarE) gene, PxαE14, was found significantly overexpressed in a field-evolved multiresistant P. xylostella population and can be dramatically induced by eight of nine tested insecticides. Results of the real-time quantitative polymerase chain reaction (RT-qPCR) showed that PxαE14 was predominantly expressed in the midgut and malpighian tubule of larvae. Knockdown of PxαE14 dramatically increased the susceptibility of the larvae to β-cypermethrin, bifenthrin, chlorpyrifos, fenvalerate, malathion, and phoxim, while overexpression of PxαE14 in Drosophila melanogaster increased the tolerance of the fruit flies to these insecticides obviously. More importantly, gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) assay showed that the recombinant PxαE14 expressed in Escherichia coli exhibited metabolic activity against the six insecticides. The homology modeling, molecular docking, and molecular dynamics simulation analyses showed that these six insecticides could stably bind to PxαE14. Taken together, these results demonstrate that constitutive and inductive overexpression of PxαE14 contributes to detoxification of multiple insecticides involved in multiresistance in P. xylostella. Our findings provide evidence for understanding the molecular mechanisms underlying the multiresistance in insect pests.

C/EBPα Regulates PxTreh1 and PxTreh2 Trehalase-Related Bt Resistance in Plutella xylostella (L.)

Simple Summary

The diamondback moth (Plutella xylostella) is a major agricultural pest of cruciferous vegetables and crops worldwide, causing economic losses of up to USD 5 billion annually. The long-term use of insecticides leads to the rapid evolution of resistance in insects, which greatly increases the difficulty of controlling pests. Trehalase regulates energy metabolism in insects by converting trehalose into two glucose molecules. The existence of trehalase is critical for insect flight and larval stress resistance. However, whether trehalase participates in the development of pesticide resistance remains unclear. In this study, we found that the activity of trehalase and the levels of gene expression in Bt-resistant and field populations of P. xylostella were significantly higher than they were in the susceptible strains. By analyzing the promoter sequences of PxTreh1 and PxTreh2, we confirmed the interaction between C/EBPα and the PxTreh2 promoter. The findings of this study suggest that C/EBPα mediates the adaptability of P. xylostella to adverse environmental stressors by regulating the expression of trehalase.

Abstract

Trehalase regulates energy metabolism in insects by converting trehalose into two glucose molecules. High amounts of trehalase are critical for insect flight and larval stress resistance. However, whether trehalase participates in the development of pesticide resistance remains unclear. In this study, we explored this phenomenon and the mechanism that underlies the regulation of Trehalase transcription. We found that overexpression of PxTreh1 and PxTreh2 induced Bacillus thuringiensis (Bt) resistance in Plutella xylostella. The promoter sequences of PxTreh1 and PxTreh2 were also cloned and identified. The dual-luciferase reporter system and RNA interference technology revealed that the expression of PxTreh1 and PxTreh2 genes is possibly regulated by the CCAAT enhancer-binding protein (C/EBPα). A yeast one-hybrid experiment confirmed the interaction between C/EBPα and the PxTreh2 promoter. The findings of this study suggest that C/EBPα mediates the adaptability of P. xylostella to adverse environmental stressors by regulating the expression of trehalase.

20E and MAPK signal pathway involved in the effect of reproduction caused by cyantraniliprole in Bactrocera dorsalis Hendel (Diptera: Tephritidae)

BACKGROUND

It is a common phenomenon that insecticides affect insect reproduction and insect hormones. After cyantraniliprole treatment, the egg production and remating behavior of female Bactrocera dorsalis were affected, a phenomenon of 'hormesis' appeared, but the change at the molecular level was unknown. Therefore, we investigated the fertility, insect hormone titers and transcription levels and used RNAi to prove the function of genes, to explore the molecular mechanism of cyantraniliprole causing reproductive changes in female B. dorsalis.

RESULTS

LC₂₀ treatment promoted egg production, while LC₅₀ treatment inhibited it. Both high and low concentrations inhibited female ovaries' development and reduced the length of the ovarian tubes. Among insect hormones, only the titer of 20-hydroxyecdysone (20E) changed significantly. According to the KEGG pathway enrichment analysis of RNA-seq, there are significant differences in insect hormone synthesis and MAPK signal pathways between treatments. Furthermore, 20E biosynthetic genes, BdVgs and BdVgR were all down-regulated, and multiple MAPK signaling pathway genes were up-regulated. Based on qRT-PCR, the expression of BdCyp307A1, BdCyp302A1, BdMEKK4 and BdMAP2K6 within 1-11 days after treatment were consistent with the change of 20E titer. The BdVg1 and BdVg2 in LC₅₀ were still suppressed, while the LC₂₀ returned to normal in 9-11 days. RNAi indicated that BdMEKK4 and BdMAP2K6 participated in the transcriptional regulation of BdCyp307A1 and BdCyp302A1, then affected the levels of BdVgs.

CONCLUSION

Cyantraniliprole affected 20E through MAPK signal pathway, causing many genes to be down-regulated during the early period but up-regulated during the late period, ultimately affecting the reproduction of B. dorsalis. © 2021 Society of Chemical Industry.

The short neuropeptide F receptor regulates olfaction-mediated foraging behavior in the oriental fruit fly Bactrocera dorsalis (Hendel)

The short neuropeptide F (sNPF) signaling system, consisting of sNPF and its receptor (sNPFR), influences many physiological processes in insects, including feeding, growth and olfactory memory. We previously showed that sNPF regulates olfactory sensitivity in the oriental fruit fly Bactrocera dorsalis (Hendel) during starvation. However, the functional analysis of sNPFR is constrained by the failure of RNA interference in this species. Here, we generated a null sNPFR mutant using the CRISPR/Cas9 system to investigate the physiological roles of this receptor in more detail. G0 adults were produced at a frequency of 60.8%, and sNPFR^-/- mutants were obtained after several generations of backcrossing followed by self-crossing among heterozygous flies. We found that the mutants were significantly less successful at foraging for certain foods and showed increased foraging latency. Electroantennogram (EAG) assays indicated that the mutants had significantly lower electrophysiological responses to three tested odorants. Furthermore, qPCR data revealed the inhibition of several olfactory receptor genes, including Orco. Immunohistochemistry showed that BdsNPFR was localized in cells under the sensillum on the antennae. Based on their shape and size, the BdsNPFR⁺ cells differ from odorant receptor neurons (ORNs), which were labeled using a Drosophila melanogaster Orco antibody. Our data suggest that sNPFR regulates olfaction-mediated foraging behavior by mediating interactions between BdsNPFR⁺ cells and selected ORNs.

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.

Comparative transcriptomic analysis reveals female-biased olfactory genes potentially involved in plant volatile-mediated oviposition behavior of Bactrocera dorsalis

Background

Olfactory systems take on important tasks to distinguish salient information from a complex olfactory environment, such as locating hosts, mating, aggression, selecting oviposition sites, and avoiding predators. The olfactory system of an adult insect consists of two pairs of main olfactory appendages on the head, the antennae, and the palps, which are covered with sensilla. Benzothiazole and 1-octen-3-ol could elicit oviposition behavior in gravid B. dorsalis are regarded as oviposition stimulants. However, the mechanism for how B. dorsalis percepts benzothiazole and 1-octen-3-ol still remains unknown.

Results

We conducted a comparative analysis of the antennal transcriptomes in different genders of B. dorsalis using Illumina RNA sequencing (RNA-seq). We identified a total of 1571 differentially expressed genes (DEGs) among the two sexes, including 450 female-biased genes and 1121 male-biased genes. Among these DEGs, we screened out 24 olfaction-related genes and validated them by qRT-PCR. The expression patterns of these genes in different body parts were further determined. In addition, we detected the expression profiles of the screened female-biased chemosensory genes in virgin and mated female flies. Furthermore, the oviposition stimulants-induced expression profilings were used to identify chemosensory genes potentially responsible for benzothiazole and 1-octen-3-ol perception in this fly.

Conclusions

The results from this study provided fundamental data of chemosensory DEGs in the B. dorsalis antenna. The odorant exposure assays we employed lay a solid foundation for the further research regarding the molecular mechanism of benzothiazole and 1-octen-3-ol mediated oviposition behavior in B. dorsalis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07325-z.

Crustacean cardioactive peptide (CCAP) of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae): Molecular characterization, distribution and its potential roles in larva-pupa ecdysis

Insects must undergo ecdysis for successful development and growth, and the crustacean cardioactive peptide (CCAP) is one of the most important hormone in this process. Here we reported a cDNA encoding for the CCAP precursor cloned from the oriental fruit fly, Bactrocera dorsalis, a most destructive insect pest of agriculture. The CCAP mature peptide (PFCNAFTGC-NH2) of B. dorsalis was generated by post-translational processing and found to be highly comparable with other insects. RT-qPCR showed that mRNA of CCAP in B. dorsalis (BdCCAP) was predominantly expressed in the central nervous system (CNS) and midgut of 3rd-instar larvae. By using immunohistochemical analysis, we also localized the endocrine cells that produce CCAP in the CNS, ring gland and midgut of 3rd-instar larvae of B. dorsalis. The synthetic CCAP mature peptide could induce the expression of mRNA of adipokinetic hormone (AKH), the metabolic neuropeptides in insects. The expression of BdCCAP mRNA in the CNS, but not in the midgut, could be upregulated in the response to the challenge of insect molting hormone, 20-hydroxyecdysone.

First Insights into the Intrapuparial Development of Bactrocera dorsalis (Hendel): Application in Predicting Emergence Time for Tephritid Fly Control

Intrapuparial development is a special pattern of metamorphosis in cyclorrhaphous flies, in which the pupa forms in an opaque, barrel-like puparium. This has been well studied in forensic insects for age estimations. In this study, the intrapuparial development of a quarantine agricultural pest, Bactrocera dorsalis (Hendel), was studied under a constant temperature of 27 ± 1 °C and 70 ± 5% relative humidity. Results showed that intrapuparial development could be divided into five stages: Larval-pupal apolysis, cryptocephalic pupa, phanerocephalic pupa, pharate adult, and emergent adult. It lays a morphology-based foundation for molecular mechanism studies and enhances the understanding of the physiological basis for changes in intrapuparial development. More importantly, the chronology of intrapuparial development can be used to predict the emergence time of tephritid flies, indicating when to spray insecticides to control these phytophagous agricultural pests. This may be an effective approach to reduce the use of insecticides and slow down the evolution of insecticidal resistance.

The South American Fruit Fly: An Important Pest Insect With RNAi-Sensitive Larval Stages

RNA interference (RNAi) technology has been used in the development of approaches for pest control. The presence of some essential genes, the so-called “core genes,” in the RNAi machinery is crucial for its efficiency and robust response in gene silencing. Thus, our study was designed to examine whether the RNAi machinery is functional in the South American (SA) fruit fly Anastrepha fraterculus (Diptera: Tephritidae) and whether the sensitivity to the uptake of double-stranded RNA (dsRNA) could generate an RNAi response in this fruit fly species. To prepare a transcriptome database of the SA fruit fly, total RNA was extracted from all the life stages for later cDNA synthesis and Illumina sequencing. After the de novo transcriptome assembly and gene annotation, the transcriptome was screened for RNAi pathway genes, as well as the duplication or loss of genes and novel target genes to dsRNA delivery bioassays. The dsRNA delivery assay by soaking was performed in larvae to evaluate the gene-silencing of V-ATPase, and the upregulation of Dicer-2 and Argonaute-2 after dsRNA delivery was analyzed to verify the activation of siRNAi machinery. We tested the stability of dsRNA using dsGFP with an in vitro incubation of larvae body fluid (hemolymph). We identified 55 genes related to the RNAi machinery with duplication and loss for some genes and selected 143 different target genes related to biological processes involved in post-embryonic growth/development and reproduction of A. fraterculus. Larvae soaked in dsRNA (dsV-ATPase) solution showed a strong knockdown of V-ATPase after 48 h, and the expression of Dicer-2 and Argonaute-2 responded with an increase upon the exposure to dsRNA. Our data demonstrated the existence of a functional RNAi machinery in the SA fruit fly, and we present an easy and robust physiological bioassay with the larval stages that can further be used for screening of target genes at in vivo organisms’ level for RNAi-based control of fruit fly pests. This is the first study that provides evidence of a functional siRNA machinery in the SA fruit fly.

Characterization and Function of Two Short Peptidoglycan Recognition Proteins Involved in the Immunity of Bactrocera dorsalis (Hendel)

Peptidoglycans (PGNs) are major bacterial components recognized by the immune systems of insects and mammals. PGN recognition proteins (PGRPs) are widely distributed and highly conserved in vertebrates and invertebrates. PGRPs are a family of pattern recognition receptors that recognize peptidoglycan and regulate immune responses. In this study, we cloned two PGRP genes (BdPGRP-SA and BdPGRP-SD) from Bactrocera dorsalis (Hendel), which encode 192 and 196 amino acid residues, respectively. Both genes were highly expressed in adults, especially in the fat body and midgut. These two genes were up-regulated when challenged by the immune triggers, PGN-EB (Escherichia coli O111:B4) and PGN-SA (Staphylococcus aureus). The suppression of transcriptional expression of either gene by RNA interference (RNAi) resulted in increased sensitivities to Gram-negative E. coli and Gram-positive S. aureus PGNs. Suppression of BdPGRP-SA and -SD expression by RNAi resulted in weak expressions of four antimicrobial peptides (AMPs) upon injected with E. coli or S. aureus. BdPGRP-SA and -SD are involved in recognizing both Gram-negative and Gram-positive bacteria independently to activate the downstream AMP’s response to bacterial infection.

The Ecdysis Triggering Hormone System, via ETH/ETHR-B, Is Essential for Successful Reproduction of a Major Pest Insect, Bactrocera dorsalis (Hendel)

Ecdysis triggering hormone (ETH), released by the Inka cells, is a master hormone in regulating the ecdysis process in insect. Here we investigated the presence and role of the ETH signaling in the female adult of the oriental fruit fly, Bactrocera dorsalis (Hendel) that is one of the most important invasive pest insects in agriculture worldwide. In the female adult, ETH was confirmed in the Inka cells at the tracheae by immunostaining and also in vitro exposure to ETH stimulated the isolated corpora allata of adult in activity. Then we prepared cDNA of females at 0, 5, 10, 15, and 20 days after adult eclosion, and RT-qPCR showed that the expression pattern of ETH and its receptor ETHR-B started from a peak at the day of adult eclosion (day 0), then dropped to basal levels and increased again between day 10 and 15 which is also the period corresponding to ovary growth. In contrast, ETHR-A was absent with Ct values of >33. The expression patterns of the ecdysteroid-producing Halloween genes Spook and Shade, and the vitellogenin genes Vg1, Vg2, and Vg3 co-occurred with peak levels at days 10–15, and also juvenile hormone acid methyltransferase (JHAMT) showed increased levels on day 15. Further in RNAi assays to better understand the role of ETH and ETHR, dsRNA was injected to adult and this led to a respective decrease in expression of 62 and 56% for ETH and ETHR-B, while ETHR-A stayed absent with Ct values of 33. In these RNAi-females, there was an apparently decreased expression for JHAMT and Vg2, together with a significant decrease of the JH titer and egg production. Injection of the JH mimetic methoprene could rescue Vg2 expression and egg production. Upstream, in dsETH/dsETHR-injected females, 20-hydroxyecdysone (20E) injection rescued the transcriptions of ETH and ETHR and also egg production. In summary, our results shed more light on the pivotal role that the ETH peptide hormone and its receptor ETHR-B play an essential role in the reproduction of the female adult of B. dorsalis, via the regulation of JH and vitellogenin, which are controlled by a pulse of 20E.

RNA-sequencing of a citrus bud-feeder, Podagricomela weisei (Coleoptera: Chrysomelidae), reveals xenobiotic metabolism/core RNAi machinery-associated genes and conserved miRNAs

The citrus leaf-mining beetle, Podagricomela weisei Heikertinger, is an important citrus pest that ingests the mesophyll and new shoots. The mechanism underlying the xenobiotic metabolism of P. weisei is not well understood, in part because of a lack of available genomic and transcriptomic data, which has hampered the development of novel pest management approaches [e.g., RNA interference (RNAi)]. In this study, we completed the deep sequencing of the P. weisei transcriptome to identify factors potentially involved in xenobiotic metabolism and the core RNAi machinery. The sequencing of the P. weisei transcriptome generated >27 million clean reads, ultimately yielding 90,410 unigenes with an N50 of 1065 bp. The unigenes were used as queries to search the Nr database, which revealed that 21,847 unigenes were homologous to known genes in various species. Transcripts encoding genes involved in xenobiotic metabolism were identified, including genes encoding cytochrome P450 monooxygenase (P450, 47 unigenes), glutathione S-transferase (GST, 12 unigenes), esterase (EST, 25 unigenes), and the ATP-binding cassette transporter (ABC transporter, 32 unigenes). A parallel sequencing of small RNAs detected 30 conserved miRNAs, with the most abundant being Pwe-miR-1-3p, with an expression level reaching 517,996 reads in the prepared library, followed by Pwe-miR-8-3p (149,402 reads). Genes encoding components of the miRNA, siRNA, and piRNA pathways were also identified, and the results indicated that P. weisei possesses only one of each gene in all three pathways. In summary, this is the first detailed analysis of the transcriptome and small RNAs of P. weisei. The datasets presented herein may form the basis for future molecular characterizations of P. weisei as well as the development of enhanced pest control strategies.

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.

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.

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.

Carboxylesterase genes in pyrethroid resistant house flies, Musca domestica

Carboxylesterases are one of the major enzyme families involved in the detoxification of pyrethroids. Up-regulation of carboxylesterase genes is thought to be a major component of insecticide resistant mechanisms in insects. Based on the house fly transcriptome and genome database, a total of 39 carboxylesterase genes of different functional clades have been identified in house flies. In this study, eleven of these genes were found to be significantly overexpressed in the resistant ALHF house fly strain compared with susceptible aabys and wild-type CS strains. Eight up-regulated carboxylesterase genes with their expression levels were further induced to a higher level in response to permethrin treatments, indicating that constitutive and inductive overexpression of carboxylesterases are co-responsible for the enhanced detoxification of insecticides. Spatial expression studies revealed these up-regulated genes to be abundantly distributed in fat bodies and genetically mapped on autosome 2 or 3 of house flies, and their expression could be regulated by factors on autosome 1, 2 and 5. Taken together, these results demonstrate that multiple carboxylesterase genes are co-upregulated in resistant house flies, providing further evidence for their involvement in the detoxification of insecticides and development of insecticide resistance.

The neuropeptides and protein hormones of the agricultural pest fruit fly Bactrocera dorsalis: What do we learn from the genome sequencing and tissue-specific transcriptomes?

Neuropeptides and protein hormones are very important signaling molecules, and are involved in the regulation and coordination of various physiological processes in invertebrates and vertebrates. Using a bioinformatics approach, we screened the recently sequenced genome and six tissue-specific transcriptome databases (central nervous system, fat body, ovary, testes, male accessory glands, antennae) of the oriental fruit fly (Bactrocera dorsalis) that is economically one of the most important pest insects of tropical and subtropical fruit. Thirty-nine candidate genes were found to encode neuropeptides or protein hormones. These include most of the known insect neuropeptides and protein hormones, with the exception of adipokinetic hormone-corazonin-related peptide, allatropin, diuretic hormone 34, diuretic hormone 45, IMFamide, inotocin, and sex peptide. Our results showed the neuropeptides and protein hormones of Diptera insects appear to have a reduced repertoire compared to some other insects. Moreover, there are also differences between B. dorsalis and the super-model of Drosophila melanogaster. Interesting features of the oriental fruit fly are the absence of genes coding for sex peptide and the presence of neuroparsin and two genes coding neuropeptide F. The majority of the identified neuropeptides and protein hormones is present in the central nervous system, with only a limited number of these in the other tissues. Moreover, we predicted their physiological functions via comparing with data of FlyBase and FlyAtlas. Taken together, owing to the large number of identified peptides, this study can be used as a reference about structure, tissue distribution and physiological functions for comparative studies in other model and important pest insects.

Functional characterization of BdB1, a well-conserved carboxylesterase among tephritid fruit flies associated with malathion resistance in Bactrocera dorsalis (Hendel)

There are many evidences that insect carboxylesterase possess important physiological roles in xenobiotic metabolism and are implicated in the detoxification of organophosphate (OP) insecticides. Despite the ongoing resistance development in the oriental fruit fly, Bactrocera dorsalis (Hendel), the molecular basis of carboxylesterase and its ability to confer OP resistance remain largely obscure. This study was initiated to provide a better understanding of carboxylesterase-mediated resistance mechanism in a tephritid pest fly. Here, we narrow this research gap by demonstrating a well-conserved esterase B1 gene, BdB1, mediates malathion resistance development via gene upregulation with the use of a laboratory selected malathion-resistant strain (MR) of B. dorsalis. No sequence mutation of BdB1 was detected between MR and the susceptible strain (MS) of B. dorsalis. BdB1 is predominantly expressed in the midgut, a key insect tissue for detoxification. As compared with transcripts in MS, BdB1 was significantly more abundant in multiple tissues in the MR. RNA interference (RNAi)-mediated knockdown of BdB1 significantly increased malathion susceptibility. Furthermore, heterologous expression along with cytotoxicity assay revealed BdB1 could probably have the function of malathion detoxification.

Ecdysis Triggering Hormone Signaling (ETH/ETHR-A) Is Required for the Larva-Larva Ecdysis in Bactrocera dorsalis (Diptera: Tephritidae)

Insects must undergo ecdysis for successful development and growth, and the ecdysis triggering hormone (ETH), released by the Inka cells, is a master hormone in this process. In this study, we determined the sequence of the ETH precursor and receptors in an agriculturally important pest insect, the oriental fruit fly Bactrocera dorsalis (Hendel). We identified two functionally distinct splice receptor isoforms: BdETH-R-A and BdETH-R-B, and when expressed in Chinese hamster ovary (CHO-WTA11) cells, they exhibited a high sensitivity to the two mature peptides BdETH1 and BdETH2. The BdETH transcript was detected in the tracheal tissue of the larvae. Inka cells were identified with immunohistochemical antibody staining against Drosophila melanogaster ETH1, and in situ hybridization with specific DNA probes. Selective RNA silencing of BdETH or BdETH-R-A, but not of BdETH-R-B, caused developmental failure at ecdysis. The dsRNA-treated larvae displayed tracheal defects and could not shed the old cuticle followed by death. Our results demonstrated that BdETH, via activation of BdETH-R-A but not ETH-R-B, plays an essential role in regulating the process of larva-larva ecdysis in B. dorsalis.

Phenotypes, antioxidant responses, and gene expression changes accompanying a sugar-only diet in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

BACKGROUND

Diet composition (yeast:carbohydrate ratio) is an important determinant of growth, development, and reproduction. Recent studies have shown that decreased yeast intake elicits numerous transcriptomic changes and enhances somatic maintenance and lifespan, which in turn reduces reproduction in various insects. However, our understanding of the responses leading to a decrease in yeast ratio to 0% is limited.

RESULTS

In the present study, we investigated the effects of a sugar-only diet (SD) on the gene expression patterns of the oriental fruit fly, Bactrocera dorsalis (Hendel), one of the most economically important pests in the family Tephritidae. RNA sequencing analyses showed that flies reared on an SD induced significant changes in the expression levels of genes associated with specific metabolic as well as cell growth and death pathways. Moreover, the observed upregulated genes in energy production and downregulated genes associated with reproduction suggested that SD affects somatic maintenance and reproduction in B. dorsalis. As expected, we observed that SD altered B. dorsalis phenotypes by significantly increasing stress (starvation and desiccation) resistance, decreasing reproduction, but did not extend lifespan compared to those that received a normal diet (ND) regime. In addition, administration of an SD resulted in a reduction in antioxidant enzyme activities and an increase in MDA concentrations, thereby suggesting that antioxidants cannot keep up with the increase in oxidative damage induced by SD regime.

CONCLUSIONS

The application of an SD diet induces changes in phenotypes, antioxidant responses, and gene expressions in B. dorsalis. Previous studies have associated extended lifespan with reduced fecundity. The current study did not observe a prolongation of lifespan in B. dorsalis, which instead incurred oxidative damage. The findings of the present study improve our understanding of the molecular, biochemical, and phenotypic response of B. dorsalis to an SD diet.

Influence of various stressors on the expression of core genes of the small interfering RNA pathway in the oriental fruit fly, Bactrocera dorsalis

RNA interference (RNAi)-based technology has emerged as a potential tool for controlling insect pests, however, previous studies found that the efficiency of RNAi in Bactrocera dorsalis was variable. In nature, insects often meet various challenges, such as pathogen infections, extreme temperatures, lack of nutrition and heavy metals. To better understand the association of the stressors with efficiency of RNAi, in the current study we tested the expression of three core genes, dicer2 (Bddcr2), r2d2 (Bdr2d2) and argonaute2 (Bdago2), of the small interfering RNA (siRNA) pathway of B. dorsalis upon various stressors. Our results showed that all three genes were upregulated by the infection of invertebrate iridescent virus 6, which suggested a function of the siRNA pathway against viral infection. The loading of FeCl₃ could also increase the expression of Bddcr2. The treatments of Escherichia coli, extremely high (40°C) and low (0°C) temperatures, as well as starvation, could negatively influence the expression of Bddcr2 and/or Bdago2. In total, our results showed that various stressors could influence the expression of core components of B. dorsalis siRNA pathway. This highlights further speculation on the RNAi efficiency upon these stressors. Considering the complexity and variation of RNAi efficiency in different conditions, these results provide initial aspects in possible environmental stressors to influence the activity of the siRNA pathway, but the real impact of RNAi efficiency posed by these stressors requires further studies.

The short neuropeptide F modulates olfactory sensitivity of Bactrocera dorsalis upon starvation

The insect short neuropeptide F (sNPF) family has been shown to modulate diverse physiological processes, such as feeding, appetitive olfactory behavior, locomotion, sleep homeostasis and hormone release. In this study, we identified the sNPF (BdsNPF) and its receptor (BdsNPFR) in an important agricultural pest, the oriental fruit fly Bactrocera dorsalis (Hendel). Afterwards, the receptor cDNA was functionally expressed in Chinese hamster ovary cell lines. Activation of BdsNPFR by sNPF peptides caused an increase in intracellular calcium ions, with a 50% effective concentration values at the nanomolar level. As indicated by qPCR, the BdsNPF and BdsNPFR transcripts were mainly detected in the central nervous system and antennae, and they showed significantly starvation-induced expression patterns. Furthermore, we found that the starved flies had an increased electroantennogram response compared to the normally fed flies. However, this enhanced olfactory sensitivity was reversed when we decreased the expression of BdsNPF by double-stranded RNA injection in adults. We concluded that sNPF plays an important role in modulating the olfactory sensitivity of B. dorsalis upon starvation. Our results will facilitate the understanding of the regulation of early olfactory processing in B. dorsalis.

Functional analysis of five trypsin-like protease genes in the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae)

Insect midgut proteases catalyze the release of free amino acids from dietary proteins and are essential for insect normal development. To date, digestive proteases as potential candidates have made great progress in pest control. To clarify the function of trypsin-like protease genes in the digestive system of Bactrocera dorsalis, a serious pest of a wide range of tropical and subtropical fruit and vegetable crops, five trypsin genes (BdTry1, BdTry2, BdTry3, BdTry4 and BdTry5) were identified from transcriptome dataset, and the effects of feeding condition on their expression levels were examined subsequently. RNA interference (RNAi) was applied to further explore their function on the growth of B. dorsalis. The results showed that all the BdTrys in starving midgut expressed at a minimal level but up-regulated upon feeding (except BdTry3). Besides, RNAi by feeding dsRNAs to larvae proved to be an effective method to cause gene silencing and the mixed dsRNAs of the five BdTrys slowed larvae growth of B. dorsalis. The current data suggest that trypsin genes are actively involved in digestion process of B. dorsalis larvae and thereafter play crucial roles in their development.

Molecular characterizations of natalisin and its roles in modulating mating in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Initially, natalisin (NTL) was identified from three holometabolous insect species, Drosophila melanogaster, Tribolium castaneum and Bombyx mori, and was documented to regulate reproductive behaviours in D. melanogaster and T. castaneum. In this study, we report the sequences of the NTL precursor and its receptor (NTLR) from an important agricultural pest, Bactrocera dorsalis (Hendel). NTLR is a typical G-protein coupled receptor and phylogenetic analysis showed that B. dorsalis NTLR was closely related to insect natalisin receptors from other species. A functional assay of NTLR transiently expressed in Chinese hamster ovary cells showed that it was activated by putative natalisin mature peptides in a concentration-dependent manner, with 50% effective concentrations (EC₅₀ ) at nanomolar or micromolar levels. As indicated by quantitative real-time PCR, both NTL and NTLR had the highest expression in the central nervous system of B. dorsalis compared with the other tested tissues. Three pairs of adult brain neurones of B. dorsalis were identified with immunohistochemical antibody staining against D. melanogaster NTL4, and in situ hybridization with specific DNA probes. Moreover, RNA interference mediated by double-stranded RNA injection in adults provided evidence for the important roles of NTL in regulating both male and female mating frequencies in this fly.

Role of a tachykinin-related peptide and its receptor in modulating the olfactory sensitivity in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Insect tachykinin-related peptide (TRP), an ortholog of tachykinin in vertebrates, has been linked with regulation of diverse physiological processes, such as olfactory perception, locomotion, aggression, lipid metabolism and myotropic activity. In this study, we investigated the function of TRP (BdTRP) and its receptor (BdTRPR) in an important agricultural pest, the oriental fruit fly Bactrocera dorsalis. BdTRPR is a typical G-protein coupled-receptor (GPCR), and it could be activated by the putative BdTRP mature peptides with the effective concentrations (EC₅₀) at the nanomolar range when expressed in Chinese hamster ovary cells. Consistent with its role as a neuromodulator, expression of BdTRP was detected in the central nervous system (CNS) of B. dorsalis, specifically in the local interneurons with cell bodies lateral to the antennal lobe. BdTRPR was found in the CNS, midgut and hindgut, but interestingly also in the antennae. To investigate the role of BdTRP and BdTRPR in olfaction behavior, adult flies were subjected to RNA interference, which led to a reduction in the antennal electrophysiological response and sensitivity to ethyl acetate in the Y-tube assay. Taken together, we demonstrate the impact of TRP/TRPR signaling on the modulation of the olfactory sensitivity in B. dorsalis. The result improve our understanding of olfactory processing in this agriculturally important pest insect.

Transcriptome analysis to identify genes for peptides and proteins involved in immunity and reproduction from male accessory glands and ejaculatory duct of Bactrocera dorsalis

In the male reproductive system of insects, the male accessory glands and ejaculatory duct (MAG/ED) are important organs and their primary function is to enhance the fertility of spermatozoa. Proteins secreted by the MAG/ED are also known to induce post-mating changes and immunity responses in the female insect. To understand the gene expression profile in the MAG/ED of the oriental fruit fly Bactrocera dorsalis (Hendel), that is an important pest in fruits, we performed an Illumina-based deep sequencing of mRNA. This yielded 54,577,630 clean reads corresponding to 4.91Gb total nucleotides that were assembled and clustered to 30,669 unigenes (average 645bp). Among them, 20,419 unigenes were functionally annotated to known proteins/peptides in Gene Orthology, Clusters of Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes pathway databases. Typically, many genes were involved in immunity and these included microbial recognition proteins and antimicrobial peptides. Subsequently, the inducible expression of these immunity-related genes was confirmed by qRT-PCR analysis when insects were challenged with immunity-inducible factors, suggesting their function in guaranteeing fertilization success. Besides, we identified some important reproductive genes such as juvenile hormone- and ecdysteroid-related genes in this de novo assembly. In conclusion, this transcriptomic sequencing of B. dorsalis MAG/ED provides insights to facilitate further functional research of reproduction, immunity and molecular evolution of reproductive proteins in this important agricultural pest.

Functional characterization of NADPH-cytochrome P450 reductase from Bactrocera dorsalis: Possible involvement in susceptibility to malathion

NADPH cytochrome P450 reductase (CPR) is essential for cytochrome P450 catalysis, which is important in the detoxification and activation of xenobiotics. In this study, two transcripts of Bactrocera dorsalis CPR (BdCPR) were cloned, and the deduced amino-acid sequence had an N-terminus membrane anchor for BdCPR-X1 and three conserved binding domains (FMN, FAD, and NADP), as well as an FAD binding motif and catalytic residues for both BdCPR-X1 and BdCPR-X2. BdCPR-X1 was detected to have the high expression levels in adults and in Malpighian tubules, fat bodies, and midguts of adults, but BdCPR-X2 expressed lowly in B. dorsalis. The levels of BdCPRs were similar in malathion-resistant strain compared to susceptible strain. However, injecting adults with double-stranded RNA against BdCPR significantly reduced the transcript levels of the mRNA, and knockdown of BdCPR increased adult susceptibility to malathion. Expressing complete BdCPR-X1 cDNA in Sf9 cells resulted in high activity determined by cytochrome c reduction and these cells had higher viability after exposure to malathion than control. The results suggest that BdCPR could affect the susceptibility of B. dorsalis to malathion and eukaryotic expression of BdCPR would lay a solid foundation for further investigation of P450 in B. dorsalis.

The Microbiome of Field-Caught and Laboratory-Adapted Australian Tephritid Fruit Fly Species with Different Host Plant Use and Specialisation

Tephritid fruit fly species display a diversity of host plant specialisation on a scale from monophagy to polyphagy. Furthermore, while some species prefer ripening fruit, a few are restricted to damaged or rotting fruit. Such a diversity of host plant use may be reflected in the microbial symbiont diversity of tephritids and their grade of dependency on their microbiomes. Here, we investigated the microbiome of six tephritid species from three genera, including species that are polyphagous pests (Bactrocera tryoni, Bactrocera neohumeralis, Bactrocera jarvisi, Ceratitis capitata) and a monophagous specialist (Bactrocera cacuminata). These were compared with the microbiome of a non-pestiferous but polyphagous tephritid species that is restricted to damaged or rotting fruit (Dirioxa pornia). The bacterial community associated with whole fruit flies was analysed by 16S ribosomal DNA (rDNA) amplicon pyrosequencing to detect potential drivers of taxonomic composition. Overall, the dominant bacterial families were Enterobacteriaceae and Acetobacteraceae (both Proteobacteria), and Streptococcaceae and Enterococcaceae (both Firmicutes). Comparisons across species and genera found different microbial composition in the three tephritid genera, but limited consistent differentiation between Bactrocera species. Within Bactrocera species, differentiation of microbial composition seemed to be influenced by the environment, possibly including their diets; beyond this, tephritid species identity or ecology also had an effect. The microbiome of D. pornia was most distinct from the other five species, which may be due to its ecologically different niche of rotting or damaged fruit, as opposed to ripening fruit favoured by the other species. Our study is the first amplicon pyrosequencing study to compare the microbiomes of tephritid species and thus delivers important information about the turnover of microbial diversity within and between fruit fly species and their potential application in pest management strategies.

Discovery of Chemosensory Genes in the Oriental Fruit Fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis, is a devastating fruit fly pest in tropical and sub-tropical countries. Like other insects, this fly uses its chemosensory system to efficiently interact with its environment. However, our understanding of the molecular components comprising B. dorsalis chemosensory system is limited. Using next generation sequencing technologies, we sequenced the transcriptome of four B. dorsalis developmental stages: egg, larva, pupa and adult chemosensory tissues. A total of 31 candidate odorant binding proteins (OBPs), 4 candidate chemosensory proteins (CSPs), 23 candidate odorant receptors (ORs), 11 candidate ionotropic receptors (IRs), 6 candidate gustatory receptors (GRs) and 3 candidate sensory neuron membrane proteins (SNMPs) were identified. The tissue distributions of the OBP and CSP transcripts were determined by RT-PCR and a subset of nine genes were further characterized. The predicted proteins from these genes shared high sequence similarity to Drosophila melanogaster pheromone binding protein related proteins (PBPRPs). Interestingly, one OBP (BdorOBP19c) was exclusively expressed in the sex pheromone glands of mature females. RT-PCR was also used to compare the expression of the candidate genes in the antennae of male and female B. dorsalis adults. These antennae-enriched OBPs, CSPs, ORs, IRs and SNMPs could play a role in the detection of pheromones and general odorants and thus could be useful target genes for the integrated pest management of B. dorsalis and other agricultural pests.

Insulin signaling pathway in the oriental fruit fly: The role of insulin receptor substrate in ovarian development

Insulin signaling pathways have integral roles in regulating organ growth and body size of insects. Here, we identified and characterized six insulin signaling pathway components-InR, IRS, PI3K92E, PI3K21B, Akt, and PDK-from Bactrocera dorsalis. Quantitative real-time polymerase chain reaction was used to establish gene expression profiles for the insulin signaling pathway components for different developmental stages and tissues, and in response to 20-hydroxyecdysone (20E) and starvation. IRS, PI3K92E, and PI3K21B were highly expressed in the head, while InR, Akt, and PDK were most abundant in Malpighian tubules. Both IRS and PI3K92E were highly expressed during the larval-pupal and pupal-adult transition, while the remaining four genes were highly expressed only during the pupal-adult transition. Following initial exposure to 20E, the expression levels of most genes were significantly decreased. However, the expression levels of IRS, PI3K92E, and PI3K21B were significantly increased at 8 and 12h post-treatment compared with the control. Moreover, we found that most insulin signaling pathway genes in B. dorsalis were up-regulated in response to starvation, but decreased when re-fed. On the contrary, transcript levels of PI3K21B decreased significantly during starvation. Furthermore, injection of IRS dsRNA into adult females significantly reduced IRS transcript levels. Suppression of IRS expression inhibited ovarian development, and the average ovary size was reduced by 33% compared with the control. This study provides new insight into the roles of insulin signaling pathway components in B. dorsalis, and demonstrates an important role for IRS in ovarian development.

Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis

The testis is a highly specialized tissue that plays a vital role in ensuring fertility by producing spermatozoa, which are transferred to the female during mating. Spermatogenesis is a complex process, resulting in the production of mature sperm, and involves significant structural and biochemical changes in the seminiferous epithelium of the adult testis. The identification of genes involved in spermatogenesis of Bactrocera dorsalis (Hendel) is critical for a better understanding of its reproductive development. In this study, we constructed a cDNA library of testes from male B. dorsalis adults at different ages, and performed de novo transcriptome sequencing to produce a comprehensive transcript data set, using Illumina sequencing technology. The analysis yielded 52 016 732 clean reads, including a total of 4.65 Gb of nucleotides. These reads were assembled into 47 677 contigs (average 443 bp) and then clustered into 30 516 unigenes (average 756 bp). Based on BLAST hits with known proteins in different databases, 20 921 unigenes were annotated with a cut-off E-value of 10(-5). The transcriptome sequences were further annotated using the Clusters of Orthologous Groups, Gene Orthology and the Kyoto Encyclopedia of Genes and Genomes databases. Functional genes involved in spermatogenesis were analysed, including cell cycle proteins, metalloproteins, actin, and ubiquitin and antihyperthermia proteins. Several testis-specific genes were also identified. The transcripts database will help us to understand the molecular mechanisms underlying spermatogenesis in B. dorsalis. Furthermore, 2913 simple sequence repeats and 151 431 single nucleotide polymorphisms were identified, which will be useful for investigating the genetic diversity of B. dorsalis in the future.

Characterizing the developmental transcriptome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae) through comparative genomic analysis with Drosophila melanogaster utilizing modENCODE datasets

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis, is an important pest of fruit and vegetable crops throughout Asia, and is considered a high risk pest for establishment in the mainland United States. It is a member of the family Tephritidae, which are the most agriculturally important family of flies, and can be considered an out-group to well-studied members of the family Drosophilidae. Despite their importance as pests and their relatedness to Drosophila, little information is present on B. dorsalis transcripts and proteins. The objective of this paper is to comprehensively characterize the transcripts present throughout the life history of B. dorsalis and functionally annotate and analyse these transcripts relative to the presence, expression, and function of orthologous sequences present in Drosophila melanogaster.

RESULTS

We present a detailed transcriptome assembly of B. dorsalis from egg through adult stages containing 20,666 transcripts across 10,799 unigene components. Utilizing data available through Flybase and the modENCODE project, we compared expression patterns of these transcripts to putative orthologs in D. melanogaster in terms of timing, abundance, and function. In addition, temporal expression patterns in B. dorsalis were characterized between stages, to establish the constitutive or stage-specific expression patterns of particular transcripts. A fully annotated transcriptome assembly is made available through NCBI, in addition to corresponding expression data.

CONCLUSIONS

Through characterizing the transcriptome of B. dorsalis through its life history and comparing the transcriptome of B. dorsalis to the model organism D. melanogaster, a database has been developed that can be used as the foundation to functional genomic research in Bactrocera flies and help identify orthologous genes between B. dorsalis and D. melanogaster. This data provides the foundation for future functional genomic research that will focus on improving our understanding of the physiology and biology of this species at the molecular level. This knowledge can also be applied towards developing improved methods for control, survey, and eradication of this important pest.

Deep sequencing of small RNA libraries reveals dynamic expression patterns of microRNAs in multiple developmental stages of Bactrocera dorsalis

In eukaryotes, microRNAs (miRNAs) are small, conserved, noncoding RNAs that have emerged as critical regulators of gene expression. The oriental fruit fly Bactrocera dorsalis is one of the most economically important fruit fly pests in East Asia and the Pacific. Although transcriptome analyses have greatly enriched our knowledge of its structural genes, little is known about post-transcriptional regulation by miRNAs in this dipteran species. In this study, small RNA libraries corresponding to four B. dorsalis developmental stages (eggs, larvae, pupae and adults) were constructed and sequenced. Approximately 30.7 million reads of 18-30 nucleotides were obtained, with 123 known miRNAs and 60 novel miRNAs identified amongst these libraries. More than half of the miRNAs were stage-specific during the four developmental stages. A set of miRNAs was found to be up- or down-regulated during development by comparison of their reads at different developmental stages. Moreover, a small part of miRNAs owned both miR-#-3p and miR-#-5p types, with enormously variable miR-#-3p/miR-#-5p ratios in the same library and amongst different developmental stages for each miRNA. Taking these findings together, the current study has uncovered a number of miRNAs and provided insights into their possible involvement in developmental regulation by expression profiling of miRNAs. Further analyses of the expression and function of these miRNAs could increase our understanding of regulatory networks in this insect and lead to novel approaches for its control.

De novo transcriptome of the Hemimetabolous German cockroach (Blattella germanica)

Background

The German cockroach, Blattella germanica, is an important insect pest that transmits various pathogens mechanically and causes severe allergic diseases. This insect has long served as a model system for studies of insect biology, physiology and ecology. However, the lack of genome or transcriptome information heavily hinder our further understanding about the German cockroach in every aspect at a molecular level and on a genome-wide scale. To explore the transcriptome and identify unique sequences of interest, we subjected the B. germanica transcriptome to massively parallel pyrosequencing and generated the first reference transcriptome for B. germanica.

Methodology/Principal Findings

A total of 1,365,609 raw reads with an average length of 529 bp were generated via pyrosequencing the mixed cDNA library from different life stages of German cockroach including maturing oothecae, nymphs, adult females and males. The raw reads were de novo assembled to 48,800 contigs and 3,961 singletons with high-quality unique sequences. These sequences were annotated and classified functionally in terms of BLAST, GO and KEGG, and the genes putatively coding detoxification enzyme systems, insecticide targets, key components in systematic RNA interference, immunity and chemoreception pathways were identified. A total of 3,601 SSRs (Simple Sequence Repeats) loci were also predicted.

Conclusions/Significance

The whole transcriptome pyrosequencing data from this study provides a usable genetic resource for future identification of potential functional genes involved in various biological processes.

Alternative splicing contributes to the coordinated regulation of ferritin subunit levels in Bactrocera dorsalis (Hendel)

A constant ratio of ferritin heavy chain homolog (HCH) and light chain homolog (LCH) subunits seems to be required to compose the ferritin heteropolymer protein in insects. However, the mechanism by which insect LCH genes regulate protein levels remains unclear. We report that alternative promoters and alternative splicing contribute to maintaining a constant ratio of the two subunits, BdFer1HCH and BdFer2LCH (ferritin 1 HCH and ferritin 2 LCH), in Bactrocera dorsalis, a notorious quarantine pest. The genes BdFer1HCH and BdFer2LCH were identified with a series of potential transcription factor binding sites and were shown to be clustered within the genome in a "head to head" fashion. Thus, we unearthed a potential post-transcriptional mechanism to regulate the levels of LCH subunits, and confirmed that the expressions of BdFer1HCH and BdFer2LCH were induced by 20-hydroecdysone, iron overload, and immune challenge.

De novo cloning and annotation of genes associated with immunity, detoxification and energy metabolism from the fat body of the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis, is a destructive pest in tropical and subtropical areas. In this study, we performed transcriptome-wide analysis of the fat body of B. dorsalis and obtained more than 59 million sequencing reads, which were assembled into 27,787 unigenes with an average length of 591 bp. Among them, 17,442 (62.8%) unigenes matched known proteins in the NCBI database. The assembled sequences were further annotated with gene ontology, cluster of orthologous group terms, and Kyoto encyclopedia of genes and genomes. In depth analysis was performed to identify genes putatively involved in immunity, detoxification, and energy metabolism. Many new genes were identified including serpins, peptidoglycan recognition proteins and defensins, which were potentially linked to immune defense. Many detoxification genes were identified, including cytochrome P450s, glutathione S-transferases and ATP-binding cassette (ABC) transporters. Many new transcripts possibly involved in energy metabolism, including fatty acid desaturases, lipases, alpha amylases, and trehalose-6-phosphate synthases, were identified. Moreover, we randomly selected some genes to examine their expression patterns in different tissues by quantitative real-time PCR, which indicated that some genes exhibited fat body-specific expression in B. dorsalis. The identification of a numerous transcripts in the fat body of B. dorsalis laid the foundation for future studies on the functions of these genes.

Characterization of Bactrocera dorsalis serine proteases and evidence for their indirect role in insecticide tolerance

The oriental fruit fly Bactrocera dorsalis (Hendel) causes devastating losses to agricultural crops world-wide and is considered to be an economically important pest. Little is known about the digestive enzymes such as serine proteases (SPs) in B. dorsalis, which are important both for energy supply and mitigation of fitness cost associated with insecticide tolerance. In this study, we identified five SP genes in the midgut of B. dorsalis, and the alignments of their deduced amino acid sequences revealed the presence of motifs conserved in the SP superfamily. Phylogenetic analyses with known SPs from other insect species suggested that three of them were trypsin-like proteases. Analyses of the expression profiles among the different developmental stages showed that all five genes were most abundant in larvae than in other stages. When larvae were continuously fed on diet containing 0.33 μg/g β-Cypermethrin, expression of all five genes were upregulated in the midgut but the larval development was delayed. Biochemical assays were consistent with the increased protease activity exhibited by SPs in the midgut after treatment with β-Cypermethrin. Taken together, these findings provide evidence for the hypothesis that enhanced SP activity may play an indirect role in relieving the toxicity stress of insecticide in B. dorsalis.