An evaluation of the suitability of COI and COII gene variation for reconstructing the phylogeny of, and identifying cryptic species in, anopheline mosquitoes (Diptera Culicidae)

Review of the Mimeticus Subgroup of Culex (Culex) (Diptera: Culicidae), With an Assessment of the Specific Status of Three Nominal Species Described From India

The Mimeticus Subgroup of the subgenus Culex, genus Culex Linnaeus, is reviewed and revised. Mitochondrial COI sequences of species of the subgroup available in GenBank were analysed to verify the identity of voucher specimens. The specific status of three species described from India, Cx. gaugleri Suman, Cx. katezari Karlekar, Andrew & Deshpande, and Cx. kodaikanalensis Suman, Veer & Chandra, was examined and found to be unsupportable. Culex gaugleri and Cx. kodaikanalensis are synonyms of Cx. jacksoni based on the similarity of DNA sequences and dubious evidence of morphological differentiation. Culex katezari is recognized as a nomen dubium because the morphological description and molecular information provided by the authors are insufficient for its identification. A number of specimens registered in GenBank were found to be misidentified. In addition, available data indicate that the Mimeticus Subgroup may include at least three unrecognized species, referred to informally as Cx. mimeticus form B, Cx. mimulus form X, and Cx. mimulus form Y.

Myzorhynchus series of Anopheles mosquitoes as potential vectors of Plasmodium bubalis in Thailand

Ungulate malaria parasites and their vectors are among the least studied when compared to other medically important species. As a result, a thorough understanding of ungulate malaria parasites, hosts, and mosquito vectors has been lacking, necessitating additional research efforts. This study aimed to identify the vector(s) of Plasmodium bubalis. A total of 187 female mosquitoes (133 Anopheles spp., 24 Culex spp., 24 Aedes spp., and 6 Mansonia spp. collected from a buffalo farm in Thailand where concurrently collected water buffalo samples were examined and we found only Anopheles spp. samples were P. bubalis positive. Molecular identification of anopheline mosquito species was conducted by sequencing of the PCR products targeting cytochrome c oxidase subunit 1 (cox1), cytochrome c oxidase subunit 2 (cox2), and internal transcribed spacer 2 (ITS2) markers. We observed 5 distinct groups of anopheline mosquitoes: Barbirostris, Hyrcanus, Ludlowae, Funestus, and Jamesii groups. The Barbirostris group (Anopheles wejchoochotei or Anopheles campestris) and the Hyrcanus group (Anopheles peditaeniatus) were positive for P. bubalis. Thus, for the first time, our study implicated these anopheline mosquito species as probable vectors of P. bubalis in Thailand.

Identification of a biological form in the Anopheles stephensi laboratory colony using the odorant-binding protein 1 intron I sequence

Background

Anopheles stephensi Listen (1901) is a major vector of malaria in Asia and has recently been found in some regions of Africa. The An. stepehnsi species complex is suspected to have three sibling species: type, intermediate, and mysorensis, each with its own vector competence to the malaria parasite and ecology. To identify the members of the species complex in our An. stephensi insectary colony, we used the morphological features of eggs and genetic markers such as AnsteObp1 (Anopheles stephensi odorant binding protein 1), mitochondrial oxidases subunit 1 and 2 (COI and COII), and nuclear internal transcribed spacer 2 locus (ITS2).

Methods

Eggs were collected from individual mosquitoes (n = 50) and counted for the number of ridges under stereomicroscope. Genomic DNA was extracted from female mosquitoes. After the amplification of partial fragments of AnsteObp1, COI, COII and ITS2 genes, the PCR products were purified and sequenced. Phylogenetic analysis was performed after aligning query sequences against the submitted sequences in GenBank using MEGA 7.

Results

The range of ridges number on each egg float was 12–13 that corresponds to the mysorensis form of An. stephensi. The generated COI, COII and ITS2 sequences showed 100%, 99.46% and 99.29% similarity with the sequences deposited for Chinese, Indian and Iranian strains of An. stephensi, respectively. All the generated AnsteObp1 intron I region sequences matched 100% with the sequences deposited for An. stephensi sibling species C (mysorensis form) from Iran and Afghanistan.

Conclusions

This manuscript precisely describes the morphological and molecular details of the ‘var mysorensis’ form of An. stephensi that could be exploited in elucidating its classification as well as in differentiation from other biotypes of the same or other anopheline species. Based on our findings, we recommend AnsteObp1 as a robust genetic marker for rapid and accurate discrimination (taxonomic identification) of the An. stephensi species complex, rather than the COI, COII, and ITS2 marker, which could only be utilized for interspecies (Anopheles) differentiation.

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 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.

Phylogeny of Anopheles (Kerteszia) (Diptera: Culicidae) Using Mitochondrial Genes

Identification of mosquito species is necessary for determining the entomological components of malaria transmission, but it can be difficult in morphologically similar species. DNA sequences are largely used as an additional tool for species recognition, including those that belong to species complexes. Kerteszia mosquitoes are vectors of human and simian malaria in the Neotropical Region, but there are few DNA sequences of Kerteszia species in public databases. In order to provide relevant information about diversity and improve knowledge in taxonomy of Kerteszia species in Peru, we sequenced part of the mitochondrial genome, including the cytochrome c oxidase I (COI) barcode region. Phylogenetic analyses structured all species of mosquitoes collected in Peru into a single clade, separate from the Brazilian species. The Peruvian clade was composed of two lineages, encompassing sequences from Anopheles (Kerteszia) boliviensis and Anopheles (Kerteszia) pholidotus. An. pholidotus sequences were recorded for the first time in Peru, whereas An. boliviensis sequences were for the first time published in the GenBank database. Sequences generated from specimens morphologically identified as Anopheles (Kerteszia) cruzii clustered into three separate clades according to the collection localities of Serra do Mar, Serra da Mantiqueira, and Serra da Cantareira, confirming An. cruzii as a species complex, composed of at least three putative species.

Description of Aedes (Hulecoeteomyia) bhutanensis n. sp. (Diptera: Culicidae) from Bhutan

A new species of the subgenus Hulecoeteomyia of the genus Aedes reared from pupae found with two larvae in a small pool on a log in mountainous forest of Bhutan is diagnosed and formally named Aedes bhutanensis Somboon & Harbach, n. sp. The specific status of the species is supported by differential morphological characters of females, larvae and pupae, and sequences of the mitochondrial COI gene. Dichotomous keys are provided for identification of the adult females and larvae of the species of the subgenus.

Effects of indoor residual spraying and outdoor larval control on Anopheles coluzzii from São Tomé and Príncipe, two islands with pre-eliminated malaria

Background

Vector control is a key component of malaria prevention. Two major vector control strategies have been implemented in São Tomé and Príncipe (STP), indoor residual spraying (IRS) and outdoor larval control using Bacillus thuringiensis israelensis (Bti). This study evaluated post-intervention effects of control strategies on vector population density, composition, and knockdown resistance mutation, and their implications for malaria epidemiology in STP.

Methods

Mosquitoes were collected by indoor and outdoor human landing catches and mosquito light traps in seven districts. Mosquito density was calculated by numbers of captured adult mosquitoes/house/working hour. Mitochondrial cytochrome c oxidase subunit I (COI) was PCR amplified and sequenced to understand the spatial–temporal population composition of malaria vector in STP. Knockdown resistance L1014F mutation was detected using allele-specific PCR. To estimate the malaria transmission risks in STP, a negative binomial regression model was constructed. The response variable was monthly incidence, and the explanatory variables were area, rainfall, entomological inoculation rate (EIR), and kdr mutation frequency.

Results

Malaria vector in STP is exophilic Anopheles coluzzii with significant population differentiation between Príncipe and São Tomé (mean F_ST = 0.16, p < 0.001). Both vector genetic diversity and knockdown resistance mutation were relatively low in Príncipe (mean of kdr frequency = 15.82%) compared to São Tomé (mean of kdr frequency = 44.77%). Annual malaria incidence rate in STP had been rapidly controlled from 37 to 2.1% by three rounds of country-wide IRS from 2004 to 2007. Long-term application of Bti since 2007 kept the mosquito density under 10 mosquitoes/house/hr/month, and malaria incidence rate under 5% after 2008, except for a rising that occurred in 2012 (incidence rate = 6.9%). Risk factors of area (São Tomé compared to Príncipe), rainfall, outdoor EIR, and kdr mutation frequency could significantly increase malaria incidence by 9.33–11.50, 1.25, 1.07, and 1.06 fold, respectively.

Conclusions

Indoor residual spraying could rapidly decrease Anopheles density and malaria incidence in STP. Outdoor larval control using Bti is a sustainable approach for controlling local vector with exophilic feature and insecticide resistance problem. Vector control interventions should be intensified especially at the north-eastern part of São Tomé to minimize impacts of outbreaks.

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.

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

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

Widespread plant specialization in the polyphagous planthopper Hyalesthes obsoletus (Cixiidae), a major vector of stolbur phytoplasma: Evidence of cryptic speciation

The stolbur phytoplasma vector Hyalesthes obsoletus is generally considered as a polyphagous species associated with numerous wild and cultivated plants. However, recent research in southeastern Europe, the distribution centre of H. obsoletus and the area of most stolbur-inflicted crop diseases, points toward specific host-plant associations of the vector, indicating specific vector-based transmission routes. Here, we study the specificity of populations associated with four host-plants using mitochondrial and nuclear genetic markers, and we evaluate the evolution of host-shifts in H. obsoletus. Host-plant use was confirmed for Convolvulus arvensis, Urtica dioica, Vitex agnus-castus and Crepis foetida. Mitochondrial genetic analysis showed sympatric occurrence of three phylogenetic lineages that were ecologically delineated by host-plant preference, but were morphologically inseparable. Nuclear data supported the existence of three genetic groups (Evanno’s ΔK(3) = 803.72) with average genetic membership probabilities > 90%. While populations associated with C. arvensis and U. dioica form a homogenous group, populations affiliated with V. agnus-castus and C. foetida constitute two independent plant-associated lineages. The geographical signal permeating the surveyed populations indicated complex diversification processes associated with host-plant selection and likely derived from post-glacial refugia in the eastern Mediterranean. This study provides evidence for cryptic species diversification within H. obsoletus sensu lato: i) consistent mitochondrial differentiation (1.1–1.5%) among host-associated populations in syntopy and in geographically distant areas, ii) nuclear genetic variance supporting mitochondrial data, and iii) average mitochondrial genetic distances among host-associated meta-populations are comparable to the most closely related, morphologically distinguishable species, i.e., Hyalesthes thracicus (2.1–3.3%).