Detection and widespread distribution of sodium channel alleles characteristic of insecticide resistance in Culex pipiens complex mosquitoes in China

Mutations and intron polymorphisms in voltage-gated sodium channel genes of different geographic populations of Culex pipiens pallens/Culex pipiens quinquefasciatus in China

BACKGROUND

Culex pipiens pallens and Culex pipiens quinquefasciatus are the dominant species of Culex mosquitoes in China and important disease vectors. Long-term use of insecticides can cause mutations in the voltage-gated sodium channel (vgsc) gene of mosquitoes, but little is known about the current status and evolutionary origins of vgsc gene in different geographic populations. Therefore, this study aimed to determine the current status of vgsc genes in Cx. p. pallens and Cx. p. quinquefasciatus in China and to investigate the evolutionary inheritance of neighboring downstream introns of the vgsc gene to determine the impact of insecticides on long-term evolution.

METHODS

Sampling was conducted from July to September 2021 in representative habitats of 22 provincial-level administrative divisions in China. Genomic DNA was extracted from 1308 mosquitoes, the IIS6 fragment of the vgsc gene on the nerve cell membrane was amplified using polymerase chain reaction, and the sequence was used to evaluate allele frequency and knockdown resistance (kdr) frequency. MEGA 11 was used to construct neighbor-joining (NJ) tree. PopART was used to build a TCS network.

RESULTS

There were 6 alleles and 6 genotypes at the L1014 locus, which included the wild-type alleles TTA/L and CTA/L and the mutant alleles TTT/F, TTC/F, TCT/S and TCA/S. The geographic populations with a kdr frequency less than 20.00% were mainly concentrated in the regions north of 38° N, and the geographic populations with a kdr frequency greater than 80.00% were concentrated in the regions south of 30° N. kdr frequency increased with decreasing latitude. And within the same latitude, the frequency of kdr in large cities is relatively high. Mutations were correlated with the number of introns. The mutant allele TCA/S has only one intron, the mutant allele TTT/F has three introns, and the wild-type allele TTA/L has 17 introns.

CONCLUSIONS

Cx. p. pallens and Cx. p. quinquefasciatus have developed resistance to insecticides in most regions of China. The neighboring downstream introns of the vgsc gene gradually decreased to one intron with the mutation of the vgsc gene. Mutations may originate from multiple mutation events rather than from a single origin, and populations lacking mutations may be genetically isolated.

Novel real-time PCR assay detects widespread distribution of knock down resistance (kdr) mutations associated with pyrethroid resistance in the mosquito, Culex quinquefasciatus, in Thailand

Susceptibility to pyrethroids in the mosquito Culex quinquefasciatus, the major vector of lymphatic filariasis, is being seriously threatened worldwide. Knockdown resistance (kdr), caused by mutations in the voltage gated sodium channel (VGSC) gene, particularly the L1014F mutation, is an important resistance mechanism. Our aim was to develop a real-time PCR with melt curve analysis to evaluate the distribution of the L1014F mutation in Cx. quinquefasciatus throughout Thailand and to determine the polymorphism pattern of a VGSC gene fragment spanning the L1014F mutation. A total of 3760 females from 18 localities across five regions of Thailand were bio-assayed by exposure to 0.05% deltamethrin WHO papers, showing mortality rates ranging from 2.4% to 83.0%. Genotyping of 753 dead and surviving mosquitoes using our novel real-time PCR assay with melt curve analysis and tetra-primer allele-specific PCR revealed the mutant F1014 allele is closely associated with the deltamethrin resistance phenotype. The L1014F mutation was found at high frequency throughout Thailand, particularly in the North. However, some survivors were homozygous for wild type L1014 allele, which were further sequenced for the IIP-IIS6 region of VGSC gene. The haplotype network of phenotypically characterized individuals indicated the presence of other possible kdr alleles/resistance mechanisms at play including two novel mutations, V978E and D992E. The finding of new putative kdr alleles and widespread distribution of the F1014 allele emphasizes the significant role of kdr mutations in pyrethroid resistance in Thai Cx. quinquefasciatus populations. Monitoring kdr variations and phenotypic resistance is critical for managing resistance in Cx. quinquefasciatus.

Mosquito population structure, pathogen surveillance and insecticide resistance monitoring in urban regions of Crete, Greece

Background

In Greece vector borne diseases (VBD) and foremost West Nile virus (WNV) pose an important threat to public health and the tourist industry, the primary sector of contribution to the national economy. The island of Crete, is one of Greece’s major tourist destinations receiving annually over 5 million tourists making regional VBD control both a public health and economic priority.

Methodology

Under the auspices of the Region of Crete, a systematic integrative surveillance network targeting mosquitoes and associated pathogens was established in Crete for the years 2018–2020. Using conventional and molecular diagnostic tools we investigated the mosquito species composition and population dynamics, pathogen infection occurrences in vector populations and in sentinel chickens, and the insecticide resistance status of the major vector species.

Principal findings

Important disease vectors were recorded across the island including Culex pipiens, Aedes albopictus, and Anopheles superpictus. Over 75% of the sampled specimens were collected in the western prefectures potentially attributed to the local precipitation patterns, with Cx. pipiens being the most dominant species. Although no pathogens (flaviviruses) were detected in the analysed mosquito specimens, chicken blood serum analyses recorded a 1.7% WNV antibody detection rate in the 2018 samples. Notably detection of the first WNV positive chicken preceded human WNV occurrence in the same region by approximately two weeks. The chitin synthase mutation I1043F (associated with high diflubenzuron resistance) was recorded at an 8% allelic frequency in Lasithi prefecture Cx. pipiens mosquitoes (sampled in 2020) for the first time in Greece. Markedly, Cx. pipiens populations in all four prefectures were found harboring the kdr mutations L1014F/C/S (associated with pyrethroid resistance) at a close to fixation rate, with mutation L1014C being the most commonly found allele (≥74% representation). Voltage gated sodium channel analyses in Ae. albopictus revealed the presence of the kdr mutations F1534C and I1532T (associated with putative mild pyrethroid resistance phenotypes) yet absence of V1016G. Allele F1534C was recorded in all prefectures (at an allelic frequency range of 25–46.6%) while I1532T was detected in populations from Chania, Rethymnon and Heraklion (at frequencies below 7.1%). Finally, no kdr mutations were detected in the Anopheles specimens included in the analyses.

Conclusions/Significance

The findings of our study are of major concern for VBD control in Crete, highlighting (i) the necessity for establishing seasonal integrated entomological/pathogen surveillance programs, supporting the design of targeted vector control responses and; ii) the need for establishing appropriate insecticide resistance management programs ensuring the efficacy and sustainable use of DFB and pyrethroid based products in vector control.

The island of Crete, is one of Greece’s major tourist destinations, receiving over 5 million tourists annually, making control of vector borne diseases (VBD) like West Nile virus (WNV) both a public health and economic priority. A comprehensive surveillance network targeting mosquitoes and associated pathogens was established in Crete for the years 2018–2020. Important disease vectors (Culex pipiens, Aedes albopictus and Anopheles superpictus) were recorded across the island. The presence of WNV antibodies in sentinel chickens was also recorded and, notably, preceded human WNV occurrence in the same region by approximately two weeks. Mutations associated with resistance to insecticides used for vector control were also detected; most importantly the chitin synthase mutation I1043F (associated with high diflubenzuron resistance) was recorded for the first time in Greece. The findings of our study are of major concern for VBD control in Greece (and Crete specifically), highlighting (i) the necessity for establishing seasonal surveillance programs and ii) the need for establishing appropriate insecticide resistance management programs for sustainable vector control.

Prevalence and molecular characterization of Wolbachia in field-collected Aedes albopictus, Anopheles sinensis, Armigeres subalbatus, Culex pipiens and Cx. tritaeniorhynchus in China

Wolbachia are maternally transmitted intracellular bacteria that can naturally and artificially infect arthropods and nematodes. Recently, they were applied to control the spread of mosquito-borne pathogens by causing cytoplasmic incompatibility (CI) between germ cells of females and males. The ability of Wolbachia to induce CI is based on the prevalence and polymorphism of Wolbachia in natural populations of mosquitoes. In this study, we screened the natural infection level and diversity of Wolbachia in field-collected mosquitoes from 25 provinces of China based on partial sequence of Wolbachia surface protein (wsp) gene and multilocus sequence typing (MLST). Among the samples, 2489 mosquitoes were captured from 24 provinces between July and September, 2014 and the remaining 1025 mosquitoes were collected month-by-month in Yangzhou, Jiangsu province between September 2013 and August 2014. Our results showed that the presence of Wolbachia was observed in mosquitoes of Aedes albopictus (97.1%, 331/341), Armigeres subalbatus (95.8%, 481/502), Culex pipiens (87.0%, 1525/1752), Cx. tritaeniorhynchus (17.1%, 14/82), but not Anopheles sinensis (n = 88). Phylogenetic analysis indicated that high polymorphism of wsp and MLST loci was observed in Ae. albopictus mosquitoes, while no or low polymorphisms were in Ar. subalbatus and Cx. pipiens mosquitoes. A total of 12 unique mutations of deduced amino acid were identified in the wsp sequences obtained in this study, including four mutations in Wolbachia supergroup A and eight mutations in supergroup B. This study revealed the prevalence and polymorphism of Wolbachia in mosquitoes in large-scale regions of China and will provide some useful information when performing Wolbachia-based mosquito biocontrol strategies in China.

The mosquitoes Aedes albopictus, Anopheles sinensis, Armigeres subalbatus, Culex pipiens and Cx. tritaeniorhynchus are native to China and the major vectors in the transmission of arboviruses, protozoans and nematodes. Recently, an innovative biocontrol strategy has been developed and evaluated based on the ability of Wolbachia to induce cytoplasmic incompatibility (CI), as well as interfere with the infection and replication of pathogens. Since the ability to induce CI largely depends on the density and diversity of Wolbachia, we investigated and characterized the natural infection of Wolbachia in above-mentioned five species of field-collected mosquitoes in 25 provinces of China. The results showed that the positive rates of Wolbachia infection were high in mosquitoes of Ae. albopictus, Ar. subalbatus and Cx. pipiens in large-scale regions of China and low in Cx. tritaeniorhynchus in Guizhou province. Phylogenetic analysis based on Wolbachia surface protein (wsp) gene and five multilocus sequence typing (MLST) loci indicated the high polymorphism of Wolbachia in Ae. albopictus, and low polymorphisms in Ar. subalbatus and Cx. pipiens. This finding contributes to the understanding of the nationwide distribution of Wolbachia and the potential application of this biocontrol strategy in China.

Microbiota and transcriptome changes of Culex pipiens pallens larvae exposed to Bacillus thuringiensis israelensis

Culex pipiens pallens is an important vector of lymphatic filariasis and epidemic encephalitis. Mosquito control is the main strategy used for the prevention of mosquito-borne diseases. Bacillus thuringiensis israelensis (Bti) is an entomopathogenic bacterium widely used in mosquito control. In this study, we profiled the microbiota and transcriptional response of the larvae of Cx. pipiens pallens exposed to different concentrations of Bti. The results demonstrated that Bti induced a significant effect on both the microbiota and gene expression of Cx. pipiens pallens. Compared to the control group, the predominant bacteria changed from Actinobacteria to Firmicutes, and with increase in the concentration of Bti, the abundance of Actinobacteria was gradually reduced. Similar changes were also detected at the genus level, where Bacillus replaced Microbacterium, becoming the predominant genus in Bti-exposed groups. Furthermore, alpha diversity analysis indicated that Bti exposure changed the diversity of the microbota, possibly because the dysbiosis caused by the Bti infection inhibits some bacteria and provides opportunities to other opportunistic taxa. Pathway analysis revealed significant enhancement for processes associated with sphingolipid metabolism, glutathione metabolism and glycerophospholipid metabolism between all Bti-exposed groups and control group. Additionally, genes associated with the Toll and Imd signaling pathway were found to be notably upregulated. Bti infection significantly changed the bacterial community of larvae of Cx. pipiens pallens.

Monitoring Insecticide Resistance and Target Site Mutations of L1014 Kdr And G119 Ace Alleles in Five Mosquito Populations in Korea

Mosquitoes are globally distributed and important vectors for the transmission of many human diseases. Mosquito control is a difficult task and the cost of preventing mosquito-borne diseases is much lower than that for curing the associated diseases. Thus, chemical control remains the most effective tool for mosquito. Due to the long-term intensive use of insecticides to control mosquito vectors, resistance to most chemical insecticides has been reported. This study aimed to investigate the relationship between insecticide resistance and target site mutation of L1014 kdr and G119 ace alleles in 5 species/species group of mosquitoes (Aedes vexans, Ae. albopictus, Anopheles spp., Culex pipiens complex, and Cx. tritaeniorhynchus) obtained from 6 collection sites. For Anopheles spp., the proportion of mosquitoes with mutated alleles in L1014 was 88.4%, homozygous resistant genotypes were observed in 46.7%, and heterozygous resistant genotypes were observed in 41.8%. For the Cx. pipiens complex and Cx. tritaeniorhynchus species, homozygous resistant genotypes were found in 25.9% and 9.8%, respectively. However, target site mutation of L1014 in the Ae. vexans nipponii and Ae. albopictus species was not observed. Anopheles spp., Cx. pipiens complex, and Cx. tritaeniorhynchus mosquitoes were resistant to deltamethrin and chlorpyriphos, whereas Ae. vexans nipponii and Ae. albopictus were clearly susceptible. We also found a correlation between the resistance phenotype and the presence of the L1014 kdr and G119 ace mutations only in the Anopheles spp. population. In this study, we suggest that insecticide resistance poses a growing threat and resistance management must be integrated into all mosquito control programs.

Susceptibility to Pyrethroids and the First Report of L1014F kdr Mutation in Culex quinquefasciatus (Diptera: Culicidae) in Colombia

The use of insecticides for the control of Aedes aegypti (L.) (Diptera: Culicidae) in Colombia has indirectly influenced the susceptibility status of Culex quinquefasciatus Say populations. We evaluated pyrethroid susceptibility in two populations of Cx. quinquefasciatus in the Atlantico Department of Colombia and its possible resistance mechanism (kdr mutation). Bottle bioassays were performed for permethrin, deltamethrin, and λ-cyhalothrin in female mosquitoes of Cx. quinquefasciatus. The resistance ratios (RRs) for KC50 and LC50 for each insecticide in the field populations examined were determined, using the Cartagena strain as the susceptible control. The L1014F kdr mutation was identified in the para gene of the voltage-gated sodium channel (vgsc), along with its allelic and genotypic frequency. Low knockdown resistance (RRKC50) to deltamethrin was found in Puerto Colombia and Soledad populations as well as low resistance to λ-cyalothrin in this latter population. Moderate knockdown resistance to permethrin was found in both populations. At 24 h post-exposure on the other hand, there was low resistance (RRLC50) to permethrin in Puerto Colombia and moderate resistance in Soledad. Moderate resistance to deltamethrin was found in Puerto Colombia and low resistance in Soledad. Low resistance to λ-cyhalothrin was seen in Puerto Colombia and moderate resistance in Soledad. Variability was found in the susceptibility to the pyrethroids in the populations of Cx. quinquefasciatus evaluated, and the L1014F kdr mutation is reported for the first time as a possible pyrethroid resistance mechanism in this species in Colombia.

Breeding Site Characteristics and Associated Factors of Culex pipiens Complex in Lhasa, Tibet, P. R. China

Characterizing the breeding sites of Culex pipiens complex is of major importance for the control of West Nile disease and other related diseases. However, little information is available about the characteristics and associated factors of the breeding sites of the Cx. pipiens complex in Lhasa, a representative high-altitude region in Southwestern China. In this study, a cross-sectional study concerning the breeding site characteristics and associated factors of the Cx. pipiens complex was carried out in Lhasa, Tibet from 2013–2016. Chi-square analysis and binary logistic regression analysis were applied to identify the key factors associated with the presence of Cx. pipiens complex larvae. Using a standard dipping method, 184 water bodies were examined and Cx. pipiens complex larvae were observed in 36 (19.57%) of them. There were significant differences in the composition of Cx. pipiens complex larvae among the breeding site stability (χ² = 19.08, p = 0.00) and presence or absence of predators (χ² = 6.986, p = 0.008). Binary logistic regression analysis indicated that breeding site stability and presence or absence of predators were significantly associated with the presence of Cx. pipiens complex larvae in Chengguan District, Lhasa. Relatively permanent water bodies such as water bodies along river fringes, ponds and puddles, and water bodies with no predators should be paid more attention for future Cx. pipiens complex larvae abatement campaigns in Lhasa, China.

Identification and classification of differentially expressed genes in pyrethroid-resistant Culex pipiens pallens

Culex pipiens pallens is an important vector that transmits Bancroftian filariasis, Japanese encephalitis and other diseases that pose a serious threat to human health. Extensive and improper use of insecticides has caused insecticide resistance in mosquitoes, which has become an important obstacle to the control of mosquito-borne diseases. It is crucial to investigate the underlying mechanism of insecticide resistance. The aims of this study were to identify genes involved in insecticide resistance based on the resistance phenotype, gene expression profile and single-nucleotide polymorphisms (SNPs) and to screen for major genes controlling insecticide resistance. Using a combination of SNP and transcriptome data, gene expression quantitative trait loci (eQTLs) were studied in deltamethrin-resistant mosquitoes. The most differentially expressed pathway in the resistant group was identified, and a regulatory network was built using these SNPs and the differentially expressed genes (DEGs) in this pathway. The major candidate genes involved in the control of insecticide resistance were analyzed by qPCR, siRNA microinjection and CDC bottle bioassays. A total of 85 DEGs that encoded putative detoxification enzymes (including 61 P450s) were identified in this pathway. The resistance regulatory network was built using SNPs, and these metabolic genes, and a major gene, CYP9AL1, were identified. The functional role of CYP9AL1 in insecticide resistance was confirmed by siRNA microinjection and CDC bottle bioassays. Using the eQTL approach, we identified important genes in pyrethroid resistance that may aid in understanding the mechanism underlying insecticide resistance and in targeting new measures for resistance monitoring and management.

A novel miRNA, miR-13664, targets CpCYP314A1 to regulate deltamethrin resistance in Culex pipiens pallens

Extensive insecticide use has led to the resistance of mosquitoes to these insecticides, posing a major barrier to mosquito control. Previous Solexa high-throughput sequencing of Culex pipiens pallens in the laboratory has revealed that the abundance of a novel microRNA (miRNA), miR-13664, was higher in a deltamethrin-sensitive (DS) strain than a deltamethrin-resistant (DR) strain. Real-time quantitative PCR revealed that the miR-13664 transcript level was lower in the DR strain than in the DS strain. MiR-13664 oversupply in the DR strain increased the susceptibility of these mosquitoes to deltamethrin, whereas inhibition of miR-13664 made the DS strain more resistant to deltamethrin. Results of bioinformatic analysis, quantitative reverse-transcriptase polymerase chain reaction, luciferase assay and miR mimic/inhibitor microinjection revealed CpCYP314A1 to be a target of miR-13664. In addition, downregulation of CpCYP314A1 expression in the DR strain reduced the resistance of mosquitoes to deltamethrin. Taken together, our results indicate that miR-13664 could regulate deltamethrin resistance by interacting with CpCYP314A1, providing new insights into mosquito resistance mechanisms.

Insecticide resistance and target site mutations (G119S ace-1 and L1014F kdr) of Culex pipiens in Morocco

Background

Control of the mosquito vector Culex pipiens with insecticides is the main way to control arboviruses that the species can transmit such as West Nile virus (WNV) and Rift Valley fever virus (RVFV). However, its efficiency has been hampered by the emergence of insecticide resistance. Little is known about the insecticide-resistance status and underlying resistance mechanisms of field-collected populations of Cx. pipiens in Morocco.

Methods

Mosquito adults from Mohammadia city in Morocco were reared from immature stages. The level of their susceptibility to insecticides was assessed using standard WHO bioassay. The two forms of the Cx. pipiens complex and their hybrids were identified by a multiplex PCR. Identified mosquitoes were then tested for the presence of the G119S ace-1 and L1014F kdr mutations using PCR-RFLP and PCR assays, respectively.

Results

WHO bioassays indicated that Cx. pipiens was resistant to all tested insecticides: lambda-cyhalothrin (49% mortality), permethrin (63% mortality), DDT (16% mortality), malation (52% mortality) and bendiocarb (39% mortality). The frequency of the 119S allele was almost identical in the pipiens form and hybrids (0.11 and 0.15, respectively) whereas it remained low in the molestus form (0.03). No significant correlation was observed between the G119S allele and the resistance phenotype to two tested insecticides (malathion and bendiocarb). The frequency of the L1014F allele was identical in the pipiens form and hybrids (0.44) whereas it was low in the molestus form (0.36) but no significant difference was detected (χ² = 1.46, df = 1, P = 0.225). The presence of the L1014F kdr mutation was significantly associated with resistance to three tested insecticides in pipiens form (P = 0.0019, P = 0.0023 and P = 0.023, respectively, to lambda-cyhalothrin, permethrin and DDT) whereas no significant correlation was observed between the L1014F kdr mutation and resistance phenotype in molestus form and hybrids to the three tested insecticides.

Conclusion

These findings showed that wild populations of Cx. pipiens have developed resistance against the main insecticide families with different modes of action: organochlorines (DDT), organophosphates (malathion), carbamates (bendiocarb), pyrethroids (lambda-cyhalothrin, permethrin). Therefore, urgent action should be taken to manage the resistance in this species to maintain the effectiveness of arbovirus control.

Analysis of population structure and insecticide resistance in mosquitoes of the genus Culex, Anopheles and Aedes from different environments of Greece with a history of mosquito borne disease transmission

Greece has been recently affected by several mosquito borne diseases with the West Nile Virus (WNV) outbreak in 2010 being one of the largest reported in Europe. Currently at the epicenter of an economic and refugee crisis and visited by over 16 million tourists a year the integrated management of diseases transmitted by mosquitoes is a public health and economic priority. Vector control programs rely mainly on insecticides, however data on insecticide resistance and the mosquito fauna is essential for successful applications. We determined the mosquito species composition and population dynamics in areas of increased vulnerability to vector borne disease transmission, as well as investigated the resistance status of major nuisance and disease vectors to insecticides. High mosquito densities were recorded in Thessaloniki and Evros, with Aedes caspius, a nuisance species, Culex pipiens, a known vector of WNV and Anopheles hyrcanus a potential vector of malaria being among the most prevalent species. Both vector species populations reached their peak in late summer. Aedes albopictus was recorded at high densities in Thessaloniki, but not in Evros. Notably, Cx. pipiens hybrids, which show an opportunistic biting behavior and are suspected to be involved in the transmission of the WNV, were recorded in considerable numbers in Thessaloniki and Attica. Culex pipiens and An. hyrcanus, but not Ae. caspius mosquitoes, showed moderate levels of resistance to deltamethrin. The presence of resistance in areas not exposed to vector control indicates that other factors could be selecting for resistance, i.e. pesticide applications for agriculture. Both L1014F and L101C kdr mutations were detected in Cx. pipiens populations. Anopheles hyrcanus resistance was not associated with mutations at the L1014 site. The Ace-1 mutations conferring insensitivity to organophosphates and carbamates were detected at low frequencies in all Cx. pipiens populations. Increased activity of P450s and esterases was found in Cx. pipiens individuals from Thessaloniki. Our study contributes evidence for sustainable and efficient vector control strategies and the prevention of disease outbreaks.

Knockdown resistance mutations predict DDT resistance and pyrethroid tolerance in the visceral leishmaniasis vector Phlebotomus argentipes

Background

Indoor residual spraying (IRS) with DDT has been the primary strategy for control of the visceral leishmaniasis (VL) vector Phlebotomus argentipes in India but efficacy may be compromised by resistance. Synthetic pyrethroids are now being introduced for IRS, but with a shared target site, the para voltage-gated sodium channel (VGSC), mutations affecting both insecticide classes could provide cross-resistance and represent a threat to sustainable IRS-based disease control.

Methodology/Principal findings

A region of the Vgsc gene was sequenced in P. argentipes from the VL hotspot of Bihar, India. Two knockdown resistance (kdr) mutations were detected at codon 1014 (L1014F and L1014S), each common in mosquitoes, but previously unknown in phlebotomines. Both kdr mutations appear largely recessive, but as homozygotes (especially 1014F/F) or as 1014F/S heterozygotes exert a strong effect on DDT resistance, and significantly predict survivorship to class II pyrethroids in short-duration bioassays. The mutations are present at high frequency in wild P. argentipes populations from Bihar, with 1014F significantly more common in higher VL areas.

Conclusions/Significance

The Vgsc mutations detected appear to be a primary mechanism underlying DDT resistance in P. argentipes and a contributory factor in reduced pyrethroid susceptibility, suggesting a potential impact if P. argentipes are subjected to suboptimal levels of pyrethroid exposure, or additional resistance mechanisms evolve. The assays to detect kdr frequency changes provide a sensitive, high-throughput monitoring tool to detecting spatial and temporal variation in resistance in P. argentipes.

Visceral leishmaniasis is a fatal disease transmitted solely by the sandfly Phlebotomus argentipes in India. For decades, indoor residual spraying targeting sandflies with DDT has been the main control tool in the region. Emergence of DDT resistance has compromised this strategy and pyrethroids are now being implemented as alternative. Here, we describe the first molecular markers in sandflies for DDT and pyrethroid resistance. Two knockdown resistance (kdr) mutations in the para voltage-gated sodium channel target site gene are found at high frequency in natural populations in Bihar, the Indian focus of VL. The kdr variants strongly predict DDT resistance in P. argentipes and survival to lower-level pyrethroid exposure. Careful monitoring of resistance, aided by these markers, and of spraying efficacy is required to monitor changes and forewarn of emergent higher-level pyrethroid resistance in the region.

Sodium Channel Mutations and Pyrethroid Resistance in Aedes aegypti

Pyrethroid insecticides are widely used to control insect pests and human disease vectors. Voltage-gated sodium channels are the primary targets of pyrethroid insecticides. Mutations in the sodium channel have been shown to be responsible for pyrethroid resistance, known as knockdown resistance (kdr), in various insects including mosquitoes. In Aedes aegypti mosquitoes, the principal urban vectors of dengue, zika, and yellow fever viruses, multiple single nucleotide polymorphisms in the sodium channel gene have been found in pyrethroid-resistant populations and some of them have been functionally confirmed to be responsible for kdr in an in vitro expression system, Xenopus oocytes. This mini-review aims to provide an update on the identification and functional characterization of pyrethroid resistance-associated sodium channel mutations from Aedes aegypti. The collection of kdr mutations not only helped us develop molecular markers for resistance monitoring, but also provided valuable information for computational molecular modeling of pyrethroid receptor sites on the sodium channel.

Diversity and frequency of kdr mutations within Anopheles sinensis populations from Guangxi, China

BACKGROUND

Anopheles sinensis is a major vector of malaria in China and its control is under great threat as the development of insecticide resistance. Voltage-gated sodium channel (VGSC) is the target of several classes of insecticides. Genetic mutations of VGSC have been documented to confer knockdown resistance (kdr) to dichlorodiphenyltrichloroethane (DDT) and pyrethroids in mosquitoes. To control this vector efficiently, it is important to know the resistance-associated genetic mutations, their distribution frequencies and genealogical relations.

METHODS

Three hundreds and thirteen (313) adults of An. sinensis collected from nine locations across Guangxi Zhuang Autonomous Region were used. The partial sequence of the An. sinensis voltage gated sodium channel gene (AS-VGSC) containing codon 1014 was sequenced. PHASE2.1 was used to construct the haplotypes of each individual, and the accuracy of haplotypes was further confirmed by clone sequencing. The genealogical relations of kdr mutations in AS-VGSC was analysed using TCS 2.1 and Network 5.0.

RESULTS

Sixteen AS-VGSC haplotypes including seven haplotypes carrying non-synonymous mutations at codon 1014, and fifty-five AS-VGSC genotypes were identified from 313 mosquitoes collected from nine geographical locations across Guangxi. The number of haplotypes in each of the nine populations ranged from 5 to 13. The frequency of haplotypes carrying kdr mutations ranged from 2.7 to 80.0 % within the nine populations, of which 1014C was unexpectedly high in the northeast of Guangxi. Genealogical analysis suggested multiple origins of kdr mutations in An. sinensis.

CONCLUSION

Diverse haplotypes of AS-VGSC are distributed in Guangxi. The presence of haplotypes carrying mutations at codon 1014 indicates a risk of pyrethroid and DDT resistance. The kdr mutations show differential distribution geographically, with high frequencies occurred in the northeast of Guangxi. Genealogical analysis suggests multiple origins of kdr mutations in An. sinensis populations in Guangxi. These findings have important practical implications for the sustainability of An. sinensis control programmes.

Seasonal dynamics of insecticide resistance, multiple resistance, and morphometric variation in field populations of Culex pipiens

Resistance to insecticides that impairs nervous transmission has been widely investigated in mosquito populations as insecticides are crucial to effective insect control. The development of insecticide resistance is also of special interest to evolutionary biologists since it represents the opportunity to observe the genetic consequences of a well-characterized alteration in the environment. Although the frequencies of resistance alleles in Culex pipiens populations against different groups of insecticides have been reported, no detailed information is available on the relative change in these allele frequencies over time. In this study, we collected mosquitoes of the Cx. pipiens complex from six locations in three seasons in the Aegean region of Turkey and examined the i) seasonal variations in resistance to four different chemical classes of insecticides, ii) seasonal fluctuations in frequencies of resistance-associated target-site mutations of the three genes (ace-1, kdr, and Rdl), and iii) potential seasonal variations in wing morphometric characters that may be modified in resistant mosquitoes. Our bioassay results indicated the presence of different levels of resistance to all tested insecticides for all three seasons in all locations. The results of the PCR-based molecular analysis revealed low frequencies of mutations in ace-1 and Rdl that are associated with resistance to malathion, bendiocarb, and dieldrin and no obvious seasonal changes. In contrast, we detected high frequencies and striking seasonal changes for two kdr mutations associated with resistance to DDT and pyrethroids. In addition, the evaluation of the field populations from all seasons in terms of the combinations of polymorphisms at four resistance-associated mutations did not reveal the presence of insects that are resistant to all pesticides. Results from the morphological analysis displayed a similar pattern for both wings and did not show a clear separation among the samples from the three different seasons. The results of this study have advanced our knowledge of the potential dynamics of insecticide resistance among populations of the Cx. pipiens complex. The implications of these results to the understanding of the evolution of insecticide resistance and the management of resistance in mosquitoes are discussed.

Comparative transcriptome analyses of deltamethrin-susceptible and -resistant Culex pipiens pallens by RNA-seq

The widespread and improper use of pyrethroid insecticides, such as deltamethrin, has resulted in the evolution of resistance in many mosquito species, including Culex pipiens pallens. With the development of high-throughput sequencing, it is possible to massively screen pyrethroid resistance-associated gene. In this study, we used Illumina-Solexa transcriptome sequencing to identify genes that are expressed differently in deltamethrin-susceptible and -resistant strains of Culex pipiens pallens as a critical knowledge base for further studies. A total of 4,961,197,620 base pairs and 55,124,418 reads were sequenced, mapped to the Culex quinquefasciatus genome and assembled into 17,679 known genes. We recorded 1826 significantly differentially expressed genes (DEGs). Among them, 1078 genes were up-regulated and 748 genes were down-regulated in the deltamethrin-resistant strain compared to -susceptible strain. These DEGs contained cytochrome P450 s, cuticle proteins, UDP-glucuronosyltransferases, lipases, serine proteases, heat shock proteins, esterases and others. Among the 1826 DEGs, we found that the transcriptional levels of CYP6AA9 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Moreover, the expression levels of the CYP6AA9 were significantly higher in the resistant strains than the susceptible strains in three different field populations. We further confirmed the association between the CYP6AA9 gene and deltamethrin resistance in mosquitoes by RNA interfering (RNAi). Altogether, we explored massive potential pyrethroid resistance-associated genes and demonstrated that CYP6AA9 participated in the pyrethroid resistance in mosquitoes.

Venom allergen 5 is Associated With Deltamethrin Resistance in Culex pipiens pallens (Diptera: Culicidae)

The mosquito, Culex pipiens pallens (L.), is an important vector of encephalitis and filariasis in northern China. The control of these mosquitoes occurs primarily via the use of pyrethroid insecticides, such as deltamethrin. The widespread and improper application of pyrethroid has resulted in the evolution of pyrethroid resistance amongst many mosquito populations, including Cx. pipiens pallens. Previous studies using high-throughput transcriptome sequencing have identified that the venom allergen 5 gene is differentially expressed between deltamethrin-susceptible and deltamethrin-resistant Cx. pipiens pallens. In this study, quantitative real-time polymerase chain reaction analyses revealed that venom allergen 5 was significantly overexpressed in adult females of both deltamethrin-resistant laboratory populations and two field populations. The transcriptional level of venom allergen 5 in the laboratory populations was elevated as the levels of deltamethrin resistance increased. Full-length cDNAs of the venom allergen 5 gene were cloned from Cx. pipiens pallens, and contained an open reading frame of 765 bp, encoding a protein with 254 amino acids. The deduced amino acid sequence shared 100% identity with the ortholog in Culex quinquefasciatus Say. The overexpression of venom allergen 5 decreased the susceptibility of mosquito cells to deltamethrin, while knockdown of this gene by RNAi increased the susceptibility of mosquitoes to deltamethrin. This study provides the first evidence of the association between the venom allergen 5 gene and deltamethrin resistance in mosquitoes.

Pyrethroid resistance in Culex pipiens mosquitoes

Mosquitoes within the Culex pipiens complex are widely distributed and important in the transmission of many human diseases. Insecticides, pyrethroids in particular, remain a mainstay for control of these important vectors. In this paper we review what is known about the levels, mechanisms and fitness costs of pyrethroid resistance in Cx. pipiens. Pyrethroid resistance in Cx. pipiens is a global problem, and resistance ratios of up to 7000-fold have been found in larvae of field collected mosquitoes. However, there is considerable variation between populations, indicating significant geographic heterogeneity of the resistance. The two major mechanisms of resistance to pyrethroids in Culex are mutations in Vssc (target site insensitivity) and overexpression of cytochrome P450(s) (increased detoxification). The most frequently reported Vssc mutation is L1014F (i.e. kdr), which has been found throughout the world. The L1014S mutation has been found in Cx. p. pallens from Japan and China, and in Cx. p. pipiens from China. The L1014C mutation has only been reported for Cx. p. pipens molestus from China and the V1016G mutation has only been reported from Saudi Arabia. Studies on the P450s of Cx. pipiens have identified several that are overexpressed (measured as transcript levels) in pyrethroid resistant strains. CYP9M10 is consistently overexpressed in pyrethroid resistant Cx. pipiens from at least seven countries, suggesting this P450 might be of global importance in resistance. Both CYP9M10-mediated pyrethroid resistance and kdr have fitness costs in the absence of insecticides under certain environmental conditions. Research needs and future directions are discussed.

A mutation in the intracellular loop III/IV of mosquito sodium channel synergizes the effect of mutations in helix IIS6 on pyrethroid resistance

Activation and inactivation of voltage-gated sodium channels are critical for proper electrical signaling in excitable cells. Pyrethroid insecticides promote activation and inhibit inactivation of sodium channels, resulting in prolonged opening of sodium channels. They preferably bind to the open state of the sodium channel by interacting with two distinct receptor sites, pyrethroid receptor sites PyR1 and PyR2, formed by the interfaces of domains II/III and I/II, respectively. Specific mutations in PyR1 or PyR2 confer pyrethroid resistance in various arthropod pests and disease vectors. Recently, a unique mutation, N(1575)Y, in the cytoplasmic loop linking domains III and IV (LIII/IV) was found to coexist with a PyR2 mutation, L(1014)F in IIS6, in pyrethroid-resistant populations of Anopheles gambiae. To examine the role of this mutation in pyrethroid resistance, N(1575)Y alone or N(1575)Y + L(1014)F were introduced into an Aedes aegypti sodium channel, AaNav1-1, and the mutants were functionally examined in Xenopus oocytes. N(1575)Y did not alter AaNav1-1 sensitivity to pyrethroids. However, the N(1575)Y + L(1014)F double mutant was more resistant to pyrethroids than the L(1014)F mutant channel. Further mutational analysis showed that N(1575)Y could also synergize the effect of L(1014)S/W, but not L(1014)G or other pyrethroid-resistant mutations in IS6 or IIS6. Computer modeling predicts that N(1575)Y allosterically alters PyR2 via a small shift of IIS6. Our findings provide the molecular basis for the coexistence of N(1575)Y with L(1014)F in pyrethroid resistance, and suggest an allosteric interaction between IIS6 and LIII/IV in the sodium channel.

Sodium channel point mutations associated with pyrethroid resistance in Chinese strains of Culex pipiens quinquefasciatus (Diptera: Culicidae)

BACKGROUND

Pesticide resistance due to sodium channel point mutations has been well documented in many mosquito species.

METHODS

We tested the resistance of six, wild, Chinese populations of the mosquito Culex pipiens quinquefasciatus to deltamethrin and cyhalothrin. The full length of the sodium channel gene was cloned and sequenced from a pooled sample of mosquitoes from each population.

RESULTS

Seven amino acid alterations were found (V250M, R436K, M943V, I973T, L1035F, L1035S and E1901D). Correlation between the frequencies of these mutations and the level of pesticide resistance (LC50) associated with them indicates that those at position L1035 (corresponding to position L1014F in the house fly, Musca domestica; GenBank Accession No.: X96668) are associated with resistance to deltamethrin and cyhalothrin. The frequency of the L1035F mutation was significantly correlated with resistance to deltamethrin (R2 = 0.536, P = 0.049) and cyhalothrin (R2 = 0.626, P = 0.030), and the combined frequency of the L1035F and L1035S mutations was significantly correlated with resistance to both deltamethrin (R2 = 0.661, P = 0.025), and cyhalothrin (R2 = 0.803, P = 0.008). None of the other mutations were correlated with either deltamethrin or cyhalothrin resistance. Interestingly, an HWE test indicated significant linkage between the M943V and I973T mutations (P < 0.01), but further research is required to determine the biological significance of this linkage.

CONCLUSIONS

Identifying these mutations may be of practical benefit to the development of pesticide resistance management programs.

A look into the Medical and Veterinary Entomology crystal ball

Medical and Veterinary Entomology (MVE) represents a leading periodical in its field and covers many aspects of the biology and control of insects, ticks, mites and other arthropods of medical and veterinary importance. Since the first issue of the journal, researchers working in both developed and developing countries have published in MVE, with direct impact on current knowledge in the field. An increasing number of articles dealing with the epidemiology and transmission of vector-borne pathogens have been published in MVE, reflecting rapid changes in vector distribution, pathogen transmission and host-arthropod interactions. This article represents a gaze into the crystal ball in which we identify areas of increasing interest, discuss the main changes that have occurred in the epidemiology of parasitic arthropods since the first issue of MVE, and predict the principal scientific topics that might arise in the next 25 years for scientists working in medical and veterinary entomology.

First molecular genotyping of voltage gated sodium channel alleles in Culex quinquefasciatus populations in Malaysia

A nationwide investigation was performed to detect the presence of 1014 mutation(s) in voltage gated sodium channel (kdr) gene of Culex quinquefasciatus from 14 residential areas across 13 states and a federal territory in Malaysia. Molecular genotyping of kdr mutation was performed via a modified three tubes allele-specific-polymerase chain reaction (AS-PCR) and direct sequencing of kdr gene. Based on the results of AS-PCR, homozygous susceptible (SS) genotype was found in nine out of 14 populations with 38 individuals from a total sample size of 140. Heterozygous (RS) genotype was most predominant (99 individuals) and distributed across all study sites. Homozygous resistance (RR) genotype was detected in Perak (one individual) and Selangor (two individuals). The resistance kdr allele frequencies ranged from 0.1 to 0.55, with the highest being detected in Cx. quinquefasciatus population from Selangor. This study has documented the first field-evolved instance of 1014F mutation in Malaysian mosquitoes and the findings of this study could be utilized in the implementation of strategic measures in vector control programs in Malaysia.

Diversity and Convergence of Sodium Channel Mutations Involved in Resistance to Pyrethroids

Pyrethroid insecticides target voltage-gated sodium channels, which are critical for electrical signaling in the nervous system. The intensive use of pyrethroids in controlling arthropod pests and disease vectors has led to many instances of pyrethroid resistance around the globe. In the past two decades, studies have identified a large number of sodium channel mutations that are associated with resistance to pyrethroids. The purpose of this review is to summarize both common and unique sodium channel mutations that have been identified in arthropod pests of importance to agriculture or human health. Identification of these mutations provides valuable molecular markers for resistance monitoring in the field and helped the discovery of the elusive pyrethroid receptor site(s) on the sodium channel.