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

COI DNA barcoding to differentiate Haemagogus janthinomys and Haemagogus capricornii (Diptera: Culicidae) mosquitoes

The genus Haemagogus (Diptera: Culicidae) includes species that are important vectors of pathogens such as the yellow fever virus. The accurate identification of these species is essential for the control of zoonoses. Females of Hg. capricornii and Hg. janthinomys are morphologically indistinguishable, which makes the use of alternative identification techniques desirable. This study aimed to obtain sequences of the mitochondrial cytochrome c oxidase I (COI) gene, in the region widely used for DNA barcoding, of Haemagogus specimens from the state of São Paulo, Brazil, to evaluate the effectiveness of these sequences in the molecular identification of the species. A total of 37 female and 2 male mosquitoes were collected in various locations in the state of São Paulo, using methods such as hand-nets, Shannon traps, CDC light traps with CO2 bait and Nasci aspirators. The sequences of a 710 bp fragment of the COI gene were amplified by PCR and sequenced. A phylogenetic tree reconstruction was conducted using the Bayesian approach implemented in MrBayes v3.2.2, providing support values for taxa where genetic clusters may indicate the presence of new or cryptic species. We obtained 39 COI sequences representing three species: Haemagogus capricornii, Haemagogus leucocelaenus, and Haemagogus janthinomys. Bayesian analysis of the sequences produced clades that corroborate the morphological identification of the species. The separation of Hg. capricornii and Hg. janthinomys received 100 % statistical support and the Hg. capricornii was very well supported (91 %). The two sequences from male specimens, morphologically identified as Hg. capricornii, were grouped in the same clade, a sister clade of Hg. janthinomys. It is important to highlight that the Hg. janthinomys were positioned in several subclades, showing a polymorphism of this species within the state, a situation not observed for Hg. capricornii. For the first time, sequences of the mtCOI gene from Hg. capricornii were obtained and related to morphologically identified specimens. COI sequences proved effective in the molecular identification of Haemagogus species. This study contributes to the expansion of the GenBank database, providing the first sequences of Hg. capricornii and new sequences for Hg. janthinomys and Hg. leucocelaenus.

High Levels of Diversity in Anopheles Subgenus Kerteszia Revealed by Species Delimitation Analyses

The Anopheles subgenus Kerteszia is a poorly understood group of mosquitoes that includes several species of medical importance. Although there are currently twelve recognized species in the subgenus, previous studies have shown that this is likely to be an underestimate of species diversity. Here, we undertake a baseline study of species delimitation using the barcode region of the mtDNA COI gene to explore species diversity among a geographically and taxonomically diverse range of Kerteszia specimens. Beginning with 10 of 12 morphologically identified Kerteszia species spanning eight countries, species delimitation analyses indicated a high degree of cryptic diversity. Overall, our analyses found support for at least 28 species clusters within the subgenus Kerteszia. The most diverse taxon was Anopheles neivai, a known malaria vector, with eight species clusters. Five other species taxa showed strong signatures of species complex structure, among them Anopheles bellator, which is also considered a malaria vector. There was some evidence for species structure within An. homunculus, although the results were equivocal across delimitation analyses. The current study, therefore, suggests that species diversity within the subgenus Kerteszia has been grossly underestimated. Further work will be required to build on this molecular characterization of species diversity and will rely on genomic level approaches and additional morphological data to test these species hypotheses.

DNA Barcoding of Morphologically Characterized Mosquitoes Belonging to the Genus Mansonia from the Atlantic Forest and Brazilian Savanna

The identification of mosquito species is necessary for determining the entomological components of disease transmission. However, identification can be difficult in species that are morphologically similar. The cytochrome c oxidase subunit I (COI) DNA barcode region is considered a valuable and reliable diagnostic tool for mosquito species recognition, including those that belong to species complexes. Mansonia mosquitoes are found in forests near swampy areas. They are nocturnal and are highly attracted to light. Hematophagous adult females exhibit aggressive biting behavior and can become infected with and transmit pathogens during their feeding, including some epizootic viruses and avian malaria. In Brazil, twelve Mansonia species have been reported. In a recent study from the São Paulo Zoo in Brazil, three morphologically distinct species were collected and identified, namely: Mansonia (Mansonia) indubitans, Ma. (Man.) pseudotitillans and Ma. (Man.) titillans. However, confirmation of these species by molecular identification was unsuccessful due to a lack of COI sequences in the GenBank database. Thus, this research aimed to describe the COI DNA barcode sequences of some morphologically characterized Mansonia (Man.) species from Brazil and to determine their utility in delimiting species collected from the Atlantic Forest and Brazilian Savanna. Accordingly, we provide tools for the genetic identification of species that play a significant role in pathogen transmission in wildlife and potentially humans. We show that the delimitation of Mansonia species via five different approaches based on COI DNA sequences (BI, NJ, ASAP, bPTP and GMYC) yield basically the same groups identified by traditional taxonomy, and we provide the identification of specimens that were previously identified only up to the subgenus level. We also provide COI sequences from two Mansonia species that were not previously available in sequence databases, Ma. wilsoni and Ma. pseudotitillans, and thus contribute to the ongoing global effort to standardize DNA barcoding as a molecular means of species identification.

Evidence of Elevational Speciation in Kerteszia cruzii (Diptera: Culicidae) in the Ribeira Valley, São Paulo, Brazil

Kerteszia cruzii [former Anopheles (Kerteszia) cruzii] is a bromeliad mosquito widespread in the Brazilian Atlantic rainforest. In South-eastern Brazil, it plays an important role in malaria transmission because it was infected with at least four Plasmodium species. There is robust evidence that Ke. cruzii is a species complex. We used single nucleotide polymorphisms (SNPs) from a nextRAD sequence (nextera-tagmented, reductively amplified DNA) to investigate the genetic structure of Ke. cruzii in the Ribeira Valley, South-eastern Brazil. Furthermore, we verified whether the genetic structure was associated with forest cover, elevation, slope, and vegetation physiognomy. Our results showed two distinct lineages in the studied region associated with elevation and isolation by distance. The first lineage included samples from coastal localities and the second comprised specimens from inland or mountain sites. At one sampling locality (Esteiro do Morro in Cananéia municipality), both lineages are sympatric. These results are in accordance with previously published data that showed elevated stratification in Ke. cruzii. However, Fst values did not indicate the existence of cryptic or sister species in Ke. cruzii in this region, we concluded that elevational speciation probably occurs, and we hypothesized that differences in population structure found might be associated with the distribution of bromeliad species.

Genetic Diversity and Molecular Characterization of Mosquitoes (Diptera: Culicidae) In North-Central Nigeria Using Ribosomal DNA ITS2 and Mitochondrial 16S-DNA Sequences

Mosquitoes are vectors of various life-threatening diseases like malaria, yellow fever, dengue fever etc. Their close proximity to human habitations allows ease for disease transmission. They have been identified by key morphological tools like their wings, legs, bristles etc. but closely related species are difficult to identify based on morphology. Molecular tools have, therefore, been employed to help with the more accurate identification. This study was aimed at identifying and characterizing different mosquito species in five different states in North-Central Nigeria using internal transcribed spacer 2 (ITS2) and mitochondrial 16S rDNA regions. Mosquito larvae were collected from stagnant water in breeding places at each collection site in North-central Nigeria. Morphological identification was carried out using standard keys. DNA extraction was performed using EZNA extraction kit. PCR amplification of ribosomal ITS2 and mitochondrial 16S-rDNA gene regions were carried out. The PCR amplicons were sequenced using primers initially used for the PCR. Sequence data were aligned in MEGA 6.0 using ClustalW multiple alignment feature and then compared with GenBank databases for similarity.  Phylogenetic analysis of DNA sequences from the ITS2 region was able to distinguish two mosquito subfamilies; Anophelinae and Culicinae as well as differentiate between and amongst Culex and Aedes species. However, it was unable to effectively distinguish between the two different species of Anopheles sequenced. Mitochondrial 16S rRNA marker was also able to distinguish the two mosquito subfamilies. It efficiently identified and differentiated Culex, Aedes and Anopheles mosquito species sequenced in this study. This study concludes that heterogeneity among Nigerian populations of Anopheles mosquitoes of may likely impact malaria vector control programs. We recommend the combination of nuclear and mitochondrial markers for effective and reliable phylogenetic study and determination of evolutionary relationship among mosquito species.

Wing geometric morphometrics for identification of mosquito species (Diptera: Culicidae) of neglected epidemiological importance

Culicidae comprises more than 3500 species, some of which are responsible for the spread of various human diseases, causing millions of deaths worldwide. Correct identification of these species is essential for the development of surveillance and control strategies. The most common method of mosquito identification is based on specific traits of the external morphology of species. However, identification of mosquitoes by morphological characters can be inaccurate or even unfeasible if the specimen is damaged or there is a lack of distinguishing features, as in the case of cryptic species complexes. Wing geometric morphometrics is a reliable, affordable tool for the identification of mosquito species, including sibling species. More importantly, it can be used in addition to both traditional morphologic identification methods as well as genetic approaches. Here, wing geometric morphometrics was used to identify sixteen mosquito species from eight genera: Aedes, Coquillettidia, Culex, Limatus, Mansonia, Psorophora, Runchomyia, and Wyeomyia. The 390 specimens used here were collected in São Paulo, Brazil using CDC traps, aspiration, and Shannon traps. Allometry was assessed by multivariate regression of the Procrustes coordinates on centroid size followed by canonical variate analysis and a pairwise cross-validated reclassification test. A Neighbor-Joining tree based on Mahalanobis distances was constructed with 1,000 bootstrap replicates using MorphoJ 1.02 and Past 2.17c. The canonical variate analysis of genera resulted in distinct clusters for Culex, Limatus, and Psorophora and partial overlapping between Aedes, Coquilettidia, and Mansonia, and between Runchomyia and Wyeomyia. Pairwise cross-validated reclassification tests indicated that genera were identified with an accuracy of at least 99% and subgenera with a mean accuracy of 96% and that in 160 of the 240 possible comparisons species were identified with an accuracy of 100%. Our results show that the eight genera in the study were correctly distinguished by wing shape, as were subgenera and most species, demonstrating that wing geometric morphometrics can be used for the identification of the mosquito species studied here.

Mitochondrial genome sequencing and phylogeny of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis (Diptera: Culicidae)

The genus Haemagogus (Diptera: Culicidae) comprises species of great epidemiological relevance, involved in transmission cycles of the Yellow fever virus and other arboviruses in South America. So far, only Haemagogus janthinomys has complete mitochondrial sequences available. Given the unavailability of information related to aspects of the evolutionary biology and molecular taxonomy of this genus, we report here, the first sequencing of the mitogenomes of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis. The mitogenomes showed an average length of 15,038 bp, average AT content of 79.3%, positive AT-skews, negative GC-skews, and comprised 37 functional subunits (13 PCGs, 22 tRNA, and 02 rRNA). The PCGs showed ATN as start codon, TAA as stop codon, and signs of purifying selection. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogenetic analyzes of Bayesian inference and Maximum Likelihood, based on concatenated sequences from all 13 PCGs, produced identical topologies and strongly supported the monophyletic relationship between the Haemagogus and Conopostegus subgenera, and corroborated with the known taxonomic classification of the evaluated taxa, based on external morphological aspects. The information produced on the mitogenomes of the Haemagogus species evaluated here may be useful in carrying out future taxonomic and evolutionary studies of the genus.