Wing geometric morphometrics and molecular assessment of members in the Albitarsis Complex from Colombia

Molecular differentiation analysis of ten putative species of Fannia (Diptera: Fanniidae) collected in carrion-baited traps from Colombia

The genus Fannia is the most representative of the Fannidae family of true flies with worldwide distribution. Some species are attracted to decomposing materials and live vertebrate animals, which makes them important in forensics, medical and veterinary fields. However, identifying Fannia species can be difficult due to the high similarity in the external morphology of females and limited descriptions and morphological keys. Herein, molecular markers could provide a complementary tool for species identification. However, molecular identification has still limited application since databases contain few data for neotropical species of Fannia. This study assessed the potential of two molecular markers, the COI-3' region and internal transcribed spacer 2 (ITS2), to differentiate 10 putative species of the genus Fannia from Colombia using distance-based and tree-based approaches. The partial ITS2 and/or COI-3' regions allowed molecular diagnosis of six species, while pairs of species Fannia colazorrensis + F. dodgei and F. laclara + F. aburrae are conflicting. Although these results might suggest that conflicting pair species are conspecific, consistent morphological differences between males do not support this hypothesis. The lack of differentiation at the nuclear and mitochondrial molecular markers for the conflicting species may be due to incomplete evolutionary lineage separation, hybridization, or introgression events. In addition, sexual selection on male morphological traits before species-specific differences in molecular markers emerge may partially explain the results. Our study provides a valuable dataset to identify and confirm some Fannia species molecularly. Further, they could be used to associate females and immature stages with their conspecifics as a baseline to deep into their biology, ecology, distribution and potential applications in forensic and medico-veterinary entomology.

Updating the bionomy and geographical distribution of Anopheles (Nyssorhynchus) albitarsis F: A vector of malaria parasites in northern South America

Anopheles albitarsis F is a putative species belonging to the Albitarsis Complex, recognized by rDNA, mtDNA, partial white gene, and microsatellites sequences. It has been reported from the island of Trinidad, Venezuela and Colombia, and incriminated as a vector of malaria parasites in the latter. This study examined mitochondrially encoded cytochrome c oxidase I (MT-CO1) sequences of An. albitarsis F from malaria-endemic areas in Colombia and Venezuela to understand its relations with other members of the Complex, revised and update the geographical distribution and bionomics of An. albitarsis F and explore hypotheses to explain its phylogenetic relationships and geographical expansion. Forty-five MT-CO1 sequences obtained in this study were analyzed to estimate genetic diversity and possible evolutionary relationships. Sequences generated 37 haplotypes clustered in a group where the genetic divergence of Venezuelan populations did not exceed 1.6% with respect to Colombian samples. Anopheles albitarsis F (π = 0.013) represented the most recent cluster located closer to An. albitarsis I (π = 0.009). Barcode gap was detected according to Albitarsis Complex lineages previously reported (threshold 0.014–0.021). Anopheles albitarsis F has a wide distribution in northern South America and might play an important role in the transmission dynamics of malaria due to its high expansion capacity. Future studies are required to establish the southern distribution of An. albitarsis F in Venezuela, and its occurrence in Guyana and Ecuador.

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

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

Morphometric variations of two populations of Anopheles albitarsis F (Diptera: Culicidae) in the Orinoquia region, Colombia

Introducción. Anopheles albitarsis F ha sido registrada con infección natural por Plasmodium falciparum. En la Orinoquía, se encuentra en simpatría con Anopheles darlingi, un vector primario de malaria.Objetivo. Determinar variaciones morfométricas en caracteres diagnósticos de ala y pata posterior de mosquitos adultos procedentes de dos poblaciones naturales de la Orinoquía, donde la malaria constituye un problema de salud pública.Materiales y métodos. Se analizaron, mediante morfometría lineal y geométrica, los patrones de manchas de la vena Costa del ala de hembras silvestres, de hembras obtenidas de series entomológicas y de isofamilias. Se estimó la forma y el tamaño alar. El segundo tarsómero posterior se analizó por morfometría lineal.Resultados. Se observaron dos patrones de manchas de la vena Costa; se encontraron diferencias estadísticamente significativas en el tamaño y forma alar (p<0.05) entre las hembras silvestres y las obtenidas en condiciones de laboratorio. Para las poblaciones estudiadas, se proponen como caracteres diagnósticos la proporción PHD/HP de 1.17-1.19 en la vena Costa y 0.38-0.54 para la proporción DSIII2/TaIII2 del segundo tarsómero posterior de la pata posterior.Conclusión. Las poblaciones estudiadas de An. albitarsis F presentaron diferencias significativas entre las hembras silvestres y las obtenidas en condiciones de laboratorio, no solo en el tamaño de las alas, sino también en la forma. Se presentan nuevos rangos para los caracteres diagnósticos de la especie.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Molecular operational taxonomic units of mosquitoes (Diptera: Culicidae) collected in high Andean mountain ecosystems of Antioquia, Colombia

Accurate taxonomic identification of highland mosquito species may be complicated because of the lack of comprehensive regional morphological keys and taxonomic specialists, particularly for mosquitoes of medical or ecological importance. We applied a multi-locus approach to explore the diversity of genera/species collected, to define the Molecular Operational Taxonomic Units (MOTUs) and to perform phylogenetic clustering. Twenty MOTUs and three single sequences were revealed from 78 concatenated cox1 + ITS2 sequences, and the species name was allocated for five of these. This study provides molecular taxonomic information of culicid fauna present in high Andean mountain ecosystems in Antioquia, Colombia. However, future morphological and integrative taxonomic studies should be conducted to achieve the specific identity of all detected MOTUs.

Potential distribution of mosquito vector species in a primary malaria endemic region of Colombia

Rapid transformation of natural ecosystems changes ecological conditions for important human disease vector species; therefore, an essential task is to identify and understand the variables that shape distributions of these species to optimize efforts toward control and mitigation. Ecological niche modeling was used to estimate the potential distribution and to assess hypotheses of niche similarity among the three main malaria vector species in northern Colombia: Anopheles nuneztovari, An. albimanus, and An. darlingi. Georeferenced point collection data and remotely sensed, fine-resolution satellite imagery were integrated across the Urabá -Bajo Cauca-Alto Sinú malaria endemic area using a maximum entropy algorithm. Results showed that An. nuneztovari has the widest geographic distribution, occupying almost the entire study region; this niche breadth is probably related to the ability of this species to colonize both, natural and disturbed environments. The model for An. darlingi showed that most suitable localities for this species in Bajo Cauca were along the Cauca and Nechí river. The riparian ecosystems in this region and the potential for rapid adaptation by this species to novel environments, may favor the establishment of populations of this species. Apparently, the three main Colombian Anopheles vector species in this endemic area do not occupy environments either with high seasonality, or with low seasonality and high NDVI values. Estimated overlap in geographic space between An. nuneztovari and An. albimanus indicated broad spatial and environmental similarity between these species. An. nuneztovari has a broader niche and potential distribution. Dispersal ability of these species and their ability to occupy diverse environmental situations may facilitate sympatry across many environmental and geographic contexts. These model results may be useful for the design and implementation of malaria species-specific vector control interventions optimized for this important malaria region.

Use of DNA barcoding to distinguish the malaria vector Anopheles neivai in Colombia

A reference 535 bp barcode sequence from a fragment of the mitochondrial gene cytochrome oxidase I (COI), acquired from specimens of An. neivai Howard, Dyar & Knab, 1913 from its type locality in Panama, was used as a tool for distinguishing this species from others in the subgenus Kerteszia. Comparisons with corresponding regions of COI between An. neivai and other species in the subgenus (An. bellator Dyar & Knab 1906, An. homunculus Komp 1937, An cruzii Dyar & Knab, 1908 and An. laneanus Corrêa & Cerqueira, 1944) produced K2P genetic distances of 8.3-12.6%, values well above those associated with intraspecific variation. In contrast, genetic distances among 55 specimens from five municipalities in the Colombian Pacific coastal state of Chocó were all within the range of 0-2.5%, with an optimized barcode threshold of 1.3%, the limit for unambiguous differentiation of An. neivai. Among specimens from the Chocó region, 18 haplotypes were detected, two of which were widely distributed over the municipalities sampled. The barcode sequence permits discrimination of An. neivai from sympatric species and indicates genetic variability within the species; aspects key to malaria surveillance and control as well as defining geographic distribution and dispersion patterns.

Phenotypic and genotypic variations among three allopatric populations of Lutzomyia umbratilis, main vector of Leishmania guyanensis

BACKGROUND

In South America, Lutzomyia umbratilis is the main vector of Leishmania guyanensis, one of the species involved in the transmission of American tegumentary leishmaniasis. In Brazil, L. umbratilis has been recorded in the Amazon region, and in the state of Pernambuco, Northeastern region, where an isolated population has been identified. This study assessed the phylogeographic structure and size and shape differences of the wing of three Brazilian populations.

METHODS

Samples of L. umbratilis were collected from Rio Preto da Eva (north of the Amazon River, Amazonas), from Manacapuru (south of the Amazon River), and from the isolated population in Recife, Pernambuco state. These samples were processed to obtain sequences of the Cytochrome Oxidase I mitochondrial gene. Geometrics morphometry analysis of the right wing shape of the three populations was made using discriminate canonical analysis.

RESULTS

Phylogenetic analysis revealed the presence of two distinct monophyletic clades: one clade comprised of the Recife and Rio Preto da Eva samples, and the other clade comprised of the Manacapuru samples. Comparing the Manacapuru population with the Recife and Rio Preto da Eva populations generated high indices of interpopulational divergence. Geometric morphometry analysis indicated two distinct groups between the studied populations. Canonical variate analysis of wing shape indicated that Rio Preto da Eva population is significantly closer to Recife population, and both populations were genetically distant from Manacapuru.

CONCLUSION

The polymorphic sites and geometric morphometry analysis indicate that the distance, lack of continuity and environmental differences have not modified the ancestral relationship between Recife and Rio Preto da Eva populations. The genetic and morphological similarities shared by the Recife and Rio Preto da Eva populations suggest that these populations are more closely related evolutionarily. These results confirm the existence of an L. umbratilis species complex in the North and Northeast regions.

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

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

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

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

Malaria vectors in ecologically heterogeneous localities of the Colombian Pacific region

The Colombian Pacific region is second nationally in number of malaria cases reported. This zone presents great ecological heterogeneity and Anopheles species diversity. However, little is known about the current spatial and temporal distribution of vector species. This study, conducted in three ecologically different localities of the Pacific region, aimed to evaluate the composition and distribution of Anopheles species and characterize transmission intensity. A total of 4,016 Anopheles mosquitoes were collected representing seven species. The composition and dominant species differed in each locality. Three species were infected with malaria parasites: Anopheles darlingi and An. calderoni were infected with Plasmodium falciparum and An. nuneztovari with Plasmodium vivax VK210 and VK247. Annual EIRs varied from 3.5–7.2 infective bites per year. These results confirm the importance of the primary vector An. nuneztovari in areas disturbed by human interventions, of An. darlingi in deforested margins of humid tropical rainforest and An. albimanus and the suspected vector An. calderoni in areas impacted by urbanization and large-scale palm oil agriculture close to the coast. This constitutes the first report in the Colombia Pacific region of naturally infected An. darlingi, and in Colombia of naturally infected An. calderoni. Further studies should evaluate the epidemiological importance of An. calderoni in the Pacific region.