Comparative transcriptome analysis of thiamethoxam susceptible and resistant Asian citrus psyllid, Diaphorina citri (Hemiptera: Liviidae), using RNA-sequencing

Gene expression changes in Maconellicoccus hirsutus in response to sublethal dose of buprofezin

The pink or hibiscus mealybug, Maconellicoccus hirsutus, is a serious pest of grapes, jute, and mesta, causing severe yield losses in India and other countries. Chemical control remains the foremost choice for farmers to manage this pest. As insecticides break down over time due to biotic and abiotic factors, insects are exposed to varying levels of these exogenous compounds. Several studies have reported that sublethal doses affect insect physiology, but only a few have examined the changes in gene expression at the molecular level. Therefore, the present study was conducted to elucidate the molecular mechanisms in M. hirsutus exposed to sublethal doses of buprofezin 25 SC. Life table analysis revealed increased fecundity in M. hirsutus exposed to the sublethal dose. A total of 1,744 differentially expressed genes were identified between the buprofezin-treated and untreated samples using transcriptome analysis. These genes were primarily associated with ribosomal proteins, proteases, cuticular proteins, and cytoskeletal structures. Ribosomes and phagosomes were the most highly enriched pathways. Interestingly, most of the DEGs were involved in restoring homeostasis rather than detoxification. To validate our RNA-sequencing results, qRT-PCR validation was performed on ten randomly selected genes. Overall, our findings provide valuable insights into intermittent changes in stress-coping genes, apart from detoxification genes.

Citrus-mediated gene silencing of cytochrome P450 suppresses insecticide resistance and increases mortality in Diaphorina citri

BACKGROUND

Asian citrus psyllid, Diaphorina citri, is a hemipteran that vectors the causal pathogen of citrus greening disease, or huanglongbing (HLB). HLB is a tree killing disease that has severely limited citrus production globally. Unfortunately, there is no cure for this disease, and mitigation depends on multiple insecticide applications to reduce vector populations. Silencing of cytochrome P450 expression associated with detoxification enzymes by RNA interference is known to increase susceptibility of D. citri to insecticides. However, dsRNA was previously introduced into psyllids by topical applications. The possible application of this technology for pest management will require effective field delivery of the dsRNA. Therefore, we evaluated a virus vector (Citrus tristeza virus; 'mild strain' T36) to deliver gene silencing directly to this sap-sucking insect via plant phloem. Citrus macrophylla plants inoculated with CTV expressing a truncated consensus sequence of CYP450 (CTV-tCYP450) constantly produced small interfering RNA in the plant phloem that targeted five cytochrome p540 (CYP450) genes in D. citri.

RESULTS

Insecticide susceptible D. citri reared on citrus infected with CTV-tCYP450 were subsequently more susceptible to imidacloprid, fenpropathrin, carbaryl, and chlorpyrifos than those reared on citrus infected with wildtype CTV or non-infected negative controls. Additionally, nymph survival and adult lifespan were significantly reduced when psyllids were reared on CTV-tCYP450 citrus plants compared with controls. Interestingly, similar results were obtained after one and two generations of rearing. Finally, field-collected psyllids from areas with known broad-spectrum insecticide resistance were rendered more susceptible to imidacloprid and fenpropathrin after feeding on CTV-tCYP450 citrus trees as compared with those reared on controls.

CONCLUSION

The integration of citrus-mediated RNA inference targeting psyllid detoxification enzymes could function as a resistance management tool and reduce insecticide input in an integrated pest management program for HLB. © 2024 Society of Chemical Industry.

Exposure to Sublethal Concentrations of Dinotefuran Induces Apoptosis in the Gut of Diaphorina citri Adults via Activating the Mitochondrial Apoptotic Pathway

Diaphorina citri is a serious citrus pest. Dinotefuran is highly insecticidal against D. citri. To analyze the sublethal effects of dinotefuran on D. citri adults, an indoor toxicity test was performed, which revealed that the lethal concentration 50 (LC50) values were 4.23 and 0.50 μg/mL for 24 and 48 h treatments, respectively. RNA-Seq led to the identification of 71 and 231 differentially expressed genes (DEGs) after dinotefuran treatments with LC20 and LC50 doses, respectively. Many of the DEGs are significantly enriched in the apoptosis pathway. Dinotefuran-induced apoptosis in the gut cells was confirmed through independent assays of 4',6-diamidino-2-phenylindole (DAPI) and TdT-mediated dUTP nick end labeling (TUNEL) staining. Increased levels of reactive oxygen species (ROS) and a loss of mitochondrial membrane potential were observed. Four caspase genes were identified, and dinotefuran treatments resulted in increased mRNA levels of DcCasp1 and DcCasp3a. These findings shed light on the sublethal effects of dinotefuran on D. citri.

Adverse effects of thiamethoxam on the behavior, biochemical responses, hepatopancreas health, transcriptome and intestinal flora of juvenile Chinese mitten crab (Eriocheir sinensis)

Frequent detection of thiamethoxam in global surface waters has provoked great concern in environmental safety, as thiamethoxam exhibits high toxicity to aquatic arthropods. However, little systematic investigation has been conducted on the chronic toxicity of thiamethoxam to crustaceans. This study exposed Eriocheir sinensis to thiamethoxam (0, 0.5, 5 and 50 μg/L) in water for 28 days. No significant difference in mortality was observed among all groups. A high concentration of thiamethoxam (50 μg/L) impaired the righting ability of E. sinensis. Thiamethoxam significantly increased antioxidant enzyme activities (superoxide dismutase, total antioxidant capacity and glutathione peroxidase) and malondialdehyde levels. Simultaneously, detoxification enzyme activities (aminopyrine N-demethylase, erythromycin N-demethylase and glutathione-S-transferase) increased under chronic thiamethoxam stress. In addition, thiamethoxam caused immune and hepatopancreas damage. Moreover, thiamethoxam induced intestinal flora dysbiosis by altering the microbiome structure. The reduced complexity of the gut microbiota further illustrated that thiamethoxam could disrupt the stability of the microbiota ecological network. The transcriptomic results revealed that the number of downregulated DEGs increased in a dose-dependent manner, and most downregulated DEGs were enriched in energy metabolism-related pathways. These results indicate that thiamethoxam can adversely affect the crab behavior, biochemistry, intestinal microflora and transcriptomic responses.

The effect of chlorogenic acid, a potential botanical insecticide, on gene transcription and protein expression of carboxylesterases in the armyworm (Mythimna separata)

Chlorogenic acid (CGA) is a potential botanical insecticide metabolite that naturally occurs in various plants. Our previous studies revealed CGA is sufficient to control the armyworm Mythimna separata. In this study, we conducted a proteomic analysis of saliva collected from M. separata following exposure to CGA and found that differentially expressed proteins (DEPs) treated with CGA for 6 h and 24 h were primarily enriched in glutathione metabolism and the pentose phosphate pathway. Notably, we observed six carboxylesterase (CarE) proteins that were enriched at both time points. Additionally, these corresponding genes were expressed at levels 5.05 to 130.25 times higher in our laboratory-selected resistance strains. We also noted a significant increase in the enzyme activity of carboxylesterase following treatments with varying CGA concentrations. Finally, we confirmed that knockdown of MsCarE14, MsCarE28, and MsCCE001h decreased the susceptibility to CGA in resistance strain, indicating three CarE genes play crucial roles in CGA detoxification. This study presents the first report on the salivary proteomics of M. separata, offering valuable insights into the role of salivary proteins. Moreover, the determination of CarE mediated susceptibility change to CGA provides new targets for agricultural pest control and highlights the potential insecticide resistance mechanism for pest resistance management.

Diapause-Linked Gene Expression Pattern and Related Candidate Duplicated Genes of the Mountain Butterfly Parnassius glacialis (Lepidoptera: Papilionidae) Revealed by Comprehensive Transcriptome Profiling

The mountain butterfly Parnassius glacialis is a representative species of the genus Parnassius, which probably originated in the high-altitude Qinhai-Tibet Plateau in the Miocene and later dispersed eastward into relatively low-altitude regions of central to eastern China. However, little is known about the molecular mechanisms underlying the long-term evolutionary adaptation to heterogeneous environmental conditions of this butterfly species. In this study, we obtained the high-throughput RNA-Seq data from twenty-four adult individuals in eight localities, covering nearly all known distributional areas in China, and firstly identified the diapause-linked gene expression pattern that is likely to correlate with local adaptation in adult P. glacialis populations. Secondly, we found a series of pathways responsible for hormone biosynthesis, energy metabolism and immune defense that also exhibited unique enrichment patterns in each group that are probably related to habitat-specific adaptability. Furthermore, we also identified a suite of duplicated genes (including two transposable elements) that are mostly co-expressed to promote the plastic responses to different environmental conditions. Together, these findings can help us to better understand this species' successful colonization to distinct geographic areas from the western to eastern areas of China, and also provide us with some insights into the evolution of diapause in mountain Parnassius butterfly species.

Functional Characterization of the Ryanodine Receptor Gene in Diaphorina citri

The Asian citrus psyllid Diaphorina citri (Hemiptera: Liviidae) is a major citrus pest spread around the world. It is also a vector of the bacterium 'Candidatus Liberibacter asiaticus', considered the cause of the fatal citrus disease huanglongbing (HLB). Insect ryanodine receptors (RyRs) are the primary target sites of diamide insecticides. In this study, full-length RyR cDNA from D. citri (named DcRyR) was isolated and identified. The 15,393 bp long open reading frame of DcRyR encoded a 5130 amino acid protein with a calculated molecular weight of 580,830 kDa. This protein had a high sequence identity (76-79%) with other insect homologs and a low sequence identity (43-46%) with mammals. An MIR domain, two RIH domains, three SPRY domains, four RyR repeat domains, an RIH-associated domain at the N-terminus, two consensus calcium-binding EF-hands, and six transmembrane domains were among the characteristics that DcRyR shared with insect and vertebrate RyRs. In expression analysis, the DcRyR gene displayed transcript abundance in all tissues and developmental stages as well as gene-differential and stage-specific patterns. In addition, diagnostic PCR experiments revealed that DcRyR had three potential alternative splice variants and that splicing events might have contributed to the various functions of DcRyR. However, diamide resistance-related amino acid residue mutations I4790M/K and G4946E were not found in DcRyR. These results can serve as the basis for further investigation into the target-based diamide pesticide resistance of D. citri.

Transcriptome analysis of Mythimna separata: De novo assembly and detection of genes related to beta-cypermethrin resistance

The oriental armyworm Mythimna separata (Walker) is a devastating pest of cereal crops mainly in Asia and Oceania and recently become resistant to beta-cypermethrin (beta-CP). However, molecular biological studies of its response to beta-CP are scarce, and related genomic information is not available. In this study, we sequenced and de novo assembled the transcriptome of beta-CP susceptible M. separata (MsS-whole, abbr. MsS-W). A total of 30,486 unigenes were obtained, with an N50 length of 506 bp. A number of 21,051 unigenes were matched to public databases, of which 10,107 were classified into 59 gene ontology subcategories, 5792 were assigned into 25 clusters of orthologous groups of proteins subgroups and 12,123 were assigned to 257 Kyoto Encyclopedia of Genes and Genomes pathways. A total of 729 simple sequence repeats were detected. In addition, a total of 323 cytochrome P450-associated sequences from nine lepidopterous species, of which 130 were from M. separata, were analyzed using the maximum likelihood method and Bayesian inference. Among the 130 cytochrome P450-associated sequences from M. separata, 60 were dropped into CYP3 clan, which is associated with metabolizing xenobiotics and plant natural compounds. Furthermore, the beta-CP susceptible (MsS-2) and resistant (MsR-2) M. separata population transcriptomes were sequenced. Certain critical genes involved in beta-CP detoxification were detected and verified by quantitative real-time polymerase chain reaction. Collectively, our results provided a basis for further studies of the molecular mechanism of insecticide resistance in M. separata.

Silencing of Glycogen Synthase Kinase 3 Significantly Inhibits Chitin and Fatty Acid Metabolism in Asian Citrus Psyllid, Diaphorina citri

Glycogen is a predominant carbohydrate reserve in various organisms, which provides energy for different life activities. Glycogen synthase kinase 3 (GSK3) is a central player that catalyzes glucose and converts it into glycogen. In this study, a GSK3 gene was identified from the D. citri genome database and named DcGSK3. A reverse transcription quantitative PCR (RT-qPCR) analysis showed that DcGSK3 was expressed at a high level in the head and egg. The silencing of DcGSK3 by RNA interference (RNAi) led to a loss-of-function phenotype. In addition, DcGSK3 knockdown decreased trehalase activity, glycogen, trehalose, glucose and free fatty acid content. Moreover, the expression levels of the genes associated with chitin and fatty acid synthesis were significantly downregulated after the silencing of DcGSK3. According to a comparative transcriptomics analysis, 991 differentially expressed genes (DEGs) were identified in dsDcGSK3 groups compared with dsGFP groups. A KEGG enrichment analysis suggested that these DEGs were primarily involved in carbon and fatty acid metabolism. The clustering analysis of DEGs further confirmed that chitin and fatty acid metabolism-related DEGs were upregulated at 24 h and were downregulated at 48 h. Our results suggest that DcGSK3 plays an important role in regulating the chitin and fatty acid metabolism of D. citri.

A Novel Reference for Bt-Resistance Mechanism in Plutella xylostella Based on Analysis of the Midgut Transcriptomes

Simple Summary

Plutella xylostella is a very serious pest to cruciferous vegetables. At present, the control methods used are mainly traditional insecticides and the cultivation of Bt crops. However, with the long-term and large-scale use of insecticides, the diamondback moth has developed strong resistance to many kinds of insecticides and Bt crops. The Cry1S1000 strain of P. xylostella used here is a strain with more than 8000 times resistance to Bt Cry1Ac protoxin. In this paper, we used transcriptome sequencing to determine the midgut transcriptome of the G88-susceptible strain, Cry1S1000-resistant strain and its corresponding toxin-induced strains to find more genes related to Bt resistance. Our results will provide a reference for optimizing the control strategy of diamondback moth resistance and improving the control efficiency of biopesticides and Bt crops.

Abstract

The diamondback moth, Plutella xylostella, is a lepidopteran insect that mainly harms cruciferous vegetables, with strong resistance to a variety of agrochemicals, including Bacillus thuringiensis (Bt) toxins. This study intended to screen genes associated with Bt resistance in P. xylostella by comparing the midgut transcriptome of Cry1Ac-susceptible and -resistant strains together with two toxin-treated strains 24 h before sampling. A total of 12 samples were analyzed by BGISEQ-500, and each sample obtained an average of 6.35 Gb data. Additionally, 3284 differentially expressed genes (DEGs) were identified in susceptible and resistant strains. Among them, five DEGs for cadherin, 14 for aminopeptidase, zero for alkaline phosphatase, 14 for ATP binding cassette transport, and five heat shock proteins were potentially involved in resistance to Cry1Ac in P. xylostella. Furthermore, DEGs associated with “binding”, “catalytic activity”, “cellular process”, “metabolic process”, and “cellular anatomical entity” were more likely to be responsible for resistance to Bt toxin. Thus, together with other omics data, our results will offer prospective genes for the development of Bt resistance, thereby providing a brand new reference for revealing the resistance mechanism to Bt of P. xylostella.