- Open Access
Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border
Malaria Journal volume 21, Article number: 91 (2022)
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.
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.
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 (FST = 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.
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.
Dengue fever and malaria continue to be critical parasitic diseases jeopardizing people’s lives in subtropical and tropical regions, primarily in the Greater Mekong Subregion (GMS, Myanmar, Vietnam, Laos, China, Thailand, and Cambodia) . Over the past 15 years, Cambodia has effectively controlled malaria through the introduction of artemisinin-based combinations, the establishment of rapid diagnostic tests (RDTs), the implementation of a village malaria health worker system, and the widespread use of long-lasting insecticidal bed nets (LLINs) [2,3,4,5].
Though its incidence declines significantly, malaria continuously poses a major burden on public health in this nation, in which there is a co-occurrence of Plasmodium vivax and Plasmodium falciparum. However, P. vivax control is obviously less effective . Cases of P. vivax have major distributions crossing six northeastern Cambodia provinces (20,954, ~ 80%), and particularly prevail in Stung Treng within the Cambodia-Laos border (~ 28% of overall vivax case identified) .
The transmission of malaria in the GMS is characterized by vector diversity and a high degree of spatial heterogeneity in the distribution pattern . Generally, Anopheles minimus, Anopheles dirus complex, and Anopheles sinensis refer to the major vector, whereas the significance pertaining to the respective species in the malaria-transmitting process changes extensively with areas . Some thorough investigations were conducted in the GMS, unfortunately, most of these researches failed to screen for “secondary” vectors (e.g., Anopheles nivipes) in terms of Plasmodium infection [10,11,12]. Anopheles nivipes was found to be Plasmodium positive according to a study in Bangladesh in 2012  and has also been reported positive for Plasmodium parasites in India and other countries [14, 15]. In Cambodia, An. nivipes was a secondary vector, as well as Anopheles philippinensis  in transmitting malaria and had been long suspected.
On the other hand, though mosquito population exhibiting various gene makeups is likely to be different in vector competence, rare data has been achieved regarding population genetics of An. nivipes in the Cambodia-Laos border. Data regarding genetic diversity and population structure can help develop available mosquito control plans [16, 17]. Owing to its significant sensitivity in the field of system and phylogenetic relationships for maternal inheritance and haploidization , Mitochondrial DNA (mtDNA) has been widely used as an effective genetic marker in describing molecular taxonomy, phylogenetic relationships, population structure, and genetic diversity in malaria vectors .
Given the existing drastic variations of malaria conditions in Cambodia, this work aimed to study the mosquito diversity and molecular phylogeny of secondary vectors i.e., An. nivipes in the Cambodia-Laos border. Siem Pang County (Stung treng Province) was taken as the investigation site and mosquito surveillance was firstly performed. The cytochrome c oxidase subunit 2 gene (COII) of mtDNA from An. nivipes samples were subsequently sequenced. Furthermore, the gene variability and phylogenetic relationship of An. nivipes along the Cambodia-Laos border were investigated and compared with other samples from India, Thailand, and Myanmar for examining their gene relations. This work largely aimed at gaining more insights into the continuous surveillance of mosquito vectors, as well as molecular phylogeny and evolution of An. nivipes in the Cambodia-Laos border by using an effective genetic marker, mtDNA-COII.
Stung treng borders Laos in the north and is located on the east bank of the Mekong River. This endemic region takes up nearly 11,000 square kilometres and has about 42,000 residents. The average population density is 3.82 people per square kilometre. It has dense forests, diversified landforms, the temperature changes greatly and the annual rainfall is 1800 mm.
Mosquito collection and identification
Adult mosquitoes were collected in two villages in Siem Pang County by overnight trapping with the battery-operated CDC light traps (Model: 1012, Origin: John W. Hock Inc, USA) in the cattle/pig pens or human rooms from 8:00 pm to 08:00 am and continued for 1–5 nights, in 2018. Four CDC light traps were operating in two cattle/pig pens and two human rooms each night. All live adult mosquitoes were killed by freezing in the refrigerator and subsequently isolated and distinguished according to subgroup, species, and sex through dissecting microscope based on standards procedures in the field office of Siem Pang County [20, 21]. All mosquitoes were initially morphologically sorted out in the field using the keys of Das et al. . Each morphologically identified specimen was kept individually in capped plastic beam capsules having silica gel and stored at 4 °C for molecular species confirmation and further processing. To avoid any deviations in the further analysis, molecular identification of An. nivipes based on COII was subsequently performed. Final species confirmation is required to have ≥ 98% sequence identity to the voucher specimens/sequences in the NCBI nt databases. To avoid the issue that COII alone did not produce significant results for the voucher sequence, the amplification and sequencing on ITS2 was also performed. As a result, ITS2 and COII database comparisons of each morphologically identified An. nivipes samples were paired to determine the species of An. nivipes (Additional file 1: Table S1).
DNA extraction and sequencing
The extraction of An. nivipes genomic DNA in individual mosquitoes was carried out following the producer’s manual (QIAamp® DNA Mini Kit, Germany). The amplification for approximately 620 bp of COII gene was carried out using primers, LYS-R (5′-ACTTGCTTTCAGTCATCTAATG-3′) and LEU-F (5′-TCTAATATGGCAGATTAGTGCA-3′). The overall PCR reaction volume reached 20 µl and the mixture of PCR reagent comprised 2 µl of DNA, 0.05 unit of Takara Taq (Dalian, China), 0.3 µM of the respective primer, 0.2 mM of dNTPs, and 2.5 µl of 10× Buffer (15 mM MgCl2, 100 mM Tris–HCl PH=8.3 as well as 500 mM KCl). Besides, the cycling parameter included 95 °C, 5 min; 95 °C/1 min, 51 °C/1 min, 72 °C/2 min for 35 cycles; with a final extension of 72 °C for 10 min. The PCR products were analysed by 1.5% agarose gel electrophoresis stained with GoldView (Solarbio, China), under UV transillumination. The sequencing reaction proceeded in both directions using an ABI Big Dye Terminator Kit (Applied Biosystems, Thermo Fisher Scientific) and was analysed through ABI Prism 3500xL Genetic Analysis Tool (Applied Biosystems, Thermo Fisher Scientific) in Shanghai (Sangon Biotech).
Sequence alignment and phylogenetic analysis
In order to analyse the sequence variations and genetic relationship of An. nivipes from Cambodia-Laos and neighbouring nations, multi-sequence alignment of COII sequence was firstly carried out after the retrieval of sequence deposited in GenBank by using the keywords “(species name) & COII” (http://www.ncbi.nlm.nih.gov/). The manual adjustment was conducted by using BioEdit V7.0.9 if required . Gaps were excluded from the analysis and characters were unweighted. Subsequently, a phylogenetic tree was generated using a Neighbour-joining algorithm bootstrapped with 1000 replicates  based on MEGA-X . The visualization of this phylogram was performed with the use of Figtree v1.4.2 . Anopheles lindesayi (GenBank Accession No. AJ620898) and Anopheles claviger (GenBank accession nos. AY129232 and DQ229313) were used as outgroup taxa to the An. nivipes, following previous studies .
The average of nucleotide differences per site (K), nucleotide diversity (π), haplotype diversity (Hd), and haplotypes (H) was calculated by using DnaSP v.5.0 . For more widely comparing haplotypes in Cambodia-Laos and other geographical regions, the existing data in GenBank from China, India, and Thailand were analysed, and a parsimony framework was subsequently carried out by using Network 4.0 .
The calculation of pairwise FST for estimating population differentiation by complying with a difference in haplotype frequency, Nei’s Nm estimated gene flow conformed to GST , as well as analysis of molecular variance (AMOVA) for determining the distribution of genetic variation in population using 1000 permutations were performed using Arlequin v.3.5 .
Moreover, a spatial analysis of molecular variance (SAMOVA 2.0) was conducted for clustering COII sequences to homogeneous populations with genetic and geographical homogeneousness . SAMOVA generates F-statistics (FCT, FST, FSC), with the use of the AMOVA approach, to K group for maximizing group variations. SAMOVA estimate was calculated for K = 2–8, with 1000 simulated annealing steps from each of 100 sets of initial starting conditions. Furthermore, isolation by distance (IBD) was examined using a nonparametric Mantel with the web-based computer program IBDWS v.3.16 .
Using the statistics D  and Fu’s Fs , the hypothesis of strict neutrality was tested in An. nivipes population based on DnaSP v.5.0. The mismatch distribution (simulated in Arlequin v.3.5) was carried out for distinguishing smooth unimodal distribution from multimodal or ragged distribution [36,37,38]. For rejecting the demographic expansion hypothesis, the difference with statistical significance in distribution under the observation and simulation was assessed with the sum of square deviations (SSD).
Surveillance of mosquito abundance and composition
A total of 1440 adult mosquitoes were collected in cattle/pig pen or human rooms through overnight trapping with the battery-operated CDC light traps (Fig. 1a). Those mosquitoes were classified into 5 genera including Anopheles, Culex, Aedes, Armigeres, Tripteroides, and 27 species. The predominant genera included Culex (9 species) and Anopheles (12 species), respectively taking up 60.07% (865/1440) and 31.25% (450/1440); Other genera including Aedes, Armigeres, and Tripteroides respectively accounted for 3.61% (52/1440), 2.71% (39/1440), and 0.07% (1/1440); In addition, undefined genera/species accounted for 2.29% (33/1400) (Fig. 1b, Additional file 2: Table S2). Among the Genus Anopheles, An. nivipes, and An. maculatus were two major species and accounted for 73.56% (331/450) and 14.22% (64/450), respectively, while Anopheles peditaeniatus, Anopheles argyropus, Anopheles tessellatus, Anopheles barbirostris, Anopheles sandaicus, Anopheles vagus, Anopheles interruptus, An. dirus, Anopheles kochi, and An. sinensis accounted for merely 0.22% (1/450), 0.22% (1/450), 1.33% (6/450), 1.11% (5/450), 0.44% (2/450), 2.00% (9/450), 0.22% (1/450), 1.33% (6/450), 4.00% (18/450), and 1.33% (6/450) (Fig. 1b, Additional file 2: Table S2).
Genetic diversity of Anopheles nivipes among populations
DNA samples were extracted from 53 mosquitos which were collected in Siem Pang and morphologically/molecularly identified as An. nivipes (Additional file 1: Table S1). The amplification of the COII sequence was achieved in these 53 specimens, representing populations from the Cambodia-Laos border. To draw a broader comparison in haplotype from Cambodia-Laos and other geographical regions, the available data in GenBank from neighbouring nations were downloaded and analysed (Additional file 3: Table S3). In general, a total of 86 An. nivipes COII sequences were generated for 10 populations, including Siem Pang (n = 53), Tripura (India; n = 4), Cheema (Punjab, India; n = 2), Bathinda (Punjab, India; n = 2), Nagaland (India; n = 6), Assam (India; n = 7), Myanmar (n = 2), Thailand (n = 2), Mizoram (India; n = 2) and Meghalaya (India; n = 6) (Table 1, Additional file 4: Fig. S1). A high level of nucleotide diversity (Pi = 0.0156), number of haplotypes (h = 47), and haplotype diversity (Hd = 0.9237) were observed in COII populations (Table 1).
Population structure and genetic differentiation
The median-joining network based on 86 COII sequences denoted the distribution pattern exhibited by 47 haplotypes in An. nivipes population. The An. nivipes populations fell into three Groups. Group 1 consisted of all the haplotypes except for Bathinda and Cheema; Group 2 consisted of haplotype from Cheema (Punjab, Northwest India); Group 3 with haplotype from Bathinda (Punjab, Northwest India). The most common haplotypes referred to H1 (n = 19), H2 (n = 13), and H6 (n = 8), as only identified in 73.58% (39/53) of Siem Pang and 50% (1/2) of Thailand. H21 (n = 1), H22 (n = 1) was only identified in Cheema; H20 (n = 1) and H23 (n = 1) were only identified in Bathinda (Fig. 2, Additional file 5: Fig. S2). In addition, a single individual denoted other 38 haplotypes but derived from haplotype 1 through a very few mutation steps. The UPGMA dendrogram based on K2P genetic distances between haplotypes indicated that H21, H22 constituted one cluster, H20 and H23 constituted another, while the other haplotypes constituted the third (Additional file 6: Fig. S3).
AMOVA analysis based on COII sequences demonstrated that most of the variances were found among group variation (86.72%) rather than within populations (4.56%) and among populations within groups (8.72%), suggesting that these populations could fall to several groups. The fixation index among groups (FCT), among populations within groups (FSC), and within populations (FST) showed statistical significance (P < 0.05) (Table 2). The neighbour-joining dendrogram indicated that the populations of Siem Pang (Stung treng, Cambodia), Tripura (Northeast India), Nagaland (Northeast India), Assam (Northeast India), Myanmar, Thailand, Mizoram (Northeast India), and Meghalaya (Northeast India) were clustered in one clade, whereas Cheema (Punjab, Northwest India) and Bathinda (Punjab, Northwest India) formed the second and third clades (Fig. 3). Anopheles lindesayi and An. claviger were introduced as outgroups which clustered into another clade.
The maximal level of genetic differentiation in accordance with the fixation index FST based on sequences analysis of COII was between Cheema and Siem Pang (FST = 0.97824, P < 0.05). Estimates of gene flow (Nm) varied extensively in populations, ranging from 0.01112 to 8.44874. The minimal was between Cheema and Siem Pang (Nm = 0.01112). The maximal was between Assam and Myanmar (Nm = 8.44874) (Additional file 7: Table S4).
Spatial genetic structure analysis
According to the COII-based SAMOVA analysis, the clustering number was between 2 and 9. High genetic differentiation between groups was detected, as the FCT value ranged from 0.74758 to 0.867 and was statistically significant. FCT value was maximal at K = 3 (0.867; P < 0.05) with Fsc (0.65541, P < 0.001). In this case, two Punjab populations (Bathinda and Cheema) became separated and assigned to two single population groups, leaving other populations, which demonstrated that there is a significant genetic barrier to restrict gene flow from Punjab to other populations. As K value increased up to K = 8, Siem Pang and Thailand remained attributed to one group, whereas the other population created 8 distinct groups (Table 3). In general, the patterns identified in SAMOVA analyses complied with the phylogenetic trees (Fig. 3). The UPGMA dendrogram based on Nei’s unbiased genetic distance between populations revealed three distinctive groups, as Bathinda and Cheema constituting two groups, and the remaining populations covered in a third group (Additional file 8: Fig. S4). A Mantel test revealed a significant correlation between geographical and genetic distances in all populations (Z = 44,983.1865, r = 0.5575, P = 0.0070), suggesting that the genetic structure observed in An. nivipes population could be partially explained by distance isolation based on COII sequence analysis (Fig. 4).
Demographic history and neutrality test
As indicated from Tajima’s D and Fu’s Fs tests based on COII, the Siem Pang population exhibited significant negativity (P < 0.05, P < 0.02), suggesting a recent population expansion or selection (Table 1). The observed smooth and unimodal mismatch distribution suggested a sudden population expansion, conforming to the mismatch distribution derived under the model of sudden expansion (Additional file 9: Fig. S5).
The spread of malaria in GMS exhibits the diversity of vector species and the high spatial heterogeneity of distribution patterns . In 2002 about 61% of Cambodia’s total land area was estimated to be forested  with over 80% of that in malaria-endemic areas . These forested areas are primarily located in provinces bordering Vietnam, Laos, and Thailand. Efficient forest malaria vectors are found, including An. dirus and An. minimus, so a high risk of malaria is imposed on people living in villages on the edge of the forest or engaged in forest activities [41,42,43,44]. Besides the main forest vector, considerable other Anopheles species appear close to human settlements. Under high coverage of insecticide-treated nets, transmission by “secondary” vectors with outdoor or early biting behaviour may be more important than transmission by primary vectors , since the secondary vectors are less anthropogenic and may be more prone to exophagous and early bites. Furthermore, the exogenous incubation period of P. vivax is highlighted to be shorter, and the secondary vector is likely to be a more effective vector as compared with P. falciparum . There is a need to study the extent to which malaria vectors in Cambodia are outdoor bites and early bites, and to assess the importance of secondary vectors.
Several studies have been conducted to screen the secondary or suspected vector in the GMS. In central Laos (Khammouane), An. nivipes was one of the predominant species and accounted for 11.55% , while in the south-eastern part of Laos (Nongceng), instead of the major vector An. minimus and An. dirus, An. nivipes is suspected to be the dominant vector, taking up over 65% . In Thailand, An. nivipes, as well as An. philippinensis, were suspected vectors and accounted for 30.18% and 27.32%, respectively . Anopheles nivipes referred to a secondary vector, as well as the An. philippinensis in Cambodia . According to St Laurent B’s research , An. nivipes referred to one of the five most prevalent species and accounted for 23.49% in Preah Vihear and 35.61% in Ratanakiri Provinces, both located in the north Cambodia and border with Stung treng Province. Likewise, in this work, An. nivipes accounted for 22.99% (331/1440) in Siem Pang County of Stung treng (Fig. 1, Additional file 2: Table S2). It is noteworthy that An. nivipes has been reported positive for Plasmodium parasites in Bangladesh, India, and other countries [13,14,15]. Besides, the vivax malaria case was largely reported in three northern provinces, Preah Vihear, Ratanakiri, and Stung treng, suggesting that An. nivipes might be a potential and important vector of P. vivax in northern Cambodia .
In mosquito barcoding studies, the most frequently used molecular markers or barcoding region is the cytochrome c oxidase subunit I (COI) of the mitochondrial genome (mtDNA), followed by the internal transcriptional spacer 2 (ITS2) of ribosomal DNA (rDNA). Since mtDNA is typically inherited maternally, any hybrid or offspring would only have the mtDNA from the maternal species . Furthermore, mitochondrial genomes evolve 5–10 times faster than nuclear genomes, making mtDNA potentially more useful than rDNA in correctly identifying recently differentiated species [51, 52]. Some researchers used COI or COII as the only marker to identify mosquito species and study their molecular evolution [16, 53,54,55,56,57,58]. Additionally, due to the typically high copy numbers, availability of conserved primer binding sequences, and ease of amplification , variation in mitochondrial DNA (mtDNA) is being used more recently to approach levels of population structure and genetic diversity within species . Nonetheless, more reliable species information can be obtained if multiple molecular markers are used simultaneously [60,61,62]. Based on mtDNA-COII sequencing in this work, the AMOVA analysis indicated that all populations could be divided into different groups (Table 2). According to phylogenetic analysis, the An. nivipes could fall to two Northwest India groups (Bathinda, Cheema) and a Northeastern India/SEA group (Figs. 2, 3). As revealed from strong genetic differentiation and limited gene flow between the Cambodia-Laos and Northwest India populations, these regions might undergo genetic isolation (Additional file 7: Tables S4). SAMOVA analysis further suggested a genetic barrier restricting gene flow between these populations (Table 3), and the patterns were also consistent with the phylogenetic trees. The Bay of Bengal, in the northeast part of the Indian Ocean, is situated between India to the west and northwest, Bangladesh to the north, and Myanmar to the east, acting as a natural barrier to interrupt the spread of malaria vector. On the contrary, the weak genetic differentiation and frequent gene flow between the Cambodia-Laos and other populations such as Northeast India, Myanmar, and Thailand suggest that these regions are not genetically isolated. Moreover, the Mantel test revealed a strong correlation of genetic and geographical distances, which demonstrated the significant isolation by the distance of An. nivipes (Fig. 4). Unimodal plots of the mismatch distribution were observed in the Cambodia-Laos population (Additional file 9: Fig. S5), as well as the Tajima’s D and Fu’s Fs test were both negative and significant, thereby demonstrating a recent expansion after a bottleneck of An. nivipes population in this border region (Table 1).
In the present study, the molecular phylogeny analysis of An. nivipes based on ITS2 was also conducted to evaluate its effectiveness compared to mtDNA-COII (Additional file 10). Specific to ITS2, no genetic barriers were identified between the Cambodia-Laos and other populations. SAMOVA analysis and Mantel test revealed no correlation between genetic and geographical distances, thereby demonstrating no isolation-by-distance was found in An. nivipes populations based on ITS2 (Additional file 10). Since ITS2 regions alone have been used in distinguishing closely related mosquito species belonging to various genera such as Anopheles , Culex , and Aedes , mtDNA-COII may be a more effective marker than rDNA-ITS2 for describing the genetic diversity and population structure of An. nivipes due to its advantages of maternal inheritance, no recombination, and high variability .
This work reported that a high diversity of mosquito vectors was found and Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. A genetic barrier limiting gene flow between Cambodia-Laos and two Northwest India populations had been found based on sequences analysis of mtDNA-COII, and a recent population expanding or selecting of An. nivipes occurred in this border, suggesting that mtDNA-COII can serve as effective markers for describing the genetic diversity and population structure of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The datasets generated and/or analysed during the current study are available in the GenBank (http://www.ncbi.nlm.nih.gov/). The raw datasets used and/or analysed during this study are available from the corresponding author upon reasonable request.
Greater Mekong subregion
Rapid diagnostic tests
Long-lasting insecticide-treated bed nets
Internal transcribed spacer 2
Cytochrome c oxidase subunit 2
Analysis of molecular variance
Spatial analysis of molecular variance
Isolation by distance
WHO. World malaria report 2016. Geneva: World Health Organization; 2016. https://www.who.int/malaria/publications/world-malaria-report-2016/report/en/. Accessed 20 Aug 2020.
Cox J, Dy Soley L, Bunkea T, Sovannaroth S, Soy Ty K, Ngak S, et al. Evaluation of community-based systems for the surveillance of day three-positive Plasmodium falciparum cases in Western Cambodia. Malar J. 2014;13:282.
Hasegawa A, Yasuoka J, Ly P, Nguon C, Jimba M. Integrating child health services into malaria control services of village malaria workers in remote Cambodia: service utilization and knowledge of malaria management of caregivers. Malar J. 2013;12:292.
Kheang ST, Duong S, Olkkonen A. Increasing access to early malaria diagnosis and prompted treatment in remote Cambodian villages. Am J Public Health. 2011;101:e6-8.
Maude RJ, Nguon C, Ly P, Bunkea T, Ngor P, de la Torre SEC, et al. Spatial and temporal epidemiology of clinical malaria in Cambodia 2004–2013. Malar J. 2014;13:385.
WHO. World malaria report 2019. Geneva: World Health Organization; 2019. https://www.who.int/malaria/publications/world-malaria-report-2019/report/en/. Accessed 20 Aug 2020.
Siv S, Roca-Feltrer A, Vinjamuri SB, Bouth DM, Lek D, Rashid MA, et al. Plasmodium vivax malaria in Cambodia. Am J Trop Med Hyg. 2016;95(6 Suppl):97–107.
WHO. Anopheline species complexes in South and South-East Asia. SEARO Technical Publication No. 57. 2007.
Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, et al. Malaria in the Greater Mekong Subregion: heterogeneity and complexity. Acta Trop. 2012;121:227–39.
Durnez L, Coosemans M. Residual transmission of malaria: an old issue for new approaches. In: Anopheles mosquitoes, new insights into malaria vectors. Manguin S (Ed). IntechOpen Publ. 2013:671–704. http://www.intechopen.com/articles/show/title/residual-transmission-of-malaria-an-old-issue-for-newapproaches. Accessed 20 Aug 2020.
Harbach RE, Gingrich JB, Pang LW. Some entomological observations on malaria transmission in a remote village in northwestern Thailand. J Am Mosq Control Assoc. 1987;3:296–301.
Suwonkerd W, Tsuda Y, Overgaard HJ, Chawprom S, Tuno N, Prajakwong S, et al. Changes in malaria vector densities over a 23 year period in Mae Hong Son Province, northern Thailand. Southeast Asian J Trop Med Public Health. 2004;35:316–24.
Alam MS, Chakma S, Khan WA, Glass GE, Mohon AN, Elahi R, et al. Diversity of anopheline species and their Plasmodium infection status in rural Bandarban, Bangladesh. Parasit Vectors. 2012;5:150.
Bhattacharyya DR, Prakash A, Sarma NP, Mohapatra PK, Singh S, Sarma DK, et al. Molecular evidence for the involvement of Anopheles nivipes (Diptera: Culicidae) in the transmission of Plasmodium falciparum in north-eastern India. Ann Trop Med Parasitol. 2010;104:331–6.
St Laurent B, Oy K, Miller B, Gasteiger EB, Lee E, Sovannaroth S, et al. Cow-baited tents are highly effective in sampling diverse Anopheles malaria vectors in Cambodia. Malar J. 2016;15:440.
Chan A, Chiang LP, Hapuarachchi HC, Tan CH, Pang SC, Lee R, et al. DNA barcoding: complementing morphological identification of mosquito species in Singapore. Parasit Vectors. 2014;7:569.
Marrelli MT, Sallum MA, Marinotti O. The second internal transcribed spacer of nuclear ribosomal DNA as a tool for Latin American anopheline taxonomy-a critical review. Mem Inst Oswaldo Cruz. 2006;101:817–32.
Avise JC. Molecular markers, natural history and evolution: natural history and evolution. Berlin: Springer Science & Business Media; 1994.
Norris DE. Genetic markers for study of the anopheline vectors of human malaria. Int J Parasitol. 2002;32:1607–15.
Perandin F, Manca N, Calderaro A, Piccolo G, Galati L, Ricci L, et al. Development of a real-time PCR assay for detection of Plasmodium falciparum, Plasmodium vivax, and Plasmodium ovale for routine clinical diagnosis. J Clin Microbiol. 2004;42:1214–9.
Stresman GH, Kamanga A, Moono P, Hamapumbu H, Mharakurwa S, Kobayashi T, et al. A method of active case detection to target reservoirs of asymptomatic malaria and gametocyte carriers in a rural area in Southern Province, Zambia. Malar J. 2010;9:265.
Das BP, Rajagopal R, Akiyama J. Pictorial key to the species of Indian Anopheline mosquitoes. Zool J Pure Appl Zool. 1990;2:131–62.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Ser. 1999;41:95–98. http://www.mbio.ncsu.edu/BioEdit/bioedit.html. Accessed 20 Aug 2020.
Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. 2004;101:11030–5.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Rambaut, A. FigTree. 2012. https://github.com/rambaut/figtree/releases. Accessed 20 Aug 2020.
Wang G, Li C, Guo X, Xing D, Dong Y, Wang Z, et al. Identifying the main mosquito species in China based on DNA barcoding. PLoS ONE. 2012;7:e47051.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA. 1973;70:3321–3.
Excofer L, Laval G, Schneider S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50.
Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11:2571–81.
Jensen JL, Bohonak AJ, Kelley ST. Isolation by distance, web service. BMC Genet. 2005;6:13.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.
Slatkin M, Hudson RR. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991;129:555–62.
Rogers AR, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992;9:552–69.
Rogers AR. Genetic evidence for a Pleistocene population explosion. Evolution. 1995;49:608–15.
SCW. The School Atlas of Cambodia. Phnom Penh, Cambodia: Save Cambodia’s Wildlife; 2006.
Guerra CA, Snow RW, Hay SI. A global assessment of closed forests, deforestation and malaria risk. Ann Trop Med Parasitol. 2006;100:189–204.
Meek SR. Vector control in some countries of Southeast Asia: comparing the vectors and the strategies. Ann Trop Med Parasitol. 1995;89:135–47.
Coosemans M, Van Bortel W. Malaria vectors in the Mekong Countries: a complex interaction between vectors, environment and human behaviour. In: proceedings of the International Conferece Hubs, Harbours and Deltas in Southeast Asia. Phnom Penh; 2006–2007. pp. 6–8.
Sanh NH, Van Dung N, Thanh NX, Trung TN, Van Co T, Cooper RD. Forest malaria in central Vietnam. Am J Trop Med Hyg. 2008;79:652–4.
Obsomer V, Dufrene M, Defourny P, Coosemans M. Anopheles species associations in Southeast Asia: indicator species and environmental influences. Parasit Vectors. 2013;6:136.
Durnez L, Mao S, Denis L, Roelants P, Sochantha T, Coosemans M. Outdoor malaria transmission in forested villages of Cambodia. Malar J. 2013;12:329.
Van Bortel W, Trung HD, le Hoi X, Van Ham N, Van Chut N, Luu ND, et al. Malaria transmission and vector behaviour in a forested malaria focus in central Vietnam and the implications for vector control. Malar J. 2010;9:373.
Kobayashi J, Nambanya S, Miyagi I, Vanachone B, Manivong K, Koubouchan T, et al. Collection of anopheline mosquitos in three villages endemic for malaria in Khammouane, Lao PDR. Southeast Asian J Trop Med Public Health. 1997;28:20.
Kobayashi J, Somboon P, Keomanila H, Inthavongsa S, Nambanya S, Inthakone S, et al. Malaria prevalence and a brief entomological survey in a village surrounded by rice fields in Khammouan province. Lao PDR Trop Med Int Health. 2000;5:17–21.
Sumarnrote A, Corbel V, Overgaard HJ, Celhay O, Marasri N, Fustec B, et al. Plasmodium infections in Anopheles mosquitoes in Ubon Ratchathani Province, Northeastern Thailand during a malaria outbreak. J Am Mosq Control Assoc. 2018;34:11–7.
Fang Y, Shi WQ, Zhang Y. Molecular phylogeny of Anopheles hyrcanus group members based on ITS2 rDNA. Parasit Vectors. 2017;10:417.
Castro JA, Picornell A, Ramon M. Mitochondrial DNA: a tool for populational genetics studies. Int Microbiol. 1998;1:327–32.
Moore WS. Inferring phylogenies from mtDNA variation: mitochondrial-gene trees versus nuclear-gene trees. Evolution. 1995;49:718–26.
Gunay F, Alten B, Simsek F, Aldemir A, Linton YM. Barcoding Turkish Culex mosquitoes to facilitate arbovirus vector incrimination studies reveals hidden diversity and new potential vectors. Acta Trop. 2015;143:112–20.
Murugan K, Vadivalagan C, Karthika P, Panneerselvam C, Paulpandi M, Subramaniam J, et al. DNA barcoding and molecular evolution of mosquito vectors of medical and veterinary importance. Parasitol Res. 2016;115:107–21.
Ashfaq M, Hebert PD, Mirza JH, Khan AM, Zafar Y, Mirza MS. Analyzing mosquito (Diptera: Culicidae) diversity in Pakistan by DNA barcoding. PLoS ONE. 2014;9:e97268.
Cywinska A, Hunter FF, Hebert PD. Identifying Canadian mosquito species through DNA barcodes. Med Vet Entomol. 2006;20:413–24.
Chen B, Harbach RE, Butlin RK. Genetic variation and population structure of the mosquito Anopheles jeyporiensis in southern China. Mol Ecol. 2004;13:3051–6.
Chen B, Pedro PM, Harbach RE, Somboon P, Walton C, Butlin RK. Mitochondrial DNA variation in the malaria vector Anopheles minimus across China, Thailand and Vietnam: evolutionary hypothesis, population structure and population history. Heredity. 2011;106:241–52.
Templeton AR. Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history. Mol Ecol. 1998;7:381–97.
Sindhania A, Das MK, Sharma G, Surendran SN, Kaushal BR, Lohani HP, et al. Molecular forms of Anopheles subpictus and Anopheles sundaicus in the Indian subcontinent. Malar J. 2020;19:417.
Yang M, Ma Y, Wu J. Mitochondrial genetic differentiation across populations of the malaria vector Anopheles lesteri from China (Diptera: Culicidae). Malar J. 2011;10:216.
Chang X, Zhong D, Lo E, Fang Q, Bonizzoni M, Wang X, et al. Landscape genetic structure and evolutionary genetics of insecticide resistance gene mutations in Anopheles sinensis. Parasit Vectors. 2016;9:228.
Toma T, Miyagi I, Crabtree MB, Miller BR. Identification of Culex vishnui subgroup (Diptera: Culicidae) mosquitoes from the Ryukyu Archipelago, Japan: development of a species-diagnostic polymerase chain reaction assay based on sequence variation in ribosomal DNA spacers. J Med Entomol. 2000;37:554–8.
Beebe NW, Whelan PI, Van den Hurk AF, Ritchie SA, Corcoran S, Cooper RD. A polymerase chain reaction-based diagnostic to identify larvae and eggs of container mosquito species from the Australian region. J Med Entomol. 2007;44:376–80.
We greatly thank the National Natural Science Foundation of China and Shanghai Science and technology innovation action plan for the study design, analysis, and interpretation of the data, as well as the writing of the manuscript and the decision to submit the manuscript for publication. We also thank the Yunnan Provincial Collaborative Innovation Center for Public Health and Disease Prevention and Control (Grant Number 2014YNPHXT03) and the China-ASEAN Public Health Cooperation Fund and Joint Control of Malaria and Dengue in the Lancang-Mekong River Area for the sample collection, data collection, and performing the experiments. We thank all the staff at the Yunnan Institute of Parasitic Diseases of China, as well as SREA SAMBO and SIEM PANG Hospital in Siem Pang County of Stung treng Province, Provincial Health Department of Stung Treng, and the National Center for Parasitology, Entomology & Malaria Control of Cambodia who participated in the sampling and data collection.
This study was supported by the National Natural Science Foundation of China (Grant Number 31601002, 81160357, 30960327, 30660160), Shanghai Science and technology innovation action plan (Grant Number 19JC1416500), Yunnan Provincial Collaborative Innovation Center for Public Health and Disease Prevention and Control (Grant Number 2014YNPHXT03) and the China-ASEAN Public Health Cooperation Fund and Joint Control of Malaria and Dengue in Lancang-Mekong River Area.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Table S1. Full list of 53 Anopheles nivipes specimens which were classified by both molecular and morphological identifications, with morphology species ID, molecular species ID based on ITS2, molecular species ID based on COII, geographical location, latitude, and longitude.
: Table S2. Genus/Species compositions of mosquitoes trapped by CDC lamp in SIEM PANG County, Stung treng Province. Adult mosquitoes were collected by using overnight trapping with the battery-operated CDC light traps. ** indicated the genus or species composition of mosquitoes which were described in percentage. (%)** indicated the genus or species composition of mosquitoes which were described in numbers.
: Table S3. COII sequences of Anopheles nivipes were downloaded from the NCBI. * indicated the longitude and latitude coordinates to the geographical center of a certain province, due to the samples being initially collected from various sampling sites in a certain province. An. niv, Anopheles nivipes.
: Figure S1. Map of the populations from different geographical regions. Populations of mitochondrial COII sequences: KH, Siem Pang County (Stung treng, Cambodia); Tri, Tripura (India); Ch, Cheema (Punjab, India); Ba, Bathinda (Punjab, India); Nag, Nagaland (India); Ass, Assam (India); MM, Myanmar; TH, Thailand; Miz, Mizoram (India); Meg, Meghalaya (India). The map was prepared by using LocaSpace Viewer.
: Figure S2. Distribution heatmap of haplotype based on COII. KH-St-Sp, Siem Pang County (Stung treng, Cambodia); MM, Myanmar; TH, Thailand; IN-Pun-Ba, Bathinda (Punjab, India); IN-Pun-Ch, Cheema (Punjab, India); IN-Tri, Tripura (India); IN-Meg, Meghalaya (India); IN-Ass, Assam (India); IN-Nag, Nagaland (India); IN-Miz, Mizoram (India). The numbers of haplotypes are shown on the right side of the figure. The color scale ranges from blue to red, showing a range from minimum number (0) to maximum numbers (19) for each haplotype.
: Figure S3. Neighbor-joining phylogenetic tree of An. nivipes haplotypes based on COII sequences from GenBank and original data in this study. Bootstrap values (1000 replicates) of Neighbor-Joining analyses are shown above/below the main lineages. Lineage designation is indicated on the right. Bars represent 2.0 substitutions per site based on COII. Different colors indicated different population groups of An. nivipes.
: Table S4. Genetic differentiation and Gene flow among the Geographic Groups based on COII. The pairwise FST values and Nm values based on the COII are shown below and above the diagonal, respectively. Characters in bold indicate the significance (P < 0.05). inf, infinite. KH_St_Sp, Siem Pang County (Stung treng, Cambodia); IN_Tri, Tripura (India); IN_Pun_Ch, Cheema (Punjab, India); IN_Pun_Ba, Bathinda (Punjab, India); IN_Nag, Nagaland (India); IN_Ass, Assam (India); MM, Myanmar; TH, Thailand; IN_Miz, Mizoram (India); IN-Meg, Meghalaya (India).
: Figure S4. Cluster analysis based on COII sequences in Anopheles nivipes populations. UPGMA dendrogram based on Nei ‘s unbiased genetic distance between different populations of An. nivipes. Bars represent 0.7 substitutions per site based on COII. KS, Siem Pang County (Stung treng, Cambodia); IN_Tri, Tripura (India); IN_Pun_Ch, Cheema (Punjab, India); IN_Pun_Ba, Bathinda (Punjab, India); IN_Nag, Nagaland (India); IN_Ass, Assam (India); MM, Myanmar; TH, Thailand; IN_Miz, Mizoram (India); IN_Meg, Meghalaya (India).
: Figure S5. Mismatch distribution graphs for Siem Pang population. The X and Y-axis show the number of pairwise differences and the frequency of the pairwise comparisons, respectively. The observed frequencies are represented by a dotted line. The frequency expected under the hypothesis of the constant population model is depicted by a solid line. (a) all populations-COII; (b) Siem Pang population-COII.
About this article
Cite this article
Zhang, Y., Zhang, C., Yang, R. et al. Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border. Malar J 21, 91 (2022). https://doi.org/10.1186/s12936-022-04121-w
- Molecular phylogeny