De novo Transcriptome Analysis of Chinese Citrus Fly, Bactrocera minax (Diptera: Tephritidae), by High-Throughput Illumina Sequencing

Trehalose-6-phosphate phosphatases are involved in trehalose synthesis and metamorphosis in Bactrocera minax

Trehalose is the principal sugar circulating in the hemolymph of insects, and trehalose synthesis is catalyzed by trehalose-6-phosphate synthase (TPS) and trehalose-6-phosphate phosphatase (TPP). Insect TPS is a fused enzyme containing both TPS domain and TPP domain. Thus, many insects do not possess TPP genes as TPSs have replaced the function of TPPs. However, TPPs are widely distributed across the dipteran insects, while the roles they play remain largely unknown. In this study, 3 TPP genes from notorious dipteran pest Bactrocera minax (BmiTPPB, BmiTPPC1, and BmiTPPC2) were identified and characterized. The different temporal-spatial expression patterns of 3 BmiTPPs implied that they exert different functions in B. minax. Recombinant BmiTPPs were heterologously expressed in yeast cells, and all purified proteins exhibited enzymatic activities, despite the remarkable disparity in performance between BmiTPPB and BmiTPPCs. RNA interference revealed that all BmiTPPs were successfully downregulated after double-stranded RNA injection, leading to decreased trehalose content and increased glucose content. Also, suppression of BmiTPPs significantly affected expression of downstream genes and increased the mortality and malformation rate. Collectively, these results indicated that all 3 BmiTPPs in B. minax are involved in trehalose synthesis and metamorphosis. Thus, these genes could be evaluated as insecticidal targets for managing B. minax, and even for other dipteran pests.

Identification and functional analysis of a chemosensory protein from Bactrocera minax (Diptera: Tephritidae)

BACKGROUND

Olfaction has an indispensable role in insect behavior, enabling location of suitable host plants and oviposition sites, finding mates and evasion of natural enemies. Chemosensory proteins (CSPs) function to screen external odorants and transport them to olfactory receptor neurons, thereby increasing the sensitivity of the olfactory system. At present, CSP genes have been identified in many insect species, but there are relatively few studies on the function of CSP, especially in Tephritidae.

RESULTS

In this study, we sequenced and analyzed 12 transcriptomes of Bactrocera minax and identified five CSP genes. The results of polymerase chain reactions with reverse transcription showed that BminCSP3 was highly expressed only in antennae. Results from competitive binding experiments showed that BminCSP3 has good binding ability to citral compared with 23 other volatile organic compounds. The docking model with citral showed hydrogen bond formation with residues (ARG97); however, no hydrogen bonds were formed in the docking of five other ligands (furfuryl alcohol, linalool, cis-3-hexenyl acetate, (R)-(+)-limonene and (+)-carvone). Electroantennogram (EAG) analyses revealed that citral was active in B. minax at the antennal level, and the EAG response value of female adults was significantly higher than that of male adults. Furthermore, the results of behavioral bioassays showed that females were significantly attracted to citral.

CONCLUSION

Our results suggest that BminCSP3 plays an important role in the recognition of citral by B. minax adults. © 2022 Society of Chemical Industry.

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.

The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock

Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.

De novo transcriptome analysis and identification of genes associated with immunity, detoxification and energy metabolism from the fat body of the tephritid gall fly, Procecidochares utilis

The fat body, a multifunctional organ analogous to the liver and fat tissue of vertebrates, plays an important role in insect life cycles. The fat body is involved in protein storage, energy metabolism, elimination of xenobiotics, and production of immunity regulator-like proteins. However, the molecular mechanism of the fat body's physiological functions in the tephritid stem gall-forming fly, Procecidochares utilis, are still unknown. In this study, we performed transcriptome analysis of the fat body of P. utilis using Illumina sequencing technology. In total, 3.71 G of clean reads were obtained and assembled into 30,559 unigenes, with an average length of 539 bp. Among those unigenes, 21,439 (70.16%) were annotated based on sequence similarity to proteins in NCBI's non-redundant protein sequence database (Nr). Sequences were also compared to NCBI's non-redundant nucleotide sequence database (Nt), a manually curated and reviewed protein sequence database (SwissProt), and KEGG and gene ontology annotations were applied to better understand the functions of these unigenes. A comparative analysis was performed to identify unigenes related to detoxification, immunity and energy metabolism. Many unigenes involved in detoxification were identified, including 50 unigenes of putative cytochrome P450s (P450s), 18 of glutathione S-transferases (GSTs), 35 of carboxylesterases (CarEs) and 26 of ATP-binding cassette (ABC) transporters. Many unigenes related to immunity were identified, including 17 putative serpin genes, five peptidoglycan recognition proteins (PGRPs) and four lysozyme genes. In addition, unigenes potentially involved in energy metabolism, including 18 lipase genes, five fatty acid synthase (FAS) genes and six elongases of very long chain fatty acid (ELOVL) genes, were identified. This transcriptome improves our genetic understanding of P. utilis and the identification of a numerous transcripts in the fat body of P. utilis offer a series of valuable molecular resources for future studies on the functions of these genes.

Identification of olfactory genes and functional analysis of BminCSP and BminOBP21 in Bactrocera minax

Insects possess highly developed olfactory systems which play pivotal roles in its ecological adaptations, host plant location, and oviposition behavior. Bactrocera minax is an oligophagous tephritid insect whose host selection, and oviposition behavior largely depend on the perception of chemical cues. However, there have been very few reports on molecular components related to the olfactory system of B. minax. Therefore, the transcriptome of B. minax were sequenced in this study, with 1 candidate chemosensory protein (CSP), 21 candidate odorant binding proteins (OBPs), 53 candidate odorant receptors (ORs), 29 candidate ionotropic receptors (IRs) and 4 candidate sensory neuron membrane proteins (SNMPs) being identified. After that, we sequenced the candidate olfactory genes and performed phylogenetic analysis. qRT-PCR was used to express and characterize 9 genes in olfactory and non-olfactory tissues. Compared with GFP-injected fly (control), dsOBP21-treated B. minax and dsCSP-treated B. minax had lower electrophysiological response to D-limonene (attractant), suggesting the potential involvement of BminOBP21 and BminCSP genes in olfactory perceptions of the fly. Our study establishes the molecular basis of olfaction, tributary for further functional analyses of chemosensory processes in B. minax.

De novo transcriptomic analysis of the alimentary tract of the tephritid gall fly, Procecidochares utilis

The tephritid gall fly, Procecidochares utilis, is an important obligate parasitic insect of the malignant weed Eupatorium adenophorum which biosynthesizes toxic secondary metabolites. Insect alimentary tracts secrete several enzymes that are used for detoxification, including cytochrome P450s, glutathione S-transferases, and carboxylesterases. To explore the adaptation of P. utilis to its toxic host plant, E. adenophorum at molecular level, we sequenced the transcriptome of the alimentary tract of P. utilis using Illumina sequencing. Sequencing and de novo assembly yielded 62,443 high-quality contigs with an average length of 604 bp that were further assembled into 45,985 unigenes with an average length of 674 bp and an N50 of 983 bp. Among the unigenes, 30,430 (66.17%) were annotated by alignment against the NCBI non-redundant protein (Nr) database, while 16,700 (36.32%), 16,267 (35.37%), and 11,530 (25.07%) were assigned functions using the Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) databases, respectively. Using the comprehensive transcriptome data set, we manually identified several important gene families likely to be involved in the detoxification of toxic compounds including 21 unigenes within the glutathione S-transferase (GST) family, 22 unigenes within the cytochrome P450 (P450) family, and 16 unigenes within the carboxylesterase (CarE) family. Quantitative PCR was used to verify eight, six, and two genes of GSTs, P450s, and CarEs, respectively, in different P. utilis tissues and at different developmental stages. The detoxification enzyme genes were mainly expressed in the foregut and midgut. Moreover, the unigenes were higher expressed in the larvae, pupae, and 3-day adults, while they were expressed at lower levels in eggs. These transcriptomic data provide a valuable molecular resource for better understanding the function of the P. utilis alimentary canal. These identified genes could be pinpoints to address the molecular mechanisms of P. utilis interacting with toxic plant host.

Transcriptomic and metabolomic profiles of Chinese citrus fly, Bactrocera minax (Diptera: Tephritidae), along with pupal development provide insight into diapause program

The Chinese citrus fly, Bactrocera minax (Enderlein), is a devastating citrus pest in Asia. This univoltine insect enters obligatory pupal diapause in each generation, while little is known about the course and the molecular mechanisms of diapause. In this study, the course of diapause was determined by measuring the respiratory rate throughout the pupal stage. In addition, the variation of transcriptomic and metabolomic profiles of pupae at five developmental stages (pre-, early-, middle-, late-, and post-diapause) were evaluated by next-generation sequencing technology and ¹H nuclear magnetic resonance spectroscopy (NMR), respectively. A total of 4,808 genes were significantly altered in ten pairwise comparisons, representing major shifts in metabolism and signal transduction as well as endocrine system and digestive system. Gene expression profiles were validated by qRT-PCR analysis. In addition, 48 metabolites were identified and quantified by ¹H NMR. Nine of which significantly contributed to the variation in the metabolomic profiles, especially proline and trehalose. Moreover, the samples collected within diapause maintenance (early-, middle-, and late-diapause) only exhibited marginal transcriptomic and metabolomic variation with each other. These findings greatly improve our understanding of B. minax diapause and lay the foundation for further pertinent studies.

RNA interference of a trehalose-6-phosphate synthase gene reveals its roles during larval-pupal metamorphosis in Bactrocera minax (Diptera: Tephritidae)

Trehalose is the major blood sugar in insects, which plays a crucial role as an instant source of energy and the starting substrate for chitin biosynthesis. In insects, trehalose is synthesized by catalysis of an important enzyme, trehalose-6-phosphate synthase (TPS). In the present study, a trehalose-6-phosphate synthase gene from Bactrocera minax (BmTPS) was cloned and characterized. BmTPS contained an open reading frame of 2445 nucleotides encoding a protein of 814 amino acids with a predicted molecular weight of 92.05kDa. BmTPS was detectable in all developmental stages of Bactrocera minax and expressed higher in the final- (third-) instar larvae. Tissue-specific expression patterns of BmTPS showed that it was mainly expressed in the fat body. The 20-hydroxyecdysone (20E) induced the expression of BmTPS and three genes in the chitin biosynthesis pathway. Moreover, injection of double-stranded RNA into third-instar larvae successfully silenced the transcription of BmTPS in B. minax, and thereby decreased the activity of TPS and trehalose content. Additionally, silencing of BmTPS inhibited the expression of three key genes in the chitin biosynthesis pathway and exhibited 52% death and abnormal phenotypes. The findings demonstrate that BmTPS is indispensable for larval-pupal metamorphosis. Besides, the establishment of RNAi experimental system in B. minax would lay a solid foundation for further investigation of molecular biology and physiology of this pest.