Knockdown resistance of Anopheles sinensis in Henan province, China

Background Vivax malaria was historically epidemic in Henan Province of China and Anopheles sinensis was the main vectors and poor farming communities bare the greatest burden of disease. Knockdown resistance in An. sinensis is one of the mechanisms of resistance against pyrethroids. In the present study, the frequency of mutations from An. sinensis was examined in Henan province, China. Methods Anopheles was collected from Kaifeng, Tongbai, Tanghe, Pingqiao, Shihe, and Yongcheng counties of Henan province in 2013. Molecular identification of Anopheles species was conducted by polymerase chain reaction (PCR) amplifying the internal transcribed spacer 2 (ITS2). Part of the IIS6 domain of the para-type sodium channel protein gene was polymerase chain reaction-amplified and directly sequenced. Frequency and geographic difference of kdr gene mutant types were analysed. Results 208 Anopheles were received molecular identification, of which 169 (81.25%) were An. sinensis, 25 (12.02%) were Anopheles yatsushiroensis, and 12 (5.77%) were Anopheles lesteri. A 325 bp fragment of the para-type sodium channel gene including position 1014 was successfully sequenced from 139 Anopheles, of which 125 (89.93%) were An. sinensis, 12 (8.63%) were An. yatsushiroensis, 2 (1.44%) were An. lesteri. The molecular analyses revealed that mutations existed at codon 1014 in An. sinensis but not in An. yatsushiroensis and An. lesteri. Frequency of kdr mutation was 73.60% (92/125) from population of An. sinensis in Henan province, of which L1014F (TTT + TTC) allele frequencies accounted for 46.40% (58/125), and was higher than that of L1014C(TGT) which accounted for 27.20% (34/125) ( χ2 = 55.423, P < 0.001). The frequency of kdr mutation in Kaifeng county was 100% (49/49), and was higher than that of 37.93% (11/29) in Tongbai, 54.55% (6/11) in Pingqiao, 50.00% (3/3) in Shihe, and 62.50% (10/16) in Yongcheng county, respectively (χ2 = 39.538, P < 0.001; χ2 = 24.298, P < 0.001; χ2 = 25.913, P < 0.001; χ2 = 20.244, P < 0.001). While 92.86% (13/14) frequency of kdr mutation was found in Tanghe county, which was higher than that in Tongbai county (χ2 = 11.550, P = 0.0018). Conclusions A high frequency of kdr gene mutations from population of An. sinensis in Henan province was found.


Background
Anopheles sinensis is the main vivax malaria vectors and poor farming communities bear the greatest burden of disease in China and other Southeast Asian countries [1][2][3][4]. Indoor residual spraying and long-lasting insecticidal nets are recommend by the World Health Organization (WHO) as effective vector control measures to prevent malaria transmission [5,6]. Nowadays, pyrethroids are emerging as the predominant insecticides for vector control because of their low toxicity to humans, high efficacy against mosquito vectors and short residual action. In the past decade, pyrethroids have become the preferred choice among the currently WHO approved compounds. 414 tonnes of pyrethroids were used annually for global vector control during the period 2000-2009 in the world. 68% (282/414 tonnes) of pyrethroid for residual spraying, 24% (100/414 tonnes) for space spraying, and the remainder for treatment of nets and larviciding [7]. The exploitation of pyrethroids in China started from 1970s, and has been used throughout the country in order to control medically and agriculturally important arthropod pests, including mosquitoes. The area treated with pyrethroids occupies more than one third of the total insecticide-treated area in China [8]. It is critical that the susceptibility of malaria vectors to pyrethroids is preserved. Indeed, it has been recommended not to use pyrethroids for indoor residual spraying where there is high coverage with treated nets [9]. Pyrethroid resistance in malaria vectors has been mostly studied in the major African malaria vector, Anopheles gambiae [10][11][12][13][14][15]. High levels of resistance to pyrethroids have been reported in An. sinensis populations from China, Korea, and Mekong region (Vietnam, Cambodia and Laos) [8,[16][17][18][19][20][21]. Resistance to insecticides can arise due to mutations in the insecticide target site (target site resistance), which in the case of pyrethroids is the para-type sodium channel gene, which is known as knockdown resistance (kdr), is caused by a single mutation in the S6 transmembrane segment of domain II in the voltage-gated sodium channel (VGSC) gene [22]. In recent years, An. sinensis in Henan Province has developed high degree of resistance to deltamethrin [23]. In the present study, the frequency of kdr mutations from An. sinensis was detected in Henan province, China.

Mosquito collection
Adult An. sinensis were captured from six different geographical sites in August 2013 in Henan province of China, including Kaifeng, Tongbai, Tanghe, Pingqiao, Shihe, and Yongcheng counties. The mosquitoes (carcasses/intact mosquitoes) were preserved individually in 1.5 ml microtubes for further molecular analyses.

DNA extraction
Each mosquito was used for DNA extraction with the MaqExtractorTm Kit (Toyobo co. Ltd). Briefly, 1) each mosquito was placed at the bottom of a 1.5 ml microtubes. 2) 750 μl lysis and binding solution and 40 μl magnetic beads were added and mixed for 10 min by vortex. 3) Supernatant was removed by magnetic capture. 4) Magnetic beads were washed three times by 900 μl washing solution and 900 μl 70% ethanol respectively. 5) 100 μl sterilized water was added and well mixed for 10 minutes. 6) Supernatant was collected by magnetic capture and place in a fresh tube. 7) Extracted DNA was stored at −20°C for PCR.

Molecular identification and detection of kdr mutation
Molecular identifications of An. sinensis species were conducted by using species-specific primers and amplification of the ITS2 [24]. To determine point mutations of the para-type sodium gene at positions 1014, a 325 bp fragment of the para-type sodium gene including position 1014 was amplified. PCR primers were designed based on the An. sinensis sequences of the DIIS6 region of the para-type sodium gene according to reference [25]. The allele-specific primers designed were: kdrF (5'-TGC CAC TCC GTG TGT TTA GA-3') and kdrR (5'-GAG CGA TGA TGA TCC GAA AT -3') in a reaction mixture (50 μl) that contained 1× Buffer, 1.5 mM of MgCl 2 , 200 μM of each dNTP, 0.5 μM of each primers and 0.625 unit of Taq DNA polymerase. The conditions of PCR were: an initial denaturation at 95°C for 5 min, followed by 35 cycles at 95°C for 30 S, 52°C for 30 S and 72°C for 30 S, and a final extension step at 72°C for 7 min. The PCR products were purified with QIAquick PCR purification kit (Qiagen) and direct sequencing was done at Sangon Biotech Inc. Primers used for sequencing were in both forward and reverse directions using the same primers.

Statistical analysis
The kdr allele frequency was calculated in each population and statistical differences among populations were examined using the χ-test.

Ethics statement
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.

Discussion
According to surveys conducted in 1930s, vivax malaria (temperate strains) was prevalent in Henan province of China [26]. There were two large epidemics of vivax malaria happened respectively in 1960s and at the beginning of 1970s. The incidence in the whole province was as high as 16,944.40/100 thousand in 1970 [27]. With active implementation of malaria control measures (integrated vector control measures and appropriate treatment of malaria cases) for more than 30 years, considerable success had been achieved and human cases infected with Plasmodium vivax have been reduced significantly in Henan Province in the end of 1980s, and malaria incidence was below 1/10,000 in most areas. By 1992, malaria had been nearly eliminated (incidence less than 1/100,000), with only 318 malaria cases observed in the province [28]. From 2000 to 2006, there was a substantial increase in malaria cases due to re-emerging vivax malaria in Huang-Huai plain. Dramatic vivax malaria resurgence appeared in Yongcheng and Xiayi county of east of Henan province in 2006, which were 307.04% and 360.94% higher than that in 2005 respectively. In Yongcheng county, malaria incidence was 1.23/10,000 in 2004 and 5.19/10,000 in 2005, which is 13 and 3.22 times higher than that of the previous year, respectively [29]. 36 malaria outbreaks and 1825 cases were found in four townships of Yongcheng county, accounted for 63.2% of total malaria cases of Henan province, the highest malaria incidence was up to 4.0% in a village with 43 malaria cases [30].
The degree of epidemicity of malaria is decided by many factors, of which vector efficiency is one of the most important ones. Ten Anopheles species were found in Henan province including An. sinensis, Anopheles lindesayi, Anopheles koreicus, Anopheles kweiyangensis, Anopheles minimus, Anopheles pattoni, Anopheles maculates, Anopheles gigas baileyi, Anopheles anthropophagus, and An. yatsushiroensis. The only two vector species of vivax malaria in Henan province are An. sinensis and An. anthropophagus [27] -An. anthropophagus and An. lesteri were same species. Anopheles yatsushiroensis has recently been declared a synonym of Anopheles pullus [31,32] and a possible vector of vivax malaria in Korea [33,34]. In this study, An. yatsushiroensis is the second most common Anopheles species in Pingqiao county (44.44%), after An. sinensis. The next step is to understand whether or not An. yatsushiroensis be a possible vector of vivax malaria and the pyrethroid susceptibility to this specie in Pingqiao.
Vector capacity of An. anthropophagus is higher than that of An. sinensis, but An. sinensis is the major transmitting vectors in the central provinces due to its widespread distribution [35,36]. The vectorial efficiency of An. sinensis increased during the warmest months from June to August of the year when people frequently sleep outdoor near their fields and unprotected by bed nets [37]. The re-emergence in the Huanghuai plain of central China, including the four provinces of Anhui, Henan, Hubei and Jiangsu were associated with the predominant vector An. sinensis [38], which also plays an important role in the maintenance of P. vivax malaria transmission [39]. Pan et al. investigated vectorial capacities of An. sinensis in 2007, the vectorial capacities of An. sinensis in Huaiyuan and Yongcheng county were 0.7740 and 0.5502, respectively [40]. The results showed that the vector capacity was about 2.3 and 1.7 times higher than 0.331 in the 1990s [41], and was 4.6 and 3.3 times higher than 0.1686 in Henan during 1996-1998 [42], respectively. It Table 2 Frequency (% (n)) of kdr alleles in three An. sinensis population from Henan province Note: a. the frequency of kdr mutation in Kaifeng county was higher than that in Tongbai, Pingqiao, Shihe, and Yongcheng county (P < 0.001); b. the frequency of kdr mutation in Tanghe county was higher than that in Tongbai county (P = 0.001).
was considered that An. sinensis was the sole potential vector of P. vivax malaria in Yongcheng city of Henan province with a 2.78-fold vectorial capacity in 2010 (0.4689) compared to 0.1686 in the 1990s [43,44]. One of the most effective measures to prevent malaria transmission relies on vector control through the use of insecticides, primarily pyrethroids. The overdoses of pyrethroids poses strong selection pressure on mosquito populations for resistance, and have quickly led to the presence and spread of insecticide-resistant mosquitoes, which have caused serious problems for malaria-controlling interventions [34]. The results showed that the KT50 of An. sinensis to deltamethrin was 1122.50, 89.65, and 960 min in Tongbai, Huaibin and Yongcheng county, respectively. The mortality rate of An. sinensis in 24 h-post-exposure to 0.05% deltamethrin was 92.08%, 77.14%, and 63.46%, with a resistance degree of M, R and R, respectively. The results showed that An. sinensis has developed high degree of resistance to deltamethrin in Henan province [23]. The same results were found in Anhui, Hubei, Jiangsu and Shandong province of central China [19,45,46], and Zhejiang, Hainan, Guangxi and Hunan of south China [47][48][49][50], which suggested that pyrethroid resistance was already widespread in natural populations of China.
Knockdown resistance, which are mutations in the paratype sodium channel gene, the target site of pyrethroids is one major resistance mechanisms, which causing a change in affinity between the insecticide and its binding site that reduces sensitivity to the insecticide. In present study, kdr mutation An. sinensis from Henan province was examined. A high frequency (100.00%) of kdr mutations was found in populations from Kaifeng county, but not in populations from Tongbai county (37.93%), with the average frequency of 73.60% (92/125) in Henan province. The previous studies reported that the frequencies of the kdr allele of An. sinensis in China ranged from <10% to >85%, indicating a similar genetic outcome under selective pressure from insecticide treatment [20,21,50,51].
In this study, molecular analysis of kdr gene revealed that mutations at codon 1014 existed only in An. sinensis, whereas no kdr mutations were observed in An. yatsushiroensis and An. lesteri. Several mutations at codon 1014 of the kdr allele, such as L1014F (Leu-to-Phe), L1014S (Leuto-Ser), and L1014C (Leuto-Cys) have been reported in many Anopheles species [52][53][54][55]. In this study, frequency of L1014F allele accounted for 46.40% (58/125), and was higher than that of L1014C which accounted for 27.20% (34/125) (χ2 = 55.423, P < 0.001). The results suggested that L1014F mutation was a major allele that showed a high allele frequency, whereas L1014C mutation was a minor allele that showed a low allele frequency within the An. sinensis populations in Henan province.
A study using stepwise multiple regression analyses in mosquito populations from central China found that both kdr mutations and monooxygenase activity were significantly associated with deltamethrin resistance, with monooxygenase activity playing a stronger role [24]. The results suggest that different mechanisms of resistance could evolve in geographically different populations.

Conclusions
The observed pyrethroid resistance and high kdr mutation frequency in populations of An. sinensis could profoundly affect the current malaria vector control programme in Henan province. The identification of widespread kdr mutations suggests that the development of reliable resistance surveillance tools is an important topic for future research. This needs an urgent call for implementing rational resistance management strategies and integrated vector control intervention.