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

Does use of domestic insecticides undermine public health control strategies?

Vector-borne diseases (VBD), particularly dengue and malaria, pose a growing threat to human health worldwide. While insecticides remain the cornerstone of vector control programmes, their efficacy is being compromised by increasing insecticide resistance in mosquito populations, leading to control failures that have significant epidemiological and socioeconomic implications. Current research has predominantly examined resistance development in the context of public health interventions and agricultural applications. However, the contribution of domestic insecticide use to resistance evolution in VBD-endemic regions remains inadequately characterised. Evidence indicates that household insecticide utilisation is extensive, with approximately 60% of residents in endemic areas regularly employing domestic insecticidal products for personal protection. This viewpoint highlights how the poorly regulated household insecticide market may significantly contribute to resistance development. Therefore, understanding the impact of domestic insecticide products and usage patterns is urgently needed to preserve the efficacy of vector control campaigns and protect public health outcomes.

Omics Approaches in Understanding Insecticide Resistance in Mosquito Vectors

In recent years, the emergence of insecticide resistance has been a major challenge to global public health. Understanding the molecular mechanisms of this phenomenon in mosquito vectors is paramount for the formulation of effective vector control strategies. This review explores the current knowledge of insecticide resistance mechanisms through omics approaches. Genomic, transcriptomic, proteomic, and metabolomics approaches have proven crucial to understand these resilient vectors. Genomic studies have identified multiple genes associated with insecticide resistance, while transcriptomics has revealed dynamic gene expression patterns in response to insecticide exposure and other environmental stimuli. Proteomics and metabolomics offer insights into protein expression and metabolic pathways involved in detoxification and resistance. Integrating omics data holds immense potential to expand our knowledge on the molecular basis of insecticide resistance in mosquitoes via information obtained from different omics platforms to understand regulatory mechanisms and differential expression of genes and protein, and to identify the transcription factors and novel molecules involved in the detoxification of insecticides. Eventually, these data will help construct predictive models, identify novel strategies, and develop targeted interventions to control vector-borne diseases.

Unveiling candidate genes for metabolic resistance to malathion in Aedes albopictus through RNA sequencing-based transcriptome profiling

Aedes albopictus, also known as the Asian tiger mosquito, is indigenous to the tropical forests of Southeast Asia. Ae. albopictus is expanding across the globe at alarming rates, raising concern over the transmission of mosquito-borne diseases, such as dengue, West Nile fever, yellow fever, and chikungunya fever. Since Ae. albopictus was reported in Houston (Harris County, Texas) in 1985, this species has rapidly expanded to at least 32 states across the United States. Public health efforts aimed at controlling Ae. albopictus, including surveillance and adulticide spraying operations, occur regularly in Harris County. Despite rotation of insecticides to mitigate the development of resistance, multiple mosquito species including Culex quinquefasciatus and Aedes aegypti in Harris County show organophosphate and pyrethroid resistance. Aedes albopictus shows relatively low resistance levels as compared to Ae. aegypti, but kdr-mutation and the expression of detoxification genes have been reported in Ae. albopictus populations elsewhere. To identify potential candidate detoxification genes contributing to metabolic resistance, we used RNA sequencing of field-collected malathion-resistant and malathion-susceptible, and laboratory-maintained susceptible colonies of Ae. albopictus by comparing the relative expression of transcripts from three major detoxification superfamilies involved in malathion resistance due to metabolic detoxification. Between these groups, we identified 12 candidate malathion resistance genes and among these, most genes correlated with metabolic detoxification of malathion, including four P450 and one alpha esterase. Our results reveal the metabolic detoxification and potential cuticular-based resistance mechanisms associated with malathion resistance in Ae. albopictus in Harris County, Texas.

Pyrethroid insecticides susceptibility profiles and evaluation of L1014F kdr mutant alleles in Culex quinquefasciatus from lymphatic filariasis endemic communities

This study investigated the dynamics in pyrethriod resistance and the presence/frequencies of L1014F knockdown resistance mutant allelles in Culex quinquefasciatus vector populations from Uruan Local Government Area of AkwaIbom State, Southern Nigeria between the months of March and November, 2021. Uruan LGA is among the endemic LGAs for lymphatic filariasis in AkwaIbomState. Female Anopheles mosquitoes from Eman Uruan, Ituk Mbang and Idu Uruan were exposed to permethrin, deltamethrin and alphacypermethrin in CDC insecticide coated bottles for susceptibility bioassay following standard protocols. The mosquitoes were obtained as aquatic forms from the study sites and reared under laboratory conditions to adults. The adult mosquitoes were used for this study. All the mosquitoes used for the insecticide susceptibility bioassay were morphologically identified. Standard Polymerase chain reaction (PCR) was used for authenticating the Culex quinquefasciatus species. A portion of the vgsc (917 bp) gene spanning the entire intron and the exon containing the L1014F mutation associated with knockdown resistance (kdr) in the vectorswere amplified using Allele-SPECIFIC POLYMERASE CHAIN REACTION (AS-PCR) in order to detect target site insensitivity in the vectors from the study sites. Results obtained revealed that vectors from all the study sites were resistant to permethrin insecticide (mortality rate: 18-23%). Suspected resistance (mortality rate: 90-93%) to deltamethrin and low resistance (mortality rate: 82-85%) to alphacypermethrin insecticides were detected. knockdown was more rapid with deltamethrin and alphacypermethrin than with permethrin across the study sites considering their KDT50 and KDT95. The frequency of the resistant phenotypes ranged from 35.14 to 55.3% across the study sites with a net of 45.1% resistant phenotype recorded in this study. The 1014F allelic frequency calculated from Hardy-Weinberg principle for vector populations across the study sites ranged from 0.500 (50.00%) to 0.7763 (77.63%). All populations witnessed significant (p < 0.05) deviations from Hardy-Weinberg equilibrium in the distribution of these alleles. The findings of this study show that there is a tendency to record an entire population of resistant vectors in this study area over time due to natural selection. The public health implication of these findings is that the use of pyrethroid based aerosols, coils, sprays, LLITNs and others for the purpose of controlling vectors of lymphatic filariasis and other diseases may be effort in futility.

Efficacy of a 'lethal house lure' against Culex quinquefasciatus from Bouaké city, Côte d'Ivoire

BACKGROUND

Eave tube technology is a novel method of insecticide application that uses an electrostatic coating system to boost insecticide efficacy against resistant mosquitoes. A series of previous experiments showed encouraging insecticidal effects against malaria vectors. This study was undertaken to assess the effects of the eave tube approach on other Culicidae, in particular Culex quinquefasciatus, under laboratory and semi-field conditions.

METHODS

Larvae of Cx. quinquefasciatus from Bouaké were collected and reared to adult stage, and World Health Organization (WHO) cylinder tests were performed to determine their resistance status. WHO standard 3-min cone bioassays were conducted using PermaNet 2.0 netting versus eave tube-treated inserts. To assess the transient exposure effect on Cx. quinquefasciatus, eave tube assay utilizing smelly socks as attractant was performed with exposure time of 30 s, 1 min, and 2 min on 10% beta-cyfluthrin-treated inserts. Residual activity of these treated inserts was then monitored over 9 months. Field tests involving release-recapture of Cx. quinquefasciatus within enclosures around experimental huts fitted with windows and untreated or insecticide-treated eave tubes were conducted to determine house entry preference and the impact of tubes on the survival of this species.

RESULTS

Bouaké Cx. quinquefasciatus displayed high resistance to three out of four classes of insecticides currently used in public health. After 3 min of exposure in cone tests, 10% beta-cyfluthrin-treated inserts induced 100% mortality in Cx. quinquefasciatus, whereas the long-lasting insecticidal net (LLIN) only killed 4.5%. With reduced exposure time on the eave tube insert, mortality was still 100% after 2 min, 88% after 1 min, and 44% after 30 s. Mortality following 1 h exposure on 10% beta-cyfluthrin-treated insert was > 80% continuously up to 7 months post-treatment. Data suggest that Cx. quinquefasciatus have a stronger preference for entering a house through the eaves than through windows. Beta-cyfluthrin-treated inserts were able to kill 51% of resistant Cx. quinquefasciatus released within the enclosure.

CONCLUSIONS

Eave tubes are a novel method for delivery of insecticide to the house. They attract nuisance host-seeking Cx. quinquefasciatus mosquitoes and are as effective in controlling them as they are against pyrethroid-resistant Anopheles gambiae, despite the high level of resistance Cx. quinquefasciatus have developed.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

First Report of Mutations Associated With Pyrethroid (L1014F) and Organophosphate (G119S) Resistance in Belgian Culex (Diptera: Culicidae) Mosquitoes

The emergence of West Nile virus and Usutu virus in Europe poses a significant risk to public health. In the absence of efficient antiviral therapy or vaccine candidates, the only strategy to control these arboviruses is to target the Culex (Diptera: Culicidae) mosquito vector. However, the selection pressure caused by exposure to insecticides for vector control or agricultural pest control can lead to insecticide resistance, thereby reducing the efficacy of insecticide-based vector control interventions. In Culex mosquitoes, two of the most common amino acid substitutions associated with insecticide resistance are the kdr L1014F in voltage gated sodium channels and G119S in acetylcholinesterase. In this study, Culex pipiens biotype pipiens, Culex torrentium, and Culex modestus were sampled from 2019 to 2021 in three distinct environmental habitats (urban, peri-urban, and agricultural) in and around the city of Leuven, Belgium. Individual mosquitoes were screened for two mutations resulting in L1014F and G119S amino acid substitutions. Both mutations were observed in Cx. pipiens and Cx. modestus but not in Cx. torrentium mosquitoes across the four collection sites. Furthermore, multi-resistance or cross-resistance in Cx. pipiens could be a threat in these areas, as both mutations were observed at low frequencies. These results provide the first report of kdr L1014F and ace-1 G119S resistance mutations in Cx. pipiens and Cx. modestus mosquitoes from Belgium, highlighting the importance of mosquito surveillance to design effective arbovirus outbreak control strategies.

Essential functions of mosquito ecdysone importers in development and reproduction

Steroid hormones control sexual maturation and reproduction in insects and humans alike. The insect steroid hormone ecdysone uses a membrane transporter named Ecdysone Importer (EcI) to enter cells and promote these physiological processes, but EcI is unexpectedly missing in mosquito genomes. Using the yellow fever mosquito Aedes aegypti, here we show that mosquitoes use alternative ecdysone importers to facilitate ecdysone-dependent development and reproduction. These transporters are also present in other insects, including fruit flies, but they are dispensable for fly development and reproduction likely due to their limited expression patterns. Our results thus indicate that differential expression of steroid hormone importers enables tissue- and stage-specific hormone responses, and some importers can obtain critical physiological functions only in certain species.

The primary insect steroid hormone ecdysone requires a membrane transporter to enter its target cells. Although an organic anion-transporting polypeptide (OATP) named Ecdysone Importer (EcI) serves this role in the fruit fly Drosophila melanogaster and most likely in other arthropod species, this highly conserved transporter is apparently missing in mosquitoes. Here we report three additional OATPs that facilitate cellular incorporation of ecdysone in Drosophila and the yellow fever mosquito Aedes aegypti. These additional ecdysone importers (EcI-2, -3, and -4) are dispensable for development and reproduction in Drosophila, consistent with the predominant role of EcI. In contrast, in Aedes, EcI-2 is indispensable for ecdysone-mediated development, whereas EcI-4 is critical for vitellogenesis induced by ecdysone in adult females. Altogether, our results indicate unique and essential functions of these additional ecdysone importers in mosquito development and reproduction, making them attractive molecular targets for species- and stage-specific control of ecdysone signaling in mosquitoes.

Emerging Mosquito Resistance to Piperonyl Butoxide-Synergized Pyrethroid Insecticide and Its Mechanism

Piperonyl butoxide (PBO)-synergized pyrethroid products are widely available for the control of pyrethroid-resistant mosquitoes. To date, no study has examined mosquito resistance after pre-exposure to PBO and subsequent enzymatic activity when exposed to PBO-synergized insecticides. We used Culex quinquefasciatus Say (Diptera: Culicidae), an important vector of arboviruses and lymphatic filariasis, as a model to examine the insecticide resistance mechanisms of mosquitoes to PBO-synergized pyrethroid using modified World Health Organization tube bioassays and biochemical analysis of metabolic enzyme expressions pre- and post-PBO exposure. Mosquito eggs and larvae were collected from three cities in Orange County in July 2020 and reared in insectary, and F0 adults were used in this study. A JHB susceptible strain was used as a control. Mosquito mortalities and metabolic enzyme expressions were examined in mosquitoes with/without pre-exposure to different PBO concentrations and exposure durations. Except for malathion, wild strain Cx quinquefasciatus mosquitoes were resistant to all insecticides tested, including PBO-synergized pyrethroids (mortality range 3.7 ± 4.7% to 66.7 ± 7.7%). Wild strain mosquitoes had elevated levels of carboxylesterase (COE, 3.8-fold) and monooxygenase (P450, 2.1-fold) but not glutathione S-transferase (GST) compared to susceptible mosquitoes. When wild strain mosquitoes were pre-exposed to 4% PBO, the 50% lethal concentration of deltamethrin was reduced from 0.22% to 0.10%, compared to 0.02% for a susceptible strain. The knockdown resistance gene mutation (L1014F) rate was 62% in wild strain mosquitoes. PBO pre-exposure suppressed P450 enzyme expression levels by 25~34% and GST by 11%, but had no impact on COE enzyme expression. Even with an optimal PBO concentration (7%) and exposure duration (3h), wild strain mosquitoes had significantly higher P450 enzyme expression levels after PBO exposure compared to the susceptible laboratory strain. These results further demonstrate other studies that PBO alone may not be enough to control highly pyrethroid-resistant mosquitoes due to multiple resistance mechanisms. Mosquito resistance to PBO-synergized insecticide should be closely monitored through a routine resistance management program for effective control of mosquitoes and the pathogens they transmit.