Identification of cytochrome P450 monooxygenase genes and their expression profiles in cyhalothrin-treated Colorado potato beetle, Leptinotarsa decemlineata

Identification and expression dynamics of CYPome across different developmental stages of Maconellicoccus hirsutus (Green)

Maconellicoccus hirsutus is a highly polyphagous insect pest, posing a substantial threat to various crop sp., especially in the tropical and sub-tropical regions of the world. While extensive physiological and biological studies have been conducted on this pest, the lack of genetic information has hindered our understanding of the molecular mechanisms underlying its growth, development, and xenobiotic metabolism. The Cytochrome P450 gene, a member of the CYP gene superfamily ubiquitous in living organisms is associated with growth, development, and the metabolism of both endogenous and exogenous substances, contributing to the insect's adaptability in diverse environments. To elucidate the specific role of the CYP450 gene family in M. hirsutus which has remained largely unexplored, a de novo transcriptome assembly of the pink mealybug was constructed. A total of 120 proteins were annotated as CYP450 genes through homology search of the predicted protein sequences across different databases. Phylogenetic studies resulted in categorizing 120 CYP450 genes into four CYP clans. A total of 22 CYP450 families and 30 subfamilies were categorized, with CYP6 forming the dominant family. The study also revealed five genes (Halloween genes) associated with the insect hormone biosynthesis pathway. Further, the expression of ten selected CYP450 genes was studied using qRT-PCR across crawler, nymph, and adult stages, and identified genes that were expressed at specific stages of the insects. Thus, the findings of this study reveal the expression dynamics and possible function of the CYP450 gene family in the growth, development, and adaptive strategies of M. hirsutus which can be further functionally validated.

Sublethal effects of halofenozide on larval development and detoxification in Phaedon brassicae (Coleoptera: Chrysomelidae)

The brassica leaf beetle, Phaedon brassicae, is a serious defoliator of cruciferous crops. Halofenozide (Hal), an ecdysone agonist, is a new class of insect growth-regulating insecticide. Our preliminary experiment revealed the outstanding larval toxicity of Hal against P. brassicae. However, the metabolic degradation of this compound in insects remains unclear. In this study, oral administration of Hal at LC10 and LC25 caused severe separation of the cuticle and epidermis, leading to larval molting failure. Sublethal dose exposure also significantly reduced the larval respiration rate as well as their pupation rates and pupal weights. Conversely, the activities of the multifunctional oxidase, carboxylesterase (CarE), and glutathione S-transferase (GST) were significantly enhanced in Hal-treated larvae. Further analysis using RNA sequencing identified 64 differentially expressed detoxifying enzyme genes, including 31 P450s, 13 GSTs, and 20 CarEs. Among the 25 upregulated P450s, 22 genes were clustered into the CYP3 clan, and the other 3 genes belonged to the CYP4 clan. Meanwhile, 3 sigma class GSTs and 7 epsilon class GSTs were dramatically increased, accounting for the majority of the upregulated GSTs. Moreover, 16 of the 18 overexpressed CarEs were clustered into the coleopteran xenobiotic-metabolizing group. These results showed the augmented expression of detoxification genes in P. brassicae after exposed to sublethal dose of Hal, and helped to better understand the potential metabolic pathways that could contribute to the reduced sensitivity to Hal in this pest. Overall, a deep insight into the detoxification mechanisms would provide practical guidance for the field management of P. brassicae.

Functional roles of two novel P450 genes in the adaptability of Conogethes punctiferalis to three commonly used pesticides

Introduction: Insect cytochrome P450 (CYP450) genes play important roles in the detoxification and metabolism of xenobiotics, such as plant allelochemicals, mycotoxins and pesticides. The polyphagous Conogethes punctiferalis is a serious economic pest of fruit trees and agricultural crops, and it shows high adaptability to different living environments. Methods: The two novel P450 genes CYP6CV1 and CYP6AB51 were identified and characterized. Quantitative real-time PCR (qRT-PCR) technology was used to study the expression patterns of the two target genes in different larval developmental stages and tissues of C. punctiferalis. Furthermore, RNA interference (RNAi) technology was used to study the potential functions of the two P450 genes by treating RNAi-silenced larvae with three commonly used pesticides. Results: The CYP6CV1 and CYP6AB51 genes were expressed throughout various C. punctiferalis larval stages and in different tissues. Their expression levels increased along with larval development, and expression levels of the two target genes in the midgut were significantly higher than in other tissues. The toxicity bioassay results showed that the LC₅₀ values of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin on C. punctiferalis larvae were 0.2028 μg/g, 0.0683 μg/g and 0.6110 mg/L, respectively. After treating with different concentrations of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin (LC₁₀, LC₃₀, LC₅₀), independently, the relative expressions of the two genes CYP6CV1 and CYP6AB51 were significantly induced. After the dsRNA injection, the expression profiles of the two CYP genes were reduced 72.91% and 70.94%, respectively, and the mortality rates of the larvae significantly increased when treated with the three insecticides independently at LC₁₀ values. Discussion: In the summary, after interfering with the CYP6CV1 and CYP6AB51 in C. punctiferalis, respectively, the sensitivity of C. punctiferalis to chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin was significantly increased, indicating that the two CYP6 genes were responsible for the adaptability of C. punctiferalis to the three chemical insecticides in C. punctiferalis. The results from this study demonstrated that CYP6CV1 and CYP6AB51 in C. punctiferalis play crucial roles in the detoxification of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin.

Leveraging RNA Interference to Impact Insecticide Resistance in the Colorado Potato Beetle, Leptinotarsa decemlineata

The Colorado potato beetle, Leptinotarsa decemlineata Say, is a potato pest that can cause important economic losses to the potato industry worldwide. Diverse strategies have been deployed to target this insect such as biological control, crop rotation, and a variety of insecticides. Regarding the latter, this pest has demonstrated impressive abilities to develop resistance against the compounds used to regulate its spread. Substantial work has been conducted to better characterize the molecular signatures underlying this resistance, with the overarching objective of leveraging this information for the development of novel approaches, including RNAi-based techniques, to limit the damage associated with this insect. This review first describes the various strategies utilized to control L. decemlineata and highlights different examples of reported cases of resistances against insecticides for this insect. The molecular leads identified as potential players modulating insecticide resistance as well as the growing interest towards the use of RNAi aimed at these leads as part of novel means to control the impact of L. decemlineata are described subsequently. Finally, select advantages and limitations of RNAi are addressed to better assess the potential of this technology in the broader context of insecticide resistance for pest management.

The effects of short-term glyphosate-based herbicide exposure on insect gene expression profiles

Glyphosate-based herbicides (GBHs) are the most frequently used herbicides worldwide. The use of GBHs is intended to tackle weeds, but GBHs have been shown to affect the life-history traits and antioxidant defense system of invertebrates found in agroecosystems. Thus far, the effects of GBHs on detoxification pathways among invertebrates have not been sufficiently investigated. We performed two different experiments-1) the direct pure glyphosate and GBH treatment, and 2) the indirect GBH experiment via food-to examine the possible effects of environmentally relevant GBH levels on the survival of the Colorado potato beetle (Leptinotarsa decemlineata) and the expression profiles of their detoxification genes. As candidate genes, we selected four cytochrome P450 (CYP), three glutathione-S-transferase (GST), and two acetylcholinesterase (AChE) genes that are known to be related to metabolic or target-site resistances in insects. We showed that environmentally relevant levels of pure glyphosate and GBH increased the probability for higher mortality in the Colorado potato beetle larvae in the direct experiment, but not in the indirect experiment. The GBHs or glyphosate did not affect the expression profiles of the studied CYP, GST, or AChE genes; however, we found a large family-level variation in expression profiles in both the direct and indirect treatment experiments. These results suggest that the genes selected for this study may not be the ones expressed in response to glyphosate or GBHs. It is also possible that the relatively short exposure time did not affect gene expression profiles, or the response may have already occurred at a shorter exposure time. Our results show that glyphosate products may affect the survival of the herbivorous insect already at lower levels, depending on their sensitivity to pesticides.

RNA interference in cytochrome P450 monooxygenase (CYP) gene results in reduced insecticide resistance in Megalurothrips usitatus Bagnall

Genes of the cytochrome P450 (CYP450) superfamily are known to be involved in the evolution of insecticide resistance. In this study, the transcriptomes of two Megalurothrips usitatus Bagnall (Thysanoptera: Thripidae) strains (resistant and susceptible) were screened for detoxification genes. MusiDN2722 encodes a protein composed of 504 amino acid residues with a relative molecular mass of 57.3 kDa. Multiple sequence alignment and phylogenetic analysis showed that MusiDN2722 is a member of the CYP450 family and has characteristics of the conserved CYP6 domain shared by typical CYP450 family members. RT-qPCR (real-time quantitative polymerase chain reaction) analysis showed that MusiDN2722 was upregulated in the acetamiprid-resistant strain compared with the susceptible strain (p < 0.05), and the relative expression level was significantly higher at 48 h after exposure than at 24 h after exposure. The interference efficiency of the injection method was higher than that of the membrane-feeding method. Silencing of MusiDN2722 through RNA interference significantly increased the sensitivity of M. usitatus to acetamiprid. Overall, this study revealed that MusiDN2722 plays a crucial role in the resistance of M. usitatus to acetamiprid. The findings will not only advance our understanding of the role of P450s in insecticide resistance but also provide a potential target for the sustainable control of destructive pests such as thrips.

Spinosad-associated modulation of select cytochrome P450s and glutathione S-transferases in the Colorado potato beetle, Leptinotarsa decemlineata

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is an insect pest that threatens potato crops. Multiple options exist to limit the impact of this pest even though insecticides remain a primary option for its control. Insecticide resistance has been reported in Colorado potato beetles and a better understanding of the molecular players underlying such process is of utmost importance to optimize the tools used to mitigate the impact of this insect. Resistance against the insecticide spinosad has been reported in this insect and this work thus aims at exploring the expression of targets previously associated with insecticide response in Colorado potato beetles exposed to this compound. Amplification and quantification of transcripts coding for cytochrome P450s and glutathione S-transferases were conducted via qRT-PCR in insects treated with varying doses of spinosad and for different time duration. This approach notably revealed differential expression of CYP6a23 and CYP12a5 in insects exposed to low doses of spinosad for 4 h as well as modulation of CYP6a13, CYP6d4, GST, GST1, and GST1-Like in insects treated with high doses of spinosad for the same duration. RNAi-based targeting of CYP4g15 and CYP6a23 was associated with marked reduction of transcript expression 7 days following dsRNA injection and reduction of the former had a marked impact on insect viability. In general, results presented here provide novel information regarding the expression of transcripts relevant to spinosad response in Colorado potato beetles and reveal a novel target to consider in the development of RNAi-based strategies aimed at this potato pest.

The role of structural variants in pest adaptation and genome evolution of the Colorado potato beetle, Leptinotarsa decemlineata (Say)

Structural variation has been associated with genetic diversity and adaptation. Despite these observations, it is not clear what their relative importance is for evolution, especially in rapidly adapting species. Here, we examine the significance of structural polymorphisms in pesticide resistance evolution of the agricultural super-pest, the Colorado potato beetle, Leptinotarsa decemlineata. By employing a parent offspring trio sequencing procedure, we develop highly contiguous reference genomes to characterize structural variation. These updated assemblies represent >100-fold improvement of contiguity and include derived pest and ancestral nonpest individuals. We identify >200,000 structural variations, which appear to be nonrandomly distributed across the genome as they co-occur with transposable elements and genes. Structural variations intersect with exons in a large proportion of gene annotations (~20%) that are associated with insecticide resistance (including cytochrome P450s), development, and transcription. To understand the role structural variations play in adaptation, we measure their allele frequencies among an additional 57 individuals using whole genome resequencing data, which represents pest and nonpest populations of North America. Incorporating multiple independent tests to detect the signature of natural selection using SNP data, we identify 14 genes that are probably under positive selection, include structural variations, and SNPs of elevated frequency within the pest lineages. Among these, three are associated with insecticide resistance based on previous research. One of these genes, CYP4g15, is coinduced during insecticide exposure with glycosyltransferase-13, which is a duplicated gene enclosed within a structural variant adjacent to the CYP4g15 genic region. These results demonstrate the significance of structural variations as a genomic feature to describe species history, genetic diversity, and adaptation.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Decreased cuticular penetration minimizes the impact of the pyrethroid insecticide λ-cyhalothrin on the insect predator Eocanthecona furcellata

The use of broad-spectrum pesticides may reduce the biological control efficacy of predatory arthropods. Hence, the risks of pesticides to predators need to be evaluated. Here, we assessed the effects of a broad spectrum pyrethroid λ-cyhalothrin on a polyphagous predatory insect Eocanthecona furcellata via contact exposure route. The recommended application rate of λ-cyhalothrin was lower than the LR₅₀ and HQ _(in-field) was equal to 0.57, indicating the risk of λ-cyhalothrin to E. furcellata was low. Dried λ-cyhalothrin residue had no effect on the mortality, body weight, protein content of cuticle, or activities of major detoxification enzymes in E. furcellata. Residual of λ-cyhalothrin was only detected in the cuticle and legs of E. furcellata with a decreasing trend as time went by and no λ-cyhalothrin was detected inside the body. Additionally, a comparative transcriptome analysis was conducted to study global changes in gene expression in E. furcellata at different time points following exposure to λ-cyhalothrin-contaminated environment. A total of 57,839 unigenes with an average length of 1044 bp and an N50 of 1820 bp were obtained. In total, 118 and 109 differentially expressed genes (DEGs) at 12 h, and 60 h were identified between two groups. The DEGs were largely enriched in functional categories related to the structural constituent of cuticle. Accordingly, multiple cuticle protein-coding genes were up-regulated at 12 h after pesticide exposure. The present study stressed the importance of evaluating the compatibility between a specific pesticide (λ-cyhalothrin) and E. furcellata via simulating the releasing predators after insecticide application. The data could help optimize the pesticide use, optimizing the ecological services of E. furcellata as a BCA, and expanding its use into more areas of agriculture.

Gut transcriptome analysis of P450 genes and cytochrome P450 reductase in three moth species feeding on gymnosperms or angiosperms

Cytochrome P450 enzymes (P450s, CYPs) are a superfamily of heme–thiolate proteins involved in the metabolism of endogenous and exogenous substances in insects. In this study, the identification of putative P450 proteins was done and the elimination of the repeated sequences resulted in 57 proteins from Gastropacha populifolia, 63 proteins from Dendrolimus punctatus, and 53 proteins from Dendrolimus tabulaeformis. The putative P450 proteins were aligned together with seven other insect species based on five conserved domains. A total of ten co-orthologous groups were identified. Interestingly, one co-orthologous gene, CYP4g15 in CYP4 clan, was identified and its 3D structure analysis showed that the highly conserved sites of the predicted motifs were close to the active sites of P450. Furthermore, this study revealed that insect CYP4g15 and two bacteria cytochrome P450 were monophyletic. This suggests that insects CYP4g15 are not only functionally conserved but also an ancient gene originating from different bacteria species.

The insecticidal effect of the botanical insecticide chlorogenic acid on Mythimna separata (Walker) is related to changes in MsCYP450 gene expression

The oriental armyworm Mythimna separata (Walker) (Lepidoptera: Noctuidae) can feed on the leaves of many crops, resulting in vast areas of damage and severe losses. Therefore, this insect has become a significant agricultural pest in north Asia. In this study, we fed 3^rd instar larvae with artificial diets containing different concentrations of chlorogenic acid and found a significant lethal effect and the mortality increased with increasing chlorogenic acid concentration. Next, we measured the sublethal effect of chlorogenic acid at LC₂₀ on the growth and development of M. separata larvae. The durations of the 4^th and 5^th instar were longer than those of the control group (prolonged by 0.8 and 0.6 days, respectively), and the 6^th instar was shorter (by 1.1 days). The total survival rate, pupation rate, eclosion rate, sex ratio, and oviposition amount in the LC₂₀ chlorogenic acid-treated group were significantly lower than those in the control group. Furthermore, transcriptome analysis of 3^rd instar larvae fed various concentrations of chlorogenic acid revealed that several MsCYP450 genes were significantly up-regulated, and this finding was further validated by qRT-PCR. In addition, various concentrations of chlorogenic acid and different treatment times significantly affected the enzyme activity of CYP450 in 3^rd instar larvae. Importantly, dietary ingestion of dsMsCYP450 significantly reduced the mRNA level of MsCYP450 genes and increased mortality in the presence of chlorogenic acid. Our results revealed that MsCYP6B6, MsCYP321A7, and MsCYP6B7-like play an essential role in the detoxification of chlorogenic acid by M. separata. This study provides evidence of control effect by botanical insecticide chlorogenic acid on M. separata, and potential detoxification mechanism mediated by P450 of botanical insecticide in arthropods.

Acute toxic effects of diclofenac exposure on freshwater crayfish (Procambarus clarkii): Insights from hepatopancreatic pathology, molecular regulation and intestinal microbiota

In this study, we exposed adult male crayfish (Procambarus clarkii) to different concentrations of diclofenac (DCF) for 96 h. In the meantime, we investigated the alternations of hepatopancreatic pathology, molecular regulation and intestinal microbiota of P. clarkii exposed to DCF. The results demonstrated DCF led to histological changes including epithelium vacuolization and tubule lumen dilatation in the hepatopancreas. Transcriptome sequencing analysis showed that 642 and 586 genes were differentially expressed in the hepatopancreas of P. clarkii exposed to 1 and 10 mg/L DCF, respectively. DCF could affect the functions of antioxidation, immunity and metabolism of hepatopancreas by inducing the abnormal expressions of immune- and redox-related genes. GO enrichment results demonstrated that 10 mg/L DCF exposure could modulate the processes of molting, amino sugar metabolism, protein hydrolysis and intracellular protein translocation of P. clarkii. Additionally, the abundances of bacterial families including Shewanellaceae, Bacteroidaceae, Vibrionaceae, Erysipelotrichaceae, Aeromonadaceae, Moraxellaceae, etc. in the intestine were significantly changed after DCF exposure, and the disruption of intestinal flora might further cause abnormal intestinal metabolism in P. clarkii. This study provides novel mechanistic insights into the toxic effects of anti-inflammatory drugs on aquatic crustaceans.

Insights into the Use of Eco-Friendly Synergists in Resistance Management of Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle (CPB), Leptinotarsa decemlineata (Say), is the most notorious insect pest of potato globally. Injudicious use of insecticides for management of this pest has resulted in resistance to all major groups of insecticides along with many human, animal health, and environmental concerns. Additionally, the input cost of insecticide development/discovery is markedly increasing because each year thousands of chemicals are produced and tested for their insecticidal properties, requiring billions of dollars. For the management of resistance in insect pests, synergists can play a pivotal role by reducing the application dose of most insecticides. These eco-friendly synergists can be classified into two types: plant-based synergists and RNAi-based synergists. The use of plant-based and RNAi-based synergists in resistance management of insect pests can give promising results with lesser environmental side effects. This review summarizes the resistance status of CPB and discusses the potential advantage of plant-based and RNAi-based synergists for CPB resistance management. It will motivate researchers to further investigate the techniques of using plant- and RNAi-based synergists in combination with insecticides.

Genome-wide identification of long non-coding (lncRNA) in Nilaparvata lugens's adaptability to resistant rice

Background

The brown planthopper (BPH), Nilaparvata lugens (Stål), is a very destructive pest that poses a major threat to rice plants worldwide. BPH and rice have developed complex feeding and defense strategies in the long-term co-evolution.

Methods

To explore the molecular mechanism of BPH's adaptation to resistant rice varieties, the lncRNA expression profiles of two virulent BPH populations were analyzed. The RNA-seq method was used to obtain the lncRNA expression data in TN1 and YHY15.

Results

In total, 3,112 highly reliable lncRNAs in TN1 and YHY15 were identified. Compared to the expression profiles between TN1 and YHY15, 157 differentially expressed lncRNAs, and 675 differentially expressed mRNAs were identified. Further analysis of the possible regulation relationships between differentially expressed lncRNAs and differentially expressed mRNAs, identified three pair antisense targets, nine pair cis-regulation targets, and 3,972 pair co-expressed targets. Function enriched found arginine and proline metabolism, glutathione metabolism, and carbon metabolism categories may significantly affect the adaptability in BPH when it is exposed to susceptible and resistant rice varieties. Altogether, it provided scientific data for the study of lncRNA regulation of brown planthopper resistance to rice. These results are helpful in the development of new control strategies for host defense against BPH and breeding rice for high yield.

Characterization of Insecticide Response-Associated Transcripts in the Colorado Potato Beetle: Relevance of Selected Cytochrome P450s and Clothianidin

Simple Summary

The Colorado potato beetle is an insect pest that can significantly harm potato crops. Various approaches are available to mitigate its damages including the use of insecticides. Unfortunately, its ability to develop resistance towards these compounds is substantial, and understanding the basis of this process is of utmost importance to design strategies to limit the impact of this insect. This work thus aims at quantifying the expression of key transcripts coding for proteins associated with insecticide resistance in Colorado potato beetles exposed to four insecticides. Significant variations were observed, notably in insects exposed to the insecticide clothianidin. Interestingly, subsequent reduction of endogenous levels of selected targets modulated by clothianidin was associated with increased insect susceptibility to this neonicotinoid. These results further highlight molecular players with potential relevance for insecticide resistance, and introduce novel targets that underlie clothianidin resistance in the Colorado potato beetle.

Abstract

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is known for its capacity to cause significant damages to potato crops worldwide. Multiple approaches have been considered to limit its spread including the use of a diverse arsenal of insecticides. Unfortunately, this insect frequently develops resistance towards these compounds. Investigating the molecular bases underlying the response of L. decemlineata against insecticides is of strong interest to ultimately devise novel and targeted approaches aimed at this pest. This work aimed to characterize, via qRT-PCR, the expression status of targets with relevance to insecticide response, including ones coding for cytochrome P450s, glutathione s-transferases, and cuticular proteins, in L. decemlineata exposed to four insecticides; chlorantraniliprole, clothianidin, imidacloprid, and spinosad. Modulation of levels associated with transcripts coding for selected cytochrome P450s was reported in insects treated with three of the four insecticides studied. Clothianidin treatment yielded the most variations in transcript levels, leading to significant changes in transcripts coding for CYP4c1, CYP4g15, CYP6a13, CYP9e2, GST, and GST-1-Like. Injection of dsRNA targeting CYP4c1 and CYP9e2 was associated with a substantial decrease in expression levels and was, in the case of the latter target, linked to a greater susceptibility of L. decemlineata towards this neonicotinoid, supporting a potential role for this target in clothianidin response. Overall, this data further highlights the differential expression of transcripts with potential relevance in insecticide response, as well as generating specific targets that warrant investigation as novel dsRNA-based approaches are developed against this insect pest.

The Leptinotarsa forkhead transcription factor O exerts a key function during larval-pupal-adult transition

Forkhead box O (FoxO) protein, a major downstream transcription factor of insulin/insulin-like growth factor signaling/target of rapamycin pathway (IIS/TOR), is involved in the regulation of larval growth and the determination of organ size. FoxO also interacts with 20-hydroxyecdysone (20E) and juvenile hormone (JH) signal transduction pathways, and hence is critical for larval development in holometabolans. However, whether FoxO plays a critical role during larval metamorphosis needs to be further determined in Leptinotarsa decemlineata. We found that 20E stimulated the expression of LdFoxO. RNA interference (RNAi)-aided knockdown of LdFoxO at the third-instar stage repressed 20E signaling and reduced larval weight. Although the resultant larvae survived through the third-fourth instar ecdysis, around 70% of the LdFoxO depleted moribund beetles developmentally arrested at prepupae stage. These LdFoxO depleted beetles were completely wrapped in the larval exuviae, gradually darkened and finally died. Moreover, approximately 12% of the LdFoxO RNAi beetles died as pharate adults. Ingestion of either 20E or JH by the LdFoxO depletion beetles excessively rescued the corresponding hormonal signals, but could not alleviate larval performance and restore defective phenotypes. Therefore, FoxO plays an important role in regulation of larval-pupal-adult transformation in L. decemlineata, in addition to mediation of IIS/TOR pathway and stimulation of ecdysteroidogenesis.

Full length sequencing reveals novel transcripts of detoxification genes along with related alternative splicing events and lncRNAs in Phyllotreta striolata

The striped flea beetle, Phyllotreta striolata (Fabricius), damages crops in the Brassicaceae. The genetic data for this pest are insufficient to reveal its insecticide resistance mechanisms or to develop molecular markers for resistance monitoring. We used PacBio Iso-Seq technology to sequence the full-length transcriptome of P. striolata. After isoform sequence clustering and removal of redundant transcripts, a total of 41,293 transcripts were obtained, and 35,640 of these were annotated in the database of gene products. Structure analysis uncovered 4,307 alternative splicing events, and 3,836 sequences were recognized as lncRNAs. Transcripts with the complete coding region of important detoxification enzymes were further classified. There were 57 transcripts of P450s distributed in CYP2, CYP3, CYP4, and Mito CYP clades, 29 transcripts of ESTs from 4 functional groups, 17 transcripts of GSTs classified into 5 families, 51 transcripts of ABCs distributed in 6 families, and 19 transcripts of UGTs. Twenty-five lncRNAs were predicted to be regulators of these detoxification genes. Full-length transcriptome sequencing is an efficient method for molecular study of P. striolata and it is also useful for gene function analysis.

Transcriptomics-Based Approach Identifies Spinosad-Associated Targets in the Colorado Potato Beetle, Leptinotarsa decemlineata

Simple Summary

The Colorado potato beetle Leptinotarsa decemlineata is a potato pest that can cause substantial damages to potato crops worldwide. Multiple approaches have been leveraged to control this pest including the use of a variety of insecticides. Resistance to different insecticides aimed at controlling this insect has been reported and much work has been conducted in recent years to elucidate the underlying molecular changes associated with insecticide resistance in L. decemlineata. However, information is sparse regarding the molecular impact associated with spinosad treatment in this insect pest. The current study thus explores transcriptional changes associated with spinosad response in L. decemlineata exposed to this compound using high-throughput sequencing. Results presented show multiple transcripts of interest that exhibit differential expression in spinosad-treated L. decemlineata and provide a preliminary footprint of transcripts affected by this insecticide in this potato pest. Select targets identified in this signature should be further explored in follow-up studies to better characterize their contribution, if any, in the process of spinosad resistance.

Abstract

The Colorado potato beetle Leptinotarsa decemlineata is an insect pest that threatens potato crops globally. The primary method to control its damage on potato plants is the use of insecticides, including imidacloprid, chlorantraniliprole and spinosad. However, insecticide resistance has been frequently observed in Colorado potato beetles. The molecular targets and the basis of resistance to imidacloprid and chlorantraniliprole have both been previously quantified. This work was undertaken with the overarching goal of better characterizing the molecular changes associated with spinosad exposure in this insect pest. Next-generation sequencing was conducted to identify transcripts that were differentially expressed between Colorado potato beetles exposed to spinosad versus control insects. Results showed several transcripts that exhibit different expression levels between the two conditions, including ones coding for venom carboxylesterase-6, chitinase 10, juvenile hormone esterase and multidrug resistance-associated protein 4. In addition, several microRNAs, such as miR-12-3p and miR-750-3p, were also modulated in the investigated conditions. Overall, this work reveals a molecular footprint underlying spinosad response in Colorado potato beetles and provides novel leads that could be targeted as part of RNAi-based approaches to control this insect pest.

Ecdysone receptor isoforms play distinct roles in larval-pupal-adult transition in Leptinotarsa decemlineata

A heterodimer of two nuclear receptors, ecdysone receptor (EcR) and ultraspiracle, mediates 20-hydroxyecdysone (20E) signaling to modulate many aspects in insect life, such as molting and metamorphosis, reproduction, diapause and innate immunity. In the present paper, we intended to determine the isoform-specific roles of EcR during larval-pupal-adult transition in the Colorado potato beetle. Double-stranded RNAs (dsRNAs) were prepared using the common (dsEcR) or isoform-specific (dsEcRA, dsEcRB1) regions of EcR as templates. Ingestion of either dsEcR or dsEcRA, rather than dsEcRB1, by the penultimate (3rd) and final (4th) instar larvae caused failure of larval-pupal and pupal-adult ecdysis. The RNA interference (RNAi) larvae remained as prepupae, or became deformed pupae and adults. Determination of messenger RNA (mRNA) levels of EcR isoforms found that LdEcRA regulates the expression of LdEcRB1. Moreover, silencing the two EcR transcripts, LdEcRA or LdEcRB1 reduced the mRNA levels of Ldspo and Ldsad, and lowered 20E titer. In contrast, the expression levels of HR3, HR4, E74 and E75 were significantly decreased in the LdEcR or LdEcRA RNAi larvae, but not in LdEcRB1 depleted specimens. Dietary supplement with 20E did not restore the expression of five 20E signaling genes (USP, HR3, HR4, E74 and E75), and only partially alleviated the pupation defects in dsEcR- or dsEcRA-fed beetles. These data suggest that EcR plays isoform-specific roles in the regulation of ecdysteroidogenesis and the transduction of 20E signal in L. decemlineata.

Microcystin-LR-induced changes of hepatopancreatic transcriptome, intestinal microbiota, and histopathology of freshwater crayfish (Procambarus clarkii)

As a hepatotoxin, microcystin-LR (MC-LR) poses a great threat to aquatic organisms. In this research, the hepatopancreatic transcriptome, intestinal microbiota, and histopathology of Procambarus clarkii (P. clarkii) in response to acute MC-LR exposure were studied. RNA-seq analysis of hepatopancreas identified 372 and 781 differentially expressed genes (DEGs) after treatment with 10 and 40 μg/L MC-LR, respectively. Among the DEGs, 23 genes were immune-related and 21 genes were redox-related. GO functional enrichment analysis revealed that MC-LR could impact nuclear-transcribed mRNA catabolic process, cobalamin- and heme-related processes, and sirohydrochlorin cobaltochelatase activity of P. clarkii. In addition, the only significantly enriched KEGG pathway induced by MC-LR was galactose metabolism pathway. Meanwhile, sequencing of the bacterial 16S rRNA gene demonstrated that MC-LR decreased bacterial richness and diversity, and altered the intestinal microbiota composition. At the phylum level, after 96 h, the abundance of Verrucomicrobia decreased after treatment with 10 and 40 μg/L MC-LR, while Firmicutes increased in the 40 μg/L MC-LR-treated group. At the genus level, the abundances of 15 genera were significantly altered after exposure to MC-LR. Our research demonstrated that MC-LR exposure caused histological alterations such as structural damage of hepatopancreas and intestines. This research provides an insight into the mechanisms associated with MC-LR toxicity in aquatic crustaceans.

Functional characterization of ultraspiracle in Leptinotarsa decemlineata using RNA interference assay

A heterodimer of ultraspiracle (USP) and ecdysone receptor (EcR) mediates 20-hydroxyecdysone (20E) signalling cascade to regulate insect moulting and metamorphosis. However, at least two questions remain to be addressed in terms of the molecular importance of USP in insect species. First, is USP involved in both regulation of ecdysteroidogenesis and mediation of 20E signalling in non-drosophilid insects, as in Drosophila melanogaster? Second, does USP play any role in larval metamorphosis except as the partner of heterodimeric receptor to activate the downstream 20E signalling genes? In this paper, we found that RNA interference (RNAi) of LdUSP in the final (fourth) instar larvae reduced the messenger RNA levels of four ecdysteroidogenesis genes (Ldspo, Ldphm, Lddib and Ldsad) and 20E titre, and repressed the expression of five 20E signal genes (EcRA, HR3, HR4, E74 and E75) in Leptinotarsa decemlineata. The LdUSP RNAi larvae remained as prepupae, with developing antennae, legs and discs of forewings and hindwings. Dietary supplement with 20E restored the expression of the five 20E signal genes, but only partially alleviated the decreased pupation rate in LdUSP RNAi beetles. Knockdown of LdUSP at the penultimate (third) instar larvae did not affect third-fourth instar moulting. However, silencing LdUSP caused similar but less severe impairments on pupation. Accordingly, we propose that USP is undoubtedly necessary for ecdysteroidogenesis, for mediation of 20E signalling and for initiation of metamorphosis in L. decemlineata.

Silencing Taiman impairs larval development in Leptinotarsa decemlineata

An exploration of novel control strategies for Leptinotarsa decemlineata is becoming more pressing given rapid evolution of insecticide resistance and rise of production loss of potato. Dietary delivery of bacterially expressed double-stranded RNA (dsRNA) is a promising alternative for management. An important first step is to uncover possible RNA-interference (RNAi)-target genes effective against both young and old larvae. Taiman (Tai) is a basic-helix-loop-helix/Per-Arnt-Sim transcription factor that is involved in the mediation of both juvenile hormone (JH) and 20-hydroxyecdysone (20E) signaling. In the present paper, we found that continuous ingestion of dsTai for three days by third (penultimate)-instar larvae caused approximately 20% larval mortality and 80% pupation failure. The larval lethality resulted from failed cuticle and tracheae shedding, which subsequently reduced foliage consumption and nutrient absorption, and depleted lipid stores. In contrast, pupation failure derived from disturbed JH and 20E signals, and disordered nutrient homeostasis including, among others, inhibition of trehalose metabolism and reduction of chitin content. Knockdown of LdTai caused similar larval lethality and pupation impairment in second and fourth (final) larval instars. Therefore, LdTai is among the most attractive candidate genes for RNAi to control L. decemlineata larvae.

Importance of Taiman in Larval-Pupal Transition in Leptinotarsa decemlineata

Insect Taiman (Tai) binds to methoprene-tolerant to form a heterodimeric complex, mediating juvenile hormone (JH) signaling to regulate larval development and to prevent premature metamorphosis. Tai also acts as a steroid receptor coactivator of 20-hydroxyecdysone (20E) receptor heterodimer, ecdysone receptor (EcR) and Ultraspiracle (USP), to control the differentiation of early germline cells and the migration of specific follicle cells and border cells in ovaries in several insect species. In holometabolous insects, however, whether Tai functions as the coactivator of EcR/USP to transduce 20E message during larval-pupal transition is unknown. In the present paper, we found that the LdTai mRNA levels were positively correlated with circulating JH and 20E titers in Leptinotarsa decemlineata; and ingestion of either JH or 20E stimulated the transcription of LdTai. Moreover, RNA interference (RNAi)-aided knockdown of LdTai at the fourth (final) instar stage repressed both JH and 20E signals, inhibited larval growth and shortened larval developing period. The knockdown caused 100% larval lethality due to failure of larval-pupal ecdysis. Under the apolysed larval cuticle, the LdTai RNAi prepupae possessed pupal thorax. In contrast, the process of tracheal ecdysis was uncompleted. Neither JH nor 20E rescued the aforementioned defectives in LdTai RNAi larvae. It appears that Tai mediates both JH and 20E signaling. Our results uncover a link between JH and 20E pathways during metamorphosis in L. decemlineata.

Global transcriptome profiling and functional analysis reveal that tissue-specific constitutive overexpression of cytochrome P450s confers tolerance to imidacloprid in palm weevils in date palm fields

Background

Cytochrome P450-dependent monooxygenases (P450s), constituting one of the largest and oldest gene superfamilies found in many organisms from bacteria to humans, play a vital role in the detoxification and inactivation of endogenous toxic compounds. The use of various insecticides has increased over the last two decades, and insects have developed resistance to most of these compounds through the detoxifying function of P450s. In this study, we focused on the red palm weevil (RPW), Rhynchophorus ferrugineus, the most devastating pest of palm trees worldwide, and demonstrated through functional analysis that upregulation of P450 gene expression has evolved as an adaptation to insecticide stress arising from exposure to the neonicotinoid-class systematic insecticide imidacloprid.

Results

Based on the RPW global transcriptome analysis, we identified 101 putative P450 genes, including 77 likely encoding protein coding genes with ubiquitous expression. A phylogenetic analysis revealed extensive functional and species-specific diversification of RPW P450s, indicating that multiple CYPs actively participated in the detoxification process. We identified highly conserved paralogs of insect P450s that likely play a role in the development of resistance to imidacloprid: Drosophila Cyp6g1 (CYP6345J1) and Bemisia tabaci CYP4C64 (CYP4LE1). We performed a toxicity bioassay and evaluated the induction of P450s, followed by the identification of overexpressed P450s, including CYP9Z82, CYP6fra5, CYP6NR1, CYP6345J1 and CYP4BD4, which confer cross-resistance to imidacloprid. In addition, under imidacloprid insecticide stress in a date palm field, we observed increased expression of various P450 genes, with CYP9Z82, CYP4BD4, CYP6NR1 and CYP6345J1 being the most upregulated detoxification genes in RPWs. Expression profiling and cluster analysis revealed P450 genes with multiple patterns of induction and differential expression. Furthermore, we used RNA interference to knock down the overexpressed P450s, after which a toxicity bioassay and quantitative expression analysis revealed likely candidates involved in metabolic resistance against imidacloprid in RPW. Ingestion of double-stranded RNA (dsRNA) successfully knocked down the expression of CYP9Z82, CYP6NR1 and CYP345J1 and demonstrated that silencing of CYP345J1 and CYP6NR1 significantly decreased the survival rate of adult RPWs treated with imidacloprid, indicating that overexpression of these two P450s may play an important role in developing tolerance to imidacloprid in a date palm field.

Conclusion

Our study provides useful background information on imidacloprid-specific induction and overexpression of P450s, which may enable the development of diagnostic tools/markers for monitoring the spread of insecticide resistant RPWs. The observed trend of increasing tolerance to imidacloprid in the date palm field therefore indicated that strategies for resistance management are urgently needed.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5837-4) contains supplementary material, which is available to authorized users.

Identification of Cuticular Protein Genes in the Colorado Potato Beetle Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Structural cuticular proteins (CPs) are the primary components of insect cuticle, linings of salivary gland, foregut, hindgut and tracheae, and midgut peritrophic membrane. Variation of CPs in insect cuticle can cause penetration resistance to insecticides. Moreover, depletion of specific CP by RNA interference may be a suitable way for the development of potential pest control traits. Leptinotarsa decemlineata (Say) CPs are poorly characterized at present, and therefore, we mined the genome and transcriptome data to better annotate and classify L. decemlineata CPs in this study, by comparison with the annotated CPs of Tribolium castaneum Browse (Coleoptera: Tenebrionidae). We identified 175 CP genes. Except one miscellaneous CP with an 18-amino acid motif, these CPs were classified into 7 families based on motifs and phylogenetic analyses (CPs with a Rebers and Riddiford motif, CPR; CPs analogous to peritrophins, CPAP3 and CPAP1; CPs with a tweedle motif, TWDL; CPs with a 44-amino acid motif, CPF; CPs that are CPF-like, CPFL; and CPs with two to three copies of C-X5-C motif, CPCFC). Leptinotarsa decemlineata CPRs could be categorized into three subfamilies: RR-1 (50), RR-2 (85), and RR-3 (2). The RR-1 proteins had an additional motif with a conserved YTADENGF sequence. The RR-2 members possessed a conserved RDGDVVKG region and three copes of G-x(3)-VV. Few genes were found in TWDL (9), CPAP1 (9), CPAP3 (8), CPF (5), CPFL (4), and CPCFC (2) families. The findings provide valuable information to explore molecular modes of penetration resistance to insecticides and to develop dsRNA-based control method in L. decemlineata.

Silencing chitin deacetylase 2 impairs larval-pupal and pupal-adult molts in Leptinotarsa decemlineata

Insect chitin deacetylases (CDAs) are carbohydrate esterases that catalyze N-deacetylation of chitin to generate chitosan, a process essential for chitin organization and compactness during the formation of extracellular chitinous structure. Here we identified two CDA2 splice variants (LdCDA2a and LdCDA2b) in Leptinotarsa decemlineata. Both splices were abundantly expressed in larval foregut, rectum, and epidermis; their levels peaked immediately before ecdysis within each instar. In vivo results revealed that the two isoforms transcriptionally responded, positively and negatively respectively, to 20-hydroxyecdysone and juvenile hormone signaling pathways. RNA interference (RNAi)-aided knockdown of the two LdCDA2 variants (hereafter LdCDA2) or LdCDA2b, rather than LdCDA2a, resulted in three negative effects. First, foliage consumption was significantly reduced, larval developing period was lengthened, and larval growth was retarded. Second, chitin contents were reduced, whereas glucose, trehalose, and glycogen contents were increased in the LdCDA2 and LdCDA2b RNAi larvae. Third, approximately 20% of LdCDA2 and LdCDA2b RNAi larvae were trapped within the exuviae and finally died. About 60% of the abnormal pupae died as pharate adults. Around 20% of the RNAi pupae emerged as deformed adults, with small size and wrinkled wings. These adults eventually died within 1 week after molting. Our results reveal that knockdown of CDA2 affects chitin accumulation. Consequently, LdCDA2 may be a potential target for control of L. decemlineata larvae.

Identification of cytochrome P450 monooxygenase genes and their expression in response to high temperature in the alligatorweed flea beetle Agasicles hygrophila (Coleoptera: Chrysomelidae)

Cytochrome P450 monooxygenases (P450s) are a large class of enzymes that play essential roles in metabolic processes such as hormone synthesis and the catabolism of toxins and other chemicals in insects. In the present study, we identified 82 P450 genes using comprehensive RNA sequencing in the flea beetle Agasicles hygrophila, and all of the sequences were validated by cloning and sequencing. Phylogenetic analysis showed that the P450 genes in A. hygrophila fell into the mitochondrial clan, CYP2 clan, CYP3 clan and CYP4 clan and were classified into 20 families and 48 subfamilies. Most A. hygrophila P450 genes had high sequence homology with those from other coleopteran insects. To understand the effects of high temperatures on the metabolic processes of female and male adults, we studied the effects of two temperature regimes (constant temperature of 28 °C for 20 h with a 4-h period of high temperatures of 30 °C and 39 °C) on the expression levels of P450 genes in A. hygrophila using RT-PCR and qRT-PCR. The results showed that there were no differences in expression in 30 P450 genes between the control and high-temperature-treated A. hygrophila adults, while 22 P450 genes showed up-regulated expression and 19 P450 genes were down-regulated in A. hygrophila female adults after high-temperature treatment. For A. hygrophila male adults exposed to high temperatures, we found that 8 P450 genes had higher expression levels and 12 P450 genes had lower expression levels under the same conditions. The P450 genes are candidates that showed significantly different expression levels after high-temperature treatments in A. hygrophila adults, and further studies are needed to determine their possible roles in metabolic processes during the response to elevated temperatures.

RNA interference against the putative insulin receptor substrate gene chico affects metamorphosis in Leptinotarsa decemlineata

It is noted that insect insulin/insulin-like growth factor/target of rapamycin signaling is critical for the regulation of metamorphosis in holometabolous insects. However, the molecular mechanism remains undetermined. Our previous findings reveal that RNA interference (RNAi)-mediated knockdown of an insulin gene (LdILP2) in Leptinotarsa decemlineata disturbs both 20-hydroxyecdysone (20E) and juvenile hormone (JH) signaling, and impairs pupation. In the present paper, we further observed that the expression of the insulin receptor substrate gene chico (Ldchico) and the phosphoinositide-3-kinase gene pi3k (Ldpi3k92E) was repressed in LdILP2 depleted larvae. Moreover, RNAi of Ldchico or Ldpi3k92E decreased food consumption, affected absorption and metabolism of amino acids and sugars, and reduced expression of several 20E (LdEcR, LdHR3 and LdE75) and JH (LdJHAMT, LdKr-h1 and LdHairy) signaling genes. As a result, larval development was postponed and larval growth was inhibited. Intriguingly, knockdown of Ldchico, rather than Ldpi3k92E, impaired larval-pupal and pupal-adult ecdysis, and specifically repressed transcription of another 20E signaling gene LdUSP. Ingestion of 20E rescued the expression of LdEcR, LdHR3 and LdE75, whereas 20E feeding restored neither the decreased LdUSP mRNA level, nor the reduced pupation and adult emergence rates in Ldchico RNAi larvae. Therefore, Chico is critical for the regulation of larval-pupal-adult transition by a PI3K-independent pathway, perhaps through activation of USP in L. decemlineata.

Cap 'n' collar C regulates genes responsible for imidacloprid resistance in the Colorado potato beetle, Leptinotarsa decemlineata

The Colorado potato beetle (CPB), Leptinotarsa decemlineata developed resistance to imidacloprid after exposure to this insecticide for multiple generations. Our previous studies showed that xenobiotic transcription factor, cap 'n' collar isoform C (CncC) regulates the expression of multiple cytochrome P450 genes, which play essential roles in resistance to plant allelochemicals and insecticides. In this study, we sought to obtain a comprehensive picture of the genes regulated by CncC in imidacloprid-resistant CPB. We performed sequencing of RNA isolated from imidacloprid-resistant CPB treated with dsRNA targeting CncC or gene coding for green fluorescent protein (control). Comparative transcriptome analysis showed that CncC regulated the expression of 1798 genes, out of which 1499 genes were downregulated in CncC knockdown beetles. Interestingly, expression of 79% of imidacloprid induced P450 genes requires CncC. We performed quantitative real-time PCR to verify the reduction in the expression of 20 genes including those coding for detoxification enzymes (P450s, glutathione S-transferases, and esterases) and ABC transporters. The genes coding for ABC transporters are induced in insecticide resistant CPB and require CncC for their expression. Knockdown of genes coding for ABC transporters simultaneously or individually caused an increase in imidacloprid-induced mortality in resistant beetles confirming their contribution to insecticide resistance. These studies identified CncC as a transcription factor involved in regulation of genes responsible for imidacloprid resistance. Small molecule inhibitors of CncC or suppression of CncC by RNAi could provide effective synergists for pest control or management of insecticide resistance.

Requirement of Leptinotarsa decemlineata gene within the 74EF puff for larval-pupal metamorphosis and appendage growth

Two Drosophila melanogaster E-twenty-six domain transcription factor isoforms (E74A and E74B) act differentially at the start of the 20-hydroxyecdysone (20E) signalling cascade to regulate larval-pupal metamorphosis. In the present paper, we identified the two isoforms (LdE74A and LdE74B) in Leptinotarsa decemlineata. During the larval development stage, the mRNA transcript levels of the two LdE74 isoforms were correlated with circulating 20E titres. In vitro midgut culture and in vivo dietary supplementation with 20E revealed that the presence of 20E induced expression peaks of both LdE74A and LdE74B, with similar patterns observed for the two isoforms. Moreover, the mRNA transcript levels of both LdE74A and LdE74B isoforms were significantly downregulated in the L. decemlineata ecdysone receptor RNA interference (RNAi) specimens, but not in the LdE75 RNAi beetles. Ingestion of 20E reduced the larval fresh weights and shortened the larval development period, irrespective of knockdown of LdE74 or not. RNAi of LdE74 did not affect 20E-induced expression of the Ecdysone induced protein 75-hormone receptor 3-fushi tarazu factor 1 (E75-HR3-FTZ-F1) transcriptional cascade. Thus, it seems that LdE74 mediates 20E signalling independent of the E75-HR3-FTZ-F1 transcriptional cascade. Furthermore, silencing of both LdE74 isoforms caused failure of ecdysis. Most of the LdE74 RNAi beetles remained as prepupae. The LdE74 RNAi prepupae exhibited adult character-like forms underneath after removal of the apolysed larval cuticle. Their appendages such as antennae, legs and wings were shorter than those of control larvae. Only a few LdE74 RNAi larvae finally became deformed pupae, with shortened antennae and legs. Therefore, LdE74 is required for larval-pupal metamorphosis and appendage growth in L. decemlineata.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Identification of two novel P450 genes and their responses to deltamethrin in the cabbage moth, Mamestra brassicae Linnaeus

Cytochrome P450 monooxygenases, found in virtually all living organisms, play an important role in the metabolism of xenobiotics such as drugs, pesticides, and plant toxins. In this study, we identified two novel cytochrome P450 genes from the cabbage moth, Mamestra brassicae Linnaeus. They were named CYP4M51 and CYP6AB56 (GenBank Accession Nos.: KX008607 and JQ901385, respectively) by the P450 Nomenclature Committee. Real-time quantitative reverse transcription revealed that CYP4M51 and CYP6AB56 were highly expressed in the fat bodies and were differentially expressed at different larval developmental stages. Expression levels of these two cytochrome P450 genes were up-regulated by deltamethrin. Analyses of their detoxification roles using RNA interference followed by a deltamethrin bioassay showed that larvae mortalities increased by 11.4% and 21.6%, respectively, after CYP4M51 and CYP6AB56 were partially silenced. These results suggest that inhibition of the novel cytochrome P450 genes CYP4M51 and CYP6AB56 could be used to increase the efficacy of cabbage moth control by deltamethrin.

Resistance irrelevant CYP417A2v2 was found degrading insecticide in Laodelphax striatellus

Cytochrome P450 monooxygenases (CYPs) usually overexpressed in resistant strain were found involved in oxidative detoxification of insecticides. In this study, an investigation was conducted to confirm if resistance irrelevant CYPs which were not overexpressed in resistant strain before, were capable of degrading insecticides. Three resistance irrelevant CYPs viz. CYP417A2v2, CYP425A1v2, and CYP4DJ1 from CYP4 family of Laodelphax striatellus were randomly selected for experiments. CYP417A2v2 and CYP425A1v2 were found expressed successfully in Sf9 cell line while CYP4DJ1 was not expressed successfully and out of two expressed CYPs, only CYP417A2v2 showed its efficient catalytic activity. For catalytic activity, three traditional model probe substrates and five insecticides were assayed. For the probe substrates screened, p-nitroanisole and ethoxycoumarin were preferentially metabolized by CYP417A2v2 (specific activity 3.76 ± 1.22 and 1.63 ± 0.37 nmol min-1 mg protein-1, respectively) and they may be potential diagnostic probes for this enzyme. Among insecticides, only imidacloprid was efficiently degraded by CYP417A2v2. Incubation of imidacloprid with CYP417A2v2 of L. striatellus and subsequent HPLC, LC-MS, and MS/MS analysis revealed the formation of imidacloprid metabolites, that is, 4' or 5'hydroxy-imidacloprid by hydroxylation. This result implies the exemption of CYPs character that it is not always, all the CYPs degrading insecticides being selected and overexpressed in resistant strains and the degrading CYPs without mutations to upregulate could be candidates during insecticide resistance evolution. This characterization of individual insect CYPs in insecticide degradation can provide insight for better understand of insecticide resistance development.

Leptinotarsa cap 'n' collar isoform C/Kelch-like ECH associated protein 1 signaling is critical for the regulation of ecdysteroidogenesis in the larvae

Drosophila cap 'n' collar isoform C (CncC) and Kelch-like ECH associated protein 1 (Keap1) regulate metamorphosis by transcriptional control of a subset of genes involved in ecdysteroidogenesis, 20-hydroxyecdysone (20E) signaling, and juvenile hormone (JH) degradation. In the present paper, we found that prothoracicotropic hormone signal was required for the activation of LdCncC and LdKeap1 in Leptinotarsa decemlineata. Moreover, RNA interference of LdCncC or LdKeap1 in the fourth-instar larvae delayed development. As a result, the treated larvae obtained heavier larval and pupal fresh weights and had larger body sizes than the controls. Furthermore, knockdown of LdCncC or LdKeap1 significantly reduced the mRNA levels of four ecdysone biosynthetic genes (Ldspo, Ldphm, Lddib and Ldsad), lowered 20E titer and decreased the transcript levels of five 20E response genes (LdEcR, LdUSP, LdE75, LdHR3 and LdFTZ-F1). However, the expression of two JH epoxide hydrolase genes and JH contents were not affected in the LdCncC and LdKeap1 RNAi larvae. Dietary supplementation with 20E shortened the developmental period to normal length, rescued the larval and pupal body mass rises, and recovered or even overcompensated the expression levels of the five 20E response genes in either LdCncC or LdKeap1 RNAi hypomorphs. Therefore, LdCncC/LdKeap1 signaling regulates several ecdysteroidogenesis genes, and consequently 20E pulse, to modulate the onset of metamorphosis in L. decemlineata.

Identification of glutathione S-transferase genes in Leptinotarsa decemlineata and their expression patterns under stress of three insecticides

Glutathione S-transferases (GSTs) is a family of multifunctional enzymes that are involved in detoxification of poisonous compounds. In the present paper, the Leptinotarsa decemlineata genome and transcriptome dataset were mined and 30 GST genes were identified. These GSTs belonged to cytosolic (29 genes) and microsomal (1 gene) classes. Among them 3 GSTs (LdGSTe2, LdGSTs4, and LdGSTo3) possessed splice variants. Of the 29 cytosolic LdGSTs, 3, 10, 5, 4, 4, and 1 members were classified as delta, epsilon, omega, sigma, theta, and zeta subclasses respectively, along with 2 unclassified genes. Phylogenetic analysis suggest that epsilon, omega and sigma subclasses appear to undergo species-specific bloom. Moreover, most epsilon, omega and sigma GSTs are tandemly arranged in three chromosome scaffolds. To find GST candidates involving in insecticide detoxification, we tested the mRNA levels of 20 GST transcripts under stress of cyhalothrin, fipronil or endosulfan. Out of them, LdGSTe2a, LdGSTe2b, LdGSTo5 and LdGSTt1 were significantly overexpressed after exposure to each of the three insecticides. Two other genes were respectively upregulated after cyhalothrin (LdGSTe10 and LdGSTu2) or endosulfan (LdGSTd1 and LdGSTu2) treatment. The diversified expression responses to insecticide exposure suggest that the LdGSTs may depend on a functionally complex system to detoxify different classes of insecticides. In addition, our findings provide a base for a better understanding of the evolution of insecticide resistance, and functional research on specific GST genes.

Transcriptome and Difference Analysis of Fenpropathrin Resistant Predatory Mite, Neoseiulus barkeri (Hughes)

Several fenpropathrin-resistant predatory mites have been reported. However, the molecular mechanism of the resistance remains unknown. In the present study, the Neoseiulus barkeri (N. barkeri) transcriptome was generated using the Illumina sequencing platform, 34,211 unigenes were obtained, and 15,987 were manually annotated. After manual annotation, attentions were attracted to resistance-related genes, such as voltage-gated sodium channel (VGSC), cytochrome P450s (P450s), and glutathione S-transferases (GSTs). A polymorphism analysis detected two point mutations (E1233G and S1282G) in the linker region between VGSC domain II and III. In addition, 43 putative P450 genes and 10 putative GST genes were identified from the transcriptome. Among them, two P450 genes, NbCYP4EV2 and NbCYP4EZ1, and four GST genes, NbGSTd01, NbGSTd02, NbGSTd03 and NbGSTm03, were remarkably overexpressed 3.64-46.69-fold in the fenpropathrin resistant strain compared to that in the susceptible strain. These results suggest that fenpropathrin resistance in N. barkeri is a complex biological process involving many genetic changes and provide new insight into the N. barkeri resistance mechanism.

A specialist herbivore pest adaptation to xenobiotics through up-regulation of multiple Cytochrome P450s

The adaptation of herbivorous insects to their host plants is hypothesized to be intimately associated with their ubiquitous development of resistance to synthetic pesticides. However, not much is known about the mechanisms underlying the relationship between detoxification of plant toxins and synthetic pesticides. To address this knowledge gap, we used specialist pest Colorado potato beetle (CPB) and its host plant, potato, as a model system. Next-generation sequencing (454 pyrosequencing) was performed to reveal the CPB transcriptome. Differential expression patterns of cytochrome P450 complement (CYPome) were analyzed between the susceptible (S) and imidacloprid resistant (R) beetles. We also evaluated the global transcriptome repertoire of CPB CYPome in response to the challenge by potato leaf allelochemicals and imidacloprid. The results showed that more than half (51.2%) of the CBP cytochrome P450 monooxygenases (P450s) that are up-regulated in the R strain are also induced by both host plant toxins and pesticide in a tissue-specific manner. These data suggest that xenobiotic adaptation in this specialist herbivore is through up-regulation of multiple P450s that are potentially involved in detoxifying both pesticide and plant allelochemicals.

Knockdown of a nutrient amino acid transporter gene LdNAT1 reduces free neutral amino acid contents and impairs Leptinotarsa decemlineata pupation

A Leptinotarsa decemlineata SLC6 NAT gene (LdNAT1) was cloned. LdNAT1 was highly expressed in the larval alimentary canal especially midgut. LdNAT1 mRNA levels were high right after the molt and low just before the molt. JH and a JH analog pyriproxyfen activated LdNAT1 expression. RNAi of an allatostatin gene LdAS-C increased JH and upregulated LdNAT1 transcription. Conversely, silencing of a JH biosynthesis gene LdJHAMT decreased JH and reduced LdNAT1 expression. Moreover, 20E and an ecdysteroid agonist halofenozide repressed LdNAT1 expression, whereas a decrease in 20E by RNAi of an ecdysteroidogenesis gene LdSHD and disruption of 20E signaling by knockdown of LdE75 and LdFTZ-F1 activated LdNAT1 expression. Thus, LdNAT1 responded to both 20E and JH. Moreover, knockdown of LdNAT1 reduced the contents of cysteine, histidine, isoleucine, leucine, methionine, phenylalanine and serine in the larval bodies and increased the contents of these amino acids in the larval feces. Furthermore, RNAi of LdNAT1 inhibited insulin/target of rapamycin pathway, lowered 20E and JH titers, reduced 20E and JH signaling, retarded larval growth and impaired pupation. These data showed that LdNAT1 was involved in the absorption of several neutral amino acids critical for larval growth and metamorphosis.

RNA interference suppression of the receptor tyrosine kinase Torso gene impaired pupation and adult emergence in Leptinotarsa decemlineata

In Drosophila melanogaster prothoracic gland (PG) cells, Torso mediates prothoracicotropic hormone (PTTH)-triggered mitogen activated protein kinase (MAPK) pathway (consisting of four core components Ras, Raf, MEK and ERK) to stimulate ecdysteroidogenesis. In this study, LdTorso, LdRas, LdRaf and LdERK were cloned in Leptinotarsa decemlineata. The four genes were highly or moderately expressed in the larval prothoracic glands. At the first- to third-instar stages, their expression levels were higher just before and right after the molt, and were lower in the mid instars. At the fourth-instar stage, their transcript levels were higher before prepupal stage. RNA interference-mediated knockdown of LdTorso delayed larval development, increased pupal weight, and impaired pupation and adult emergence. Moreover, knockdown of LdTorso decreased the mRNA levels of LdRas, LdRaf and LdERK, repressed the transcription of two ecdysteroidogenesis genes (LdPHM and LdDIB), lowered 20E titer, and downregulated the expression of several 20E-response genes (LdEcR, LdUSP, LdHR3 and LdFTZ-F1). Furthermore, silencing of LdTorso induced the expression of a JH biosynthesis gene LdJHAMT, increased JH titer, and activated the transcription of a JH early-inducible gene LdKr-h1. Thus, our results suggest that Torso transduces PTTH-triggered MAPK signal to regulate ecdysteroidogenesis in the PGs in a non-drosophiline insect.

CHARACTERIZATION AND FUNCTIONAL STUDY OF A PUTATIVE JUVENILE HORMONE DIOL KINASE IN THE COLORADO POTATO BEETLE Leptinotarsa decemlineata (Say)

Juvenile hormone diol kinase (JHDK) is an enzyme involved in JH degradation. In the present article, a putative JHDK cDNA (LdJHDK) was cloned from the Colorado potato beetle Leptinotarsa decemlineata. The cDNA consists of 814 bp, containing a 555 bp open reading frame encoding a 184 amino acid protein. LdJHDK reveals a high degree of identity to the previously reported insect JHDKs. It possesses three conserved purine nucleotide-binding elements, and contains three EF-hand motifs (helix-loop-helix structural domains). LdJHDK mRNA was mainly detected in hindgut and Malpighian tubules. Besides, a trace amount of LdJHDK mRNA was also found in thoracic muscles, brain-corpora cardiaca-corpora allata complex, foregut, midgut, ventral ganglia, fat body, epidermis, and hemocytes. Moreover, LdJHDK was expressed throughout all developmental stages. Within the first, second, and third larval instar, the expression levels of LdJHDK were higher just before and right after the molt, and were lower in the intermediate instar. In the fourth larval instar, the highest peak of LdJHDK occurred 56 h after ecdysis. Ingestion of double-stranded RNA (dsRNA) against LdJHDK successfully knocked down the target gene, increased JH titer, and significantly upregulated LdKr-h1 mRNA level. Knockdown of LdJHDK significantly impaired adult emergence. Thus, we provide a line of experimental evidence in L. decemlineata to support that LdJHDK encodes function protein involved in JH degradation.

Elevated expression of esterase and cytochrome P450 are related with lambda-cyhalothrin resistance and lead to cross resistance in Aphis glycines Matsumura

A resistant strain of the Aphis glycines Matsumura (CRR) has developed 76.67-fold resistance to lambda-cyhalothrin compared with the susceptible (CSS) strain. Synergists piperonyl butoxide (PBO), S,S,S-Tributyltrithiophosphate (DEF) and triphenyl phosphate (TPP) dramatically increased the toxicity of lambda-cyhalothrin to the resistant strain. Bioassay results indicated that the CRR strain had developed high levels of cross-resistance to chlorpyrifos (11.66-fold), acephate (8.20-fold), cypermethrin (53.24-fold), esfenvalerate (13.83-fold), cyfluthrin (9.64-fold), carbofuran (14.60-fold), methomyl (9.32-fold) and bifenthrin (4.81-fold), but did not have cross-resistance to chlorfenapyr, imidacloprid, diafenthiuron, abamectin. The transcriptional levels of CYP6A2-like, CYP6A14-like and cytochrome b-c1 complex subunit 9-like increased significantly in the resistant strain than that in the susceptible. Similar trend were observed in the transcripts and DNA copy number of CarE and E4 esterase. Overall, these results demonstrate that increased esterase hydrolysis activity, combined with elevated cytochrome P450 monooxygenase detoxicatication, plays an important role in the high levels of lambda-cyhalothrin resistance and can cause cross-resistance to other insecticides in the CRR strain.

A putative Δ1-pyrroline-5-carboxylate synthetase involved in the biosynthesis of proline and arginine in Leptinotarsa decemlineata

Delta 1-pyrroline-5-carboxylate synthetase (P5CS) catalyzes the conversion of glutamate (Glu) to Glu semialdehyde (GSA). GSA spontaneously cyclizes to form P5C. P5C is then reduced to proline (Pro) or is converted to ornithine, the intermediate for arginine (Arg) biosynthesis. In the present study, a full-length Ldp5cs complementary DNA was cloned from the Colorado potato beetle Leptinotarsa decemlineata, a notorious insect defoliator of potato in most potato-growing regions of the world. Ldp5cs encodes a 792-amino-acid protein which shares high identity to homologues from other insect species. Quantitative reverse transcription polymerase chain reaction revealed that Ldp5cs was ubiquitously expressed in the eggs, first to fourth-instar larvae, wandering larvae, pupae and sexually mature adults. In the adults, Ldp5cs mRNA levels were higher in the fat body, foregut, midgut and hindgut, moderate in the ventral ganglion, lower in the thorax muscles, epidermis and Malpighian tubules. Two double-stranded RNAs (dsRNAs) (dsLdp5cs1 and dsLdp5cs2) targeting Ldp5cs were constructed and bacterially expressed. Ingestion during 3 consecutive days of dsLdp5cs1 or dsLdp5cs2 successfully silenced Ldp5cs, significantly reduced the contents of Pro and Arg in the hemolymph, decreased flight speed and shortened flight distance of the resulting adults. Furthermore, knocking down Ldp5cs significantly increased adult mortality. Thus, our results suggest that identified Ldp5cs encodes a functional P5CS enzyme that is involved in the biosynthesis of Pro and Arg in L. decemlineata.

Involvement of FTZ-F1 in the regulation of pupation in Leptinotarsa decemlineata (Say)

During the final instar larvae of holometabolous insects, a pulse of 20-hydroxyecdysone (20E) and a drop in juvenile hormone (JH) trigger larval-pupal metamorphosis. In this study, two LdFTZ-F1 cDNAs (LdFTZ-F1-1 and LdFTZ-F1-2) were cloned in Leptinotarsa decemlineata. Both LdFTZ-F1-1 and LdFTZ-F1-2 were highly expressed just before or right after each molt, similar to the expression pattern of an ecdysteroidogenesis gene LdSHD. Ingestion of an ecdysteroid agonist halofenozide (Hal) enhanced LdFTZ-F1-1 and LdFTZ-F1-2 expression in the final larval instar. Conversely, a decrease in 20E by feeding a double-stranded RNA (dsRNA) against LdSHD repressed the expression. Moreover, Hal rescued the expression levels in LdSHD-silenced larvae. Thus, 20E peaks seem to induce the transcription of LdFTZ-F1s. Furthermore, ingesting dsLdFTZ-F1 from a common fragment of LdFTZ-F1-1 and LdFTZ-F1-2 successfully knocked down both LdFTZ-F1s, and impaired pupation. Finally, knocking down LdFTZ-F1s significantly repressed the transcription of three ecdysteroidogenesis genes, lowered 20E titer, and reduced the expression of two 20E receptor genes. Silencing LdFTZ-F1s also induced the expression of a JH biosynthesis gene, increased JH titer, but decreased the mRNA level of a JH early-inducible gene. Thus, LdFTZ-F1s are involved in the regulation of pupation by modulating 20E and JH titers and mediating their signaling pathways.

The P450 enzyme Shade mediates the hydroxylation of ecdysone to 20-hydroxyecdysone in the Colorado potato beetle, Leptinotarsa decemlineata

Ecdysone 20-monooxygenase (E20MO), a cytochrome P450 monooxygenase (CYP314A1), catalyses the conversion of ecdysone (E) to 20-hydroxyecdysone (20E). We report here the cloning and characterization of the Halloween gene Shade (Shd) encoding E20MO in the Colorado potato beetle, Leptinotarsa decemlineata. LdSHD has five conserved motifs typical of insect P450s, ie the Helix-C, Helix-I, Helix-K, PxxFxPE/DRF (PERF) and heme-binding motifs. LdShd was expressed in developing eggs, the first to fourth instars, wandering larvae, pupae and adults, with statistically significant fluctuations. Its mRNA was ubiquitously distributed in the head, thorax and abdomen. The recombinant LdSHD protein expressed in Spodoptera frugiperda 9 (Sf9) cells catalysed the conversion of E to 20E. Dietary introduction of double-stranded RNA (dsRNA) of LdShd into the second instar larvae successfully knocked down the LdShd expression level, decreased the mRNA level of the ecdysone receptor (LdEcR) gene, caused larval lethality, delayed development and affected pupation. Moreover, ingestion of LdShd-dsRNA by the fourth instars also down-regulated LdShd and LdEcR expression, reduced the 20E titre, and negatively influenced pupation. Introduction of 20E and a nonsteroidal ecdysteroid agonist halofenozide into the LdShd-dsRNA-ingested second instars, and of halofenozide into the LdShd-dsRNA-ingested fourth instars almost completely relieved the negative effects on larval performance. Thus, LdSHD functions to regulate metamorphotic processes by converting E to 20E in a coleopteran insect species Le. decemlineata.

De novo transcriptome analysis of Perna viridis highlights tissue-specific patterns for environmental studies

BACKGROUND

The tropical green-lipped mussel Perna viridis is a common biomonitor throughout the Indo-Pacific region that is used for environmental monitoring and ecotoxicological investigations. However, there is limited molecular data available regarding this species. We sought to establish a global transcriptome database from the tissues of adductor muscle, gills and the hepatopancreas of P. viridis in an effort to advance our understanding of the molecular aspects involved during specific toxicity responses in this sentinel species.

RESULTS

Illumina sequencing results yielded 544,272,542 high-quality filtered reads. After de novo assembly using Trinity, 233,257 contigs were generated with an average length of 1,264 bp and an N50 length of 2,868 bp; 192,879 assembled transcripts and 150,111 assembled unigenes were obtained after clustering. A total of 93,668 assembled transcripts (66,692 assembled genes) with putative functions for protein domains were predicted based on InterProScan analysis. Based on similarity searches, 44,713 assembled transcripts and 25,319 assembled unigenes were annotated with at least one BLAST hit. A total of 21,262 assembled transcripts (11,947 assembled genes) were annotated with at least one well-defined Gene Ontology (GO) and 5,131 assembled transcripts (3,181 assembled unigenes) were assigned to 329 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The quantity of assembled unigenes and transcripts obtained from male and female mussels were similar but varied among the three studied tissues, with the highest numbers recorded in the gills, followed by the hepatopancreas, and then the adductor muscle. Multivariate analyses revealed strong tissue-specific patterns among the three different tissues, but not between sexes in terms of expression profiles for annotated genes in various GO terms, and genes associated with stress responses and degradation of xenobiotics. The expression profiles of certain selected genes in each tissue type were further validated using real-time quantitative polymerase chain reaction assays and a similar tissue-specific trend was seen.

CONCLUSIONS

The extensive sequence data generated from this study will provide a valuable molecular resource for facilitating environmental studies with P. viridis, and highlight the importance of tissue-specific approaches in the future.