Molecular comparison of topotypic specimens confirms Anopheles (Nyssorhynchus) dunhami Causey (Diptera: Culicidae) in the Colombian Amazon

The presence of Anopheles (Nyssorhynchus) dunhami Causey in Colombia (Department of Amazonas) is confirmed for the first time through direct comparison of mtDNA cytochrome c oxidase I (COI) barcodes and nuclear rDNA second internal transcribed spacer (ITS2) sequences with topotypic specimens of An. dunhami from Tefé, Brazil. An. dunhami was identified through retrospective correlation of DNA sequences following misidentification as Anopheles nuneztovari s.l. using available morphological keys for Colombian mosquitoes. That An. dunhami occurs in Colombia and also possibly throughout the Amazon Basin, is of importance to vector control programs, as this non-vector species is morphologically similar to known malaria vectors including An. nuneztovari, Anopheles oswaldoi and Anopheles trinkae. Species identification of An. dunhami and differentiation from these closely related species are highly robust using either DNA ITS2 sequences or COI DNA barcode. DNA methods are advocated for future differentiation of these often sympatric taxa in South America.     

Molecular taxonomy of members of the Anopheles hyrcanus group from Thailand and Indonesia

During studies of malaria vectors in Indonesia and Thailand, several specimens identified by field staff as members of the Anopheles barbirostris group (Diptera: Culicidae) were found to belong to the Anopheles hyrcanus group, as shown by marked differences in the size of the nuclear rDNA second internal transcribed spacer (ITS2) between the barbirostris (~1500 bp) and hyrcanus (~600 bp) groups. Identification of the species concerned required a more detailed study of ITS2 sequences and subunit I of the mitochondrial DNA cytochrome oxidase gene (COI). A phylogenetic analysis, based on Bayesian methods, revealed that the hyrcanus group specimens comprised five distinct clades, two of which corresponded with known species, Anopheles peditaeniatus and Anopheles sinensis. The remaining specimens formed three additional clades, for which there are no similar sequences in GenBank and which cannot be linked to previously described species. The misidentification of hyrcanus group species has important implications for malaria vector control; more comprehensive studies employing gene sequences are required to clarify the number of species in the group, their distribution and vector status.     

Intragenomic rDNA ITS2 variation in the neotropical Anopheles (Nyssorhynchus) albitarsis complex (Diptera: Culicidae)

We cloned and sequenced the rDNA internal transcribed spacer 2 (ITS2) of 4 species belonging to the neotropical Anopheles (Nyssorhynchus) albitarsis complex, that is, A. albitarsis; A. albitarsis B; Anopheles marajoara, a proven malaria vector; and Anopheles deaneorum, a suspected vector. Even though the ITS2 sequences of these species were very similar (< or =1.17% divergence), we found differences suitable for species identification and intragenomic variation of possible consequence in phylogenetic reconstruction. Variation came from 2 microsatellite regions and a number of indels and base substitutions. The existence of partially correlated subsets of clones in A. albitarsis is hypothesized either to be separate rDNA loci or to be semi-independently evolving portions of a single rDNA locus. No differences were found between males and females, suggesting that similar rDNA arrays exist on both the X and Y chromosomes. In addition, highly variant clones, possibly pseudogenes, were found in A. marajoara from Venezuela.     

Molecular evidence for a single taxon, Anopheles nuneztovari s.l., from two endemic malaria regions in Colombia

To elucidate the Anopheles nuneztovari s.l. taxonomic status at a microgeographic level in four malaria endemic localities from Antioquia and Córdoba, Colombia, fragments of the cytochrome oxidase subunit I (COI) and the white gene were used. The COI analysis showed low genetic differentiation with fixation index (F(ST)) levels between -0.02-0.137 and Nm values between 3-∞, indicating the presence of high gene flow among An. nuneztovari s.l. populations from the four localities. The COI network showed a single most common haplotype, type 1 (n = 55), present in all localities, as the likely ancestral haplotype. Analysis of the white gene showed that An. nuneztovari s.l. populations from both departments grouped with haplotypes 19 and 20, which are part of lineage 3 reported previously. The results of the present study suggest that An. nuneztovari s.l. is a single taxon in the area of the present study.     

Sequence analysis of the ribosomal DNA internal transcribed spacer 2 from populations of Anopheles nuneztovari (Diptera: Culicidae)

Sequence variation of the ribosomal DNA internal transcribed spacer 2 (ITS2) was examined for populations of the malaria vector Anopheles nuneztovari collected in Colombia, Venezuela, Bolivia, Suriname, and Brazil. Mosquitoes from Colombia and Venezuela had identical ITS2 sequences and were distinguished from sequences in other populations by three insertion/deletion events (indels) and by one transversion. The length of the ITS2 was 363-369 bp, and it had a G+C content of 55.3%- 55.7%. Variation in the length of the ITS2 between and within populations was due to indels in simple repeats. ITS2 consensus sequences were similar or identical for samples from the following three groups: (1) Colombia, Bolivia, and Venezuela; (2) Suriname and northern Brazil; and (3) eastern and central Brazil. The presence of two different consensus sequences from a single location near Manaus, Brazil, suggests that populations from eastern Brazil and those from Suriname converge in this region of the Amazon Basin. These data show that putative cryptic species of An. nuneztovari are distinguished by very minor differences in DNA sequence of the ITS2 region.      

Molecular Phylogeny of Neotropical Anopheles (Nyssorhynchus) Albitarsis Species Complex (Diptera: Culicidae)

A phylogeny was reconstructed for four species belonging to the Neotropical Anopheles (Nyssorhynchus) albitarsis complex using partial sequences from the mitochondrial cytochrome oxidase I (COI) and NADH dehydrogenase 4 (ND4) genes and the ribosomal DNA ITS2 and D2 expansion region of the 28S subunit. The basis for initial characterization of each member of the complex was by correlated random amplification of polymorphic DNA-polymerase chain reaction (RAPD-PCR) markers. Analyses were carried out with and without an outgroup (An.(Nys.) argyritarsis Robineau-Desvoidy) by using maximum parsimony, maximum likelihood, and Bayesian methods. A total evidence approach without the outgroup, using separate models for "fast" (COI and ND4 position 3) and "slow" (rDNA ITS2 and D2, and COI and ND4 position 1) partitions, gave the best supported topology, showing close relationships of An. albitarsis Lynch-Arribálzaga to An. albitarsis B and An. marajoara Galv?o & Damasceno to An. deaneorum Rosa-Freitas. Analyses with the outgroup included showed poorer support, possibly because of a long branch attraction effect caused by a divergent outgroup, which caused one of the An. marajoara specimens to cluster with An. deaneorum in some analyses. The relationship of the above-mentioned result to a separately proposed hypothesis suggesting a fifth species in the complex is discussed.     

Cryptic Species in the Anopheles (Nyssorhynchus) albitarsis (Diptera: Culicidae) Complex: Incongruence Between Random Amplified Polymorphic DNA-Polymerase Chain Reaction Identification and Analysis of Mitochondrial DNA COI Gene Sequences

Random amplified polymorphic DNA (RAPD) diagnostic bands are one tool used to differentiate cryptic mosquito species in the Anopheles albitarsis Complex. Monophyly of four species (A. albitarsis Lynch-Arribálzaga, A. albitarsis B, A. deaneorum Rosa-Freitas, and A. marajoara Galv?o & Damasceno) currently identified with the RAPD technique was assessed using sequences of the cytochrome oxidase I (COI) mitochondrial DNA (mtDNA) gene. Maximum parsimony, maximum likelihood, and Bayesian analyses support monophyly for A. albitarsis s.s., A. albitarsis B, and A. deaneorum. Anopheles marajoara, as identified by RAPD banding patterns, was either polyphyletic or paraphyletic in all phylogenetic analyses. The phylogenetic pattern and within-species genetic distances observed in A. marajoara suggest the existence of a previously unidentified species (species E) in northern Brazil and Venezuela. Diagnostic RAPD bands were unable to distinguish between A. marajoara and species E, probably because of the low number of correlated bands used to identify species and weaknesses of the RAPD technique, in particular, violations of the untested assumption of homology of comigrating bands. A. marajoara (even without species E) is paraphyletic with respect to A. deaneorum; if A. deaneorum is a separate species from A. marajoara, then A. marajoara may consist of two or more species in Amazonian Brazil. Based on mtDNA COI sequences, there are at least four phylogenetic species within the Albitarsis Complex: A. albitarsis s.s., A. albitarsis B, A. marajoara, and species E; the species status of A. deaneorum is ambiguous.     

Molecular systematics of the Philippine malaria vector Anopheles flavirostris

Allozyme and molecular sequence data from the malaria vector Anopheles flavirostris (Ludlow) (Diptera: Culicidae) were analysed from 34 sites throughout the Philippines, including the type locality, to test the hypothesis that this taxon is a single panmictic species. A finer-scaled allozyme study, of mainly Luzon samples, revealed no fixed genetic differences in sympatric sites and only low levels of variation. We obtained data from partial sequences for the internal transcribed spacer 2 (ITS2) (483 bp), the third domain (D3) (330 bp) of the 28S ribosomal DNA subunit and cytochrome c oxidase subunit I (COI) of mitochondrial DNA (261 bp). No sequence variation was observed for ITS2, only a one base pair difference was observed between Philippine and Indonesian D3 sequences and An. flavirostris sequences were unique, confirming their diagnostic value for this taxon. Sixteen COI haplotypes were identified, giving 25 parsimony informative sites. Neighbour-Joining, Maximum Parsimony, Maximum Likelihood and Bayesian phylogenetic analysis of COI sequences for An. flavirostris and outgroup taxa revealed strong branch support for the monophyly of An. flavirostris, thus confirming that Philippine populations of this taxon comprise a single separate species within the Minimus Subgroup of the Funestus Group. Variation in the behaviour of An. flavirostris is likely to be intraspecific rather than interspecific in origin.     

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.     

Evidence for a conspecific relationship between two morphologically and cytologically different forms of Korean Anopheles pullus mosquito

Despite the medical importance of anopheline mosquitoes as vectors of Korean vivax malaria, differentiation between Korean anopheline mosquitoes by traditional morphological taxonomic criteria is difficult. An. yatsushiroensis is the second most common Anopheles mosquito species in Korea and a possible vector of Korean vivax malaria together with An. sinensis, the predominant anopheline species. Recently, An. yatsushiroensis has been declared a synonym of An. pullus, based on comparisons of egg morphology and adult progeny, although they differ in ecology and morphology. To verify the species status of these two ambiguous forms, we established isofemale lines of Korean An. pullus and An. yatsushiroensis (An. pullus form yatsushiroensis) mosquitoes and investigated their genetic relationship by metaphase karyotype analysis, comparing the DNA sequences of rDNA internal transcribed spacer 2 (ITS2) and mitochondrial cytochrome c oxidase subunits I (COI) and II (COII), and by hybridization experiments. Two isofemale lines had differently shaped X and Y chromosomes. However reciprocal crosses between them yielded viable progeny with completely synaptic salivary gland polytene chromosomes. DNA analyses also strongly supported their conspecificity. The two strains also showed great sequence similarity in the ITS2, COI and COII regions (variation rate = 0.0 to 0.8%). Based on these findings, we conclude that the two forms, though differing distinctly in morphological, cytological and ecological traits, remain interfertile.     

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.     

我国赫坎按蚊复合体成员种的rDNA-ITS2区序列差异及系统发育分析

测定了我国赫坎按蚊复合体 9成员种的核糖体DNA第二内转录间隔区 (rDNA ITS2 )序列 ,根据序列差异分析各蚊种间的系统发育关系。结果显示 :( 1 )ITS2区序列最长的是中华按蚊 ( 4 6 8bp) ,最短的是克劳按蚊和赫坎按蚊 ( 4 36bp) ;GC含量为 4 4 9%～ 4 6 8% ;( 2 )发现该复合体 4成员种的ITS2区序列存在种内个体间差异 ,幅度为 0～ 3 8% ,明显小于种间差异 ;( 3)将各蚊种的ITS2区序列进行同源排序比较 ,发现其变异大多是简单重复单元的拷贝数不同 ;种间差异性最大的是克劳按蚊与嗜人按蚊( 32 3% ) ,最小的是贵阳按蚊与凉山按蚊 ( 9 0 % )平均差异率为 2 2 3% ;( 4 )根据ITS2区序列特征 ,用 3种方法构建的树状图拟合一致。以上结果表明赫坎按蚊复合体各成员种rDNA ITS2序列在种内非常保守 ,以种间序列差异分析为基础的分子鉴别技术是甄别蚊种分类地位混淆和错误的有效方法。     

Confirmation of Anopheles (Anopheles) calderoni Wilkerson, 1991 (Diptera: Culicidae) in Colombia and Ecuador through molecular and morphological correlation with topotypic material

The morphologically similar taxa Anopheles calderoni, Anopheles punctimacula, Anopheles malefactor and Anopheles guarao are commonly misidentified. Isofamilies collected in Valle de Cauca, Colombia, showed morphological characters most similar to An. calderoni, a species which has never previously been reported in Colombia. Although discontinuity of the postsubcostal pale spots on the costa (C) and first radial (R1) wing veins is purportedly diagnostic for An. calderoni, the degree of overlap of the distal postsubcostal spot on C and R1 were variable in Colombian specimens (0.003-0.024). In addition, in 98.2% of larvae, seta 1-X was located off the saddle and seta 3-C had 4-7 branches in 86.7% of specimens examined. Correlation of DNA sequences of the second internal transcribed spacer and mtDNA cytochrome c oxidase subunit I gene (COI) barcodes (658 bp of the COI gene) generated from Colombian progeny material and wild-caught mosquitoes from Ecuador with those from the Peruvian type series of An. calderoni confirmed new country records. DNA barcodes generated for the closely related taxa, An. malefactor and An. punctimacula are also presented for the first time. Examination of museum specimens at the University of the Valle, Colombia, revealed the presence of An. calderoni in inland localities across Colombia and at elevations up to 1113 m.     

Population analysis using the nuclear white gene detects Pliocene/Pleistocene lineage divergence within Anopheles nuneztovari in South America

Anopheles (Nyssorhynchus) nuneztovari Gabaldón (Diptera: Culicidae), a locally important malaria vector in some regions of South America, has been hypothesized to consist of at least two cryptic incipient species. We investigated its phylogeographic structure in several South American localities to determine the number of lineages and levels of divergence using the nuclear white gene, a marker that detected two recently diverged genotypes in the primary Neotropical malaria vector Anopheles darlingi Root. In An. nuneztovari, five distinct lineages (1-5) were elucidated: (1) populations from northeastern and central Amazonia; (2) populations from Venezuela east and west of the Andes; (3) populations from Colombia and Venezuela west of the Andes; (4) southeastern and western Amazonian Brazil populations, and (5) southeastern and western Amazonian Brazil and Bolivian populations. There was a large amount of genetic differentiation among these lineages. The deepest and earliest divergence was found between lineage 3 and lineages 1, 2 and 4, which probably accounts for the detection of lineage 3 in some earlier studies. The multiple lineages within Amazonia are partially congruent with previous mtDNA and ITS2 data, but were undetected in many earlier studies, probably because of their recent (Pleistocene) divergence and the differential mutation rates of the markers. The estimates for the five lineages, interpreted as recently evolved or incipient species, date to the Pleistocene and Pliocene. We hypothesize that the diversification in An. nuneztovari is the result of an interaction between the Miocene/Pliocene marine incursion and Pleistocene climatic changes leading to refugial isolation. The identification of cryptic lineages in An. nuneztovari could have a significant impact on local vector control measures.     

Analysis of the evolutionary forces shaping mitochondrial genomes of a Neotropical malaria vector complex

Many vectors of human malaria belong to complexes of morphologically indistinguishable cryptic species. Here we report the analysis of the newly sequenced complete mitochondrial DNA molecules from six recognized or putative species of one such group, the Neotropical Anopheles albitarsis complex. The molecular evolution of these genomes had been driven by purifying selection, particularly strongly acting on the RNA genes. Directional mutation pressure associated with the strand-asynchronous asymmetric mtDNA replication mechanism may have shaped a pronounced DNA strand asymmetry in the nucleotide composition in these and other Anopheles species. The distribution of sequence polymorphism, coupled with the conflicting phylogenetic trees inferred from the mitochondrial DNA and from the published white gene fragment sequences, indicates that the evolution of the complex may have involved ancient mtDNA introgression. Six protein coding genes (nad5, nad4, cox3, atp6, cox1 and nad2) have high levels of sequence divergence and are likely informative for population genetics studies. Finally, the extent of the mitochondrial DNA variation within the complex supports the notion that the complex consists of a larger number of species than until recently believed.     

Presence of Anopheles culicifacies B in Cambodia established by the PCR-RFLP assay developed for the identification of Anopheles minimus species A and C and four related species

Abstract A polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) assay developed for identification of five species of the Anopheles minimus Theobald group and a related mosquito species of the Myzomyia Series (Diptera: Culicidae) was applied to morphologically identified adult female specimens collected in Ratanakiri Province, north‐eastern Cambodia. In addition to finding An. aconitus Dönitz, An. minimus species A and An. pampanai Büttiker & Beales, some specimens showed a new restriction banding pattern. Siblings of specimens that exhibited this new PCR‐RFLP pattern were morphologically identified as An. culicifacies James sensu lato. Based on nucleotide sequences of the ribonuclear DNA internal transcribed spacer 2 region (ITS2) and the mitochondrial cytochrome oxidase I gene (COI), these specimens were recognized as An. culicifacies species B (sensu Green & Miles, 1980 ), the first confirmed record of the An. culicifacies complex from Cambodia. This study shows that the PCR‐RFLP assay can detect species not included in the initial set‐up and is capable of identifying at least seven species of the Myzomyia Series, allowing better definition of those malaria vector and non‐vector anophelines in South‐east Asia.     

Infectivity of two strains of Plasmodium vivax to Anopheles albitarsis mosquitoes from Colombia

Anopheles albitarsis obtained from Villavicencio, Colombia, were colonized in the laboratory using force-mating techniques. Laboratory reared mosquitoes were allowed to feed on Aotus monkeys infected with the Salvador II or the Rio Meta strains of Plasmodium vivax from El Salvador and Colombia, respectively. In comparison with other species, the An. albitarsis were less susceptible than Anopheles freeborni, Anopheles culicifacies and strains of Anopheles albimanus from El Salvador, Panama and Colombia and more susceptible than a strain of An. albimanus from Haiti.     

DNA Barcodes indicate members of the Anopheles fluviatilis (Diptera: Culicidae) species complex to be conspecific in India

Anopheles fluviatilis, a major vector of malaria in India has been described as a complex of three sibling species members, named as S, T and U, based on variations in chromosomal inversions. Also, ribosomal DNA markers (repetitive Internal Transcribed Spacer 2 (ITS2) and 28S D3 region) were described to differentiate these three sibling species members. However, controversies prevail on the genetic isolation status of these cryptic species. Hence, we evaluated this taxonomic incongruence employing DNA barcoding, the well established methodology for species identification, using 60 An. fluviatilis sensu lato specimens, collected from two malaria endemic eastern states of India. These specimens were also subjected to sibling species characterization by ITS2 and D3 DNA markers. The former marker identified 31 specimens among these as An. fluviatilis S and 21 as An. fluviatilis T. Eight specimens amplified DNA fragments specific for both S and T. The D3 marker characterized 39 specimens belonging to species S and 21 to species T. Neither marker identified species U. Neighbor Joining analysis of mitochondrial cytochrome c oxidase gene 1 sequences (the DNA barcode) categorized all the 60 specimens into a single operational taxonomic unit, their Kimura 2 parameter (K2P) genetic variability being only 0.8%. The genetic differentiation (F_ST) and gene flow (N_m) estimates were 0.00799 and 62.07, respectively, indicating these two ‘species’ (S & T) as genetically con‐specific intermixing populations with negligible genetic differentiation. Earlier investigations have refuted the existence of species U. Also, this study demonstrated that An. fluviatilis and the closely related An. minimus could be taxonomically differentiated by the DNA Barcode approach (K2P = 5.0%).     

High levels of population structure caused by habitat islands in the malarial vector Anopheles scanloni

The genetic structure of four populations of the malarial vector Anopheles scanloni in Thailand was studied using mitochondrial DNA sequences. Four highly divergent lineages were observed, all with signals of population expansion. Since An. scanloni is restricted to 'islands' of limestone karst habitat, this suggests there is a metapopulation-type dynamic in this species, with restricted gene flow, extinctions and drift all contributing to lineage divergence. Historical environmental change and marine transgressions may also have contributed to population extinction, expansion and divergence. Although there is some current gene flow inferred between nearby populations, it is extremely restricted between the northern and southern populations, which also differed by one fixed polymorphism at the ITS2 rDNA locus. Crossing experiments showed no post-mating barriers existing between the north and the south, but the lack of gene flow between these populations could ultimately result in speciation and has implications for malaria control strategies.     

Phylogeny of anopheline (Diptera: Culicidae) species in southern Africa,based on nuclear and mitochondrial genes

A phylogeny of anthropophilic and zoophilic anopheline mosquito species was constructed, using the nuclear internal transcribed spacer 2 (ITS2) and mitochondrial cytochrome oxidase subunit I (COI) genes. The ITS2 alignment, typically difficult due to its noncoding nature and large size variations, was aided by using predicted secondary structure, making this phylogenetically useful gene more amenable to investigation. This phylogeny is unique in explicitly including zoophilic, non‐vector anopheline species in order to illustrate their relationships to malaria vectors. Two new, cryptic species, Anopheles funestus‐like and Anopheles rivulorum‐like, were found to be present in Zambia for the first time. Sequences from the D3 region of the 28S rDNA suggest that the Zambian An. funestus‐like may be a hybrid or geographical variant of An. funestus‐like, previously reported in Malawi. This is the first report of An. rivulorum‐like sympatric with An. rivulorum (Leeson), suggesting that these are separate species rather than geographic variants.