Skip to main content
Download PDF

Insecticide resistance of Anopheles sinensis after elimination of malaria in Henan Province, China

Parasites & Vectors volume 16, Article number: 180 (2023) Cite this article

Abstract

Background

Historically, malaria due to Plasmodium vivax has been epidemic in Henan Province, China, with Anopheles sinensis as the main vector. The most effective measures to prevent malaria transmission are based on vector control through the use of insecticides. However, insecticides exert a strong selective pressure on mosquito populations for insecticide resistance. The aim of this study was to investigate the susceptibility profile and population genetic characteristics of An. sinensis to provide basic data and scientific guidance for the study of resistance mechanisms and the control of An. sinensis in Henan Province.

Methods

Adult Anopheles mosquitoes were collected at sites near local farmers' sheepfolds, pigsties and/or cowsheds located in Pingqiao, Xiangfu, Xiangcheng and Tanghe counties/districts of Henan Province during July–September 2021 for insecticide susceptibility testing. Molecular identification of collected mosquitoes as belonging to genus Anopheles was by PCR, and the frequencies of mutations in the knockdown resistance gene (kdr) and acetylcholinesterase-1 gene (ace-1) were detected using gene amplification. The mitochondrial DNA cytochrome oxidase subunit I (COI) gene was amplified in deltamethrin-resistant and deltamethrin-sensitive mosquitoes to analyze the genetic evolutionary relationship.

Results

A total of 1409 Anopheles mosquitoes were identified by molecular identification, of which 1334 (94.68%) were An. sinensis, 28 (1.99%) were An. yatsushiroensis, 43 (3.05%) were An. anthropophagus and four (0.28%) were An. belenrae. The 24-h mortality rates of An. sinensis in Pingqiao, Tanghe, Xiangcheng and Xiangfu counties/districts exposed to deltamethrin were 85.85%, 25.38%, 29.73% and 7.66%, respectively; to beta-cyfluthrin, 36.24%, 70.91%, 34.33% and 3.28%, respectively; to propoxur, 68.39%, 80.60%, 37.62% and 9.29%, respectively; and to malathion, 97.43%, 97.67%, 99.21% and 64.23%, respectively. One mutation, G119S, was detected in the ace-1 gene. The frequencies of the main genotypes were 84.21% of specimens collected in Xiangfu (G/S), 90.63% of speciments collected in Xiangcheng (G/G) and 2.44% of speciments collected in Tanghe (S/S). Significantly higher G119S allele frequencies were observed in both propoxur- and malathion-resistant mosquitoes than in their sensitive counterparts in the Tanghe population (P < 0.05). Three mutations, L1014F (41.38%), L1014C (9.15%) and L1014W (0.12%), were detected in the kdr gene. The genotypes with the highest frequency in the populations of An. sinensis in Xiangfu and Tanghe were the mutant TTT (F/F) and wild-type TTG (L/L), at 67.86% (57/84) and 74.29% (52/70), respectively. In Pingqiao and Xiangfu, higher frequencies of the L1014F allele and lower frequencies of the L1014C allele were observed in mosquitoes resistant for beta-cyfluthrin than in those which were sensitive for this insecticide (P < 0.05). The results of Tajima's D and of Fu and Li's D and F were not significantly negative (P > 0.10), and each haplotype was interlaced and did not form two distinct branches.

Conclusions

High resistance to pyrethroids and propoxur was observed at four sites, but the resistance to malathion varied according to the location. Anopheles belenrae and the L1014W (TGG) mutation in An. sinensis were first discovered in Henan Province. The deltamethrin-resistant and deltamethrin-sensitive mosquito populations showed no genetic differentiation. The generation of resistance might be the result of the combination of multiple factors.

Graphical Abstract

Background

Malaria, caused by infection with Plasmodium spp., is one of the most significant life-threatening infectious diseases in humans worldwide. According to the WHO Malaria Report 2022, the total number of deaths due to malaria worldwide reached 619,000 in 2021, with approximately 247 million reported malaria cases [1]. Anopheles sinensis is the predominant Anopheles species in Henan Province, China and the main vector of the malaria parasite Plasmodium vivax [2]. In 2006, there was an outbreak of malaria caused by P. vivax in Yongcheng county, Henan Province [3]. Malaria vectors are known to still exist in the original malaria-endemic areas and if there are also imported malaria cases, the risk of a malaria outbreak is extremely high. It has become an important measure to control the population density of An. sinensis and other malaria vectors and prevent secondary transmission caused by imported cases [4, 5]. Indigenous malaria transmission has been effectively controlled, and China officially obtained the WHO certification for malaria elimination in 2021. With the development of global trade and the transnational economy, infections imported from overseas have become the main source of malaria cases in Henan Province. Imported malaria has introduced new challenges to achieving the overall elimination of malaria [6,7,8].

Insecticides have been the most important measure to control the number of mosquitoes. Because of their broad-spectrum, efficiency and durable characteristics, insecticides are widely used across the globe [9]. Five major categories of insecticides have been applied: pyrethroids, organochlorines, carbamates, neonicotinoids and organophosphates [10]. Pyrethroids are emerging as the predominant insecticides for vector control because of their short residual action, high efficacy against mosquito vectors and low toxicity to humans [11]. The main targets of insecticides on mosquito vectors are the sodium ion channels (SC), acetylcholinesterase (AchE) and the gamma-aminobutyric (GABA) receptor-chloride channel complex [12]. Among these, mutations in the genes associated with the sodium channel in the axon produce knockdown resistance (kdr). The mutation at codon 1014 of the knockdown gene (kdr) in the IIS6 region of the sodium channel is also known as the kdr mutation [13, 14]. Four types of kdr mutations have been identified: (i) the L1014C (TTG → TGT) mutation leads to a change from leucine to cysteine substitution; (ii) the L1014W (TTG → TGG) mutation leads to a leucine to tryptophan substitution; and (iii) two L1014F (TTG → TTT, and TTG → TTC) mutations lead to a leucine to phenylalanine substitution. A mutation at codon 119 of the acetylcholinesterase-1 gene (ace-1) leads to a single amino acid substitution of glycine to serine in the binding pocket of acetylcholinesterase, conferring resistance to organophosphates and carbamates [15].

Understanding the genetic differentiation of the population helps to infer the current situation and possible diffusion path of insecticide resistance, thus guiding the use of insecticides [16]. Population genetics, based on the characteristics of population genetic structure, analyzes the population evolutionary process, studies the genetic structure and change dynamics of organisms and explains the differences between populations. As an important component of the mitochondrial gene, the cytochrome oxidase subunit I (COI) gene has rich variation and is often used for phylogenetic and population genetic diversity analysis. The mitochondrial DNA (mtDNA)-COI gene sequence has been used in population genetic analyses due to its simple structure, fast evolutionary rate and extremely low probability of recombination [17,18,19].

Although the insecticide resistance of An. sinensis in Henan Province has been studied, little resistance monitoring data from recent years are available [20,21,22]. Therefore, using previous monitoring data, we have studied the insecticide resistance of An. sinensis in four regions of Henan Province, with the aim to determine the insecticide resistance status of An. sinensis in Henan Province. To further understand the current situation of target resistance and population genetic characteristics of An. sinensis in Henan Province, we selected samples of An. sinensis to detect the mutation and frequency of the kdr and ace-1 genes, and carried out population genetic analysis of sensitive and resistant mosquitoes based on the mtDNA-COI gene, providing a scientific basis for guiding the control measures required for An. sinensis and basic data for further study of the population genetic characteristics of An. sinensis in China.

Methods

Study sites and mosquito collection

This study on the insecticide resistance of mosquitoes was carried out during July–September 2021 in four districts/counties of Henan Province, namely Pingqiao, Xiangfu, Xiangcheng and Tanghe. Adult Anopheles mosquitoes were collected using a Maxttrac Mosquito Trap (Shanghai Daolan Industrial Co., Ltd, China) near local farmers' sheepfolds, pigsties and/or cowsheds. Collected mosquitoes were transported to the insectary of the Henan Provincial Center for Disease Control and Prevention (CDC) where they were acclimatized under standard conditions (27 ± 1 °C, and 80% relative humidity [RH]) and provided with 8% glucose solution until the following day when they were used in the experiments. All collected mosquitoes were morphologically identified to their respective complexes and groups [23].

Insecticide susceptibility test

In accordance with the standard WHO insecticide susceptibility tube-test procedures [24], adult female mosquitoes were used to conduct bioassays to determine susceptibility to four insecticides (supplied by China CDC): malathion (5%), deltamethrin (0.05%), beta-cyfluthrin (0.15%) and propoxur (0.1%). As controls, papers treated with olive oil, silicone oil, silicone oil and olive oil, respectively, were used. For each insecticide, 120 female mosquitoes were tested, with 20 mosquitoes per tube and two control tubes without insecticide. The number of knockdown mosquitoes was recorded during the exposure period at 10, 15, 20, 30, 40, 50 and 60 min. After 60 min of exposure, all mosquitoes were transferred to holding tubes, provided with 8% sugar water and placed under conditions of 26 °C and 80% RH. At the end of this 24-h recovery period, the number of dead mosquitoes was counted and recorded. After the bioassays, all mosquitoes were preserved individually in 1.5-ml Eppendorf tubes for further molecular analyses.

DNA extraction and molecular identification

The mosquitoes collected in the field and morphologically identified as An. sinensis were tested for insecticide resistance and then their identity molecularly confirmed. One leg of each mosquito was used for DNA extraction with the Animal Tissue Genomic DNA Extraction Kit (Xi’an Tianlong Science & Technology Co., Ltd, Xi’an City, China). More precisely, one leg was placed in 100-μl aliquot of tissue lysis solution and fully ground using an electric tissue grinder. Then, 200 μl of tissue lysis solution and 40 µl of protein K solution were added and mixed by vortexing. An automatic nucleic acid extractor (Qiagen, Hilden, Germany) was used for automatic DNA extraction. Molecular identification of An. sinensis was performed by using species-specific primers and amplification of internal transcribed spacer 2 (ITS2) [25]. The samples identified as An. sinensis were again used for subsequent gene amplification and data analysis.

Detection of ace-1 and kdr mutations

To determine point mutations of the ace-1 gene at codon 119, a 193-bp fragment was amplified in An. sinensis using primers ace-1-F (467F; 5’-GTGCGACCATGTGGAACC-3’) and ace-1-R (660R; 5’-ACCACGATCACGTTCTCCTC-3’) [26]. Amplification reactions were performed in a 25-µl reaction volume containing 2 µl of the DNA sample, 12 µl of a 2 × Go Taq Green Master Mix (Promega Inc., Madison, WI, USA), 2 µl of target primers and 9 µl of ddH2O. PCR cycling consisted of 1 cycle at 95 °C for 3 min, followed by 35 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 45 s, with a final extension at 72 °C for 6 min. Sequencing was performed by Shanghai DNA Bio Technologies Co., Ltd. (Shanghai, China). Genotyping of point mutations of the kdr gene at codon 1014 was conducted as previously described [11].

COI gene amplification of deltamethrin-resistant and -sensitive mosquitoes

To clarify the genetic evolutionary relationship between resistant and sensitive mosquitoes of An. sinensis, the mtDNA-COI gene was amplified in deltamethrin-resistant and deltamethrin-sensitive An. sinensis using the forward primer GGTCAACAAATCATAAAGATATTGG and the reverse primer TAAACTTCAGGGTGACCAAAAAATCA, following the methods previously described by Liang et al. [27]. All PCR products were analyzed by 2% agarose gel electrophoresis and the purified products were sequenced using an ABI 3730XL automated sequencer (Applied Biosystems®, Thermo Fisher Scientific, Waltham, MA, USA.

Statistical analysis

The mortality rates of the test and control samples were calculated after the bioassay for each insecticide group. If the control mortality was > 20%, the tests were discarded. When control mortality was 5% and 20%, then the observed mortality was corrected using Abbott’s formula [28]. Mosquito resistance status was interpreted in accordance with WHO guidelines as: (i) confirmed resistance (CR), mortality rate < 90%; (ii) possible resistance (PR), mortality rate between 90% and 97%; and (iii) susceptibility (S), mortality rate ≥ 98% [24]. Statistical analyses were performed using SPSS version 21.0 software (IBM Corp., Armonk, NY, USA). The mortality rates of An. sinensis at the survey sites, ace-1 and kdr mutation rates and gene mutation rates of the different populations were calculated using the Chi-square (χ2) test. Nucleotide polymorphism, natural selection and population differentiation were described according to He et al. [29]. In addition, we constructed a haplotype network diagram among the different geographical populations based on the median connection method, with the aim to analyze the genetic relationship among different geographical populations using PopArt software. A P value of < 0.05 was considered to indicate significance.

Results

Molecular identification of Anopheles species and insecticide susceptibility test

The molecular identification of all samples used for insecticide resistance determination was carried out based on the length of the electrophoresis fragments of the PCR product. A total of 1409 Anopheles mosquitoes were identified by species-specific PCR, of which 1334 (94.68%) were An. sinensis, 28 (1.99%) were An. yatsushiroensis, 43 (3.05%) were An. anthropophagus and four (0.28%) were An. belenrae (Table 1). Analysis of the insecticide resistance of the 1334 An. sinensis mosquitoes showed that the 24-h mortality rates in mosquitoes exposed to deltamethrin, beta-cyfluthrin and propoxur ranged from 3.28% to approximately 85.85%, which reached the CR level. The 24-h mortality rate of An. sinensis to malathion in Pingqiao and Tanghe was 97.43% and 97.67%, respectively, which reached the PR level. The 24-h mortality rates for malathion in Xiangfu and Xiangcheng were 64.23% (CR) and 99.21% (S), respectively. The difference between the knockdown rate and 24-h mortality of the different insecticides at the four sampling locations was statistically significant (χ2 = 32.635 to approx. 143.433; P < 0.05). For beta-cyfluthrin, deltamethrin, propoxur and malathion, the median knockdown time for 50% of the tested mosquitoes (KT50) ranged from 29.81 min to approximately 315.73 min in Pingqiao, Tanghe, Xiangcheng and Xiangfu. For beta-cyfluthrin and deltamethrin, the KT50 could not be calculated due to the low knockdown rate in Xiangfu (Table 2).

Table 1 Sampling locations and molecular identification of Anopheles species collected in Henan Province, China
Full size table
Table 2 Adjusted mortality rates and knockdown times of Anopheles sinensis exposed to four insecticides in Henan Province
Full size table

Association between ace-1 genotype and insecticide resistance

Mosquitoes molecularly identified as An. sinensis were used for amplification of ace-1. A total of 315 mosquitoes, including propoxur- and malathion-resistant and -sensitive mosquitoes, were used in this study. After splicing the sequencing results, a 193-bp-long gene sequence was obtained. Three genotypes of the ace-1 gene in different populations of An. sinensis were detected through sequence alignment: GGC (G/G), AGC (S/S) and GGC/AGC (G/S). Among these, G/G is the homozygous wild type, S/S is a homozygous mutant and G/S is a heterozygous mutant. Three genotypes were detected in three populations: Pingqiao, Tanghe and Xiangfu. Only two genotypes were detected in Xiangcheng, with no mosquito with genotype S/S found at this location. The main genotype of the Pingqiao and Xiangfu populations was G/S, with frequencies of 68.67% (57/83) and 84.21% (80/95), respectively. The Tanghe and Xiangcheng populations were genotype-based, with G/G having the highest frequency in both populations at 73.17% (30/41) and 90.63% (87/96), respectively, and genotype S/S only found in the Tanghe population, with a frequency of only 2.44% (1/41). Two alleles were found in the four sampling populations of An. sinensis in Henan Province: wild-type G119G and mutant-type G119S, with a total mutation rate of 33.65%. Significantly higher G119S allele frequencies were observed in both propoxur- and malathion-resistant mosquitoes than in propoxur- and malathion-sensitive mosquitoes in the Tanghe population (P < 0.05), but this distribution was not observed in the other three populations. Table 3 presents detailed data on the relationship between the ace-1 mutant gene and insecticide resistance in An. sinensis.

Table 3 Relationship between the ace-1 mutant gene and insecticide resistance in An. sinensis
Full size table

Association between kdr genotype and insecticide resistance

A 325-bp fragment that included position 1014 was successfully sequenced from 339 deltamethrin and beta-cyfluthrin-resistant and deltamethrin and beta-cyfluthrin-sensitive An. sinensis. The mutation frequency at codon 1014 of the kdr gene at the four collection sites was 68.83% (Pingqiao), 18.04% (Tanghe), 17.70% (Xiangcheng) and 98.05% (Xiangfu). Ten genotypes of the kdr gene in the different populations of An. sinensis were detected through sequence alignment: TTG (L/L), TTT (F/F), TGT (C/C), TTG/TGT (L/C), TGG/TTT (W/F), TTT/TGT (F/C), TGT/TTC (C/F), TTG/TTT (L/F), TTG/TTC (L/F) and TTT/TTC (F/F). Among these, TTG is the wild type; the other nine are mutations, with TTT and TGT being homozygous mutations, and the remaining seven genotypes being heterozygous mutations. The highest number of genotypes was detected in the Pingqiao sample, including nine genotypes. The genotypes with the highest frequency in the Pingqiao and Xiangfu populations of An. sinensis were the mutant TTT (F/F) genotype, at 32.18% (28/87) and 67.86% (57/84), respectively. The highest frequency of the wild-type TTG (L/L) genotype was found in the Tanghe and Xiangcheng populations, at 74.29% (52/70) and 67.35% (66/98), respectively (Additional file 1: Table S1). As shown in Table 4, the total mutation rate of the kdr gene was 50.65%, with the L1014F mutation having the highest frequency (41.38%); the frequency of the L1014C mutation was 9.15%, and the L1014W mutation had the lowest frequency (0.12%). A significantly higher frequency of the L1014F allele was observed in beta-cyfluthrin-resistant mosquitoes than in their sensitive counterparts in the Pingqiao, Tanghe and Xiangfu populations, but a significantly lower frequency of the L1014C allele was observed in beta-cyfluthrin-resistant mosquitoes than in their sensitive counterparts in the Pingqiao and Xiangfu populations (P < 0.05).

Table 4 Relationship between the kdr mutant gene and insecticide resistance in An. sinensis
Full size table

Genetic diversity analysis of the COI gene of deltamethrin-resistant and deltamethrin-sensitive mosquitoes

Genetic diversity analysis of 58 COI gene sequences of mosquitoes collected in the 2020–2021 study period showed that there were 59 mutation sites, including 27 singleton variable sites and 32 parsimony informative sites. As shown in Table 5, the total haploid diversity (Hd ± standard deviation) was 0.953 ± 0.022 according to the analysis of the DNA diversity index in 2021, and the haploid diversity index (Hd) ± SD was 1.00 in 2020. The total nucleotide diversity was 0.00910 (in 2021) and 0.00864 (in 2020). The results of Tajima's D and of Fu and Li's D and F were not significantly negative from 2020–2021 (P > 0.10). The genetic distance between resistant and sensitive mosquito populations was 0.009, but their intraspecific genetic distances were 0.009 and 0.01, indicating a close genetic relationship between the two populations. The genetic differentiation index (Fst) between the two populations was not significantly negative (P > 0.10), indicating no significant genetic differentiation between the two populations (Additional file 2: Table S2).

Table 5 Genetic diversity analysis of mosquitoes collected in Xiangcheng from 2020 to 2021
Full size table

A total of 36 haplotypes were detected in 52 mtDNA-COI gene sequences of mosquitoes collected in 2021, including 32 unique haplotypes and four shared haplotypes. Among the 32 unique haplotypes, 18 were from resistant mosquitoes and 14 were from sensitive mosquitoes. The haplotype network diagram shows that shared haplotype Hap_8 accounted for the largest number (6 resistant mosquitoes, 5 sensitive mosquitoes) of mosquitoes, and haplotypes Hap_6, Hap_12 and Hap_ 21 were shared by two populations, but there were only two. As shown in Fig. 1, each haplotype was interlaced and did not form two distinct branches.

Fig. 1
figure 1

Haplotype network of deltamethrin-resistant and deltamethrin-sensitive mosquitoes. Circles represent a haplotype, black dots represent the assumed intermediate haplotype and the horizontal line represents a mutation step. The size of the circle is proportional to the frequency of the haplotype

Full size image

Discussion

Molecular monitoring and bioassays of insecticide resistance status in local malaria vectors are important tools for managing and controlling resistance vectors [30]. A previous study demonstrated the initial development of resistance in An. sinensis to deltamethrin and cyfluthrin in Nanyang City [31]. The results of Hu et al.’s study indicate that resistance of An. sinensis to deltamethrin and beta-cyfluthrin in Puyang City occurred in 2018 and 2020 [20]. The results of resistance monitoring in our study showed that the resistance of An. sinensis to pyrethroid insecticides in Henan Province is increasing, with the level of resistance to the propoxur insecticide already reached in the four studied regions; however, the resistance to malathion was different in terms of sensitivity to resistance. The results for malathion are similar to those reported in previous studies [20, 21], possibly due to the lower use of malathion in Henan Province in recent decades [32]. Given this present-day sensitivity of mosquitoes to malathion, alternative insecticides can be used to prolong the development of mosquito resistance to malathion.

In our study, only one mutation type, G119S, was found. This result is in line with observations reported in published studies suggesting that the G119S mutation is widely distributed in An. sinensis [33]. Compared to other mosquitoes, this result is also consistent with the ace gene mutations in Anopheles gambiae and Culex pipiens [34]. Qin et al. reported a modest to high (45–75%) ace-1 mutation frequency in An. sinensis populations [33]. Chang et al. demonstrated that the frequencies of G119S mutations in the Anhui (China) and Yunnan (China) populations were 58.9% and 38.5%, respectively [35]. These results suggest that the G119S mutation rate of An. sinensis is high and that this mutation is widely distributed. A previous study revealed that the G119S allele frequency of a homozygous mutation was 100% in mosquitoes from the two An. sinensis populations, indicating that it is not involved in malathion resistance [36]. The gene mutation rate of the Xiangfu population in Henan Province was as high as 54.74%, indicating that the high frequency of the G119S mutation might be the result of the long-term use of organophosphorus insecticides for agricultural pest control [27]. The high frequency of the G119S mutation indicated that these An. sinensis populations have a high risk of developing resistance to organophosphorus and carbamate insecticides. Other researchers found that the G285A and F350Y mutations on the ace gene are related to resistance to organic phosphorus, and that the F105S, G285A and F305Y mutations are related to resistance to carbamate [37]. Higher G119S allele frequencies were observed in both propoxur- and malathion-resistant mosquitoes than in propoxur- and malathion-sensitive mosquitoes in the Tanghe population, but such a development was not observed in the other three sites, suggesting that the production of resistance might be the result of the joint action of multiple factors.

Kdr was first identified in the house-fly Musca domestica L. [38]. The L1014F substitution has been reported in pyrethroid-resistant pest species, including An. gambiae and Cx. p. pallens [39, 40]. The L1014C mutation, with TTG (Leu) being replaced with TGT (Cys), is a newly reported mutation [41]. The L1014F and L1014C mutations were detected in the Republic of Korea, where the TTC L1014F mutation was observed for the first time [42]. Previous studies revealed that the frequencies of the L1014F allele were significantly associated with deltamethrin-resistant and DDT-resistant phenotypes, but that this was not the case for the L1014C allele. In Kaifeng (Henan, China), the difference in kdr mutation frequency between surviving and dead mosquitoes was statistically significant [36, 43]. However, Sun et al. reported that the difference in kdr mutation frequency between surviving and dead mosquitoes was not statistically significant in Yingjiang (Yunnan, China) and Suining (Jiangsu, China) [43].

Three types of kdr mutant alleles were found in the present study, including L1014F, L1014C and L1014W. In contrast to the situation with An. sinensis in Sichuan, China, where the frequency of L1014F was observed to be higher than that of L1014L. L1014L was found to be the predominant resistance allele in Henan, but L1014C had lower frequencies than that L1014F and L1014L, which is similar to the results of the present study [33]. L1014F and L1014C mutations were found in Henan Province in a previous study, and the L1014W mutation was first found in the present study [11]. In the present study, the proportion of the L1014W mutation was very small (0.12%) and only distributed in the Pingqiao population. Tan et al. reported that the presence of the L1014W mutation in 52 An. sinensis in Guangxi, based on their kdr gene mutation study, accounting for a small proportion (1%) [44]. This result is similar to that of the current study. Higher L1014F and lower L1014C allele frequencies were observed in beta-cyfluthrin-resistant mosquitoes than in beta-cyfluthrin-sensitive mosquitoes, suggesting that the relationship between higher L1014F and lower L1014C allele frequencies and insecticide resistance requires further in-depth analysis.

The evolution of insect resistance is based on the genetics of insecticide-resistant populations. The Hd was > 0.95, and nucleotide diversity was > 0.008 from 2020–2021, indicating that the overall genetic diversity of An. sinensis in Henan Province is high. The values of Tajima's D and of Fu and Li's D and F suggest that the populations conform to the neutral selection hypothesis and are not subject to obvious selection pressure during the evolution process. Yang et al. studied the population genetic variation and population structure characteristics of five species groups of An. sinensis in nine sampling sites in Yunnan Province [45]. These authors reported gene exchange between the four species groups of An. sinensis in Yunnan Province, except for the YU population group, and that genetic differentiation was not apparent [45]. The values of the Fst and the haplotype network diagram indicated that there was no genetic differentiation between deltamethrin-resistant and deltamethrin-sensitive mosquito populations, possibly because the mosquitoes in this study were all adult mosquitoes collected in the field, and they lived in the same ecological environment and faced the same pressure of insecticide selection. External environmental factors might neutralize the degree of genetic differentiation, so genetic differentiation between the two populations was not observed. The results of this study show that the emergence of mosquito resistance might be the result of multiple mechanisms.

The advantages of adult female mosquitoes directly captured in the field for the study of insecticide resistance are that fewer facilities are required, and the age distribution of the vectors is representative of the wild vector population at a given time and location. The limitations of this study include failing to accurately assess adult female mosquitoes’ blood-sucking history, age and survival status, any of which might affect the test results. No further study on mutations at other sites besides the kdr and ace-1 genes and no comprehensive study of resistance mechanisms have been carried out from the perspective of enzymology and proteomics.

Conclusions

Four Anopheles species were identified in the present study: An. sinensis, An. yatsushiroensis, An. anthropophagus and An. belenrae. Anopheles belenrae was discovered for the first time in Henan Province. A high resistance to pyrethroids and propoxur was noted in the mosquitoes collected in Henan Province, but resistance to malathion varied according to the sampling location. The L1014W (TGG) mutation was detected for the first time in this study. Only one mutation, G119S, was found in the ace-1 gene. Mutations in the kdr and ace-1 genes play a role in the generation of insecticide resistance in An. sinensis. The genetic diversity of deltamethrin-resistant and deltamethrin-sensitive mosquitoes showed no genetic differentiation. The generation of resistance might be the result of the combination of multiple factors. Further experiments are needed to collect comprehensive research data on resistance mechanisms from the perspectives of enzymology and proteomics. These findings are a step towards providing missing data on the insecticide resistance of An. sinensis in Henan Province and provide the scientific basis for guiding the control of An. sinensis.

Availability of data and materials

The datasets used during the current study are available from the corresponding author on reasonable request.

Abbreviations

ace-1 :

Acetylcholinesterase-1

CDC:

Center for Disease Control and Prevention

COI:

Cytochrome oxidase subunit I

CR:

Confirmed resistance

GABA:

Gamma-aminobutyric acid

kdr :

Knockdown resistance

mtDNA:

Mitochondria DNA

PR:

Possible resistance

S:

Susceptible

SC:

Sodium ion channels

References

  1. WHO. World malaria report 2022. 2022. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2022. Accessed 7 Feb 2023.

  2. Zhang HW, Su YP, Xu BL. A review of malaria control and research in Henan Province of China. J Pathog Biol. 2006;1:64–6.

    CAS  Google Scholar 

  3. Zhang HW, Su YP, Zhou GC, Liu Y, Xu BL. Epidemiological survey of malaria in Shangqiu City, Henan Province in 2000–2005. Chin Trop Dis. 2006;7:1157–8.

    Google Scholar 

  4. Chu WD, Zhang YJ. Imported malaria epidemic situation and it’s surveillance in China. J Trop Med. 2018;18:127–30.

    Google Scholar 

  5. Chen HY, Li KL, Wang XH, Yang XY, Lin Y, Cai F, et al. First identification of kdr allele F1534S in VGSC gene and its association with resistance to pyrethroid insecticides in Aedes albopictus populations from Haikou City, Hainan Island, China. Infect Dis Poverty. 2016;5:31–8.

    Article  PubMed  PubMed Central  Google Scholar 

  6. National Health Commission of the People’s Republic of China. China successfully passed the World Health Organization certification for malaria elimination. 2021. http://www.nhc.gov.cn/jkj/s5874/202106/cffeb8ba568c4939a7e72668dd446d4c.shtml. Accessed 8 Feb 2023.

  7. Feng J, Zhang L, Tu H, Zhou SS, Xia ZG. From elimination to post elimination: epidemic characteristics, challenges and strategies to prevent re transmission of imported malaria in China. Chin Trop Dis. 2021;21:5–10.

    Google Scholar 

  8. Yang CY, Qian D, Lu DL, Liu Y, Zhou RM, Li SH, et al. Epidemic status of malaria and progress of malaria elimination in Henan Province, 2018. Chin J Schisto Control. 2020;32:298–300.

    CAS  Google Scholar 

  9. Jiang L, Zhou YM, Gao YF. Major advance in mosquito resistance research. Port Health Control. 2019;24:59–62.

    Google Scholar 

  10. WHO. Prequalified vector control products. 2023. https://extranet.who.int/pqweb/vector-control-products/prequalified-product-list. Accessed 31 Mar 2023.

  11. Zhang HW, Liu Y, Hu T, Zhou RM, Chen JS, Qian D, et al. Knockdown resistance of Anopheles sinensis in Henan province, China. Malar J. 2015;14:137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Lin KM, Yang YC, Li J. Progress of insecticide resistance of the main malaria vectors in China. Chin Trop Dis. 2019;19:584–90.

    Google Scholar 

  13. Chang XL, Zhong DB, Eugenia L, Fang Q, Bonizzoni M, Wang XM, et al. Landscape genetic structure and evolutionary genetics of insecticide resistance gene mutations in Anopheles sinensis. Parasit Vectors. 2016;9:228.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Tan WL, Li CX, Lv RC, Dong YD, Guo XX, Xing D, et al. The polymorphism and geographical distribution of knockdown resistance of adult Anopheles sinensis populations in eastern China. Malar J. 2019;18:164.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Essandoh J, Yawson A, Weetman D. Acetylcholinesterase (Ace-1) target site mutation 119S is strongly diagnostic of carbamate and organophosphate resistance in Anopheles gambiae s.s. and Anopheles coluzzii across southern Ghana. Malar J. 2013;12:404.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Yang FL. Taxonomy of the Anopheles hyrcanus group from southern China and population structure of Anopheles sinensis in Yunnan (Diptera: Culicidae: Anopheles). Chongquing: Chongqing Normal University; 2014. p. 1–64.

  17. Wang XT, Zhang YJ, He X, Mei T, Chen B. Identification, characteristics and distribution of microsatellites in the whole genome of Anopheles sinensis (Diptera: Culicidae). Acta Entomol Sin. 2016;59:1058–68.

    Google Scholar 

  18. Li J, Zhang XF, Xu LL, Shen YY, Li XX, Wang ZY. Genetic structure and Wolbachia infection in geographical populations of Monolepta hieroglyphica (Coleoptera: Chrysomelidae) in South China. Acta Entomol Sin. 2021;64:730–42.

    Google Scholar 

  19. Guo XX, Cheng P, Liu LJ, Zhang CX, Wang HF, Wang HW, et al. Analysis of population genetic diversity of mosquitoes from Shandong Province based on mitochondrial DNA cytochrome oxidase subunit Igene fragment. Chin J Schisto Control. 2018;30:37–41.

    Google Scholar 

  20. Hu YB, He ZQ, Liu Y, Qian D, Yang CY, Lu DL, et al. Susceptibility of Anopheles sinensis to insecticides in Puyang City, Henan Province. Chin J Schisto Control. 2021;33:501–4+9.

    CAS  Google Scholar 

  21. Liu Y, Chen JS, Zhou RM, Qian D, Chen QW, Xu BL, et al. Investigation on the sensitivity of Anopheles sinensis to insecticide. Chin J Parasitol Parasit Dis. 2012;30:309–11.

    CAS  Google Scholar 

  22. Zhang YQ, Zhao Q, Liu JQ, Guo XS, Tang ZQ, Zhao XD. Resistance of Anopheles sinensis to three kinds of insecticide in Xinyang city. Chin J Hyg Insect Equip. 2014;20:368–70.

    CAS  Google Scholar 

  23. Wu GL. Human parasitology. 4th ed. Beijing: People’s Health Publishing House; 2013.

    Google Scholar 

  24. WHO. Test procedures for insecticide resistance monitoring in malaria vector mosquitoes. Second edition. 2018. https://who.int/iris/bitstream/handle/10665/250677/9789241511575-eng.pdf. Accessed 31 Mar 2023.

  25. Yajun M, Xu J. The Hyrcanus group of Anopheles (Anopheles) in China (Diptera: Culicidae): species discrimination and phylogenetic relationships inferred by ribosomal DNA internal transcribed spacer 2 sequences. J Med Entomol. 2005;42:610–9.

    PubMed  Google Scholar 

  26. Chang XL, Zhong DB, Fang Q, Hartsel J, Zhou GF, Shi LN, et al. 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. 2014;8:e2889.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Liang QG. Detection of mutations in kdr and ace-1 genes of resistance targets and population genetics of Anopheles sinensis in Guizhou Province, China. Guizhou: Guizhou Medical University; 2019. p. 1–50.

  28. Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–7.

    Article  CAS  Google Scholar 

  29. He ZQ, Zhang QQ, Wang D, Hu YB, Zhou RM, Qian D, et al. Genetic polymorphism of circumsporozoite protein of Plasmodium falciparum among Chinese migrant workers returning from Africa to Henan Province. Malar J. 2022;21:248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zhong DB, Aung PL, Mya MM, Wang XM, Qin Q, Soe MT, et al. Community structure and insecticide resistance of malaria vectors in northern-central Myanmar. Parasit Vectors. 2022;15:155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Cui SL, Li F, Chen GY, Zhang HP, Wang TZ. Sensitivity of Anopheles sinensis to insecticides. Chin J Hyg Insect Equip. 2013;19:223–4.

    CAS  Google Scholar 

  32. Bo YX, Tang ZQ, Zhang YQ, Li SJ. Insecticide resistance of Culex pipiens pallens in Henan Province. Chin J Hyg Insect Equip. 2010;16:214–6.

    CAS  Google Scholar 

  33. Qian WP, Liu N, Yang Y, Liu J, He JH, Chen ZH, et al. A survey of insecticide resistance-conferring mutations in multiple targets in Anopheles sinensis populations across Sichuan, China. Parasit Vectors. 2021;14:169.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, et al. Insecticide resistance in mosquito vectors. Nature. 2003;423:136-7. https://doi.org/10.1038/423136b.

  35. Qin Q, Li YJ, Zhang DB, Zhou N, Chang XL, Li CY, et al. Insecticide resistance of Anopheles sinensis and An. vagus in Hainan Island, a malaria-endemic area of China. Parasit Vectors. 2014;7:92.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chen SX, Qin Q, Zhong DB, Fang X, He HJ, Wang LJ, et al. Insecticide resistance status and mechanisms of Anopheles sinensis (Diptera: Culicidae) in Wenzhou, an important coastal port city in China. J Med Entomol. 2019;56:803–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Vaughan A, Rocheleau T, Ffrench-Constant R. Site-directed mutagenesis of an acetylcholinesterase gene from the yellow fever mosquito Aedes aegypti confers insecticide insensitivity. Exp Parasitol. 1997;87:237–44.

    Article  CAS  PubMed  Google Scholar 

  38. Milani R. Comportamento mendeliano della resistenza alla azione jabbattante del DDT: correlazione tran abbattimento emortalia in Musca domestica L. Riv Parasitol. 1954;15:513–42.

    Google Scholar 

  39. Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, et al. Molecular characterization of pyrethoid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s.. Insect Mol Biol. 1998;7:179–84.

    Article  CAS  PubMed  Google Scholar 

  40. Martinez-Torres D, Foster SP, Field LM, Devonshire AL, Williamson MS. A sodium channel point mutation is associated with resistance to DDT and pyrethroid insecticides in the peach-potato aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae). Insect Mol Biol. 1999;8:339–46.

    Article  CAS  PubMed  Google Scholar 

  41. Tan WL, Wang ZM, Li CX, Chu HL, Xu Y, Dong YD, et al. First report on co-occurrence knockdown resistance mutations and susceptibility to beta-cypermethrin in Anopheles sinensis from Jiangsu Province, China. PLoS ONE. 2012;7:e29242.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kang S, Jung J, Lee S, Hwang H, Kim W. The polymorphism and the geographical distribution of the knockdown resistance (kdr) of Anopheles sinensis in the Republic of Korea. Malar J. 2012;11:151.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Sun YW. Study on deltamethrin resistance mechanism of Anopheles sinensis in central and southern of China. Beijing: Chinese Center for Disease Control and Prevention; 2016. p. 1–78.

  44. Tan WL, Li CX, Wang ZM, Liu MD, Dong YD, Feng XY, et al. First detection of multiple knockdown resistance (kdr)-like mutations in voltage-gated sodium channel using three new genotyping methods in Anopheles sinensis from Guangxi Province, China. J Med Entomol. 2012;49:1012–20.

    Article  CAS  PubMed  Google Scholar 

  45. Yang FL, Li XD, Yan ZT, Fu WB, Chen B. Genetic variation and population structure of Anopheles sinensis ( Diptera: Culicidae) in Yunnan. Acta Entomol Sin. 2015;35:5449–57.

    Google Scholar 

Download references

Acknowledgements

We wish to thank all the participants for volunteering their time to participate in this study.

Funding

This work was supported by the Henan Medical Science and Technology Research Program (LHGJ20210145, LHGJ20220178).

Author information

Author notes

  1. Zhi-quan He and Ya-bo Hu contributed equally to this work.

Authors and Affiliations

  1. Department of Parasite Disease Control and Prevention, Henan Center for Disease Control and Prevention, Zhengzhou, China

    Zhi-quan He, Dan Wang, Cheng-yun Yang, Dan Qian, Rui-min Zhou, De-ling Lu, Su-hua Li, Ying Liu & Hong-wei Zhang

  2. Henan Provincial Medical Key Laboratory of Parasitic Pathogen and Vector, No. 105 South Agricultural Road, Zhengdong New District, Zhengzhou, 450016, China

    Zhi-quan He, Dan Wang, Cheng-yun Yang, Dan Qian, Rui-min Zhou, De-ling Lu, Su-hua Li, Ying Liu & Hong-wei Zhang

  3. Zhengzhou Second Hospital, No. 90 Hanghai Middle Road, Erqi District, Zhengzhou, 450016, China

    Ya-bo Hu

  4. College of Public Health, Zhengzhou University, No.100 Science Avenue, High-Tech District, Zhengzhou, 450016, China

    Yu-ting Liu

Authors
  1. Zhi-quan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ya-bo Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Dan Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yu-ting Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Cheng-yun Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Dan Qian
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Rui-min Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. De-ling Lu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Su-hua Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Ying Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Hong-wei Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ZQH, YBH, YL, and HWZ conceived the study. DW and YTL analyzed the data. YL, SHL, DQ and CYY participated in the study. ZQH drafted the manuscript. DLL, RMZ and HWZ revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ying Liu or Hong-wei Zhang.

Ethics declarations

Ethics approval and consent to participate

No specific permits were required for the described field studies. For mosquito collection, oral consent was obtained from field and house owners in each location. These locations were not protected land, and the field studies did not involve endangered or protected species.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Distribution and frequency of 10 genotypes at 1014 in the kdr gene in Anopheles sinensis.

Additional file 2: Table S2.

Genetic distance and genetic differentiation index of the two populations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, Zq., Hu, Yb., Wang, D. et al. Insecticide resistance of Anopheles sinensis after elimination of malaria in Henan Province, China. Parasites Vectors 16, 180 (2023). https://doi.org/10.1186/s13071-023-05796-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-023-05796-z

Keywords

Parasites & Vectors

ISSN: 1756-3305

Contact us