- 1Institute of Disinfection and Pest Control, Beijing Center for Disease Prevention and Control, Beijing, China
- 2State Key Laboratory of Integrated Management of Pest Insects and Rodents, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
Background: Aedes (Stegomyia) albopictus (Skuse) is an invasive and widespread mosquito species that can transmit dengue, chikungunya, yellow fever, and Zika viruses. Its control heavily relies on the use of insecticides. However, the efficacy of the insecticide-based intervention is threatened by the increasing development of resistance to available insecticides. Understanding the current status and potential mechanisms of insecticide resistance is an important prerequisite for devising strategies to maintain the sustainability of vector control programs. In this study, we investigated the current status and probable candidate detoxification genes associated with insecticide resistance in the Asian tiger mosquito in Beijing, the capital city of China.
Methods: Bioassays were conducted on three field populations of Ae. albopictus collected from urban communities in Beijing by exposure to diagnostic doses of permethrin, deltamethrin, malathion, and propoxur. Differentially expressed genes (DEGs) associated with insecticide resistance were screened by transcriptomic analysis using Illumina RNA sequencing data (RNA-seq) from 12 independent RNA libraries constructed from female strains of the three field populations and one susceptible strain.
Results: The bioassay results indicated that all the three field populations were resistant to propoxur (carbamate), deltamethrin, and permethrin (pyrethroids), but susceptible to malathion (organophosphate). Eighteen (18) cytochrome P450s (P450s), five (5) glutathione S-transferases (GSTs), four (4) carboxy/cholinesterases (CCEs), eight (8) UDP-glycosyltransferases (UGTs), and three (3) ATP-binding cassette transporters (ABCs) were found to be significantly overexpressed in the three field populations relative to the susceptible strain via transcriptomic analysis.
Conclusion: This study demonstrates that the Ae. albopictus field populations in Beijing exhibit multiple phenotypic resistance to commonly used pyrethroids and carbamate. The identification of a number of DEGs associated with insecticide resistance indicates that the mechanisms underlying resistance in field populations are complicated, and detoxifying enzymes may play important roles. The multiple resistance status detected in the three field populations suggests that resistance management strategies such as insecticide rotation and non-chemical-based measures should be implemented in order to sustain effective control of the disease vector and vector-borne diseases.
1 Introduction
The Asian tiger mosquito Aedes albopictus (Skuse) is a highly invasive disease vector. It has been documented that this vector is able to transmit dengue, chikungunya, yellow fever, and Zika viruses through biting humans (Paixão et al., 2018) and act as a bridge vector for zoonotic pathogens to humans due to the wide host range from primary mammalian hosts to birds, reptiles, and amphibians. With limited antiviral drugs and vaccines, control of Aedes vectors is the primary practice to prevent dengue, Zika, and chikungunya transmission. Chemical insecticides have been used successfully for many years to control Aedes vectors and vector-borne diseases, due to their high efficacy and speed of action. However, the sustainability of insecticide-based vector control faces a major obstacle: the ever-increasing resistance to the few number of insecticides registered for vector control (Brogdon and McAllister, 1998).
Resistance to insecticides in Aedes vectors has become a global problem. Previous studies have shown that Aedes mosquitoes have developed resistance to the three classes of commonly used insecticides (carbamates, organophosphates, and pyrethroids) in addition to the abandoned organochlorides worldwide (Whalon et al., 2015; Moyes et al., 2017; Dusfour et al., 2015; Dusfour et al., 2019; Zulfa et al., 2022), including in China (Liu et al., 2020; Wei et al., 2021; Li et al., 2018; Li et al., 2021).
Resistance mechanisms in Aedes are multifaceted, involving target site insensitivity, metabolic detoxification, and cuticular modifications, each contributing differently depending on environmental pressures and genetic factors. One well-established mechanism is target site insensitivity, particularly knockdown resistance (kdr) mutations in the voltage-gated sodium channel (VGSC). A common mutation such as V1016G reduces the effectiveness of pyrethroids, one of the most widely used insecticide classes (Kasai et al., 2014). Metabolic resistance, though less explored in Ae. albopictus compared to other mosquito species such as Ae. aegypti (Schluep and Buckner, 2021) and Anopheles gambiae (Gueye et al., 2020), involves the overexpression of detoxifying enzymes such as cytochrome P450 monooxygenases, glutathione S-transferases (GSTs), and carboxylesterases. These enzymes degrade or sequester insecticides before they can reach their target sites. Recent studies suggest that specific P450 genes, such as those in the CYP9J subfamily, may play a role in pyrethroid resistance in Ae. aegypti (Dusfour et al., 2015). Cuticular resistance, involving exoskeleton changes, reduces insecticide absorption mainly by slowing its penetration due to the upregulation of cuticular protein gene expression or altered chitin biosynthesis of chitin or hydrocarbon in Ae. aegypti and An. gambiae (Jacobs et al., 2023; Balabanidou et al., 2019). The limited knowledge of metabolic resistance, especially in terms of specific enzyme contributions, hinders our understanding of how Ae. albopictus is adapting to insecticide pressures. Further research is needed to identify the mechanisms occurring in insecticide-resistant Ae. albopictus populations, which is crucial for developing effective control strategies.
Beijing is the capital city of China with a population of approximately 22 million. Thanks to the Healthy Beijing Program, mosquito control campaigns have been carried out since 2016 during the rainy seasons by means of larval source management and insecticide application (residual and space sprays). As a consequence, resistance to insecticides in the target mosquitoes is gradually increasing (unpublished data). Understanding the current status and mechanisms of insecticide resistance is an important prerequisite for the subsequent development of strategies to slow the development of resistance and keep the control program sustainable. As efforts in this regard, we previously completed an extensive survey on target site mutation-mediated resistance and found that multiple mutations were present in the voltage-gated sodium channel (VGSC) gene in the Ae. albopictus populations in Beijing (Zhou et al., 2019a). However, information on the resistance status is limited, and other possible resistance mechanisms remain unknown in Ae. albopictus field populations in Beijing. In this study, we attempted to identify genes that may contribute to metabolic resistance by genome-wide transcriptome analysis on three field populations of Ae. albopictus sampled from urban communities of Beijing.
2 Materials and methods
2.1 Mosquito samples
Larvae of Ae. albopictus were collected from three sampling sites in Beijing from July to August in 2019. Heping Street community (DT) in Chaoyang District and Ganjiakou community (GJK) in Haidian District were residential communities, and mosquito collection in both communities was authorized with permission from the households. JK population mosquitoes were collected from the courtyard garden of Beijing Center for Disease Prevention and Control (Beijing CDC) in Dongcheng District, surrounded by a park, residential buildings, a school, and a stadium, with permission of the authorized representative. Surroundings (within 100 m radius) of the sampling points were thoroughly searched for possible breeding sites of Aedes mosquitoes. All identified artificial and natural water bodies (such as tree holes, rock pools, leaf axils, bamboo stumps, and small ponds) and containers (such as bottles, tires, cans, paper cups, barrels, foam boxes, and pots) likely to harbor mosquito larvae were visually searched thoroughly for the presence of Aedes mosquito larvae. Aedes larvae were sampled using the dipping method. At least five breeding containers were used for collection at each sample site, and approximately 300 larvae were collected overall for each site/population. The susceptible strain (SS) of Ae. albopictus from Beijing CDC was kept in a laboratory environment without exposure to any insecticide since 1984 and showed no tolerance to different insecticides. All samples were maintained at 25°C–27°C, 60%–80% relative humidity, and a 14:10 h light:dark photoperiod. Larvae were fed with mouse feed, and adults were provided with a 10% sucrose solution and mouse for blood meal. Ae. albopictus specimens were identified morphologically and confirmed molecularly by the PCR method based on the rDNA-ITS sequence (Higa et al., 2010). Insecticide tolerance bioassays and transcriptome analyses were completed within the third generation (F3) after mosquito collection.
2.2 Insecticide resistance bioassays
All collected mosquito larvae were reared to adult mosquitoes in the laboratory, and bioassays were completed as soon as possible. Non-blood fed 3- to 5-day-old female adults of the first to third generation (F1 to F3) were used directly for bioassays without extra insecticide pre-treatment. To prepare the insecticide-impregnated paper, the insecticide was diluted to the diagnostic concentration using a mixture of liquid paraffin and ether (1:2 v/v). Next, 2 mL of the prepared insecticide solution was applied evenly onto a 15 cm × 12 cm filter paper. After allowing the ether to evaporate completely, the impregnated paper may be stored in airtight, light-protected containers or used immediately for bioassay procedures. Adults were exposed to 0.03% deltamethrin-, 0.4% permethrin-, 0.5% malathion-, and 0.05% propoxur-impregnated filter membranes according to WHO tube bioassay protocols (WHO, 2013), and the discriminating (or diagnostic) concentrations were most frequently used for Ae. albopictus in China (Wu et al., 2022). Meanwhile, the susceptible strain mosquito was assayed by control filter (only with the solvent and without the active ingredient) in parallel. All the test and control filters were supplied by the Division of Vector Control, Chinese Center for Disease Control and Prevention. At least triple replicates were performed per insecticide, with no less than 30 adults per replicate. After 1-h exposure, mosquitoes were transferred to a recovery tube and maintained on 10% sucrose solution for 24 h when the mortality rates were recorded. If mortality in the control exceeded 5%, Abbott’s correction was applied. According to WHO guidelines, the resistance status of mosquito populations was scored as “susceptible” if mortality at 24 h after exposure was ≥98%, “possibly resistant” if mortality ranged between 90% and 97%, and “resistant” if mortality was ≤90% (WHO, 2013; Chang et al., 2014).
2.3 RNA extraction, RNA-seq library preparation, sequencing, and functional annotation
Twelve (three field populations and one susceptible strain, with three replicates/pools for each population or strain) pools of 30 unfed female adults were stored in liquid nitrogen for subsequent RNA extraction. All the mosquitoes used in subsequent transcriptomic sequencing were not challenged with insecticides. Total RNA was extracted from each pool using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s protocol. Recombinant DNase I (Takara, Japan) was used to remove potential genomic DNA. RNA integrity was evaluated using NanoDrop 2000 (Thermo Fisher Scientific, Delaware, ME, United States). RNA quality was analyzed on an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, United States).
The RNA-seq libraries were constructed with 500 μg of starting total RNA with the Illumina TruSeq Stranded mRNA Sample Preparation Kit (Illumina, San Diego, CA) following the kit protocol. The resulting libraries were sequenced on an Illumina HiSeq 2500 instrument to generate strand-specific, paired-end reads of length 125 bp (HiSeq SBS Kit v4 sequencing reagents). Library construction, sequencing, and read trimming were performed at the high-throughput genomics and bioinformatic analysis platform of Majorbio (Majorbio, Shanghai, China). According to FastQC version 0.11.4, no reads were tagged as poor quality. All the reads generated in Illumina-Solexa sequencing were deposited in the NCBI Sequence Read Archive (SRA) database with BioProject accession number PRJNA1110376.
The software Fastx_toolkit_0.0.14 was used to analyze the raw transcriptome data to obtain clean data, which was further aligned with the reference genome (GCF_006496715.1) using the software HISAT2 (2.2.1). The returned SAM format files were assembled using Samtools (1.9) and StringTie (2.1.2) to obtain the whole transcriptome. The assembled transcripts were annotated against NR, SWISS-Prot, Pfam, COG, GO, and KEGG.
2.4 Analysis of differentially expressed genes (DEGs) and functional enrichment
Gene expression level analysis was carried out using the software RSEM (1.3.3) with the transcripts per million reads (TPM) index. Expression differences of genes among different samples were analyzed using DESeq2, where the laboratory susceptible strain (SS) was used as the control and the three field-resistant populations (DT, GJK, and JK) were used as the treatment groups. The criteria for significant differences in gene expression were set as |log2Ratio| ≥ 1 and p-value < 0.05, where p-value was corrected by false discovery rate (FDR) ≤ 0.001. If there was more than one transcript for a gene, the longest transcript was used to calculate its expression level and coverage.
KEGG functional enrichment analysis of significantly differentially expressed genes in field populations was performed by Fisher’s exact test using the bioinformatics cloud platform developed by Majorbio (https://cloud.majorbio.com/). To exclude false positives for significant enrichment of KEGG Pathway, the Benjamini–Hochberg (BH) method (Benjamini et al., 2001) was used for correction with p-value ≤ 0.05.
In order to screen for DEGs probably involved in metabolic resistance, the following analysis pipeline was applied. First, structural domain screening was performed among the DEGs using HMMER (3.2.1) and the Pfam database, with the PF00067 domain for cytochrome P450s, PF00043 and PF02798 domains for GSTs, the PF00135 domain for CCEs, the PF00201 domain for UGTs, and the PF00005 domain for ABCs, and a total of five metabolic resistance-related gene families were searched. Second, the retrieved genes were filtered for a number of amino acid residues based on the protein sequences obtained from ORF translation using the awk command in Linux OS, and CYP450s: 450 aa–600 aa; GSTs: 200 aa–300 aa; CCEs: 500 aa–700 aa; UGTs: 500 aa–700 aa; CCEs: 500 aa–700 aa. 700 aa; UGTs: 500 aa–700 aa; ABC-transporters: 500 aa–2,000 aa corresponding genes. Then, the genes obtained from step 2 were filtered according to their expression levels using Python scripts, and genes with TPM less than 10 in all four populations were excluded.
3 Results
3.1 Phenotypic resistance
Female strains aged 3–5 days from the three field populations (GJK, DT, and JK) and the laboratory susceptible strain (SS) were bioassayed to evaluate their susceptibility to four commonly used insecticides at discriminating concentrations. In all the controls, no knock-down and dead mosquito was observed after the 1-h exposure to solvent-treated filter membrane and 24-h normal feeding. For 0.5% malathion, possible resistance (96.85% mortality) was only detected in GJK, while other two populations exhibited a mortality of more that 98%. For 0.05% propoxur, a mortality of 74.02%, 57.85%, and 50.76% was observed in DT, GJK, and JK, respectively. The three field populations exhibited significant resistance, with the mortality ranging from 61.96% to 80.13% and from 68.32% to 86.79% after exposure to 0.4% permethrin and 0.03% deltamethrin, respectively. Overall, the highest resistance level against the tested insecticides was detected in JK, except for malathion. These results indicated that the three Ae. albopictus field populations were resistant to permethrin, propoxur, and deltamethrin, but remain largely susceptible to malathion (Table 1).
Table 1. Resistance status of Ae. albopictus populations from Beijing to four insecticides by WHO tube bioassays.
3.2 Transcriptomes obtained from four Ae. albopictus populations
Sequencing the 12 transcriptomes produced 622,916,642 paired-end raw reads in total. After analysis of base composition and sequence quality using the software Fastx_toolkit_0.0.14, 618,559,408 clean reads were obtained. The error rates were below 0.03%, and Q30 percentages were > 93% (Supplementary Material 1). The clean reads reached 83.13%–85.06% coverage when mapped to the Ae. albopictus reference genome (GCF_006496715.1) (Supplementary Material 2).
3.3 Screening of differentially expressed genes associated with insecticide resistance
By comparing the transcriptomes, 2,063, 1,674, and 2,340 significantly upregulated (Ur) genes and 1,503, 1,424, and 1,745 significantly downregulated (Dr) genes were observed in the field-resistant DT, GJK, and JK populations, respectively, relative to the laboratory susceptible strain (SS) (Figure 1A). Considering the difference of the genetic background among the four populations, we focused on the shared subsets of genes that were differentially transcribed in all the three pairwise comparisons. In this manner, a total of 586 genes were significantly upregulated and 512 genes were significantly downregulated in the three field populations as compared with the SS strain (Figure 1B). Furthermore, KEGG pathway enrichment analyses on these commonly up- or downregulated genes in the three field populations were especially performed to further explore their potential biological functions and pathways, and the top 20 enrichment results are presented in Figure 1C. The enrichment results indicated that metabolism of xenobiotics by cytochrome P450, chemical carcinogenesis, and drug metabolism–cytochrome P450 pathways were the common significantly enriched pathways in the three field-resistant populations (Figure 2).
Figure 1. Identification of differentially expressed genes in three field populations (DT, GJK, and JK) compared to the laboratory susceptible strain. (A) Volcano plots (volcano plot) of differential gene expression between each of the field populations and the laboratory-sensitive populations (SS); red points represent significantly upregulated (Ur) expression genes [log2 (fold change) ≥ 1 and p-value ≤ 0.05], and blue points represent significantly downregulated (Dr) expression genes [log2 (fold change) ≤ −1 and p-value ≤ 0.05]. (B) Venn analysis of significantly differentially expressed genes in three field-resistant populations; the left panel shows significantly upregulated differentially expressed genes in field-resistant populations, and the right panel shows significantly downregulated differentially expressed genes in field-resistant populations. (C) The left side illustrates the enrichment results of genes that are commonly upregulated in three field populations, while the right side depicts the enrichment results of genes that are commonly downregulated in those same populations. The horizontal axis denotes the KEGG pathway, while the vertical axis signifies the Rich factor. The gradient of the bar colors indicates the significance of the enrichment, with a Padjust < 0.001 marked as ***, Padjust < 0.01 marked as **, and Padjust < 0.05 marked as *.
Figure 2. KEGG pathway enrichment analysis of significantly upregulated gene expressions in three field populations relative to the laboratory susceptible strain. The vertical axis represents the KEGG pathway, and the horizontal axis represents the Rich factor [ratio of (number of genes/transcripts enriched in the pathway) to (number of annotated genes/transcripts)]. The size of the dots indicates the number of genes in the pathway, and the color of the dots corresponds to different Padjust ranges (the smaller the Padjust, the higher the enrichment). Only the KEGG pathway with Padjust ≤ 0.05 and the top 20 enrichment are shown.
Focusing on xenobiotics detoxification-related genes, we found that a total of 18 P450, 5 GST, 4 CCE, 8 UGT, and 3 ABC genes were significantly overexpressed in the three field populations relative to the susceptible strain (Table 2). Among the P450s, there are 2 P450s from the CYP4 family, 6 from the CYP6 family, and 10 from the CYP9 family (9 from the CYP9J subfamily). Specifically, JK, DT, and GJK overexpressed 33, 26, and 17 detoxification-related genes, respectively. Notably, 3 P450 (CYP9J18, CYP9J32, and CYP9J74), 2 GST (GSTe2a and GSTI1), 1 CCE (COesterase6a), and 3 UGT (UGT36D1a, UGT36D1b, and UGT49C1a) genes were concurrently overexpressed in the three field-resistant populations. In addition to the above-mentioned, 8 P450s (CYP6AG17, CYP6N4, CYP6Z8, CYP6Z23, CYP9J19, CYP9J26, CYP9J65, and CYP9J68), 1 GSTs (GSTX2b), 2 CCEs (COesterase5a and COesterase6b), 3 UGTs (UGT35E2, UGT49B1, and UGT49C1b), and 3 ABCs (ABCG1, ABCG4, and ABCG1-X1) were overexpressed in two field populations.
Table 2. Xenobiotics detoxification genes that are up-regulated in the three field Ae. Albopictus populations relative to the susceptible strain.
4 Discussion
Insecticide resistance in Ae. albopictus is now threatening the global fight against arboviral human diseases such as dengue, yellow fever, chikungunya, and Zika (Asgarian et al., 2023). In China, insecticide resistance in Ae. albopictus is rapidly emerging in most areas (Li et al., 2018; Gao et al., 2018; Su et al., 2019; Yuan et al., 2023; Zhao et al., 2023). Therefore, monitoring and management of insecticide resistance is becoming a crucial element in vector control. The purpose of this study was to evaluate the insecticide resistance status and screen for genes potentially involved in resistance in Ae. albopictus field populations in the urban areas of Beijing.
Bioassay results reveal that the three examined populations have developed varying degrees of resistance to the commonly used pyrethroid and carbamate insecticides, and the extent of resistance is related to the intensity of insecticide usage against different Ae. albopictus populations (Table 1). For example, in JK, where insecticide sprays were applied three times each month from July to September during the years 2016–2019, Ae. albopictus developed the highest resistance to deltamethrin, permethrin, and propoxur. By contrast, no significant resistance to malathion was detected in all the three populations, probably because organophosphate insecticides have been used less frequently. The observations that these populations were susceptible to malathion, but resistant to propoxur, may suggest that propoxur resistance is attributed to enhanced detoxification by enzymes that selectively catalyze propoxur rather than mutations in their shared target acetylcholinesterase.
Similarly, Ae. albopictus populations in the urban areas of Guangzhou in southern China have developed resistance to deltamethrin, permethrin, propoxur, and bendiocarb, but they have remained susceptible to malathion (Li et al., 2018; Su et al., 2019); field populations of Ae. albopictus from Zhejiang Province in eastern China were widely resistant to pyrethroids but sensitive to organophosphate insecticides (Chang et al., 2014). A recent study also reported that seven populations from Guangyuan City of Sichuan Province in western China are susceptible to malathion but resistant to beta-cypermethrin and deltamethrin (Zhao et al., 2024). The above findings from different provinces of China suggest that most urban Ae. albopictus mosquitoes may be subjected to similar selective pressures of insecticides and consequently exhibit similar profiles of insecticide resistance. Although the annual most active period of Ae. albopictus in Beijing is relatively short (from July to August), the large-scale mosquito control campaigns in urban areas have imposed continuing pressure to select for insecticide-resistant populations of Ae. albopictus. The prevalence of phenotypic resistance in Ae. albopictus populations in Beijing strongly suggests that reduced use or rotation of insecticides should be considered in future practices of vector control and highlights the necessity to understand the genetic mechanism underlying multiple insecticide resistance.
It has been extensively accepted that enhanced metabolic detoxification of insecticides is a major mechanism in insecticide resistance in mosquitoes (Liu, 2015). To identify genes that are associated with insecticide resistance, we sequenced 12 Ae. albopictus transcriptomes. A total of 586 upregulated and 512 downregulated genes were commonly observed in the DK, GJK, and JK populations as compared with the SS strain (Figure 1). Insecticide resistance often incurs a metabolic cost, so downregulated genes may reflect suppression of non-essential pathways to conserve resources. KEGG pathway enrichment analysis of DEGs revealed that several detoxification genes were upregulated in the field-resistant populations. There are 186 P450, 64 CCE, 32 GST, and 71 ABC transporter genes in the Ae. albopictus genome (Chen et al., 2015), of which 18 P450, 5 GST, 4 CCE, 8 UGT, and 3 ABC genes were significantly expressed in the three field populations relative to the susceptible strain, indicating that detoxification-mediated mechanisms may play crucial roles in the phenotypic resistance to pyrethroids (Marcombe et al., 2024) and carbamates in Ae. albopictus.
P450s are well recognized for their roles in insecticide metabolism and resistance. Of the 18 P450 genes that are overexpressed in at least one of the three field populations, 9 genes are members of the CYP9J subfamily (Table 2). P450 genes belonging to the CYP9J subfamily (CYP9J6, CYP9J8, CYP9J9, CYP9J10, CYP9J19, CYP9J22, CYP9J23, CYP9J24, CYP9J26, CYP9J27, CYP9J28, CYP9J31, and CYP9J32) were also reported to be overexpressed in several pyrethroid-resistant populations of Ae. aegypti (Strode et al., 2008; Saavedra-Rodriguez et al., 2012; Stevenson et al., 2012; Rahman et al., 2021; Bariami et al., 2012). Moreover, functional studies in vitro demonstrate that some CYP9J members (e.g., CYP9J24, CYP9J26, CYP9J28, and CYP9J32) are capable of metabolizing pyrethroids permethrin or deltamethrin (Stevenson et al., 2012). Moreover, we observed two members (CYP6Z8 and CYP6Z23) of the CYP6Z subfamily over-transcribed in two (DT and JK) of the three field populations. The involvement of CYP6Z in insecticide resistance in other mosquito species was also reported. For example, CYP6Z8 and CYP6Z9 were found to play a pivotal role in pyrethroid clearance in Ae. aegypti (Bariami et al., 2012; Chandor-Proust et al., 2013); three CYP6Z genes in Ae. albopictus were involved in the resistance to deltamethrin in vivo (Zou et al., 2024). In addition, upregulation of CYP6AG17 was also observed in DT and JK populations in this study. Three CYP6AG (Dusfour et al., 2015; Dusfour et al., 2019; Zulfa et al., 2022) homologous genes were associated with deltamethrin resistance in field Ae. aegypti (Saavedra-Rodriguez et al., 2019). These observations coupled with the mRNA levels (Table 2) make us hypothesize that CYP6Z8, CYP6Z23, and CYP9J18 may be key players in P450-mediated insecticide resistance in Ae. albopictus urban populations in Beijing.
In this study, five glutathione S-transferases (GSTs) were observed to be overexpressed in the 3 field populations (Table 2). GSTs may contribute to the development of resistance to all main classes of insecticides via direct metabolism and/or sequestration of insecticides or via indirectly providing protection against oxidative stress induced by insecticide exposure (Pavlidi et al., 2018; Zhai et al., 2024). GSTs were found to be related with pyrethroid resistance in other species of mosquitoes such as Ae. aegypti (Strode et al., 2008; Lumjuan et al., 2011; Seixas et al., 2017; Leong et al., 2019), An. sinensis (Zhu et al., 2014), and An. gambiae (Bonizzoni et al., 2015). GST activities were enhanced in most Ae. albopictus populations from Sri Lanka resistant to pyrethroids (NWNP et al., 2021). Notably, the abundances of GSTe2a and GSTI1 in the three field populations were significantly higher than those of the susceptible strain (Table 4), indicating that these two GSTs may play crucial roles in resistance.
The CCE family contains both catalytically active enzymes and non-catalytic proteins and can act by rapid binding or sequestration to affect the interactions between insecticides and targets (Mourya et al., 1993; Karunaratne et al., 1995; Safi et al., 2019). In the current study, four CCEs were observed to be overexpressed in at least one of the three field populations (Table 2). Notably, COesterase6a had the highest abundance and was commonly overexpressed in all the three resistant populations, highlighting its importance in resistance. Interestingly, the homologous COesterase5a and COesterase6a genes were reported to undergo an independent amplification respectively and spread between two countries in field-resistant Ae. albopictus populations (Grigoraki et al., 2017).
In addition to the above-mentioned detoxification enzymes, an increasing number of studies have indicated the involvement of UGTs and ABC transporters in insecticide resistance (Bariami et al., 2012; Dermauw and Van Leeuwen, 2014; Lv et al., 2016; He et al., 2019; Rault et al., 2019; Zhou et al., 2019b; Bruckmueller and Cascorbi, 2021; Israni et al., 2022; Chen et al., 2023). In this study, we found that 8 UGTs and 3 ABCs were upregulated in the three resistant populations relative to the susceptible strain of Ae. albopictus (Table 2). Further research is required to clarify the role of these DEGs in resistance.
This study has some limitations. First, the diagnostic dosage method could only determine whether the field population are resistant or susceptible to the insecticides based on survival rates at the discriminating dose, without quantifying the resistance levels or detecting low levels of resistance or early-stage resistance development. Second, the differences in some genes expressions observed between the susceptible strain and field populations could be due to the differences in the genetic background and/or living (rearing) conditions rather than insecticide resistance per se. Including a field-susceptible strain would help address these potential confounding factors and provide a more robust analysis of the genetic basis of insecticide resistance in Ae. albopictus. This approach would provide a more accurate baseline for comparison. However, to our knowledge, it is quite difficult to find a susceptible strain sourced from the field due to the massive insecticide use historically. Another limitation of the study was the lack of the supporting biochemical data or synergist such as piperonyl butoxide (PBO) and S, S, S-tributyl phosphorothioate (DEF) assay. However, our research group is currently conducting in vivo and in vitro functional validations on three overexpressed P450 genes, aiming to clarify their contribution to insecticide tolerance.
In summary, the availability of Ae. albopictus genomes (Chen et al., 2015; Palatini et al., 2020) has greatly benefited transcriptome analyses. The high-quality transcriptomes obtained in the study allow us to make a list of potential genes involved in metabolic resistance to insecticides in the field-resistant populations. Our results support the notion that the three major detoxification enzyme families (P450s, CCEs, and GSTs) together with ABCs and UGTs are involved in insecticide resistance in Aedes mosquitoes (Kasai et al., 2014; Chen et al., 2015; Strode et al., 2008; Zhu et al., 2014; Lv et al., 2016; Liu et al., 2011; Simma et al., 2019; Lien et al., 2019; Messenger et al., 2021). The exact contribution of the above-mentioned detoxification-related genes to resistance will require further work. Such effort is underway in our lab.
5 Conclusion
Bioassay results demonstrate that all the examined Ae. albopictus field populations in Beijing have developed multiple phenotypic resistance to commonly used pyrethroids and carbamate. The identification of several DEGs associated with insecticide resistance indicates that the mechanisms underlying resistance in field populations are complicated. The observation that multiple detoxification-related genes are upregulated in all the three resistance populations strongly suggests that enhanced detoxification mediated by increased expression of detoxifying enzymes may play important roles in the development of multiple resistance. The significant resistance status detected in all the three populations highlights that resistance management strategies such as insecticide rotation and/or non-chemical-based measures should be implemented in order to sustain effective control of the disease vector and vector-borne diseases.
Data availability statement
All the reads generated in transcriptome sequencing were deposited in the NCBI Sequence Read Archive (SRA) database with BioProject accession number PRJNA1110376.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
XZ: funding acquisition, investigation, methodology, writing–original draft, and writing–review and editing. JL: investigation, methodology, project administration, and writing–review and editing. RN: data curation, methodology, software, and writing–original draft. XQ: conceptualization, funding acquisition, methodology, and writing–review and editing. YZ: resources, supervision, validation, and writing–review and editing. YT: project administration, supervision, validation, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by grants from the Beijing Municipal Health Commission High-level Public Health Technical Talent Construction Program (Discipline backbone-01-010, to XZ), the Beijing Natural Science Foundation (52222020, to XQ), and the State Key Laboratory of Integrated Management of Pest Insects and Rodents (IPM2305, to XZ).
Acknowledgments
The authors appreciate David R. Nelson of the University of Tennessee Health Science Center for his great help in the nomenclature of P450 genes in Aedes albopictus.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2024.1498313/full#supplementary-material
References
Asgarian T. S., Vatandoost H., Hanafi-Bojd A. A., Nikpoor F. (2023). Worldwide status of insecticide resistance of Aedes aegypti and Ae. Albopictus, vectors of arboviruses of chikungunya, dengue, zika and yellow fever. J. Arthropod Borne Dis. 17, 1–27. doi:10.18502/jad.v17i1.13198
Balabanidou V., Kefi M., Aivaliotis M., Koidou V., Girotti J. R., Mijailovsky S. J., et al. (2019). Mosquitoes cloak their legs to resist insecticides. Proc. Biol. Sci. 286, 20191091. doi:10.1098/rspb.2019.1091
Bariami V., Jones C. M., Poupardin R., Vontas J., Ranson H. (2012). Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti. PLoS Negl. Trop. Dis. 6, e1692. doi:10.1371/journal.pntd.0001692
Benjamini Y., Drai D., Elmer G., Kafkafi N., Golani I. (2001). Controlling the false discovery rate in behavior genetics research. Behav. Brain Res. 125, 279–284. doi:10.1016/s0166-4328(01)00297-2
Bonizzoni M., Ochomo E., Dunn W. A., Britton M., Afrane Y., Zhou G., et al. (2015). RNA-seq analyses of changes in the Anopheles gambiae transcriptome associated with resistance to pyrethroids in Kenya: identification of candidate-resistance genes and candidate-resistance SNPs. Parasit. Vectors 8, 474. doi:10.1186/s13071-015-1083-z
Brogdon W. G., McAllister J. C. (1998). Insecticide resistance and vector control. Emerg. Infect. Dis. 4 (4), 605–613. doi:10.3201/eid0404.980410
Bruckmueller H., Cascorbi I. (2021). ABCB1, ABCG2, ABCC1, ABCC2, and ABCC3 drug transporter polymorphisms and their impact on drug bioavailability: what is our current understanding? Expert Opin. Drug Metab. Toxicol. 17, 369–396. doi:10.1080/17425255.2021.1876661
Chandor-Proust A., Bibby J., Régent-Kloeckner M., Roux J., Guittard-Crilat E., Poupardin R., et al. (2013). The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling. Biochem. J. 455, 75–85. doi:10.1042/BJ20130577
Chang X., Zhong D., Fang Q., Hartsel J., Zhou G., Shi L., et al. (2014). Multiple resistances and complex mechanisms of Anopheles sinensis mosquito: a major obstacle to mosquito-borne diseases control and elimination in China. PLoS Negl. Trop. Dis. 8, e2889. doi:10.1371/journal.pntd.0002889
Chen M. L., Zhang S. X., Guo P. Y., Qin Q. S., Meng L. W., Yuan G. R., et al. (2023). Identification and characterization of UDP-glycosyltransferase genes and the potential role in response to insecticides exposure in Bactrocera dorsalis. Pest Manag. Sci. 79, 666–677. doi:10.1002/ps.7234
Chen X. G., Jiang X., Gu J., Xu M., Wu Y., Deng Y., et al. (2015). Genome sequence of the Asian Tiger mosquito, Aedes albopictus, reveals insights into its biology, genetics, and evolution. Proc. Natl. Acad. Sci. U. S. A. 112, E5907–E5915. doi:10.1073/pnas.1516410112
Dermauw W., Van Leeuwen T. (2014). The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. Insect Biochem. Mol. Biol. 45, 89–110. doi:10.1016/j.ibmb.2013.11.001
Dusfour I., Vontas J., David J. P., Weetman D., Fonseca D. M., Corbel V., et al. (2019). Management of insecticide resistance in the major Aedes vectors of arboviruses: Advances and challenges. PLoS Negl. Trop. Dis. 13, e0007615. doi:10.1371/journal.pntd.0007615
Dusfour I., Zorrilla P., Guidez A., Issaly J., Girod R., Guillaumot L., et al. (2015). Deltamethrin resistance mechanisms in Aedes aegypti populations from three French overseas territories worldwide. PLoS Negl. Trop. Dis. 9, e0004226. doi:10.1371/journal.pntd.0004226
Gao J. P., Chen H. M., Shi H., Peng H., Ma Y. J. (2018). Correlation between adult pyrethroid resistance and knockdown resistance (kdr) mutations in Aedes albopictus (Diptera: Culicidae) field populations in China. Infect. Dis. Poverty 7, 86. doi:10.1186/s40249-018-0471-y
Grigoraki L., Pipini D., Labbé P., Chaskopoulou A., Weill M., Vontas J. (2017). Carboxylesterase gene amplifications associated with insecticide resistance in Aedes albopictus: geographical distribution and evolutionary origin. PLoS Negl. Trop. Dis. 11, e0005533. doi:10.1371/journal.pntd.0005533
Gueye O. K., Tchouakui M., Dia A. K., Faye M. B., Ahmed A. A., Wondji M. J., et al. (2020). Insecticide resistance profiling of Anopheles coluzzii and Anopheles gambiae populations in the southern Senegal: role of target sites and metabolic resistance mechanisms. Genes (Basel) 11, 1403. doi:10.3390/genes11121403
He Q., Yan Z., Si F., Zhou Y., Fu W., Chen B. (2019). ATP-binding cassette (ABC) transporter genes involved in pyrethroid resistance in the malaria vector Anopheles sinensis: genome-wide identification, characteristics, phylogenetics, and expression profile. Int. J. Mol. Sci. 20, 1409. doi:10.3390/ijms20061409
Higa Y., Toma T., Tsuda Y., Miyagi I. (2010). A multiplex PCR-based molecular identification of five morphologically related, medically important subgenus Stegomyia mosquitoes from the genus Aedes (Diptera: Culicidae) found in the Ryukyu Archipelago, Japan. Jpn. J. Infect. Dis. 63, 312–316. doi:10.7883/yoken.63.312
Israni B., Luck K., Römhild S. C. W., Raguschke B., Wielsch N., Hupfer Y., et al. (2022). Alternative transcript splicing regulates UDP-glucosyltransferase-catalyzed detoxification of DIMBOA in the fall armyworm (Spodoptera frugiperda). Sci. Rep. 12, 10343. doi:10.1038/s41598-022-14551-w
Jacobs E., Chrissian C., Rankin-Turner S., Wear M., Camacho E., Scott J. G., et al. (2023). Cuticular profiling of insecticide resistant Aedes aegypti. Sci. Rep. 13, 10154. doi:10.1038/s41598-023-36926-3
Karunaratne S. H., Hemingway J., Jayawardena K. G., Dassanayaka V., Vaughan A. (1995). Kinetic and molecular differences in the amplified and non-amplified esterases from insecticide-resistant and susceptible Culex quinquefasciatus mosquitoes. J. Biol. Chem. 270, 31124–31128. doi:10.1074/jbc.270.52.31124
Kasai S., Komagata O., Itokawa K., Shono T., Ng L. C., Kobayashi M., et al. (2014). Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism. PLoS Negl. Trop. Dis. 8, e2948. doi:10.1371/journal.pntd.0002948
Leong C. S., Vythilingam I., Liew J. W., Wong M. L., Wan-Yusoff W. S., Lau Y. L. (2019). Enzymatic and molecular characterization of insecticide resistance mechanisms in field populations of Aedes aegypti from Selangor, Malaysia. Parasit. Vectors 12, 236. doi:10.1186/s13071-019-3472-1
Li Y., Xu J., Zhong D., Zhang H., Yang W., Zhou G., et al. (2018). Evidence for multiple-insecticide resistance in urban Aedes albopictus populations in southern China. Parasit. Vectors 11, 4. doi:10.1186/s13071-017-2581-y
Li Y., Zhou G., Zhong D., Wang X., Hemming-Schroeder E., David R. E., et al. (2021). Widespread multiple insecticide resistance in the major dengue vector Aedes albopictus in Hainan Province, China. Pest Manag. Sci. 77, 1945–1953. doi:10.1002/ps.6222
Lien N. T. K., Ngoc N. T. H., Lan N. N., Hien N. T., Tung N. V., Ngan N. T. T., et al. (2019). Transcriptome sequencing and analysis of changes associated with insecticide resistance in the dengue mosquito (Aedes aegypti) in vietnam. Am. J. Trop. Med. Hyg. 100, 1240–1248. doi:10.4269/ajtmh.18-0607
Liu H., Liu L., Cheng P., Yang L., Chen J., Lu Y., et al. (2020). Bionomics and insecticide resistance of Aedes albopictus in Shandong, a high latitude and high-risk dengue transmission area in China. Parasit. Vectors 13, 11. doi:10.1186/s13071-020-3880-2
Liu N. (2015). Insecticide resistance in mosquitoes: impact, mechanisms, and research directions. Annu. Rev. Entomol. 60, 537–559. doi:10.1146/annurev-ento-010814-020828
Liu N., Li T., Reid W. R., Yang T., Zhang L. (2011). Multiple Cytochrome P450 genes: their constitutive overexpression and permethrin induction in insecticide resistant mosquitoes, Culex quinquefasciatus. PLoS One 6, e23403. doi:10.1371/journal.pone.0023403
Lumjuan N., Rajatileka S., Changsom D., Wicheer J., Leelapat P., Prapanthadara L. A., et al. (2011). The role of the Aedes aegypti Epsilon glutathione transferases in conferring resistance to DDT and pyrethroid insecticides. Insect Biochem. Mol. Biol. 41, 203–209. doi:10.1016/j.ibmb.2010.12.005
Lv Y., Wang W., Hong S., Lei Z., Fang F., Guo Q., et al. (2016). Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq. Mol. Genet. Genomics 291, 309–321. doi:10.1007/s00438-015-1109-4
Marcombe S., Doeurk B., Thammavong P., Veseli T., Heafield C., Mills M. A., et al. (2024). Metabolic resistance and not voltage-gated sodium channel gene mutation is associated with pyrethroid resistance of Aedes albopictus (skuse, 1894) from Cambodia. Insects 15, 358. doi:10.3390/insects15050358
Messenger L. A., Impoinvil L. M., Derilus D., Yewhalaw D., Irish S., Lenhart A. (2021). A whole transcriptomic approach provides novel insights into the molecular basis of organophosphate and pyrethroid resistance in Anopheles arabiensis from Ethiopia. Insect Biochem. Mol. Biol. 139, 103655. doi:10.1016/j.ibmb.2021.103655
Mourya D. T., Hemingway J., Leake C. J. (1993). Changes in enzyme titres with age in four geographical strains of Aedes aegypti and their association with insecticide resistance. Med. Vet. Entomol. 7, 11–16. doi:10.1111/j.1365-2915.1993.tb00645.x
Moyes C. L., Vontas J., Martins A. J., Ng L. C., Koou S. Y., Dusfour I., et al. (2017). Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLoS Negl. Trop. Dis. 11, e0005625. doi:10.1371/journal.pntd.0005625
Nwnp N., Silva Wapp De, Weeraratne T. C., Parakrama Karunaratne S. H. P. (2021). Kdr type mutations and enhanced GST based insecticide resistance in dengue vector mosquitoes Aedes aegypti and Aedes albopictus. Int. J. Trop. Insect Sci. 41, 409–417. doi:10.1007/s42690-020-00219-3
Paixão E. S., Teixeira M. G., Rodrigues L. C. (2018). Zika, chikungunya and dengue: the causes and threats of new and re-emerging arboviral diseases. BMJ Glob. Health 3 (Suppl. 1), e000530. doi:10.1136/bmjgh-2017-000530
Palatini U., Masri R. A., Cosme L. V., Koren S., Thibaud-Nissen F., Biedler J. K., et al. (2020). Improved reference genome of the arboviral vector Aedes albopictus. Genome Biol. 21, 215. doi:10.1186/s13059-020-02141-w
Pavlidi N., Vontas J., Van Leeuwen T. (2018). The role of glutathione S-transferases (GSTs) in insecticide resistance in crop pests and disease vectors. Curr. Opin. Insect Sci. 27, 97–102. doi:10.1016/j.cois.2018.04.007
Rahman R. U., Souza B., Uddin I., Carrara L., Brito L. P., Costa M. M., et al. (2021). Insecticide resistance and underlying targets-site and metabolic mechanisms in Aedes aegypti and Aedes albopictus from Lahore, Pakistan. Sci. Rep. 11, 4555. doi:10.1038/s41598-021-83465-w
Rault L. C., Johnson E. J., O'Neal S. T., Chen R., McComic S. E., Swale D. R., et al. (2019). Age- and sex-related ABC transporter expression in pyrethroid-susceptible and -resistant Aedes aegypti. Sci. Rep. 9, 19551. doi:10.1038/s41598-019-56134-2
Saavedra-Rodriguez K., Campbell C. L., Lenhart A., Penilla P., Lozano-Fuentes S., Black W. C. (2019). Exome-wide association of deltamethrin resistance in Aedes aegypti from Mexico. Insect Mol. Biol. 28, 591–604. doi:10.1111/imb.12575
Saavedra-Rodriguez K., Suarez A. F., Salas I. F., Strode C., Ranson H., Hemingway J., et al. (2012). Transcription of detoxification genes after permethrin selection in the mosquito Aedes aegypti. Insect Mol. Biol. 21, 61–77. doi:10.1111/j.1365-2583.2011.01113.x
Safi N. H. Z., Ahmadi A. A., Nahzat S., Warusavithana S., Safi N., Valadan R., et al. (2019). Status of insecticide resistance and its biochemical and molecular mechanisms in Anopheles stephensi (Diptera: Culicidae) from Afghanistan. Malar. J. 18, 249. doi:10.1186/s12936-019-2884-x
Schluep S. M., Buckner E. A. (2021). Metabolic resistance in permethrin-resistant Florida Aedes aegypti (Diptera: Culicidae). Insects 12, 866. doi:10.3390/insects12100866
Seixas G., Grigoraki L., Weetman D., Vicente J. L., Silva A. C., Pinto J., et al. (2017). Insecticide resistance is mediated by multiple mechanisms in recently introduced Aedes aegypti from Madeira Island (Portugal). PLoS Negl. Trop. Dis. 11, e0005799. doi:10.1371/journal.pntd.0005799
Simma E. A., Dermauw W., Balabanidou V., Snoeck S., Bryon A., Clark R. M., et al. (2019). 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. Pest Manag. Sci. 75, 1808–1818. doi:10.1002/ps.5374
Stevenson B. J., Pignatelli P., Nikou D., Paine M. J. (2012). Pinpointing P450s associated with pyrethroid metabolism in the dengue vector, Aedes aegypti: developing new tools to combat insecticide resistance. PLoS Negl. Trop. Dis. 6, e1595. doi:10.1371/journal.pntd.0001595
Strode C., Wondji C. S., David J. P., Hawkes N. J., Lumjuan N., Nelson D. R., et al. (2008). Genomic analysis of detoxification genes in the mosquito Aedes aegypti. Insect Biochem. Mol. Biol. 38, 113–123. doi:10.1016/j.ibmb.2007.09.007
Su X., Guo Y., Deng J., Xu J., Zhou G., Zhou T., et al. (2019). Fast emerging insecticide resistance in Aedes albopictus in Guangzhou, China: alarm to the dengue epidemic. PLoS Negl. Trop. Dis. 13, e0007665. doi:10.1371/journal.pntd.0007665
Wei Y., Zheng X., He S., Xin X., Zhang J., Hu K., et al. (2021). Insecticide susceptibility status and knockdown resistance (kdr) mutation in Aedes albopictus in China. Parasit. Vectors 14, 609. doi:10.1186/s13071-021-05095-5
Whalon M. E., Mota-Sanchez D., Hollingworth R. M. (2015). Arthropod pesticide resistance database. Available at: http://www.pesticideresistance.org (Accessed August 1, 2024).
WHO (2013). Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Second edition. Geneva (Switzerland): World Health Organization.
Wu Y., Wang J., Liu Q., Li T., Luo M., Gong Z. (2022). Practice of integrated vector surveillance of arthropod vectors, pathogens and reservoir hosts to monitor the occurrence of tropical vector-borne diseases in 2020 in Zhejiang Province, China. Front. Vet. Sci. 9, 1003550. doi:10.3389/fvets.2022.1003550
Yuan H., Shan W., Zhang Y., Yan H., Li Y., Zhou Q., et al. (2023). High frequency of Voltage gated sodium channel (VGSC) gene mutations in Aedes albopictus (Diptera: Culicidae) suggest rapid insecticide resistance evolution in Shanghai, China. PLoS Negl. Trop. Dis. 17, e0011399. doi:10.1371/journal.pntd.0011399
Zhai X. D., Wang S. H., Ma M., Pan D., Wang J. J., Wei D. (2024). Suppressing the expression of glutathione S-transferase gene GSTd10 increases the sensitivity of Zeugodacus cucurbitae against β-cypermethrin. Insect Mol. Biol. 33, 218–227. doi:10.1111/imb.12892
Zhao C., Zhou X., Xue C., Lun X., Li W., Liu X., et al. (2023). Knockdown resistance mutations distribution and characteristics of Aedes albopictus field populations within eleven dengue local epidemic provinces in China. Front. Cell Infect. Microbiol. 12, 981702. doi:10.3389/fcimb.2022.981702
Zhao Q., Jia Y., Lu X., Liu Y., Yin Z., Zhang Y., et al. (2024). Insecticide susceptibility and kdr mutations in Aedes albopictus collected from seven districts of Guangyuan city, northern Sichuan, China. J. Am. Mosq. Control Assoc. 40, 20–25. doi:10.2987/23-7155
Zhou X., Yang C., Liu N., Li M., Tong Y., Zeng X., et al. (2019a). Knockdown resistance (kdr) mutations within seventeen field populations of Aedes albopictus from Beijing China: first report of a novel V1016G mutation and evolutionary origins of kdr haplotypes. Parasit. Vectors 12, 180. doi:10.1186/s13071-019-3423-x
Zhou Y., Fu W. B., Si F. L., Yan Z. T., Zhang Y. J., He Q. Y., et al. (2019b). UDP-glycosyltransferase genes and their association and mutations associated with pyrethroid resistance in Anopheles sinensis (Diptera: Culicidae). Malar. J. 18, 62. doi:10.1186/s12936-019-2705-2
Zhu G., Zhong D., Cao J., Zhou H., Li J., Liu Y., et al. (2014). Transcriptome profiling of pyrethroid resistant and susceptible mosquitoes in the malaria vector, Anopheles sinensis. BMC Genomics 15, 448. doi:10.1186/1471-2164-15-448
Zou F., Guo Q., Shen B., Lv Y. (2024). Pyrethroid susceptibility status and functional analysis of cytochrome P450 CYP6 subfamily genes in field Aedes albopictus. Jiangsu, China. J. Asia-Pac Entomol. 27, 102210. doi:10.1016/j.aspen.2024.102210
Keywords: transcriptome, insecticide resistance, cytochrome P450, metabolic resistance, Aedes albopictus
Citation: Zhou X, Li J, Ni R, Qiu X, Zhang Y and Tong Y (2024) Insecticide resistance in the field populations of the Asian tiger mosquito Aedes albopictus in Beijing: resistance status and associated detoxification genes. Front. Physiol. 15:1498313. doi: 10.3389/fphys.2024.1498313
Received: 26 September 2024; Accepted: 14 November 2024;
Published: 18 December 2024.
Edited by:
Paul O. Mireji, Kenya Agricultural and Livestock Research Organization, KenyaReviewed by:
Yong Huang, Anhui Agricultural University, ChinaJames Mutunga, Penn State Harrisburg, United States
Copyright © 2024 Zhou, Li, Ni, Qiu, Zhang and Tong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Zhang, emhhbmd5Y2RjQDEyNi5jb20=; Ying Tong, dG9uZ3lpbmc5NkAxMjYuY29t
†These authors have contributed equally to this work and share first authorship