Insect ATP-Binding Cassette (ABC) Transporters: Roles in Xenobiotic Detoxification and Bt Insecticidal Activity

Transfer and bioaccumulation of pesticides in terrestrial arthropods and food webs: State of knowledge and perspectives for research

Arthropods represent an entry point for pesticide transfers in terrestrial food webs, and pesticide accumulation in upper chain organisms, such as predators can have cascading consequences on ecosystems. However, the mechanisms driving pesticide transfer and bioaccumulation in food webs remain poorly understood. Here we review the literature on pesticide transfers mediated by terrestrial arthropods in food webs. The transfer of pesticides and their potential for bioaccumulation and biomagnification are related to the chemical properties and toxicokinetic of the substances, the resistance and detoxification abilities of the contaminated organisms, as well as by their effects on organisms' life history traits. We further identify four critical areas in which knowledge gain would improve future predictions of pesticides impacts on terrestrial food webs. First, efforts should be made regarding the effects of co-formulants and pesticides mixtures that are currently understudied. Second, progress in the sensitivity of analytical methods would allow the detection of low concentrations of pesticides in small individual arthropods. Quantifying pesticides in arthropods preys, their predators, and arthropods or vertebrates at higher trophic level would bring crucial insights into the bioaccumulation and biomagnification potential of pesticides in real-world terrestrial food webs. Finally, quantifying the influence of the trophic structure and complexity of communities on the transfer of pesticides could address several important sources of variability in bioaccumulation and biomagnification across species and food webs. This narrative review will inspire future studies aiming to quantify pesticide transfers in terrestrial food webs to better capture their ecological consequences in natural and cultivated landscapes.

The role of ATP-binding cassette transporters in arthropod pesticide toxicity and resistance

Pesticide resistance in arthropods threatens agricultural productivity and the control of vector-borne diseases. The ATP-binding cassette (ABC) transporters have emerged as important factors in the toxicity of synthetic pesticides, as well as for Bacillus thuringiensis insecticidal Cry protein binding. Depending on the localization of expression, both higher and lower expression of ABCs have been linked with pesticide resistance. The recent development of genetic-based approaches such as RNAi and CRISPR/Cas9 gene editing in nonmodel species, has greatly contributed to unveil their functional importance in pesticide toxicity and resistance. Using these tools, we are now poised to further unravel the molecular genetic mechanisms of gene regulation uncovering more elusive regulatory resistance genes.

Emerging Role of ABC Transporters in Glia Cells in Health and Diseases of the Central Nervous System

ATP-binding cassette (ABC) transporters play a crucial role for the efflux of a wide range of substrates across different cellular membranes. In the central nervous system (CNS), ABC transporters have recently gathered significant attention due to their pivotal involvement in brain physiology and neurodegenerative disorders, such as Alzheimer's disease (AD). Glial cells are fundamental for normal CNS function and engage with several ABC transporters in different ways. Here, we specifically highlight ABC transporters involved in the maintenance of brain homeostasis and their implications in its metabolic regulation. We also show new aspects related to ABC transporter function found in less recognized diseases, such as Huntington's disease (HD) and experimental autoimmune encephalomyelitis (EAE), as a model for multiple sclerosis (MS). Understanding both their impact on the physiological regulation of the CNS and their roles in brain diseases holds promise for uncovering new therapeutic options. Further investigations and preclinical studies are warranted to elucidate the complex interplay between glial ABC transporters and physiological brain functions, potentially leading to effective therapeutic interventions also for rare CNS disorders.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein

BACKGROUND

Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood.

RESULTS

To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes.

CONCLUSIONS

This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Chemical tolerance related to the ABC transporter gene and DNA methylation in cladocera (Daphnia magna)

We performed multigenerational tests to clarify the chemical tolerance mechanisms of a nontarget aquatic organism, Daphnia magna. We continuously exposed D. magna to a carbamate insecticide (pirimicarb) at lethal or sublethal concentrations (0, 3.8, 7.5, and 15 μg/L) for 15 generations (F0-F14). We then determined the 48 h-EC50 values and mRNA expression levels of acetylcholinesterase, glutathione S-transferase, and ATP (Adenosine triphosphate)-binding cassette transporter (ABCt) in neonates (<24 h old) from F0, F4, F9, and F14. To ascertain the effects of DNA methylation on pirimicarb sensitivity, we measured 5-methylcytosine levels (DNA methylation levels) in neonates of parents in the last generation (F14). In addition, we cultured groups exposed to 0 and 7.5 μg/L (the latter of which acquired chemical tolerance to pirimicarb) with or without 5-azacytidine (de-methylating agent) and determined methylation levels and 48 h-EC50 values in neonates (<24 h old) from the treated parents. The EC50 values (30.3-31.6 μg/L) in F14 of the 7.5 and 15 μg/L groups were approximately two times higher than that in the control (16.0 μg/L). A linear mixed model analysis showed that EC50 and ABCt mRNA levels were significantly increased with generational alterations; further analysis showed that the ABCt mRNA level was positively related to the EC50 . Therefore, ABCt may be associated with altered pirimicarb sensitivity. In addition, the EC50 value and DNA methylation levels in pirimicarb-tolerant clones decreased after exposure to 5-azacytidine, suggesting that DNA methylation contributes to chemical tolerance. These findings improved our knowledge regarding the acquisition of chemical tolerance in aquatic organisms.

The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics

BACKGROUND

There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches.

METHODS

Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime.

RESULTS

Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures.

CONCLUSION

This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

Population genomics of Agrotis segetum provide insights into the local adaptive evolution of agricultural pests

BACKGROUND

The adaptive mechanisms of agricultural pests are the key to understanding the evolution of the pests and to developing new control strategies. However, there are few studies on the genetic basis of adaptations of agricultural pests. The turnip moth, Agrotis segetum (Lepidoptera: Noctuidae) is an important underground pest that affects a wide range of host plants and has a strong capacity to adapt to new environments. It is thus a good model for studying the adaptive evolution of pest species.

RESULTS

We assembled a high-quality reference genome of A. segetum using PacBio reads. Then, we constructed a variation map of A. segetum by resequencing 98 individuals collected from six natural populations in China. The analysis of the population structure showed that all individuals were divided into four well-differentiated populations, corresponding to their geographical distribution. Selective sweep analysis and environmental association studies showed that candidate genes associated with local adaptation were functionally correlated with detoxification metabolism and glucose metabolism.

CONCLUSIONS

Our study of A. segetum has provided insights into the genetic mechanisms of local adaptation and evolution; it has also produced genetic resources for developing new pest management strategies.

Individual transmembrane domains of SfABCC2 from Spodoptera frugiperda do not serve as functional Cry1F receptors

The fall armyworm (Spodoptera frugiperda) is a major global pest causing severe damage to various crops, especially corn. Transgenic corn producing the Cry1F pesticidal protein from the bacterium Bacillus thuringiensis (Cry1F corn) showed effectiveness in controlling this pest until S. frugiperda populations at locations in North and South America evolved practical resistance. The mechanism for practical resistance involved disruptive mutations in an ATP binding cassette transporter subfamily C2 gene (SfABCC2), which serves as a functional Cry1F receptor in the midgut cells of susceptible S. frugiperda. The SfABCC2 protein contains two transmembrane domains (TMD1 and TMD2), each with a cytosolic nucleotide (ATP) binding domain (NBD1 and NBD2, respectively). Previous reports have demonstrated that disruptive mutations in TMD2 were linked with resistance to Cry1F, yet whether the complete SfABCC2 structure is needed for receptor functionality or if a single TMD-NBD protein can serve as functional Cry1F receptor remains unknown. In the present study, we separately expressed TMD1 and TMD2 with their corresponding NBDs in cultured insect cells and tested their Cry1F receptor functionality. Our results show that the complete SfABCC2 structure is required for Cry1F receptor functionality. Moreover, binding competition assays revealed that Cry1F specifically bound to SfABCC2, whereas neither SfTMD1-NBD1 nor SfTMD2-NBD2 exhibited any significant binding. These results provide insights into the molecular mechanism of Cry1F recognition by SfABCC2 in S. frugiperda, which could facilitate the development of more effective insecticidal proteins.

Gene expression plasticity facilitates different host feeding in Ips sexdentatus (Coleoptera: Curculionidae: Scolytinae)

Host shift is ecologically advantageous and a crucial driver for herbivore insect speciation. Insects on the non-native host obtain enemy-free space and confront reduced competition, but they must adapt to survive. Such signatures of adaptations can often be detected at the gene expression level. It is astonishing how bark beetles cope with distinct chemical environments while feeding on various conifers. Hence, we aim to disentangle the six-toothed bark beetle (Ips sexdentatus) response against two different conifer defences upon host shift (Scots pine to Norway spruce). We conducted bioassay and metabolomic analysis followed by RNA-seq experiments to comprehend the beetle's ability to surpass two different terpene-based conifer defence systems. Beetle growth rate and fecundity were increased when reared exclusively on spruce logs (alternative host) compared to pine logs (native host). Comparative gene expression analysis identified differentially expressed genes (DEGs) related to digestion, detoxification, transporter activity, growth, signalling, and stress response in the spruce-feeding beetle gut. Transporter genes were highly abundant during spruce feeding, suggesting they could play a role in pumping a wide variety of endogenous and xenobiotic compounds or allelochemicals out. Trehalose transporter (TRET) is also up-regulated in the spruce-fed beetle gut to maintain homeostasis and stress tolerance. RT-qPCR and enzymatic assays further corroborated some of our findings. Taken together, the transcriptional plasticity of key physiological genes plays a crucial role after the host shift and provides vital clues for the adaptive potential of bark beetles on different conifer hosts.

Differential gene expression and microRNA profile in corpora allata-corpora cardiaca of Aedes aegypti mosquitoes with weak juvenile hormone signalling

The corpora allata-corpora cardiaca (CA-CC) is an endocrine gland complex that regulates mosquito development and reproduction through the synthesis of juvenile hormone (JH). Epoxidase (Epox) is a key enzyme in the production of JH. We recently utilized CRISPR/Cas9 to establish an epoxidase-deficient (epox-/-) Aedes aegypti line. The CA from epox-/- mutants do not synthesize epoxidated JH III but methyl farneosate (MF), a weak agonist of the JH receptor, and therefore have reduced JH signalling. Illumina sequencing was used to examine the differences in gene expression between the CA-CC from wild type (WT) and epox-/- adult female mosquitoes. From 18,034 identified genes, 317 were significantly differentially expressed. These genes are involved in many biological processes, including the regulation of cell proliferation and apoptosis, energy metabolism, and nutritional uptake. In addition, the same CA-CC samples were also used to examine the microRNA (miRNA) profiles of epox-/- and WT mosquitoes. A total of 197 miRNAs were detected, 24 of which were differentially regulated in epox-/- mutants. miRNA binding sites for these particular miRNAs were identified using an in silico approach; they target a total of 101 differentially expressed genes. Our results suggest that a lack of epoxidase, besides affecting JH synthesis, results in the diminishing of JH signalling that have significant effects on Ae. aegypti CA-CC transcriptome profiles, as well as its miRNA repertoire.

Knockdown of the ABCG23 Gene Disrupts the Development and Lipid Accumulation of Panonychus citri (Acari/Tetranychidae)

Panonychus citri is a worldwide citrus pest that is currently controlled through the use of insecticides. However, alternative strategies are required to manage P. citri. Recent studies suggest that the ATP-binding cassette (ABC) transporter G subfamily plays a crucial role in transporting cuticular lipids, which are essential for the insect's barrier function against microbial penetration. Therefore, investigating the potential of the ABC transporter G subfamily as a control measure for P. citri could be a promising approach. Based on the genome database, the gene was cloned, and the transcriptional response of ABCG23 for the different developmental stages of P. citri and under spirobudiclofen stress was investigated. Our results showed that the expression level of ABCG23 was significantly lower in adult females exposed to treatment compared to the control and was higher in females than males. The knockdown of ABCG23 using RNAi led to a decrease in the survival rate, fecundity, and TG contents of P. citri. Additionally, a lethal phenotype was characterized by body wrinkling and darkening. These results indicate that ABCG23 may be involved in cuticular lipid transportation and have adverse effects on the development and reproduction of P. citri, providing insight into the discovery of new targets for pest management based on the insect cuticle's penetration barrier function.

Analysis of Genes Associated with Feeding Preference and Detoxification in Various Developmental Stages of Aglais urticae

Herbivorous insects and host plants have developed a close and complex relationship over a long period of co-evolution. Some plants provide nutrients for insects, but plants' secondary metabolites also influence their growth and development. Urtica cannabina roots and leaves are poisonous, yet Aglais urticae larvae feed on them, so we aimed to clarify the mechanism enabling this interaction. At present, studies on the detoxification mechanism of the A. urticae are rare. In our study, first, we used the A. urticae larval odor selection behavior bioassay and choice feeding preference assay to analyze the feeding preferences of A. urticae on its host plant, U. cannabina. Next, we used transcriptome sequencing to obtain the unigenes annotated and classified by various databases, such as KEGG and GO. In this study, we found that U. cannabina could attract A. urticae larvae to feed via scent, and the feeding preference assay confirmed that larvae preferred U. cannabina leaves over three other plants: Cirsium japonicum, Cannabis sativa, and Arctium lappa. The activity of detoxifying enzymes GST and CarE changed in larvae that had consumed U. cannabina. Furthermore, through transcriptomic sequencing analysis, 77,624 unigenes were assembled from raw reads. The numbers of differentially expressed genes were calculated using pairwise comparisons of all life stages; the expression of detoxification enzyme genes was substantially higher in larvae than in the pupal and adult stages. Finally, we identified and summarized 34 genes associated with detoxification enzymes, such as UDP-glucose 4-epimerase gene, 5 Glutathione S-transferase genes, 4 Carboxylesterase genes, 4 Cytochrome P450 genes, 10 ATP-binding cassette genes, 4 Superoxide dismutase, and Peroxidase. Moreover, we identified 28 genes associated with the development of A. urticae. The qRT-PCR results were nearly consistent with the transcriptomic data, showing an increased expression level of four genes in larvae. Taken together, this study examines the correlation between A. urticae and host plants U. cannabina, uncovering a pronounced preference for A. urticae larvae toward host plants. Consistent with RNA-seq, we investigated the mechanism of A. urticae's interaction with host plants and identified detoxification-related genes. The present study provides theoretical support for studying insect adaptation mechanisms and biological control.

Overexpression of ABCB transporter genes confer multiple insecticide tolerances in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

Bactrocera dorsalis is a notable invasive pest that has developed resistance to several commonly used insecticides in the field, such as avermectin, beta-cypermethrin and malathion. Investigating the mechanisms of insecticide resistance in this pest is of paramount importance for ensuring its effective control. The ATP-binding cassette transporter subfamily B (ABCB) genes, responsible for encoding transmembrane efflux transporters, represent a potential source of insecticide detoxification activity or transportation that remains largely unexplored in B. dorsalis. In this study, seven BdABCB genes were identified and comprehensive analyzed based on the latest genome and transcriptome dataset. Subsequently, we characterized the expression profiles of these genes across different development stages and tissues, as well as under different insecticide exposures. The results showed that the BdABCB genes were expressed at all stages in B. dorsalis, with BdABCB2 and BdABCB7 being highly expressed in the pupal stage, while BdABCB5 and BdABCB6 were highly expressed in the larval stage. Besides, the BdABCBs were highly expressed in the detoxification metabolic tissues. Among them, BdABCB5 and BdABCB6 were significantly overexpressed in the midgut and Malpighian tubules, respectively. Furthermore, with the exception of BdABCB6, the expression levels of the other six BdABCBs were significantly up-regulated following induction with avermectin, beta-cypermethrin and malathion. Six BdABCBs (BdABCB1-5 and BdABCB7) were knocked down by RNA interference, and the interference efficiencies were 46.58%, 39.50%, 45.60%, 33.74%, 66.37% and 63.83%, respectively. After injecting dsBdABCBs, the mortality of flies increased by 25.23% to 39.67% compared to the control upon exposure to the three insecticides. These results suggested that BdABCBs play crucial roles in the detoxification or tolerance of B. dorsalis to multiple insecticides.

Time-series analysis of transcriptomic changes due to permethrin exposure reveals that Aedes aegypti undergoes detoxification metabolism over 24 h

Insecticide resistance is a multifaceted response and an issue across taxa. Aedes aegypti, the mosquito that vectors Zika, dengue, chikungunya, and yellow fever, demonstrates high levels of pyrethroid resistance across the globe, presenting a challenge to public health officials. To examine the transcriptomic shifts across time after exposure to permethrin, a 3'Tag-Seq analysis was employed on samples 6, 10, and 24 h after exposure along with controls. Differential expression analysis revealed significant shifts in detoxifying enzymes and various energy-producing metabolic processes. These findings indicate significant alterations in gene expression associated with key energy mobilization pathways within the system. These changes encompass a coordinated response involving lipolysis, beta-oxidation, and the citric acid cycle, required for the production of energetic molecules such as ATP, NADH, NADPH, and FADH. These findings highlight a complex interplay of metabolic processes that may have broader implications for understanding insect physiology and response to environmental stimuli. Among the upregulated detoxifying enzymes are cytochrome P450s, glutathione s-transferases and peroxidases, and ATP-binding cassette transporters. Additionally, eight heat shock genes or genes with heat shock domains exhibit the highest fold change across time. Twenty-four hours after exposure, samples indicate a global downregulation of these processes, though principal component analysis suggests lasting signatures of the response. Understanding the recovery response to insecticide exposure provides information on possible new genetic and synergist targets to explore.

Assessment of the role of an ABCC transporter TuMRP1 in the toxicity of abamectin to Tetranychus urticae

The rapid evolution of pest resistance threatens the sustainable utilization of bioinsecticides such as abamectin, and so deciphering the molecular mechanisms affecting toxicity and resistance is essential for their long-term application. Historical studies of abamectin resistance in arthropods have mainly focused on mechanisms involving the glutamate-gated chloride channel (GluCl) targets, with the role of metabolic processes less clear. The two-spotted spider mite, Tetranychus urticae, is a generalist herbivore notorious for rapidly developing resistance to pesticides worldwide, and abamectin has been widely used for its control in the field. After reanalyzing previous transcriptome and RNA-seq data, we here identified an ABC transporter subfamily C gene in T. urticae named multidrug resistance-associated protein 1 (TuMRP1), whose expression differed between susceptible and resistant populations. Synergism bioassays with the inhibitor MK-571, the existence of a genetic association between TuMRP1 expression and susceptibility to abamectin, and the effect of RNA interference mediated silencing of TuMRP1 were all consistent with a direct role of this transporter protein in the toxicity of abamectin. Although ABC transporters are often involved in removing insecticidal compounds from cells, our data suggest either an alternative role for these proteins in the mechanism of action of abamectin or highlight an indirect association between their expression and abamectin toxicity.

Natural Coumarin Shows Toxicity to Spodoptera litura by Inhibiting Detoxification Enzymes and Glycometabolism

Coumarin and its derivatives are plant-derived compounds that exhibit potent insecticidal properties. In this study, we found that natural coumarin significantly inhibited the growth and development of Spodoptera litura larvae through toxicological assay. By transcriptomic sequencing, 80 and 45 differentially expressed genes (DEGs) related to detoxification were identified from 0 to 24 h and 24 to 48 h in S. litura after coumarin treatment, respectively. Enzyme activity analysis showed that CYP450 and acetylcholinesterase (AChE) activities significantly decreased at 48 h after coumarin treatment, while glutathione S-transferases (GST) activity increased at 24 h. Silencing of SlCYP324A16 gene by RNA interference significantly increased S. litura larval mortality and decreased individual weight after treatment with coumarin. Additionally, the expression levels of DEGs involved in glycolysis and tricarboxylic acid (TCA) cycle were inhibited at 24 h after coumarin treatment, while their expression levels were upregulated at 48 h. Furthermore, metabonomics analysis identified 391 differential metabolites involved in purine metabolism, amino acid metabolism, and TCA cycle from 0 to 24 h after treated with coumarin and 352 differential metabolites associated with ATP-binding cassette (ABC) transporters and amino acid metabolism. These results provide an in-depth understanding of the toxicological mechanism of coumarin on S. litura.

The silkworm (Bombyx mori) gut microbiota is involved in metabolic detoxification by glucosylation of plant toxins

Herbivores have evolved the ability to detoxify feed components through different mechanisms. The oligophagous silkworm feeds on Cudrania tricuspidata leaves (CTLs) instead of mulberry leaves for the purpose of producing special, high-quality silk. However, CTL-fed silkworms are found to have smaller bodies, slower growth and lower silk production than those fed mulberry leaves. Here, we show that the high content of prenylated isoflavones (PIFs) that occurred in CTLs is converted into glycosylated derivatives (GPIFs) in silkworm faeces through the silkworm gut microbiota, and this biotransformation is the key process in PIFs detoxification because GPIFs are found to be much less toxic, as revealed both in vitro and in vivo. Additionally, adding Bacillus subtilis as a probiotic to remodel the gut microbiota could beneficially promote silkworm growth and development. Consequently, this study provides meaningful guidance for silk production by improving the adaptability of CTL-fed silkworms.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Identification and Characterization of Glutathione S-transferase Genes in Spodoptera frugiperda (Lepidoptera: Noctuidae) under Insecticides Stress

Insect glutathione S-transferases (GSTs) serve critical roles in insecticides and other forms of xenobiotic chemical detoxification. The fall armyworm, Spodoptera frugiperda (J. E. Smith), is a major agricultural pest in several countries, especially Egypt. This is the first study to identify and characterize GST genes in S. frugiperda under insecticidal stress. The present work evaluated the toxicity of emamectin benzoate (EBZ) and chlorantraniliprole (CHP) against the third-instar larvae of S. frugiperda using the leaf disk method. The LC₅₀ values of EBZ and CHP were 0.029 and 1.250 mg/L after 24 h of exposure. Moreover, we identified 31 GST genes, including 28 cytosolic and 3 microsomal SfGSTs from a transcriptome analysis and the genome data of S. frugiperda. Depending on the phylogenetic analysis, sfGSTs were divided into six classes (delta, epsilon, omega, sigma, theta, and microsomal). Furthermore, we investigated the mRNA levels of 28 GST genes using qRT-PCR under EBZ and CHP stress in the third-instar larvae of S. frugiperda. Interestingly, SfGSTe10 and SfGSTe13 stood out with the highest expression after the EBZ and CHP treatments. Finally, a molecular docking model was constructed between EBZ and CHP using the most upregulated genes (SfGSTe10 and SfGSTe13) and the least upregulated genes (SfGSTs1 and SfGSTe2) of S. frugiperda larvae. The molecular docking study showed EBZ and CHP have a high binding affinity with SfGSTe10, with docking energy values of -24.41 and -26.72 kcal/mol, respectively, and sfGSTe13, with docking energy values of -26.85 and -26.78 kcal/mol, respectively. Our findings are important for understanding the role of GSTs in S. frugiperda regarding detoxification processes for EBZ and CHP.

Chlorbenzuron caused growth arrest through interference of glycolysis and energy metabolism in Hyphantria cunea (Lepidoptera: Erebidae) larvae

Chlorbenzuron is a kind of benzoylphenylureas (BPUs), which plays a broad role in insect growth regulators (IGRs), with an inhibitory effect on chitin biosynthesis. However, BPUs how to regulate glycolysis and insect growth remains largely unclear. Here, we investigated the effects of chlorbenzuron on growth, nutritional indices, glycolysis, and carbohydrate homeostasis in Hyphantria cunea, a destructive and highly polyphagous forest pest, to elucidate the action mechanism of chlorbenzuron from the perspective of energy metabolism. The results showed that chlorbenzuron dramatically restrained the growth and nutritional indices of H. cunea larvae and resulted in lethality. Meanwhile, we confirmed that chlorbenzuron significantly decreased carbohydrate levels, adenosine triphosphate (ATP), and pyruvic acid (PA) in H. cunea larvae. Further studies indicated that chlorbenzuron caused a significant enhancement in the enzyme activities and mRNA expressions of hexokinase (HK), phosphofructokinase (PFK), and pyruvate kinase (PK), resulting in increased glycolytic flux. Expressions of genes involved in the AMP-activated protein kinase (AMPK) signaling pathway were also upregulated. Moreover, chlorbenzuron had remarkable impacts on H. cunea larvae from the perspective of metabolite enrichment, including the tricarboxylic acid (TCA) cycle and glycolysis, indicating an energy metabolism disorder in larvae. The findings provide a novel insight into the molecular mechanism by which chlorbenzuron abnormally promotes glycolysis and eventually interferes with insect growth and nutritional indices.

Expression profile of the entire detoxification gene inventory of the western honeybee, Apis mellifera across life stages

The western honeybee, Apis mellifera, is a managed pollinator of many crops and potentially exposed to a wide range of foreign compounds, including pesticides throughout its life cycle. Honeybees as well as other insects recruit molecular defense mechanisms to facilitate the detoxification of xenobiotic compounds. The inventory of detoxification genes (DETOXome) is comprised of five protein superfamilies: cytochrome P450 monooxygenases (P450), carboxylesterases, glutathione S-transferases (GST), UDP-glycosyl transferases (UGT) and ATP-binding cassette (ABC) transporters. Here we characterized the gene expression profile of the entire honeybee DETOXome by analyzing 47 transcriptomes across the honeybee life cycle, including different larval instars, pupae, and adults. All life stages were well separated by principal component analysis, and K-means clustering revealed distinct temporal patterns of gene expression. Indeed, >50% of the honeybee detoxification gene inventory is found in one cluster and follows strikingly similar expression profiles, i.e., increased expression during larval development, followed by a sharp decline after pupation and a steep increase again in adults. This cluster includes 29 P450 genes dominated by CYP3 and CYP4 clan members, 15 ABC transporter genes mostly belonging to the ABCC subfamily and 13 carboxylesterase genes including almost all members involved in dietary/detox and hormone/semiochemical processing. RT-qPCR analysis of selected detoxification genes from all families revealed high expression levels in various tissues, especially Malpighian tubules, fatbody and midgut, supporting the view that these tissues are essential for metabolic clearance of environmental toxins and pollutants in honeybees. Our study is meant to spark further research on the molecular basis of detoxification in this critical pollinator to better understand and evaluate negative impacts from potentially toxic substances. Additionally, the entire gene set of 47 transcriptomes collected and analyzed provides a valuable resource for future honeybee research across different disciplines.

Preliminary Analysis of Transcriptome Response of Dioryctria sylvestrella (Lepidoptera: Pyralidae) Larvae Infected with Beauveria bassiana under Short-Term Starvation

The Dioryctria genus contains several destructive borer pests that are found in coniferous forests in the Northern Hemisphere. Beauveria bassiana spore powder was tested as a new method of pest control. In this study, Dioryctria sylvestrella (Lepidoptera: Pyralidae) was used as the object. A transcriptome analysis was performed on a freshly caught group, a fasting treatment control group, and a treatment group inoculated with a wild B. bassiana strain, SBM-03. Under the conditions of 72-h fasting and a low temperature of 16 ± 1 °C, (i) in the control group, 13,135 of 16,969 genes were downregulated. However, in the treatment group, 14,558 of 16,665 genes were upregulated. (ii) In the control group, the expression of most genes in the upstream and midstream of the Toll and IMD pathways was downregulated, but 13 of the 21 antimicrobial peptides were still upregulated. In the treatment group, the gene expression of almost all antimicrobial peptides was increased. Several AMPs, including cecropin, gloverin, and gallerimycin, may have a specific inhibitory effect on B. bassiana. (iii) In the treatment group, one gene in the glutathione S-transferase system and four genes in the cytochrome P450 enzyme family were upregulated, with a sharp rise in those that were upregulated significantly. In addition, most genes of the peroxidase and catalase families, but none of the superoxide dismutase family were upregulated significantly. Through innovative fasting and lower temperature control, we have a certain understanding of the specific defense mechanism by which D. sylvestrella larvae may resist B. bassiana in the pre-wintering period. This study paves the way for improving the toxicity of B. bassiana to Dioryctria spp.

Identification and characterization of ABC proteins in an important rice insect pest, Cnaphalocrocis medinalis unveil their response to Cry1C toxin

Rice leaffolder (Cnaphalocrocis medinalis) is an important insect pest in paddy fields. Due to their essential role in the physiology and insecticidal resistance, ATP-binding cassette (ABC) proteins were studied in many insects. In this study, we identified the ABC proteins in C. medinalis through genomic data and analyzed their molecular characteristics. A total of 37 sequences with nucleotide-binding domain (NBD) were identified as ABC proteins and belonged to eight families (ABCA-ABCH). Four structure styles of ABC proteins were found in C. medinalis, including full structure, half structure, single structure, and ABC2 structure. In addition to these structures, TMD-NBD-TMD, NBD-TMD-NBD, and NBD-TMD-NBD-NBD were found in C. medinalis ABC proteins. Docking studies suggested that in addition to the soluble ABC proteins, other ABC proteins including ABCC4, ABCH1, ABCG3, ABCB5, ABCG1, ABCC7, ABCB3, ABCA3, and ABCC5 binding with Cry1C had higher weighted scores. The upregulation of ABCB1 and downregulation of ABCB3, ABCC1, ABCC7, ABCG1, ABCG3, and ABCG6 were associated with the C. medinalis response to Cry1C toxin. Collectively, these results help elucidate the molecular characteristics of C. medinalis ABC proteins, pave the way for further functional studies of C. medinalis ABC proteins, including their interaction with Cry1C toxin, and provide potential insecticide targets.

Computational Analysis of the Ligand-Binding Sites of the Molecular Chaperone OppA from Yersinia pseudotuberculosis

The function of chaperones is to correct or degrade misfolded proteins inside the cell. Classic molecular chaperones such as GroEL and DnaK have not been found in the periplasm of Yersinia pseudotuberculosis. Some periplasmic substrate-binding proteins could be bifunctional, such as OppA. Using bioinformatic tools, we try to elucidate the nature of the interactions between OppA and ligands from four proteins with different oligomeric states. Using the crystal structure of the proteins Mal12 alpha-glucosidase from Saccharomyces cerevisiae S288C, LDH rabbit muscle lactate dehydrogenase, EcoRI endonuclease from Escherichia coli and THG Geotrichum candidum lipase, a hundred models were obtained in total, including five different ligands from each enzyme with five conformations of each ligand. The best values for Mal12 stem from ligands 4 and 5, with conformation 5 for both; for LDH, ligands 1 and 4, with conformations 2 and 4, respectively; for EcoRI, ligands 3 and 5, with conformation 1 for both; and for THG, ligands 2 and 3, with conformation 1 for both. The interactions were analyzed with LigProt, and the length of the hydrogen bridges has an average of 2.8 to 3.0 Å. The interaction within the OppA pocket is energetically favored due to the formation of hydrogen bonds both of OppA and of the selected enzymes. The Asp 419 residue is important in these junctions.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Transcriptional response and functional analysis of ATP-binding cassette transporters to tannic acid in pea aphid, Acyrthosiphon pisum (Harris)

Although tannins are widely distributed in broad beans and alfalfa, the pea aphid (Acyrthosiphon pisum) can still destroy them. The ATP binding cassette (ABC) transporters participate in the metabolism of plant secondary metabolites and pesticides in insects. However, whether ABC transporter genes play a role in the metabolism of tannins in the pea aphid is unclear. Here, we found that verapamil (an ABC transporter inhibitor) significantly increased the mortality of tannic acid to pea aphid, which indicated that ABC transporter gene was related to the metabolism of tannic acid by pea aphid. Then, we identified 54 putative ABC transporter genes from the genome database of A. pisum. These genes were divided into eight subfamilies, ApABCA to ApABCH, of which subfamily G has the largest number of genes with 19, followed by the subfamily C with 14. RT-qPCR results show that the expression levels of ApABCA2, ApABCC7, ApABCG2, and ApABCG3 were highly expressed in the first instar, while those of ApABCA3, ApABCG6, ApABCG7, ApABCH3, and ApABCH4 were highly expressed in adults. Furthermore, transcription levels of many ABC transporter genes were induced by tannic acid. Especially, ApABCG17 and ApABCH2 were obviously induced after being exposed to tannic acid. Meanwhile, knockdown of ApABCG17 by RNA interference resulted in increased sensitivity of pea aphid to tannic acid. These results suggest that ApABCG17 may be involved in tannic acid metabolism in pea aphid. This study will help us to understand the mechanism of tannic acid metabolism in pea aphid, and provides a basis for further research on the physiological function of ABC transporter genes in pea aphid.

Cyclosporin A acts as a novel insecticide against Cry1Ac-susceptible and -resistant Helicoverpa armigera

Cotton bollworm (Helicoverpa armigera) is an economically important pest, which is difficult to manage due to its biological and ecological traits, and resistance to most insecticides. Alternative compounds for the sustainable management of H. armigera are needed. As a fungal metabolite, Cyclosporin A (CsA) has not been applied in agriculture pests. Here, CsA was evaluated as a propective insecticide for H. armigera. The results showed that CsA displayed high insecticidal activity against both Cry1Ac-susceptible and -resistant populations of H. armigera. Moreover, lower concentrations of CsA had clear effects, including significantly reduced pupal weight, pupation rate, emergence rate, ovary size, female fecundity and egg hatchability. Further study confirmed that CsA suppressed calcineurin activity and the subsequent expression of endogenous antimicrobial peptide genes (APMs), leading to impaired immunity, ultimately resulting in delayed development and increased mortality. Thus, CsA treatment could control the cotton bollworm population and even showed efficacy against those with Bt resistance. In addition, the morphological changes observed in insects fed CsA with lower concentrations provide insight into insect immunity, regulation of growth and development, regulation of body color, ovary development and sexual selection under external pressure. Overall, our study provides information on biological control potential of Cry1Ac-susceptible and -resistant populations of H. armigera to develop novel bioinsecticides.

Comparative genomic analysis of ABC transporter genes in Tenebrio molitor and four other tenebrionid beetles (Coleoptera: Tenebrionidea)

ATP-binding cassette (ABC) transporters, one of the largest transmembrane protein families, transport a diverse number of substate across membranes. Details of their diverse physiological functions have not been established. Here, we identified 87 ABC transporter genes in the genomes of Tenebrio molitor along with those from Asbolus verrucosus (104), Hycleus cichorii (65), and Hycleus phaleratus (80). Combining these genes (336 in total) with genes reported in Tribolium castaneum (73), we analyzed the phylogeny of ABC transporter genes in all five Tenebrionids. They are assigned into eight subfamilies (ABCA-H). In comparison to other species, the ABCC subfamily in this group of beetles appears expanded. The expression profiles of the T. molitor genes at different life stages and in various tissues were also investigated using transcriptomic analysis. Most of them display developmental specific expression patterns, suggesting to us their possible roles in development. Most of them are highly expressed in detoxification-related tissues including gut and Malpighian tubule, from which we infer their roles in insecticide resistance. We detected specific or abundant expressions of many ABC transporter genes in various tissues such as salivary gland, ovary, testis, and antenna. This new information helps generate new hypotheses on their biological significance within tissues.

Cyclosporin A as a Potential Insecticide to Control the Asian Corn Borer Ostrinia furnacalis Guenée (Lepidoptera: Pyralidae)

The long-term use of chemical insecticides has caused serious problems of insect resistance and environmental pollution; new insecticides are needed to solve this problem. Cyclosporin A (CsA) is a polypeptide produced by many fungi, which is used to prevent or treat immune rejection during organ transplantation. However, little is known about the utility of CsA as an insecticide. Therefore, this study evaluated the insecticidal activity of CsA using Ostrinia furnacalis as a model. The results demonstrated that CsA was toxic to O. furnacalis with LC₅₀ values of 113.02 μg/g and 198.70 μg/g for newly hatched neonates and newly molted third-instar larvae, respectively. Furthermore, CsA treatment had sublethal effects on the development of O. furnacalis, and significantly reduced the fecundity of adults; this suggests that CsA has great potential to suppress O. furnacalis populations. Further analysis revealed that CsA suppressed calcineurin activity in larvae. CsA had independent or synergistic toxic effects on O. furnacalis when combined with β-cypermethrin, indoxacarb, emamectin benzoate, azadirachtin, and the Bacillus thuringiensis toxin Cry1Ac, which suggests that CsA can help prevent or manage resistance. Our study provides detailed information on the potential of CsA as an insecticide for controlling lepidopterans.

A single transcription factor facilitates an insect host combating Bacillus thuringiensis infection while maintaining fitness

Maintaining fitness during pathogen infection is vital for host survival as an excessive response can be as detrimental as the infection itself. Fitness costs are frequently associated with insect hosts countering the toxic effect of the entomopathogenic bacterium Bacillus thuringiensis (Bt), which delay the evolution of resistance to this pathogen. The insect pest Plutella xylostella has evolved a mechanism to resist Bt toxins without incurring significant fitness costs. Here, we reveal that non-phosphorylated and phosphorylated forms of a MAPK-modulated transcription factor fushi tarazu factor 1 (FTZ-F1) can respectively orchestrate down-regulation of Bt Cry1Ac toxin receptors and up-regulation of non-receptor paralogs via two distinct binding sites, thereby presenting Bt toxin resistance without growth penalty. Our findings reveal how host organisms can co-opt a master molecular switch to overcome pathogen invasion with low cost, and contribute to understanding the underlying mechanism of growth-defense tradeoffs during host-pathogen interactions in P. xylostella.

Metabolization and sequestration of plant specialized metabolites in insect herbivores: Current and emerging approaches

Herbivorous insects encounter diverse plant specialized metabolites (PSMs) in their diet, that have deterrent, anti-nutritional, or toxic properties. Understanding how they cope with PSMs is crucial to understand their biology, population dynamics, and evolution. This review summarizes current and emerging cutting-edge methods that can be used to characterize the metabolic fate of PSMs, from ingestion to excretion or sequestration. It further emphasizes a workflow that enables not only to study PSM metabolism at different scales, but also to tackle and validate the genetic and biochemical mechanisms involved in PSM resistance by herbivores. This review thus aims at facilitating research on PSM-mediated plant-herbivore interactions.

High-quality genomes reveal significant genetic divergence and cryptic speciation in the model organism Folsomia candida (Collembola)

The collembolan Folsomia candida Willem, 1902, is widely distributed throughout the world and has been frequently used as a test organism in soil ecology and ecotoxicology studies. However, it is questioned as an ideal "standard" because of differences in reproductive modes and cryptic genetic diversity between strains from various geographical origins. In this study, we obtained two high-quality chromosome-level genomes of F. candida, for a parthenogenetic strain (named as FCDK, 219.08 Mb, 25,139 protein-coding genes) and a sexual strain (named as FCSH, 153.09 Mb, 21,609 protein-coding genes), reannotated the genome of the parthenogenetic strain reported by Faddeeva-Vakhrusheva et al. in 2017 (named as FCBL, 221.7 Mb, 25,980 protein-coding genes), and conducted comparative genomic analyses of three strains. High genome similarities between FCDK and FCBL on synteny, genome architecture, mitochondrial and nuclear gene sequences support they are conspecific. The seven chromosomes of FCDK are each 25-54% larger than the corresponding chromosomes of FCSH, showing obvious repetitive element expansions and large-scale inversions and translocations but no whole-genome duplication. The strain-specific genes, expanded gene families and genes in nonsyntenic chromosomal regions identified in FCDK are highly related to the broader environmental adaptation of parthenogenetic strains. In addition, FCDK has fewer strain-specific microRNAs than FCSH, and their mitochondrial and nuclear genes have diverged greatly. In conclusion, FCDK/FCBL and FCSH have accumulated independent genetic changes and evolved into distinct species since 10 Mya. Our work provides important genomic resources for studying the mechanisms of rapidly cryptic speciation and soil arthropod adaptation to soil ecosystems.

Novel-miR-310 mediated response mechanism to Cry1Ac protoxin in Plutella xylostella (L.)

The diamondback moth (DBM), Plutella xylostella (L.), has evolved resistance to multiple insecticides including Bacillus thuringiensis (Bt). ATP-binding cassette (ABC) transporters are a class of transmembrane protein families, involved in multiple physiological processes and pesticide resistances in insects. However, the role and regulatory mechanism of ABC transporter in mediating the response to Bt Cry1Ac toxin remain unclear. Here, we characterized a MAPK signaling pathway-enriched ABCG subfamily gene PxABCG20 from DBM, and found it was differentially expressed in the Cry1Ac-resistant and Cry1Ac-susceptible strains. RNAi knockdown of PxABCG20 increased the tolerance of DBM to Cry1Ac protoxin. To explore the regulatory mechanism of PxABCG20 expression, we predicted the potential miRNAs targeting PxABCG20 using two target prediction algorithms. Luciferase reporter assay confirmed that novel-miR-310 was able to down-regulate PxABCG20 expression in HEK293T cells. Furthermore, injection of novel-miR-310 agomir markedly inhibited PxABCG20 expression, resulting in increased tolerance to Cry1Ac protoxin in susceptible strain, while injection of novel-miR-310 antagomir markedly induced the expression of PxABCG20, leading to decreased tolerance to Cry1Ac protoxin. Our work provides theoretical basis for exploring novel targets for the DBM response to Cry1Ac toxin and expands the understanding of miRNA role in mediating the susceptibility of insect pest to Cry1Ac toxin.

Searching and identifying pigmentation genes from Neocaridina denticulate sinensis via comparison of transcriptome in different color strains

Aquaria species are characterized by their amazing colors and patterns. Research on the breeding molecular genetics of ornamental shrimps is surprisingly limited. We conducted a transcriptome analysis to investigate the expression of encoding genes in the integument of the strains Neocaridina denticulate sinensis. After assembled and filtered, 19,992 unigenes were annotated by aligning with public functional databases (NR, Swiss-Prot, KEGG, COG). 14,915 unigenes with significantly different expressions were found by comparing three strains integument transcriptomes. Ribosomal protein genes, ABC transporter families, calmodulin, carotenoid proteins and crustacyanin may play roles in the cytological process of pigment migration and chromatophore maintenance. Numerous color genes associated with multiple pathways including melanin, ommochrome and pteridines pathways were identified. The expression patterns of 25 candidate genes were analysis by qPCR in red, yellow, transparent and glass strains. The qPCR results in red, yellow and transparent were consistent with the level of RPKM values in the transcriptomes. The above results will advance our knowledge of integument color varieties in N. denticulate sinensis and help the genetic selection of crustaceans with consumer-favored colors. Furthermore, it also provides some candidate pigmentation genes to investigate the correlation between coloration and sympatric speciation in crustaceans.

Differential expression profiling of Oxycarenus laetus Kirby (Hemiptera: Lygaeidae) upon exposure to gossypol

BACKGROUND

Gossypium hirsutum seeds are rich in gossypol. In addition to its diverse beneficial properties, it is a known anti-fertility inducing agent in humans. Oxycarenus laetus feeds on the cottonseeds and yet its courtship, mating and reproduction is unaffected.

METHODS AND RESULTS

In this study, we performed a transcriptome profiling of O. laetus fed on Abutilon indicum (AB-no gossypol), G. hirsutum (GH-natural gossypol) and 1400 ppm commercial gossypol-soaked GH seeds (GHGO). Illumina NextSeq-500 paired-end 75 bp reads were generated and de novo assembled (48,214 genes) to identify the differentially expressed transcripts (DET) between the samples. Gene enrichment, KEGG pathway and cluster profiling of the DETs resulted in the identification of vital genes involved in the detoxification, pheromone biosynthesis, cuticle protein in the GHGO sample. Cyp4C1, Cyp6a13, Cyp6a14, Cyp4g15, Cyp4em8, Cyp303a1 were the detoxification related genes identified. Similarly, SDR dehydrogenase family 11 and fatty acid synthase in pheromone biosynthesis and cuticle proteins (RR1 and RR2) coding transcripts were found to be differentially expressed.

CONCLUSION

This is the first study to report the expression of genes induced by gossypol in O. laetus. Based on the findings from the DET analysis, we conclude that the detoxification related genes of gossypol treated samples were affected.

Genome-wide identification and comparative analysis of Cry toxin receptor families in 7 insect species with a focus on Spodoptera litura

Cadherin, aminopeptidase N (APN) and alkaline phosphatase (ALP) have been characterized as Cry receptors. In this study, comparative genomic analysis of the 3 receptor families was performed in 7 insects. ALPs and APNs are divided into three and eight clades in phylogenetic trees, respectively. ALPs in clade 3 and APNs in clade 1 contain multiple paralogs within each species and most paralogs are located closely in chromosomes. Drosophila melanogaster has expanded APNs in clade 5 and were lowly expressed in midgut. Cadherins are divided into 16 clades; they may diverge before holometabolous insect speciation except for BtR and Cad89D-like clades. Eight insects from different orders containing BtR orthologs are sensitive to Cry1A or Cry3A, while five species without BtR are insensitive to both toxins. Most APNs in clade 1, several ALPs in clade 3, BtR and Cad89D-like genes were highly or moderately expressed in larval midgut of Spodoptera litura and the other six species, and several members in these clades have been identified as Cry receptors. Expressions of putative S. litura Cry receptors in the midgut after exposing to Bt toxins were also analyzed.

Map-based cloning and functional analysis revealed ABCC2 is responsible for Cry1Ac toxin resistance in Bombyx mori

Bt toxins are parasporal crystals produced by Bacillus thuringiensis (Bt). They have specific killing activity against various insects and have been widely used to control agricultural pests. However, their widespread use has developed the resistance of many target insects. To maintain the sustainable use of Bt products, the resistance mechanism of insects to Bt toxins must be fully clarified. In this study, Bt-resistant and Bt-susceptible silkworm strains were used to construct genetic populations, and the genetic pattern of silkworm resistance to Cry1Ac toxin was determined. Sequence-tagged site molecular marker technology was used to finely map the resistance gene and to draw a molecular genetic linkage map, and the two closest markers were T1590 and T1581, indicating the resistance gene located in the 155 kb genetic region. After analyzing the sequence of the predicted gene in the genetic region, an ATP binding cassette transporter (ABCC2) was identified as the candidate gene. Molecular modeling and protein-protein docking result showed that a tyrosine insertion in the mutant ABCC2 might be responsible for the interaction between Cry1Ac and ABCC2. Moreover, CRISPR/Cas9-mediated genome editing technology was used to knockout ABCC2 gene. The homozygous mutant ABCC2 silkworm was resistant to Cry1Ac toxin, which indicated ABCC2 is the key gene that controls silkworm resistance to Cry1Ac toxin. The results have laid the foundation for elucidating the molecular resistance mechanism of silkworms to Cry1Ac toxin and could provide a theoretical basis for the biological control of lepidopteran pests.

Transcriptional Analysis of Cotton Bollworm Strains with Different Genetic Mechanisms of Resistance and Their Response to Bacillus thuringiensis Cry1Ac Toxin

Transgenic crops producing Bacillus thuringiensis (Bt) insecticidal proteins are grown widely for pest control, but the evolution of resistance in target pests could reduce their efficacy. Mutations in genes encoding cadherin, ABC transporter or tetraspanin were linked with resistance to Cry1Ac in several lepidopteran insects, including the cotton bollworm (Helicoverpa armigera), a worldwide agricultural pest. However, the detailed molecular mechanisms by which these mutations confer insect resistance to Cry1Ac remain largely unknown. In this study, we analyzed the midgut transcriptomes of a susceptible SCD strain and three SCD-derived Cry1Ac-resistant strains of H. armigera (SCD-r1, with a naturally occurring deletion mutation of cadherin; SCD-KI, with a knock-in T92C point mutation in tetraspanin; and C2/3-KO, with both ABCC2 and ABCC3 knocked out). Evaluation of midgut transcript profiles of the four strains without Cry1Ac exposure identified many constitutively differentially expressed genes (DEGs) in the resistant SCD-r1 (n = 1355), SCD-KI (n = 1254) and C2/3-KO (n = 2055) strains. Analysis of DEGs in the midguts of each strain after Cry1Ac exposure revealed similar patterns of response to Cry1Ac in the SCD and SCD-r1 strains, but unique responses in the SCD-KI and C2/3-KO strains. Expression of midgut epithelium healing and defense-related genes was strongly induced by Cry1Ac intoxication in the SCD and SCD-r1 strains, while immune-related pattern recognition receptor and effector genes were highly expressed in the SCD-KI strain after Cry1Ac exposure. This study advances our knowledge of the transcriptomic basis for insect resistance to Bt toxins and provides a valuable resource for further molecular characterization of insect response to Cry1Ac toxin in H. armigera and other pest species.

Genomic features of the polyphagous cotton leafworm Spodoptera littoralis

BACKGROUND

The cotton leafworm, Spodoptera littoralis, is a highly polyphagous pest of many cultivated plants and crops in Africa and Europe. The genome of this pest will help us to further understand the molecular mechanisms of polyphagy.

RESULTS

Herein, the high-quality genome of S. littoralis was obtained by Pacific Bioscience (PacBio) sequencing. The assembled genome size of S. littoralis is 436.55 Mb with a scaffold N50 of 6.09 Mb, consisting of 17,207 annotated protein-coding genes. Phylogenetic analysis shows that S. littoralis and its sibling species S. litura diverged about 5.44 million years ago. Expanded gene families were mainly involved in metabolic detoxification and tolerance to toxic xenobiotics based on GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis. Comparative genomics analysis showed that gene families involved in detoxification and chemosensation were significantly expanded in S. littoralis, representing genetic characteristics related to polyphagy and an extensive host range.

CONCLUSIONS

We assembled and annotated the reference genome of S. littoralis, and revealed that this pest has the genetic features of strong detoxification capacity, consistent with it being a significant risk to a wide range of host crops. These data resources will provide support for risk assessment and early warning monitoring of major polyphagous agricultural pests.

Identification and Validation of ATP-Binding Cassette Transporters Involved in the Detoxification of Abamectin in Rice Stem Borer, Chilo suppressalis

Chilo suppressalis has developed high levels of resistance to abamectin in many areas of China, while the underline resistance mechanisms are largely unclear. ATP-binding cassette (ABC) transporters function in transporting a large diversity of substrates including insecticides and play important roles in the detoxification metabolism of insects. In this study, synergism bioassay revealed that the ABC transporters were involved in the detoxification of C. suppressalis to abamectin. Six ABC transporter genes were upregulated in C. suppressalis after abamectin exposure, among which five genes CsABCC8, CsABCE1, CsABCF1, CsABCF2, and CsABCH1 were induced in the detoxification-related tissues. In addition, the five ABC transporters were recombinantly expressed in Sf9 cells, and the cytotoxicity assay showed that the viabilities of cells expressing CsABCC8 or CsABCH1 were significantly increased when compared with the viabilities of cells expressing EGFP after abamectin, chlorantraniliprole, cyantraniliprole, fipronil, and chlorpyrifos treatment, respectively. Overexpression of CsABCE1 significantly increased the viabilities of cells to abamectin, chlorantraniliprole, deltamethrin, and indoxacarb exposure, respectively. These results suggested that CsABCC8, CsABCE1, and CsABCH1 might participate in the detoxification and transport of abamectin and several other classes of insecticides in C. suppressalis. Our study provides valuable insights into the transport-related detoxification mechanisms in C. suppressalis and other insects.

The Mycotoxin Beauvericin Exhibits Immunostimulatory Effects on Dendritic Cells via Activating the TLR4 Signaling Pathway

Beauvericin (BEA), a mycotoxin of the enniatin family produced by various toxigenic fungi, has been attributed multiple biological activities such as anti-cancer, anti-inflammatory, and anti-microbial functions. However, effects of BEA on dendritic cells remain unknown so far. Here, we identified effects of BEA on murine granulocyte–macrophage colony-stimulating factor (GM-CSF)-cultured bone marrow derived dendritic cells (BMDCs) and the underlying molecular mechanisms. BEA potently activates BMDCs as signified by elevated IL-12 and CD86 expression. Multiplex immunoassays performed on myeloid differentiation primary response 88 (MyD88) and toll/interleukin-1 receptor (TIR) domain containing adaptor inducing interferon beta (TRIF) single or double deficient BMDCs indicate that BEA induces inflammatory cytokine and chemokine production in a MyD88/TRIF dependent manner. Furthermore, we found that BEA was not able to induce IL-12 or IFNβ production in Toll-like receptor 4 (Tlr4)-deficient BMDCs, whereas induction of these cytokines was not compromised in Tlr3/7/9 deficient BMDCs. This suggests that TLR4 might be the functional target of BEA on BMDCs. Consistently, in luciferase reporter assays BEA stimulation significantly promotes NF-κB activation in mTLR4/CD14/MD2 overexpressing but not control HEK-293 cells. RNA-sequencing analyses further confirmed that BEA induces transcriptional changes associated with the TLR4 signaling pathway. Together, these results identify TLR4 as a cellular BEA sensor and define BEA as a potent activator of BMDCs, implying that this compound can be exploited as a promising candidate structure for vaccine adjuvants or cancer immunotherapies.

Gut Metagenomic Profiling of Gossypol Induced Oxycarenus laetus (Hemiptera: Lygaeidae) Reveals Gossypol Tolerating Bacterial Species

Oxycarenus laetus is a cotton pest that primarily feeds on seeds that are rich in gossypol. Though gossypol is toxic to general herbivores, O. laetus does not show such complications and instead grows and reproduces well on cotton plants compared to its other hosts. In this study, we have fed O. laetus with natural and induced gossypol-based diets to explore the difference in its gut microbiota. We performed NGS 16S rRNA amplicon sequencing on the Illumina MiSeq platform and analyzed the data using the QIIME2 pipeline supplemented with Greengenes and EZBioCloud reference databases. We also used culture-based methods to identify a few abundant gut bacteria present in O. laetus. Enterococcus faecalis, Wolbachia bourtzisii, Wolbachia pipientis, Corynebacterium glyciniphilum, Staphylococcus sciuri, and Kocuria rosea were some of the major species that formed the core gut microbiome of O. laetus. We have also observed that some species were present only in the sample with the highest concentration of gossypol, signifying that they might have the potential to degrade gossypol.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12088-021-00964-0.

MAP4K4 controlled transcription factor POUM1 regulates PxABCG1 expression influencing Cry1Ac resistance in Plutella xylostella (L.)

Deciphering the molecular mechanisms of insect resistance to Bacillus thuringiensis (Bt) based biotechnology products including Bt sprays and Bt crops is critical for the long-term application of Bt technology. Previously, we established that down-regulation of the ABC transporter gene PxABCG1, trans-regulated by the MAPK signaling pathway, contributed to high-level resistance to Bt Cry1Ac toxin in diamondback moth, Plutella xylostella (L.). However, the underlying transcriptional regulatory mechanism was unknown. Herein, we identified putative binding sites (PBSs) of the transcription factor (TF) POUM1 in the PxABCG1 promoter and used a dual-luciferase reporter assay (DLRA) and yeast one-hybrid (Y1H) assay to reveal that POUM1 activates PxABCG1 via interaction with one of these sites. The expression of POUM1 was significantly decreased in the midgut tissue of Cry1Ac-resistant P. xylostella strains compared to a Cry1Ac-susceptible P. xylostella strain. Silencing of POUM1 expression resulted in reduced expression of the PxABCG1 gene and an increase in larval tolerance to Bt Cry1Ac toxin in the Cry1Ac-susceptible P. xylostella strain. Furthermore, silencing of PxMAP4K4 expression increased the expression of both POUM1 and PxABCG1 genes in the Cry1Ac-resistant P. xylostella strain. These results indicate that the POUM1 induces PxABCG1 expression, while the activated MAPK cascade represses PxABCG1 expression thus reducing Cry1Ac susceptibility in P. xylostella. This result deepens our understanding of the transcriptional regulatory mechanism of midgut Cry receptor genes and the molecular basis of the evolution of Bt resistance in insects.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

Chromosome-Scale Genome of Masked Palm Civet (Paguma larvata) Shows Genomic Signatures of Its Biological Characteristics and Evolution

The masked palm civet (Paguma larvata) is a small carnivore with distinct biological characteristics, that likes an omnivorous diet and also serves as a vector of pathogens. Although this species is not an endangered animal, its population is reportedly declining. Since the severe acute respiratory syndrome (SARS) epidemic in 2003, the public has been particularly concerned about this species. Here, we present the first genome of the P. larvata, comprising 22 chromosomes assembled using single-tube long fragment read (stLFR) and Hi-C technologies. The genome length is 2.41 Gb with a scaffold N50 of 105.6 Mb. We identified the 107.13 Mb X chromosome and one 1.34 Mb Y-linked scaffold and validated them by resequencing 45 P. larvata individuals. We predicted 18,340 protein-coding genes, among which 18,333 genes were functionally annotated. Interestingly, several biological pathways related to immune defenses were found to be significantly expanded. Also, more than 40% of the enriched pathways on the positively selected genes (PSGs) were identified to be closely related to immunity and survival. These enriched gene families were inferred to be essential for the P. larvata for defense against the pathogens. However, we did not find a direct genomic basis for its adaptation to omnivorous diet despite multiple attempts of comparative genomic analysis. In addition, we evaluated the susceptibility of the P. larvata to the SARS-CoV-2 by screening the RNA expression of the ACE2 and TMPRSS2/TMPRSS4 genes in 16 organs. Finally, we explored the genome-wide heterozygosity and compared it with other animals to evaluate the population status of this species. Taken together, this chromosome-scale genome of the P. larvata provides a necessary resource and insights for understanding the genetic basis of its biological characteristics, evolution, and disease transmission control.

Identification and Functional Analysis of Differentially Expressed Genes in Myzus persicae (Hemiptera: Aphididae) in Response to Trans-anethole

Plant essential oils, with high bioactivity and biodegradability, provide promising alternatives to synthetic pesticides for pest control. Trans-anethole is the major component of essential oil from star anise, Illicium verum Hook. The compound has a strong contact toxicity against the green peach aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae), which is a major insect pest of many vegetables and crops. However, little information is known about how M. persicae responds to trans-anethole at the molecular level. We conducted a comparative transcriptome analysis of M. persicae in response to a LD50 dose of trans-anethole. A total of 559 differentially expressed genes were detected in the treated individuals, with 318 genes up-regulated, and 241 genes down-regulated. Gene ontology (GO) analysis revealed that these genes were classified into different biological processes and pathways. We also found that genes encoding ATP-binding cassette (ABC) transporters, DnaJ, and cuticle proteins were dramatically up-regulated in response to trans-anethole. To study the function of these genes, we performed RNA interference (RNAi) analysis. Knockdown of an ABC transporter gene (ABCG4) and a DnaJ gene (DnaJC1) resulted in a significantly increased mortality rate in M. persicae following trans-anethole exposure, indicating the involvement of these two genes in the toxicity response to trans-anethole. The findings provide new insights into the mechanisms of M. persicae in coping with plant essential oils.

Engineering of Cry3Bb1 provides mechanistic insights toward countering western corn rootworm resistance

•

Cry3Bb1 engineering for receptor retargeting to counter WCR resistance described

•

Phage display screens against BBMV or recombinant insecticidal-protein receptors

•

77 gut binding peptides selected and engineered into 284 new variants of Cry3Bb1

•

112 variants were active against susceptible but not resistant WCR

•

Engineering of Cry3Bb1 Domain II loops 1 and 2 disrupted insecticidal activity

•

Sites for Cry3Bb1 engineering and implications for Cry3Bb1 resistance discussed

The western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae), is an economically important pest of corn (maize) in North America and Europe. Current management practices for WCR involve transgenic expression of insecticidal proteins to minimize larval feeding damage to corn roots. The evolution of resistant WCR populations to transgenic corn expressing insecticidal proteins (e.g. Cry3Bb1, Gpp34Ab1/Tpp35Ab1) necessitates efforts to discover and deploy new modes of action for WCR control. Here, we tested the hypothesis that the addition of short peptides selected for binding to the WCR gut would restore insecticidal activity of Cry3Bb1 to resistant insects. Phage display technology coupled with deep sequencing was used to identify peptides selected for binding to WCR brush border membrane vesicles and to recombinant putative receptors aminopeptidase and cadherin. The binding and specificity of selected peptides was confirmed by ELISA and pull-down assays, and candidate gut surface binding partners were identified. Although production of 284 novel Cry3Bb1 variants with these peptides did not restore activity against resistant WCR in artificial diet bioassays, 112 variants were active against susceptible insects. These results provided insights for the mechanism of Cry3Bb1 activity and toward engineering a new mode-of-action via receptor re-targeting in the context of protein structure and function.

ABCC2 is a functional receptor of Bacillus thuringiensis Cry1Ca in Spodoptera litura

Spodoptera litura is a serious polyphagous pest in the whole world, which has developed resistance to most conventional insecticides and even some Bacillus thuringiensis (Bt) toxins. Cry1Ca has excellent insecticide activity against S. litura with potential application to control S. litura and delay the development of insect resistance. However, the mode of action of Cry1Ca in S. litura is poorly understood. Here, Cry1Ca-binding proteins were identified from S. litura by using pull down assays and liquid chromatography-tandem mass spectrometry (LC-MS/MS). The results indicated that aminopeptidase-N (APN), ATP binding cassette subfamily C member 2 (ABCC2), polycalin, actin and V-type proton ATPase subunit A may bind with Cry1Ca. Further study confirmed that ABCC2 fragment expressed in vitro can bind to Cry1Ca as demonstrated by Ligand blot and homologous competition experiments. The over-expression of endogenous SlABCC2 in Sf9 cells increased Cry1Ca cytotoxicity. Correspondingly, the vivo loss of function analyses by SlABCC2 small interfering RNAs (siRNAs) in S. litura larvae decreased the toxicity of Cry1Ca to larvae. Altogether, these results show that ABCC2 of S. litura is a functional receptor that is involved in the action mode of Cry1Ca.

Silencing of cytochrome P450 gene CYP321A1 effects tannin detoxification and metabolism in Spodoptera litura

Cytochrome P450 monooxygenase (P450 or CYP) plays an important role in the metabolism of insecticides and plant allelochemicals by insects. CYP321B1, a novel Spodoptera litura P450 gene, was identified and characterized. CYP321B1 contains a 1488 bp open reading frame (ORF) that encodes a 495 amino acid protein. In fourth instar larvae, the highest CYP321B1 expression levels were found in the midgut and fat body. In the tannin feeding test, tannin can significantly induce the expression of CYP321B1 in the midgut and fat body of 4th instar larvae. To verify the function of CYP321B1, RNA interference and metabolome analysis were performed. The results showed that silencing CYP321B1 significantly reduced the rate of weight gain under tannin induction. Metabolome analysis showed silencing affected 47 different metabolites, mainly involved in secondary metabolite biosynthesis and amino acid metabolism, including amino acids, lipid fatty acids, organic acids and their derivatives. Henoxyacetic acid and cysteamine are the most highly regulated metabolites, respectively. These findings demonstrate that CYP321B1 plays an important role in tannin detoxification and metabolism. Functional knowledge about metabolite detoxification genes in this major herbivorous insect pest can provide new insights into this biological process and provide new targets for agricultural pest control.