Molecular identification of mosquito species

Revolutionizing Malaria Vector Control: The Importance of Accurate Species Identification through Enhanced Molecular Capacity

Many factors, such as the resistance to pesticides and a lack of knowledge of the morphology and molecular structure of malaria vectors, have made it more challenging to eradicate malaria in numerous malaria-endemic areas of the globe. The primary goal of this review is to discuss malaria vector control methods and the significance of identifying species in vector control initiatives. This was accomplished by reviewing methods of molecular identification of malaria vectors and genetic marker classification in relation to their use for species identification. Due to its specificity and consistency, molecular identification is preferred over morphological identification of malaria vectors. Enhanced molecular capacity for species identification will improve mosquito characterization, leading to accurate control strategies/treatment targeting specific mosquito species, and thus will contribute to malaria eradication. It is crucial for disease epidemiology and surveillance to accurately identify the Plasmodium spp. that are causing malaria in patients. The capacity for disease surveillance will be significantly increased by the development of more accurate, precise, automated, and high-throughput diagnostic techniques. In conclusion, although morphological identification is quick and achievable at a reduced cost, molecular identification is preferred for specificity and sensitivity. To achieve the targeted malaria elimination goal, proper identification of vectors using accurate techniques for effective control measures should be prioritized.

Ecological Niche Modelling Approaches: Challenges and Applications in Vector-Borne Diseases

Vector-borne diseases (VBDs) pose a major threat to human and animal health, with more than 80% of the global population being at risk of acquiring at least one major VBD. Being profoundly affected by the ongoing climate change and anthropogenic disturbances, modelling approaches become an essential tool to assess and compare multiple scenarios (past, present and future), and further the geographic risk of transmission of VBDs. Ecological niche modelling (ENM) is rapidly becoming the gold-standard method for this task. The purpose of this overview is to provide an insight of the use of ENM to assess the geographic risk of transmission of VBDs. We have summarised some fundamental concepts and common approaches to ENM of VBDS, and then focused with a critical view on a number of crucial issues which are often disregarded when modelling the niches of VBDs. Furthermore, we have briefly presented what we consider the most relevant uses of ENM when dealing with VBDs. Niche modelling of VBDs is far from being simple, and there is still a long way to improve. Therefore, this overview is expected to be a useful benchmark for niche modelling of VBDs in future research.

Morphological and Molecular Characterization Using Genitalia and CoxI Barcode Sequence Analysis of Afrotropical Mosquitoes with Arbovirus Vector Potential

Potential arboviral Afrotropical mosquito vectors are underrepresented in public databases of CoxI barcode sequences. Furthermore, available CoxI sequences for many species are often not associated with voucher specimens to match the corresponding fine morphological characterization of specimens. Hence, this study focused on the characterization of Culicine mosquitoes from South Africa, Mozambique, and Angola and their classification using a complementary approach including a morphological analysis of specimens’ genitalia and phylogenetic study based on the analysis of CoxI barcode sequences using maximum likelihood and Bayesian phylogenetic inference methods, alongside Median-Joining Network and PCOORD analyses. Overall, 800 mosquitoes (652 males and 148 females) from 67 species, were analyzed. Genitalia from 663 specimens allowed the identification of 55 species of 10 genera. A total of 247 CoxI partial gene sequences corresponding to 65 species were obtained, 11 of which (Aedes capensis, Ae. mucidus, Culex andersoni, Cx. telesilla, Cx. inconspicuosus, Eretmapodites subsimplicipes, Er. quinquevittatus, Ficalbia uniformis, Mimomyia hispida, Uranotaenia alboabdominalis, and Ur. mashonaensis) are, to the best of our knowledge, provided here for the first time. The presence of Cx. pipiens ecotypes molestus and pipiens and their hybrids, as well as Cx. infula, is newly reported in the Afrotropical region. The rates of correct sequence identification using BOLD and BLASTn (≥95% identity) were 64% and 53%, respectively. Phylogenetic analysis revealed that, except for subgenus Eumelanomyia of Culex, there was support for tribes Aedini, Culicini, Ficalbiini, and Mansoniini. A divergence >2% was observed in conspecific sequences, e.g., Aedeomyia africana, Ae. cumminsii, Ae. unilineatus, Ae. metallicus, Ae. furcifer, Ae. caballus, and Mansonia uniformis. Conversely, sequences from groups and species complexes, namely, Ae. simpsoni, Ae. mcintoshi, Cx. bitaeniorhynchus, Cx. simpsoni, and Cx. pipiens were insufficiently separated. A contribution has been made to the barcode library of Afrotropical mosquitoes with associated genitalia morphological identifications.

Detection and persistence of environmental DNA (eDNA) of the different developmental stages of a vector mosquito, Culex pipiens pallens

Preventing mosquito-borne infectious diseases requires that vector mosquitoes are monitored and controlled. Targeting immature mosquitoes (eggs, larvae, and pupae), which have less mobility than adults, is an effective management approach. However, conducting these surveys is often difficult due to the limitations of morphological classification and survey costs. The application of environmental DNA (eDNA) analysis can solve these issues because it allows easy estimation of species distribution and morphology-independent species identification. Although a few previous studies have reported mosquito eDNA detection, there is a gap in knowledge regarding the dynamics related to the persistence of immature mosquito eDNA. We used Culex pipiens pallens, a vector of West Nile fever, as a model species. First, we developed a species-specific detection assay and confirmed its specificity using in silico and in vitro tests. Next, we conducted laboratory experiments using breeding tanks. Water samples were collected at each developmental stage. In addition, water samples were collected daily until the seventh day after emergence from the pupae. We quantified eDNA using real-time PCR with the developed assay to investigate the dynamics of mosquito eDNA. The specificity of the developed assay was confirmed by in silico and in vitro tests. Mosquito eDNA was detected at all developmental stages and detected up to seven days after emergence of pupae. In particular, high concentrations of eDNA were detected immediately after hatching from eggs and after emergence from pupae. Highly frequent positive eDNA signals were continuously detected between egg hatching and pupa hatching. Mosquito eDNA was detected immediately after the eggs were introduced, and eDNA-positive detections continued until pupae emergence, suggesting that eDNA analysis is useful for monitoring mosquito larvae. In the future, monitoring immature mosquitoes using eDNA analysis will contribute to prevent mosquito-borne infectious diseases.

Development of microsatellite markers for population genetics of biting midges and a potential tool for species identification of Culicoides sonorensis Wirth & Jones

Background

Proper vector surveillance relies on the ability to identify species of interest accurately and efficiently, though this can be difficult in groups containing cryptic species. Culicoides Latreille is a genus of small biting flies responsible for the transmission of numerous pathogens to a multitude of vertebrates. Regarding pathogen transmission, the C. variipennis species complex is of particular interest in North America. Of the six species within this group, only C. sonorensis Wirth & Jones is a proven vector of bluetongue virus and epizootic hemorrhagic disease virus. Unfortunately, subtle morphological differences, cryptic species, and mitonuclear discordance make species identification in the C. variipennis complex challenging. Recently, single-nucleotide polymorphism (SNP) analysis enabled discrimination between the species of this group; however, this demanding approach is not practical for vector surveillance.

Methods

The aim of the current study was to develop a reliable and affordable way of distinguishing between the species within the C. variipennis complex, especially C. sonorensis. Twenty-five putative microsatellite markers were identified using the C. sonorensis genome and tested for amplification within five species of the C. variipennis complex. Machine learning was then used to determine which markers best explain the genetic differentiation between species. This led to the development of a subset of four and seven markers, which were also tested for species differentiation.

Results

A total of 21 microsatellite markers were successfully amplified in the species tested. Clustering analyses of all of these markers recovered the same species-level identification as the previous SNP data. Additionally, the subset of seven markers was equally capable of accurately distinguishing between the members of the C. variipennis complex as the 21 microsatellite markers. Finally, one microsatellite marker (C508) was found to be species-specific, only amplifying in the vector species C. sonorensis among the samples tested.

Conclusions

These microsatellites provide an affordable way to distinguish between the sibling species of the C. variipennis complex and could lead to a better understanding of the species dynamics within this group. Additionally, after further testing, marker C508 may allow for the identification of C. sonorensis with a single-tube assay, potentially providing a powerful new tool for vector surveillance in North America.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05189-8.

Identification of mosquitoes (Diptera: Culicidae) from Mexico State, Mexico using morphology and COI DNA barcoding

Mosquitoes are commonly identified to species level using morphological traits, but complementary methods for identification are often necessary when specimens are collected as immature stages, stored inadequately, or when delineation of species complexes is problematic. DNA-barcoding using the mitochondrial cytochrome c oxidase subunit 1 (COI) gene is one such tool used for the morphological identification of species. A comprehensive entomological survey of mosquito species in Mexico State identified by COI DNA barcoding and morphology is documented in this paper. Specimens were collected from all the physiographic provinces in Mexico State between 2017 and 2019. Overall, 2,218 specimens were collected from 157 localities representing both subfamilies Anophelinae and Culicinae. A species checklist that consists of 6 tribes, 10 genera, 20 subgenera, and 51 species, 35 of which are new records for Mexico State, is provided. Three hundred and forty-two COI sequences of 46 species were analysed. Mean intraspecific and interspecific distances ranged between 0% to 3.9% and from 1.2% to 25.3%, respectively. All species groups were supported by high bootstraps values in a Neighbour-Joining analysis, and new COI sequences were generated for eight species: Aedes chionotum Zavortink, Ae. vargasi Schick, Ae. gabriel Schick, Ae. guerrero Berlin, Ae. ramirezi Vargas and Downs, Haemagogus mesodentatus Komp and Kumm, Culex restrictor Dyar and Knab, and Uranotaenia geometrica Theobald. This study provides a detailed inventory of the Culicidae from Mexico State and discusses the utility of DNA barcoding as a complementary tool for accurate mosquito species identification in Mexico.

Integrated systematics of Anopheles subpictus (Diptera: Culicidae) in the Oriental Region, with emphasis on forms in Thailand and Sulawesi, Indonesia

The Anopheles subpictus complex consists of four species informally designated, based on fixed inversions of polytene chromosomes and morphology, as species A, B, C and D in India. However, recent studies revealed the presence of only species A and B in Sri Lanka. Little is known about the specific identity of the taxon in other countries in Asia. This paper reports the results of a molecular and morphological study of An. subpictus in Thailand and South Sulawesi, Indonesia. The maxillary palpi of most females from Thailand have the apical pale band longer than the subapical dark band, seta 7-I of pupae branched and short, and eggs with 18-25 float ridges. These characters do not agree with those described for species A, B, C and D in India. The females of An. subpictus from South Sulawesi usually have the subapical dark band of the maxillary palpus equal in length to the apical pale band. Phylogenetic analyses of sequences of the internal transcribed spacer 2 (ITS2) region of rDNA and the cytochrome c oxidase subunit I (COI) gene of mtDNA of specimens from Thailand, and South Sulawesi, and from various localities in GenBank, were conducted. ITS2 sequences of specimens from all localities in Thailand were identical, except for a small divergence in specimen from Phang Nga Province. Three distinct COI clades were detected in specimens from Chiang Mai Province in northern Thailand. However, crossing experiments between the three clades revealed no genetic incompatibility, suggesting that they were conspecific. ITS2 and COI sequences of most specimens from Thailand fell in clades other than those of An. subpictus species A and B and An. subpictus from Indonesia (East Nusa Tenggara, Java, South Sulawesi) and the Philippines. ITS2 sequences from South Sulawesi and East Nusa Tenggara were very similar, and fell in a clade consisting of specimen from Phang Nga in southern Thailand and sequences of some specimens from Cambodia and Vietnam, but their COI sequences were distinct. DNA sequences and morphological differences suggest the presence of two species within An. subpictus in Thailand, and more than one species in Indonesia.

Molecular and morphological evidence for sibling species within Anopheles (Anopheles) lindesayi Giles (Diptera: Culicidae) in Bhutan

This paper reports the results of a comparative molecular and morphological study of An. lindesayi collected from various districts of Bhutan and An. l. cameronensis from Thailand, compared with GenBank accessions and publications for An. l. japonicus from Japan, South Korea and China, An. l. pleccau from Taiwan, and An. lindesayi from India. Phylogenetic analyses based on ribosomal (ITS2) and mitochondrial (COI) DNA sequences using the Maximum Likelihood method revealed five genetically distinct clades (A, B, C, D and E) in Bhutan. Specimens in Clade A correspond to the original description of An. lindesayi, particularly in wing markings, the pattern of basal pale scales on the hindfemur and the single seta 4-C of larvae, and their COI sequences were closely related to one Indian sequence. Larvae of Clades B, C, D and E are similar in having seta 4-C branched rather than single. The adults of Clades C, D and E (B not available) are distinguishable from those of Clade A and other subspecies. Specimens of Clade C are unique in having a long pale spot on wing vein R and the subcosta, scattered pale scales on several veins and a dark spot at the tip of vein R₂. The adults of Clades D and E are similar in having a dark spot at the tip of vein R₂ and no scattered pale scales on all other veins. We provisionally recognize mosquitoes of Clades A, B, C, D and E as species A, B, C, D and E, respectively, of the Lindesayi Complex. Species A is An. lindesayi sensu stricto and the others are unnamed species. Concomitantly, the previous concept of the "Lindesayi Complex", which included An. lindesayi, An. menglangensis, An. nilgiricus and An. wellingtonianus, is now recognized as the Lindesayi Subgroup of the Lindesayi Group (Anopheles Series, subgenus Anopheles) with the five sibling species of An. lindesayi comprising a more apposite Lindesayi Complex within the subgroup.

Molecular identification of mosquitoes of the Anopheles maculatus group of subgenus Cellia (Diptera: Culicidae) in the Indonesian Archipelago

This study reports the molecular differentiation of females of Anopheles maculatus s.l. collected in eight localities on five islands in the Indonesian Archipelago: Hargowilis and Hargotirto villages of Central Java Province, North Kalimantan Province, Sabang off the northern tip of Sumatra Province, Sumba Island of East Nusa Tenggara Province and Sulawesi Province. Analyses based on rDNA (ITS2 and D3) and mtDNA (COII) sequences revealed the presence of An. greeni for the first time in North Kalimantan, and at least one novel (previously unrecognized) species of the Maculatus Group in Central Java (Hargowilis). Despite the similarity of rDNA markers of specimens of An. maculatus s.l. from Central Java and Sulawesi, their COII sequences are highly divergent (3.3%), which might indicate the presence of a further new species. Specimens of An. maculatus s.l. from the other localities had identical rDNA sequences to most An. maculatus s.s. from mainland Southeast Asia, but moderate divergence in their COII sequences (1.2-2.1%). The latter might indicate there are further novel species within the Maculatus Complex. However, as the divergence at COII may be the result of geographical structuring within species related to the historical biogeography of the region, further studies are needed to shed light on this possibility.

Is Genetic Continuity Between Anopheles sinensis (Diptera: Culicidae) and its Sibling Species Due to Gene Introgression or Incomplete Speciation?

The Anopheles mosquito Hyrcanus Group is widely distributed geographically across both Palearctic and Oriental regions and comprises 26 valid species. Although the species Anopheles sinensis Wiedemann (1828) is the most common in China and has a low potential vector rank, it has nevertheless long been thought to be an important natural malaria vector within the middle and lower reaches of the Yangtze River. A number of previous research studies have found evidence to support the occurrence of natural hybridization between An. sinensis and Anopheles kleini Rueda, 2005 (a competent malaria vector). We, therefore, collected a sample series of An. sinensis and morphologically similar species across China and undertook ribosomal and mitochondrial DNA analyses in order to assess genetic differentiation (Fst) and gene flow (Nm) amongst different groups. This enabled us to evaluate divergence times between morphologically similar species using the cytochrome oxidase I (COI) gene. The results of this study reveal significant genetic similarities between An. sinensis, An. kleini, and Anopheles belenrae Rueda, 2005 and therefore imply that correct molecular identifications will require additional molecular markers. As results also reveal the presence of gene flow between these three species, their taxonomic status will require further work. Data suggest that An. kleini is the most basal of the three species, while An. sinensis and An. belenrae share the closest genetic relationship.

Molecular Differentiation of the African Yellow Fever Vector Aedes bromeliae (Diptera: Culicidae) from Its Sympatric Non-vector Sister Species, Aedes lilii

INTRODUCTION

Yellow fever continues to be a problem in sub-Saharan Africa with repeated epidemics occurring. The mosquito Aedes bromeliae is a major vector of yellow fever, but it cannot be readily differentiated from its non-vector zoophilic sister species Ae. lilii using morphological characters. Genetic differences have been reported between anthropophilic Ae. bromeliae and zoophilic Ae. lilii and between forest and domestic populations. However, due to the application of different molecular markers and non-overlapping populations employed in previous studies, interpretation of species delimitation is unclear.

METHODOLOGY/PRINCIPLE FINDINGS

DNA sequences were generated from specimens of Ae. simpsoni s.l. from the Republic of Benin, Tanzania and Uganda for two nuclear genes apolipophorin 2 (apoLp2) and cytochrome p450 (CYPJ92), the ribosomal internal transcribed spacer region (ITS) and the mitochondrial cytochrome c oxidase (COI) barcoding region. Nuclear genes apoLp2 and CYPJ92 were unable to differentiate between species Ae. bromeliae and Ae. lilii due to ancestral lineage sorting, while ITS sequence data provided clear topological separation on a phylogeny. The standard COI barcoding region was shown to be subject to species introgression and unable to clearly distinguish the two taxa. Here we present a reliable direct PCR-based method for differentiation of the vector species Ae. bromeliae from its isomorphic, sympatric and non-biomedically important sister taxon, Ae. lilii, based on the ITS region. Using molecular species verification, we describe novel immature habitats for Ae. lilii and report both sympatric and allopatric populations. Whereas only Ae. lilii is found in the Republic of Benin and only Ae. bromeliae in Tanzania, both species are sympatric in Uganda.

CONCLUSIONS/SIGNIFICANCE

Our accurate identification method will allow informed distribution and detailed ecological studies that will facilitate assessment of arboviral disease risk and development of future targeted vector control.

Deep Interisland Genetic Divergence in the Macaronesian Endemic Mosquito Ochlerotatus eatoni (Diptera: Culicidae), Indication of Cryptic Species

Ochlerotatus eatoni (Edwards, 1916) is a species endemic to Canary and Madeira Islands that, based on morphology, is considered to be single species. Mitochondrial 16S rRNA and cytochrome oxidase I (COI) sequence data demonstrate that the populations from Tenerife and Madeira Islands are highly differentiated (F(ST) = 0.93). The phylogenetic analysis also separates the two populations into two highly distinct groups. The sharp mitochondrial genetic differentiation between islands is congruent with the published nuclear (allozyme) data. However, mtDNA data did not reveal any significant genetic differentiation within islands. Extreme interisland genetic divergence, but lack of morphological variation, is indicative of the existence of cryptic species. I suggest the elevation of populations to at least incipient species status, designating the populations from Tenerife and Madeira Islands as Oc. eatoni. hewitti and Oc. eatoni. krimbasi, respectively.

Phylogenetic analysis of the GST family in Anopheles (Nyssorhynchus) darlingi

Anopheles darlingi Root, 1926 and Anopheles gambiae (Diptera: Culicidae) are the most important human malaria vectors in South America and Africa, respectively. The two species are estimated to have diverged 100 million years ago. Studies on the phylogenetics and evolution of gene sequences, such as glutathione S-transferase (GST) in disease-transmitting mosquitoes are scarce. The sigma class GST (KC890767) from the transcriptome of An. darlingi captured in the Brazilian Amazon was studied by in silico hybridization, and mapped to chromosome 3 of An. gambiae. The sigma class GST of An. darlingi was used for phylogenetic analyses to understand the GST base composition of the most recent common ancestor between An. darlingi, Anopheles gambiae, Aedes aegypti and Culex quinquefasciatus. The GST (KC890767) of An. darlingi was studied to generate the main divergence branches using a Neighbor-Joining and bootstrapping approaches to confirm confidence levels on the tree nodes that separate the An. darlingi and other mosquito species. The results showed divergence between An. gambiae, Ae. Aegypti, Cx. quinquefasciatus, and Phlebotomus papatasi as outgroup, and the homology relationship between sigma class GST of An. darlingi and GSTS1_1 gene of An. gambiae was valuable for phylogenetic and evolutionary studies.

Analyzing mosquito (Diptera: culicidae) diversity in Pakistan by DNA barcoding

BACKGROUND

Although they are important disease vectors mosquito biodiversity in Pakistan is poorly known. Recent epidemics of dengue fever have revealed the need for more detailed understanding of the diversity and distributions of mosquito species in this region. DNA barcoding improves the accuracy of mosquito inventories because morphological differences between many species are subtle, leading to misidentifications.

METHODOLOGY/PRINCIPAL FINDINGS

Sequence variation in the barcode region of the mitochondrial COI gene was used to identify mosquito species, reveal genetic diversity, and map the distribution of the dengue-vector species in Pakistan. Analysis of 1684 mosquitoes from 491 sites in Punjab and Khyber Pakhtunkhwa during 2010-2013 revealed 32 species with the assemblage dominated by Culex quinquefasciatus (61% of the collection). The genus Aedes (Stegomyia) comprised 15% of the specimens, and was represented by six taxa with the two dengue vector species, Ae. albopictus and Ae. aegypti, dominant and broadly distributed. Anopheles made up another 6% of the catch with An. subpictus dominating. Barcode sequence divergence in conspecific specimens ranged from 0-2.4%, while congeneric species showed from 2.3-17.8% divergence. A global haplotype analysis of disease-vectors showed the presence of multiple haplotypes, although a single haplotype of each dengue-vector species was dominant in most countries. Geographic distribution of Ae. aegypti and Ae. albopictus showed the later species was dominant and found in both rural and urban environments.

CONCLUSIONS

As the first DNA-based analysis of mosquitoes in Pakistan, this study has begun the construction of a barcode reference library for the mosquitoes of this region. Levels of genetic diversity varied among species. Because of its capacity to differentiate species, even those with subtle morphological differences, DNA barcoding aids accurate tracking of vector populations.

Evidence to support natural hybridization between Anopheles sinensis and Anopheles kleini (Diptera: Culicidae): possibly a significant mechanism for gene introgression in sympatric populations

Background

Malaria caused by Plasmodium vivax is still a public health problem in the Republic of Korea (ROK), particularly regarding the recent re-emergence of this malarial species near the demilitarized zone in northwestern Paju City, Gyeonggi-do Province. Currently, at least 4 species (An. kleini, An. pullus, An. belenrae and An. lesteri) of the Hyrcanus Group are reported as possible natural vectors of vivax malaria in the ROK, and An. sinensis, which is the most dominant species, has long been incriminated as an important natural vector of this P. vivax. However, An. sinensis was ranked recently as a low potential vector. According to the discovery of natural hybrids between An. sinensis (a low potential vector for P. vivax) and An. kleini (a high potential vector for P. vivax) in Paju City, intensive investigation of this phenomenon is warranted under laboratory conditions.

Methods

Mosquitoes were collected during 2010-2012 from Paju City, ROK. Hybridization experiments used iso-female line colonies of these anophelines together with DNA analysis of ribosomal DNA [second internal transcribed spacer (ITS2)] and mitochondrial DNA [cytochrome c oxidase subunit I (COI)] of the parental colonies, F₁-hybrids and repeated backcross progenies were performed intensively by using a PCR-based assay and pyrosequencing technology.

Results

The results from hybridization experiments and molecular investigations revealed that the mitochondrial COI gene was introgressed from An. sinensis into An. kleini. The An. sinensis progenies obtained from consecutive repeated backcrosses in both directions, i.e., F₂-₁₁ progeny [(An. sinensis x An. kleini) x An. sinensis] and F₃-₅ progeny [(An. kleini x An. sinensis) x An. kleini] provided good supportive evidence.

Conclusions

This study revealed introgression of the mitochondrial COI gene between An. sinensis and An. kleini through consecutive repeated backcrosses under laboratory conditions. This new body of knowledge will be emphasized in reliable promising strategies in order to replace the population of An. kleini as a high potential vector for P. vivax, with that of a low potential vector, An. sinensis, through the mechanism of gene introgression in nature.

Intragenomic variation in the second internal transcribed spacer of the ribosomal DNA of species of the genera Culex and Lutzia (Diptera: Culicidae)

Culex is the largest genus of Culicini and includes vectors of several arboviruses and filarial worms. Many species of Culex are morphologically similar, which makes their identification difficult, particularly when using female specimens. To aid evolutionary studies and species distinction, molecular techniques are often used. Sequences of the second internal transcribed spacer (ITS2) of ribosomal DNA (rDNA) from 16 species of the genus Culex and one of Lutzia were used to assess their genomic variability and to verify their applicability in the phylogenetic analysis of the group. The distance matrix (uncorrected p-distance) that was obtained revealed intragenomic and intraspecific variation. Because of the intragenomic variability, we selected ITS2 copies for use in distance analyses based on their secondary structures. Neighbour-joining topology was obtained with an uncorrected p-distance. Despite the heterogeneity observed, individuals of the same species were grouped together and correlated with the current, morphology-based classification, thereby showing that ITS2 is an appropriate marker to be used in the taxonomy of Culex.

The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis

Background

This is the second in a series of three articles documenting the geographical distribution of 41 dominant vector species (DVS) of human malaria. The first paper addressed the DVS of the Americas and the third will consider those of the Asian Pacific Region. Here, the DVS of Africa, Europe and the Middle East are discussed. The continent of Africa experiences the bulk of the global malaria burden due in part to the presence of the An. gambiae complex. Anopheles gambiae is one of four DVS within the An. gambiae complex, the others being An. arabiensis and the coastal An. merus and An. melas. There are a further three, highly anthropophilic DVS in Africa, An. funestus, An. moucheti and An. nili. Conversely, across Europe and the Middle East, malaria transmission is low and frequently absent, despite the presence of six DVS. To help control malaria in Africa and the Middle East, or to identify the risk of its re-emergence in Europe, the contemporary distribution and bionomics of the relevant DVS are needed.

Results

A contemporary database of occurrence data, compiled from the formal literature and other relevant resources, resulted in the collation of information for seven DVS from 44 countries in Africa containing 4234 geo-referenced, independent sites. In Europe and the Middle East, six DVS were identified from 2784 geo-referenced sites across 49 countries. These occurrence data were combined with expert opinion ranges and a suite of environmental and climatic variables of relevance to anopheline ecology to produce predictive distribution maps using the Boosted Regression Tree (BRT) method.

Conclusions

The predicted geographic extent for the following DVS (or species/suspected species complex*) is provided for Africa: Anopheles (Cellia) arabiensis, An. (Cel.) funestus*, An. (Cel.) gambiae, An. (Cel.) melas, An. (Cel.) merus, An. (Cel.) moucheti and An. (Cel.) nili*, and in the European and Middle Eastern Region: An. (Anopheles) atroparvus, An. (Ano.) labranchiae, An. (Ano.) messeae, An. (Ano.) sacharovi, An. (Cel.) sergentii and An. (Cel.) superpictus*. These maps are presented alongside a bionomics summary for each species relevant to its control.

Mitochondrial DNA variation in the malaria vector Anopheles minimus across China, Thailand and Vietnam: evolutionary hypothesis, population structure and population history

The effects of Pleistocene environmental fluctuations on the distribution and diversity of organisms in Southeast Asia are much less well known than in Europe and North America. In these regions, the combination of palaeoenvironmental reconstruction and inferences about population history from genetic data has been very powerful. In Southeast Asia, mosquitoes are good candidates for the genetic approach, with the added benefit that understanding the relative contributions of historical and current processes to population structure can inform management of vector species. Genetic variation among populations of Anopheles minimus was examined using 144 mtDNA COII sequences from 23 sites in China, Thailand and Vietnam. Haplotype diversity was high, with two distinct lineages that have a sequence divergence of over 2% and exhibit different geographical distributions. We compare alternative hypotheses concerning the origin of this pattern. The observed data deviate from the expectations based on a single-panmictic population with or without growth, or a stable but spatially structured population. However, they can be readily accommodated by a model of past fragmentation into eastern and western refugia, followed by growth and range expansion. This is consistent with the palaeoenvironmental reconstructions currently available for the region.

Resolving genetic diversity in Australasian Culex mosquitoes: incongruence between the mitochondrial cytochrome c oxidase I and nuclear acetylcholine esterase 2

Insects that vector pathogens are under constant surveillance in Australasia although the repertoire of genetic markers to distinguish what are often cryptic mosquito species remains limited. We present a comparative assessment of the second exon-intron region of the acetylcholine esterase 2 gene (ace-2) and the mitochondrial DNA cytochrome c oxidase I (COI) using two closely related Australasia mosquitoes Culex annulirostris and Culex palpalis. The COI revealed eight divergent lineages of which four were confirmed with the ace-2. We dissect out the nuclear chromosomal haplotypes of the ace-2 as well as the exon-intron regions by assessing the protein's tertiary structure to reveal a hypervariable 5'-exon that forms part of an external protein loop and displays a higher polymorphic rate than the intron. We retrace the evolutionary history of these mosquitoes by phylogenetic inference and by testing different evolutionary hypotheses. We conclude that DNA barcoding using COI may overestimate the diversity of Culex mosquitoes in Australasia and should be applied cautiously with support from the nuclear DNA such as the ace-2. Together the COI and ace-2 provide robust evidence for distinct cryptic Culex lineages--one of which correlates exactly with the southern limit of Japanese encephalitis virus activity in Australasia.

Complex population history of two Anopheles dirus mosquito species in Southeast Asia suggests the influence of Pleistocene climate change rather than human-mediated effects

Anopheles dirus and Anopheles baimaii are closely related species which feed on primates, particularly humans, and transmit malaria in the tropical forests of mainland Southeast Asia. Here, we report an in-depth phylogeographic picture based on 269 individuals from 21 populations from mainland Southeast Asia. Analysis of 1537 bp of mtDNA sequence revealed that the population history of A. baimaii is far more complex than previously thought. An old expansion (pre-300 kyr BP) was inferred in northern India/Bangladesh with a wave of south-eastwards expansion arriving at the Thai border (ca 135-173 kyr BP) followed by leptokurtic dispersal very recently (ca 16 kyr BP) into peninsular Thailand. The long and complex population history of these anthropophilic species suggests their expansions are not in response to the relatively recent (ca 40 kyr BP) human expansions in mainland Southeast Asia but, rather, fit well with our understanding of Pleistocene climatic change there.

Classification and identification of mosquito species using artificial neural networks

An artificial neural network method is presented for classification and identification of Anopheles mosquito species based on the internal transcribed spacer2 (ITS2) data of ribosomal DNA string. The method is implemented in two different multi-layered feed-forward neural network model forms, namely, multi-input single-output neural network (MISONN) and multi-input multi-output neural network (MIMONN). A number of data sequences in varying sizes of different Anopheline malarial vectors and their corresponding species coding are employed to develop the neural network models. The classification efficiency of the network models for untrained data sequences is evaluated in terms of quantitative performance criteria. The results demonstrate the efficiency of the neural network models to extract the genetic information in ITS2 sequences and to adapt to new data. The method of MISONN is found to exhibit superior performance over MIMONN in distinguishing and identification of the mosquito vectors.

Bionomics, taxonomy, and distribution of the major malaria vector taxa of Anopheles subgenus Cellia in Southeast Asia: an updated review

There is high diversity of Anopheles mosquitoes in Southeast Asia and the main vectors of malaria belong to complexes or groups of species that are difficult or impossible to distinguish due to overlapping morphological characteristics. Recent advances in molecular systematics have provided simple and reliable methods for unambiguous species identification. This review summarizes the latest information on the seven taxonomic groups that include principal malaria vectors in Southeast Asia, i.e. the Minimus, Fluviatilis, Culicifacies, Dirus, Leucosphyrus, and Sundaicus Complexes, and the Maculatus Group. Main issues still to be resolved are highlighted. The growing knowledge on malaria vectors in Southeast Asia has implications for vector control programs, the success of which is highly dependant on precise information about the biology and behavior of the vector species. Acquisition of this information, and consequently the application of appropriate, sustainable control measures, depends on our ability to accurately identify the specific vectors.

Molecular characterization and species identification of the Anopheles dirus and An. minimus complexes in north-east India using r-DNA ITS-2

The sibling species composition of the Anopheles minimus and Anopheles dirus complexes is poorly known in the highly malarious north-eastern region of India where these two vector taxa are accountable for most of the malaria transmission among 30.7 million inhabitants. Prevalent members of these two complexes in this part of India were identified using sequences for the second internal transcribed spacer (ITS2) of ribosomal DNA. Anopheles baimaii (species D) of the An. dirus complex and An. minimus s.s. (species A) of the An. minimus complex were detected in Arunachal Pradesh, Assam, Meghalaya and Nagaland states. No intraspecific variation was observed in the ITS2 sequence (479bp) of An. baimaii whereas a single substitution was detected in the ITS2 sequence (372bp) of An. minimus from Nagaland state.

Genetic diversity and molecular identification of mosquito species in the Anopheles maculatus group using the ITS2 region of rDNA

The species diversity and genetic structure of mosquitoes belonging to the Anopheles maculatus group in Southeast Asia were investigated using the internal transcribed spacer 2 (ITS2) of ribosomal DNA (rDNA). A molecular phylogeny indicates the presence of at least one hitherto unrecognised species. Mosquitoes of chromosomal form K from eastern Thailand have a unique ITS2 sequence that is 3.7% divergent from the next most closely related taxon (An. sawadwongporni) in the group. In the context of negligible intraspecific variation at ITS2, this suggests that chromosomal form K is most probably a distinct species. Although An. maculatus sensu stricto from northern Thailand and southern Thailand/peninsular Malaysia differ from each other in chromosomal banding pattern and vectorial capacity, no intraspecific variation was observed in the ITS2 sequences of this species over this entire geographic area despite an extensive survey. A PCR-based identification method was developed to distinguish five species of the group (An. maculatus, An. dravidicus, An. pseudowillmori, An. sawadwongporni and chromosomal form K) to assist field-based studies in northwestern Thailand. Sequences from 187 mosquitoes (mostly An. maculatus and An. sawadwongporni) revealed no intraspecific variation in specimens from Thailand, Cambodia, mainland China, Malaysia, Taiwan and Vietnam, suggesting that this identification method will be widely applicable in Southeast Asia. The lack of detectable genetic structure also suggests that populations of these species are either connected by gene flow and/or share a recent common history.

Sequence analysis of the rDNA internal transcribed spacer 2 and polymerase chain reaction identification of Anopheles fluminensis (Diptera: Culicidae: Anopheles) in Bolivia

Anopheles fluminensis Root is a member of the Arribalzagia Series in the subgenus Anopheles. We report the first record of this species in the department of Cochabamba, Bolivia. This species was sampled from two locations in the foothills of the eastern Andes Mountains within the Chapare Valley. Larvae were collected in fast-flowing, shaded streams at the edges of rocky pools. We provide the first sequence data for the rDNA of An. fluminensis, a partial sequence of the 5.8S and the internal transcribed spacer 2 (ITS2). The ITS2 of An. fluminensis, sequenced from two individuals at one site, was at least 596 bp, had 56.5% GC, and included three large repeats (approximately equal to 125 bp each). We describe a polymerase chain reaction protocol and species-specific primers for identifying this species in the Chapare Valley, Bolivia.

Identification of Anopheles (Nyssorhynchus) albitarsis complex species (Diptera: Culicidae) using rDNA internal transcribed spacer 2-based polymerase chain reaction primes

Anopheles (Nyssorhynchus) marajoara is a proven primary vector of malaria parasites in Northeast Brazil, and An. deaneorum is a suspected vector in Western Brazil. Both are members of the morphologically similar Albitarsis Complex, which also includes An. albitarsis and an undescribed species, An. albitarsis "B". These four species were recognized and can be identified using random amplified polymorphic DNA (RAPD) markers, but various other methodologies also point to multiple species under the name An. albitarsis. We describe here a technique for identification of these species employing polymerase chain reaction (PCR) primers based on ribosomal DNA internal transcribed spacer 2 (rDNA ITS2) sequence. Since this method is based on known sequence it is simpler than the sometimes problematical RAPD-PCR. Primers were tested on samples previously identified using RAPD markers with complete correlation.

Crossing experiments of Anopheles minimus species C and putative species E

In the Anopheles minimus complex, 2 sibling species (A and C) are generally accepted. Recently, a 3rd species, provisionally designated An. minimus species E, has been described from the Ryukyu Archipelago, Japan, based on crossing experiments (A and E), DNA analysis, mitotic karyotypes, and some morphological characteristics. The present study reports the results of crossing experiments between species C and putative species E. Hybridization between the progeny of An. minimus species C from Thailand and putative species E from Japan revealed postzygotic genetic incompatibility. Although F1 hybrid progeny were obtained from both directions of crosses, the hybrid males from C female x E male crosses were completely sterile, with atrophied testes and accessory glands. In addition, the external terminalia of all of these males never completely rotated and the males failed to copulate by artificial mating. In E female x C male crosses, the hybrid males showed partially sterile testes in which most spermatozoa were abnormal (enlarged head) and inactive, and they had very little success in inseminating females. The salivary gland polytene chromosomes of F1 hybrid larvae from species C female x species E male showed a fixed heterozygous inversion on the 3L arm. Those F1 hybrids from species E female x species C male showed partial asynapsis on identified arms (2R and 3L) and a fixed heterozygous inversion on the 3R arm. When the F1 hybrid females from both directions of crosses were backcrossed with either C or E males, they produced male progeny with abnormal spermatozoa. Study of mating behavior in a 30 x 30 x 30-cm cage showed that the C males failed to mate with either C or E females, indicating that species C males cannot breed in confined spaces (lack stenogamy). Putative species E males had little success in inseminating species C females. This study provides strong evidence of genetic incompatibility between An. minimus species C and putative species E, supporting previous data that species E is a distinct species in the An. minimus complex.

Dentification of four vectors of human Plasmodium spp. by multiplex PCR: Anopheles rangeli, An. strodei, An. triannulatus, and An. trinkae (Diptera: Culicidae: Nyssorhynchus)

One of the major obstacles for studies of the biology, ecology, and behavior of Neotropical vectors of human Plasmodium has been the lack of reliable and efficient means of identifying many species. Although the subgenus Nyssorhynchus includes most species responsible for human transmission in South America, there are no polymerase chain reaction (PCR)-based techniques for identifying members of this taxon. We describe the first multiplex PCR for identifying four species in the subgenus Nyssorhynchus that are vectors of human Plasmodium spp. Four species specific primers, together with a universal primer that anneals to the 5.8S rDNA region, produce amplicons of the internal transcribed spacer two with base pair sizes of 131,308,371, and 441 for An. triannulatus, An. trinkae, An. strodei, and An. rangeli, respectively.

Ribosomal DNA ITS2 sequences differentiate six species in the Anopheles crucians complex (Diptera: Culicidae)

Anopheles crucians Wiedemann (sensu lato) was investigated for the presence of cryptic species using rDNA ITS2 sequences. This complex of species presently contains the named species An. crucians, An. bradleyi King, and An. georgianus King. Adult female mosquitoes were collected at 28 sites in Alabama, Florida, Georgia, North Carolina, Mississippi, and Louisiana, resulting in 245 progeny broods. Species were identified using preliminary morphological characters, and the internal transcribed spacer two (ITS2) was amplified from all broods. The result was five distinct sizes of amplification product, and based on morphological characters, one of the size classes was suspected to consist of two species. All six putative species were then sequenced: five directly, and the sixth, because of extreme intragenomic (each individual with many variants) size variability, cloned. The ITS2 sequences were markedly distinct for all six species. Species designations and ITS2 sequence lengths (base pairs in parentheses) were A (461), B (1,000+), C (204), D (293), E (195), and An. bradleyi (208). Species B showed both large intraspecific and intragenomic sequence variability and is distinguished by having the longest ITS2 found so far in an Anopheles. Based on these data, we found that all species could be identified with polymerase chain reaction (PCR) using a mixture of four primers in a single reaction. Members of this complex were often found in sympatry, with the adults of five species collected at a single site in central Florida.

Molecular systematics of Anopheles: from subgenera to subpopulations

The century-old discovery of the role of Anopheles in human malaria transmission precipitated intense study of this genus at the alpha taxonomy level, but until recently little attention was focused on the systematics of this group. The application of molecular approaches to systematic problems ranging from subgeneric relationships to relationships at and below the species level is helping to address questions such as anopheline phylogenetics and biogeography, the nature of species boundaries, and the forces that have structured genetic variation within species. Current knowledge in these areas is reviewed, with an emphasis on the Anopheles gambiae model. The recent publication of the genome of this anopheline mosquito will have a profound impact on inquiries at all taxonomic levels, supplying better tools for estimating phylogeny and population structure in the short term, and ultimately allowing the identification of genes and/or regulatory networks underlying ecological differentiation, speciation, and vectorial capacity.

Genetic population structure and introgression in Anopheles dirus mosquitoes in South-east Asia

Genetic structure and species relationships were studied in three closely related mosquito species, Anopheles dirus A, C and D in Thailand using 11 microsatellite loci and compared with previous mitochondrial DNA (mtDNA) data on the same populations. All three species were well differentiated from each other at the microsatellite loci. Given the almost complete absence of mtDNA differentiation between An. dirus A and D, this endorses the previous suggestion of mtDNA introgression between these species. The high degree of differentiation between the northern and southern population of An. dirus C (RST = 0.401), in agreement with mtDNA data, is suggestive of incipient species. The lack of genetic structure indicated by microsatellites in four populations of An. dirus A across northern Thailand also concurs with mtDNA data. However, in An. dirus D a limited but significant level of structure was detected by microsatellites over ~400 km in northern Thailand, whereas the mtDNA detected no population differentiation over a much larger area (>1200 km). There is prior evidence for population expansion in the mtDNA. If this is due to a selective sweep originating in An. dirus D, the microsatellite data may indicate greater barriers to gene flow within An. dirus D than in species A. Alternatively, there may have been historical introgression of mtDNA and subsequent demographic expansion which occurred first in An. dirus D so enabling it to accumulate some population differentiation. In the latter case the lack of migration-drift equilibrium precludes the inference of absolute or relative values of gene flow in An. dirus A and D.

Population structure and population history of Anopheles dirus mosquitoes in Southeast Asia

Separating the confounding effects of long-term population history from gene flow can be difficult. Here, we address the question of what inferences about gene flow can be made from mitochondrial sequence data in three closely related species of mosquitoes, Anopheles dirus species A, C, and D, from southeast Asia. A total of 84 sequences of 923 bp of the mitochondrial cytochrome oxidase I gene were obtained from 14 populations in Thailand, Myanmar, and Bangladesh. The genealogy of sequences obtained from two populations of AN: dirus C indicates no contemporary gene flow between them. The F(ST) value of 0.421 therefore probably represents a recent common history, perhaps involving colonization events. Anopheles dirus A and D are parapatric, yet no differentiation was seen either within or between species. The starlike genealogy of their haplotypes, smooth unimodal mismatch distributions, and excess of low frequency mutations indicate population expansion in An. dirus A and D. This, rather than widespread gene flow, explains their low within-species F(ST) values (0.018 and 0.022). The greater genetic diversity of An. dirus D suggests that expansion occurred first in species D and subsequently in species A. The current geographical separation and low hybrid fitness of these species also argue against ongoing interspecific gene flow. They suggest instead either historical introgression of mtDNA from An. dirus D into species A followed by independent range expansions, or a selective sweep of mtDNA that originated in An. dirus D. While not excluding contemporary gene flow, historical population processes are sufficient to explain the data in An. dirus A and D. The genealogical relationships between haplotypes could not be used to make inferences of gene flow because of extensive homoplasy due to hypervariable sites and possibly also recombination. However, it is concluded that this approach, rather than the use of fixation indices, is required in the future to understand contemporary gene flow in these mosquitoes. The implications of these results for understanding gene flow in another important and comparable group of malaria vector mosquitoes in Africa, the An. gambiae complex, are also discussed.