Landscape genetic structure and evolutionary genetics of insecticide resistance gene mutations in Anopheles sinensis

Wolbachia and mosquitoes: Exploring transmission modes and coevolutionary dynamics in Shandong Province, China

Vector-borne diseases leave a large footprint on global health. Notable culprits include West Nile virus (WNV), St. Louis encephalitis virus (SLEV), and Japanese encephalitis virus (JEV), all transmitted by Culex mosquitoes. Chemical insecticides have been widely used to reduce the spread of mosquito-borne diseases. Still, mosquitoes are becoming more and more resistant to most chemical insecticides which cause particular harm to the ecology. Wolbachia belongs to the family Ehrlichiaceae in the order Rickettsiales and is a matrilineally inherited endosymbiont present in 60% of insects in nature. Wolbachia is capable of inducing a wide range of reproductive abnormalities in its hosts, such as cytoplasmic incompatibility, and can alter mosquito resistance to pathogen infection. Wolbachia has been proposed as a biological alternative to chemical vector control, and specific research progress and effectiveness have been achieved. Despite the importance of Wolbachia, this strategy has not been tested in Culex pipiens pallens, the most prevalent mosquito species in Shandong Province, China. Little is known about how the mass release of Wolbachia-infected mosquitoes may impact the genetic structure of Culex pipiens pallens, and how the symbiotic bacterium Wolbachia interacts with mitochondria during host mosquito transmission. Based on the population genetic structure of Culex pipiens pallens in Shandong Province, this study investigated the infection rate and infection type of Wolbachia in Shandong Province and jointly analysed the evolutionary relationship between the host mosquito and the symbiotic bacterium Wolbachia. Our study showed that Wolbachia naturally infected by Culex pipiens pallens in Shandong Province was less homologous to Wolbachia infected by Aedes albopictus released from mosquito factory in Guangzhou. Our results also show that Culex pipiens pallens is undergoing demographic expansion in Shandong Province. The overall Wolbachia infection rate of Culex pipiens pallens was 92.8%, and a total of 15 WSP haplotypes were detected. We found that the genetic diversity of Wolbachia was low in Culex pipiens pallens from Shandong Province, and the mosquitoes were infected only with type B Wolbachia. Visualizing the relationship between Culex pipiens pallens and Wolbachia using a tanglegram revealed patterns of widespread associations. A specific coevolutionary relationship exists between the host mosquito and Wolbachia. Knowledge of this mosquito-Wolbachia relationship will provide essential scientific information required for Wolbachia-based vector control approaches in Shandong Province and will lead to a better understanding of the diversity and evolution of Wolbachia for its utility as a biocontrol agent.

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.

Susceptibility Status of Several Field-Collected German Cockroaches (Blattella germanica) to a Pyrethroid Insecticide and Molecular Detection of Knockdown Resistance (kdr)

Background

High frequency of insecticide used to control German cockroaches may lead to insecticide resistance development. We aimed to compare the level of insecticide resistance and heterogeneity in insecticide resistance of field-collected German cockroaches in eight selected zones of Mashhad City, Razavi Khorasan Province, Iran.

Methods

The present study was conducted from October 2019 to May 2021. Adult and nymphs of German cockroaches were collected from infested restaurants or hotels in eight zones of Mashhad City and then colonized in an insectarium. The cockroaches were subjected to bioassay against cypermethrin insecticide. In addition, the genomic DNA of each cockroach population was analyzed in the region where the kdr mutations reside in the German cockroach.

Results

The LT50 values against cypermethrin in the susceptible strain (SS) and eight field-collected strains so varied from 17.52 to 95.36 min. The resistance ratio of the strains was also different in response to the insecticide. The similarity of multiple sequence alignments at the amino acid level was 97.5%-100%. An exon 20 transversion mutation and a nonsynonymous substitution were found. The L1014F substitution was detected in 83.4% of the cockroach samples.

Conclusion

The field-collected strains were resistant to cypermethrin at different levels. Furthermore, the molecular study confirms the heterogeneity in the level of resistance among eight strains. Therefore, eight strains might have a different history of insecticide treatment.

Vector Biology and Integrated Management of Malaria Vectors in China

Malaria is an infectious disease caused by Plasmodium parasites, transmitted by Anopheles sinensis, Anopheles lesteri, Anopheles minimus, and Anopheles dirus in China. In 2021, the disease was eliminated in China after more than 70 years of efforts implementing an integrated mosquito management strategy. This strategy comprised indoor residual spray, insecticide-treated bed nets, irrigation management, and rice-fish coculture based on an understanding of taxonomic status and ecological behaviors of vector species, in conjunction with mass drug administration and promotion of public education. However, China still faces postelimination challenges, including the importation of approximately 2,000-4,000 cases of malaria into the country each year, as well as widespread resistance to pyrethroid insecticides in An. sinensis; these challenges require long-term vector surveillance to understand the distribution, population density, and development of resistance in vector mosquitoes to prevent local epidemics caused by imported malaria cases.

Genetic diversity and wing geometric morphometrics among four populations of Aedes aegypti (Diptera: Culicidae) from Benin

BACKGROUND

The impact of the arbovirus vector Aedes aegypti is of major concern for global public health as the viruses that it transmits affect millions of people each year worldwide. Originating in Africa, Ae. aegypti has now spread throughout much of the world. While the genetic makeup of Ae. aegypti in the New World has been extensively studied, there is limited knowledge on its genetic diversity in Africa, particularly at a microgeographical level.

METHODS

We investigated mitochondrial cytochrome oxidase I of four Ae. aegypti populations from Benin and employed wing morphometric analyses as a cost-effective and reliable tool to explore population structure. Our sampling encompassed various areas of Benin, from the southern to the northern borders of the country, and included urban, semi-urban, and sylvatic sites.

RESULTS

We observed a notable level of genetic diversity (haplotype diversity of 0.8333) and nucleotide diversity (0.00421986), and identified seven distinct haplotypes. Sylvatic and semi-urban sites exhibited a greater number of haplotypes compared to urban sites. Utilizing 18 wing landmarks, we calculated the centroid size, which revealed significant variation among the three landscape types. However, principal component analysis, employed to assess wing shape variation, did not demonstrate significant differences between populations based on landscape type.

CONCLUSIONS

Our findings indicate substantial genetic and morphological diversity among Ae. aegypti populations in Benin, and provide insight into important biological characteristics of these populations with respect to their potential to transmit viruses. To the best of our knowledge, this is the first study undertaken in Africa to integrate genetics with morphology to analyse the population structure of the major arbovirus vector Ae. aegypti.

High frequency of Voltage-gated sodium channel (VGSC) gene mutations in Aedes albopictus (Diptera: Culicidae) suggest rapid insecticide resistance evolution in Shanghai, China

BACKGROUND

Dengue fever is an infectious disease that is imported into Shanghai, China and requires prevention and control measures. Controlling the vector Aedes albopictus through insecticide use is a key approach to dengue control. However, the rapid evolution of insecticide resistance in Ae. albopictus has raised concerns about the failure of dengue control efforts. Knockdown resistance (kdr) caused by point mutations in the voltage-gated sodium channel (VGSC) gene is a primary mechanism of pyrethroid resistance. In this study, we investigated the kdr mutations of Ae. albopictus in Shanghai and evaluated the trend in its evolution.

METHODOLOGY/PRINCIPAL FINDINGS

We collected 17 populations of Ae. albopictus from 15 districts in Shanghai in 2020, extracted genomic DNA from individual mosquitoes, and amplified Domain II, III, and IV in VGSC using PCR. Following sequencing, we obtained 658 VGSC sequences. We detected the nonsynonymous mutations V1016G, I1532T, and F1534S/C/I, among which V1016G and F1534C/I were reported in Shanghai for the first time and F1534I was a novel mutant allele in Ae. albopictus. The overall mutation frequency was 84.65%, with individual mutation frequencies ranging from 46.81% to 100%, excluding the Fengxian District population, which had a frequency of 0%. The V1016G and I1532T mutation types accounted for 7.14% and 3.42%, respectively. The mutant allele at codon 1534 accounted for 63.98% of all mutations, including TCC/S (62.77%), TGC/C (1.06%), and ATC/I (0.15%). We identified and classified five intron types in Domain III by length, including A (83 bp, 12.07%), B (68 bp, 87.30%), C (80 bp, 0.16%), D (72 bp, 0.16%), and E (70 bp, 0.31%). Individuals with intron B had a significant mutation tendency at codon 1534 relative to intron A (chi-square test, p < 0.0001). We found no correlation between mutation frequency and the amount of pyrethroid used (Pearson correlation, p = 0.4755).

CONCLUSIONS/SIGNIFICANCE

In recent years, kdr mutations in the Ae. albopictus population in Shanghai have rapidly evolved, as evidenced by an increase in mutation types and significantly increased mutation frequency. The F1534I/ATC mutant allele was found to be a novel mutation, F1534C/TGC was reported for the first time in Shanghai, and intron B in Domain III was significantly associated with mutation frequency at codon 1534. Continuous monitoring of resistance changes and strict regulation of insecticide use are required.

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

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.

Genetic structure of Aedes albopictus (Diptera: Culicidae) populations in China and relationship with the knockdown resistance mutations

BACKGROUND

Mosquito control is needed to prevent dengue fever, which is mainly spread by Aedes albopictus in China. Application of insecticides is one of the main mosquito control methods; however, this approach can fail due to the knockdown resistance (kdr) gene mutation that causes decreased sensitivity to insecticides in Ae. albopictus. The kdr mutation patterns among different regions in China differ significantly. However, the underlying mechanism and factors that influence kdr mutation remain unclear. To explore the potential influence of genetic background on the development of insecticide resistance in Ae. albopictus, we analyzed the genetic structure of Ae. albopictus populations in China and its correlation with major kdr mutations.

METHODS

We collected Ae. albopictus from 17 sites in 11 provinces (municipalities) across China from 2016 to 2021 and extracted the genomic DNA from individual adult mosquitoes. We selected eight microsatellite loci for genotyping, and based on microsatellite scores, we estimated intraspecific genetic diversity, population structure, and effective population size. The association between the intrapopulation genetic variation and F1534 mutation rate was evaluated by the Pearson correlation coefficient.

RESULTS

Based on variation analysis of the microsatellite loci of 453 mosquitoes representing 17 populations throughout China, more than 90% of the variation occurred within individuals, whereas only about 9% of the variation occurred among populations, indicating that field populations of Ae. albopictus are highly polymorphic. The northern populations tended to belong to gene pool I (BJFT 60.4%, SXXA 58.4%, SDJN 56.1%, SXYC 46.8%), the eastern populations tended to belong to pool III (SH 49.5%, JZHZ 48.1%), and the southern populations tended to belong to three different gene pools. Moreover, we observed that the greater the fixation index (F_ST), the lower the wild-type frequency of F1534 of VSGC.

CONCLUSIONS

The degree of genetic differentiation among Ae. albopictus populations in China was low. These populations were divided into three gene pools, in which the northern and eastern pools are relatively homogeneous, while the southern gene pool is heterogeneous. The possible correlation between its genetic variations and kdr mutations is also noteworthy.

Evaluation of the evolutionary genetics and population structure of Culex pipiens pallens in Shandong province, China based on knockdown resistance (kdr) mutations and the mtDNA-COI gene

BACKGROUND

Mosquitoes are important vectors for a range of diseases, contributing to high rates of morbidity and mortality in the human population. Culex pipiens pallens is dominant species of Culex mosquito in northern China and a major vector for both West Nile virus and Bancroftian filariasis. Insecticide application were largely applied to control the mosquito-mediated spread of these diseases, contributing to increasing rates of resistance in the mosquito population. The voltage-gated sodium channel (Vgsc) gene is the target site of pyrethroids, and mutations in this gene cause knockdown resistance (kdr). While these kdr mutations are known to be critical to pyrethroid resistance, their evolutionary origins remain poorly understood. Clarifying the origins of these mutations is potential to guide further vector control and disease prevention efforts. Accordingly, the present study was designed to study the evolutionary genetics of kdr mutations and their association with the population structure of Cx. p. pallens in Shandong province, China.

METHODS

Adult Culex females were collected from Shandong province and subjected to morphological identification under a dissection microscope. Genomic DNA were extracted from the collected mosquitoes, the Vgsc gene were amplified via PCR and sequenced to assess kdr allele frequencies, intron polymorphisms, and kdr codon evolution. In addition, population genetic diversity and related population characteristics were assessed by amplifying and sequencing the mitochondrial cytochrome C oxidase I (COI) gene.

RESULTS

Totally, 263 Cx. p. pallens specimens were used for DNA barcoding and sequencing analyses to assess kdr allele frequencies in nine Culex populations. The kdr codon L1014 in the Vgsc gene identified two non-synonymous mutations (L1014F and L1014S) in the analyzed population. These mutations were present in the eastern hilly area and west plain region of Shandong Province. However, only L1014F mutation was detected in the southern mountainous area and Dongying city of Shandong Province, where the mutation frequency was low. Compared to other cities, population in Qingdao revealed significant genetic differentiation. Spatial kdr mutation patterns are likely attributable to some combination of prolonged insecticide-mediated selection coupled with the genetic isolation of these mosquito populations.

CONCLUSIONS

These data suggest that multiple kdr alleles associated with insecticide resistance are present within the Cx. p. pallens populations of Shandong Province, China. The geographical distributions of kdr mutations in this province are likely that the result of prolonged and extensive insecticide application in agricultural contexts together with frequent mosquito population migrations. In contrast, the low-frequency kdr mutation detected in central Shandong Province populations may originate from the limited selection pressure in this area and the relative genetic isolation. Overall, the study compares the genetic patterns revealed by a functional gene with a neutral marker and demonstrates the combined impact of demographic and selection factors on population structure.

Impact of MiraNet® long-lasting insecticidal net against Anopheles arabiensis wild population of Northern Tanzania

Despite high levels of pyrethroid resistance reported in malaria vectors, long-lasting insecticidal nets (LNs) still play a key role in controlling malaria transmission. This study tested the efficacy of MiraNet®, a pyrethroid-based LN against a wild population of Anopheles arabiensis in northern Tanzania. DuraNet® was used as a positive control in this evaluation. Standard WHO laboratory bioefficacy evaluations of MiraNet and DuraNet that were unwashed or had been washed 20 times indicated optimal knockdown and mortality for both net types against a susceptible strain of Anopheles gambiae s.s. Standard experimental hut evaluations were conducted to evaluate the efficacy of both nets against a wild population of An. arabiensis. The killing effect of MiraNet was 54.5% for unwashed and 50% for 20 times washed while DuraNet achieved 44.4% mortality for unwashed and 47.4% for 20 times washed against wild An. arabiensis. Both DuraNet and MiraNet exhibited significantly higher killing effects (> 44.4%). There was no significant difference in deterrence or induced exophily detected between the treatment arms for either net. Additionally, there were no adverse effects reported among hut sleepers. The results of this study indicate that the pyrethroid net MiraNet can be used effectively against wild populations of An. gambiae s.l. of low to moderate resistant levels from Northern Tanzania.

Genetic Diversity and Population Structure of Anopheles funestus in Western Kenya Based on Mitochondrial DNA Marker COII

The mitochondrial marker, COII, was employed to assess the genetic structure and diversity of Anopheles funestus, a very important malaria vector in Africa that adapt and colonize different ecological niches in western Kenya. Mosquitoes were collected using mechanical aspirators in four areas (Bungoma, Port Victoria, Kombewa, and Migori) in western Kenya. Following morphological identification, PCR was used to confirm the species. The COII gene was amplified, sequenced, and analyzed to determine genetic diversity and population structure. A total of 126 (Port Victoria-38, Migori-38, Bungoma-22, and Kombewa-28) sequences of COII were used for population genetic analysis. Anopheles funestus had a high haplotype diversity (Hd = 0.97 to 0.98) but low nucleotide diversity (Π = 0.004 to 0.005). The neutrality test revealed negative Tajima's D and Fs values indicating an excess of low-frequency variation. This could be attributed to either population expansion or negative selection pressure across all the populations. No genetic or structural differentiation (Fst = -0.01) and a high level of gene flow (Gamma St, Nm = 17.99 to 35.22) were observed among the populations. Population expansion suggests the high adaptability of this species to various ecological requirements, hence sustaining its vectorial capacity and malaria transmission.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

In the course of malaria elimination in the Greater Mekong Subregion (GMS), malaria epidemiology has experienced drastic spatiotemporal changes with residual transmission concentrated along international borders and the rising predominance of Plasmodium vivax. The emergence of Plasmodium falciparum parasites resistant to artemisinin and partner drugs renders artemisinin-based combination therapies less effective while the potential spread of multidrug-resistant parasites elicits concern. Vector behavioral changes and insecticide resistance have reduced the effectiveness of core vector control measures. In recognition of these problems, the Southeast Asian International Center of Excellence for Malaria Research (ICEMR) has been conducting multidisciplinary research to determine how human migration, antimalarial drug resistance, vector behavior, and insecticide resistance sustain malaria transmission at international borders. These efforts allow us to comprehensively understand the ecology of border malaria transmission and develop population genomics tools to identify and track parasite introduction. In addition to employing in vivo, in vitro, and molecular approaches to monitor the emergence and spread of drug-resistant parasites, we also use genomic and genetic methods to reveal novel mechanisms of antimalarial drug resistance of parasites. We also use omics and population genetics approaches to study insecticide resistance in malaria vectors and identify changes in mosquito community structure, vectorial potential, and seasonal dynamics. Collectively, the scientific findings from the ICEMR research activities offer a systematic view of the factors sustaining residual malaria transmission and identify potential solutions to these problems to accelerate malaria elimination in the GMS.

Population genetic structure and evolutionary genetics of Anopheles sinensis based on knockdown resistance (kdr) mutations and mtDNA-COII gene in China-Laos, Thailand-Laos, and Cambodia-Laos borders

BACKGROUND

Vector control is still a pivotal method for preventing malaria, and its potency is weakened by the increasing resistance of vectors to chemical insecticides. As the most abundant and vital malaria vector in Southeast Asia, the chemical insecticide resistance status in Anopheles sinensis remains elusive in Laos, which makes it imperative to evaluate the true nature of chemical insecticide resistance-associated genetic mutations in An. sinensis in Laos.

METHODS

Adult An. sinensis were collected from three border regions in Laos. DNA was extracted from individual mosquitoes. PCR amplification and DNA sequencing of a fragment containing codon 1014 of the voltage-gated sodium channel (vgsc) gene were completed to study the kdr allele frequency distribution, kdr intron polymorphism, population genetic diversity, and the evolutionary status of the kdr codon. The mitochondrial cytochrome c oxidase subunit II gene (COII) was amplified and sequenced to examine population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

Nine wild kdr haplotypes of the vgsc gene were detected in this study, and eight of them, namely 1014L1, 1014L2, 1014L4, 1014L7, 1014L9, 1014L10, 1014L11, and 1014L21, were discovered in the China-Laos border (northern Laos), while 1014L3 was only detected in the Thailand-Laos border (northwestern Laos) and Cambodia-Laos border (southern Laos). The newly identified haplotype, 1014L21, was uniquely distributed in the China-Laos border and was not identified in other countries. Based on sequence analysis of the mitochondrial COII genes, significant genetic differentiation and limited gene flow were detected between the China-Laos and Cambodia-Laos An. sinensis populations, which suggested that those two regions were genetically isolated. The distinct distribution of the kdr haplotype frequencies is probably the result of geographical isolation in mosquito populations.

CONCLUSIONS

Lack of kdr mutations in the vgsc gene was probably due to genetic isolation and the absence of intense selection pressure in the three border regions of Laos. This study reveals that pyrethroid-based chemical insecticides are still appropriate for battling An. sinensis in parts of Laos, and routine monitoring of chemical insecticide resistance should be continuously implemented and focused on more restricted areas as part of chemical insecticide resistance management.

Community structure and insecticide resistance of malaria vectors in northern-central Myanmar

BACKGROUND

Myanmar is one of the six countries in the Greater Mekong Subregion (GMS) of Southeast Asia. Malaria vectors comprise many Anopheles species, which vary in abundance and importance in malaria transmission among different geographical locations in the GMS. Information about the species composition, abundance, and insecticide resistance status of vectorial systems in Myanmar is scarce, hindering our efforts to effectively control malaria vectors in this region.

METHODS

During October and November 2019, larvae and adult females of Anopheles mosquitoes were collected in three sentinel villages of Banmauk township in northern Myanmar. Adult female mosquitoes collected by cow-baited tent collection (CBTC) and adults reared from field-collected larvae (RFCL) were used to determine mortality rates and knockdown resistance (kdr) against deltamethrin using the standard WHO susceptibility test. Molecular species identification was performed by multiplex PCR and ITS2 PCR, followed by DNA sequencing. The kdr mutation at position 1014 of the voltage-gated sodium channel gene was genotyped by DNA sequencing for all Anopheles species tested.

RESULTS

A total of 1596 Anopheles mosquitoes from seven morphologically identified species groups were bioassayed. Confirmed resistance to deltamethrin was detected in the populations of An. barbirostris (s.l.), An. hyrcanus (s.l.), and An. vagus, while possible resistance was detected in An. annularis (s.l.), An. minimus, and An. tessellatus. Anopheles kochi was found susceptible to deltamethrin. Compared to adults collected by CBTC, female adults from RFCL had significantly lower mortality rates in the four species complexes. A total of 1638 individuals from 22 Anopheles species were molecularly identified, with the four most common species being An. dissidens (20.5%) of the Barbirostris group, An. peditaeniatus (19.4%) of the Hyrcanus group, An. aconitus (13.4%) of the Funestus group, and An. nivipes (11.5%) of the Annularis group. The kdr mutation L1014F was only detected in the homozygous state in two An. subpictus (s.l.) specimens and in a heterozygous state in one An. culicifacies (s.l.) specimen.

CONCLUSIONS

This study provides updated information about malaria vector species composition and insecticide resistance status in northern Myanmar. The confirmed deltamethrin resistance in multiple species groups constitutes a significant threat to malaria vector control. The lack or low frequency of target-site resistance mutations suggests that other mechanisms are involved in resistance. Continual monitoring of the insecticide resistance of malaria vectors is required for effective vector control and insecticide resistance management.

Spatial heterogeneity of knockdown resistance mutations in the dengue vector Aedes albopictus in Guangzhou, China

Background

The city of Guangzhou has been the epicenter of dengue fever in China since the 1990s, with Aedesalbopictus being the primary vector. The main method used to control vectors and prevent dengue fever has been the application of chemical insecticides; however, this control strategy has resulted in the development of resistance to these insecticides in mosquitoes. Here we report our investigation of the patterns of knockdown resistance (kdr) mutations in 15 field populations of Ae.albopictus collected from 11 districts in Guangzhou.

Results

Four mutant alleles (V1016G, F1534S, F1534C, F1534L) were detected in domain II and III of the voltage-gated sodium channel (VGSC) gene. Various allele frequencies of kdr mutations were observed (3.1–25.9% for V1016G, 22.6–85.5% for F1534S, 0–29.0% for F1534L, 0.6–54.2% for F1534C). Seven kdr haplotypes (VF, VS, VL, VC, GF, GC, GS) were identified; the highest frequency of haplotypes was found for the single mutant haplotype VS (50.8%), followed by the wild-type VF haplotype (21.7%) and the single mutant haplotype VC (11.9%). Of the three double mutant haplotypes, GF was the most frequent (8.8%), followed by GC (1.2%) and GS (0.8%). Aedesalbopictus showed spatial heterogeneity in deltamethrin resistance in populations collected in Guangzhou. We also observed significant differences in haplotype frequency. The frequency of the VC haplotype was significantly higher in high-risk dengue areas than in low-risk ones.

Conclusions

The kdr allele V1016G was discovered for the first time in Guangzhou. Genetic isolation in mosquito populations and long-term insecticide selection seem to be responsible for the persistent, patchy distribution of kdr mutant alleles. The small-scale spatial heterogeneity in the distribution and frequency of kdr mutations may have important implications for vector control operations and insecticide resistance management strategies.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05241-7.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

Identification and phylogenetic analysis of voltage-gated sodium channel haplotypes in the malaria vector Anopheles sinensis using a high-throughput amplicon sequencing approach

Background

Anopheles sinensis is a dominant vector for malaria transmission in Asian countries. Voltage-gated sodium channel (VGSC) mutation-mediated knock-down resistance (kdr) has developed in many A. sinensis populations because of intensive and long-term use of pyrethroids. Our previous study showed that multiple mutations at position 1014 of the VGSC were heterogeneously distributed in A. sinensis populations across Sichuan, China.

Methods

To understand resistance genotypes at the haplotype level and reconstruct the phylogenetic relationship of VGSC haplotypes, a cost-effective next-generation sequencing (NGS)-based amplicon sequencing approach was established to clarify haplotypes containing codon 1014 of the VGSC gene from a total of 446 adults collected in 12 locations of Sichuan, China.

Results

Nineteen (19) haplotypes were identified, including 11 wild 1014L, 6 resistance 1014F, and 2 resistance 1014C haplotypes. We found that resistance haplotypes of A. sinensis VGSC were widely distributed at frequencies ranging from 3.67 to 92.61%. The frequencies of the 1014C haplotype in the southeast of Sichuan (Luzhou, Guangan, and Suining) were relatively higher than those in other sampling locations. Phylogenetic analyses support that kdr-type mutation at position 1014 is not singly originated and resistance 1014C haplotypes evolve from TTT-encoding 1014F.

Conclusions

A cost-effective next-generation sequencing (NGS)-based amplicon sequencing approach has been established in this study. The data revealed the patchy distribution of VGSC resistance haplotypes with overall high frequencies in Sichuan, China. Phylogenetic analyses support multiple origins and sequential evolution (1014L → 1014F → 1014C) for kdr-type mutations in A. sinensis.

Graphical abstract

ASDB: A comprehensive omics database for Anopheles sinensis

Anopheles sinensis is a key disease vector for human malaria and parasitic diseases such as malayan filariasis, and it is considered to be one of the most important malaria vectors in China and Southeast Asia. As high-throughput sequencing and assembly technology are widely used in An. sinensis, a lot of omics data have been generated, and abundant genome, mRNA transcriptome, miRNA transcriptome and resequencing results have been accumulated. In addition, lots of valuable morphological images and publications have been produced with the in-depth studies on An. sinensis. However, the increased quantity, variety, and structure complexity of the omics data create inconveniences for researchers to use and manage this information. We have built an An. sinensis omics database (ASDB, http://asdb.jungleran.com/) - a comprehensive and integrated database to promote scientific research on An. sinensis. Docker was used to deploy a development environment and Drupal to build ASDB. ASDB provides a Blast tool to do sequence alignment of genome sequence, gene sequence and protein sequence of An. sinensis. It also offers JBrowse (a next-generation genome visualization and analysis web platform) to facilitate researchers visualize the gene structure, non-coding RNA (include miRNA, snRNA, tRNA and so on) structure and genomic variation sites as desired. ASDB has integrated various latest omics data of An. Sinensis, including de novo genome and its annotation data, genome variation data (such as SNP and InDel), transcriptome and its expression value, miRNA expression value and miRNA-mRNA interaction, metagenomes. The database has also included the morphological images of different developmental stages and tissues, and important literatures associated with An. sinensis. ASDB provides a user-friendly search and displays pages. The integration of these resources will contribute to the study of basic biology and functional genome of An. sinensis.

Influence of a Major Mountainous Landscape Barrier (Mount Cameroon) on the Spread of Metabolic (GSTe2) and Target-Site (Rdl) Resistance Alleles in the African Malaria Vector Anopheles funestus

Increased levels of insecticide resistance in major malaria vectors such as Anopheles funestus threaten the effectiveness of insecticide-based control programmes. Understanding the landscape features impacting the spread of resistance makers is necessary to design suitable resistance management strategies. Here, we examined the influence of the highest mountain in West Africa (Mount Cameroon; 4095 m elevation) on the spread of metabolic and target-site resistance alleles in An. funestus populations. Vector composition varied across the four localities surveyed along the altitudinal cline with major vectors exhibiting high parity rate (80.5%). Plasmodium infection rates ranged from 0.79% (An. melas) to 4.67% (An. funestus). High frequencies of GSTe2^R (67–81%) and Rdl^R (49–90%) resistance alleles were observed in An. funestus throughout the study area, with GSTe2^R frequency increasing with altitude, whereas the opposite is observed for Rdl^R. Patterns of genetic diversity and population structure analyses revealed high levels of polymorphisms with 12 and 16 haplotypes respectively for GSTe2 and Rdl. However, the reduced diversity patterns of resistance allele carriers revealed signatures of positive selection on the two genes across the study area irrespective of the altitude. Despite slight variations associated with the altitude, the spread of resistance alleles suggest that control strategies could be implemented against malaria vectors across mountainous landscapes.

Genetic structure of deltamethrin-resistant populations of Triatoma infestans (Hemiptera: Reduviidae) in the Gran Chaco

The genetic structure of natural populations offers insight into the complexities of their dynamics, information that can be relevant to vector control strategies. Microsatellites are useful neutral markers to investigate the genetic structure and gene flow in Triatoma infestans, one of the main vectors of Chagas disease in South America. Recently, a heterogeneous pyrethroid-resistant hotspot was found in the Argentine Gran Chaco, characterized by the highest levels of deltamethrin resistance found at the present time. We applied population genetics analyses to microsatellite and village data and search for associations between the genetic variability and the heterogeneous toxicological pattern previously found. We genotyped 10 microsatellite loci in 67 T. infestans from 6 villages with no, low, and high pyrethroid resistance. The most genetically diverse populations were those susceptible or with low values of resistance. In contrast, high-resistance populations had lower herozygosity and some monomorphic loci. A negative association was found between variability and resistant ratios. Global and pairwise F_STs indicated significant differentiation between populations. The only susceptible population was discriminated in all the performed studies. Low-resistance populations were also differentiated by a discriminant analysis of principal components (DAPC) and were composed mostly by the same two genetic clusters according to STRUCTURE Bayesian algorithm. Individuals from the high-resistance populations were overlapped in the DAPC and shared significant proportions of a genetic cluster. These observations suggest that the resistant populations might have a common origin, although more genetic markers and samples are required to test this hypothesis more rigorously.

Tracing temporal and geographic distribution of resistance to pyrethroids in the arboviral vector Aedes albopictus

BACKGROUND

The arboviral vector Aedes albopictus became established on all continents except Antarctica in the past 50 years. A consequence of its rapid global invasion is the transmission of diseases previously confined to the tropics and subtropics occurring in temperate regions of the world, including the re-emergence of chikungunya and dengue in Europe. Application of pyrethroids is among the most widely-used interventions for vector control, especially in the presence of an arboviral outbreak. Studies are emerging that reveal phenotypic resistance and monitor mutations at the target site, the para sodium channel gene, primarily on a local scale.

METHODS

A total of 512 Ae. albopictus mosquitoes from twelve geographic sites, including those from the native home range and invaded areas, were sampled between 2011 and 2018, and were analyzed at five codons of the para sodium channel gene with mutations predictive of resistance phenotype. Additionally, to test for the origin of unique kdr mutations in Mexico, we analyzed the genetic connectivity of southern Mexico mosquitoes with mosquitoes from home range, the Reunion Island, America and Europe.

RESULTS

We detected mutations at all tested positions of the para sodium channel gene, with heterozygotes predominating and rare instance of double mutants. We observed an increase in the distribution and frequency of F1534C/L/S mutations in the ancestral China population and populations in the Mediterranean Greece, the appearance of the V1016G/I mutations as early as 2011 in Italy and mutations at position 410 and 989 in Mexico. The analyses of the distribution pattern of kdr alleles and haplotype network analyses showed evidence for multiple origins of all kdr mutations.

CONCLUSIONS

Here we provide the most-up-to-date survey on the geographic and temporal distribution of pyrethroid-predictive mutations in Ae. albopictus by combining kdr genotyping on current and historical samples with published data. While we confirm low levels of pyrethroid resistance in most analyzed samples, we find increasing frequencies of F1534C/S and V1016G in China and Greece or Italy, respectively. The observed patterns of kdr allele distribution support the hypothesis that on site emergence of resistance has contributed more than spread of resistance through mosquito migration/invasions to the current widespread of kdr alleles, emphasizing the importance of local surveillance programs and resistance management.

Culex quinquefasciatus Egg Membrane Alteration and Ovicidal Activity of Cipadessa baccifera (Roth) Plant Extracts Compared to Synthetic Insect Growth Regulators

Background

Insecticide resistance among mosquito vectors for synthetic insecticides still remains a major problem for control efforts. This study assessed the ovicidal potential of crude solvent extracts from the medicinal plant Cipadessa baccifera comparatively to standard registered synthetic insect growth regulators (IGR) on freshly laid eggs of Culex quinquefasciatus.

Method

Five plant extracts were prepared using different solvents. The batches of eggs were exposed to different concentrations of each solvent extract comparatively to synthetic IGR. The hatched eggs of Cx. quinquefasciatus were subjected to different concentrations. The first instars that emerged from the eggs were counted daily. The egg hatching inhibition was observed 24, 48 and 72 hrs post treatment. The desiccation median time (DT₅₀ and DT₉₀) was calculated.

Results

The percent egg hatching inhibition was inversely proportional to the concentration of extracts. The morphological damage to the eggs was observed. Among five solvent extracts, acetone extracts showed the highest ovicidal activity. The changes in eggshell morphology were observed. The maximum ovicidal activity was observed in acetone extracts with DT₅₀ value of 1.70 hrs (0.91–2.22). The methanol plant extract using gas chromatography-mass spectrometry identified 14 compounds.

Conclusion

These results suggest that the acetone extracts of C. baccifera have the potential to be used as an ovicidal agent for controlling mosquito populations in aquatic stages. The biodegradability of the extracts has the advantage of being eco-friendly.

Phylogeny of certain members of Hyrcanus group (Diptera: Culicidae) in China based on mitochondrial genome fragments

Background

Species of the Anopheles hyrcanus group are widely distributed in Palearctic and Oriental regions and some of them are important malaria vectors. The cryptic species of An. hyrcanus group was almost impossible to identify based only on their morphology. The phylogenetic relationship of An. hyrcanus group was also not clear.

Methods

Five members of An. hyrcanus group were identified by rDNA ITS2 sequencing as An. yatsushiroensis, An. belenrae, An. kleini, An. lesteri and An. sineroides. The mitochondrial genome fragments were sequenced and annotated using the mitochondrial genome of An. sinensis as reference. Based on the four segments and Joint Data sequences of these species, and other four anopheline species downloaded from GenBank, intraspecific as well as interspecific genetic distances were calculated and the phylogenetic trees were reconstructed by the methods of neighbor joining, maximum parsimony, minimum evolution and maximum likelihood.

Findings

Four parts of mitochondrial genomes, which were partial fragments COI + tRNA + COII (F5), ATP6 + COIII(F7 + F8), ND1(F19) and lrRNA (F21), were obtained. All fragments were connected as one sequence (referred as Joint Data), which had a total length of 3393 bp. All fragment sequences were highly conservative within species, with the maximum p distance (0.026) calculated by F19 of An. belenrae. The pairwise interspecific p distance calculated by each fragment showed minor or even no difference among An. sinensis, An. kleini and An. belenrae. However, interspecific p distances calculated by the Joint Data sequence ranged from 0.004 (An. belenrae vs An. kleini) to 0.089 (An. sineroides vs An. minimus), and the p distances of the six members of An. hyrcanus group were all less than 0.029. The phylogenetic tree showed two major clades: all subgenus Anopheles species (including six members of An. hyrcanus group, An. atroparvus and An. quadrimaculatus A) and subgenus Cellia (including An. dirus and An. minimus). The An. hyrcanus group was divided into two clusters as ((An. lesteri, An. sineroides) An. yatsushiroensis) and ((An. belenrae, An. sinensis) An. kleini)).

Conclusions

The An. hyrcanus group in this study could be divided into two clusters, in one of which An. belenrae, An. sinensis and An. kleini were most closely related. More molecular markers would make greater contribution to phylogenetic analysis.

Insecticide Resistance Status and Mechanisms of Anopheles sinensis (Diptera: Culicidae) in Wenzhou, an Important Coastal Port City in China

Although scaled-up interventions and effective control efforts have drastically reduced malaria morbidity and mortality, malaria remains a serious threat to public health worldwide. Anopheles sinensis Wiedemann 1828 is a historically important vector of Plasmodium vivax (Haemosporida: Plasmodiidae) malaria in China. Insecticide resistance has become a major obstacle to vector-borne disease control. However, little is known about the insecticide resistance of An. sinensis in Wenzhou, an important coastal port city in Zhejiang province, China. The aim of this study was to examine insecticide resistance and mechanisms in An. sinensis field mosquito populations. Evidence of multiple insecticide resistance was found in An. sinensis adult female populations. Medium to high frequencies of target site kdr together with fixed ace-1 mutations was detected in both the Ruian and Yongjia populations. Both populations showed an association between kdr L1014 mutation and resistance phenotype when tested against deltamethrin and DDT. Significantly different metabolic enzyme activities were found between the susceptible laboratory strain and field-collected mosquitoes from both Ruian and Yongjia. Both field collected An. sinensis populations exhibited significantly higher P450 enzyme activity compared with the laboratory strain, while the field-collected resistant mosquitoes exhibited various GST and COE enzyme activities. These results indicate multiple resistance mechanisms in An. sinensis field populations. Effective implementation of insecticide resistance management strategies is urgently needed. The data collected in this study will be valuable for modeling insecticide resistance spread and vector-control interventions.

Target-site mutations (AChE-G119S and kdr) in Guangxi Anopheles sinensis populations along the China-Vietnam border

Background

In South Asia, the epidemiology of malaria is complex, and transmission mainly occurs in remote areas near international borders. Vector control has been implemented as a key strategy in malaria prevention for decades. A rising threat to the efficacy of vector control efforts is the development of insecticide resistance, thus it is important to monitor the type and frequency of insecticide resistant alleles in the disease vectors such as An. sinensis along the China-Vietnam border. Such information is needed to synthesize effective malaria vector control strategies.

Methods

A total of 208 adults of An. sinensis, collected from seven sites in southwest Guangxi along the China-Vietnam border, were inspected for the resistance-conferring G119S mutation in acetylcholinesterase (AChE) by PCR-RFLP (polymerase chain reaction restriction fragment length polymorphism) and kdr mutations in the voltage-gated sodium channel (VGSC) by sequencing. In addition, the evolutionary origin of An. sinensis vgsc gene haplotypes was analyzed using Network 5.0.

Results

The frequencies of mutant 119S of AChE were between 0.61–0.85 in the seven An. sinensis populations. No susceptible homozygote (119GG) was detected in three of the seven sites (DXEC, LZSK and FCGDX). Very low frequencies of kdr (0.00–0.01) were detected in the seven populations, with most individuals being susceptible homozygote (1014LL). The 1014F mutation was detected only in the southeast part (FCGDX) at a low frequency of 0.03. The 1014S mutation was distributed in six of the seven populations with frequencies ranging from 0.04 to 0.08, but absent in JXXW. Diverse haplotypes of 1014L and 1014S were found in An. sinensis along the China-Vietnam border, while only one 1014F haplotype was detected in this study. Consistent with a previous report, resistant 1014S haplotypes did not have a single origin.

Conclusions

The G119S mutation of AChE was present at high frequencies (0.61–0.85) in the An. sinensis populations along the China-Vietnam border, suggesting that the vector control authorities should be cautious when considering carbamates and organophosphates as chemicals for vector control. The low frequencies (0.00–0.11) of kdr in these populations suggest that pyrethroids remain suitable for use against An. sinensis in these regions.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

Phylogeography of Diaphorina citri (Hemiptera: Liviidae) and its primary endosymbiont, 'Candidatus Carsonella ruddii': an evolutionary approach to host-endosymbiont interaction

BACKGROUND

In insects, little is known about the co-evolution between their primary endosymbionts and hosts at the intraspecific level. This study examined co-diversification between the notorious agricultural pest Diaphorina citri and its primary endosymbionts (P-endosymbiont), 'Candidatus Carsonella ruddii' at the population level.

RESULTS

Maximum likelihood, haplotype network, principal components and Bayesian clustering identified three lineages for D. citri and its P-endosymbiont: a Western clade containing individuals from Pakistan, Bhutan (Phuentsholing), Vietnam (Son La), USA, Myanmar and China (Ruili, Yunnan); a Central clade, with accessions originating from Southwest China, Bhutan (Tsirang) and Bangladesh; and an Eastern clade containing individuals from Southeast Asia, and East and South China. A more diverse genetic structure was apparent in the host mitochondrial DNA than their P-endosymbionts; however, the two sets of data were strongly congruent.

CONCLUSION

This study provides evidence for the co-diversification of D. citri and its P-endosymbiont during the migration from South Asia to East and Southeast Asia. We also suggest that the P-endosymbiont may facilitate investigations into the genealogy and migration history of the host. The biogeography of D. citri and its P-endosymbiont indicated that D. citri colonized and underwent a secondary dispersal from South Asia to East and Southeast Asia. © 2018 Society of Chemical Industry.

Sequencing and analysis of the complete mitochondrial genome in Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis (Diptera: Culicidae) is a primary vector of Plasmodium vivax and Brugia malayi in most regions of China. In addition, its phylogenetic relationship with the cryptic species of the Hyrcanus Group is complex and remains unresolved. Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, of which mitochondrial genome data of An. sinensis is not available.

Methods

An. sinensis samples was collected from Shandong, China, and identified by molecular marker. Genomic DNA was extracted, followed by the Illumina sequencing. Two complete mitochondrial genomes were assembled and annotated using the mitochondrial genome of An. gambiae as reference. The mitochondrial genomes sequences of the 28 known Anopheles species were aligned and reconstructed phylogenetic tree by Maximum Likelihood (ML) method.

Findings

The length of complete mitochondrial genomes of An. sinensis was 15,076 bp and 15,138 bp, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and an AT-rich control region. As in other insects, most mitochondrial genes are encoded on the J strand, except for ND5, ND4, ND4L, ND1, two rRNA and eight tRNA genes, which are encoded on the N strand. The bootstrap value was set as 1000 in ML analyses. The topologies restored phylogenetic affinity within subfamily Anophelinae. The ML tree showed four major clades, corresponding to the subgenera Cellia, Anopheles, Nyssorhynchus and Kerteszia of the genus Anopheles.

Conclusions

The complete mitochondrial genomes of An. sinensis were obtained. The number, order and transcription direction of An. sinensis mitochondrial genes were the same as in other species of family Culicidae.

Electronic supplementary material

The online version of this article (10.1186/s40249-017-0362-7) contains supplementary material, which is available to authorized users.

Spotlight on Plasmodium falciparum evolutionary system in the southeastern Atlantic forest

Abstract Malaria elimination is now set to occur in Brazil until 2030. While this achievement is feasible, as it is for other endemic regions worldwide, it is important to recognize resistance of parasites and vectors against anti-malarial interventions. Resistance against drugs and insecticides can lead to discontinuities of malaria transmission, known as residual malaria transmission. Herein, we described a novel phenomenon that is occurring in a residual malaria transmission scenario in the southeastern Atlantic forest. This novel phenomenon does not belong to what is known and therefore we decided to explain it based on an evolutionary perspective. Although it shall not be viewed as a threat to public health, the phenomenon has important aspects that should be highlighted. Specifically, it represents an adaptation of P. falciparum among vectors and hosts in the southeastern Atlantic forest. Knowledge about this phenomenon could be of importance, including to the on-going malaria elimination programs.

2b-RAD genotyping for population genomic studies of Chagas disease vectors: Rhodnius ecuadoriensis in Ecuador

BACKGROUND

Rhodnius ecuadoriensis is the main triatomine vector of Chagas disease, American trypanosomiasis, in Southern Ecuador and Northern Peru. Genomic approaches and next generation sequencing technologies have become powerful tools for investigating population diversity and structure which is a key consideration for vector control. Here we assess the effectiveness of three different 2b restriction site-associated DNA (2b-RAD) genotyping strategies in R. ecuadoriensis to provide sufficient genomic resolution to tease apart microevolutionary processes and undertake some pilot population genomic analyses.

METHODOLOGY/PRINCIPAL FINDINGS

The 2b-RAD protocol was carried out in-house at a non-specialized laboratory using 20 R. ecuadoriensis adults collected from the central coast and southern Andean region of Ecuador, from June 2006 to July 2013. 2b-RAD sequencing data was performed on an Illumina MiSeq instrument and analyzed with the STACKS de novo pipeline for loci assembly and Single Nucleotide Polymorphism (SNP) discovery. Preliminary population genomic analyses (global AMOVA and Bayesian clustering) were implemented. Our results showed that the 2b-RAD genotyping protocol is effective for R. ecuadoriensis and likely for other triatomine species. However, only BcgI and CspCI restriction enzymes provided a number of markers suitable for population genomic analysis at the read depth we generated. Our preliminary genomic analyses detected a signal of genetic structuring across the study area.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that 2b-RAD genotyping is both a cost effective and methodologically simple approach for generating high resolution genomic data for Chagas disease vectors with the power to distinguish between different vector populations at epidemiologically relevant scales. As such, 2b-RAD represents a powerful tool in the hands of medical entomologists with limited access to specialized molecular biological equipment.

Insecticide resistance in malaria vectors along the Thailand-Myanmar border

BACKGROUND

There is a paucity of data about the susceptibility status of malaria vectors to Public Health insecticides along the Thailand-Myanmar border. This lack of data is a limitation to guide malaria vector-control in this region. The aim of this study was to assess the susceptibility status of malaria vectors to deltamethrin, permethrin and DDT and to validate a simple molecular assay for the detection of knock-down resistance (kdr) mutations in the study area.

METHODS

Anopheles mosquitoes were collected in four sentinel villages during August and November 2014 and July 2015 using human landing catch and cow bait collection methods. WHO susceptibility tests were carried out to measure the mortality and knock-down rates of female mosquitoes to deltamethrin (0.05%), permethrin (0.75%) and DDT (4%). DNA sequencing of a fragment of the voltage-gated sodium channel gene was carried out to identify knock-down resistance (kdr) mutations at position 1014 in mosquitoes surviving exposure to insecticides.

RESULTS

A total of 6295 Anopheles belonging to ten different species were bioassayed. Resistance or suspected resistance to pyrethroids was detected in An. barbirostris (s.l.) (72 and 84% mortality to deltamethrin (n = 504) and permethrin (n = 493) respectively), An. hyrcanus (s.l.) (33 and 48% mortality to deltamethrin (n = 172) and permethrin (n = 154), respectively), An. jamesii (87% mortality to deltamethrin, n = 111), An. maculatus (s.l.) (85 and 97% mortality to deltamethrin (n = 280) and permethrin (n = 264), respectively), An. minimus (s.l.) (92% mortality, n = 370) and An. vagus (75 and 95% mortality to deltamethrin (n =148) and permethrin (n = 178), respectively). Resistance or suspected resistance to DDT was detected in An. barbirostris (s.l.) (74% mortality, n = 435), An. hyrcanus (s.l.) (57% mortality, n = 91) and An. vagus (97% mortality, n = 133). The L1014S kdr mutation at both heterozygous and homozygous state was detected only in An. peditaeniatus (Hyrcanus Group).

CONCLUSION

Resistance to pyrethroids is present along the Thailand-Myanmar border, and it represents a threat for malaria vector control. Further investigations are needed to better understand the molecular basis of insecticide resistance in malaria vectors in this area.

Genetic diversity and population structure of the primary malaria vector Anopheles sinensis (Diptera: Culicidae) in China inferred by cox1 gene

Background

Anopheles sinensis is a primary vector for Plasmodium vivax malaria in most regions of China. A comprehensive understanding of genetic variation and structure of the mosquito would be of benefit to the vector control and in a further attempt to contribute to malaria elimination in China. However, there is only inadequate population genetic data pertaining to An. sinensis currently.

Methods

Genetic variations and structure among populations of An. sinensis was examined and analyzed based on the nucleotide sequences of a 662 nt variable region of the mitochondrial cox1 gene among 15 populations from 20 collection sites in China.

Results

A total of 453 individuals in 15 populations were analyzed. The cox1 gene sequences were aligned, and 247 haplotypes were detected, 41 of these shared between populations. The range of haplotype diversity was from 0.709 (Yunnan) to 0.998 (Anhui). The genealogic network showed that the haplotypes were divided into two clusters, cluster I was at a high level of homoplasy, while cluster II included almost all individuals from the Yunnan population. The Yunnan population displayed a significantly high level of genetic differentiation (0.452−0.622) and a restricted gene flow with other populations. The pairwise F_ST values among other populations were lower. The AMOVA result showed that the percentage of variation within populations (83.83%) was higher than that among populations (16.17%). Mantel test suggested that geographical distance did not significantly contribute to the genetic differentiation (R² = 0.0125, P = 0.59). Neutral test and mismatch analysis results showed that the An. sinensis population has undergone demographic expansions.

Conclusions

Anopheles sinensis populations showed high genetic polymorphism by cox1 gene. The weak genetic structure may be a consequence of low genetic differentiation and high gene flow among populations, except the Yunnan samples. The Yunnan population was isolated from the other populations, gene flow limited by geographical distance and barriers. These findings will provide a theoretical basis for vector surveillance and vector control in China.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2013-z) contains supplementary material, which is available to authorized users.

A neural network prediction of environmental determinants of Anopheles sinensis knockdown resistance mutation to pyrethroids in China

Selection pressure caused by long-term intensive use of insecticides is the key driving force in resistance development. Additional parameters such as environmental conditions may affect both the mosquito response to insecticides and the selection of resistance mechanisms. In this context, we analyzed the environmental determinants of kdr prevalence in Anopheles sinensis across China. We collected kdr frequency from 48 sites across central and southern China, together with key environmental factors including long-term climatic data, topographic features, main crops, and land cover types. Trend surface analysis found that the distribution of kdr frequency can be partitioned into three regions, namely central China (kdr frequency >80%), western China (kdr frequency varies from 0% to 60%), and southern China (kdr frequency <10%). Seven predictor variables were selected based on a radial basis function neural network model. A multilayer perceptron (MLP) network model revealed that the number of crops in a year was the most important predictor for the kdr mutation rate. Topography, long-term mean climate and land cover all contributed to the kdr mutation rate. The observed mean kdr frequency was 53.0% and the MLP network model-predicted mean was 52.6%, a 0.1% relative error. Predicted kdr frequencies closely matched the observed values. The model explained 92% of the total variance in kdr frequency. The results indicated that kdr was associated with the intensity of pesticide usage. Crop cultivation information, together with environmental factors, may well predict the spatial heterogeneity of kdr mutations in An. sinensis in China.

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.