Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae)

Characteristics of Trypsin genes and their roles in insecticide resistance based on omics and functional analyses in the malaria vector Anopheles sinensis

Trypsin is one of the most diverse and widely studied protease hydrolases. However, the diversity and characteristics of the Trypsin superfamily of genes have not been well understood, and their role in insecticide resistance is yet to be investigated. In this study, a total of 342 Trypsin genes were identified and classified into seven families based on homology, characteristic domains and phylogenetics in Anopheles sinensis, and the LY-Domain and CLECT-Domain families are specific to the species. Four Trypsin genes, (Astry2b, Astry43a, Astry90, Astry113c) were identified to be associated with pyrethroid resistance based on transcriptome analyses of three field resistant populations and qRT-PCR validation, and the knock-down of these genes significantly decrease the pyrethroid resistance of Anopheles sinensis based on RNAi. The activity of Astry43a can be reduced by five selected insecticides (indoxacarb, DDT, temephos, imidacloprid and deltamethrin); and however, the Astry43a could not directly metabolize these five insecticides, like the trypsin NYD-Tr did in earlier reports. This study provides the overall information frame of Trypsin genes, and proposes the role of Trypsin genes to insecticide resistance. Further researches are necessary to investigate the metabolism function of these trypsins to insecticides.

The role of P450 enzymes in malaria and other vector-borne infectious diseases

Vector-borne infectious diseases are still an important global health problem. Malaria is the most important among them, mainly pediatric, life-threatening disease. Malaria and other vector-borne disorders caused by parasites, bacteria, and viruses have a strong impact on public health and significant economic costs. Most vector-borne diseases could be prevented by vector control, with attention to the ecological and biodiversity conservation aspects. Chemical control with pesticides and insecticides is widely used as a measure of prevention although increasing resistance to insecticides is a serious issue in vector control. Metabolic resistance is the most common mechanism and poses a big challenge. Insect enzyme systems, including monooxygenase CYP P450 enzymes, are employed by vectors mainly to metabolize insecticides thus causing resistance. The discovery and application of natural specific inhibitors/blockers of vector P450 enzymes as synergists for commonly used pesticides will contribute to the "greening" of insecticides. Besides vector CYPs, host CYP enzymes could also be exploited to fight against vector-borne diseases: using mostly their detoxifying properties and involvement in the immune response. Here, we review published research data on P450 enzymes from all players in vector-borne infections, that is, pathogens, vectors, and hosts, regarding the potential role of CYPs in disease. We discuss strategies on how to exploit cytochromes P450 in vector-borne disease control.

Transcriptomic Analysis Reveals the Detoxification Mechanism of Chilo suppressalis in Response to the Novel Pesticide Cyproflanilide

Chilo suppressalis is one of the most damaging rice pests in China’s rice-growing regions. Chemical pesticides are the primary method for pest control; the excessive use of insecticides has resulted in pesticide resistance. C. suppressalis is highly susceptible to cyproflanilide, a novel pesticide with high efficacy. However, the acute toxicity and detoxification mechanisms remain unclear. We carried out a bioassay experiment with C. suppressalis larvae and found that the LD10, LD30 and LD50 of cyproflanilide for 3rd instar larvae was 1.7 ng/per larvae, 6.62 ng/per larvae and 16.92 ng/per larvae, respectively. Moreover, our field trial results showed that cyproflanilide had a 91.24% control efficiency against C. suppressalis. We investigated the effect of cyproflanilide (LD30) treatment on the transcriptome profiles of C. suppressalis larvae and found that 483 genes were up-regulated and 305 genes were down-regulated in response to cyproflanilide exposure, with significantly higher CYP4G90 and CYP4AU10 expression in the treatment group. The RNA interference knockdown of CYP4G90 and CYP4AU10 increased mortality by 20% and 18%, respectively, compared to the control. Our results indicate that cyproflanilide has effective insecticidal toxicological activity, and that the CYP4G90 and CYP4AU10 genes are involved in detoxification metabolism. These findings provide an insight into the toxicological basis of cyproflanilide and the means to develop efficient resistance management tools for C. suppressalis.

Identification and Functional Characterization of CYP4D2 Putatively Associated with β-Cypermethrin Detoxification in Phortica okadai

Phortica okadai, a polyphagous pest, serves as a vector for Thelazia callipaeda in China. Currently, there are no effective control strategies for this vector. Agricultural pest control may cause P. okadai to become a threat due to the development of pesticide resistance. Cytochrome P450s (CYP450) plays a significant role in detoxifying xenobiotics in insects. In this study, we performed RNA sequencing of P. okadai exposed to β-cypermethrin for 0 and 1 h and then gene cloning of the five up-regulated CYP450 genes. Three CYP450 genes were successfully cloned, and their expression patterns in different developmental stages and in different tissues were analyzed by RT-qPCR. Pocyp4d2 was observed to have the highest expression in the midgut (fold change 2.82 for Pocyp4d2, 2.62 for Pocyp49a1, and 1.77 for Pocyp28d2). Functional analysis was carried out according to overexpression in S2 cells from the pfastbac1 vector and RNAi with siRNA. The results of the CCK8 assay indicated that the overexpression of the recombinant protein PoCYP4D2 suppressed the decrease in S2 cell viability due to β-cypermethrin. The expression levels of PoCYP4D2 decreased significantly, and the mortality rates increased from 6.25% to 15.0% at 3 h and from 15.0% to 27.5% at 6 h after Pocyp4d2-siRNA injection. These results suggest that Pocyp4d2 may be an essential key gene in the metabolism of β-cypermethrin in P. okadai. This study constitutes a foundation to explore further the functions of P. okadai CYP450 genes in insecticide metabolism.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

Glutathione S-transferase (GST) genes and their function associated with pyrethroid resistance in the malaria vector Anopheles sinensis

BACKGROUND

Glutathione S-transferases (GSTs), a multifunctional protein family, are involved in insecticide resistance. However, a systematic analysis of GSTs in Anopheles sinensis, an important vector for malaria transmission, is lacking. In this study, we investigated the diversity and characteristics of GST genes, and analyzed their expression patterns and functions associated with insecticide resistance in this species.

RESULTS

We identified 32 putative cytosolic and three putative microsomal GST genes in the An. sinensis genome. Transcriptome analysis showed that GSTs were highly expressed in larvae, and mainly expressed in the antennae, midgut and Malpighian tubules of adults. In addition, we found that GSTd2 and GSTe2 were significantly upregulated in four An. sinensis pyrethroid-resistant field populations. Furthermore, silencing of GSTd2 and GSTe2 significantly increased the susceptibility of An. sinensis to deltamethrin, and recombinant GSTd2 and GSTe2 exhibited high enzymatic activity in the metabolism of 1-chloro-2,4-dinitrobenzene and dichlorodiphenyltrichloroethane (DDT).

CONCLUSION

These results showed that GSTs are involved in the development of insecticide resistance in An. sinensis through transcriptional overexpression and enzymatic metabolization, facilitating our understanding of insecticide resistance in insects. © 2022 Society of Chemical Industry.

Genome-Wide Identification of P450 Genes in Chironomid Propsilocerus akamusi Reveals Candidate Genes Involved in Gut Microbiota-Mediated Detoxification of Chlorpyrifos

Chironomids commonly dominate macroinvertebrate assemblages in aquatic habitats and these non-biting midges are known for their ability to tolerate contaminants. Studies regarding the interplay between gut microbiota and host detoxification ability is currently a point of interest. Cytochrome P450s (P450s) are critical metabolic enzymes in which a subset is involved in xenobiotic detoxification. In this study, we first conducted an integrated global investigation of P450s based on the whole genomic sequence of Propsilocerus akamusi and retrieved a series of 64 P450 genes which were further classified into 4 clans and 25 families on the basis of phylogenetic relationships. With assistance of RNA-Seq and RT-qPCR validation, the expression profile of screened PaP450s in guts was compared between chlorpyrifos-challenged larvae with deficient gut microbiota (GD) and those with a conventional gut community (CV). An increasing prevalence of chlorpyrifos from sublethal to lethal dosages induced a greater mortality rate of individuals coupled with remarkable downregulation of 14 P450s in GD larval guts when compared to CV ones. Moreover, it turned out that the decreased level of PaCYP3998B1 and PaCYP3987D1 might imply impaired host endogenous detoxification capability potentiated by gut dysbiosis, reflected by a remarkably severe mortality in GD larvae treated with lethal chlorpyrifos. Collectively, our study unveiled candidate P450 genes that might be mediated by gut symbionts in chlorpyrifos-challenged P. akamusi larvae, possibly facilitating further understanding of the detoxified mechanism that chironomids might employ to alleviate poisonousness.

A Blueprint of Microstructures and Stage-Specific Transcriptome Dynamics of Cuticle Formation in Bombyx mori

Insect cuticle is critical for the environmental adaptability and insecticide resistance of insects. However, there is no clear understanding of the structure and protein components of the cuticle during each developmental stage of holometabolous insects, and knowledge about the protein components within each layer is vague. We conducted serial sectioning, cuticular structure analysis, and transcriptome sequencing of the larval, pupal, and adult cuticles of Bombyx mori. The deposition processes of epicuticle, exocuticle, and endocuticle during larval, pupal, and adult cuticle formation were similar. Transcriptome analysis showed that these cuticle formations share 74% of the expressed cuticular protein (CP) genes and 20 other structural protein genes, such as larval serum protein and prisilkin. There are seven, six, and eleven stage-specific expressed CP genes in larval, pupal, and adult cuticles, respectively. The types and levels of CP genes may be the key determinants of the properties of each cuticular layer. For example, the CPs of the RR-2 protein family with high contents of histidine (His) are more essential for the exocuticle. Functional analysis suggested that BmorCPAP1-H is involved in cuticle formation. This study not only offers an in-depth understanding of cuticle morphology and protein components but also facilitates the elucidation of molecular mechanisms underlying cuticle formation in future studies.

Comparative transcriptome and RNA interference reveal CYP6DC1 and CYP380C47 related to lambda-cyhalothrin resistance in Rhopalosiphum padi

The bird-cherry-oat aphid, Rhopalosiphum padi, is a serious agricultural pest of Triticeae crops, and pyrethroids are the most widely used chemical pesticides for the control of the aphid. Our previous studies found that some R. padi field populations have developed resistance against pyrethroids; an M918L target-site mutation of the voltage gated sodium channel was present in the pyrethroid resistant individuals, while the high-level resistance to lambda-cyhalothrin revealed the presence of other mechanisms in the pest. Here, we conducted genome-wide transcriptional analysis for the lambda-cyhalothrin susceptible (SS) and resistant (LC-RR) strains of R. padi. Results indicated that 2457 genes were differently expressed between the SS and LC-RR strains. In the LC-RR, a total of 1265 and 1192 genes were up- and down-regulated, respectively. KEGG analysis implicated enrichment of P450 involved in insecticide metabolic pathways in the resistant transcriptome. qRT-PCR results confirmed that two P450 genes (CYP6DC1 and CYP380C47) were significantly overexpressed in the LC-RR individuals. Furthermore, RNA interference (RNAi) of CYP6DC1 or CYP380C47 significantly increased mortality of R. padi exposure to lambda-cyhalothrin. These results suggest that the overexpression of CYP6DC1 and CYP380C47 contributed to the lambda-cyhalothrin resistance in the pest. This study provides knowledge for further analyzing the molecular mechanism of resistance to pyrethroids in R. padi.

Genetic variation at the Cyp6m2 putative insecticide resistance locus in Anopheles gambiae and Anopheles coluzzii

Background

The emergence of insecticide resistance is a major threat to malaria control programmes in Africa, with many different factors contributing to insecticide resistance in its vectors, Anopheles mosquitoes. CYP6M2 has previously been recognized as an important candidate in cytochrome P450-mediated detoxification in Anopheles. As it has been implicated in resistance against pyrethroids, organochlorines and carbamates, its broad metabolic activity makes it a potential agent in insecticide cross-resistance. Currently, allelic variation within the Cyp6m2 gene remains unknown.

Methods

Here, Illumina whole-genome sequence data from Phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) was used to examine genetic variation in the Cyp6m2 gene across 16 populations in 13 countries comprising Anopheles gambiae and Anopheles coluzzii mosquitoes. To identify whether these alleles show evidence of selection either through potentially modified enzymatic function or by being linked to variants that change the transcriptional profile of the gene, hierarchical clustering of haplotypes, linkage disequilibrium, median joining networks and extended haplotype homozygosity analyses were performed.

Results

Fifteen missense biallelic substitutions at high frequency (defined as > 5% frequency in one or more populations) are found, which fall into five distinct haplotype groups that carry the main high frequency variants: A13T, D65A, E328Q, Y347F, I359V and A468S. Despite consistent reports of Cyp6m2 upregulation and metabolic activity in insecticide resistant Anophelines, no evidence of directional selection is found occurring on these variants or on the haplotype clusters in which they are found.

Conclusion

These results imply that emerging resistance associated with Cyp6m2 is potentially driven by distant regulatory loci such as transcriptional factors rather than by its missense variants, or that other genes are playing a more significant role in conferring metabolic resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03757-4.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

MiR-4448 is involved in deltamethrin resistance by targeting CYP4H31 in Culex pipiens pallens

Background

Culex pipiens (Cx. pipiens) complex, which acts as a vector of viruses and is widespread and abundant worldwide, including West Nile virus, Japanese encephalitis virus, and Sindbis virus, can cause serious vector-borne diseases affecting human health. Unfortunately, mosquitoes have developed deltamethrin resistance because of its long-term overuse, representing a major challenge to mosquito control. Understanding the molecular regulatory mechanisms of resistance is vital to control mosquitoes. MicroRNAs (miRNAs) are short non-coding RNAs that have been demonstrated to be important regulators of gene expression across a wide variety of organisms, which might function in mosquito deltamethrin resistance. In the present study, we aimed to investigate the regulatory functions of miR-4448 and CYP4H31 in the formation of insecticidal resistance in mosquito Culex pipiens pallens.

Methods

We used quantitative real-time reverse transcription PCR to measure miR-4448 and CYP4H31 (encoding a cytochrome P450) expression levels. The regulatory functions of miR-4448 and CYP4H31 were assessed using dual-luciferase reporter assays. Then, oral feeding, RNA interference, and the American Centers for Disease Control and Prevention bottle bioassay were used to determine miR-4448’s association with deltamethrin resistance by targeting CYP4H31in vivo. Cell Counting Kit-8 (CCK-8) was also used to detect the viability of pIB/V5-His-CYP4H31-transfected C6/36 cells after deltamethrin treatment in vitro.

Results

MiR-4448 was downregulated in the deltamethrin-resistant strain (DR strain), whereas CYP4H31 was downregulated in deltamethrin-susceptible strain. CYP4H31 expression was downregulated by miR-4448 recognizing and binding to its 3′ untranslated region. Functional verification experiments showed that miR-4448 overexpression resulted in lower expression of CYP4H31. The mortality of miR-4448 mimic-injected DR strain mosquitoes was higher than that of the controls. CCK-8 assays showed that CYP4H31 decreased cellular resistance to deltamethrin in vitro and the mortality of the DR strain increased when CYP4H31 was knocked down in vivo.

Conclusions

In mosquitoes, miR-4448 participates in deltamethrin resistance by targeting CYP4H31. The results of the present study increase our understanding of deltamethrin resistance mechanisms.

Transcriptome profiling and functional analysis suggest that the constitutive overexpression of four cytochrome P450s confers resistance to abamectin in Tetranychus urticae from China

BACKGROUND

The two-spotted spider mite Tetranychus urticae is a polyphagous and cosmopolitan pest that has developed high resistance to abamectin, making it difficult to control. Although 'target resistance' related to glutamate-gated chloride channel mutations was found in T. urticae field populations in China, other resistance mechanisms appear to be involved. Here, we conducted genome-wide transcriptome profiling using RNA-sequencing of two abamectin-resistant populations (NB-ZJ and SY-BJ) and one susceptible strain (Lab-SS) to identify differentially expressed genes that might contribute to the resistance of T. urticae to abamectin in China.

RESULTS

Our experiments showed that abamectin resistance was synergized by piperonyl butoxide (PBO) and triphenyl phosphate (TPP), with synergistic ratios (SR) of 2.95-fold and 2.21-fold for PBO and 3.55-fold and 2.84-fold for TPP in NB-ZJ and SY-BJ populations, respectively. Transcriptome data and quantitative real-time PCR (qRT-PCR) revealed that seven detoxification enzyme genes were overexpressed in the two resistant populations. Furthermore, functional analysis by RNA interference (RNAi) indicated that the mortality caused by abamectin was significantly increased by the separate silencing of the P450 genes CYP389C10, CYP392D8, CYP392A11, and CYP392A12.

CONCLUSION

qRT-PCR expression and RNAi data suggest that the overexpression of P450 genes CYP389C10, CYP392D8, CYP392A11, and CYP392A12 may be involved in the abamectin-resistance of field populations of T. urticae in China. This knowledge could facilitate the elucidation of resistance mechanisms and the development of resistance management of T. urticae field populations. © 2020 Society of Chemical Industry.

Citronellal perception and transmission by Anopheles gambiae s.s. (Diptera: Culicidae) females

Anopheles gambiae s.s. is a key vector of Plasmodium parasites. Repellents, which may be a promising alternative to pesticides used to control malaria mosquitoes. Although citronellal is a known mosquito repellent, its repellency characteristics are largely unknown. Determining the specific odorant-binding proteins (OBPs) and odorant receptors (ORs) that detect and transfer the citronellal molecule in A. gambiae s.s. will help to define the mode of action of this compound. In this research, we assessed the repellent activity of citronellal in A. gambiae s.s. using a Y-tube olfactory meter, screened candidate citronellal-binding OBPs and ORs using reverse molecular docking, clarified the binding properties of predicted proteins for citronellal using fluorescence competition binding assay. Results showed that citronellal had a dosage effect on repelling A. gambiae s.s.. The 50% repellent rate was determined to be 4.02 nmol. Results of simulated molecular docking showed that the only proteins that bound tightly with citronellal were AgamOBP4 and AgamORC7. Fluorescence competitive binding assays confirmed the simulations. This research determined that citronellal was captured by AgamOBP4 and transmitted to AgamORC7 in A. gambiae s.s.. Our study will be beneficial in the further understanding the repellent mechanism of citronellal against A. gambiae s.s..

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Bactrocera dorsalis (Hendel) and screening of candidate P450 genes associated with malathion resistance

BACKGROUND

Bactrocera dorsalis (Hendel) is a notorious agricultural pest worldwide, and its resistance to insecticides is a major obstacle in successful control. Cytochrome P450s (P450s) are major metabolic enzymes associated with insecticide resistance. The genome of B. dorsalis was sequenced recently, allowing an integrated genome-wide analysis of P450 genes (P450s) and the analysis of correlations between these genes and insecticide resistance in this pest.

RESULTS

Totally, 101 P450s were identified in the B. dorsalis genome and classified into four clans, 25 families and 57 subfamilies. Quantitative reverse transcription polymerase chain reaction results showed that most of these genes were highly expressed in adults (46) and in metabolic tissues, including the fatbody (63), midgut (61) and Malphagian tubules (66). In a malathion-resistant strain, 13 and 9 genes were significantly upregulated and downregulated, respectively, compared with a susceptible strain, and these genes were screened as candidate genes associated with malathion resistance.

CONCLUSION

This study provides useful information for understanding the evolution and potential functions of P450s in B. dorsalis, and the results lay the foundation for further studies on the correlations between P450s and malathion resistance in B. dorsalis. © 2020 Society of Chemical Industry.

Deciphering pyrethroid resistance in Cx. pipiens (L): Implications of cytochrome P450; expression profiling and regulatory microRNA

Over the past decades, the extensive use of pyrethroids insecticides for vector control has resulted in the development of insecticide resistance. Cytochrome P450 has been recognized to play a critical role in the metabolic detoxification of insecticides. In the current study, Culex pipiens mosquitoes were collected from Giza Governorate in Egypt and tested for insecticide susceptibility against deltamethrin. First detection of Knockdown resistance gene (Kdr) mutations in field collected mosquitoes was performed. Activities of cytochrome oxidase P450 detoxification enzyme that synchronized with the resistance development, was assessed. Expression profiles of cytochrome P450s and their putative corresponding regulating miRNAs, which was previously reported in Cx. pipiens pallens were evaluated in pyrethroid resistant field-collected Cx. pipiens using RT-qPCR and stem-loop RT-qPCR, respectively. Specific stem-loop reverse transcription primers and forward primers were designed for miRNAs profiling. Our results elucidated the pyrethroid resistance development and revealed its relation to the metabolic and target site modification mechanisms with a first report of L1014F-kdr mutation detection. RT-qPCR results have showed an up-regulation in the expression of the studied P450 transcripts. Negative correlations were found between the expression of P450s and their regulatory miRNAs except for CYP9J35, where positive correlation was found with its corresponding miR-13. Interestingly, our data was the first to detect negative correlation between miR-285 and its putative CYP6Cp1 target gene. These findings highlighted the significance of identifying P450 gene along with regulatory miRNAs as a key mechanism implicated in pyrethroid resistance in field Culex vector population. The elucidation of this mechanism would shed light on the development of insecticide resistance and would help in shaping strategies to combat such vectors.

Mitochondrial genes associated with pyrethroid resistance revealed by mitochondrial genome and transcriptome analyses in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Insecticides are still the main method of mosquito control, but mosquito resistance presents a large obstacle. The function of mitochondrial genes in the evolution of insecticide resistance is still poorly understood. Pyrethroid is the most commonly used insecticide, and Anopheles sinensis is an important malaria vector in China and Southeast Asia. In this study, we investigated the mitochondrial genes associated with pyrethroid resistance through their genetic and expression variation based on analyses of transcriptomes and 36 individuals with resequencing in three geographical populations in China.

RESULTS

The nucleotide diversity (Pi) in 18 resistant individuals was much lower than that in 18 susceptible individuals, which suggests that some sites experienced purifying selection subject to pyrethroid stress. Ka/Ks and amino acid analyses showed that ND4 experienced positive selection and had 23 amino acid mutations due to pyrethroid stress. These mutations might change the ND4 structure and function and thus alter the efficiency of the respiratory chain. ND5 was significantly upregulated, and ATP8 was significantly downregulated in these three pyrethroid resistant populations, which suggests that these two genes function in the production and maintenance of pyrethroid resistance. There are differences in mitochondrial genes involved in pyrethroid resistance among these three populations.

CONCLUSION

This is the first study to reveal the association of mitochondrial genes in the evolution of insecticide resistance through amino acid mutation and expression patterns and can help us further understand insecticide resistance mechanisms. © 2019 Society of Chemical Industry.

The cytochrome P450 CYP389C16 contributes to the cross-resistance between cyflumetofen and pyridaben in Tetranychus cinnabarinus (Boisduval)

BACKGROUND

Acaricide resistance is a serious problem in spider mites. Cyflumetofen is a new complex II inhibitor, whereas pyridaben acts at complex I and has been used for decades. Although cross-resistance between cyflumetofen and pyridaben has been observed in Tetranychus cinnabarinus, the specific mechanisms at play have not yet been investigated.

RESULTS

Investigation into the cross-resistance mechanisms identified five P450s, among which CYP389C16 was evaluated as the most likely candidate conferring cross-resistance. Knockdown of CYP389C16 expression via RNA interference diminished the level of cross-resistance in the cyflumetofen-resistant strain. In addition, recombinant CYP389C16 (40 pmol) effectively metabolized 25.0 ± 0.7% of cyflumetofen, 39.7 ± 1.0% of pyridaben, and 69.3 ± 3.3% of AB-1 (active de-esterified metabolite of cyflumetofen) within 2 h. In addition, hydroxylation metabolite of AB-1 was identified by HPLC-MS/MS.

CONCLUSIONS

The study reveals that overexpressed CYP389C16 is involved in the cross-resistance between cyflumetofen and pyridaben in T. cinnabarinus. © 2019 Society of Chemical Industry.

Comparative analyses of simple sequence repeats (SSRs) in 23 mosquito species genomes: Identification, characterization and distribution (Diptera: Culicidae)

Simple sequence repeats (SSRs) exist in both eukaryotic and prokaryotic genomes and are the most popular genetic markers, but the SSRs of mosquito genomes are still not well understood. In this study, we identified and analyzed the SSRs in 23 mosquito species using Drosophila melanogaster as reference at the whole-genome level. The results show that SSR numbers (33 076-560 175/genome) and genome sizes (574.57-1342.21 Mb) are significantly positively correlated (R² = 0.8992, P < 0.01), but the correlation in individual species varies in these mosquito species. In six types of SSR, mono- to trinucleotide SSRs are dominant with cumulative percentages of 95.14%-99.00% and densities of 195.65/Mb-787.51/Mb, whereas tetra- to hexanucleotide SSRs are rare with 1.12%-4.22% and 3.76/Mb-40.23/Mb. The (A/T)n, (AC/GT)n and (AGC/GCT)n are the most frequent motifs in mononucleotide, dinucleotide and trinucleotide SSRs, respectively, and the motif frequencies of tetra- to hexanucleotide SSRs appear to be species-specific. The 10-20 bp length of SSRs are dominant with the number of 110 561 ± 93 482 and the frequency of 87.25% ± 5.73% on average, and the number and frequency decline with the increase of length. Most SSRs (83.34% ± 7.72%) are located in intergenic regions, followed by intron regions (11.59% ± 5.59%), exon regions (3.74% ± 1.95%), and untranslated regions (1.32% ± 1.39%). The mono-, di- and trinucleotide SSRs are the main SSRs in both gene regions (98.55% ± 0.85%) and exon regions (99.27% ± 0.52%). An average of 42.52% of total genes contains SSRs, and the preference for SSR occurrence in different gene subcategories are species-specific. The study provides useful insights into the SSR diversity, characteristics and distribution in 23 mosquito species of genomes.

HSP superfamily of genes in the malaria vector Anopheles sinensis: diversity, phylogenetics and association with pyrethroid resistance

BACKGROUND

Heat shock proteins (HSPs) are molecular chaperones that are involved in many normal cellular processes and various kinds of environmental stress. There is still no report regarding the diversity and phylogenetics research of HSP superfamily of genes at whole genome level in insects, and the HSP gene association with pyrethroid resistance is also not well known. The present study investigated the diversity, classification, scaffold location, characteristics, and phylogenetics of the superfamily of genes in Anopheles sinensis genome, and the HSP genes associated with pyrethroid resistance.

METHODS

The present study identified the HSP genes in the An. sinensis genome, analysed their characteristics, and deduced phylogenetic relationships of all HSPs in An. sinensis, Anopheles gambiae, Culex quinquefasciatus and Aedes aegypti by bioinformatic methods. Importantly, the present study screened the HSPs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and looked over the HSP gene expression pattern for the first time in An. sinensis on the time-scale post insecticide treatment with RT-qPCR.

RESULTS

There are 72 HSP genes in An. sinensis genome, and they are classified into five families and 11 subfamilies based on their molecular weight, homology and phylogenetics. Both RNA-seq and qPCR analysis revealed that the expression of AsHSP90AB, AsHSP70-2 and AsHSP21.7 are significantly upregulated in at least one field pyrethroid-resistant population. Eleven genes are significantly upregulated in different period after pyrethroid exposure. The HSP90, sHSP and HSP70 families are proposed to be involved in pyrethroid stress response based in expression analyses of three field pyrethroid-resistant populations, and expression pattern on the time scale post insecticide treatment. The AsHSP90AB gene is proposed to be the essential HSP gene for pyrethroid stress response in An. sinensis.

CONCLUSIONS

This study provides the information frame for HSP superfamily of genes, and lays an important basis for the better understanding and further research of HSP function in insect adaptability to diverse environments.

ATP-Binding Cassette (ABC) Transporter Genes Involved in Pyrethroid Resistance in the Malaria Vector Anopheles sinensis: Genome-Wide Identification, Characteristics, Phylogenetics, and Expression Profile

background: The ATP-binding cassette (ABC) transporters family is one of the largest families of membrane proteins existing in all living organisms. Pyrethroid resistance has become the largest unique obstacle for mosquito control worldwide. ABC transporters are thought to be associated with pyrethroid resistance in some agricultural pests, but little information is known for mosquitoes. Herein, we investigated the diversity, location, characteristics, phylogenetics, and evolution of ABC transporter family of genes in the Anopheles sinensis genome, and identified the ABC transporter genes associated with pyrethroid resistance through expression profiles using RNA-seq and qPCR. Results: 61 ABC transporter genes are identified and divided into eight subfamilies (ABCA-H), located on 22 different scaffolds. Phylogenetic and evolution analyses with ABC transporters of A. gambiae, Drosophila melanogaster, and Homo sapiens suggest that the ABCD, ABCG, and ABCH subfamilies are monophyly, and that the ABCC and ABCG subfamilies have experienced a gene duplication event. Both RNA-seq and qPCR analyses show that the AsABCG28 gene is uniquely significantly upregulated gene in all three field pyrethroid-resistant populations (Anhui, Chongqing, and Yunnan provinces) in comparison with a laboratory-susceptible strain from Jiangsu province. The AsABCG28 is significantly upregulated at 12-h and 24-h after deltamethrin exposure in three-day-old female adults. Conclusion: This study provides the information frame for ABC transporter subfamily of genes, and lays an important basis for the better understanding and further research of ABC transporter function in insecticide toxification. The AsABCG28 gene is associated with pyrethroid detoxification, and it functions at later period in the detoxification process for xenobiotics transportation.

Benzo(a)pyrene degradation by cytochrome P450 hydroxylase and the functional metabolism network of Bacillus thuringiensis

The relationship between benzo(a)pyrene biodegradation and certain target biomolecules has been investigated. To regulate the degradation process, the associated metabolism network must be clarified. To this end, benzo(a)pyrene degradation, carbon substrate metabolism and exometabolomic mechanism of Bacillus thuringiensis were analyzed. Benzo(a)pyrene was degraded through hydroxylation catalyzed by cytochrome P450 hydroxylase. After the treatment of 0.5 mg L^-1 of benzo(a)pyrene by 0.2 g L^-1 of cells for 9 d, biosorption and degradation efficiencies were measured at approximately 90% and 80%, respectively. During this process, phospholipid synthesis, glycogen, asparagine, arginine, itaconate and xylose metabolism were significantly downregulated, while glycolysis, pentose phosphate pathway, citrate cycle, amino sugar and nucleotide sugar metabolism were significantly upregulated. These findings offer insight into the biotransformation regulation of polycyclic aromatic hydrocarbons.

UDP-glycosyltransferase genes and their association and mutations associated with pyrethroid resistance in Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

UDP-glycosyltransferase (UGT) is an important biotransformation superfamily of enzymes. They catalyze the transfer of glycosyl residues from activated nucleotide sugars to acceptor hydrophobic molecules, and function in several physiological processes, including detoxification, olfaction, cuticle formation, pigmentation. The diversity, classification, scaffold location, characteristics, phylogenetics, and evolution of the superfamily of genes at whole genome level, and their association and mutations associated with pyrethroid resistance are still little known.

METHODS

The present study identified UGT genes in Anopheles sinensis genome, classified UGT genes in An. sinensis, Anopheles gambiae, Aedes aegypti and Drosophila melanogaster genomes, and analysed the scaffold location, characteristics, phylogenetics, and evolution of An. sinensis UGT genes using bioinformatics methods. The present study also identified the UGTs associated with pyrethroid resistance using three field pyrethroid-resistant populations with RNA-seq and RT-qPCR, and the mutations associated with pyrethroid resistance with genome re-sequencing in An. sinensis.

RESULTS

There are 30 putative UGTs in An. sinensis genome, which are classified into 12 families (UGT301, UGT302, UGT306, UGT308, UGT309, UGT310, UGT313, UGT314, UGT315, UGT36, UGT49, UGT50) and further into 23 sub-families. The UGT308 is significantly expanded in gene number compared with other families. A total of 119 UGTs from An. sinensis, An. gambiae, Aedes aegypti and Drosophila melanogaster genomes are classified into 19 families, of which seven are specific for three mosquito species and seven are specific for Drosophila melanogaster. The UGT308 and UGT302 are proposed to main families involved in pyrethroid resistance. The AsUGT308D3 is proposed to be the essential UGT gene for the participation in biotransformation in pyrethroid detoxification process, which is possibly regulated by eight SNPs in its 3' flanking region. The UGT302A3 is also associated with pyrethroid resistance, and four amino acid mutations in its coding sequences might enhance its catalytic activity and further result in higher insecticide resistance.

CONCLUSIONS

This study provides the diversity, phylogenetics and evolution of UGT genes, and potential UGT members and mutations involved in pyrethroid resistance in An. sinensis, and lays an important basis for the better understanding and further research on UGT function in defense against insecticide stress.

Complete mitogenome of Anopheles sinensis and mitochondrial insertion segments in the nuclear genomes of 19 mosquito species

Anopheles sinensis is a major malarial vector in China and Southeast Asia. The mitochondria is involved in many important biological functions. Nuclear mitochondrial DNA segments (NUMTs) are common in eukaryotic organisms, but their characteristics are poorly understood. We sequenced and analyzed the complete mitochondrial (mt) genome of An. sinensis. The mt genome is 15,418 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs and a large non-coding region. Its gene arrangement is similar to previously published mosquito mt genomes. We identified and analyzed the NUMTs of 19 mosquito species with both nuclear genomes and mt genome sequences. The number, total length and density of NUMTs are significantly correlated with genome size. About half of NUMTs are short (< 200 bp), but larger genomes can house longer NUMTs. NUMTs may help explain genome size expansion in mosquitoes. The expansion due to mitochondrial insertion segments is variable in different insect groups. PCGs are transferred to nuclear genomes at a higher frequency in mosquitoes, but NUMT origination is more different than in mammals. Larger-sized nuclear genomes have longer mt genome sequences in both mosquitoes and mammals. The study provides a foundation for the functional research of mitochondrial genes in An. sinensis and helps us understand the characteristics and origin of NUMTs and the potential contribution to genome expansion.

The Effect of Permethrin Resistance on Aedes aegypti Transcriptome Following Ingestion of Zika Virus Infected Blood

Aedes aegypti (L.) is the primary vector of many emerging arboviruses. Insecticide resistance among mosquito populations is a consequence of the application of insecticides for mosquito control. We used RNA-sequencing to compare transcriptomes between permethrin resistant and susceptible strains of Florida Ae. aegypti in response to Zika virus infection. A total of 2459 transcripts were expressed at significantly different levels between resistant and susceptible Ae. aegypti. Gene ontology analysis placed these genes into seven categories of biological processes. The 863 transcripts were expressed at significantly different levels between the two mosquito strains (up/down regulated) more than 2-fold. Quantitative real-time PCR analysis was used to validate the Zika-infection response. Our results suggested a highly overexpressed P450, with AAEL014617 and AAEL006798 as potential candidates for the molecular mechanism of permethrin resistance in Ae. aegypti. Our findings indicated that most detoxification enzymes and immune system enzymes altered their gene expression between the two strains of Ae. aegypti in response to Zika virus infection. Understanding the interactions of arboviruses with resistant mosquito vectors at the molecular level allows for the possible development of new approaches in mitigating arbovirus transmission. This information sheds light on Zika-induced changes in insecticide resistant Ae. aegypti with implications for mosquito control strategies.