Identification and classification of detoxification enzymes from Culex quinquefasciatus (Diptera: Culicidae)

Potential key genes involved in metabolic resistance to malathion in the southern house mosquito, Culex quinquefasciatus, and functional validation of CYP325BC1 and CYP9M12 as candidate genes using RNA interference

BACKGROUND

Metabolic detoxification is one of the major mechanisms contributing to the development of resistance in mosquitoes, including the southern house mosquito, Culex quinquefasciatus. The three major detoxification supergene families, cytochrome P450s, glutathione S-transferases and general esterases, have been demonstrated to play an important role in metabolic resistance. In this study, we performed differential gene expression analysis based on high-throughput transcriptome sequencing on samples from four experimental groups to give insight into key genes involved in metabolic resistance to malathion in Cx. quinquefasciatus. We conducted a whole transcriptome analysis of field captured wild Cx. quinquefasciatus from Harris County (WI), Texas and a malathion susceptible laboratory-maintained Sebring colony (CO) to investigate metabolic insecticide resistance. Field captured mosquitoes were also phenotypically classified into the malathion resistant and malathion susceptible groups following a mortality response measure conducted using a Centers for Disease Control and Prevention (CDC) bottle assay. The live (MR) and dead (MS) specimens from the bottle assay, along with an unselected WI sample and a CO sample were processed for total RNA extraction and subjected to whole-transcriptome sequencing.

RESULTS

We demonstrated that the genes coding for detoxification enzymes, particularly cytochrome P450s, were highly up-regulated in the MR group compared to the MS group with similar up-regulation observed in the WI group compared to the CO group. A total of 1,438 genes were differentially expressed in comparison between MR and MS group, including 614 up-regulated genes and 824 down-regulated genes. Additionally, 1,871 genes were differentially expressed in comparison between WI and CO group, including 1,083 up-regulated genes and 788 down-regulated genes. Further analysis on differentially expressed genes from three major detoxification supergene families in both comparisons resulted in 16 detoxification genes as candidates potentially associated with metabolic resistance to malathion. Knockdown of CYP325BC1 and CYP9M12 using RNA interference on the laboratory-maintained Sebring strain significantly increased the mortality of Cx. quinquefasciatus after exposure to malathion.

CONCLUSION

We generated substantial transcriptomic evidence on metabolic detoxification of malathion in Cx. quinquefasciatus. We also validated the functional roles of two candidate P450 genes identified through DGE analysis. Our results are the first to demonstrate that knockdown of CYP325BC1 and CYP9M12 both significantly increased malathion susceptibility in Cx. quinquefasciatus, indicating involvement of these two genes in metabolic resistance to malathion.

Comparative genomic analysis of carboxylesterase genes in Tenebrio molitor and other four tenebrionids

Carboxylesterases (COEs) have various functions in wide taxons of organisms. In insects, COEs are important enzymes involved in the hydrolysis of a variety of ester-containing xenobiotics, neural signal transmission, pheromone degradation, and reproductive development. Understanding the diversity of COEs is basic to illustrate their functions. In this study, we identified 53, 105, 37, and 39 COEs from the genomes of Tenebrio molitor, Asbolus verucosus, Hycleus cichorii, and H. phaleratus in the superfamily of Tenebrionidea, respectively. Phylogenetic analysis showed that 234 COEs from these four species and those reported in Tribolium castaneum (63) could be divided into 12 clades and three major classes. The α-esterases significantly expanded in T. molitor, A. verucosus, and T. castaneum compared to dipteran and hymenopteran insects. In T. molitor, most COEs showed tissue and stage-specific but not a sex-biased expression. Our results provide insights into the diversity and evolutionary characteristics of COEs in tenebrionids, and lay a foundation for the functional characterization of COEs in the yellow mealworm.

Comparative Molecular Description of a Novel GST Gene in Culex pipiens (Diptera: Culicidae)

Repeated exposure to insecticides, particularly pyrethroids and organophosphates, has resulted in the development of insecticide resistance in the mosquito Culex pipiens, a primary disease vector. Glutathione S-transferase (GST) is involved in the phase II detoxification of numerous xenobiotics, including insecticides. In this study, a GST gene (CPIJ002678) was amplified, sequenced, and used in comprehensive molecular analyses ending up in development of a rapid assay to distinguish more tolerant individuals from susceptible Culex pipiens using the Restriction Fragment Length Polymorphism (RFLP) technique. Field collected Culex pipiens strains from untreated areas, organophosphates-treated areas and a lab strain reared for many generations, all were used in CDC bottle bioassays to evaluate the susceptibility status of the studied individuals to malathion insecticide. Interestingly, both field sites collected groups showed high levels of resistance at the malathion diagnostic time. Gene amplification, and bidirectional direct sequencing results were analyzed. Compared with the reference genome sequence, the pairwise alignment of the amplified sequences showed 96.6% similarity to the reference sequence in the GenBank database. The confirmed gene sequences were assembled and aligned using various bioinformatic softwares. The assembled contigs were used in NEBcutter V2.0 for constructing restriction maps and checked for the availability of differences (if present) between susceptible and more tolerant strains. Specific molecular RFLP markers were successfully recognized to differentiate the more tolerant from the susceptible Culex pipiens phenotypes.

Insecticide resistance, fitness and susceptibility to Zika infection of an interbred Aedes aegypti population from Rio de Janeiro, Brazil

Background

Aedes aegypti is a vector of high relevance, since it transmits several arboviruses, including dengue, chikungunya and Zika. Studies on vector biology are usually conducted with laboratory strains presenting a divergent genetic composition from field populations. This may impair vector control policies that were based on laboratory observations employing only long maintained laboratory strains. In the present study we characterized a laboratory strain interbreed with Ae. aegypti collected from five different localities in Rio de Janeiro (Aedes Rio), for insecticide resistance (IR), IR mechanisms, fitness and Zika virus infection.

Methods

We compared the recently established Aedes Rio with the laboratory reference strain Rockefeller. Insecticide resistance (deltamethrin, malathion and temephos), activity of metabolic resistance enzymes and kdr mutation frequency were determined. Some life table parameters (longevity, blood-feeding, number and egg viability) and Zika virus susceptibility was also determined.

Results

Aedes Rio showed resistance to deltamethrin (resistance ratio, RR₅₀ = 32.6) and temephos (RR₅₀ = 7.0) and elevated activity of glutathione S-transferase (GST) and esterases (α-EST and pNPA-EST), but not acetylcholinesterase (AChE). In total, 92.1% of males genotyped for kdr presented a “resistant” genotype. Weekly blood-fed females from both strains, presented reduced mortality compared to sucrose-fed mosquitoes; however, Aedes Rio blood-fed females did not live as long (mean lifespan: Rockefeller = 70 ± 3.07; Aedes Rio = 53.5 ± 2.16 days). There were no differences between strains in relation to blood-feeding and number of eggs, but Aedes Rio eggs presented reduced viability (mean hatch: Rockefeller = 77.79 ± 1.4%; Aedes Rio = 58.57 ± 1.77%). Zika virus infection (plaque-forming unit, PFU) was similar in both strains (mean PFU ± SE: Aedes Rio: 4.53 × 10⁴ ± 1.14 × 10⁴ PFU; Rockefeller: 2.02 × 10⁴ ± 0.71 × 10⁴ PFU).

Conclusion

Selected conditions in the field, such as IR mechanisms, may result in pleiotropic effects that interfere in general physiology of the insect. Therefore, it is important to well characterize field populations to be tested in parallel with laboratory reference strains. This practice would improve the significance of laboratory tests for vector control methods.

Transcriptomic analysis of insecticide resistance in the lymphatic filariasis vector Culex quinquefasciatus

Culex quinquefasciatus plays an important role in transmission of vector-borne diseases of public health importance, including lymphatic filariasis (LF), as well as many arboviral diseases. Currently, efforts to tackle C. quinquefasciatus vectored diseases are based on either mass drug administration (MDA) for LF, or insecticide-based interventions. Widespread and intensive insecticide usage has resulted in increased resistance in mosquito vectors, including C. quinquefasciatus. Herein, the transcriptome profile of Ugandan bendiocarb-resistant C. quinquefasciatus was explored to identify candidate genes associated with insecticide resistance. High levels of insecticide resistance were observed for five out of six insecticides tested, with the lowest mortality (0.97%) reported to permethrin, while for DDT, lambdacyhalothrin, bendiocarb and deltamethrin the mortality rate ranged from 1.63-3.29%. Resistance to bendiocarb in exposed mosquitoes was marked, with 2.04% mortality following 1 h exposure and 58.02% after 4 h. Genotyping of the G119S Ace-1 target site mutation detected a highly significant association (p < 0.0001; OR = 25) between resistance and Ace1-119S. However, synergist assays using the P450 inhibitor PBO, or the esterase inhibitor TPP resulted in markedly increased mortality (to ≈80%), suggesting a role of metabolic resistance in the resistance phenotype. Using a novel, custom 60 K whole-transcriptome microarray 16 genes significantly overexpressed in resistant mosquitoes were detected, with the P450 Cyp6z18 showing the highest differential gene expression (>8-fold increase vs unexposed controls). These results provide evidence that bendiocarb resistance in Ugandan C. quinquefasciatus is mediated by both target-site mechanisms and over-expression of detoxification enzymes.

Identification and classification of differentially expressed genes in pyrethroid-resistant Culex pipiens pallens

Culex pipiens pallens is an important vector that transmits Bancroftian filariasis, Japanese encephalitis and other diseases that pose a serious threat to human health. Extensive and improper use of insecticides has caused insecticide resistance in mosquitoes, which has become an important obstacle to the control of mosquito-borne diseases. It is crucial to investigate the underlying mechanism of insecticide resistance. The aims of this study were to identify genes involved in insecticide resistance based on the resistance phenotype, gene expression profile and single-nucleotide polymorphisms (SNPs) and to screen for major genes controlling insecticide resistance. Using a combination of SNP and transcriptome data, gene expression quantitative trait loci (eQTLs) were studied in deltamethrin-resistant mosquitoes. The most differentially expressed pathway in the resistant group was identified, and a regulatory network was built using these SNPs and the differentially expressed genes (DEGs) in this pathway. The major candidate genes involved in the control of insecticide resistance were analyzed by qPCR, siRNA microinjection and CDC bottle bioassays. A total of 85 DEGs that encoded putative detoxification enzymes (including 61 P450s) were identified in this pathway. The resistance regulatory network was built using SNPs, and these metabolic genes, and a major gene, CYP9AL1, were identified. The functional role of CYP9AL1 in insecticide resistance was confirmed by siRNA microinjection and CDC bottle bioassays. Using the eQTL approach, we identified important genes in pyrethroid resistance that may aid in understanding the mechanism underlying insecticide resistance and in targeting new measures for resistance monitoring and management.

Identification of carboxylesterase genes associated with pyrethroid resistance in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Carboxylesterases (CCEs) are one of three large detoxification enzyme families. Some CCEs are active on synthetic insecticides with ester structures. Anopheles sinensis is an important malaria vector in eastern Asia. This study identified and characterized the CCE genes in the A. sinensis genome and determined CCE genes associated with pyrethroid resistance using RNA sequencing (RNA-seq) and quantitative reverse transcription - polymerase chain reaction (qRT-PCR), in A. sinensis from Anhui, Chongqing, and Yunnan in China.

RESULTS

Fifty-seven putative CCEs were identified and placed into three classes, 12 subfamilies and 14 clades through phylogenetic and homology analyses. Exon sizes ranged from 31 to 4317 bp, with 49 CCEs having two to five exons and eight having six to 11 exons. A total of 183 introns were recognized with sizes ranging from 31 to 4317 bp. The 57 CCEs were located on 14 scaffolds, with 70% located on four scaffolds. The alpha-esterase subfamily was significantly expanded compared with that of Anopheles gambiae. In a pyrethroid-resistant strain, RNA-seq detected five upregulated CCE genes and qRT-PCR detected 12 upregulated CCE genes. The α-esterase 10 (AsAe10) and acetylcholinesterase 1 (AsAce1) genes were the main CCE genes associated with pyrethroid resistance.

CONCLUSION

This information will be useful for further study of the CCE gene family and pyrethroid resistance mechanisms mediated by CCEs. © 2017 Society of Chemical Industry.

Transcription factor cap n collar C regulates multiple cytochrome P450 genes conferring adaptation to potato plant allelochemicals and resistance to imidacloprid in Leptinotarsa decemlineata (Say)

Colorado potato beetle (CPB), Leptinotarsa decemlineata is a notorious pest of potato. Co-evolution with Solanaceae plants containing high levels of toxins (glycoalkaloids) helped this insect to develop an efficient detoxification system and resist almost every chemical insecticide introduced for its control. Even though the cross-resistance between plant allelochemicals and insecticides is well acknowledged, the underlying molecular mechanisms are not understood. Here, we investigated the molecular mechanisms involved in detoxification of potato plant allelochemicals and imidacloprid resistance in the field-collected CPB. Our results showed that the imidacloprid-resistant beetles employ metabolic detoxification of both potato plant allelochemicals and imidacloprid by upregulation of common cytochrome P450 genes. RNAi aided knockdown identified four cytochromes P450 genes (CYP6BJ^a/b, CYP6BJ1v1, CYP9Z25, and CYP9Z29) that are required for defense against both natural and synthetic chemicals. These P450 genes are regulated by the xenobiotic transcription factors Cap n Collar C, CncC and muscle aponeurosis fibromatosis, Maf. Studies on the CYP9Z25 promoter using the luciferase reporter assay identified two binding sites (i.e. GCAGAAT and GTACTGA) for CncC and Maf. Overall, these data showed that CPB employs the metabolic resistance mediated through xenobiotic transcription factors CncC and Maf to regulate multiple P450 genes and detoxify both imidacloprid and potato plant allelochemicals.

Rhodnius prolixus supergene families of enzymes potentially associated with insecticide resistance

Chagas disease or American trypanosomiasis, is a potentially life-threatening illness caused by the protozoan parasite, Trypanosoma cruzi. Once known as an endemic health problem of poor rural populations in Latin American countries, it has now spread worldwide. The parasite is transmitted by triatomine bugs, of which Rhodnius prolixus (Hemiptera, Reduviidae, Triatominae) is one of the vectors and a model organism. This species occurs mainly in Central and South American countries where the disease is endemic. Disease prevention focuses on vector control programs that, in general, rely intensely on insecticide use. However, the massive use of chemical insecticides can lead to resistance. One of the major mechanisms is known as metabolic resistance that is associated with an increase in the expression or activity of detoxification genes. Three of the enzyme families that are involved in this process - carboxylesterases (CCE), glutathione s-transferases (GST) and cytochrome P450s (CYP) - are analyzed in the R. prolixus genome. A similar set of detoxification genes to those of the Hemipteran Acyrthosiphon pisum but smaller than in most dipteran species was found in R. prolixus genome. All major CCE classes (43 genes found) are present but the pheromone/hormone processing class had fewer genes than usual. One main expansion was detected on the detoxification/dietary class. The phosphotriesterase family, recently associated with insecticide resistance, was also represented with one gene. One microsomal GST gene was found and the cytosolic GST gene count (14 genes) is extremely low when compared to the other hemipteran species with sequenced genomes. However, this is similar to Apis mellifera, a species known for its deficit in detoxification genes. In R. prolixus 88 CYP genes were found, with representatives in the four clans (CYP2, CYP3, CYP4 and mitochondrial) usually found in insects. R. prolixus seems to have smaller species-specific expansions of CYP genes than mosquitoes and beetles, among others. The number of R. prolixus CYP genes is similar to the hemipteran Ac. pisum, although with a bigger expansion in CYP3 and CYP4 clans, along with several gene fragments, mostly in CYP4 clan. Eleven founding members of new families were detected, consisting of ten genes in the CYP3 clan and 1 gene in the CYP4 clan. Members of these clans were proposed to have important detoxification roles in insects. The identification of CCE, GST and CYP genes is of utmost importance for directing detoxification studies on triatomines that can help insecticide management strategies in control programs.

Genomic Analysis of Detoxification Supergene Families in the Mosquito Anopheles sinensis

Anopheles sinensis is an important malaria vector in China and other Southeast Asian countries, and the emergence of insecticide resistance in this mosquito poses a serious threat to the efficacy of malaria control programs. The recently published An. sinensis genome and transcriptome provide an opportunity to understand the molecular mechanisms of insecticide resistance. Analysis of the An. sinensis genome revealed 174 detoxification genes, including 93 cytochrome P450s (P450s), 31 glutathione-S-transferases (GSTs), and 50 choline/carboxylesterases (CCEs). The gene number was similar to that in An. gambiae, but represented a decrease of 29% and 42% compared with Aedes aegypti and Culex quinquefasciatus, respectively. The considerable contraction in gene number in Anopheles mosquitoes mainly occurred in two detoxification supergene families, P450s and CCEs. The available An. sinensis transcriptome was also re-analyzed to further identify key resistance-associated detoxification genes. Among 174 detoxification genes, 124 (71%) were detected. Several candidate genes overexpressed in a deltamethrin-resistant strain (DR-strain) were identified as belonging to the CYP4 or CYP6 family of P450s and the Delta GST class. These generated data provide a basis for identifying the resistance-associated genes of An. sinensis at the molecular level.