Review of genetic diversity in malaria vectors (Culicidae: Anophelinae)

Integrating wing morphometrics and mitochondrial DNA analysis to assess the filaria vector Mansonia uniformis (Diptera: Culicidae) populations in Thailand

Mansonia uniformis (Diptera: Culicidae) is recognized as a vector of Brugia malayi and has been reported to transmit Wuchereria bancrofti, both causing lymphatic filariasis in humans. This study employed geometric morphometrics (GM) to investigate wing shape variation and analyzed genetic diversity through cytochrome c oxidase subunit 1 (COI) gene analyses in Ma. uniformis populations across Thailand. Wing GM analyses indicated significant differences in wing shape based on Mahalanobis distances among nearly all population pairs (p < 0.05), with no significant correlation between wing shape and geographic distance (r = 0.210, p > 0.05). Genetic analyses identified 63 haplotypes and 49 polymorphic sites, with the overall population exhibiting a nucleotide diversity of 0.006 (± 0.001) and a haplotype diversity of 0.912 (± 0.017). Deviations from neutrality, as indicated by Tajima's D and Fu's FS tests for the overall Ma. uniformis populations in Thailand, were statistically significant and negative, suggesting population expansion (both p < 0.05). Analysis of molecular variance revealed no significant genetic structure when all populations were categorized based on collection sites and geographic regions. However, significant differences in FST values were observed between some populations. These findings enhance our understanding of the geographical and genetic factors influencing Ma. uniformis populations, which are crucial for developing effective control strategies in Thailand.

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.

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.

Effect of land cover and landscape fragmentation on anopheline mosquito abundance and diversity in an important Colombian malaria endemic region

Landscape structure influences the distribution and abundance of anopheline mosquitoes and has an indirect impact on malaria transmission. This work aimed to determine the effect of land cover and landscape fragmentation on anopheline mosquito abundance and diversity in an important Colombian malaria endemic area, the Bajo Cauca region. Diversity indices were calculated for Anopheles mosquitoes collected in various localities of the region. Land cover types were characterized using orthorectified aerial photographs to estimate landscape metrics. The relationship between landscape fragmentation and species diversity was evaluated by regression analysis. The correlation between species abundance and land cover types was determined using canonical correspondence analyses. Results showed a statistically significant tendency for a lower diversity of the Anopheles community in landscapes with higher patch number, patch density and effective mesh size. For most species, there was evidence of a significant relationship between species abundance and land covers modified by anthropic activities which generate forest loss. These results indicate that activities that modify the landscape structure and land cover composition generate changes that affect the spatial distribution and composition of epidemiologically-important Anopheles species, which may impact malaria distribution in a region. This information is useful to guide control interventions that promote unfavorable landscapes for malaria vector propagation.

Molecular characterization and phylogenetic analysis of anopheline (Anophelinae: Culicidae) mosquitoes of the Oriental and Afrotropical Zoogeographic zones in Saudi Arabia

The Kingdom of Saudi Arabia (KSA) has a diverse fauna due to its peculiar position bordering the Afrotropical, Oriental and Palaearctic zoogeographic zones. The present study reports the phylogenetics of five mosquito species belonging to five series of Anopheles (Cellia) . We collected mosquito larvae from eastern, western and southwestern regions of KSA. The sampled mosquitoes were morphologically identified using the pictorial keys of mosquitoes and characterized by using single and multi-locus analysis of -internal transcribed spacer 2 (ITS2) region and cytochrome oxidase c subunit I (COI). Based on the morphological and molecular data, five species were recognized, like An. stephensi (Neocellia) (Oriental), An. arabiensis (Pyretophorus) (Afrotropical), An. dthali (Myzomyia) (Oriental and Palaearctic), An. cinereus (Paramyzomyia) and An. rhodesiensis rupicola (Neomyzomyia) (Oriental and Palaearctic). The phylogenetic analysis showed that An. stephensi is a monophyletic species with different ecotypes found in different geographic regions. Comprehensive phylogenetics and population genetics studies are crucial for a better understanding of the role of these five mosquito species in malarial transmission across various zoogeographic zones of different ecological and demographic characteristics.

Diversity and Molecular Characterization of Mosquitoes (Diptera: Culicidae) in Selected Ecological Regions in Kenya

Mosquitoes play a predominant role as leading agents in the spread of vector-borne diseases and the consequent mortality in humans. Despite reports on increase of new and recurrent mosquito borne-disease outbreaks such as chikungunya, dengue fever and Rift Valley fever in Kenya, little is known about the genetic characteristics and diversity of the vector species that have been incriminated in transmission of disease pathogens. In this study,  mosquito species were collected from Kisumu city, Kilifi town and Nairobi city and we determined their genetic diversity and phylogenetic relationships. PCR was used to amplify the partial cytochrome oxidase subunit 1 (CO1) gene of mosquito samples. Molecular-genetic and phylogenetic analysis of the partial cytochrome oxidase subunit 1 (CO1) gene were employed to identify their relationship with known mosquito species. Fourteen (14) haplotypes belonging to genus Aedes, nine (9) haplotypes belonging to genus Anopheles and twelve (12) haplotypes belonging to genus Culex were identified in this study. Findings from this study revealed a potentially new haplotype belonging to Anopheles genus and reported the first molecular characterization of Aedes cumminsii in Kenya. Sequence results revealed variation in mosquito species from Kilifi, Kisumu and Nairobi. Since vector competence varies greatly across species as well as species-complexes and is strongly associated with specific behavioural adaptations, proper species identification is important for vector control programs.

Microevolution of medically important mosquitoes - A review

This review intends to discuss central issues regarding the microevolution of mosquito (Culicidae) vectors of several pathogens and how this process impacts vector biology, disease transmission, and vector control attempts. On the microevolutionary context, it comparatively discusses the current knowledge on the population genetics of representatives of the genera Aedes, Anopheles and Culex, and comments on insecticide resistance of culicids. It also discusses other biological aspects of culicids that are not usually addressed in microevolutionary studies, such as vectorial competence, endosymbiosis, and wing morphology. One conclusion is that mosquitoes are highly genetically variable, adaptable, fast evolving, and have versatile vectorial competence. Unveiling microevolutionary patterns is fundamental for the design and maintenance of all control programs. Sampling methods for assessing microevolution must be standardized and must follow meaningful guidelines, such as those of "landscape genetics". A good understanding of microevolution requires more than a collection of case studies on population genetics and resistance. Future research could deal not only with the microevolution sensu stricto, but also with evolutionarily meaningful issues, such as inheritable characters, epigenetics, physiological cost-free plasticity, vector immunity, symbiosis, pathogen-mosquito co-evolution and environmental variables. A genotyping panel for seeking adaptive phenotypes as part of the standardization of population genetics methods is proposed. The investigative paradigm should not only be retrospective but also prospective, despite the unpredictability of evolution. If we integrate all suggestions to tackle mosquito evolution, a global revolution to counter vector-borne diseases can be provoked.

Molecular taxonomy and evolutionary relationships in the Oswaldoi-Konderi complex (Anophelinae: Anopheles: Nyssorhynchus) from the Brazilian Amazon region

Recent studies have shown that Anopheles oswaldoi sensu lato comprises a cryptic species complex in South America. Anopheles konderi, which was previously raised to synonymy with An. oswaldoi, has also been suggested to form a species complex. An. oswaldoi has been incriminated as a malaria vector in some areas of the Brazilian Amazon, Colombia, Peru and Venezuela, but was not recognized as a vector in the remaining regions in its geographic distribution. The role of An. konderi as a malaria vector is unknown or has been misattributed to An. oswaldoi. The focus of this study was molecular identification to infer the evolutionary relationships and preliminarily delimit the geographic distribution of the members of these complexes in the Brazilian Amazon region. The specimens were sampled from 18 localities belonging to five states in the Brazilian Amazon and sequenced for two molecular markers: the DNA barcode region (COI gene of mitochondrial DNA) and Internal Transcribed Spacer 2 (ITS2 ribosomal DNA). COI (83 sequences) and ITS2 (27 sequences) datasets generated 43 and 10 haplotypes, respectively. Haplotype networks and phylogenetic analyses generated with the barcode region (COI gene) recovered five groups corresponding to An. oswaldoi s.s., An. oswaldoi B, An. oswaldoi A, An. konderi and An. sp. nr. konderi; all pairwise genetic distances were greater than 3%. The group represented by An. oswaldoi A exhibited three strongly supported lineages. The molecular dating indicated that the diversification process in these complexes started approximately 2.8 Mya, in the Pliocene. These findings confirm five very closely related species and present new records for these species in the Brazilian Amazon region. The paraphyly observed for the An. oswaldoi complex suggests that An. oswaldoi and An. konderi complexes may comprise a unique species complex named Oswaldoi-Konderi. Anopheles oswaldoi B may be a potential malaria vector in the extreme north of the Brazilian Amazon, whereas evidence of sympatry for the remaining species in other parts of the Brazilian Amazon (Acre, Amazonas, Pará and Rondônia) precluded identification of probable vectors in those areas.

Cryptic diversity in an Atlantic Forest malaria vector from the mountains of South-East Brazil

BACKGROUND

Anopheles (Kerteszia) cruzii is the primary vector of human and simian malarias in Brazilian regions covered by the Atlantic Rainforest. Previous studies found that An. cruzii presents high levels of behavioural, chromosomal and molecular polymorphisms, which led to the hypothesis that it may be a complex of cryptic species. Here, An. cruzii specimens were collected in five sites in South-East Brazil located at different altitudes on the inner and coastal slopes of two mountain ranges covered by Atlantic Rainforest, known as Serra do Mar and Serra da Mantiqueria. Partial sequences for two genes (Clock and cpr) were generated and compared with previously published sequences from Florianópolis (southern Brazil). Genetic diversity was analysed with estimates of population structure (F _ST ) and haplotype phylogenetic trees in order to understand how many species of the complex may occur in this biome and how populations across the species distribution are related.

RESULTS

The sequences from specimens collected at sites located on the lower coastal slopes of Serra do Mar (Guapimirim, Tinguá and Sana) clustered together in the phylogenetic analysis, while the major haplotypes from sites located on higher altitude and at the continental side of the same mountains (Bocaina) clustered with those from Serra da Mantiqueira (Itatiaia), an inner mountain range. These two An. cruzii lineages showed statistically significant genetic differentiation and fixed characters, and have high F _ST values typical of between species comparisons. Finally, in Bocaina, where the two lineages occur in sympatry, we found deviations from Hardy-Weinberg equilibrium due to a deficit of heterozygotes, indicating partial reproductive isolation. These results strongly suggest that at least two distinct lineages of An. cruzii (provisorily named "Group 1" and "Group 2") occur in the mountains of South-East Brazil.

CONCLUSIONS

At least two genetically distinct An. cruzii lineages occur in the Atlantic Forest covered mountains of South-East Brazil. The co-occurrence of distinct lineages of An. cruzii (possibly incipient species) in those mountains is an interesting biological phenomenon and may have important implications for malaria prevalence, Plasmodium transmission dynamics and control.

Genetic structure and phenotypic variation of Anopheles darlingi in northwest Colombia

This study evaluated the influence of environmental heterogeneity on Anopheles darlingi genetic and morphometric traits at a microgeographic level. Specimens of An. darlingi collected from multiple municipalities in the Colombian malaria endemic region Urabá-Bajo Cauca and Alto Sinú (UCS) were analyzed using 13 microsatellite loci. Spatial genetic structure, population variation and wing geometric morphometric analyses were performed. Microsatellite results showed low genetic differentiation and high gene flow among populations; four highly admixed subpopulations were detected with no particular association to the municipalities. Wing geometric morphometrics analysis showed a subtle but significant difference in wing shape for El Bagre vs. Mutatá populations, possibly influenced by geographical distance. Discrimination among populations in the morphospace showed a slight separation of the Tierralta population. There was no significant correlation between the genetic and geographic or genetic and environmental distances. We hypothesize that environmental heterogeneity in the UCS region does not reach a threshold to affect population structure of An. darlingi. Another possibility is that microsatellites are not sensitive enough to detect existing structure. It remains to be determined which local factors govern phenotypic variation among these populations and how, or whether these may affect mosquito biology and transmission capacity.

Phylogenetic Relationships Among Malaria Vectors and Closely Related Species in Thailand Using Multilocus DNA Sequences

The evolutionary and taxonomic status is important for understanding speciation events and phylogenetic relationships between closely related vector and nonvector species. This information is useful for targeting important disease vector species groups for the development of novel genetic-based vector and pathogen control methods. In this study, different phylogenetic analyses were performed to reconstruct phylogenetic trees for the primary malaria vectors in Thailand based on sequence information of 4 DNA fragments from the nuclear and mitochondrial regions. The primary Anopheles species in the subgenus Cellia involved in malaria transmission in Thailand separate clearly into 3 distinct clades: the Leucosphyrus group, Minimus subgroup, and Maculatus group. The phylogenetic trees based on different reconstructed algorithms and different gene regions provided congruent phylogenetic status of the mosquito species studied. The phylogenetic relationships of malaria vector species examined followed similar patterns based on morphological characters. An estimate of the divergence time among the Anopheles species infers that they were present during the Eocene and Miocene periods (>41 million years ago). Congruent phylogenetic analysis of malaria vectors is presented with different algorithms and gene regions. The nuclear TOLL6 fragment appears useful for molecular phylogenetic, species DNA barcode, and Anopheles population genetic analyses.

Advancing vector biology research: a community survey for future directions, research applications and infrastructure requirements

Vector-borne pathogens impact public health, animal production, and animal welfare. Research on arthropod vectors such as mosquitoes, ticks, sandflies, and midges which transmit pathogens to humans and economically important animals is crucial for development of new control measures that target transmission by the vector. While insecticides are an important part of this arsenal, appearance of resistance mechanisms is increasingly common. Novel tools for genetic manipulation of vectors, use of Wolbachia endosymbiotic bacteria, and other biological control mechanisms to prevent pathogen transmission have led to promising new intervention strategies, adding to strong interest in vector biology and genetics as well as vector-pathogen interactions. Vector research is therefore at a crucial juncture, and strategic decisions on future research directions and research infrastructure investment should be informed by the research community. A survey initiated by the European Horizon 2020 INFRAVEC-2 consortium set out to canvass priorities in the vector biology research community and to determine key activities that are needed for researchers to efficiently study vectors, vector-pathogen interactions, as well as access the structures and services that allow such activities to be carried out. We summarize the most important findings of the survey which in particular reflect the priorities of researchers in European countries, and which will be of use to stakeholders that include researchers, government, and research organizations.

Genome-Wide Divergence in the West-African Malaria Vector Anopheles melas

Anopheles melas is a member of the recently diverged An. gambiae species complex, a model for speciation studies, and is a locally important malaria vector along the West-African coast where it breeds in brackish water. A recent population genetic study of An. melas revealed species-level genetic differentiation between three population clusters. An. melas West extends from The Gambia to the village of Tiko, Cameroon. The other mainland cluster, An. melas South, extends from the southern Cameroonian village of Ipono to Angola. Bioko Island, Equatorial Guinea An. melas populations are genetically isolated from mainland populations. To examine how genetic differentiation between these An. melas forms is distributed across their genomes, we conducted a genome-wide analysis of genetic differentiation and selection using whole genome sequencing data of pooled individuals (Pool-seq) from a representative population of each cluster. The An. melas forms exhibit high levels of genetic differentiation throughout their genomes, including the presence of numerous fixed differences between clusters. Although the level of divergence between the clusters is on a par with that of other species within the An. gambiae complex, patterns of genome-wide divergence and diversity do not provide evidence for the presence of pre- and/or postmating isolating mechanisms in the form of speciation islands. These results are consistent with an allopatric divergence process with little or no introgression.

Evidence of new species for malaria vector Anopheles nuneztovari sensu lato in the Brazilian Amazon region

BACKGROUND

Anopheles nuneztovari sensu lato comprises cryptic species in northern South America, and the Brazilian populations encompass distinct genetic lineages within the Brazilian Amazon region. This study investigated, based on two molecular markers, whether these lineages might actually deserve species status.

METHODS

Specimens were collected in five localities of the Brazilian Amazon, including Manaus, Careiro Castanho and Autazes, in the State of Amazonas; Tucuruí, in the State of Pará; and Abacate da Pedreira, in the State of Amapá, and analysed for the COI gene (Barcode region) and 12 microsatellite loci. Phylogenetic analyses were performed using the maximum likelihood (ML) approach. Intra and inter samples genetic diversity were estimated using population genetics analyses, and the genetic groups were identified by means of the ML, Bayesian and factorial correspondence analyses and the Bayesian analysis of population structure.

RESULTS

The Barcode region dataset (N = 103) generated 27 haplotypes. The haplotype network suggested three lineages. The ML tree retrieved five monophyletic groups. Group I clustered all specimens from Manaus and Careiro Castanho, the majority of Autazes and a few from Abacate da Pedreira. Group II clustered most of the specimens from Abacate da Pedreira and a few from Autazes and Tucuruí. Group III clustered only specimens from Tucuruí (lineage III), strongly supported (97 %). Groups IV and V clustered specimens of A. nuneztovari s.s. and A. dunhami, strongly (98 %) and weakly (70 %) supported, respectively. In the second phylogenetic analysis, the sequences from GenBank, identified as A. goeldii, clustered to groups I and II, but not to group III. Genetic distances (Kimura-2 parameters) among the groups ranged from 1.60 % (between I and II) to 2.32 % (between I and III). Microsatellite data revealed very high intra-population genetic variability. Genetic distances showed the highest and significant values (P = 0.005) between Tucuruí and all the other samples, and between Abacate da Pedreira and all the other samples. Genetic distances, Bayesian (Structure and BAPS) analyses and FCA suggested three distinct biological groups, supporting the barcode region results.

CONCLUSIONS

The two markers revealed three genetic lineages for A. nuneztovari s.l. in the Brazilian Amazon region. Lineages I and II may represent genetically distinct groups or species within A. goeldii. Lineage III may represent a new species, distinct from the A. goeldii group, and may be the most ancestral in the Brazilian Amazon. They may have differences in Plasmodium susceptibility and should therefore be investigated further.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Mitochondrial DNA-based genetic diversity of Anopheles nivipes in North East India

Malaria is a major public health problem in north-east India mainly transmitted by Anopheles baimai and An. minimus while Anopheles nivipes plays an important supportive role. The genetic diversity of An. nivipes in north-east India was investigated by employing two mitochondrial DNA markers namely NADH dehydrogenase 5 (ND5) and cytochrome oxidase sub unit II (COII). High genetic diversity in An. nivipes was observed by the detection of 16 haplotypes among 30 sequences of ND5 gene and 29 haplotypes among 29 COII sequences. Anopheles nivipes of north-east India was significantly differentiated genetically with that of neighboring South-east Asian An. nivipes as revealed by pairwise F_ST values of 0.127 (p < 0.01) and 0.044 (p < 0.001) for ND5 and COII genes, respectively, suggesting geographical barriers to gene flow in this species between the two geographical areas resulting in significant population structuring.

Molecular taxonomy provides new insights into anopheles species of the neotropical arribalzagia series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated An. apicimacula MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.

Genetic diversity and population structure in the Leishmania guyanensis vector Lutzomyia anduzei (Diptera, Psychodidae) from the Brazilian Amazon

Lutzomyia (Nyssomyia) anduzei has been recognized as a secondary vector of Leishmania guyanensis in the Brazilian Amazon region. Since L. anduzei is anthropophilic, co-exists in areas of high leishmaniasis transmission and has been found infected with L. guyanensis, the understanding of the vector population structure and of the process responsible for it is paramount to the vector management and control efforts. In this study we analyzed 74 and 67 sequences of the COI and Cytb loci, respectively, from mitochondrial DNA, aiming to estimate the intra-population genetic variability and population structure in six L. anduzei samples from the Brazilian Amazon region. For COI, we found 58 haplotypes, low to high (FST=0.0310-0.4128) and significant (P=0.0033) genetic structure, and reduced gene flow among populations. The haplotype network yielded many reticulations that likely resulted from hypervariability in the locus. For Cytb, we observed 27 haplotypes, low to moderate (FST=0.0077-0.1954) and nonsignificant (P>0.05) genetic structure for the majority of comparisons and extensive gene flow among populations, in line with the haplotypes network data. AMOVA analysis indicated that most of the variation occurred within populations (83.41%, 90.94%); nevertheless, there were significant differences (ΦST=0.0906-0.1659; P=0.00098; P=0.00000) among them for both loci. The Mantel test showed that the genetic structure is not associated to an isolation-by-distance (IBD) model in either of both loci. These data suggest that L. anduzei is genetically very diverse. The genetic structure lacking IBD may be due to adaptation to local habitats and the low dispersal capacity of the sandflies, and both could lead to population fragmentation and geographic isolation. These findings have important implications for epidemiology, surveillance and vector control and may be a first step in understanding the evolutionary history of this species.

Altitudinal population structure and microevolution of the malaria vector Anopheles cruzii (Diptera: Culicidae)

BACKGROUND

In Brazil, the autochthonous transmission of extra-Amazonian malaria occurs mainly in areas of the southeastern coastal Atlantic Forest, where Anopheles cruzii is the primary vector. In these locations, the population density of the mosquito varies with altitude (5-263 m above sea level), prompting us to hypothesise that gene flow is also unevenly distributed. Describing the micro-geographical and temporal biological variability of this species may be a key to understanding the dispersion of malaria in the region. We explored the homogeneity of the An. cruzii population across its altitudinal range of distribution using wing shape and mtDNA gene analysis. We also assessed the stability of wing geometry over time.

METHODS

Larvae were sampled from lowland (5-20 m) and hilltop (81-263 m) areas in a primary Atlantic Forest region, in the municipality of Cananéia (State of São Paulo, Brazil). The right wings of males and females were analysed by standard geometric morphometrics. Eighteen landmarks were digitised for each individual and a discriminant analysis was used to compare samples from the hilltop and lowland. A 400-bp DNA fragment of the mitochondrial cytochrome oxidase gene subunit I (CO-I) was PCR-amplified and sequenced.

RESULTS

Wing shapes were distinct between lowland and hilltop population samples. Results of cross-validated tests based on Mahalanobis distances showed that the individuals from both micro-environments were correctly reclassified in a range of 54-96%. The wings of hilltop individuals were larger. The CO-I gene was highly polymorphic (haplotypic diversity = 0.98) and altitudinally structured (Фst = 0.085 and Jaccard = 0.033). We found 60 different haplotypes but only two were shared by the lowland and hilltop populations. Wing shape changed over the brief study period (2009-2013).

CONCLUSIONS

Wing geometry and CO-I gene analysis indicated that An. cruzii is vertically structured. Wing shape varied rapidly, but altitude structure was maintained. Future investigations should identify the biotic/abiotic causes of these patterns and their implications in the local epidemiology of malaria.

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.

Molecular characterization, biological forms and sporozoite rate of Anopheles stephensi in southern Iran

OBJECTIVE

To identify the biological forms, sporozoite rate and molecular characterization of the Anopheles stephensi (An. stephensi) in Hormozgan and Sistan-Baluchistan provinces, the most important malarious areas in Iran.

METHODS

Wild live An. stephensi samples were collected from different malarious areas in southern Iran. The biological forms were identified based on number of egg-ridges. Molecular characterization of biological forms was verified by analysis of the mitochondrial cytochrome oxidase subunit I and II (mtDNA-COI/COII). The Plasmodium infection was examined in the wild female specimens by species-specific nested-PCR method.

RESULTS

Results showed that all three biological forms including mysorensis, intermediate and type are present in the study areas. Molecular investigations revealed no genetic variation between mtDNA COI/COII sequences of the biological forms and no Plasmodium parasites was detected in the collected mosquito samples.

CONCLUSIONS

Presence of three biological forms with identical sequences showed that the known biological forms belong to a single taxon and the various vectorial capacities reported for these forms are more likely corresponded to other epidemiological factors than to the morphotype of the populations. Lack of malaria parasite infection in An. stephensi, the most important vector of malaria, may be partly due to the success and achievement of ongoing active malaria control program in the region.

Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation

Anopheles albimanus is a major malaria mosquito vector in Colombia. In the present study, wing variability (size and shape) in An. albimanus populations from Colombian Maracaibo and Chocó bio-geographical eco-regions and the relationship of these phenotypic traits with environmental factors were evaluated. Microsatellite and morphometric data facilitated a comparison of the genetic and phenetic structure of this species. Wing size was influenced by elevation and relative humidity, whereas wing shape was affected by these two variables and also by rainfall, latitude, temperature and eco-region. Significant differences in mean shape between populations and eco-regions were detected, but they were smaller than those at the intra-population level. Correct assignment based on wing shape was low at the population level (<58%) and only slightly higher (>70%) at the eco-regional level, supporting the low population structure inferred from microsatellite data. Wing size was similar among populations with no significant differences between eco-regions. Population relationships in the genetic tree did not agree with those from the morphometric data; however, both datasets consistently reinforced a panmictic population of An. albimanus. Overall, site-specific population differentiation is not strongly supported by wing traits or genotypic data. We hypothesize that the metapopulation structure of An. albimanus throughout these Colombian eco-regions is favoring plasticity in wing traits, a relevant characteristic of species living under variable environmental conditions and colonizing new habitats.

Functional constraints and evolutionary dynamics of the repeats in the rDNA internal transcribed spacer 2 of members of the Anopheles barbirostris group

BACKGROUND

The Anopheles barbirostris group is widely distributed in Southeast Asia. Although seven species have been formally described, a molecular analysis of the rDNA ITS2 and the mitochondrial cytochrome oxidase I gene suggests that the group includes species that are morphologically very similar or identical.We have previously shown that species in the Anopheles barbirostris Subgroup have an exceptionally large ITS2 (>1.5 kb), greater than in any other Anopheline group. However, the molecular processes responsible for generating such a large ITS2 have not previously been explored.

METHODS

To determine the processes by which this large ITS2 is generated, we examined the sequence and secondary structure of the ITS2 of 51 specimens from five species of the Anopheles barbirostris Subgroup. These include the anthropophilic species An. campestris and three morphospecies of the Barbirostris Complex: An. vanderwulpi, An. barbirostris I and III, together with a previously undescribed member of this group (Clade IV).

RESULTS AND CONCLUSIONS

All the specimens were found to have an ITS2 greater than 1.5 kb in length. The possibility that the spacer sequences amplified were pseudogenes was examined and discarded. The large size of ITS2 in the species studied is due to the presence of internal repeats of approximately 110 bp in length, confined to the central region of the spacer. Repeats varied markedly between the species examined, with respect to their organization, number and sequence similarity. The nucleotide diversity increased in direct relation to size variation and the presence of non-repeated elements.A secondary structure analysis showed that the repeats form hairpin structures with a wide range of free energy values. These hairpin structures are known to facilitate the subsequent processing of mature rRNA. An analysis of the repeats from the different species suggests they originate from a common ancestor, with the repeats appearing before speciation of the Barbirostris Group.

Needs for monitoring mosquito transmission of malaria in a pre-elimination world

As global efforts to eliminate malaria intensify, accurate information on vector populations and transmission dynamics is critical for directing control efforts, developing new control tools, and predicting the effects of these interventions under various conditions. Currently available sampling tools for mosquito population monitoring suffer from well-recognized limitations. As reported in this workshop summary, a recent gathering of medical entomologists, modelers, and malaria experts reviewed these issues and agreed that efforts are needed to improve methods to monitor key transmission parameters. Identified needs include standardized methods for sampling of both mosquito adults and larvae, improved tools for mosquito species identification and age-grading, and a better means for determining the entomological inoculation rate.

Novel genetic diversity within Anopheles punctimacula s.l.: phylogenetic discrepancy between the Barcode cytochrome c oxidase I (COI) gene and the rDNA second internal transcribed spacer (ITS2)

Anopheles punctimacula s.l. is a regional malaria vector in parts of Central America, but its role in transmission is controversial due to its unresolved taxonomic status. Two cryptic species, An. malefactor and An. calderoni, have been previously confused with this taxon, and evidence for further genetic differentiation has been proposed. In the present study we collected and morphologically identified adult female mosquitoes of An. punctimacula s.l. from 10 localities across Panama and one in Costa Rica. DNA sequences from three molecular regions, the three prime end of the mitochondrial cytochrome c oxidase I gene (3' COI), the Barcode region in the five prime end of the COI (5' COI), and the rDNA second internal transcribed spacer (ITS2) were used to test the hypothesis of new molecular lineages within An. punctimacula s.l. Phylogenetic analyses using the 3' COI depicted six highly supported molecular lineages (A-F), none of which was An. malefactor. In contrast, phylogenetic inference with the 5' COI demonstrated paraphyly. Tree topologies based on the combined COI regions and ITS2 sequence data supported the same six lineages as the 3' COI alone. As a whole this evidence suggests that An. punctimacula s.l. comprises two geographically isolated lineages, but it is not clear whether these are true species. The phylogenetic structure of the An. punctimacula cluster as well as that of other unknown lineages (C type I vs C type II; D vs E) appears to be driven by geographic partition, because members of these assemblages did not overlap spatially. We report An. malefactor for the first time in Costa Rica, but our data do not support the presence of An. calderoni in Panama.

Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

BACKGROUND

Anopheles (Kerteszia) cruzii (Diptera: Culicidae) is a primary vector of human and simian malaria parasites in southern and southeastern Brazil. Earlier studies using chromosome inversions, isoenzymes and a number of molecular markers have suggested that An. cruzii is a species complex.

RESULTS

In this study, a multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was carried out to investigate in more detail the genetic differentiation between the An. cruzii populations from Florianópolis-Santa Catarina (southern Brazil) and Itatiaia-Rio de Janeiro States (southeastern Brazil). The analyses were performed first comparing Florianópolis and Itatiaia, and then comparing the two putative sympatric incipient species from Itatiaia (Itatiaia A and Itatiaia B). The analysis revealed high FST values between Florianópolis and Itatiaia (considering Itatiaia A and B together) and also between the sympatric Itatiaia A and Itatiaia B, irrespective of their function. Also, using the IM program, no strong indication of migration was found between Florianópolis and Itatiaia (considering Itatiaia A and B together) using all loci together, but between Itatiaia A and Itatiaia B, the results show evidence of migration only in the direction of Itatiaia B.

CONCLUSIONS

The results of the multilocus analysis indicate that Florianópolis and Itatiaia represent different species of the An. cruzii complex that diverged around 0.6 Mya, and also that the Itatiaia sample is composed of two sympatric incipient species A and B, which diverged around 0.2 Mya. Asymmetric introgression was found between the latter two species despite strong divergence in some loci.

Genetic analysis of invasive Aedes albopictus populations in Los Angeles County, California and its potential public health impact

The Asian tiger mosquito, Aedes albopictus, is an anthropophilic aggressive daytime-biting nuisance and an efficient vector of certain arboviruses and filarial nematodes. Over the last 30 years, this species has spread rapidly through human travel and commerce from its native tropical forests of Asia to every continent except Antarctica. In 2011, a population of Asian tiger mosquito (Aedes albopictus) was discovered in Los Angeles (LA) County, California. To determine the probable origin of this invasive species, the genetic structure of the population was compared against 11 populations from the United States and abroad, as well as preserved specimens from a 2001 introduction into California using the mitochondrial cytochrome c oxidase 1 (CO1) gene. A total of 66 haplotypes were detected among samples and were divided into three main groups. Aedes albopictus collected in 2001 and 2011 from LA County were genetically related and similar to those from Asia but distinct from those collected in the eastern and southeastern United States. In view of the high genetic similarities between the 2001 and 2011 LA samples, it is possible that the 2011 population represents in part the descendants of the 2001 introduction. There remains an imperative need for improved surveillance and control strategies for this species.

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

BACKGROUND

The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated.

METHODS

The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches.

RESULTS

Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis.

CONCLUSION

Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.

An epidemiological overview of malaria in Bangladesh

Bangladesh is one of the four major malaria-endemic countries in South-East Asia having approximately 34% of its population at risk of malaria. This paper aims at providing an overview of the malaria situation in this country. Relevant information was retrieved from published articles and reports in PubMed and Google Scholar. Malaria in Bangladesh is concentrated in 13 districts with a prevalence ranging between 3.1% and 36%, and is mostly caused by Plasmodium falciparum. Geographical conditions pose a potential risk for Plasmodium knowlesi malaria. Resistance to a number of drugs previously recommended for treatment has been reported. Low socio-economic status, poor schooling and close proximity to water bodies and forest areas comprise important risk factors. Despite the significant steps in Long Lasting Insecticide Net (LLIN)/Insecticide Treated Net (ITN) coverage in Bangladesh, there are still many challenges including the extension of malaria support to the remote areas of Bangladesh, where malaria prevalence is higher, and further improvements in the field of referral system and treatment.

Lutzomyia umbratilis, the main vector of Leishmania guyanensis, represents a novel species complex?

BACKGROUND

Lutzomyia umbratilis is an important Leishmania guyanensis vector in South America. Previous studies have suggested differences in the vector competence between L. umbratilis populations situated on opposite banks of the Amazonas and Negro Rivers in the central Amazonian Brazil region, likely indicating a species complex. However, few studies have been performed on these populations and the taxonomic status of L. umbratilis remains unclear.

METHODOLOGY/PRINCIPAL FINDINGS

Phylogeographic structure was estimated for six L. umbratilis samples from the central Amazonian region in Brazil by analyzing mtDNA using 1181 bp of the COI gene to assess whether the populations on opposite banks of these rivers consist of incipient or distinct species. The genetic diversity was fairly high and the results revealed two distinct clades ( = lineages) with 1% sequence divergence. Clade I consisted of four samples from the left bank of the Amazonas and Negro Rivers, whereas clade II comprised two samples from the right bank of Negro River. No haplotypes were shared between samples of two clades. Samples within clades exhibited low to moderate genetic differentiation (F(ST) = -0.0390-0.1841), whereas samples between clades exhibited very high differentiation (F(ST) = 0.7100-0.8497) and fixed differences. These lineages have diverged approximately 0.22 Mya in the middle Pleistocene. Demographic expansion was detected for the lineages I and II approximately 30,448 and 15,859 years ago, respectively, in the late Pleistocene.

CONCLUSIONS/SIGNIFICANCE

The two genetic lineages may represent an advanced speciation stage suggestive of incipient or distinct species within L. umbratilis. These findings suggest that the Amazonas and Negro Rivers may be acting as effective barriers, thus preventing gene flow between populations on opposite sides. Such findings have important implications for epidemiological studies, especially those related to vector competence and anthropophily, and for vector control strategies. In addition, L. umbratilis represents an interesting example in speciation studies.