Evaluation of endogenous references for gene expression profiling in different tissues of the oriental fruit fly Bactrocera dorsalis (Diptera: Tephritidae)

AANAT1 regulates insect midgut detoxification through the ROS/CncC pathway

Insecticide resistance has been a problem in both the agricultural pests and vectors. Revealing the detoxification mechanisms may help to better manage insect pests. Here, we showed that arylalkylamine N-acetyltransferase 1 (AANAT1) regulates intestinal detoxification process through modulation of reactive oxygen species (ROS)-activated transcription factors cap"n"collar isoform-C (CncC): muscle aponeurosis fibromatosis (Maf) pathway in both the oriental fruit fly, Bactrocera dorsalis, and the arbovirus vector, Aedes aegypti. Knockout/knockdown of AANAT1 led to accumulation of biogenic amines, which induced a decreased in the gut ROS level. The reduced midgut ROS levels resulted in decreased expression of CncC and Maf, leading to lower expression level of detoxification genes. AANAT1 knockout/knockdown insects were more susceptible to insecticide treatments. Our study reveals that normal functionality of AANAT1 is important for the regulation of gut detoxification pathways, providing insights into the mechanism underlying the gut defense against xenobiotics in metazoans.

Characterization of an agmatine N-acetyltransferase from Bactrocera dorsalis that modulates ovary development

Agmatine N-acetyltransferase (AgmNAT), which catalyzes the formation of N-acetylagmatine from acetyl-CoA and agmatine, is a member of the GCN5-related N-acetyltransferase family. So far, knowledge of the physiological roles of AgmNAT in insects is limited. Here, we identified one gene encoding protein homologous to that of Drosophila AgmNAT using sequence information from an activity-verified Drosophila AgmNAT in a BLAST search of the Bactrocera dorsalis genome. We expressed and purified B. dorsalis AgmNAT in Escherichia coli and used the purified enzyme to define the substrate specificity for acyl-CoA and amine substrates. Our application of the screening strategy to BdorAgmNAT led to the identification of agmatine as the best amine substrate for this enzyme, with the highest kcat/Km value. We successfully obtained a BdorAgmNAT knockout strain based on a wild-type strain (WT) using the CRISPR/Cas9 technique. The ovary development of the BdorAgmNAT knockout mutants was delayed for 10 days compared with the WT specimens. Moreover, mutants had a much smaller mature ovary size and laid far fewer eggs than WT. Loss of function of BdorAgmNAT caused by RNAi with mature WT females did not affect their fecundity. These findings indicate that BdorAgmNAT is critical for oogenesis. Our data provide the first evidence for AgmNAT in regulating ovary development.

Identification of antennal alternative splicing by combining genome and full-length transcriptome analysis in Bactrocera dorsalis

Alternative splicing is an essential post-transcriptional regulatory mechanism that diversifies gene function by generating multiple protein isoforms from a single gene and act as a crucial role in insect environmental adaptation. Olfaction, a key sense for insect adaptation, relies heavily on the antennae, which are the primary olfactory organs expressing most of the olfactory genes. Despite the extensive annotation of olfactory genes within insect antennal tissues facilitated by high-throughput sequencing technology advancements, systematic analyses of alternative splicing are still relatively less. In this study, we focused on the oriental fruit fly (Bactrocera dorsalis), a significant pest of fruit crops. We performed a detailed analysis of alternative splicing in its antennae by utilizing the full-length transcriptome of its antennal tissue and the insect's genome. The results revealed 8600 non-redundant full-length transcripts identified in the oriental fruit fly antennal full-length transcriptome, spanning 4,145 gene loci. Over 40% of these loci exhibited multiple isoforms. Among these, 161 genes showed sex-biased isoform switching, involving seven different types of alternative splicing. Notably, events involving alternative transcription start sites (ATSS) and alternative transcription termination sites (ATTS) were the most common. Of all the genes undergoing ATSS and ATTS alternative splicing between male and female, 32 genes were alternatively spliced in protein coding regions, potentially affecting protein function. These genes were categorized based on the length of the sex-biased isoforms, with the highest difference in isoform fraction (dIF) associated with the ATSS type, including genes such as BdorABCA13, BdorCAT2, and BdorTSN3. Additionally, transcription factor binding sites for doublesex were identified upstream of both BdorABCA13 and BdorCAT2. Besides being expressed in the antennal tissues, BdorABCA13 and BdorCAT2 are also expressed in the mouthparts, legs, and genitalia of both female and male adults, suggesting their functional diversity. This study reveals alternative splicing events in the antennae of Bactrophora dorsalis from two aspects: odorant receptor genes and other types of genes expressed in the antennae. This study not only provides a research foundation for understanding the regulation of gene function by alternative splicing in the oriental fruit fly but also offers new insights for utilizing olfaction-based behavioral manipulation techniques to manage this pest.

RNAi-mediated pest control targeting the Troponin I (wupA) gene in sweet potato weevil, Cylas formicarius

The sweet potato weevil (Cylas formicarius) is a critical pest producing enormous global losses in sweet potato crops. Traditional pest management approaches for sweet potato weevil, primarily using chemical pesticides, causes pollution, food safety issues, and harming natural enemies. While RNA interference (RNAi) is a promising environmentally friendly approach to pest control, its efficacy in controlling the sweet potato weevil has not been extensively studied. In this study, we selected a potential target for controlling C. formicarius, the Troponin I gene (wupA), which is essential for musculature composition and crucial for fundamental life activities. We determined that wupA is abundantly expressed throughout all developmental stages of the sweet potato weevil. We evaluated the efficiency of double-stranded RNAs in silencing the wupA gene via microinjection and oral feeding of sweet potato weevil larvae at different ages. Our findings demonstrate that both approaches significantly reduced the expression of wupA and produced high mortality. Moreover, the 1st instar larvae administered dswupA exhibited significant growth inhibition. We assessed the toxicity of dswupA on the no-target insect silkworm and assessed its safety. Our study indicates that wupA knockdown can inhibit the growth and development of C. formicarius and offer a potential target gene for environmentally friendly control.

Distribution analysis of TRH in Bactrocera dorsalis using a CRISPR/Cas9-mediated reporter knock-in strain

Although the study of many genes and their protein products is limited by the availability of high-quality antibodies, this problem could be solved by fusing a tag/reporter to an endogenous gene using a gene-editing approach. The type II bacterial CRISPR/Cas system has been demonstrated to be an efficient gene-targeting technology for many insects, including the oriental fruit fly Bactrocera dorsalis. However, knocking in, an important editing method of the CRISPR/Cas9 system, has lagged in its application in insects. Here, we describe a highly efficient homology-directed genome editing system for B. dorsalis that incorporates coinjection of embryos with Cas9 protein, guide RNA and a short single-stranded oligodeoxynucleotide donor. This one-step procedure generates flies carrying V5 tag (42 bp) in the BdorTRH gene. In insects, as in other invertebrates and in vertebrates, the neuronal tryptophan hydroxylase (TRH) gene encodes the rate-limiting enzyme for serotonin biosynthesis in the central nervous system. Using V5 monoclonal antibody, the distribution of TRH in B. dorsalis at different developmental stages was uncovered. Our results will facilitate the generation of insects carrying precise DNA inserts in endogenous genes and will lay foundation for the investigation of the neural mechanisms underlying the serotonin-mediated behaviour of B. dorsalis.

Selection and Validation of Reference Genes for Gene Expression in Bactericera gobica Loginova under Different Insecticide Stresses

Insecticide resistance has long been a problem in crop pest control. Bactericera gobica is a major pest on the well-known medicinal plants Lycium barbarum L. Investigating insecticide resistance mechanisms of B. gobica will help to identify pesticide reduction strategies to control the pest. Gene expression normalization by RT-qPCR requires the selection and validation of appropriate reference genes (RGs). Here, 15 candidate RGs were selected from transcriptome data of B. gobica. Their expression stability was evaluated with five algorithms (Delta Ct, GeNorm, Normfinder, BestKeeper and RefFinder) for sample types differing in response to five insecticide stresses and in four other experimental conditions. Our results indicated that the RGs RPL10 + RPS15 for Imidacloprid and Abamectin; RPL10 + AK for Thiamethoxam; RPL32 + RPL10 for λ-cyhalothrin; RPL10 + RPL8 for Matrine; and EF2 + RPL32 under different insecticide stresses were the most suitable RGs for RT-qPCR normalization. EF1α + RPL8, EF1α + β-actin, β-actin + EF2 and β-actin + RPS15 were the optimal combination of RGs under odor stimulation, temperature, developmental stages and both sexes, respectively. Overall, EF2 and RPL8 were the two most stable RGs in all conditions, while α-TUB and RPL32 were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target cytochrome P450 CYP6a1 gene between adult and nymph stages and under imidacloprid stress. The results of CYP6a1 expression were consistent with transcriptome data. This study is the first research on the most stable RG selection in B. gobica nymphs exposed to different insecticides, which will contribute to further research on insecticide resistance mechanisms in B. gobica.

Mitochondrial coding genes mediate insecticide tolerance in the oriental fruit fly, Bactrocera dorsalis (Hendel)

The oriental fruit fly, Bactrocera dorsalis (Hendel), an invasive insect pest infesting fruits and vegetables, possesses a remarkable capacity for environmental adaptation. The investigation of behind mechanisms of the stress adaptability in B. dorsalis holds significantly practical relevance. Previous studies on the molecular mechanism underlying stress resistance in B. dorsalis have predominantly focused on nuclear-coding genes, with limited exploration on organelle-coding genes. In this study, we assessed alterations in the mitochondrial physiological parameters of B. dorsalis under exposure to malathion, avermectin, and beta-cypermethrin at LD50 dosages. The results showed that all three insecticides were capable of reducing mitochondrial complex IV activity and ATP content. Expression patterns of mitochondrial coding genes across different developmental stages, tissues and insecticide exposures were analyzed by RT-qPCR. The results revealed that these mitochondrial coding genes were expressed in various tissues and at different developmental stages. Particularly noteworthy, atp6, cox2, and cytb exhibited substantial up-regulation in response to malathion and avermectin treatment. Furthermore, RNAi-mediated knockdown of atp6 and cox2 resulted in the increased toxicity of malathion and avermectin against B. dorsalis, and cox2 silencing was also associated with the decreased complex IV activity. These findings suggest that atp6 and cox2 most likely play pivotal roles in mediating tolerance or resistance to malathion and avermectin in B. dorsalis. Our results provide novel insights into the role of mitochondrial coding genes in conferring tolerance to insecticides in B. dorsalis, with practical implications for controlling this pest in the field.

MFS Transporter Bdorwp Does Not Affect Antennal Electrophysiology but Regulates Reproductive Behaviors in Bactrocera dorsalis

Developing behavioral modifying chemicals through molecular targets is a promising way to improve semiochemical-based technology for pest management. Identifying molecular targets that affect insect behavior largely relies on functional genetic techniques such as deletions, insertions, and substitutions. Selectable markers have thus been developed to increase the efficiency of screening for successful editing events. However, the effect of selectable markers on relevant phenotypic traits needs to be considered. In this study, we cloned the wp gene ofBactrocera dorsalis. Knocking out Bdorwp causes white pupae phenotypes. Reproductive behaviors in both males and females were strongly regulated by Bdorwp. Remarkably, Bdorwp did not affect the antennal electrophysiology response to 63 chemical components with various structures. It is recommended to indirectly apply Bdorwp as a selectable marker in functional gene research on behavioral modifying chemicals. Moreover, Bdorwp could also be a potential molecular target for developing new insecticides for tephritid species control.

A ferritin protein is involved in the development and reproduction of the whitefly, Bemisia tabaci

Ferritins are conserved iron-binding proteins that exist in most living organisms and play an essential role in the maintenance of cellular iron homeostasis. Although ferritin has been studied in many species, little is known about its role in the whitefly, Bemisia tabaci. In this study, we identified an iron-binding protein from B. tabaci and named it BtabFer1. The full-length cDNA of BtabFer1 is 1,043 bp and encodes a protein consisting of 224 amino acids with a deduced molecular weight of 25.26 kDa, and phylogenetic analysis shows that BtabFer1 is conserved among Hemiptera insects. The expression levels of BtabFer1 in different developmental stages and tissues were analyzed by real-time PCR, and results showed that BtabFer1 was ubiquitously expressed at all developmental stages and in all examined tissues. The RNAi-mediated knockdown of BtabFer1 caused a significant reduction in survival rate, egg production, and egg hatching rate of whiteflies. Knockdown of BtabFer1 also inhibited the transcription of genes in the juvenile hormone signal transduction pathway. Taken together, these results suggest that BtabFer1 plays a critical role in the development and reproduction of whiteflies. This study can broaden our understanding of ferritin in insect fecundity and development, as well as provide baseline data for future studies.

Identification and functional analysis of isopentenyl pyrophosphate isomerase genes in the whiteflies Bemisia tabaci (Hemiptera: Aleyrodidae)

The juvenile hormone (JH) plays a vital role in the regulation of a number of physiological processes, including development, reproduction, and ovarian maturation. Isopentenyl pyrophosphate isomerase (IPPI) is a key enzyme in the biosynthetic pathway of JH. In this study, we identified an isopentenyl pyrophosphate isomerase protein from Bemisia tabaci and named it BtabIPPI. The open reading frame (ORF) of BtabIPPI is 768 bp and encodes a protein of 255 amino acids that contains a conserved domain of the Nudix family. The temporal and spatial expression profiles showed that BtabIPPI was highly expressed in the female adults.RNA interference (RNAi)-mediated silencing of BtabIPPI reduced JH titers and the relative expression of vitellogenin receptor (VgR) and JH signaling pathway genes, resulting in a dramatic reduction in fecundity and hatchability. These results indicate that the BtabIPPI gene plays an important role in the female fecundity of B. tabaci. This study will broaden our understanding of the function of IPPI in regulating insect reproduction and provide a theoretical basis for targeting IPPI for pest control in the future.

High-Sucrose Diet Exposure on Larvae Contributes to Adult Fecundity and Insecticide Tolerance in the Oriental Fruit Fly, Bactrocera dorsalis (Hendel)

Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) is one of the broad host ranges and economically-important insect pests in tropical and subtropical areas. A wide range of hosts means they have strong adaptation ability to changes in dietary macronutrients (e.g., sucrose and protein). However, the effects of dietary conditions on the phenotypes and genotypes of B. dorsalis are still unclear. In this study, we aimed to investigate the effects of larval dietary sucrose on the life history traits and stress tolerance of B. dorsalis, and its defense response at the molecular level. The results showed that low-sucrose (LS) induced decreased body size, shortened developmental duration, and enhanced sensitivity to beta-cypermethrin. Otherwise, high-sucrose (HS) diet increased developmental duration, adult fecundity, and tolerance to malathion. Based on transcriptome data, 258 and 904 differentially expressed genes (DEGs) were identified in the NS (control) versus LS groups, and NS versus HS groups, respectively. These yielded DEGs were relevant to multiple specific metabolisms, hormone synthesis and signaling, and immune-related pathways. Our study will provide biological and molecular perspective to understand phenotypic adjustments to diets and the strong host adaptability in oriental fruit flies.

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA-gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA-gene pairs were identified from the miRNA-mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation-reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.

A Hepatocyte Nuclear Factor BtabHNF4 Mediates Desiccation Tolerance and Fecundity in Whitefly (Bemisia tabaci)

Hepatocyte nuclear factor 4 (HNF4) is essential for glucose homeostasis and lipid metabolism in insects. However, little is known about the role of HNF4 in whiteflies. In the present study, we identified a hepatocyte nuclear factor protein from Bemsia tabaci (Diptera: Drosophilidae) and named it BtabHNF4. The full-length of BtabHNF4 was 3,006 bp, encoding a sequence of 434 amino acids that contains a conserved zinc-finger DNA-binding domain (DBD) and a well-conserved ligand-binding domain (LBD). The temporal and spatial expression showed that BtabHNF4 was highly expressed in the female adult stage and abdominal tissues of B. tabaci. A leaf-mediated RNA interference method was used to explore the function of BtabHNF4 in whiteflies. Our results showed that the knockdown of BtabHNF4 influences the desiccation tolerance, egg production, and egg hatching rate of whiteflies. Additionally, BtabHNF4 silencing significantly inhibited the expression level of vitellogenin. These results expand the function of HNF4 and pave the way for understanding the molecular mechanisms of HNF4 in regulating multiple physiological processes.

Protein feeding mediates sex pheromone biosynthesis in an insect

Protein feeding is critical for male reproductive success in many insect species. However, how protein affects the reproduction remains largely unknown. Using Bactrocera dorsalis as the study model, we investigated how protein feeding regulated sex pheromone synthesis. We show that protein ingestion is essential for sex pheromone synthesis in male. While protein feeding or deprivation did not affect Bacillus abundance, transcriptome analysis revealed that sarcosine dehydrogenase (Sardh) in protein-fed males regulates the biosynthesis of sex pheromones by increasing glycine and threonine (sex pheromone precursors) contents. RNAi-mediated loss-of-function of Sardh decreases glycine, threonine, and sex pheromone contents and results in decreased mating ability in males. The study links male feeding behavior with discrete patterns of gene expression that plays role in sex pheromone synthesis, which in turn translates to successful copulatory behavior of the males.

Chromosome-level genome assembly reveals potential epigenetic mechanisms of the thermal tolerance in the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis (Hendel), has very strong ecological adaptability and phenotypic plasticity. Here, the genome of B. dorsalis was assembled into 549.45 Mb sequences with a contig N50 length of 12.81 Mb. Among, 95.67 % assembled genome sequences were anchored on six chromosomes with an N50 length of 94.63 Mb. According to the basic characteristics of the sex chromosomes of Tephritidae, the X chromosome of B. dorsalis was identified. Significant gene expansions were detected in several important gene families related to adaptability. In particular, we annotated 50 histone modification enzymes (HMEs) in this genome. A comparative transcriptome analysis indicated that 12 HME genes were differentially expressed in two thermo-tolerant strains (heat and cold). Interestingly, four and seven of the 12 HME genes responded to heat shock or cold hardening, respectively. These evidences suggested that the histone modification as an epigenetic modification may be involved in the thermal tolerance of B. dorsalis, but with different regulation mechanisms in thermal acclimation and hardening. The high quality genome of B. dorsalis provides an invaluable resource for further functional genomic study. Moreover, comparative genomic analysis will shed insights on revealing the mechanisms of adaptive evolution in this fly.

BdcSP10 is a prophenoloxidase-activating protease in Bactrocera dorsalis

Insects rely on a powerful and efficient innate immune system against microbial invaders. One of the most important immune processes is the melanization reaction, in which eumelanin is synthesized and deposited on the physically injured site or the surface of invading pathogens. The melanization reaction is mediated by prophenoloxidase (PPO), which is synthesized as an inactive zymogen and requires proteolytic activation through a clip serine protease cascade. This cascade has been characterized in several Lepidoptera insect species, but it is less understood in most Diptera insects. Here, with the means of reverse genetics and biochemistry, we characterized the function of a clip serine protease BdcSP10 from the oriental fruit fly Bactrocera dorsalis (Hendel), a significant agriculture pest to a broad variety of fruit and vegetable crops. BdcSP10 knockdown inhibited the melanization reaction and rendered adult flies more vulnerable to pathogenic infections. In addition, purified and activated BdcSP10 proteases promoted the melanization reaction in larval hemolymph and directly cleaved and activated purified PPO1 and PPO2 in vitro. Taken together, we identified BdcSP10 as a PPO-activating protease and validated its important role in the defense against microbial infection in B. dorsalis. This work broadens the understanding of the activation mechanism of the melanization reaction in Diptera insects.

Transcriptome-Based Selection and Validation of Reference Genes for Gene Expression in Goji Fruit Fly (Neoceratitis asiatica Becker) under Developmental Stages and Five Abiotic Stresses

Goji fruit fly, Neoceratitis asiatica, is a major pest on the well-known medicinal plant Lycium barbarum. Dissecting molecular mechanisms of infestation and host selection of N. asiatica will contribute to the determination of best management practices for pest fly control. Gene expression normalization by Real-time quantitative PCR (qPCR) requires the selection and validation of appropriate reference genes (RGs). Hence, 15 candidate RGs were selected from transcriptome data of N. asiatica. Their expression stability was evaluated with five algorithms (∆Ct, Normfinder, GeNorm, BestKeeper, and RefFinder) for sample types differing in the developmental stage, sex, tissue type, and in response to five different abiotic stresses. Our results indicated that the RGs β-Actin + GST for sex, RPL32 + EF1α for tissue type, RPS13+ EF1α for developmental stages along with odor stimulation, color induction, and starvation-refeeding stresses, EF1α + GAPDH under insecticide stress, RPS13 + RPS18 under temperature stress, respectively, were selected as the most suitable RGs for qPCR normalization. Overall, RPS18 and EF1α were the two most stable RGs in all conditions, while RPS15 and EF1β were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target odorant-binding protein OBP56a gene in male and female antennae, different tissues, and under odor stimulation. The results of OBP56a expression were consistent with transcriptome data. Our study is the first research on the most stable RGs selection in N. asiatica, which will facilitate further studies on the mechanisms of host selection and insecticide resistance in N. asiatica.

Knockdown of a Nicotinic Acetylcholine Receptor Subunit Gene Bdorβ1 Decreases Susceptibility to Oxa-Bridged trans- instead of cis-Nitromethylene Neonicotinoid Insecticides in Bactrocera dorsalis

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that mediate the fast action of acetylcholine in synaptic cholinergic transmissions. Insect nAChRs are the target of several classes of insecticides. Here, the full-length cDNA encoding a nAChR beta1 subunit (Bdorβ1) was identified and characterized from a destructive pest, Bactrocera dorsalis. The amino acid sequence of Bdorβ1 shows high identities to other insect nAChRs β1 subunits. Double injection of dsBdorβ1 reduced the expression of Bdorβ1 and in turn significantly decreased susceptibility to oxa-bridged trans- instead of cis-nitromethylene neonicotinoids. Our results support the involvement of Bdorβ1 in the susceptibility of B. dorsalis to oxa-bridged trans- instead of cis-nitromethylene neonicotinoids and imply that these two classes of neonicotinoids might be acting at different nAChR subtypes.

Gut bacteria induce oviposition preference through ovipositor recognition in fruit fly

Gut bacteria play important roles in insect life cycle, and various routes can be used by insects to effectively transmit their gut bacteria. However, it is unclear if the gut bacteria can spread by actively attracting their insect hosts, and the recognition mechanisms of host insects are poorly understood. Here, we explore chemical interactions between Bactrocera dorsalis and its gut bacterium Citrobacter sp. (CF-BD). We found that CF-BD could affect the development of host ovaries and could be vertically transmitted via host oviposition. CF-BD could attract B. dorsalis to lay eggs by producing 3-hexenyl acetate (3-HA) in fruits that were hosts of B. dorsalis. Furthermore, we found that B. dorsalis could directly recognize CF-BD in fruits with their ovipositors in which olfactory genes were expressed to bind 3-HA. This work reports an important mechanism concerning the active spread of gut bacteria in their host insects.

The modulatory effects of biogenic amines on male mating performance in Bactrocera dorsalis

In insects, the emergence of mating behavior requires the interplay among sex-determination hierarchy mechanisms that regulate sex-specific differentiation, perception and integration of different sensory cues, and precisely patterned behavioral outputs. Biogenic amines, including octopamine (OA), dopamine (DA), tyramine (TA), serotonin and histamine, have been identified and proposed as putative neurotransmitters, neurohormones and/or neuromodulators in the central nervous system of insects to influence multiple physiologies and behaviors. The current study provides the physiological roles and pharmacology of these biogenic amines in the mating performance of Bactrocera dorsalis. Silencing gene expressions coding for biosynthetic enzymes of DA and serotonin in male flies could decrease mating rates, while OA, TA and histamine had no such effects on mating. Furthermore, injection of DA or the DA receptor antagonist chlorpromazine could affect mating rate, as well as injection of serotonin. Pharmacological treatments with other biogenic amines or their receptor antagonists in male flies have no roles in regulating mating performance. We conclude that DA and its receptors are involved in regulating male mating behaviors in B. dorsalis, while changes in serotonin levels in male flies could also affect mating rates. In the current study, the modulatory effects of these biogenic amines on mating performance were investigated, and these results will be helpful in providing a new strategy for controlling B. dorsalis.

Serotonin modulates insect gut bacterial community homeostasis

Background

Metazoan guts are in permanent contact with microbial communities. However, the host mechanisms that have developed to manage the dynamic changes of these microorganisms and maintain homeostasis remain largely unknown.

Results

Serotonin (5-hydroxytryptamine [5-HT]) was found to modulate gut microbiome homeostasis via regulation of a dual oxidase (Duox) gene expression in both Bactrocera dorsalis and Aedes aegypti. The knockdown of the peripheral 5-HT biosynthetic gene phenylalanine hydroxylase (TPH) increased the expression of Duox and the activity of reactive oxygen species, leading to a decrease in the gut microbiome load. Moreover, the TPH knockdown reduced the relative abundance of the bacterial genera Serratia and Providencia, including the opportunistic pathogens, S. marcescens and P. alcalifaciens in B. dorsalis. Treatment with 5-hydroxytryptophan, a precursor of 5-HT synthesis, fully rescued the TPH knockdown-induced phenotype.

Conclusions

The findings reveal the important contribution of 5-HT in regulating gut homeostasis, providing new insights into gut–microbe interactions in metazoans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01319-x.

Selection of Reference Genes for RT-qPCR Analysis of Wing Dimorphism in English Grain Aphid, Sitobion avenae (Hemiptera: Aphididae)

The English grain aphid, Sitobion avenae (Fabricius), exhibits classic and dramatic phenotypic plasticity in wing development. Both genetic and environmental inputs contribute to the wing polyphenism in aphids, an extreme form of phenotypic plasticity in which a single genotype produces discrete winged and wingless morphs. Validated reference genes are needed to accurately normalize temporal and spatial variation in gene expression estimates by RT-qPCR. In this research, the stability of 11 candidate reference genes selected from S. avenae transcriptomes was evaluated under an array of abiotic and biotic conditions relevant to wing development. RefFinder, a comprehensive software integrating rankings from delta Ct, BestKeeper, NormFinder, and geNorm, offered a series of reference genes for every experimental condition. Overall, helicase (HEL) and ubiquitin ribosomal protein S27A fusion protein (RpS27) are suited for most of the conditions examined in this study, although exceptions do exist. Specifically, NADH dehydrogenase (Ap-NADH) and 28S ribosomal RNA (28S) are recommended for insecticide and antibiotic treatments, while ribosomal RNA L14 (RPL14) and 18S ribosomal RNA (18S) are selected for density treatment, respectively. This study provides a suite of reference genes to investigate the wing polyphenism in S. avenae, and is important for application of RT-qPCR in future experiments of novel tactics to control aphids.

Selection and Validation of Reference Genes for RT-qPCR Normalization in Bradysia odoriphaga (Diptera: Sciaridae) Under Insecticides Stress

Bradysia odoriphaga (Diptera: Sciaridae) is the most serious root maggot pest which causes substantial damage to the Chinese chive. Organophosphate (OP) and neonicotinoid insecticides are widely used chemical pesticides and play important roles in controlling B. odoriphaga. However, a strong selection pressure following repeated pesticide applications has led to the development of resistant populations of this insect. To understand the insecticide resistance mechanism in B. odoriphaga, gene expression analysis might be required. Appropriate reference gene selection is a critical prerequisite for gene expression studies, as the expression stability of reference genes can be affected by experimental conditions, resulting in biased or erroneous results. The present study shows the expression profile of nine commonly used reference genes [elongation factor 1α (EF-1α), actin2 (ACT), elongation factor 2α (EF-2α), glucose-6-phosphate dehydrogenase (G6PDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein S3 (RPS3), ubiquitin-conjugating enzyme (UBC), and α-tubulin (TUB)] was systematically analyzed under insecticide stress. Moreover, we also evaluated their expression stability in other experimental conditions, including developmental stages, sexes, and tissues. Five programs (NormFinder, geNorm, BestKeeper, RefFinder, and ΔCt) were used to validate the suitability of candidate reference genes. The results revealed that the most appropriate sets of reference genes were RPL10 and ACT across phoxim; ACT and TUB across chlorpyrifos and chlorfluazuron; EF1α and TUB across imidacloprid; EF1α and EF2α across developmental stages; RPL10 and TUB across larvae; EF1α and ACT across tissues, and ACT and G6PDH across sex. These results will facilitate the standardization of RT-qPCR and contribute to further research on B. odoriphaga gene function under insecticides stress.

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.

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.

Evaluation of Reference Genes for Quantitative Reverse Transcription Polymerase Chain Reaction in Bactrocera dorsalis (Diptera: Tephritidae) Subjected to Various Phytosanitary Treatments

Simple Summary

In this study, seven internal reference genes (G6PDH, GAPDH, RPL-32, Rpl-13, Rps-3, α-tub, and 18S) of Bactroceradorsalis under different quarantine treatments (heat treatment, cold treatment, methyl bromide fumigation, and irradiation) were screened. Finally, the most stable internal reference gene was selected, which laid a foundation for the further study of its resistance mechanisms to some abiotic stresses.

Abstract

Bactrocera dorsalis is a major pest that causes serious damage to many fruits. Although phytosanitary treatment methods have been developed for Bactrocera control, there is a lack of information related to the gene expression pattern of B. dorsalis subjected to phytosanitary treatment conditions. Prior to quantitative reverse transcription polymerase chain reaction analysis of the most stable reference genes in B. dorsalis (Diptera: Tephritidae), B. dorsalis third-instar larvae were exposed to various phytosanitary treatments; seven candidate reference genes (18S, G6PDH, GAPDH, RPL-13, RPL-32, RPS-3, and α-Tub) were amplified and their expression stabilities were evaluated using geNorm, NormFinder, BestKeeper, and RefFinder algorithms. Different reference genes were found under different stress conditions. G6PDH was the most stable gene after heat treatment. After cold treatment, α-Tub exhibited the highest expression stability. G6PDH expression stability was the highest after fumigation with methyl bromide. RPL-32 showed the highest expression stability after irradiation treatment. Collectively, RefFinder analysis results revealed G6PDH and RPL-32 as the most suitable genes for analyzing phytosanitary treatment in B. dorsalis. This study provides an experimental basis for further gene expression analyses in B. dorsalis subjected to various phytosanitary treatments, which can aid in the development of novel phytosanitary treatments against insect pests.

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.

Expression stability of candidate RT-qPCR housekeeping genes in Spodoptera frugiperda (Lepidoptera: Noctuidae)

Reverse-transcription quantitative polymerase chain reaction (RT-qPCR) is commonly used to quantify gene expression. For normalization, the expression of each gene is compared with a reference "housekeeping" gene that is stably expressed under relevant stress. Unfortunately, there have been no reports on the stability of such reference genes under various treatments of the Spodoptera frugiperda. In this study, we used five tools (RefFinder, GeNorm, NormFinder, BestKeeper, and ΔCt methods) to evaluate the stability of 12 candidate reference genes (RPS18, β-tubulin, GAPDH, RPS7, RPS15, RPL7, RPL32, Actin-5C, EF1-α, EF1-γ, RPL27, and ACE) in different instars, tissues, and treatments (high and low temperature, UV-A, and emamectin benzoate). Several ribosomal proteins (RPS7, RPS15, RPL32, RPS18, and RPL7), GAPDH, Actin-5C, and β-tubulin, were relatively stable, suggesting that they are ideal housekeeping genes for various treatments. ACE was extremely unstable under various experimental treatments, rendering it unsuitable as an internal reference. This study identified the reference housekeeping genes stably expressed by S. frugiperda under different treatments, thus setting a foundation for further exploration of the physiological and biochemical mechanisms.

Molecular Characterization and Transcriptional Expression Analysis of ABC Transporter H Subfamily Genes in the Oriental Fruit Fly

The oriental fruit fly, Bactrocera dorsalis Hendel (Diptera: Tephretidae), is a serious pest of fruits and vegetables and has developed high levels of insecticide resistance. ATP-binding cassette transporter genes (ABC transporters) are involved in mediating the energy-driven transport of many substances across membranes and are closely associated with development and insecticide detoxification. In this study, three ABC transporters in the H subfamily were identified, and the possible roles of these genes in B. dorsalis are discussed. Bioinformatics analysis revealed that those genes are conserved, typical of half-transporters. The expression profiles of BdABCH genes (BdABCHs) in the developmental stages, tissues, and following insecticide exposure, extreme temperature, warm- and cold-acclimated strain, starvation, and desiccation stress were determined by quantitative real-time PCR. Expression of BdABCHs can be detected in various tissues and in different developmental stages. They were most highly expressed in the hindgut and in newly emerged adults. The mRNA levels of BdABCHs in males (including most tissues and body segments) were higher than in females. The expression of BdABCH1 was significantly upregulated 3.8-fold in the cold-acclimated strain, and was significantly upregulated by 1.9-, 3.8- and 4.1-fold in the 0°C, starvation, and desiccation treatments, respectively. Treatment with malathion and avermectin at LD20 and LD30 concentrations produced no obvious changes in the levels of BdABCHs. BdABCHs may be involved in the transport of related hormones during eclosion, as well as water and inorganic salts. BdABCH1 also demonstrated that it is related to the ability to cope with adverse environments.

The adipokinetic hormone signaling system regulates the sensitivity of Bactrocera dorsalis to malathion

The neuropeptide adipokinetic hormone (AKH) binds to the AKH receptor (AKHR) to regulate carbohydrate and lipid metabolism. It also participates in the insect anti-stress response. We used RT-qPCR to detect the expression levels of 39 neuropeptides in malathion-susceptible (MS) and malathion-resistant (MR) strains of Bactrocera dorsalis. AKH and AKHR were highly expressed in the MR strain. Using a malathion bioassay and RNA interference (RNAi), we demonstrated that AKHR is involved in the susceptibility of B. dorsalis to malathion. We found significantly reduced expression of two detoxification enzyme genes (glutathione-S-transferase, GST and α-esterase, CarE) after AKHR RNAi. Based on our previous data, GSTd10 and CarE6 participate the direct metabolism of malathion in this fly, which is also verified by a malathion metabolism assay by HPLC using the crude enzymes in the current study. These results suggest that AKHR plays an important role in affecting malathion susceptibility via detoxification enzyme genes.

Genome-wide identification of long non-coding RNAs (lncRNAs) associated with malathion resistance in Bactrocera dorsalis

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in the regulation of biological processes and have been identified in many species including insects. However, the association between lncRNAs and pesticide resistance in insect species such as Bactrocera dorsalis is unknown.

RESULTS

RNA-seq was performed on malathion resistant (MR1) and susceptible (MS) strains of B. dorsalis and a total of 6171 lncRNAs transcripts were identified. These included 3728 lincRNAs, 653 antisense lncRNAs, 1402 intronic lncRNAs, and 388 sense lncRNAs. A total of 40 and 52 upregulated lncRNAs were found in females and males of the MR1 strain compared to 54 and 49 in the same sexes of the MS strain, respectively. Twenty-seven of these lncRNAs showed the same trend of expression in both females and males in the MR1 strain, in which 15 lncRNAs were upregulated and 12 were downregulated. RT-qPCR results indicated that the differentially expressed lncRNAs were associated with malathion resistance. The lnc15010.10 and lnc3774.2 were highly expressed in the cuticle of the MR1 strain, indicating that these two lncRNAs may be related to malathion resistance. RNAi of lnc3774.2 and a bioassay showed that malathion resistance was possibly influenced by changes in the B. dorsalis cuticle.

CONCLUSION

LncRNAs of B. dorsalis potentially related to the malathion resistance were identified. Two lncRNAs appear to influence malathion resistance via modulating the structure, or components, of the cuticle. © 2021 Society of Chemical Industry.

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.

Molecular Characteristics of Fat Body Protein 1 in the Oriental Fruit Fly, Bactrocera dorsalis

Simple Summary

Bactrocera dorsalis fat body protein 1 (Bdfbp1) cDNA was cloned. The deduced amino acid sequence contains three motifs: hemocyanin N, high molecular weight glutenin (gultenin hmw), and hemocyanin C from N to C termini. The glutenin hmw allows Bdfbp1 to fold into a compact form for storage. Bdfbp1 was highly expressed in the late third instar larvae and day 0 pupae. This suggests that Bdfbp1 is stored during larval stages as a storage protein for construction of adult tissues during pupal stages, and may be associated with adult eclosion.

Abstract

Bactrocera dorsails fat body protein 1 (Bdfbp1) cDNA was cloned (GenBank accession no. MT514270), and the complete 3,749-bp cDNA encoded a 1,152-amino acid protein. The phylogenetic relationship of dipteran fbp1s was analyzed. The sequence XP_028900815 from the insect genome project for Zeugodacus cucurbitae (LOC105219342) was proposed that two fbp1 genes were present in the sequence. The developmental transcriptional expression profiles were determined. In the larval stages, Bdfbp1 mRNA had significantly higher expression in the late third instar larvae compared with first, second, and early third instar larvae. In the pupal stages, the highest expression of Bdfbp1 mRNA was found in the newly pupated pupae and then decreased with age. In the fat body of female adults, Bdfbp1 was highly expressed in newly emerged samples and decreased rapidly over the following three days. In the fat body of male adults, Bdfbp1 was highly expressed in newly eclosed samples. RNAi treatment decreased the expression level of Bdfbp1 without statistical difference. However, RNAi treatment significantly decreased the rate of eclosion. These results suggest that Bdfbp1 may function as a storage protein and be associated with adult eclosion.

Selection of Reference Genes for Normalization of Gene Expression in Thermobia domestica (Insecta: Zygentoma: Lepismatidae)

Zygentoma occupies a key evolutionary position for understanding the evolution of insect metamorphosis but has received little attention in terms of genetic analysis. To develop functional genomic studies in this insect, we evaluated five candidate internal reference genes for quantitative RT-PCR (qPCR) studies from Thermobia domestica, a representative species of Zygentoma, including Actin 5C (Actin5C), Elongation factor-1 alpha (EF1A), Ribosome protein S26 (RPS26), Ribosome protein L32 (RPL32), and Superoxide dismutase 2 (SOD2), at different developmental stages, in various body parts, and under dsRNA microinjection and starvation stresses, using four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) and a comparative algorithm (RefFinder). Specific suitable reference genes were recommended across specific experimental conditions, and the combination of RPS26 and RPL32 was appropriate for all tested samples. Employing our selected reference gene combination, we investigated the gene expression pattern of Myoglianin (Myo), a crucial gene-regulating insect metamorphosis, in ametabolous T. domestica, and demonstrated the efficiency of RNA interference (RNAi) in firebrat nymphs. This study provides a basis for reliable quantitative studies of genes and greatly benefits evolutionary and functional genomics studies in Zygentoma.

Selection and validation of experimental condition-specific reference genes for qRT-PCR in Metopolophium dirhodum (Walker) (Hemiptera: Aphididae)

Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) is one of the most common aphid pests of winter cereals. To facilitate accurate gene expression analyses with qRT-PCR assays, the expression stability of candidate reference genes under specific experimental conditions must be verified before they can be used to normalize target gene expression levels. In this study, 10 candidate reference genes in M. dirhodum were analyzed by qRT-PCR under various experimental conditions. Their expression stability was evaluated with delta Ct, BestKeeper, geNorm, and NormFinder methods, and the final stability ranking was determined with RefFinder. The results indicate that the most appropriate sets of internal controls were SDHB and RPL8 across geographic population; RPL8, Actin, and GAPDH across developmental stage; SDHB and NADH across body part; RPL8 and Actin across wing dimorphism and temperature; RPL4 and EF1A across starvation stress; AK and RPL4 across insecticide treatments; RPL8 and NADH across antibiotic treatments; RPL8, RPL4, Actin, and NADH across all samples. The results of this study provide useful insights for establishing a standardized qRT-PCR procedure for M. dirhodum and may be relevant for identifying appropriate reference genes for molecular analyses of related insects.

The Bactrocera dorsalis caspase-1 gene is expressed throughout development and required for female fertility

BACKGROUND

The oriental fruit fly Bactrocera dorsalis is one of the most destructive pests of fruits and vegetables. The sterile insect technique (SIT) is an effective and environmentally friendly approach to the control of tephritid fruit flies. The pro-apoptotic gene head involution defective (hid) has been used as an effective lethal effector in SIT. It initiates an interaction cascade including activation of caspase-like proteases. However, the biological role of caspase activity in tephritid fruit flies has yet to be explored.

RESULTS

In this study, the B. dorsalis caspase-1 gene (Bdcp-1) was cloned and characterized. Sequence comparison showed that Bdcp-1 protein shared highly homology with Drosophila effector caspases Drice and Dcp-1. It is predicted to contain a short pro-domain because two proteolytic cleavage sites (Asp¹⁶ and Asp²²³ ) are present. Expression patterns indicated that Bdcp-1 is highly transcribed in embryos and expression was upregulated during metamorphosis and upon ultraviolet irradiation. RNA interference showed that Bdcp-1 is essential for ovarian development and female fertility. For example, knockdown of Bdcp-1 caused transcriptional downregulation of expression of the yolk protein-1 gene (Bdyp-1) and delayed ovarian development. The percentage of spawning females and female fecundity were significantly reduced.

CONCLUSION

This study illustrates the function of the Bdcp-1 gene and provides an attractive method to develop a biological way to control the oriental fruit fly through the control of caspases. © 2020 Society of Chemical Industry.

Development of Sarcophaga dux (diptera: Sarcophagidae) at constant temperatures and differential gene expression for age estimation of the pupae

Sarcophaga dux (Diptera: Sarcophagidae) is a necrophagous flesh fly species with potential forensic value for estimating minimum postmortem interval (PMI_min). The basic developmental data and precise age estimates of the pupae are significant for PMI_min estimation in forensic investigations. In the present study, we investigated the development data of that species at seven constant temperatures varying from 16 °C to 34 °C, including body length changes of the larve, developmental duration and accumulated degree hours of the preadults. Several reference genes for relative quantification of the differentially expressed genes (DEGs) were firstly selected and evaluated in the pupae of different ages under different temperatures. The DEGs of the insects during the pupal period at different constant temperatures (34, 25 and 16 °C) were further analyzed for more precise age estimation. The results showed that the developmental durations of the preadults at 16, 19, 22, 25, 28, 31 and 34 °C were 1478.6 ± 18.3 h, 726.1 ± 15.8 h, 538.5 ± 0.9 h, 394.1 ± 9.5 h, 375.6 ± 10.8 h, 284.1 ± 7.3 h, and 252.5 ± 6.1 h, respectively. The developmental threshold temperature the flies was 12.27 ± 0.35 °C, and the thermal summation constant was 5341.71 ± 249.29° hours. The most reliable reference genes during the pupal period at different temperatures were found: GST1 and 18S rRNA for the 34 °C group, GST1 and RPL49 for 25 °C, and 18S rRNA and 28S rRNA for 16 °C. The four differential expression genes (Hsp60, A-alpha, ARP, and RPL8) have the potential to be used for more precise age estimation of pupal S. dux. This work provides important basic developmental data and a more precise age estimation method for pupal S. dux, and improves the value of this species for PMI_min estimation in forensic investigations.

Evaluation of suitable reference genes for normalization of quantitative real-time PCR analysis in rice plants under Xanthomonas oryzae pv. oryzae--infection and melatonin supplementation

Exogenous melatonin (MT) was found to be an interesting tool for enhancing the resistance of rice to Xanthomonasoryzaepv. oryzae (Xoo)-caused bacterial blight (BB). However, the accurate comparison of the expression levels across samples was a challenging task. In this work, the stability of 10 common used housekeeping genes under Xoo-infection and MT supplementation in rice was analyzed using quantitative real-time PCR (qRT-PCR), and algorithms geNorm, NormFinder and BestKeeper. Our results indicated that most reference genes remained stable in Xoo-infected rice plants, while a number of reference genes were affected by MT supplementation. Among all studied genes, the transcript levels of 18S(18S ribosomal RNA) and UBC (Ubiquitin-conjugating enzyme E2) remained unaltered by Xoo infection, while UBC and UBQ5(Ubiquitin 5) were the most stable genes when examining simultaneous Xoo-infection and MT supplementation, demonstrating that UBC is a suitable reference gene for qRT-PCR data normalization in rice under Xoo-infection and MT supplementation.

The RNA Virome and Its Dynamics in an Invasive Fruit Fly, Bactrocera dorsalis, Imply Interactions Between Host and Viruses

The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important invasive agricultural insect pest with a wide host range, and has spread around the world over the last century. This evolutionary trait may have arisen primarily from interactions between B. dorsalis and other invertebrates that share the same ecological niches. The invasive behavior of B. dorsalis also frequently exposes them to diverse species of viruses. Thereby, RNA viromes may be useful microbial markers to understand the ecological evolution of B. dorsalis as well as to investigate virus-host interactions. Here, we reported eight novel RNA viruses in B. dorsalis of a lab colony, including four positive-strand RNA viruses, two negative-strand RNA viruses, and two double-stranded RNA viruses using high-throughput sequencing technology. Analysis of the virus-derived small RNAs suggested that most of these viruses may be active and trigger the host antiviral RNAi responses. The viruses were also detected in various geographical populations of B. dorsalis, implying that there is a strong association between the viromes and host. In addition, these viruses infected specific fly tissues, predominately the central nervous system and gut. Furthermore, we explored the dynamics of the viruses when hosts were exposed to short- or long-term stressors, which showed that titers of some viruses were responsively altered in the stressed B. dorsalis. The discovery of these viruses may enrich our understanding of the species diversity of RNA viruses and also provide information on viruses in association with host adaptation in insects.

Overexpression of Glutathione S-Transferase Genes in Field λ-Cyhalothrin-Resistant Population of Cydia pomonella: Reference Gene Selection and Expression Analysis

Analysis of the glutathione S-transferase (GST) gene expression in an insecticide-resistant strain of Cydia pomonella using real-time quantitative polymerase chain reaction is a key step toward more mechanism studies that require suitable reference genes with stable expression. Here, nine commonly used reference genes were selected, and their expression stabilities were analyzed. Results showed that EF-1α was the most stable reference gene in all of the experimental sets. The combinations of EF-1α and 18S, EF-1α and RPL12, and EF-1α and GAPDH were sufficient for normalization of gene expression analysis accurately in developmental stages, tissues, and larvae exposed to sublethal dose of λ-cyhalothrin, respectively. Additionally, the suitability of particular reference genes was verified by analyzing the spatiotemporal and insecticide-induced expression profiles of CpGSTe3, CpGSTd3, and CpGSTd4, which were overexpressed in a λ-cyhalothrin-resistant population from northeast China. These genes were used to confer the practicability of reference genes chosen in this study.

Divergent molecular evolution in glutathione S-transferase conferring malathion resistance in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Insect glutathione S-transferases (GSTs) are important in insecticide detoxification and Insect-specific GSTs, Epsilon and Delta, have largely expanded in insects. In this study, we functionally expressed and characterized an epsilon class GST gene (BdGSTe8), predominant in the adult Malpighian tubules of Bactrocera dorsalis. This gene may be associated with malathion resistance based on transcriptional studies of resistant and susceptible strains. RNA interference-mediated knockdown of this gene significantly recovered malathion susceptibility in the adults of a malathion-resistant strain, and overexpression of BdGSTe8 enhanced resistance in transgenic Drosophila. Analysis of BdGSTe8 polymorphism showed that several point mutations may be associated with metabolic resistance to malathion. A cytotoxicity assay in Escherichia coli indicated that both of the recombinant BdGSTe8 proteins may play a functional role in protecting cells from toxicity. The allele of BdGSTe8-B conferred higher levels of malathion detoxification capability. Liquid chromatography and ultra-performance liquid chromatography-tandem mass spectrometry analysis showed that the BdGSTe8-A allele did not metabolize malathion directly. However, the BdGSTe8-B allele was involved in the direct metabolism of malathion, which was caused by a mutation in V128A. Further analysis of the sequence suggests that BdGSTe8 evolved rapidly. It maybe play the role of a backup gene and could become a new gene in the future in order to retain the ability of detoxification of malathion, which was driven by positive selection. These results suggest that divergent molecular evolution in BdGSTe8 has played a role in metabolic resistance to malathion in B. dorsalis.

Selection of reference genes for tissue/organ samples of adults of Eucryptorrhynchus scrobiculatus

Eucryptorrhynchus scrobiculatus is an important wood-boring pest of Ailanthus altissima in China, where it causes a large number of these trees to weaken or even die. To identify genes related to economic traits or specific cellular processes in E. scrobiculatus, gene expression in multiple tissue/organ samples is commonly surveyed, and reference genes are required in this process as a control for normalization. In the present study, 18 candidate reference genes from E. scrobiculatus were identified, and the expression levels of these reference genes were estimated through quantitative real-time PCR. Differences in expression levels were analyzed with four algorithms (geNorm, NormFinder, BestKeeper and delta Ct) and comprehensively with RefFinder. With the most stable levels of expression in different tissues, RPL13, RPS3 and RPL36 were determined to be suitable for use as candidate reference genes. Moreover, the expression profile of one target gene (glycoside hydrolase family 45, GH45) confirmed the reliability of the selected candidate reference genes. This study provides the first set of suitable candidate reference genes for gene expression studies in E. scrobiculatus, and the findings will facilitate subsequent transcriptomics and functional gene research on this pest.

Intestinal responses of the oriental fruit fly Bactrocera dorsalis (Hendel) after ingestion of an entomopathogenic bacterium strain

BACKGROUND

Bactrocera dorsalis (Hendel) is the main fruit fly pest of tropical and subtropical countries. The application of insecticides to manage this pest has led to serious resistance problems; therefore, new ways to control B. dorsalis are required. Pathogenic bacteria are sources of biocontrol agents for pest management.

RESULTS

We determined that a pathogenic bacterial strain, Serratia marcescens PS-1, isolated from a moribund striped flea beetle (Phyllotreta striolata), was lethal to B. dorsalis adults following ingestion. Histological analyses revealed that PS-1 damaged the intestinal epithelium, resulting in cell death within 24 h. We then generated a gut transcriptomic data set using RNA-Seq at two time points (6 and 24 h) after PS-1 infection. We found that genes encoding the peritrophic matrix constituent were down-regulated, whereas genes involved in lipid and glycan metabolism, and renewal of the gut epithelium, along with genes encoding digestive enzymes and stress response factors, were up-regulated. In addition, 14 cecropin genes were identified and cloned from B. dorsalis. To our knowledge, the number of cecropins identified in the present study is greater than that reported in the insects of earlier studies. Moreover, some of the cecropins identified were significantly down-regulated after PS-1 treatment.

CONCLUSION

Our findings provide new insights into the insect gut response to pathogenic bacterial invasion and may aid the development of new strategies for the biological control of B. dorsalis. © 2019 Society of Chemical Industry.

Differential Gene Expression for Age Estimation of Forensically Important Sarcophaga peregrina (Diptera: Sarcophagidae) Intrapuparial

Sarcophaga peregrina is an important flesh fly species for estimating the minimum postmortem interval (PMImin) in forensic entomology. The accurate determination of the developmental age is a crucial task for using necrophagous sarcophagids to estimate PMImin. During larval development, the age determination is straight forward by the morphological changes and variation of length, weight, and width; however, the age estimation of sarcophagid intrapuparial is more difficult due to anatomical and morphological changes not being visible. The analysis of differentially expressed genes (DEGs) during sarcophagid metamorphosis is a potential method for age estimation of intrapuparial. In the present study, real-time quantitative polymerase chain reaction (RT-qPCR) was used to analyze the differential gene expression level of S. peregrina intrapuparial in different constant temperatures (35°C, 25°C, and 15°C). In addition, the appropriate reference genes of S. peregrina were selected in the intrapuparial and at different temperatures to obtain reliable and valid gene expression profiles. The results indicated that two candidate genes (18S rRNA and 28S rRNA) were the most reliable reference genes, and four DEGs (Hsp90, A-alpha, AFP, AFBP) have the potential to be used to more accuracy estimate the age of S. peregrina intrapuparial.

Selection of valid reference genes for quantitative real-time PCR in Cotesia chilonis (Hymenoptera: Braconidae) exposed to different temperatures

In quantitative real-time PCR (qRT-PCR), data are normalized using reference genes, which helps to control for internal differences and reduce error among samples. In this study, the expression profiles of eight candidate housekeeping genes, 18S ribosomal (18S rRNA), elongation factor (EF1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein L17 (RPL17), histone 3 (H3), arginine kinase (AK), amd β-Actin (ACTB), were evaluated in the parasitic wasp Cotesia chilonis in response to different temperatures. Specifically, the performance and stabilities of these genes were compared in adult wasps maintained in a growth condition at 27°C (normal storage conditions) and in adults obtained from pupae refrigerated at 4°C for five days (cold storage conditions). Data were analyzed using the ΔCt method, BestKeeper, NormFinder, and geNorm. The optimal numbers and stabilities of reference genes varied between the two temperature treatments (27°C and 4°C). In samples stored at normal developmental temperature (27°C), the requirement for normalization in response to low temperature exposures was three genes (18S, H3, AK), whereas normalization in response to high temperature exposures required only two reference genes (GAPDH, ACTB). In samples stored at cold temperature (4°C), for low temperature exposures two reference genes (RPL17, RPL10) were required for standardization, while following high temperature exposures three reference genes (18S, H3, ACTB) were needed. This study strengthens understanding of the selection of reference genes before qRT-PCR analysis in C. chilonis. The reference genes identified here will facilitate further investigations of the biological characteristics of this important parasitoid.

Efficiency of RNA interference is improved by knockdown of dsRNA nucleases in tephritid fruit flies

RNA interference (RNAi) in insects is routinely used to ascertain gene function, but also has potential as a technology to control pest species. For some insects, such as beetles, ingestion of small quantities of double-stranded RNA (dsRNA) is able to knock down a targeted gene's expression. However, in other species, ingestion of dsRNA can be ineffective owing to the presence of nucleases within the gut, which degrade dsRNA before it reaches target cells. In this study, we observed that nucleases within the gut of the Queensland fruit fly (Bactrocera tryoni) rapidly degrade dsRNA and reduce RNAi efficacy. By complexing dsRNA with liposomes within the adult insect's diet, RNAi-mediated knockdown of a melanin synthesis gene, yellow, was improved significantly, resulting in strong RNAi phenotypes. RNAi efficiency was also enhanced by feeding both larvae and adults for several days on dsRNAs that targeted two different dsRNase gene transcripts. Co-delivery of both dsRNase-specific dsRNAs and yellow dsRNA resulted in almost complete knockdown of the yellow transcripts. These findings show that the use of liposomes or co-feeding of nuclease-specific dsRNAs significantly improves RNAi inhibition of gene expression in B. tryoni and could be a useful strategy to improve RNAi-based control in other insect species.

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.

CRISPR/Cas9-mediated knockout of the eye pigmentation gene white leads to alterations in colour of head spots in the oriental fruit fly, Bactrocera dorsalis

The intensely studied white gene is widely used as a genetic marker in Drosophila melanogaster. Here, we cloned and characterized the white gene in an important pest of the fruit industry, Bactrocera dorsalis, to understand its functional role in pigmentation. We obtained BdWhite knockout strains, based on the wild-type strain, using the CRISPR/Cas9 genome editing system, and found that mutants lost pigmentation in the compound eye and their black head spots. We then examined differences in the expression levels of genes associated with melanin pigmentation between mutants and the wild-type strain using quantitative reverse transcription PCR. We found that transcription levels of the Bd-yellow1 were lower in the head of mutants than in the wild-type strain, and there were no significant differences in expression of the other six genes between mutants and the wild type. Since yellow is critical for melanin biosynthesis (Heinze et al., Scientific Reports. 2017;7:4582), the lower levels of expression of Bd-yellow1 in mutants led to reduced dark pigmentation in head spots. Our results provide the first evidence, to our knowledge, that white may play a functional role in cuticle pigmentation by affecting the expression of yellow.

Evaluation of Appropriate Reference Genes For Investigating Gene Expression in Chlorops oryzae (Diptera: Chloropidae)

Reverse transcription quantitative polymerase chain reaction (PCR) has become an invaluable technique for analyzing gene expression in many insects. However, this approach requires the use of stable reference genes to normalize the data. Chlorops oryzae causes significant economic damage to rice crops throughout Asia. The lack of suitable reference genes has hindered research on the molecular mechanisms underlying many physiological processes of this species. In this study, we used quantitative real-time PCR to evaluate the expression of eight C. oryzae housekeeping genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH), β-actin (βACT), beta-tubulin (βTUB), Delta Elongation factor-1 (EF1δ), ribosomal protein S11 (RPS11), RPS15, C-terminal-Binding Protein (CtBP), and ribosomal protein 49 (RP49) in different developmental stages and tissues in both larvae and adults. We analyzed the data with four different software packages: geNorm, NormFinder, BestKeeper, and RefFinder and compared the results obtained with each method. The results indicate that PRS15 and RP49 can be used as stable reference genes for quantifying gene expression in different developmental stages and larval tissues. GAPDH and βACT, which have been considered stable reference genes by previous studies, were the least stable of the candidate genes with respect to larval tissues. GAPDH was, however, the most stable reference gene for adult tissues. We verified the candidate reference genes identified and found that the expression levels of Cadherins (Cads) changed when different reference genes were used to normalize gene expression. This study provides a valuable foundation for future research on gene function, and investigating the molecular basis of physiological processes, in C. oryzae.