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

The Dual Resistance Mechanism of CYP325G4 and CYP6AA9 in Culex pipiens pallens Legs According to Transcriptome and Proteome Analysis

Mosquitoes within the Culex pipiens complex play a crucial role in human disease transmission. Insecticides, especially pyrethroids, are used to control these vectors. Mosquito legs are the main entry point and barrier for insecticides to gain their neuronal targets. However, the resistance mechanism in mosquito legs is unclear. Herein, we employed transcriptomic analyses and isobaric tags for relative and absolute quantitation techniques to investigate the resistance mechanism, focusing on Cx. pipiens legs. We discovered 2346 differentially expressed genes (DEGs) between deltamethrin-resistant (DR) and deltamethrin-sensitive (DS) mosquito legs, including 41 cytochrome P450 genes. In the same comparison, we identified 228 differentially expressed proteins (DEPs), including six cytochrome P450 proteins. Combined transcriptome and proteome analysis revealed only two upregulated P450 genes, CYP325G4 and CYP6AA9. The main clusters of DEGs and DEPs were associated with metabolic processes, such as cytochrome P450-mediated metabolism of drugs and xenobiotics. Transcription analysis revealed high CYP325G4 and CYP6AA9 expression in the DR strain at 72 h posteclosion compared with that in the DS strain, particularly in the legs. Mosquitoes knocked down for CYP325G4 were more sensitive to deltamethrin than the controls. CYP325G4 knockdown reduced the expression of several chlorinated hydrocarbon (CHC)-related genes, which altered the cuticle thickness and structure. Conversely, CYP6AA9 knockdown increased CHC gene expression without altering cuticle thickness and structure. P450 activity analysis demonstrated that CYP325G4 and CYP6AA9 contributed to metabolic resistance in the midgut and legs. This study identified CYP325G4 as a novel mosquito deltamethrin resistance factor, being involved in both metabolic and cuticular resistance mechanisms. The previously identified CYP6AA9 was investigated for its involvement in metabolic resistance and potential cuticular resistance in mosquito legs. These findings enhance our comprehension of resistance mechanisms, identifying P450s as promising targets for the future management of mosquito vector resistance, and laying a theoretical groundwork for mosquito resistance management.

Knockdown of Thitarodes host genes influences dimorphic transition of Ophiocordyceps sinensis in the host hemolymph

The Chinese cordyceps, a unique parasitic complex of Thitarodes/Hepialus ghost moths and Ophiocordyceps sinensis fungus in the Tibetan Plateau, is a highly valuable biological resource for medicine and health foods in Asian countries. Efficient system for artificial cultivation of Chinese cordyceps relies on understanding the gene functions involved in the induction of growing blastospores into hyphae in the larval hemolymph of insect host, during O. sinensis infection. Transcriptome analysis and ribonucleic acid interference (RNA interference) method were employed to identify the key differentially expressed genes and to demonstrate their functions in Thitarodes xiaojinensis. Key larval genes critical for O. sinensis blastospore development or filamentation were identified. Nine of the 20 top upregulated genes encoded cuticles proteins, indicating that these proteins highly activated when the larval hemolymph was full of blastospores. Small interfering RNA (siRNA) knockdown of five larval genes such as Flightin, larval cuticle protein LCP-30, 26-hydroxylase (CYP18A1), cuticle protein 18.6, isoform B, and probable chitinase 3 significantly stimulated the dimorphic transition from blastospores to prehyphae in O. sinensis in the larval hemolymph after 120 h after injection. The expressions of these genes determined by quantitative real-time PCR were suppressed in various levels from 38.64% to 91.54%, compared to the controls. These results demonstrated that injection of the siRNAs of key upregulated genes into the larval hemolymph containing high load of blastospores caused the gene silence in T. xiaojinensis larvae and induced the fungal transition from blastospores to prehyphae, providing novel knowledge on the regulation of O. sinensis fungal dimorphism by Thitarodes host and cues for further study of Thitarodes biology and commercial cultivation of Chinese cordyceps.

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Identification of insect cuticular protein genes LCP17 and SgAbd5 from Helicoverpa armigera and evaluation their roles in fenvalerate resistance

Insect cuticular protein (ICP) plays an important role in insect growth and development. However, research on the role of ICP in insecticide resistance is very limited. In this study, insect cuticular protein genes LCP17 and SgAbd5 were cloned and characterized in Helicoverpa armigera based on previous transcriptome data. The functions of LCP17 and SgAbd5 genes in fenvalerate resistance were assessed by RNA interference (RNAi), and their response to fenvalerate was further detected. The results showed that LCP17 and SgAbd5 were overexpressed in the fenvalerate-resistant strain comparing with a susceptible strain. The open reading frames of LCP17 and SgAbd5 genes were 423 bp and 369 bp, encoding 141 and 123 amino acids, respectively. LCP17 and SgAbd5 genes were highly expressed in the larval stage, but less expressed in the adult and pupal stages. The expression level of LCP17 and SgAbd5 genes increased significantly after fenvalerate treatment at 24 h. When the cotton bollworms larvae were exposed to fenvalerate at LD50 level, RNAi-mediated silencing of LCP17 and SgAbd5 genes increased the mortality from 50.68% to 68.67% and 63.89%, respectively; the mortality increased to even higher level, which was 73.61%, when these two genes were co-silenced. Moreover, silencing of these two genes caused the cuticle lamellar structure to become loose, which led to increased penetration of fenvalerate into the larvae. The results suggested that LCP17 and SgAbd5 may be involved in the resistance of cotton bollworm to fenvalerate, and LCP17 and SgAbd5 could serve as potential targets for H. armigera control.

A cuticular protein, BgCPLCP1, contributes to insecticide resistance by thickening the cockroach endocuticle

Overuse of insecticides has led to severe environmental problems. Insect cuticle, which consists mainly of chitin, proteins and a thin outer lipid layer, serves multiple functions. Its prominent role is as a physical barrier that impedes the penetration of xenobiotics, including insecticides. Blattella germanica (L.) is a major worldwide indoor pest that causes allergic disease and asthma. Extensive use of pyrethroid insecticides, including β-cypermethrin, has selected for the rapid and independent evolution of resistance in cockroach populations on a global scale. We demonstrated that BgCPLCP1, the first CPLCP (cuticular proteins of low complexity with a highly repetitive proline-rich region) family cuticular protein in order Blattodea, contributes to insecticide penetration resistance. Silencing BgCPLCP1 resulted in 85.0 %-85.7 % and 81.0 %-82.0 % thinner cuticle (and especially thinner endocuticle) in the insecticide-susceptible (S) and β-cypermethrin-resistant (R) strains, respectively. The thinner and more permeable cuticles resulted in 14.4 % and 20.0 % lower survival of β-cypermethrin-treated S- and R-strain cockroaches, respectively. This study advances our understanding of cuticular penetration resistance in insects and opens opportunities for the development of new efficiently and environmentally friendly insecticides targeting the CPLCP family of cuticular proteins.

Transcriptomic investigation of the molecular mechanisms underlying resistance to the neonicotinoid thiamethoxam and the pyrethroid lambda-cyhalothrin in Euschistus heros (Hemiptera: Pentatomidae)

BACKGROUND

Laboratory-selected resistant strains of Euschistus heros to thiamethoxam (NEO) and lambda-cyhalothrin (PYR) were recently reported in Brazil. However, the mechanisms conferring resistance to these insecticides in E. heros remain unresolved. We utilized comparative transcriptome profiling and single nucleotide polymorphism (SNP) calling of susceptible and resistant strains of E. heros to investigate the molecular mechanism(s) underlying resistance.

RESULTS

The E. heros transcriptome was assembled, generating 91 673 transcripts with a mean length of 720 bp and N50 of 1795 bp. Comparative gene expression analysis between the susceptible (SUS) and NEO strains identified 215 significantly differentially expressed (DE) transcripts. DE transcripts associated with the xenobiotic metabolism were all up-regulated in the NEO strain. The comparative analysis of the SUS and PYR strains identified 204 DE transcripts, including an esterase (esterase FE4), a glutathione-S-transferase, an ABC transporter (ABCC1) and aquaporins that were up-regulated in the PYR strain. We identified 9588 and 15 043 nonsynonymous SNPs in the PYR and NEO strains. One of the SNPs (D70N) detected in the NEO strain occurs in a subunit (α5) of the nAChRs, the target site of neonicotinoid insecticides. Nevertheless, this residue position in α5 is not conserved among insects.

CONCLUSIONS

Neonicotinoid and pyrethroid resistance in laboratory-selected E. heros is associated with a potential metabolic resistance mechanism by the overexpression of proteins commonly involved in the three phases of xenobiotic metabolism. Together these findings provide insight into the potential basis of resistance in E. heros and will inform the development and implementation of resistance management strategies against this important pest. © 2023 Society of Chemical Industry.

Development of Sustainable Insecticide Candidates for Protecting Pollinators: Insight into the Bioactivities, Selective Mechanism of Action and QSAR of Natural Coumarin Derivatives against Aphids

Plants employ abundant toxic secondary metabolites to withstand insect attack, while pollinators can tolerate some natural defensive compounds. Coumarins, as promising green alternatives to chemical insecticides, possess wide application prospects in the crop protection field. Herein, the bioactivities of 30 natural coumarin derivatives against Aphis gossypii were assessed and revealed that 6-methylcoumarin exhibited potent aphicidal activity against aphids but displayed no toxicity to honeybees. Additionally, using biochemical, bioinformatic, and molecular assays, we confirmed that the action mode of 6-methylcoumarin against aphids was by inhibiting acetylcholinesterase (AChE). Meanwhile, functional assays revealed that the difference in action site, which located in Lys585 in aphid AChE (equivalent to Val548 in honeybee AChE), was the principal reason for 6-methylcoumarin being toxic to aphids but safe to pollinators. This action site was further validated by mutagenesis data, which uncovered how 6-methylcoumarin was unique selective to the aphid over honeybee or mammalian AChE. Furthermore, a 2D-QSAR model was established, revealing that the central structural feature was H3m, which offers guidance for the future design of more potent coumarin compounds. This work provides a sustainable strategy to take advantage of coumarin analogues for pest management while protecting nontarget pollinators.

Identification and functional analysis of ABC transporter genes related to deltamethrin resistance in Culex pipiens pallens

BACKGROUND

Pathogens that reproduce or develop in mosquitoes can transmit several diseases, endanger human health, and overwhelm health systems. Synthetic pyrethroids are the most widely used insecticides against adult mosquitoes, but their widespread use has led to resistance. The adenosine triphosphate (ATP)-binding cassette (ABC) transporters are involved in the resistance monitoring of insects, but their role and underlying mechanisms in insecticide resistance have not been fully elucidated. In the present study, we identified ABC transporter genes in Culex pipiens and investigated their role in the development of insecticide resistance.

RESULTS

We identified 63 ABC transporter genes in Cx. pipiens and classified them as per the ABC transporter subfamilies. We also performed phylogenetic analysis. The knockdown rate of the mosquitoes orally fed with the ABC transporter inhibitor verapamil increased after deltamethrin treatment compared with that of the control group. Several genes from the ABCB, ABCC, and ABCG subfamilies were highly expressed in resistant mosquitoes. Immunofluorescence analysis revealed that ABCG6032427 was expressed in the head, chest, abdomen, wings, and legs, and the expression was the highest in the legs. Subsequently, knockdown of ABCG6032427 using RNA interference (RNAi) increased the sensitivity of the mosquitoes to deltamethrin, and scanning and transmission electron microscopy revealed that ABCG6032427 knockdown reduced cuticle thickness and the cuticle became loose and irregular.

CONCLUSIONS

ABCG6032427 may modulate cuticle thickness and structure, thus play an important role in the development of deltamethrin resistance in mosquitoes. Thus, it could be a potential target for deltamethrin resistance management in Cx. pipiens. © 2023 Society of Chemical Industry.

Cuticular competing endogenous RNAs regulate insecticide penetration and resistance in a major agricultural pest

BACKGROUND

The continuously developing pesticide resistance is a great threat to agriculture and human health. Understanding the mechanisms of insecticide resistance is a key step in dealing with the phenomenon. Insect cuticle is recently documented to delay xenobiotic penetration which breaks the previous stereotype that cuticle is useless in insecticide resistance, while the underlying mechanism remains scarce.

RESULTS

Here, we find the integument contributes over 40.0% to insecticide resistance via different insecticide delivery strategies in oriental fruit fly. A negative relationship exists between cuticle thickening and insecticide penetration in resistant/susceptible, also in field strains of oriental fruit fly which is a reason for integument-mediated resistance. Our investigations uncover a regulator of insecticide penetration that miR-994 mimic treatment causes cuticle thinning and increases susceptibility to malathion, whereas miR-994 inhibitor results in opposite phenotypes. The target of miR-994 is a most abundant cuticle protein (CPCFC) in resistant/susceptible integument expression profile, which possesses capability of chitin-binding and influences the cuticle thickness-mediated insecticide penetration. Our analyses find an upstream transcriptional regulatory signal of miR-994 cascade, long noncoding RNA (lnc19419), that indirectly upregulates CPCFC in cuticle of the resistant strain by sponging miR-994. Thus, we elucidate the mechanism of cuticular competing endogenous RNAs for regulating insecticide penetration and demonstrate it also exists in field strain of oriental fruit fly.

CONCLUSIONS

We unveil a regulatory axis of lnc19419 ~ miR-994 ~ CPCFC on the cuticle thickness that leads to insecticide penetration resistance. These findings indicate that competing endogenous RNAs regulate insecticide resistance by modulating the cuticle thickness and provide insight into the resistance mechanism in insects.

A chromosome-level genome assembly of Stenchaetothrips biformis and comparative genomic analysis highlights distinct host adaptations among thrips

Insects have a limited host range due to genomic adaptation. Thysanoptera, commonly known as thrips, occupies distinct feeding habitats, but there is a lack of comparative genomic analyses and limited genomic resources available. In this study, the chromosome-level genome of Stenchaetothrips biformis, an oligophagous pest of rice, is assembled using multiple sequencing technologies, including PacBio, Illumina short-reads, and Hi-C technology. A 338.86 Mb genome is obtained, consisting of 1269 contigs with a contig N50 size of 381 kb and a scaffold N50 size of 18.21 Mb. Thereafter, 17,167 protein-coding genes and 36.25% repetitive elements are annotated. Comparative genomic analyses with two other polyphagous thrips, revealing contracted chemosensory-related and expanded stress response and detoxification gene families in S. biformis, potentially facilitating rice adaptation. In the polyphagous thrips species Frankliniella occidentalis and Thrips palmi, expanded gene families are enriched in metabolism of aromatic and anthocyanin-containing compounds, immunity against viruses, and detoxification enzymes. These expansion gene families play crucial roles not only in adapting to hosts but also in development of pesticide resistance, as evidenced by transcriptome results after insecticides treatment. This study provides a chromosome-level genome assembly and lays the foundation for further studies on thrips evolution and pest management.

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.

Fox transcription factor AccGRF1 in response to glyphosate stress in Apis cerana cerana

Glyphosate is an herbicide commonly used in agriculture, and its widespread use has adversely affected the survival of nontarget organisms. Among these organisms, bees in particular are important pollinators, and declining bee populations have severely affected crop yields around the world. However, the molecular mechanism by which glyphosate harms bees remains unclear. In our experiment, we screened and cloned a glyphosate-induced gene in Apis cerana cerana (A. c. cerana) and named glyphosate response factor 1 (AccGRF1). Sequence analysis showed that AccGRF1 contains a winged-helix DNA binding domain, which suggests that it belongs to the Forkhead box (Fox) protein family. qRT-PCR and heterologous expression in Escherichia coli and yeast showed that AccGRF1 can respond to glyphosate and oxidative stress. After AccGRF1 knockdown by means of RNA interference (RNAi), the resistance of A. c. cerana to glyphosate stress improved. The results suggested that AccGRF1 is involved in A. c. cerana glyphosate stress tolerance. This study reveals the functions of Fox transcription factors in response to glyphosate stress and provides molecular insights into the regulation of glyphosate responses in honeybees.

Systematic identification and functional analysis of long noncoding RNAs involved in indoxacarb resistance in Spodoptera litura

Long noncoding RNAs (lncRNAs) are noncoding transcripts that are more than 200 nucleotides long. They play essential roles in regulating a variety of biological processes in many species, including insects, and some lncRNAs have been found to be associated with insecticide resistance. However, the characteristics and biological functions of lncRNAs involved in indoxacarb resistance are unknown in Spodoptera litura. We performed RNA sequencing in the SS, InRS, and FInRS of S. litura and identified 11 978 lncRNAs, including 3 136 intergenic lncRNAs, 7 393 intronic lncRNAs, and 1 449 anti-sense lncRNAs. Compared with the SS, 51 lncRNAs were upregulated and 134 lncRNAs were downregulated in the two resistant strains, and 908 differentially expressed mRNAs were predicted as the target genes of the 185 differentially expressed lncRNAs. Further analysis showed that 112 of differentially expressed lncRNAs may be associated with indoxacarb resistance by regulating the expression of 14 P450s, seven CCEs, one GST, six UGTs, five ABC transporters, and 24 cuticle protein genes, and 79 of differentially expressed lncRNAs may regulate the expression of 14 detoxification genes and 19 cuticle protein genes to participate in indoxacarb resistance by sponging 10 microRNAs. Interestingly, 47 of differentially expressed lncRNAs may mediate indoxacarb resistance through both lncRNA-mRNA and lncRNA-miRNA-mRNA regulatory pathways. Furthermore, quantitative PCR, RNA interference, and indoxacarb bioassay analyses indicated that overexpressed LNC_004867 and LNC_006576 were involved in indoxacarb resistance. This study provides comprehensive information for lncRNAs of S. litura, and presents evidence that lncRNAs have key roles in conferring insecticide resistance in S. litura.

Cuticle Modifications and Over-Expression of the Chitin-Synthase Gene in Diflubenzuron-Resistant Phenotype

Insecticide resistance is a major threat challenging the control of harmful insect species. The study of resistant phenotypes is, therefore, pivotal to understand molecular mechanisms underpinning insecticide resistance and plan effective control and resistance management strategies. Here, we further analysed the diflubenzuron (DFB)-resistant phenotype due to the point-mutation I1043M in the chitin-synthase 1 gene (chs1) in the mosquito Culex pipiens. By comparing susceptible and resistant strains of Cx. pipiens through DFB bioassays, molecular analyses and scanning electron microscopy, we showed that the I1043M-resistant mosquitoes have: (i) a striking level of DFB resistance (i.e., resistance ratio: 9006); (ii) a constitutive 11-fold over-expression of the chs1 gene; (iii) enhanced cuticle thickness and cuticular chitin content. Culex pipiens is one of the most important vector species in Europe and the rapid spread of DFB resistance can threaten its control. Our results, by adding new data about the DFB-resistant phenotype, provide important information for the control and management of insecticide resistance.

Identification of the cuticle protein AccCPR2 gene in Apis cerana cerana and its response to environmental stress

Cuticular proteins (CPs) are known to play important roles in insect development and defence responses. The loss of CP genes can lead to changes in insect morphology and sensitivity to the external environment. In this study, we identified the AccCPR2 gene, which belongs to the CPR family (including the R&R consensus motif) of CPs, and explored its function in the response of Apis cerana cerana to adverse external stresses. Our results demonstrated that AccCPR2 was highly expressed in the late pupal stage and epidermis, and the expression of AccCPR2 may be induced or inhibited under different stressors. RNA interference experiments showed that knockdown of AccCPR2 reduced the activity of antioxidant enzymes, led to the accumulation of oxidative damage and suppressed the expression of several antioxidant genes. In addition, knockdown of AccCPR2 also reduced the pesticide resistance of A. cerana cerana. The overexpression of AccCPR2 in a prokaryotic system further confirmed its role in resistance to various stresses. In summary, AccCPR2 may play pivotal roles in the normal development and environmental stress response of A. cerana cerana. This study also enriched the theoretical knowledge of the resistance biology of bees.

Cuticular protein genes involve heat acclimation of insect larvae under global warming

Cuticular proteins (CPs) play important roles in insect growth and development. However, it is unknown whether CPs are related to heat tolerance. Cnaphalocrocis medinalis, a serious pest of rice, occurs in summer and exhibits strong adaptability to high temperature, but the underlying mechanism is unclear. Here, the role of CP genes in heat acclimation was studied. Heat tolerance of the heat-acclimated larvae was significantly stronger than the unacclimated larvae. The cuticular protein content in the heat-acclimated larvae was higher than that of the unacclimated larvae. 191 presumed CP genes of C. medinalis (CmCPs) were identified. Expression patterns of 14 CmCPs were different between the heat acclimated (S39) and unacclimated (S27) larvae under heat stress. CmCPs were specifically expressed in epidermis and the head except CmCPR20 mainly expressed in Malpighian tubules. CmCPR20 was upregulated in S39 while downregulated in S27, but CmTweedle1 and CmCPG1 were upregulated in S27 and downregulated in S39. RNAi CmTweedle1 or CmCPG1 remarkably decreased heat tolerance and cuticular protein content of the heat-acclimated larvae but not the unacclimated larvae. RNAi CmCPR20 decreased heat tolerance and cuticular protein content of the unacclimated larvae but not the heat-acclimated larvae. CmTweedle1 and CmCPG1 genes involve heat acclimation of C. medinalis.

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

BACKGROUND

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

METHODS AND RESULTS

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

CONCLUSION

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

Distinct molecular impact patterns of abamectin on Apis mellifera ligustica and Apis cerana cerana

The effects of insecticides on bee health are a topic of intensive research. Although abamectin is toxic to bees, the molecular impact of abamectin needs to be clarified. Here, we found that Apis cerana cerana exhibited a higher mortality rate when exposed to abamectin than Apis mellifera ligustica. In addition, A. cerana cerana had markedly higher numbers of differentially expressed genes (DEGs), differentially expressed proteins (DEPs) and differentially expressed metabolites (DEMs) than A. mellifera ligustica during exposure to abamectin. These results indicate that abamectin exposure exerts stronger effects on A. cerana cerana than on A. mellifera ligustica. In addition, six DEGs, two DEPs and two DEMs overlapped between the two bee species under abamectin exposure; however, some genes or proteins from the zinc finger protein, superoxide dismutase and peroxiredoxin families and the energy metabolism pathway were only unregulated in A. cerana cerana, which indicates a significant difference in the impact of abamectin on the two bee species. Despite these differences, several of the same gene families, such as heat shock proteins, cytochrome P450, odorant-binding proteins and cuticle proteins, and pathways, including the carbohydrate metabolism, immune system, lipid metabolism, amino acid metabolism, sensory system, locomotion and development pathways, were influenced by abamectin exposure in both A. cerana cerana and A. mellifera ligustica. Together, our results indicate that abamectin causes adverse effects on bees and thus poses a risk to bee populations and that abamectin exposure affects A. cerana cerana more strongly than A. mellifera ligustica. These findings improve our understanding of the behavioural and physiological effects of abamectin on bees.

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.

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

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

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.

Integration of whole genome sequencing and transcriptomics reveals a complex picture of the reestablishment of insecticide resistance in the major malaria vector Anopheles coluzzii

Insecticide resistance is a major threat to gains in malaria control, which have been stalling and potentially reversing since 2015. Studies into the causal mechanisms of insecticide resistance are painting an increasingly complicated picture, underlining the need to design and implement targeted studies on this phenotype. In this study, we compare three populations of the major malaria vector An. coluzzii: a susceptible and two resistant colonies with the same genetic background. The original colonised resistant population rapidly lost resistance over a 6-month period, a subset of this population was reselected with pyrethroids, and a third population of this colony that did not lose resistance was also available. The original resistant, susceptible and re-selected colonies were subject to RNAseq and whole genome sequencing, which identified a number of changes across the transcriptome and genome linked with resistance. Firstly, an increase in the expression of genes within the oxidative phosphorylation pathway were seen in both resistant populations compared to the susceptible control; this translated phenotypically through an increased respiratory rate, indicating that elevated metabolism is linked directly with resistance. Genome sequencing highlighted several blocks clearly associated with resistance, including the 2Rb inversion. Finally, changes in the microbiome profile were seen, indicating that the microbial composition may play a role in the resistance phenotype. Taken together, this study reveals a highly complicated phenotype in which multiple transcriptomic, genomic and microbiome changes combine to result in insecticide resistance.

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.

MiR-279-3p regulates deltamethrin resistance through CYP325BB1 in Culex pipiens pallens

Background

The overuse of insecticides to control insect vectors has promoted extensive insecticide resistance in mosquitoes. In this study, the functions of microRNA (miR)-279-3p and its target CYP325BB1 in the regulation of deltamethrin resistance in Culex pipiens pallens was investigated.

Methods

Quantitative real-time reverse transcription PCR was used to detect the expression levels of miR-279-3p and CYP325BB1. Then, the dual-luciferase reporter assay system, RNA interference, CDC bottle bioassay and Cell Counting Kit-8 (CCK-8) assay were used to explore the roles of these molecules in deltamethrin resistance both in vivo and in vitro.

Results

The expression patterns of miR-279-3p and CYP325BB1 were compared between deltamethrin-sensitive (DS-strain) and deltamethrin-resistant (DR-strain) mosquitoes. Luciferase activity was downregulated by miR-279-3p, the effect of which was ablated by a mutation of the putative binding site for CYP325BB1. In DR-strain mosquitoes, the expression of miR-279-3p was increased by microinjection and oral feeding of miR-279-3p agomir (mimic). CYP325BB1 mRNA levels were downregulated, which resulted in a higher mortality of the mosquitoes in miR-279-3p mimic-treated groups. In the DS-strain mosquitoes, microinjection of a miR-279-3p inhibitor decreased miR-279-3p expression, whereas the expression of CYP325BB1 was increased; the mortality of these mosquitoes decreased significantly. In addition, overexpression of pIB/V5-His-CYP325BB1 changed the sensitivity of C6/36 cells to deltamethrin in vitro. Also in DR-strain mosquitoes, downregulation of CYP325BB1 expression by microinjection of si-CYP325BB1 increased mosquito mortality in vivo.

Conclusions

These findings provide empirical evidence of the involvement of miRNAs in the regulation of insecticide resistance and indicate that miR-279-3p suppresses the expression of CYP325BB1, which in turn decreases deltamethrin resistance, resulting in increased mosquito mortality. Taken together, the results provide important information for use in the development of future mosquito control strategies.

Graphical abstract

Identification and expression analysis of cuticular protein genes in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae)

Structural cuticular proteins (CPs) are major components of the insect cuticle, and they play critical roles in insect development and insecticide resistance. Here, a total of 196 CP genes were successfully annotated in the Plutella xylostella genome. On the basis of motif analysis, these CPs were classified into 10 different families, including 122 CPR, 12 CPAP1, 8 CPAP3, 9 CPLCP, 2 Tweedle, 1 CPF, 1 CPFL, 1 CPCFC, 17 CPG and 2 18 aa proteins, and the remaining 21 unclassified CPs were classed as cuticular proteins hypothetical (CPH). A phylogenetic analysis of CPs from different insects revealed species-specific clades of RR-1 and RR-2 genes, suggesting that CP gene duplication might occur independently among insect taxa, while we also found that some other CPs (such as CPAP1 and CPAP3) had a closer relationship based on their conserved domain architecture. Using available RNAseq libraries, the expression profiles of the CPs were analyzed over the four developmental stages of the insect (i.e., egg, larva, pupa, and adult), revealing stage-specific expression patterns for the CPs. In a chlorpyrifos resistant strain, 18 CP genes were found to be more than two-fold upregulated compared to the susceptible control strain, and qRT-PCR analysis showed that these CP genes were overexpressed after exposure to chlorpyrifos, suggesting a potential role in the molecular mechanism of insecticide resistance in P. xylostella. This study provides the tools and molecular basis to study the role of CPs in the post-embryonal development and the mechanisms of insecticide resistance of P. xylostella.

Behavioral and ultrastructural effects of novaluron on Aedes aegypti larvae

Chitin synthesis inhibitors (CSI) are supposed to inhibit formation of chitin microfibrils in newly synthesized cuticle during molting process. Conversely, there has been comparatively few data on morphological effects of CSI on non-target insect organs. In this work, the effects of the CSI novaluron on behavior and midgut of A. aegypti were evaluated. Toxicity bioassays revealed that novaluron is toxic to A. aegypti larva with LC₅₀ = 18.57 mg L^-1 when exposed in aqueous solution for 24 h. Novaluron treated larvae were less active and spent more time resting compared to the control group. Histopathology showed that midguts of novaluron-treated larvae had cytoplasm vacuolization and damaged brush border. Cytotoxic effects in midguts of treated larvae induced necrosis, autophagy and damage to mitochondria. Despite being chitin synthesis inhibitor, novaluron did not induce alterations in the integument of A. aegypti larvae. Fluorescence microscopy revealed that the number of digestive cells were higher in novaluron-treated larvae than in control, in response to digestive cell apoptosis. The present study highlights the importance of novaluron against A. aegypti larvae by causing injuries to non-target organs, altering behaviors, inducing cell death and inhibiting cell proliferation.

Transcriptomic analysis reveals pronounced changes in gene expression due to sub-lethal pyrethroid exposure and ageing in insecticide resistance Anopheles coluzzii

BACKGROUND

Malaria control is heavily reliant on the use of insecticides that target and kill the adult female Anopheline vector. The intensive use of insecticides of the pyrethroid class has led to widespread resistance in mosquito populations. The intensity of pyrethroid resistance in some settings in Africa means mosquitoes can contact bednets treated with this insecticide class multiple times with minimal mortality effects. Furthermore, both ageing and diel cycle have been shown to have large impacts on the resistance phenotype. Together, these traits may affect other aspects of vector biology controlling the vectorial capacity or fitness of the mosquito.

RESULTS

Here we show that sublethal exposure of a highly resistant Anopheles coluzzii population originally from Burkina Faso to the pyrethroid deltamethrin results in large and sustained changes to transcript expression. We identify five clear patterns in the data showing changes to transcripts relating to: DNA repair, respiration, translation, behaviour and oxioreductase processes. Further, we highlight differential regulation of transcripts from detoxification families previously linked with insecticide resistance, in addition to clear down-regulation of the oxidative phosphorylation pathway both indicative of changes in metabolism post-exposure. Finally, we show that both ageing and diel cycle have major effects on known insecticide resistance related transcripts.

CONCLUSION

Sub-lethal pyrethroid exposure, ageing and the diel cycle results in large-scale changes in the transcriptome of the major malaria vector Anopheles coluzzii. Our data strongly supports further phenotypic studies on how transcriptional changes such as reduced expression of the oxidative phosphorylation pathway or pyrethroid induced changes to redox state might impact key mosquito traits, such as vectorial capacity and life history traits.

Genome-wide identification of long non-coding RNAs (lncRNAs) associated with malathion resistance in Bactrocera dorsalis

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in the regulation of biological processes and have been identified in many species including insects. However, the association between lncRNAs and pesticide resistance in insect species such as Bactrocera dorsalis is unknown.

RESULTS

RNA-seq was performed on malathion resistant (MR1) and susceptible (MS) strains of B. dorsalis and a total of 6171 lncRNAs transcripts were identified. These included 3728 lincRNAs, 653 antisense lncRNAs, 1402 intronic lncRNAs, and 388 sense lncRNAs. A total of 40 and 52 upregulated lncRNAs were found in females and males of the MR1 strain compared to 54 and 49 in the same sexes of the MS strain, respectively. Twenty-seven of these lncRNAs showed the same trend of expression in both females and males in the MR1 strain, in which 15 lncRNAs were upregulated and 12 were downregulated. RT-qPCR results indicated that the differentially expressed lncRNAs were associated with malathion resistance. The lnc15010.10 and lnc3774.2 were highly expressed in the cuticle of the MR1 strain, indicating that these two lncRNAs may be related to malathion resistance. RNAi of lnc3774.2 and a bioassay showed that malathion resistance was possibly influenced by changes in the B. dorsalis cuticle.

CONCLUSION

LncRNAs of B. dorsalis potentially related to the malathion resistance were identified. Two lncRNAs appear to influence malathion resistance via modulating the structure, or components, of the cuticle. © 2021 Society of Chemical Industry.

High insecticide resistance mediated by different mechanisms in Culex quinquefasciatus populations from the city of Yaoundé, Cameroon

Culex mosquitoes particularly Culex quinquefasciatus are important arboviral and filariasis vectors, however despite this important epidemiological role, there is still a paucity of data on their bionomics. The present study was undertaken to assess the insecticide resistance status of Cx. quinquefasciatus populations from four districts of Yaoundé (Cameroon). All Culex quinquefasciatus populations except one displayed high resistance to bendiocarb and malathion with mortalities ranging from 0 to 89% while high resistance intensity against both permethrin and deltamethrin was recorded. Molecular analyses revealed high frequencies of the ACE-1 G119S mutation (ranging from 0 to 33%) and kdr L1014F allele (ranging from 55 to 74%) in all Cx. quinquefasciatus populations. Significant overexpression was detected for cytochrome P450s genes CYP6AA7 and CYP6Z10, as well as for Esterase A and Esterase B genes. The total cuticular hydrocarbon content, a proxy of cuticular resistance, was significantly increased (compared to the S-lab strain) in one population. The study confirms strong insecticide resistance mediated by different mechanisms in Cx. quinquefasciatus populations from the city of Yaoundé. The expansion of insecticide resistance in Culex populations could affect the effectiveness of current vector control measures and stress the need for the implementation of integrated vector control strategies in urban settings.

Phenotypic Insecticide Resistance in Anopheles gambiae (Diptera: Culicidae): Specific Characterization of Underlying Resistance Mechanisms Still Matters

An effective control of malaria vectors requires an extensive knowledge of mechanisms underlying the resistance-phenotypes developed by these vectors against insecticides. We investigated Anopheles gambiae mosquitoes from Benin and Togo for their intensity of insecticide resistance and we discussed the involvement of genotyped mechanisms in the resistance-phenotypes observed. Three- to five-day-old adult mosquitoes emerged from field and laboratory An. gambiae larvae were assayed using WHO tube intensity tests against various doses of deltamethrin: 1× (0.05%); 2× (0.1%); 5× (0.25%); 7.5× (0.375%) and those of pirimiphos-methyl: 0.5× (0.125%); 1× (0.25%). Members of An. gambiae complex were screened in field populations using polymerase chain reaction (PCR) assays. The presence of kdrR(1014F/1014S) and ace-1R(119S) mutations was also investigated using TaqMan and PCR-RFLP techniques, respectively. Anopheles gambiae from field were very resistant to deltamethrin, whereas KisKdr and AcerKdrKis strains displayed 100% mortality rates at 2× the diagnostic dose. In contrast, the field mosquitoes displayed a low resistance-intensity against 1× the diagnostic dose of pirimiphos-methyl, whereas AcerKis and AcerKdrKis strains showed susceptibility at 0.5× the diagnostic dose. Anopheles gambiae s.s., Anopheles coluzzii, and Anopheles arabiensis were identified. Allelic frequencies of kdrR (1014F) and ace-1R (119S) mutations in the field populations varied from 0.65 to 1 and 0 to 0.84, respectively. The field An. gambiae displayed high-resistance levels against deltamethrin and pirimiphos-methyl when compared with those of the laboratory An. gambiae-resistant strains. These results exhibit the complexity of underlying insecticide resistance mechanisms in these field malaria vectors.

RNA-seq Analysis Reveals Mitochondrial and Cuticular Protein Genes Are Associated with Phosphine Resistance in the Rusty Grain Beetle (Coleoptera:Laemophloeidae)

The rusty grain beetle, Cryptolestes ferrugineus (Stephens), is a serious pest of stored grain, which has developed high levels of resistance to phosphine. In this study, five geographically distant populations of C. ferrugineus had been collected in China, specifically in granaries where phosphine fumigant is used for pest control, and they showed a high resistance ratio up to 1,907 (LC50 = 21.0 mg/liter). Then, a reference transcriptome was constructed to use as a basis for investigating the molecular mechanisms of phosphine resistance in this species, which consisted of 47,006 unigenes with a mean length of 1,090. Subsequently, the RNA-Seq analysis of individuals from the most susceptible and resistant populations led to the identification of 54 genes that are differentially expressed. GO and KEGG analysis demonstrated that genes associated with mitochondrial and respiration functions were significantly enriched. Also, the 'structural constituent of cuticle' term was annotated in the GO enrichment analysis and further qRT-PCR confirmed that the expression levels of nine cuticular protein genes were significantly increased in the resistant population. In conclusion, we present here a transcriptome-wide overview of gene expression changes between resistant and susceptible populations of C. ferrugineus, and this in turn documents that mitochondria and cuticular protein genes may play together a crucial role in phosphine resistance. Further gene function analysis should enable the provision of advice to expedite resistance management decisions.

Transcriptomic response of Anopheles gambiae sensu stricto mosquito larvae to Curry tree (Murraya koenigii) phytochemicals

BACKGROUND

Insect growth regulators (IGRs) can control insect vector populations by disrupting growth and development in juvenile stages of the vectors. We previously identified and described the curry tree (Murraya koenigii (L.) Spreng) phytochemical leaf extract composition (neplanocin A, 3-(1-naphthyl)-L-alanine, lumiflavine, terezine C, agelaspongin and murrayazolinol), which disrupted growth and development in Anopheles gambiae sensu stricto mosquito larvae by inducing morphogenetic abnormalities, reducing locomotion and delaying pupation in the mosquito. Here, we attempted to establish the transcriptional process in the larvae that underpins these phenotypes in the mosquito.

METHODS

We first exposed third-fourth instar larvae of the mosquito to the leaf extract and consequently the inherent phytochemicals (and corresponding non-exposed controls) in two independent biological replicates. We collected the larvae for our experiments sampled 24 h before peak pupation, which was 7 and 18 days post-exposure for controls and exposed larvae, respectively. The differences in duration to peak pupation were due to extract-induced growth delay in the larvae. The two study groups (exposed vs control) were consequently not age-matched. We then sequentially (i) isolated RNA (whole larvae) from each replicate treatment, (ii) sequenced the RNA on Illumina HiSeq platform, (iii) performed differential bioinformatics analyses between libraries (exposed vs control) and (iv) independently validated the transcriptome expression profiles through RT-qPCR.

RESULTS

Our analyses revealed significant induction of transcripts predominantly associated with hard cuticular proteins, juvenile hormone esterases, immunity and detoxification in the larvae samples exposed to the extract relative to the non-exposed control samples. Our analysis also revealed alteration of pathways functionally associated with putrescine metabolism and structural constituents of the cuticle in the extract-exposed larvae relative to the non-exposed control, putatively linked to the exoskeleton and immune response in the larvae. The extract-exposed larvae also appeared to have suppressed pathways functionally associated with molting, cell division and growth in the larvae. However, given the age mismatch between the extract-exposed and non-exposed larvae, we can attribute the modulation of innate immune, detoxification, cuticular and associated transcripts and pathways we observed to effects of age differences among the larvae samples (exposed vs control) and to exposures of the larvae to the extract.

CONCLUSIONS

The exposure treatment appears to disrupt cuticular development, immune response and oxidative stress pathways in Anopheles gambiae s.s larvae. These pathways can potentially be targeted in development of more efficacious curry tree phytochemical-based IGRs against An. gambiae s.s mosquito larvae.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

Transcriptomic and metabolomic landscape of the molecular effects of glyphosate commercial formulation on Apis mellifera ligustica and Apis cerana cerana

Understanding the causes of the decline in bee population has attracted intensive attention worldwide. The indiscriminate use of agrochemicals is a persistent problem due to their physiological and behavioural damage to bees. Glyphosate and its commercial formulation stand out due to their wide use in agricultural areas and non-crop areas, such as parks, railroads, roadsides, industrial sites, and recreational and residential areas, but the mode of action of glyphosate on bees at the molecular level remains largely unelucidated. Here, we found that the numbers of differentially expressed genes and metabolites under glyphosate commercial formulation (GCF) stress were significantly higher in Apis cerana cerana than in Apis mellifera ligustica. Despite these differences, the number of differentially expressed transcripts increased following an increase in the GCF treatment time in both A. cerana cerana and A. mellifera ligustica. GCF exerted adverse impacts on the immune system, digestive system, nervous system, amino acid metabolism, carbohydrate metabolism, growth and development of both bee species by influencing their key genes and metabolites to some extent. The expression of many genes involved in immunity, agrochemical detoxification and resistance, such as antimicrobial peptides, cuticle proteins and cytochrome P450 families, was upregulated by GCF in both bee species. Collectively, our results indicate that both A. cerana cerana and A. mellifera ligustica strive to mitigate the pernicious effects caused by GCF by regulating detoxification and immune systems. Moreover, A. cerana cerana might be better able to withstand the toxic effects of GCF with lower fitness costs than A. mellifera ligustica. Our work will contribute to elucidating the deleterious physiological and behavioural impacts of GCF on bees.

Ovarian morphological features and proteome reveal fecundity fitness disadvantages in β-cypermethrin-resistant strains of Blattella germanica (L.) (Blattodea: Blattellidae)

To evaluate whether the development of β-cypermethrin resistance in Blattella germanica (L.) (Blattaria: Blattellidae) affects the fecundity fitness of this insect and to determine the underlying mechanism, we compared fecundity differences between β-cypermethrin-resistant (R) and sensitive (S) strains of B. germanica, observed the physiological structural changes of ovaries from an visual perspective, and analyzed differences in the ovarian proteome using proteomic methods. The results showed that, compared with the S strain of B. germanica, the R strain of B. germanica had a significantly higher ootheca shedding rate, a significantly lower number of hatched and surviving nymphs, a significantly higher female proportion in the population and defective ovarian development. Ovarian proteomic analysis showed a total of 64 differentially expressed proteins in the R strain, including 18 upregulated proteins and 46 downregulated proteins. Twenty-four significantly differentially expressed proteins were further studied, and 14 were successfully identified, which were mainly classified into the following categories: immunity-related proteins, development-related proteins, structural proteins, energy metabolism-related proteins and proteins with unknown functions. The differential expression of these proteins reflects the overall changes in cell structure and metabolism associated with β-cypermethrin resistance and explains the possible molecular mechanism of fecundity fitness disadvantages. In summary, β-cypermethrin resistance can cause fecundity fitness disadvantages in B. germanica. The metabolic deviations needed to overcome the adverse effects of insecticides may result in an energy exchange that affects energy allocation and, ultimately, the basic needs of the insect. The fitness cost due to insecticide resistance is critical to the delay of the evolution of resistance.

Transcription factor FTZ-F1 regulates mosquito cuticular protein CPLCG5 conferring resistance to pyrethroids in Culex pipiens pallens

Background

Culex pipiens pallens poses a serious threat to human health because of its widespread distribution, high carrier capacity for several arboviruses, frequent human-biting, and growth in urban environments. Pyrethroid insecticides have been mainly used to control adult Cx. pipiens pallens during outbreaks of mosquito-borne diseases. Unfortunately, mosquitoes have developed resistance, rendering the insecticides ineffective. Cuticular resistance is the primary mechanism of pyrethroid resistance. Previously, we revealed that cuticular protein of low complexity CPLCG5 is a major cuticular protein associated with deltamethrin resistance in Cx. pipiens pallens, which is enriched in the cuticle of mosquitoes’ legs and participates in pyrethroid resistance by forming a rigid matrix. However, the regulatory mechanisms of its transcription remain unknown.

Results

First, qRT-PCR analysis revealed that the expression of FTZ-F1 (encoding Fushi tarazu-Factor 1) was ~ 1.8-fold higher in the deltamethrin-resistant (DR) than deltamethrin-susceptible (DS) strains at 24 h post-eclosion (PE) and ~ 2.2-fold higher in the DR strain than in the DS strain at 48 h PE. CPLCG5 and FTZ-F1 were co-expressed in the legs, indicating that they might play an essential role in the legs. Dual luciferase reporter assays and EMSA (electrophoretic mobility shift experiments) revealed that FTZ-F1 regulates the transcription of CPLCG5 by binding to the FTZ-F1 response element (− 870/− 864). Lastly, knockdown of FTZ-F1 not only affected CPLCG5 expression but also altered the cuticle thickness and structure of the legs, increasing the susceptibility of the mosquitoes to deltamethrin in vivo.

Conclusions

The results revealed that FTZ-F1 regulates the expression of CPLCG5 by binding to the CPLCG5 promoter region, altering cuticle thickness and structure, and increasing the susceptibility of mosquitoes to deltamethrin in vivo. This study revealed part of the mechanism of cuticular resistance, providing a deeper understanding of insecticide resistance.

Involvement of integument-rich CYP4G19 in hydrocarbon biosynthesis and cuticular penetration resistance in Blattella germanica (L.)

BACKGROUND

Cuticle penetration plays an important role as a mechanism of insecticide resistance, but the underlying molecular mechanism remains poorly understood. In Blattella germanica, the cytochrome P450 gene, CYP4G19, is overexpressed in a pyrethroid-resistant strain. Here, we investigated whether CYP4G19 is involved in the biosynthesis of hydrocarbons and further contributes to cuticular penetration resistance in B. germanica.

RESULTS

Compared with the susceptible strain, pyrethroid-resistant cockroaches showed lower cuticular permeability with Eosin Y staining. Removal of epicuticular lipids, mainly nonpolar hydrocarbons, with a hexane wash intensified the cuticular permeability and decreased the resistance index of the resistant strain. CYP4G19 was predominately expressed in the abdominal integument and could be upregulated by desiccation stress or short exposure to beta-cypermethrin. Overexpression of CYP4G19 in the resistant strain was positively correlated with a higher level of cuticular hydrocarbons (CHCs). RNAi-mediated knockdown of CYP4G19 significantly decreased its expression and caused a reduction in CHCs. Meanwhile, CYP4G19 suppression resulted in a non-uniform array of the lipid layer, enhanced cuticle permeability, and compromised insecticide tolerance.

CONCLUSION

Our findings confirm that CYP4G19 is involved in hydrocarbon production and appears to contribute to hydrocarbon-based penetration resistance in B. germanica. This study highlights the lipid-based penetration resistance, advancing our understanding of the molecular mechanism underlying P450-mediated cuticular penetration resistance in insects. © 2019 Society of Chemical Industry.

Transcriptome Analysis and Identification of Insecticide Tolerance-Related Genes after Exposure to Insecticide in Sitobion avenae

Aphids cause serious losses to the production of wheat. The grain aphid, Sitobion avenae, which is the dominant species of aphid in all wheat regions of China, is resistant to a variety of insecticides, including imidacloprid and chlorpyrifos. However, the resistance and mechanism of insecticide tolerance of S. avenae are still unclear. Therefore, this study employed transcriptome analysis to compare the expression patterns of stress response genes under imidacloprid and chlorpyrifos treatment for 15 min, 3 h, and 36 h of exposure. S. avenae adult transcriptome was assembled and characterized first, after which samples treated with insecticides for different lengths of time were compared with control samples, which revealed 60–2267 differentially expressed unigenes (DEUs). Among these DEUs, 31–790 unigenes were classified into 66–786 categories of gene ontology (GO) functional groups, and 24–760 DEUs could be mapped into 54–268 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Finally, 11 insecticide-tolerance-related unigenes were chosen to confirm the relative expression by quantitative real-time polymerase chain reaction (qRT-PCR) in each treatment. Most of the results between qRT-PCR and RNA sequencing (RNA-Seq) are well-established. The results presented herein will facilitate molecular research investigating insecticide resistance in S. avenae, as well as in other wheat aphids.

Mosquitoes cloak their legs to resist insecticides

Malaria incidence has halved since the year 2000, with 80% of the reduction attributable to the use of insecticides. However, insecticide resistance is now widespread, is rapidly increasing in spectrum and intensity across Africa, and may be contributing to the increase of malaria incidence in 2018. The role of detoxification enzymes and target site mutations has been documented in the major malaria vector Anopheles gambiae; however, the emergence of striking resistant phenotypes suggests the occurrence of additional mechanisms. By comparing legs, the most relevant insect tissue for insecticide uptake, we show that resistant mosquitoes largely remodel their leg cuticles via enhanced deposition of cuticular proteins and chitin, corroborating a leg-thickening phenotype. Moreover, we show that resistant female mosquitoes seal their leg cuticles with higher total and different relative amounts of cuticular hydrocarbons, compared with susceptible ones. The structural and functional alterations in Anopheles female mosquito legs are associated with a reduced uptake of insecticides, substantially contributing to the resistance phenotype.

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.

Culex species diversity, susceptibility to insecticides and role as potential vector of Lymphatic filariasis in the city of Yaoundé, Cameroon

BACKGROUND

Culex species are widespread across Cameroon and responsible for high burden of nuisance in most urban settings. However, despite their high nuisance, they remain less studied compared to anophelines. The present study aimed to assess Culex species distribution, susceptibility to insecticide, bionomics and role in Lymphatic Filariasis (LF) transmission in the city of Yaoundé.

METHODS

Mosquito collections were conducted from March to December 2017 using Centre for Disease Control light traps (CDC-LT), human landing catches (HLC) and larval collections. Mosquitoes were identified using morphological identification keys. Mosquitoes from the Culex pipiens complex were further identified using Polymerase Chain Reaction (PCR) to assess the presence of sibling species. Bioassays were conducted with 2-5 day-old unfed females to assess mosquito susceptibility to DDT, permethrin, deltamethrin and bendiocarb following WHO guidelines. Dead, control and surviving mosquitoes from bioassays were screened by PCR to detect the presence of knockdown resistance (kdr) alleles. Pools of mosquitoes were examined by PCR to detect the presence of Wuchereria bancrofti.

RESULTS

A total of 197,956 mosquitoes belonging to thirteen species were collected. The density of mosquito collected varied according to the collection methods, districts and seasons. Culex quinquefasciatus emerged as the most abundant and the only species of the Culex pipiens complex in Yaoundé. Culex species were found breeding in different types of breeding sites including polluted and unpolluted sites. All Culex species including Cx antennatus, Cx duttoni, Cx perfuscus and Cx tigripes were found to be highly resistant to permethrin, deltamethrin and DDT. Culex quinquefasciatus was also found to be resistant to bendiocarb. A high frequency of the West Africa kdr allele was recorded in resistant Cx. quinquefasciatus. Out of the 247 pooled samples of 25 Culex spp. examined for the presence of Wuchereria bancrofti, none was found infected.

CONCLUSION

The study confirms the high adaptation of Culex species particularly Culex quinquefasciatus to the urban environment and no implication of this species in the transmission of LF in Yaoundé Cameroon. Culex species predominance in urban settings highlight potential transmission risk of West Nile and rift valley fever in Yaoundé.

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

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

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.