Transcriptional analysis of insecticide resistance in Anopheles stephensi using cross-species microarray hybridization

Genome-wide annotation and comparative analysis revealed conserved cuticular protein evolution among non-biting midges with varied environmental adaptability

Chironomidae, non-biting midges, a diverse and abundant insect group in global aquatic ecosystems, represent an exceptional model for investigating genetic adaptability mechanisms in aquatic insects due to their extensive species diversity and resilience to various environmental conditions. The cuticle in insects acts as the primary defense against ecological pressures. Cuticular Proteins (CPs) determine cuticle characteristics, facilitating adaptation to diverse challenges. However, systematic annotation of CP genes has only been conducted for one Chironomidae species, Propsilocerus akamusi, by our team. In this study, we expanded this annotation by identifying CP genes in eight additional Chironomidae species, covering all Chironomidae species with available genome data. We identified a total of 889 CP genes, neatly categorized into nine CP families: 215 CPR RR1 genes, 272 CPR RR2 genes, 23 CPR RR3 genes, 21 CPF genes, 16 CPLCA genes, 19 CPLCG genes, 28 CPLCP genes, 77 CPAP genes, and 37 Tweedle genes. Subsequently, we conducted a comprehensive phylogenetic analysis of CPs within the Chironomidae family. This expanded annotation of CP genes across diverse Chironomidae species significantly contributes to our understanding of their remarkable adaptability.

Gut microbiome helps honeybee (Apis mellifera) resist the stress of toxic nectar plant (Bidens pilosa) exposure: Evidence for survival and immunity

Honeybee (Apis mellifera) ingestion of toxic nectar plants can threaten their health and survival. However, little is known about how to help honeybees mitigate the effects of toxic nectar plant poisoning. We exposed honeybees to different concentrations of Bidens pilosa flower extracts and found that B. pilosa exposure significantly reduced honeybee survival in a dose-dependent manner. By measuring changes in detoxification and antioxidant enzymes and the gut microbiome, we found that superoxide dismutase, glutathione-S-transferase and carboxylesterase activities were significantly activated with increasing concentrations of B. pilosa and that different concentrations of B. pilosa exposure changed the structure of the honeybee gut microbiome, causing a significant reduction in the abundance of Bartonella (p < 0.001) and an increase in Lactobacillus. Importantly, by using Germ-Free bees, we found that colonization by the gut microbes Bartonella apis and Apilactobacillus kunkeei (original classification as Lactobacillus kunkeei) significantly increased the resistance of honeybees to B. pilosa and significantly upregulated bee-associated immune genes. These results suggest that honeybee detoxification systems possess a level of resistance to the toxic nectar plant B. pilosa and that the gut microbes B. apis and A. kunkeei may augment resistance to B. pilosa stress by improving host immunity.

Identifying potential insecticide resistance markers through genomic-level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.

Identification of cuticular protein genes and analysis of their roles in phosphine resistance of the rusty grain beetle Cryptolestes ferrugineus

The rusty grain beetle, Cryptolestes ferrugineus (Stephens) is one of the most economically important stored grain pests, and it has evolved the high resistance to phosphine. Cuticular proteins (CPs) are the major structural components of insect cuticle, and previous studies have confirmed that CPs were involved in insecticide resistance. However, the CPs of C. ferrugineus are still poorly characterized, and thus we conducted transcriptome-wide identification of CP genes and analyze their possible relationships with phosphine resistance in this pest. In this study, a total of 122 putative CPs were annotated in the C. ferrugineus transcriptome data by blasting with the known CPs of Tribolium castaneum. The analysis of conserved motifs revealed these CPs of C. ferrugineus belonging to 9 different families, including 87 CPR, 13 CPAP1, 7 CPAP3, 3 Tweedle, 1 CPLCA, 1 CPLCG, 5 CPLCP, 2 CPCFC, and 3 CPFL proteins. The further phylogenetic analysis showed the different evolutionary patterns of CPs. Namely, we found some CPs (CPR family) formed species-specific protein clusters, indicating these CPs might occur independently among insect taxa, and while some other CPs (CPAP1 and CPAP3 family) shared a closer correlation based on the architecture of protein domains. Subsequently, the previous RNA-seq data were applied to establish the expression profiles of CPs in a phosphine susceptible and resistant populations of C. ferrugineus, and a large amount of CP genes were found to be over-expressed in resistant insects. Lastly, an up-regulated CP gene (CPR family) was selected for the further functional analysis, and after this gene was silenced via RNA interference (RNAi), the sensitivity to phosphine was significantly enhanced in C. ferrugineus. In conclusion, the present results provided us an overview of C. ferrugineus CPs, and which suggested that the CPs might play the critical roles in phosphine resistance.

Genome-wide analyses of Glutathione S-transferase gene family and expression profiling under deltamethrin exposure in non-biting midge Propsilocerus akamusi

Glutathione S-transferases (GSTs) are major enzymes in detoxification phase II, and have been functioned in resistance to various insecticides or oxidative stress. Herein, we selected the non-biting midge, Propsilocerus akamusi, widespread in Asian aquatic ecosystems, to uncover the gene location, structure, and phylogenetics relationship of GSTs at genome scale first time. Thirty-three cytosolic and four microsomal GST genes were identified and located on the four chromosomes. The cytosolic GSTs involved in the eight subclasses and five GST genes were unclassified. The expansion of GST genes in P. akamusi experienced duplication events on the delta, theta, xi, iota, and unclassified subclasses. The RNA-Seq analyses and RT-qPCR validation showed that the expression of PaGSTt2 gene is significantly elevated, with deltamethrin concentration increasing. The tertiary structure of PaGSTt2 enzyme was reconstructed, which was different from the other theta gene in the active site. In addition, the GST genes of six chironomids were first described based on the assembled genomes to explore the difference of those in the adaptation to kinds of environments. The GST frame for P. akmusi and its expression profiles provide valuable resources to understand their role in insecticide resistance of this species, as well as those of other biting midges.

Detect Cytochrome C Oxidase- and Glutathione-S-Transferase-Mediated Detoxification in a Permethrin-Resistant Population of Lygus lineolaris

Frequent sprays on cotton prompted resistance development in the tarnished plant bug (TPB). Knowledge of global gene regulation is highly desirable to better understand resistance mechanisms and develop molecular tools for monitoring and managing resistance. Novel microarray expressions of 6688 genes showed 3080 significantly up- or down-regulated genes in permethrin-treated TPBs. Among the 1543 up-regulated genes, 255 code for 39 different enzymes, and 15 of these participate in important pathways and metabolic detoxification. Oxidase is the most abundant and over-expressed enzyme. Others included dehydrogenases, synthases, reductases, and transferases. Pathway analysis revealed several oxidative phosphorylations associated with 37 oxidases and 23 reductases. One glutathione-S-transferase (GST LL_2285) participated in three pathways, including drug and xenobiotics metabolisms and pesticide detoxification. Therefore, a novel resistance mechanism of over-expressions of oxidases, along with a GST gene, was revealed in permethrin-treated TPB. Reductases, dehydrogenases, and others may also indirectly contribute to permethrin detoxification, while two common detoxification enzymes, P450 and esterase, played less role in the degradation of permethrin since none was associated with the detoxification pathway. Another potential novel finding from this study and our previous studies confirmed multiple/cross resistances in the same TPB population with a particular set of genes for different insecticide classes.

Over-expression of CP9 and CP83 increases whitefly cell cuticle thickness leading to imidacloprid resistance

Cuticular proteins (CPs) play an important role in protecting insects from adverse environmental conditions, like neonicotinoid insecticides, which are heavily used for numerous pests and caused environmental problems and public health concerns worldwide. However, the relationship between CPs and insecticides resistance in Bemisia tabaci, a serious and developed high insecticide resistance, is lacking. In this study, 125 CPs genes were identified in B. tabaci. Further phylogenetic tree showed the RR-2-type genes formed large gene groups in B. tabaci. Transcriptional expression levels of CPs genes at different developmental stages revealed that some CPs genes may play a specific role in insect development. The TEM results indicated that the cuticle thickness of susceptible strain was thinner than imidacloprid-resistance strain. Furthermore, 16 CPs genes (5 in RR-1 subfamily, 7 in RR-2 subfamily, 3 in CPAP3 subfamily and 1 in CPCFC subfamily) were activated in response to imidacloprid. And RNAi results indicated that CP9 and CP83 involved in imidacloprid resistance. In conclusion, this study was the first time to establish a basic information framework and evolutionary relationship between CPs and imidacloprid resistance in B. tabaci, which provides a basis for proposing integrated pest management strategies.

Genome-wide annotation of cuticular protein genes in non-biting midge Propsilocerus akamusi and transcriptome analysis of their response to heavy metal pollution

The insect cuticle is a sophisticated chitin-protein extracellular structure for mutable functions. The cuticles varied their structures and properties in different species, and the same species but in different regions or at different stages, to fill the requirements of different functions. The alteration of cuticle structures may also be induced due to challenges by some environmental crises, such as pollution exposures. The physical properties of the cuticle were determined by the cuticle proteins (CPs) they contain. The cuticle proteins are large protein groups in all insects, which are commonly divided into different families according to their conserved protein sequence motifs. Although Chironomidae is an abundant and universal insect in global aquatic ecosystems and a popular model for aquatic toxicology, no systematic annotation of CPs was done for any species in Chironomidae before. In this work, we annotated the CP genes of Propsilocerus akamusi, the most abundant Chironomidae species in Asia. A total of 160 CP genes were identified, and 97 of them could be well classified into eight CP families: 76 CPR genes can be subdivided into three groups (further divided into three subgroups: 36 RR1 genes, 37 RR2 genes, and 3 RR3 genes), 2 CPF genes, 3 CPLCA genes, 1 CPLCG gene, 8 CPAP genes, and 3 Tweedle genes. Additionally, we analyzed the response of P. akamusi CP genes at expression level to Cu exposure, which is related to the high heavy metal tolerance and the earlier onset of pupariation in heavy metal polluted water.

Sublethal toxicity, transgenerational effects, and transcriptome expression of the neonicotinoid pesticide cycloxaprid on demographic fitness of Coccinella septempunctata

Evaluating side effects of new neonicotinoids in terms of sublethal doses and transcriptome expression is a crucial but challenging part of integrated pest management (IPM) approaches. To this end, a study of lethal and sublethal effects on Coccinella septempunctata larvae was conducted, and an age-stage, two-sex life table procedure was performed to investigate life-table parameters. Cycloxaprid (CYC) was shown to have adverse effects on survival, development, total longevity, reproductive capacity, and predation ability in C. septempunctata. In addition, demographic growth parameters of the F1 generation such as net reproductive rate, and the intrinsic and finite rates of increase were significantly decreased under sublethal dosage LR₃₀ (1.91 g ai/hm²). These results demonstrated that the population growth of C. septempunctata was impacted by a sublethal dosage of CYC. For transcriptome expression, 544 up- and 338 down-regulated significantly differentially expressed genes (DEGs), were observed between LR₃₀ treatment and control groups. Moreover, pathways related to metabolism of retinol, carcinogenesis, biosynthesis of steroid hormone, P450 metabolism, and metabolism of xenobiotics were identified in KEGG pathway analysis. Ten DEGs were chosen and confirmed with quantitative real-time PCR analysis. Based on these findings, CYC should be considered as a component of IPM strategies in the field.

CPR63 promotes pyrethroid resistance by increasing cuticle thickness in Culex pipiens pallens

The cuticle protein (CP) encoded by CPR63 plays a role in deltamethrin resistance in Culex pipiens pallens. Herein, we investigated the distribution of CPR63 transcripts in this organism and observed high expression levels in legs and wings. Furthermore, expression of CPR63 in the legs of deltamethrin-resistant (DR) strains was 2.17-fold higher than in deltamethrin-susceptible (DS) strains. Cuticle analysis of small interfering RNA (siRNA) groups by scanning electron microscopy (SEM) revealed a significantly thinner cuticle of the tarsi in the siCPR63 group than in the siNC (negative control siRNA) group. Transmission electron microscopy (TEM) revealed that the exocuticle and endocuticle thickness of the tarsi were significantly thinner, which contributes the thinner procuticle of tarsi in the siCPR63 group than in the siNC group. Our results suggested that CPR63 might contribute to the resistance phenotype by thickening the cuticle and thereby possibly increasing the tolerance of mosquitoes to deltamethrin.

Transcriptomic analysis of resistance and short-term induction response to pyrethroids, in Anopheles coluzzii legs

BACKGROUND

Insecticide-treated bed nets and indoor residual spraying comprise the major control measures against Anopheles gambiae sl, the dominant vector in sub-Saharan Africa. The primary site of contact with insecticide is through the mosquitoes' legs, which represents the first barrier insecticides have to bypass to reach their neuronal targets. Proteomic changes and leg cuticle modifications have been associated with insecticide resistance that may reduce the rate of penetration of insecticides. Here, we performed a multiple transcriptomic analyses focusing on An. coluzzii legs.

RESULTS

Firstly, leg-specific enrichment analysis identified 359 genes including the pyrethroid-binder SAP2 and 2 other chemosensory proteins, along with 4 ABCG transporters previously shown to be leg enriched. Enrichment of gene families included those involved in detecting chemical stimuli, including gustatory and ionotropic receptors and genes implicated in hydrocarbon-synthesis. Subsequently, we compared transcript expression in the legs of a highly resistant strain (VK7-HR) to both a strain with very similar genetic background which has reverted to susceptibility after several generations without insecticide pressure (VK7-LR) and a lab susceptible population (NG). Two hundred thirty-two differentially expressed genes (73 up-regulated and 159 down-regulated) were identified in the resistant strain when compared to the two susceptible counterparts, indicating an over-expression of phase I detoxification enzymes and cuticular proteins, with decrease in hormone-related metabolic processes in legs from the insecticide resistant population. Finally, we analysed the short-term effect of pyrethroid exposure on An. coluzzii legs, comparing legs of 1 h-deltamethrin-exposed An. coluzzii (VK7-IN) to those of unexposed mosquitoes (VK7-HR) and identified 348 up-regulated genes including those encoding for GPCRs, ABC transporters, odorant-binding proteins and members of the divergent salivary gland protein family.

CONCLUSIONS

The data on An. coluzzii leg-specific transcriptome provides valuable insights into the first line of defense in pyrethroid resistant and short-term deltamethrin-exposed mosquitoes. Our results suggest that xenobiotic detoxification is likely occurring in legs, while the enrichment of sensory proteins, ABCG transporters and cuticular genes is also evident. Constitutive resistance is primarily associated with elevated levels of detoxification and cuticular genes, while short-term insecticide-induced tolerance is linked with overexpression of transporters, GPCRs and GPCR-related genes, sensory/binding and salivary gland proteins.

Genomic and Transcriptomic Analysis Reveals Cuticular Protein Genes Responding to Different Insecticides in Fall Armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a serious pest of crucial crops causing great threats to the food security of the world. It has evolved resistance to various insecticides, while the underlying molecular mechanisms remain largely unknown. Cuticular proteins (CPs), as primary components in cuticle, play an important role in insects' protection against environmental stresses. Few of them have been documented as participating in insecticide resistance in several insect species. In order to explore whether CP genes of the FAW exhibit a functional role in responding to insecticides stress, a total of 206 CPs, classified into eight families, were identified from the genome of the FAW through a homology-based approach coupled with manual efforts. The temporal expression profiles of all identified CP genes across developmental stages and their responses to 23 different insecticides were analyzed using the RNA-seq data. Expression profiling indicated that most of the CP genes displayed stage-specific expression patterns. It was found that the expression of 51 CP genes significantly changed after 48 h exposure to 17 different insecticides. The expression of eight CP genes responding to four insecticides were confirmed by RT-PCR analysis. The results showed that their overall expression profiles were consistent with RNA-seq analysis. The findings provide a basis for further functional investigation of CPs implied in insecticide stress in FAW.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

Comparative proteomics reveals mechanisms that underlie insecticide resistance in Culex pipiens pallens Coquillett

Mosquito control based on chemical insecticides is considered as an important element of the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of insecticide resistance of important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. In contrast to target site resistance, other mechanisms are far from being fully understood. Global protein profiles among cypermethrin-resistant, propoxur-resistant, dimethyl-dichloro-vinyl-phosphate-resistant and susceptible strain of Culex pipiens pallens were obtained and proteomic differences were evaluated by using isobaric tags for relative and absolute quantification labeling coupled with liquid chromatography/tandem mass spectrometric analysis. A susceptible strain of Culex pipiens pallens showed elevated resistance levels after 25 generations of insecticide selection, through iTRAQ data analysis detected 2,502 proteins, of which 1,513 were differentially expressed in insecticide-selected strains compared to the susceptible strain. Finally, midgut differential protein expression profiles were analyzed, and 62 proteins were selected for verification of differential expression using iTRAQ and parallel reaction monitoring strategy, respectively. iTRAQ profiles of adaptation selection to three insecticide strains combined with midgut profiles revealed that multiple insecticide resistance mechanisms operate simultaneously in resistant insects of Culex pipiens pallens. Significant molecular resources were developed for Culex pipiens pallens, potential candidates were involved in metabolic resistance and reducing penetration or sequestering insecticide. Future research that is targeted towards RNA interference of the identified metabolic targets, such as cuticular proteins, cytochrome P450s, glutathione S-transferases and ribosomal proteins proteins and biological pathways (drug metabolism—cytochrome P450, metabolism of xenobiotics by cytochrome P450, oxidative phosphorylation, ribosome) could lay the foundation for a better understanding of the genetic basis of insecticide resistance in Culex pipiens pallens.

Global protein profiles were compared among a susceptible strain of Cx. pipiens pallens and strains that were cypermethrin-resistant, propoxur-resistant, and dimethyl-dichloro-vinyl-phosphate-resistant after 25 generations of selection by distinct chemical insecticide families, multiple mechanisms were found to operate simultaneously in resistant mosquitoes of Cx. pipiens pallens, including mechanisms to lower penetration of or sequester the insecticide or to increase biodegradation of the insecticide via subtle alterations in either the cuticular protein levels or the activities of detoxification enzymes (P450s and glutathione S-transferases).

Anopheles metabolic proteins in malaria transmission, prevention and control: a review

The increasing resistance to currently available insecticides in the malaria vector, Anopheles mosquitoes, hampers their use as an effective vector control strategy for the prevention of malaria transmission. Therefore, there is need for new insecticides and/or alternative vector control strategies, the development of which relies on the identification of possible targets in Anopheles. Some known and promising targets for the prevention or control of malaria transmission exist among Anopheles metabolic proteins. This review aims to elucidate the current and potential contribution of Anopheles metabolic proteins to malaria transmission and control. Highlighted are the roles of metabolic proteins as insecticide targets, in blood digestion and immune response as well as their contribution to insecticide resistance and Plasmodium parasite development. Furthermore, strategies by which these metabolic proteins can be utilized for vector control are described. Inhibitors of Anopheles metabolic proteins that are designed based on target specificity can yield insecticides with no significant toxicity to non-target species. These metabolic modulators combined with each other or with synergists, sterilants, and transmission-blocking agents in a single product, can yield potent malaria intervention strategies. These combinations can provide multiple means of controlling the vector. Also, they can help to slow down the development of insecticide resistance. Moreover, some metabolic proteins can be modulated for mosquito population replacement or suppression strategies, which will significantly help to curb malaria transmission.

Effects of insecticide resistance and exposure on Plasmodium development in Anopheles mosquitoes

The spread of insecticide resistance in anopheline mosquitoes is a serious threat to the success of malaria control and prospects of elimination, but the potential impact(s) of insecticide resistance or sublethal insecticide exposure on Plasmodium-Anopheles interactions are poorly understood. Only a few studies have attempted to investigate such interactions, despite their clear epidemiological significance for malaria transmission. This short review provides an update on our understanding of the interactions between insecticide resistance and exposure and Plasmodium development, focusing on the mechanisms which might underpin any interactions, and identifying some key knowledge gaps.

Diverse patterns of constitutive and inducible overexpression of detoxifying enzyme genes among resistant Aphis glycines populations

Understanding the mechanisms of pyrethroid resistance is essential to the effective management of pesticide resistance in Aphis glycines Matsumura (Hemiptera: Aphididae). We mined putative detoxifying enzyme genes in the draft genome sequence of A. glycines for cytochrome oxidase P450 (CYP), glutathione-S-transferase (GST) and esterases (E4 and carboxylesterases-CES). Aphids from clonal populations resistant to pyrethroids from three sites in Minnesota, USA, were screened against a diagnostic LC99 concentration of either λ-cyhalothrin or bifenthrin and detoxifying enzyme genes expression in survivors was analyzed by qPCR. Their expression profiles were compared relative to a susceptible clonal population. We found 61 CYP (40 full-length), seven GST (all full-length), seven E4 (five full-length) and three CES (two full-length) genes, including 24 possible pseudogenes. The detoxifying enzymes had different expression profiles across resistant aphid populations, possibly reflecting differences in the genetic background and pyrethroid selection pressures as the number of constitutively overexpressed detoxifying enzyme genes was correlated with the level of resistance. Our findings will strengthen the understanding of the pyrethroid resistance mechanisms in A. glycines.

Larval exposure to a pyrethroid insecticide and competition for food modulate the melanisation and antibacterial responses of adult Anopheles gambiae

The insecticides we use for agriculture and for vector control often arrive in water bodies, where mosquito larvae may be exposed to them. Not only will they then likely affect the development of the larvae, but their effects may carry over to the adults, potentially affecting their capacity at transmitting infectious diseases. Such an impact may be expected to be more severe when mosquitoes are undernourished. In this study, we investigated whether exposing larvae of the mosquito Anopheles gambiae to a sub-lethal dose of permethrin (a pyrethroid) and forcing them to compete for food would affect the immune response of the adults. We found that a low dose of permethrin increased the degree to which individually reared larvae melanised a negatively charged Sephadex bead and slowed the replication of injected Escherichia coli. However, if mosquitoes had been reared in groups of three (and thus had been forced to compete for food) permethrin had less impact on the efficacy of the immune responses. Our results show how larval stressors can affect the immune response of adults, and that the outcome of exposure to insecticides strongly depends on environmental conditions.

Insecticide resistance genes affect Culex quinquefasciatus vector competence for West Nile virus

Insecticide resistance has been reported to impact the interactions between mosquitoes and the pathogens they transmit. However, the effect on vector competence for arboviruses still remained to be investigated. We examined the influence of two insecticide resistance mechanisms on vector competence of the mosquito Culex quinquefasciatus for two arboviruses, Rift Valley Fever virus (RVFV) and West Nile virus (WNV). Three Cx. quinquefasciatus lines sharing a common genetic background were used: two insecticide-resistant lines, one homozygous for amplification of the Ester² locus (SA2), the other homozygous for the acetylcholinesterase ace-1 G119S mutation (SR) and the insecticide-susceptible reference line Slab. Statistical analyses revealed no significant effect of insecticide-resistant mechanisms on vector competence for RVFV. However, both insecticide resistance mechanisms significantly influenced the outcome of WNV infections by increasing the dissemination of WNV in the mosquito body, therefore leading to an increase in transmission efficiency by resistant mosquitoes. These results showed that insecticide resistance mechanisms enhanced vector competence for WNV and may have a significant impact on transmission dynamics of arboviruses. Our findings highlight the importance of understanding the impacts of insecticide resistance on the vectorial capacity parameters to assess the overall consequence on transmission.

Insecticidal Activity and Insecticidal Mechanism of Total Saponins from Camellia oleifera

Chemical pesticides are commonly used during the cultivation of agricultural products to control pests and diseases. Excessive use of traditional pesticides can cause environmental and human health risks. There are ongoing searches for new plant-derived pesticides to reduce the use of chemical pesticides. In this study, tea saponin extracts of different purities were extracted from Camellia oleifera seeds using AB-8 macroporous resin and gradient elution with ethanol. The insecticidal effects of the tea saponin extracts were evaluated by contact toxicity tests and stomach toxicity tests using the lepidopteran pest of tea plantation, Ectropis obliqua. The total saponins extracted using 70% ethanol showed strong contact toxicity (LC₅₀ = 8.459 mg/L) and stomach toxicity (LC₅₀ = 22.395 mg/L). In-depth mechanistic studies demonstrated that tea saponins can disrupt the waxy layer of the epidermis, causing serious loss of water, and can penetrate the inside of the intestine of E. obliqua. After consumption of the tea saponins, the intestinal villi were shortened and the cavities of the intestinal wall were disrupted, which resulted in larval death. This study highlights the potential of tea saponins as a natural, plant-derived pesticide for the management of plant pests.

Genome-wide gene expression profiling reveals that cuticle alterations and P450 detoxification are associated with deltamethrin and DDT resistance in Anopheles arabiensis populations from Ethiopia

BACKGROUND

Vector control is the main intervention in malaria control and elimination strategies. However, the development of insecticide resistance is one of the major challenges for controlling malaria vectors. Anopheles arabiensis populations in Ethiopia showed resistance against both DDT and the pyrethroid deltamethrin. Although an L1014F target-site resistance mutation was present in the voltage gated sodium channel of investigated populations, the levels of resistance indicated the presence of additional resistance mechanisms. In this study, we used genome-wide transcriptome profiling by RNAseq to assess differentially expressed genes between three deltamethrin and DDT resistant An. arabiensis field populations - Asendabo, Chewaka and Tolay - and two susceptible strains - Sekoru and Mozambique.

RESULTS

Both RNAseq analysis and RT-qPCR showed that a glutathione-S-transferase, gstd3, and a cytochrome P450 monooxygenase, cyp6p4, were significantly overexpressed in the group of resistant populations compared to the susceptible strains, suggesting that the enzymes they encode play a key role in metabolic resistance against deltamethrin or DDT. Furthermore, a gene ontology enrichment analysis showed that expression changes of cuticle related genes were strongly associated with insecticide resistance. Although this did not translate in increased thickness of the procuticle, a higher cuticular hydrocarbon content was observed in a resistant population.

CONCLUSION

Our transcriptome sequencing of deltamethrin and DDT resistant An. arabiensis populations from Ethiopia suggests non-target site resistance mechanisms and paves the way for further investigation of the role of cuticle composition in insecticide resistance of malaria vectors. © 2019 Society of Chemical Industry.

Anopheles aquasalis transcriptome reveals autophagic responses to Plasmodium vivax midgut invasion

Background

Elimination of malaria depends on mastering transmission and understanding the biological basis of Plasmodium infection in the vector. The first mosquito organ to interact with the parasite is the midgut and its transcriptomic characterization during infection can reveal effective antiplasmodial responses able to limit the survival of the parasite. The vector response to Plasmodium vivax is not fully characterized, and its specificities when compared with other malaria parasites can be of fundamental interest for specific control measures.

Methods

Experimental infections were performed using a membrane-feeding device. Three groups were used: P. vivax-blood-fed, blood-fed on inactivated gametocytes, and unfed mosquitoes. Twenty-four hours after feeding, the mosquitoes were dissected and the midgut collected for transcriptomic analysis using RNAseq. Nine cDNA libraries were generated and sequenced on an Illumina HiSeq2500. Readings were checked for quality control and analysed using the Trinity platform for de novo transcriptome assembly. Transcript quantification was performed and the transcriptome was functionally annotated. Differential expression gene analysis was carried out. The role of the identified mechanisms was further explored using functional approaches.

Results

Forty-nine genes were identified as being differentially expressed with P. vivax infection: 34 were upregulated and 15 were downregulated. Half of the P. vivax-related differentially expressed genes could be related to autophagy; therefore, the effect of the known inhibitor (wortmannin) and activator (spermidine) was tested on the infection outcome. Autophagic activation significantly reduced the intensity and prevalence of infection. This was associated with transcription alterations of the autophagy regulating genes Beclin, DRAM and Apg8.

Conclusions

Our data indicate that P. vivax invasion of An. aquasalis midgut epithelium triggers an autophagic response and its activation reduces infection. This suggests a novel mechanism that mosquitoes can use to fight Plasmodium infection.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3506-8) contains supplementary material, which is available to authorized users.

Identification of detoxification genes in imidacloprid-resistant Asian citrus psyllid (Hemiptera: Lividae) and their expression patterns under stress of eight insecticides

BACKGROUND

The Asian citrus psyllid, Diaphorina citri, is one of the major pests in citrus-growing areas around the world. The application of insecticides is the most effective method to reduce the population of D. citri. However, D. citri has developed resistance to multiple classes of insecticides. Understanding resistance mechanisms is crucial to the management of D. citri. In this study, molecular assays were performed to characterize imidacloprid resistance mechanisms.

RESULTS

Based on the D. citri transcriptome database and other known insect resistance genes, 16 cytochrome P450, eight glutathione-S-transferase and six esterase genes were selected for cloning and sequencing. The gene expression analysis of 30 detoxification genes demonstrated that the relative expression of CYP4g15, CYP303A1, CYP4C62, CYP6BD5, GSTS1 and EST-6 were moderately high (>5-fold increase) in the imidacloprid-resistant strain. Feeding of double-stranded RNA (dsRNA) reduced the expression of the six genes (46.7%-72.1%) and resulted in significant adult mortality (65.62%-82.76%). We also determined the ability of different insecticides to induce the six selected genes. The expression of CYP4C62 and GSTS1 genes were the most significantly upregulated in adults treated with all insecticides, except for chlorfenapyr. In chlorfenapyr-treated D. citri, expression of CYP4g15 and CYP303A1 were the most highly induced.

CONCLUSION

Overexpressed detoxification genes were associated with imidacloprid resistance, as confirmed by RNA interference feeding tests. The induction of the six selected genes when exposed to different insecticides supported the hypothesis that they were involved in the metabolism of the tested insecticides. © 2018 Society of Chemical Industry.

Targeting detoxification genes by phloem-mediated RNAi: A new approach for controlling phloem-feeding insect pests

Many phloem-feeding insects are considered severe pests of agriculture and are controlled mainly by chemical insecticides. Continued extensive use of these inputs is environmentally undesirable, and also leads to the development of insecticide resistance. Here, we used a plant-mediated RNA interference (RNAi) approach, to develop a new control strategy for phloem-feeding insects. The approach aims to silence "key" detoxification genes, involved in the insect's ability to neutralize defensive and toxic plant chemistry. We targeted a glutathione S-transferase (GST) gene, BtGSTs5, in the phloem-feeding whitefly Bemisia tabaci, a devastating global agricultural pest. We report three major findings. First, significant down regulation of the BtGSTs5 gene was obtained in the gut of B. tabaci when the insects were fed on Arabidopsis thaliana transgenic plants expressing dsRNA against BtGSTs5 under a phloem-specific promoter. This brings evidence that phloem-feeding insects can be efficiently targeted by plant-mediated RNAi. Second, in-silico and in-vitro analyses indicated that the BtGSTs5 enzyme can accept as substrates, hydrolyzed aliphatic- and indolic-glucosinolates, and produce their corresponding detoxified conjugates. Third, performance assays suggested that the BtGSTs5 gene silencing prolongs the developmental period of B. tabaci nymphs. Taken together, these findings suggest that BtGSTs5 is likely to play an important role in enabling B. tabaci to successfully feed on glucosinolate-producing plants. Targeting the gene by RNAi in Brassicaceae cropping systems, will likely not eliminate the pest populations from the fields but will significantly reduce their success over the growing season, support prominent activity of natural enemies, eventually allowing the establishment of stable and sustainable agroecosystem.

Insect cuticle: a critical determinant of insecticide resistance

Intense use of insecticides has resulted in the selection of extreme levels of resistance in insect populations. Therefore understanding the molecular basis of insecticide resistance mechanisms becomes critical. Penetration resistance refers to modifications in the cuticle that will eventually slow down the penetration of insecticide molecules within insects' body. So far, two mechanisms of penetration resistance have been described, the cuticle thickening and the altering of cuticle composition. Cuticular modifications are attributed to the over-expression of diversified genes or proteins, which belong to structural components (cuticular proteins mainly), enzymes that catalyze enzymatic reactions (CYP4G16 and laccase 2) or ABC transporters that promote cuticular translocation. In the present review we summarize recent studies and discuss future perspectives.

Genome-wide annotation of cuticular proteins in the oriental fruit fly (Bactrocera dorsalis), changes during pupariation and expression analysis of CPAP3 protein genes in response to environmental stresses

Cuticular proteins (CPs) are essential components of the insect cuticle as they create a structural and protective shield and may have a role in insect development. In this paper, we studied the CPs in the oriental fruit fly (Bactrocera dorsalis), one of the most economically important pests in the Tephritidae family around the world. The availability of a complete genome sequence (NCBI Assembly: ASM78921v2) allowed the identification of 164 CP genes in B. dorsalis. Comparative analysis of the CPs in B. dorsalis with those in the model insect Drosophila melanogaster and the closely related Ceratitis capitata, and CPs from mosquitoes, Lepidoptera, Hymenoptera and Coleoptera identified Diptera-specific genes and cuticle development patterns. Analysis of their evolutionary relationship revealed that some CP families had evolved according to the phylogeny of the different insect species, while others shared a closer relationship based on domain architecture. Subsequently, transcriptome analysis showed that while most of the CPs (60-100% of the family members) are expressed in the epidermis, some were also present in internal organs such as the fat body and the reproductive organs. Furthermore, the study of the expression profiles throughout development revealed a profound change in the expression of CPs during the formation of the puparium (pupariation). Further analysis of the expression profiles of the CPAP3 genes under various environmental stresses revealed them to be involved in the response to pesticides and arid and extreme temperatures conditions. In conclusion, the data provide a particular overview of CPs and their evolutionary and transcriptional dynamics, and in turn they lay a molecular foundation to explore their roles in the unique developmental process of insect metamorphosis and stress responses.

3,4-Dihydroxyphenylacetaldehyde synthase and cuticle formation in insects

Cuticle is the most important structure that protects mosquitoes and other insect species from adverse environmental conditions and infections of microorganism. The physiology and biochemistry of insect cuticle formation have been studied for many years and our understanding of cuticle formation and hardening has increased considerably. This is especially true for flexible cuticle. The recent discovery of a novel enzyme that catalyzes the production of 3,4-dihydroxyphenylacetaldehyde (DOPAL) in insects provides intriguing insights concerning the flexible cuticle formation in insects. For convenience, the enzyme that catalyzes the production DOPAL from l-dopa is named DOPAL synthase. In this mini-review, we summarize the biochemical pathways of cuticle formation and hardening in general and discuss DOPAL synthase-mediated protein crosslinking in insect flexible cuticle in particular.

Gut Microbiota Mediate Insecticide Resistance in the Diamondback Moth, Plutella xylostella (L.)

The development of insecticide resistance in insect pests is a worldwide concern and elucidating the underlying mechanisms is critical for effective crop protection. Recent studies have indicated potential links between insect gut microbiota and insecticide resistance and these may apply to the diamondback moth, Plutella xylostella (L.), a globally and economically important pest of cruciferous crops. We isolated Enterococcus sp. (Firmicutes), Enterobacter sp. (Proteobacteria), and Serratia sp. (Proteobacteria) from the guts of P. xylostella and analyzed the effects on, and underlying mechanisms of insecticide resistance. Enterococcus sp. enhanced resistance to the widely used insecticide, chlorpyrifos, in P. xylostella, while in contrast, Serratia sp. decreased resistance and Enterobacter sp. and all strains of heat-killed bacteria had no effect. Importantly, the direct degradation of chlorpyrifos in vitro was consistent among the three strains of bacteria. We found that Enterococcus sp., vitamin C, and acetylsalicylic acid enhanced insecticide resistance in P. xylostella and had similar effects on expression of P. xylostella antimicrobial peptides. Expression of cecropin was down-regulated by the two compounds, while gloverin was up-regulated. Bacteria that were not associated with insecticide resistance induced contrasting gene expression profiles to Enterococcus sp. and the compounds. Our studies confirmed that gut bacteria play an important role in P. xylostella insecticide resistance, but the main mechanism is not direct detoxification of insecticides by gut bacteria. We also suggest that the influence of gut bacteria on insecticide resistance may depend on effects on the immune system. Our work advances understanding of the evolution of insecticide resistance in this key pest and highlights directions for research into insecticide resistance in other insect pest species.

ASGDB: a specialised genomic resource for interpreting Anopheles sinensis insecticide resistance

BACKGROUND

Anopheles sinensis is an important malaria vector in Southeast Asia. The widespread emergence of insecticide resistance in this mosquito species poses a serious threat to the efficacy of malaria control measures, particularly in China. Recently, the whole-genome sequencing and de novo assembly of An. sinensis (China strain) has been finished. A series of insecticide-resistant studies in An. sinensis have also been reported. There is a growing need to integrate these valuable data to provide a comprehensive database for further studies on insecticide-resistant management of An. sinensis.

RESULTS

A bioinformatics database named An. sinensis genome database (ASGDB) was built. In addition to being a searchable database of published An. sinensis genome sequences and annotation, ASGDB provides in-depth analytical platforms for further understanding of the genomic and genetic data, including visualization of genomic data, orthologous relationship analysis, GO analysis, pathway analysis, expression analysis and resistance-related gene analysis. Moreover, ASGDB provides a panoramic view of insecticide resistance studies in An. sinensis in China. In total, 551 insecticide-resistant phenotypic and genotypic reports on An. sinensis distributed in Chinese malaria-endemic areas since the mid-1980s have been collected, manually edited in the same format and integrated into OpenLayers map-based interface, which allows the international community to assess and exploit the high volume of scattered data much easier. The database has been given the URL: http://www.asgdb.org /.

CONCLUSIONS

ASGDB was built to help users mine data from the genome sequence of An. sinensis easily and effectively, especially with its advantages in insecticide resistance surveillance and control.

The Evolution and Metamorphosis of Arthropod Proteomics and Genomics

This article presents an overview of the development of techniques for analyzing cuticular proteins (CPs), their transcripts, and their genes over the past 50 years based primarily on experience in the laboratory of J.H. Willis. It emphasizes changes in the kind of data that can be gathered and how such data provided insights into the molecular underpinnings of insect metamorphosis and cuticle structure. It describes the techniques that allowed visualization of the location of CPs at both the anatomical and intracuticular levels and measurement of the appearance and deployment of transcripts from CP genes as well as what was learned from genomic and transcriptomic data. Most of the early work was done with the cecropia silkmoth, Hyalophora cecropia, and later work was with Anopheles gambiae.

Culex pipiens pallens cuticular protein CPLCG5 participates in pyrethroid resistance by forming a rigid matrix

Background

Chemical insecticides have hugely reduced the prevalence of vector-borne diseases around the world, but resistance threatens their continued effectiveness. Despite its importance, cuticle resistance is an under-studied area, and exploring the detailed molecular basis of resistance is critical for implementing suitable resistance management strategies.

Methods

We performed western blotting of cuticular protein CPLCG5 in deltamethrin-susceptible (DS) and laboratory-produced deltamethrin-resistant (DR) strains of Culex pipiens pallens. Immunofluorescence assays using a polyclonal antibody to locate cuticular CPLCG5 in mosquitoes. EM immunohistochemical analysis of the femur segment was used to compare the cuticle in control and CPLCG5-deficient siRNA experimental groups.

Results

The gene CPLCG5 encodes a cuticle protein that plays an important role in pyrethroid resistance. Based on a prior study, we found that expression of CPLCG5 was higher in the resistant (DR) strain than the susceptible (DS) strain. CPLCG5 transcripts were abundant in white pupae and 1-day-old adults, but expression was dramatically decreased in 3-day-old adults, then remained stable thereafter. Western blotting revealed that the CPLCG5 protein was ~2.2-fold higher in the legs of the DR strain than the DS strain. Immunofluorescence assays revealed CPLCG5 expression in the head, thorax, abdomen, wing, and leg, and expression most abundant in the leg and wing. EM immunohistochemical analysis suggested that the exocuticle thickness of the femur was significantly thinner in the CPLCG5-deficient siCPLCG5 strain (0.717 ± 0.110 μm) than the siNC strain (0.946 ± 0.126 μm). Depletion of CPLCG5 by RNA interference resulted in unorganised laminae and a thinner cuticle.

Conclusions

The results suggest CPLCG5 participates in pyrethroid resistance by forming a rigid matrix and increasing the thickness of the cuticle.

Electronic supplementary material

The online version of this article (doi: 10.1186/s13071-017-2567-9) contains supplementary material, which is available to authorized users.

Contributions of cuticle permeability and enzyme detoxification to pyrethroid resistance in the major malaria vector Anopheles gambiae

To tackle the problem of insecticide resistance, all resistance mechanisms need to be studied. This study investigated the involvement of the cuticle in pyrethroid resistance in a strain of Anopheles gambiae, MRS, free of kdr mutations. Bioassays revealed MRS to be resistant to pyrethroids and DDT, indicated by increasing knockdown times and resistance ratios. Moreover, biochemical analysis indicated that metabolic resistance based on enhanced CYP450 activity may also play a role. Insecticide penetration assays showed that there were significantly lower amounts of insecticide in the MRS strain than in the susceptible control. Analysis of the levels of the selected transcripts by qPCR showed that CYP6M2, a major pyrethroid metaboliser, CYP4G16, a gene implicated in resistance via its contribution to the biosynthesis of elevated epicuticular hydrocarbons that delay insecticide uptake, and the cuticle genes CPAP3-E and CPLCX1 were upregulated after insecticide exposure. Other metabolic (CYP6P3, GSTe2) and cuticle (CPLCG3, CPRs) genes were also constitutively upregulated. Microscopic analysis showed that the cuticle layers of the MRS strain were significantly thicker than those of the susceptible strain. This study allowed us to assess the contribution made by the cuticle and metabolic mechanisms to pyrethroid resistance in Anopheles gambiae without target-site mutations.

Rapid evolution of female-biased genes among four species of Anopheles malaria mosquitoes

Understanding how phenotypic differences between males and females arise from the sex-biased expression of nearly identical genomes can reveal important insights into the biology and evolution of a species. Among Anopheles mosquito species, these phenotypic differences include vectorial capacity, as it is only females that blood feed and thus transmit human malaria. Here, we use RNA-seq data from multiple tissues of four vector species spanning the Anopheles phylogeny to explore the genomic and evolutionary properties of sex-biased genes. We find that, in these mosquitoes, in contrast to what has been found in many other organisms, female-biased genes are more rapidly evolving in sequence, expression, and genic turnover than male-biased genes. Our results suggest that this atypical pattern may be due to the combination of sex-specific life history challenges encountered by females, such as blood feeding. Furthermore, female propensity to mate only once in nature in male swarms likely diminishes sexual selection of post-reproductive traits related to sperm competition among males. We also develop a comparative framework to systematically explore tissue- and sex-specific splicing to document its conservation throughout the genus and identify a set of candidate genes for future functional analyses of sex-specific isoform usage. Finally, our data reveal that the deficit of male-biased genes on the X Chromosomes in Anopheles is a conserved feature in this genus and can be directly attributed to chromosome-wide transcriptional regulation that de-masculinizes the X in male reproductive tissues.

piRNA-3312: A Putative Role for Pyrethroid Resistance in Culex pipiens pallens (Diptera: Culicidae)

Piwi-interacting RNAs (piRNAs) are a newly identified class of small noncoding RNAs. They are associated with chromatin organization, messenger RNA stability, and genome structure. Although the overexpression of piRNA-3312 in deltamethrin-susceptible (DS) strain of Culex pipiens pallens (L.) was observed in our previous large-scale transcriptome data, the roles of piRNA in insecticide resistance have not been clearly defined. The aim of the present study was to investigate how piRNA-3312 is involved in insecticide resistance. The lower expression level of piRNA-3312 in deltamethrin-resistant (DR) strain of Cx. pipiens pallens was confirmed by quantitative real time polymerase chain reaction (qRT-PCR). Overexpression of piRNA-3312 in the DR strain made the mosquitoes more sensitive to deltamethrin, whereas inhibiting the expression of piRNA-3312 in the DS strain made the mosquitoes more resistant to deltamethrin. Piwi-interacting RNA-3312 was also found to bind 3' UTR (Untranslated Regions) of gut esterase 1 gene and could induce its degradation. In addition, knockdown of gut esterase 1 gene increased the sensitivity of DR strain to deltamethrin. In conclusion, we found that piRNA-3312 targeted the gut esterase 1 gene to negatively regulate the insecticide resistance. This finding facilitates the understanding of various functions of piRNAs and their association with insecticide resistance.

Insecticide resistance is mediated by multiple mechanisms in recently introduced Aedes aegypti from Madeira Island (Portugal)

BACKGROUND

Aedes aegypti is a major mosquito vector of arboviruses, including dengue, chikungunya and Zika. In 2005, Ae. aegypti was identified for the first time in Madeira Island. Despite an initial insecticide-based vector control program, the species expanded throughout the Southern coast of the island, suggesting the presence of insecticide resistance. Here, we characterized the insecticide resistance status and the underlying mechanisms of two populations of Ae. aegypti from Madeira Island, Funchal and Paúl do Mar.

METHODOLOGY/PRINCIPAL FINDINGS

WHO susceptibility bioassays indicated resistance to cyfluthrin, permethrin, fenitrothion and bendiocarb. Use of synergists significantly increased mortality rates, and biochemical assays indicated elevated activities of detoxification enzymes, suggesting the importance of metabolic resistance. Microarray-based transcriptome analysis detected significant upregulation in both populations of nine cytochrome P450 oxidase genes (including four known pyrethroid metabolizing enzymes), the organophosphate metabolizer CCEae3a, Glutathione-S-transferases, and multiple putative cuticle proteins. Genotyping of knockdown resistance loci linked to pyrethroid resistance revealed fixation of the 1534C mutation, and presence with moderate frequencies of the V1016I mutation in each population.

CONCLUSIONS/SIGNIFICANCE

Significant resistance to three major insecticide classes (pyrethroid, carbamate and organophosphate) is present in Ae. aegypti from Madeira Island, and appears to be mediated by multiple mechanisms. Implementation of appropriate resistance management strategies including rotation of insecticides with alternative modes of action, and methods other than chemical-based vector control are strongly advised to delay or reverse the spread of resistance and achieve efficient control.

Gene expression modulation of ABC transporter genes in response to permethrin in adults of the mosquito malaria vector Anopheles stephensi

Living organisms have evolved an array of genes coding for detoxifying enzymes and efflux protein pumps, to cope with endogenous and xenobiotic toxic compounds. The study of the genes activated during toxic exposure is relevant to the area of arthropod vector control, since these genes are one of the targets upon which natural selection acts for the evolution of insecticide resistance. ATP-binding cassette (ABC) transporters participate to insecticide detoxification acting as efflux pumps, that reduce the intracellular concentration of toxic compounds, or of their metabolic derivatives. Here we analyzed the modulation of the expression of six genes coding for ABC transporters, after the exposure of adult females and males of the mosquito Anopheles stephensi, a major malaria vector in Asia, to permethrin. Male and female mosquitoes were exposed to insecticide for one hour, then the expression profiles of the ABC transporter genes AnstABCB2, AnstABCB3, AnstABCB4, AnstABCBmember6, AnstABCC11, and AnstABCG4 were analysed after one and 24h. Our results showed that three genes (AnstABCB2, AnstABCBmember6, AnstABCG4) were up-regulated in both sexes; two of these (AnstABCBmember6 and AnstABCG4) have previously been shown to be up-regulated also in larval stages of An. stephensi, supporting a role for these genes in permethrin defence in larvae as well as in adults. Finally, the same ABC transporter genes were activated both in females and males; however, the timing of gene induction was different, with a prompter induction in females than in males.

The choreography of the chemical defensome response to insecticide stress: insights into the Anopheles stephensi transcriptome using RNA-Seq

Animals respond to chemical stress with an array of gene families and pathways termed “chemical defensome”. In arthropods, despite many defensome genes have been detected, how their activation is arranged during toxic exposure remains poorly understood. Here, we sequenced the transcriptome of Anopheles stephensi larvae exposed for six, 24 and 48 hours to the LD₅₀ dose of the insecticide permethrin to monitor transcriptional changes of defensome genes across time. A total of 177 genes involved in insecticide defense were differentially expressed (DE) in at least one time-point, including genes encoding for Phase 0, I, II, III and antioxidant enzymes and for Heat Shock and Cuticular Proteins. Three major patterns emerged throughout time. First, most of DE genes were down-regulated at all time-points, suggesting a reallocation of energetic resources during insecticide stress. Second, single genes and clusters of genes turn off and on from six to 48 hours of treatment, showing a modulated response across time. Third, the number of up-regulated genes peaked at six hours and then decreased during exposure. Our results give a first picture of how defensome gene families respond against toxicants and provide a valuable resource for understanding how defensome genes work together during insecticide stress.

Cytochrome P450 associated with insecticide resistance catalyzes cuticular hydrocarbon production in Anopheles gambiae

The role of cuticle changes in insecticide resistance in the major malaria vector Anopheles gambiae was assessed. The rate of internalization of (14)C deltamethrin was significantly slower in a resistant strain than in a susceptible strain. Topical application of an acetone insecticide formulation to circumvent lipid-based uptake barriers decreased the resistance ratio by ∼50%. Cuticle analysis by electron microscopy and characterization of lipid extracts indicated that resistant mosquitoes had a thicker epicuticular layer and a significant increase in cuticular hydrocarbon (CHC) content (∼29%). However, the CHC profile and relative distribution were similar in resistant and susceptible insects. The cellular localization and in vitro activity of two P450 enzymes, CYP4G16 and CYP4G17, whose genes are frequently overexpressed in resistant Anopheles mosquitoes, were analyzed. These enzymes are potential orthologs of the CYP4G1/2 enzymes that catalyze the final step of CHC biosynthesis in Drosophila and Musca domestica, respectively. Immunostaining indicated that both CYP4G16 and CYP4G17 are highly abundant in oenocytes, the insect cell type thought to secrete hydrocarbons. However, an intriguing difference was indicated; CYP4G17 occurs throughout the cell, as expected for a microsomal P450, but CYP4G16 localizes to the periphery of the cell and lies on the cytoplasmic side of the cell membrane, a unique position for a P450 enzyme. CYP4G16 and CYP4G17 were functionally expressed in insect cells. CYP4G16 produced hydrocarbons from a C18 aldehyde substrate and thus has bona fide decarbonylase activity similar to that of dmCYP4G1/2. The data support the hypothesis that the coevolution of multiple mechanisms, including cuticular barriers, has occurred in highly pyrethroid-resistant An gambiae.

Delayed mortality effects cut the malaria transmission potential of insecticide-resistant mosquitoes

Malaria transmission has been substantially reduced across Africa through the distribution of long-lasting insecticidal nets (LLINs). However, the emergence of insecticide resistance within mosquito vectors risks jeopardizing the future efficacy of this control strategy. The severity of this threat is uncertain because the consequences of resistance for mosquito fitness are poorly understood: while resistant mosquitoes are no longer immediately killed upon contact with LLINs, their transmission potential may be curtailed because of longer-term fitness costs that persist beyond the first 24 h after exposure. Here, we used a Bayesian state-space model to quantify the immediate (within 24 h of exposure) and delayed (>24 h after exposure) impact of insecticides on daily survival and malaria transmission potential of moderately and highly resistant laboratory populations of the major African malaria vector Anopheles gambiae Contact with LLINs reduced the immediate survival of moderately and highly resistant An. gambiae strains by 60-100% and 3-61%, respectively, and delayed mortality impacts occurring beyond the first 24 h after exposure further reduced their overall life spans by nearly one-half. In total, insecticide exposure was predicted to reduce the lifetime malaria transmission potential of insecticide-resistant vectors by two-thirds, with delayed effects accounting for at least one-half of this reduction. The existence of substantial, previously unreported, delayed mortality effects within highly resistant malaria vectors following exposure to insecticides does not diminish the threat of growing resistance, but posits an explanation for the apparent paradox of continued LLIN effectiveness in the presence of high insecticide resistance.

Constitutive overexpression of cytochrome P450 associated with imidacloprid resistance in Laodelphax striatellus (Fallén)

BACKGROUND

Imidacloprid is a principal insecticide for controlling rice planthoppers worldwide. Resistance to imidacloprid has been reported in a field population of Laodelphax striatellus. The present work was conducted to study the molecular mechanisms of imidacloprid resistance.

RESULTS

An imidacloprid-resistant strain was produced by selecting a field population with imidacloprid for 24 generations. Piperonyl butoxide (PBO) showed a 1.70-fold synergistic effect. Enzyme activity assays were conducted, and cytochrome P450 monooxygenase showed 1.88-fold activity. The mRNA expression levels of 57 P450 genes were compared. Four CYP genes were found to be overexpressed and significantly different to the susceptible strain. Four strains were selected with imidacloprid for a short period, and the expression levels of ten identified detoxification genes were then compared. Only CYP353D1v2 overexpressed and was significantly different to the susceptible strain. Strong correlation was found between CYP353D1v2 expression levels and imidacloprid treatments. Additionally, gene-silencing RNAi via dsRNA feeding showed that depressing the expression of CYP353D1v2 could significantly enhance the sensitivity of L. striatellus to imidacloprid.

CONCLUSION

Constitutive overexpression of four CYP genes was associated with imidacloprid resistance in long-term selection, and expression of CYP353D1v2 with imidacloprid resistance in short-term selection in L. striatellus.

Transcriptome mining: Multigene panel to test delousing drug response in the sea louse Caligus rogercresseyi

Controlling infestations of copepodid ectoparasites in the salmon industry is increasingly problematic given higher instances of drug resistance or loss of sensitivity. Despite the importance of this issue, the molecular mechanisms and genes implicated in resistance/susceptibility are only scarcely understood. The objective of the present study was to identify and evaluate the expression levels of candidate genes associated with delousing drug response in the sea louse Caligus rogercresseyi. From RNA-seq data obtained for adult male and female sea lice, 62.48 M reads were assembled in 70,349 high-quality contigs. BLASTX analysis against UniprotKB/Swiss-Prot and the ESTs available for crustaceans in the NCBI database identified 870 transcripts previously related to genes associated with delousing drug response. Furthermore, 14 candidate genes were validated through RT-qPCR and were evaluated with deltamethrin and azamethiphos bioassays. The results evidenced an overregulation of genes involved in ion transport in salmon lice treated with deltamethrin, while those treated with azamethiphos evidenced an overregulation of genes such as cytochrome P450, Carboxylesterase, and acetylcholine receptors. The present study provides a multigene panel to test delousing drug response to pyrethroids and organophosphates in a highly prevalent pathogen of the Chilean salmon industry.

The cuticle proteins: a putative role for deltamethrin resistance in Culex pipiens pallens

Insecticide resistance has been a major public health challenge. It is impendent to study the mechanism on insecticide resistance. In our previous study, 14 differentially accumulated insect cuticle proteins (ICPs) based on insecticide resistance proteomes and transcriptomes were found in the deltamethrin-resistant (DR) and -susceptible (DS) strains of Culex pipiens pallens. To investigate if these ICPs are associated with deltamethrin resistance, different transcriptional levels of the 14 ICPs were detected in the DS and DR strains from laboratory and field populations by using quantitative real-time polymerase chain reaction (qRT-PCR). The expression levels of the 14 ICPs were also measured after short-term exposure of the DS strain to deltamethrin. The full-length complementary DNA (cDNA) of CpCPLCG5 gene, which encodes one of the 14 ICPs, was cloned from Cx. pipiens pallens. Homology analysis and phylogenetic analysis were carried out with some other insects. Furthermore, small interfering RNA (siRNA) was used to knockdown the expression level of CpCPLCG5 gene for characterizing its contribution to deltamethrin resistance. The results showed that the expression level of CpCPLCG5 gene was higher in DR strain than in DS strain both in laboratory and field populations while the other 13 ICPs were downregulated. The full-length cDNA of CpCPLCG5 gene was 732 bp, with the ORF of 390 bp and deduced 129 amino acids (GenBank/KF723314,2013). Knockdown of CpCPLCG5 gene increased the susceptibility of the DR strain while the expression level of the other 13 ICPs elevated. Our findings indicate that the cuticle proteins are associated with deltamethrin resistance in Cx. pipiens pallens.

Identification of proteins associated with pyrethroid resistance by iTRAQ-based quantitative proteomic analysis in Culex pipiens pallens

Background

Mosquito control based on chemical insecticides is considered as an important element in the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of pyrethroid resistance in important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. To date, the mechanisms of pyrethroid resistance are still unclear. Recent advances in proteomic techniques can facilitate to identify pyrethroid resistance-associated proteins at a large-scale for improving our understanding of resistance mechanisms, and more importantly, for seeking some genetic markers used for monitoring and predicting the development of resistance.

Methods

We performed a quantitative proteomic analysis between a deltamethrin-susceptible strain and a deltamethrin-resistant strain of laboratory population of Culex pipiens pallens using isobaric tags for relative and absolute quantitation (iTRAQ) analysis. Gene Ontology (GO) analysis was used to find the relative processes that these differentially expressed proteins were involved in. One differentially expressed protein was chosen to confirm by Western blot in the laboratory and field populations of Cx. pipiens pallens.

Results

We identified 30 differentially expressed proteins assigned into 10 different categories, including oxidoreductase activity, transporter activity, catalytic activity, structural constituent of cuticle and hypothetical proteins. GO analysis revealed that 25 proteins were sub-categorized into 35 hierarchically-structured GO classifications. Western blot results showed that CYP6AA9 as one of the up-regulated proteins was confirmed to be overexpressed in the deltamethrin-resistant strains compared with the deltamethrin-susceptible strains both in the laboratory and field populations.

Conclusions

This is the first study to use modern proteomic tools for identifying pyrethroid resistance-related proteins in Cx. pipiens. The present study brought to light many proteins that were not previously thought to be associated with pyrethroid resistance, which further expands our understanding of pyrethroid resistance mechanisms. CYP6AA9 was overexpressed in the deltamethrin-resistant strains, indicating that CYP6AA9 may be involved in pyrethroid resistance and may be used as a potential genetic marker to monitor and predict the pyrethroid resistance level of field populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0709-5) contains supplementary material, which is available to authorized users.

Mosquito genomics. Highly evolvable malaria vectors: the genomes of 16 Anopheles mosquitoes

Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. To investigate the genomic basis of vectorial capacity and explore new avenues for vector control, we sequenced the genomes of 16 anopheline mosquito species from diverse locations spanning ~100 million years of evolution. Comparative analyses show faster rates of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses, relative to Drosophila. Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts.

Comparative microarray analyses of adult female midgut tissues from feeding Rhipicephalus species

The cattle tick, Rhipicephalus microplus, has a debilitating effect on the livestock industry worldwide, owing to its being a vector of the causative agents of bovine babesiosis and anaplasmosis. In South Africa, co-infestation with R. microplus and R. decoloratus, a common vector species on local livestock, occurs widely in the northern and eastern parts of the country. An alternative to chemical control methods is sought in the form of a tick vaccine to control these tick species. However, sequence information and transcriptional data for R. decoloratus is currently lacking. Therefore, this study aimed at identifying genes that are shared between midgut tissues of feeding adult female R. microplus and R. decoloratus ticks. In this regard, a custom oligonucleotide microarray comprising of 13,477 R. microplus sequences was used for transcriptional profiling and 2476 genes were found to be shared between these Rhipicephalus species. In addition, 136 transcripts were found to be more abundantly expressed in R. decoloratus and 1084 in R. microplus. Chi-square analysis revealed that genes involved in lipid transport and metabolism are significantly overrepresented in R. microplus and R. decoloratus. This study is the first transcriptional profiling of R. decoloratus and is an additional resource that can be evaluated further in future studies for possible tick control.

Impact of agriculture on the selection of insecticide resistance in the malaria vector Anopheles gambiae: a multigenerational study in controlled conditions

Background

Resistance of mosquitoes to insecticides is mainly attributed to their adaptation to vector control interventions. Although pesticides used in agriculture have been frequently mentioned as an additional force driving the selection of resistance, only a few studies were dedicated to validate this hypothesis and characterise the underlying mechanisms. While insecticide resistance is rising dramatically in Africa, deciphering how agriculture affects resistance is crucial for improving resistance management strategies. In this context, the multigenerational effect of agricultural pollutants on the selection of insecticide resistance was examined in Anopheles gambiae.

Methods

An urban Tanzanian An. gambiae population displaying a low resistance level was used as a parental strain for a selection experiment across 20 generations. At each generation larvae were selected with a mixture containing pesticides and herbicides classically used in agriculture in Africa. The resistance levels of adults to deltamethrin, DDT and bendiocarb were compared between the selected and non-selected strains across the selection process together with the frequency of kdr mutations. A microarray approach was used for pinpointing transcription level variations selected by the agricultural pesticide mixture at the adult stage.

Results

A gradual increase of adult resistance to all insecticides was observed across the selection process. The frequency of the L1014S kdr mutation rose from 1.6% to 12.5% after 20 generations of selection. Microarray analysis identified 90 transcripts over-transcribed in the selected strain as compared to the parental and the non-selected strains. Genes encoding cuticle proteins, detoxification enzymes, proteins linked to neurotransmitter activity and transcription regulators were mainly affected. RT-qPCR transcription profiling of candidate genes across multiple generations supported their link with insecticide resistance.

Conclusions

This study confirms the potency of agriculture in selecting for insecticide resistance in malaria vectors. We demonstrated that the recurrent exposure of larvae to agricultural pollutants can select for resistance mechanisms to vector control insecticides at the adult stage. Our data suggest that in addition to selected target-site resistance mutations, agricultural pollutants may also favor cuticle, metabolic and synaptic transmission-based resistance mechanisms. These results emphasize the need for integrated resistance management strategies taking into account agriculture activities.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0480-z) contains supplementary material, which is available to authorized users.

Biochemical and molecular characterisation and cross-resistance in field and laboratory chlorpyrifos-resistant strains of Laodelphax striatellus (Hemiptera: Delphacidae) from eastern China

BACKGROUND

Laboratory selection is often employed in resistance mechanism studies because field-derived populations commonly do not have high enough resistance for such studies. In the present study, a field-collected Laodelphax striatellus population from eastern China was laboratory selected for chlorpyrifos resistance and susceptibility, and the developed strains, along with a field population, were studied for cross-resistance and resistance mechanisms at biochemical and molecular levels.

RESULTS

A 158.58-fold chlorpyrifos-resistant strain (JH-chl) and a chlorpyrifos-susceptible strain (JHS) were established after laboratory selection of 25 generations. Cross-resistance to deltamethrin, diazinon, methomyl, carbosulfan, acephate and imidacloprid were detected in JH-chl and a field-collected strain (JHF). Synergism and enzyme activity data suggested potential involvement of P450s and esterases in JH-chl as well as AChE alteration. Furthermore, CYP6AY3v2, CYP306A2v2, CYP353D1v2 and LSCE36 genes were significantly overexpressed in JH-chl (6.87-12.14-fold). Feeding of dsRNAs reduced the expression of the four target genes (35.6-56.8%) and caused significant adult mortality (75.21-88.45%), implying resistance reduction. However, mechanism(s) conferring chlorpyrifos resistance in JHF were unclear.

CONCLUSION

In contrast to previous reports, multiple overexpressed detoxification genes were potentially associated with chlorpyrifos resistance, as confirmed by RNAi feeding tests. Chlorpyrifos resistance exhibits cross-resistance with insecticides in the same and different classes.

Genetic basis of pyrethroid resistance in a population of Anopheles arabiensis, the primary malaria vector in Lower Moshi, north-eastern Tanzania

BACKGROUND

Pyrethroid resistance has been slower to emerge in Anopheles arabiensis than in An. gambiae s.s and An. funestus and, consequently, studies are only just beginning to unravel the genes involved. Permethrin resistance in An. arabiensis in Lower Moshi, Tanzania has been linked to elevated levels of both P450 monooxygenases and β-esterases. We have conducted a gene expression study to identify specific genes linked with metabolic resistance in the Lower Moshi An. arabiensis population.

METHODS

Microarray experiments employing an An. gambiae whole genome expression chip were performed on An. arabiensis, using interwoven loop designs. Permethrin-exposed survivors were compared to three separate unexposed mosquitoes from the same or a nearby population. A subsection of detoxification genes were chosen for subsequent quantitative real-time PCR (qRT-PCR).

RESULTS

Microarray analysis revealed significant over expression of 87 probes and under expression of 85 probes (in pairwise comparisons between permethrin survivors and unexposed sympatric and allopatric samples from Dar es Salaam (controls). For qRT-PCR we targeted over expressed ABC transporter genes (ABC '2060'), a glutathione-S-transferase, P450s and esterases. Design of efficient, specific primers was successful for ABC '2060'and two P450s (CYP6P3, CYP6M2). For the CYP4G16 gene, we used the primers that were previously used in a microarray study of An. arabiensis from Zanzibar islands. Over expression of CYP4G16 and ABC '2060' was detected though with contrasting patterns in pairwise comparisons between survivors and controls. CYP4G16 was only up regulated in survivors, whereas ABC '2060' was similar in survivors and controls but over expressed in Lower Moshi samples compared to the Dar es Salaam samples. Increased transcription of CYP4G16 and ABC '2060' are linked directly and indirectly respectively, with permethrin resistance in Lower Moshi An. arabiensis.

CONCLUSIONS

Increased transcription of a P450 (CYP4G16) and an ABC transporter (ABC 2060) are linked directly and indirectly respectively, with permethrin resistance in Lower Moshi An. arabiensis. Our study provides replication of CYP4G16 as a candidate gene for pyrethroid resistance in An. arabiensis, although its role may not be in detoxification, and requires further investigation.

Cuticular differences associated with aridity acclimation in African malaria vectors carrying alternative arrangements of inversion 2La

BACKGROUND

Principal malaria vectors in Africa, An. gambiae and An. coluzzii, share an inversion polymorphism on the left arm of chromosome 2 (2La/2L+a) that is distributed non-randomly in the environment. Genomic sequencing studies support the role of strong natural selection in maintaining steep clines in 2La inversion frequency along environmental gradients of aridity, and physiological studies have directly implicated 2La in heat and desiccation tolerance, but the precise genetic basis and the underlying behavioral and physiological mechanisms remain unknown. As the insect cuticle is the primary barrier to water loss, differences in cuticle thickness and/or epicuticular waterproofing associated with alternative 2La arrangements might help explain differences in desiccation resistance.

METHODS

To test that hypothesis, two subcolonies of both An. gambiae and An. coluzzii were established that were fixed for alternative 2La arrangements (2La or 2L+a) on an otherwise homosequential and shared genetic background. Adult mosquitoes reared under controlled environmental conditions (benign or arid) for eight days post-eclosion were collected and analyzed. Measurements of cuticle thickness were made based on scanning electron microscopy, and cuticular hydrocarbon (CHC) composition was evaluated by gas chromatography-mass spectrometry.

RESULTS

After removing the allometric effects of body weight, differences in mean cuticle thickness were found between alternative 2La karyotypes, but not between alternative environments. Moreover, the thicker cuticle of the An. coluzzii 2La karyotype was contrary to the known higher rate of water loss of this karyotype relative to 2L+a. On the other hand, quantitative differences in individual CHCs and overall CHC profiles between alternative karyotypes and environmental conditions were consistent with expectation based on previous physiological studies.

CONCLUSIONS

Our results suggest that alternative arrangements of the 2La inversion are associated with differences in cuticle thickness and CHC composition, but that only CHC composition appears to be relevant for desiccation resistance. Differences in the CHC composition were consistent with previous findings of a lower rate of water loss for the 2L+a karyotype at eight days post-eclosion, suggesting that CHC composition is an important strategy for maintaining water balance in this genetic background, but not for 2La. Despite a higher rate of water loss at eight days, higher body water content of the 2La karyotype confers a level of desiccation resistance equivalent to that of the 2L+a karyotype.