Mitochondrial genetic differentiation across populations of the malaria vector Anopheles lesteri from China (Diptera: Culicidae)

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.

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.

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.

Detection of alleles associated with resistance to chemical insecticide in the malaria vector Anopheles arabiensis in Santiago, Cabo Verde

BACKGROUND

Mosquitoes of the Anopheles gambiae complex are the main malaria vectors worldwide. Due to the lack of a vaccine to prevent malaria, the principal way to reduce the impact of this disease relies on the use of chemical insecticides to control its vectors. However, the intensive use of such compounds has led to the emergence of insecticide resistance in several Anopheles populations in Africa. This study aimed to investigate the presence of resistance alleles in an Anopheles arabiensis population from the City of Praia, capital of the Archipelago Cabo Verde, one of the countries on the World Health Organization list of countries that are on a path to eliminate local transmission of malaria.

METHODS

Larvae from the Anopheles genus were collected using a one-pint dipper in three areas of City of Praia. Larvae were fed and maintained until the emergence of adult mosquitoes, and these were morphologically identified. In addition, molecular identification was performed using IGS markers and all An. arabiensis samples were subjected to PCR to screen for mutations associated to resistance in the Ace-1, Na_v and GSTE2 genes.

RESULTS

From a total of 440 mosquitoes collected, 52.3% were morphologically identified as An. gambiae sensu lato (s.l.) and 46.7% as Anopheles pretoriensis. The molecular identification showed that 100% of the An. gambiae s.l. were An. arabiensis. The mutations G119S in the Ace-1 gene and L119F in the GSTE2 gene were screened but not found in any sample. However, sequencing analysis for GSTE2 revealed the presence of 37 haplotypes, 16 polymorphic sites and a high genetic diversity (π = 2.67). The L1014S mutation in the Na_v (voltage-gated sodium channel gene) was detected at a frequency of 7.3%.

CONCLUSION

This is the first study to investigate the circulation of insecticide resistance alleles in An. arabiensis from Cabo Verde. The circulation of the L1014S allele in the population of An. arabiensis in the city of Praia suggests that pyrethroid resistance may arise, be quickly selected, and may affect the process of malaria elimination in Cabo Verde. Molecular monitoring of resistance should continue in order to guide the development of strategies to be used in vector control in the study region.

Mitochondrial variation in subpopulations of Anopheles balabacensis Baisas in Sabah, Malaysia (Diptera: Culicidae)

Anopheles balabacensis, the primary vector of Plasmodium knowlesi in Sabah, Malaysia, is both zoophilic and anthropophilic, feeding on macaques as well as humans. It is the dominant Anopheles species found in Kudat Division where it is responsible for all the cases of P. knowlesi. However there is a paucity of basic biological and ecological information on this vector. We investigated the genetic variation of this species using the sequences of cox1 (1,383 bp) and cox2 (685 bp) to gain an insight into the population genetics and inter-population gene flow in Sabah. A total of 71 An. balabacensis were collected from seven districts constituting 14 subpopulations. A total of 17, 10 and 25 haplotypes were detected in the subpopulations respectively using the cox1, cox2 and the combined sequence. Some of the haplotypes were common among the subpopulations due to gene flow occurring between them. AMOVA showed that the genetic variation was high within subpopulations as compared to between subpopulations. Mantel test results showed that the variation between subpopulations was not due to the geographical distance between them. Furthermore, Tajima's D and Fu's Fs tests showed that An. balabacensis in Sabah is experiencing population expansion and growth. High gene flow between the subpopulations was indicated by the low genetic distance and high gene diversity in the cox1, cox2 and the combined sequence. However the population at Lipasu Lama appeared to be isolated possibly due to its higher altitude at 873 m above sea level.

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.

Phylogenetic relationships and evolutionary history of the greater horseshoe bat, Rhinolophus ferrumequinum, in Northeast Asia

The greater horseshoe bat, Rhinolophus ferrumequinum, is an important model organism for studies on chiropteran phylogeographic patterns. Previous studies revealed the population history of R. ferrumequinum from Europe and most Asian regions, yet there continue to be arguments about their evolutionary process in Northeast Asia. In this study, we obtained mitochondrial DNA cyt b and D-loop data of R. ferrumequinum from Northeast China, South Korea and Japan to clarify their phylogenetic relationships and evolutionary process. Our results indicate a highly supported monophyletic group of Northeast Asian greater horseshoe bats, in which Japanese populations formed a single clade and clustered into the mixed branches of Northeast Chinese and South Korean populations. We infer that R. ferrumequinum in Northeast Asia originated in Northeast China and South Korea during a cold glacial period, while some ancestors likely arrived in Japan by flying or land bridge and subsequently adapted to the local environment. Consequently, during the warm Eemian interglaciation, the Korea Strait, between Japan and South Korea, became a geographical barrier to Japanese and inland populations, while the Changbai Mountains, between China and North Korea, did not play a significant role as a barrier between Northeast China and South Korea populations.

Historical survey of the kdr mutations in the populations of Anopheles sinensis in China in 1996-2014

Background

Anopheles sinensis has become an important malaria vector in China. The long-term extensive utilization of pyrethroids for ITNs and IRS for mosquito control in the last three decades has resulted in the occurrence of resistant An. sinensis populations in many regions. Knockdown resistance (kdr), caused by point mutations in the VGSC gene, is one of the mechanisms that confer resistance to DDT and pyrethroids. Recently, several investigations revealed the kdr occurrence in some An. sinensis populations, however, no kdr data were available earlier than 2009. A survey tracking the dynamics of the kdr mutations in past decades would provide invaluable information to understand how the kdr alleles spread in mosquito populations temporally and spatially.

Methods

A survey was conducted on the kdr alleles at condon 1014 of the VGSC gene and their distributions in 733 specimens of An. sinensis and 232 specimens of the other eight member species of the Anopheles hyrcanus group that were collected from 17 provinces in China in 1996–2014.

Results

A total of three kdr alleles, TTT (F), TTG (F) and TGT (C) were detected, and TGT (C) and TTT (F) were already present in the specimens from Jiangsu and Shandong as early as 1997. The TTT (F) was the most frequent mutant allele, and largely distributed in central China, namely Shandong, Jiangsu, Anhui, Henan, Shanghai, Jiangxi and Hubei. When data were analysed in three time intervals, 1996–2001, 2005–2009, 2010–2014, the prevalence of kdr alleles increased progressively over time in the populations in central China. In contrast, the kdr alleles were less frequent in the samples from other regions, especially in Yunnan and Hainan, despite the documented presence of pyrethroid resistant populations in those regions. Interestingly, no mutant alleles were detected in all 232 specimens of eight other species in the An. hyrcanus group.

Conclusion

The survey revealed that the kdr occurrence and accumulation in the An. sinensis populations were more frequent in central China than in the other regions, suggesting that the kdr mutations may contribute significantly to the pyrethroid resistance in the mosquitoes in central China.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0644-0) contains supplementary material, which is available to authorized users.

High morphological and genetic variabilities of Ochlerotatus scapularis, a potential vector of filarias and arboviruses

Background

Ochlerotatus scapularis is a potential vector of filarias and arboviruses in the Neotropics. This species was once typically associated with sylvatic environments; however, cases of synanthropy and urbanization of this species have been increasingly reported in southeast Brazil. Despite the medical relevance of Oc. scapularis, its populational variability is not yet known. To our knowledge, this is the first report describing the morphological and genetic variabilities of this species.

Methods

Population samples were characterized using the cytochrome oxidase subunit I (COI) mitochondrial gene and wing geometrics. Adult mosquitoes were collected from five sampling sites from remnants of the Atlantic forest embedded in the urban or rural areas of southeast Brazil.

Results

In the 130 individuals analyzed, 46 COI haplotypes were detected. Haplotype diversity was high and ranged from 0.66 to 0.97. Six haplotypes were present in 61% of the individuals, whereas the remaining haplotypes were less frequent (39%). Wing shape was also highly polymorphic. Differentiation of populations across sampling sites according to genetic distances (F_st = −0.009 to 0.060) and morphological distances (Q_st = 0.47) indicated that populations were not identical. No correlations were noted for phenetic and genetic diversities (p = 0.19) or for genetic or phenetic distances with geographical distances (p = 0.2 and p = 0.18, respectively).

Conclusions

Our study results suggest that Oc. scapularis has a rich genetic patrimony, even though its habitat is fragmented. Implications of such genetic richness with respect to vectorial competence, plasticity, and ability to exploit urbanized areas need to be further investigated.

Electronic supplementary material

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

Distinct population structure for co-occurring Anopheles goeldii and Anopheles triannulatus in Amazonian Brazil

To evaluate whether environmental heterogeneity contributes to the genetic heterogeneity in Anopheles triannulatus, larval habitat characteristics across the Brazilian states of Roraima and Pará and genetic sequences were examined. A comparison with Anopheles goeldii was utilised to determine whether high genetic diversity was unique to An. triannulatus. Student t test and analysis of variance found no differences in habitat characteristics between the species. Analysis of population structure of An. triannulatus and An. goeldii revealed distinct demographic histories in a largely overlapping geographic range. Cytochrome oxidase I sequence parsimony networks found geographic clustering for both species; however nuclear marker networks depicted An. triannulatus with a more complex history of fragmentation, secondary contact and recent divergence. Evidence of Pleistocene expansions suggests both species are more likely to be genetically structured by geographic and ecological barriers than demography. We hypothesise that niche partitioning is a driving force for diversity, particularly in An. triannulatus.

Diversity of Anopheles mosquitoes in Binh Phuoc and Dak Nong Provinces of Vietnam and their relation to disease

Background

Human malaria is still a burden in Dak Nong and Binh Phuoc Provinces in south-central Vietnam that border Cambodia. Several Anopheles species that transmit human malarial Plasmodium may also transmit Wuchereria bancrofti, the nematode that causes Bancroftian lymphatic filariasis. The objective of this study was to investigate the role of Anopheles species in the transmission of these two pathogens in the two highly malaria endemic provinces of Vietnam.

Methods

Anopheles mosquitoes were collected in Dak Nong and Binh Phuoc Provinces in November and December of 2010 and 2011. Human landing catches, paired collections on human and buffalo, and resting captures were made with mouth aspirators. Collections were also made with light traps. Morphological and PCR-based methods were used to identify the species. Real-time PCR was used to detect Plasmodium species and W. bancrofti in individual mosquitoes.

Results

Twenty-four Anopheles species were identified among 797 captured mosquitoes. Anopheles dirus was found in both provinces and was the predominant species in Binh Phuoc Province; An. maculatus was the most prevalent species in Dak Nong Province. Anopheles minimus was collected only in Binh Phuoc Province. Some specimens of An. minimus and An. pampanai were misidentified based on morphology. Four specimens of An. scanloni were identified, and this is the first report of this species of the Dirus Complex in Vietnam. Two females, one An. dirus and one An. pampanai, collected in Binh Phuoc Province were infected with P. vivax, for an overall infection rate of 0.41% (2/486): 0.28% for An. dirus (1/361) and 20% for An. pampanai (1/5). No mosquitoes were found to be infected with P. falciparum, P. knowlesi or W. bancrofti in either province.

Conclusion

A diversity of Anopheles species occurs in Dak Nong and Binh Phuoc Provinces of Vietnam, several of which are considered to be actual and potential vectors of malarial protozoa and microfilariae. It is highly likely that two of the species, An. dirus and An. pampanai, are active in malaria transmission based on the detection of P. vivax in females of these species. This is the first report of An. scanloni in Vietnam.

Taxonomic assessment of Anopheles crawfordi and An. dangi of the Hyrcanus Group of subgenus Anopheles in Vietnam

Anopheles dangi, introduced as a new species of the Hyrcanus Group of subgenus Anopheles in an illustrated dichotomous key for the identification of the Anopheles mosquitoes of Vietnam published in 1987, was distinguished from Anopheles crawfordi based on the presence of a humeral pale spot on the base of the costal vein of the wing. However, this character has been known to occur occasionally in An. crawfordi. To determine whether An. dangi is distinct from An. crawfordi, we analyzed nucleotide sequences of the COI, COII and Cyt-b genes of mtDNA and the D3 gene of rDNA obtained from specimens collected in south-central Vietnam that were identified as An. dangi and An. crawfordi based on the presence or absence, respectively, of a humeral pale spot. Maximum Likelihood and Bayesian analyses of the sequences showed a low mean genetic distance of 0.004 for specimens identified as An. crawfordi and 0.008 for those identified as An. dangi. The mean genetic distance between the two nominal species was 0.006, compared with 0.077 for any group versus the outgroup taxa Anopheles dirus and Anopheles minimus, and the specimens of the two forms clustered in a single strongly supported clade. Consequently, An. dangi is merely a morphological variant of An. crawfordi and is deemed to be a synonym of that nominal species.

Genetic variations of ND5 gene of mtDNA in populations of Anopheles sinensis (Diptera: Culicidae) malaria vector in China

Background

Anopheles sinensis is a principal vector for Plasmodium vivax malaria in most parts of China. Understanding of genetic structure and genetic differentiation of the mosquito should contribute to the vector control and malaria elimination in China.

Methods

The present study investigated the genetic structure of An. sinensis populations using a 729 bp fragment of mtDNA ND5 among 10 populations collected from seven provinces in China.

Results

ND5 was polymorphic by single mutations within three groups of An. sinensis that were collected from 10 different geographic populations in China. Out of 140 specimens collected from 10 representative sites, 84 haplotypes and 71 variable positions were determined. The overall level of genetic differentiation of An. sinensis varied from low to moderate across China and with a F_ST range of 0.00065 – 0.341. Genealogy analysis clustered the populations of An. sinensis into three main clusters. Each cluster shared one main haplotype. Pairwise variations within populations were higher (68.68%) than among populations (31.32%) and with high fixation index (F_ST = 0.313). The results of the present study support population growth and expansion in the An. sinensis populations from China. Three clusters of An. sinensis populations were detected in this study with each displaying different proportion patterns over seven Chinese provinces. No correlation between genetic and geographic distance was detected in overall populations of An. sinensis (R² = 0.058; P = 0.301).

Conclusions

The results indicate that the ND5 gene of mtDNA is highly polymorphic in An. sinensis and has moderate genetic variability in the populations of this mosquito in China. Demographic and spatial results support evidence of expansion in An. sinensis populations.