Random amplified polymorphic DNA analysis of Anopheles nuneztovari (Diptera: Culicidae) from Western and northeastern Colombia

GENETIC DIVERSITY OF Lutzomyia (Nyssomyia) intermedia IN AN ENDEMIC AREA OF AMERICAN CUTANEOUS LEISHMANIASIS, STATE OF PARANÁ, BRAZIL

Lutzomyia intermedia (Diptera: Psychodidae) features as one of the main vectors that are involved in the transmission of American cutaneous leishmaniasis (ACL) in the Neotropical region. However, genetic studies involving this taxon are still incipient and important for understanding the level of variability of different populations, their role, and implications as vectors. The aim of this study was to determine the level of genetic diversity of L. intermedia present in the Ribeira River Valley, an area of ACL transmission in the state of Paraná, Brazil, through the Random amplified polymorphic DNA (RAPD). Two municipalities were chosen to collect sand flies: Cerro Azul (new transmission area of the ACL) and Adrianópolis (endemic area of the ACL). The insects were captured in the house, in the peridomicile and in the wild (forest). Two of the used markers made it possible to estimate the polymorphism of the studied populations, resulting in 40 genotypes, most of them from peridomicile. The dendrogram generated by the analysis with the primer A10 showed different degrees of similarity, suggesting that there may be gene flow in the studied populations. The Principal Coordinate Analysis (PCO) with the A2 primer, was useful in grouping L. intermedia according to its ecological and geographical origin. There was no distinction between the lineages composing the L. intermedia complex. The results of this study, with the record of great genotypic diversity in L. intermedia, may contribute to explain the maintenance of the life cycle of Leishmania braziliensis (Kinetoplastida: Trypanosomatidae) in the region.

Potential distribution of main malaria vector species in the endemic Colombian Pacific region

OBJECTIVE

To assess the existing fundamental niche, potential distribution and degree of niche overlap for the three main Colombian malaria vectors Anopheles albimanus, Anopheles darlingi and Anopheles nuneztovari in the major malaria endemic Pacific region.

METHODS

We used models based on presence records and Normalised Difference Vegetation Index (NDVI) data, created using the maximum entropy algorithm.

RESULTS

The three vector species occupied heterogeneous environments, and their NDVI values differed. Anopheles albimanus had the largest niche amplitude and was distributed mainly on coastal areas. Environmentally suitable areas for An. albimanus and An. nuneztovari were the dry forest of inter-Andean Valleys in south-western Colombia, as confirmed for An. albimanus during model validation. There was a slight degree of niche overlap between An. darlingi and An. nuneztovari, and the species co-occurred in humid forests, predominantly in riparian zones of the San Juan and Atrato rivers.

CONCLUSION

The information obtained may be used for the implementation of vector control interventions in selected priority areas to reduce malaria risk in this region while optimising resources.

Monitoring of the Apple Fruit Moth: Detection of Genetic Variation and Structure Applying a Novel Multiplex Set of 19 STR Markers

The apple fruit moth Argyresthia conjugella (Lepidoptera, Yponomeutidae) is a seed predator of rowan (Sorbus aucuparia) and is distributed in Europe and Asia. In Fennoscandia (Finland, Norway and Sweden), rowan fruit production is low every 2–4 years, and apple (Malus domestica) functions as an alternative host, resulting in economic loss in apple crops in inter-mast years. We have used Illumina MiSeq sequencing to identify a set of 19 novel tetra-nucleotide short tandem repeats (STRs) in Argyresthia conjugella. Such motifs are recommended for genetic monitoring, which may help to determine the eco-evolutionary processes acting on this pest insect. The 19 STRs were optimized and amplified into five multiplex PCR reactions. We tested individuals collected from Norway and Sweden (n = 64), and detected very high genetic variation (average 13.6 alleles, He = 0.75) compared to most other Lepidoptera species studied so far. Spatial genetic differentiation was low and gene flow was high in the test populations, although two non-spatial clusters could be detected. We conclude that this set of genetic markers may be a useful resource for population genetic monitoring of this economical important insect species.

Population dynamics of Anopheles nuneztovari in Colombia

Anopheles nuneztovari is an important Colombian malaria vector widespread on both sides of the Andean Mountains, presenting morphological, behavioral and genetic heterogeneity throughout the country. The aim of this study was to evaluate whether the population structure and distribution of An. nuneztovari in Colombia are associated with ecological and physical barriers present in a heterogeneous landscape. Further, differences in behavior were addressed. A total of 5392 specimens of An. nuneztovari were collected. Mitochondrial and nuclear marker analyses detected subdivision among the northwest-west, northeast and east populations. For both markers, isolation by distance (~53%) and isolation by resistance (>30%) were determinants of population genetic differentiation. This suggests that physical barriers, geographical distance and ecological differences on both sides of the Andean Mountains promoted the genetic differentiation and population subdivision of An. nuneztovari in Colombia. This species showed the highest biting activity after 20:00h; indoor and outdoor preferences were found in all localities. These results indicated that the most effective interventions for controlling vector populations on both sides of the Andes need to be region-specific.

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.

Malaria vector species in Colombia: a review

Here we present a comprehensive review of the literature on the vectorial importance of the major Anopheles malaria vectors in Colombia. We provide basic information on the geographical distribution, altitudinal range, immature habitats, adult behaviour, feeding preferences and anthropophily, endophily and infectivity rates. We additionally review information on the life cycle, longevity and population fluctuation of Colombian Anopheles species. Emphasis was placed on the primary vectors that have been epidemiologically incriminated in malaria transmission: Anopheles darlingi, Anopheles albimanus and Anopheles nuneztovari. The role of a selection of local, regional or secondary vectors (e.g., Anopheles pseudopunctipennis and Anopheles neivai) is also discussed. We highlight the importance of combining biological, morphological and molecular data for the correct taxonomical determination of a given species, particularly for members of the species complexes. We likewise emphasise the importance of studying the bionomics of primary and secondary vectors along with an examination of the local conditions affecting the transmission of malaria. The presence and spread of the major vectors and the emergence of secondary species capable of transmitting human Plasmodia are of great interest. When selecting control measures, the anopheline diversity in the region must be considered. Variation in macroclimate conditions over a species' geographical range must be well understood and targeted to plan effective control measures based on the population dynamics of the local Anopheles species.

Population structure analyses and demographic history of the malaria vector Anopheles albimanus from the Caribbean and the Pacific regions of Colombia

BACKGROUND

Anopheles albimanus is an important malaria vector in some areas throughout its distribution in the Caribbean and the Pacific regions of Colombia, covering three biogeographic zones of the neotropical region, Maracaibo, Magdalena and Chocó.

METHODS

This study was conducted to estimate intra-population genetic diversity, genetic differentiation and demographic history of An. albimanus populations because knowledge of vector population structure is a useful tool to guide malaria control programmes. Analyses were based on mtDNA COI gene sequences and four microsatellite loci of individuals collected in eight populations from the Caribbean and the Pacific regions of Colombia.

RESULTS

Two distinctive groups were consistently detected corresponding to COI haplotypes from each region. A star-shaped statistical parsimony network, significant and unimodal mismatch distribution, and significant negative neutrality tests together suggest a past demographic expansion or a selective sweep in An. albimanus from the Caribbean coast approximately 21,994 years ago during the late Pleistocene. Overall moderate to low genetic differentiation was observed between populations within each region. However, a significant level of differentiation among the populations closer to Buenaventura in the Pacific region was observed. The isolation by distance model best explained genetic differentiation among the Caribbean region localities: Los Achiotes, Santa Rosa de Lima and Moñitos, but it could not explain the genetic differentiation observed between Turbo (Magdalena providence), and the Pacific region localities (Nuquí, Buenaventura, Tumaco). The patterns of differentiation in the populations from the different biogeographic provinces could not be entirely attributed to isolation by distance.

CONCLUSION

The data provide evidence for limited past gene flow between the Caribbean and the Pacific regions, as estimated by mtDNA sequences and current gene flow patterns among An. albimanus populations as measured by MS loci which may be mainly influenced by semi-permeable natural barriers in each biogeographical region that lead to the genetic differences and effective population sizes detected. The relatively high genetic differentiation in the port city of Buenaventura may be the result of specific ecological conditions, human migration and activities and/or differences in effective population sizes. This knowledge could serve to evaluate and coordinate vector control strategies in these regions of Colombia.

Resurrection of Anopheles goeldii from synonymy with Anopheles nuneztovari (Diptera, Culicidae) and a new record for Anopheles dunhami in the Brazilian Amazon

Nucleotide sequences of the internal transcribed spacer 2 (ITS2) rDNA and partial sequences of the cytochrome coxidase subunit I (COI) mtDNA and white gene nDNA were obtained from specimens of Anopheles nuneztovari A collected in Macapá (state of Amapá), Obidos, Prainha and Almeirim (state of Pará), Itacoatiara and Parintins (state of Amazonas), Brazil, and compared with previously published sequences of A. nuneztovari s.l. Results of the Bayesian phylogenetic analyses performed using either COI or combined ITS2, COI and white gene sequences suggest that An. nuneztovari B/C is distinct from specimens obtained in the Amazonas/Solimões River basin. Anopheles goeldii, currently in synonymy with An. nuneztovari, was described from individuals collected in Belterra (= Fordlândia) in the Tapajós River, state of Pará, Southern Amazonas River. Morphological comparisons of the characteristics of the male genitalia indicated that An. nuneztovari A and An. goeldii are similar but distinct from An. nuneztovariB/C by the apex of the aedeagus. In considering the results of the phylogenetic analyses and morphological comparisons, An. goeldii is resurrected from synonymy with An. nuneztovari. Additionally, Anopheles dunhamiis reported for the first time in Parintins. This species can be distinguished from An. goeldiiby characters of the male genitalia and molecular data.

Morphological analysis of three populations of Anopheles (Nyssorhynchus) nuneztovari Gabaldón (Diptera: Culicidae) from Colombia

Based on the results of comparative analyses of 1,039 specimens of several progenies of Anopheles nuneztovarifrom three localities in Colombia, eight costal wing spot patterns were observed. Patterns I and III were the most frequent: 77.96% and 11.36%, respectively. Using the diagnostic characters ratio of the length of the basal dark area of hind tarsomere II/length of hind tarsomere II, ratio of the length of the humeral pale spot/length of the pre-humeral dark spot, and the ratio of the length of the subcostal pale spot/length of the distal sector dark spot (DS-III2/Ta-III2, HP/PHD, SCP/DSD) approximately 5% of the adult females were misidentified as a species of Nyssorhynchus, different from An. nuneztovari. Approximately 5% of the specimens showed DS-III2/Ta-III2 ratio less than 0.25 (range 0.21 - 0.24), and among them 3.34% shared a HP/PHD ratio less than 1.50. Consequently, 1.52% of An. nuneztovari individuals can be misidentified as Anopheles oswaldoi. In those specimens with the DS-III2/Ta-III2 ratios higher than 0.25, 34.45% displayed SCP/DSD values greater than 0.50 and of these, 3.65% displayed HP/PHD values greater than 1.8. This combination of characters could lead one to misidentify samples of An. nuneztovari as Anopheles rangeli. Similarly, 2.43% of the females could be identified erroneously as either Anopheles aquasalis or Anopheles benarrochi. Individuals with a HP/PHD ratio greater than 2.0, could be misidentified as Anopheles trinkae, Anopheles strodei or Anopheles evansae. A distinct combination of diagnostic characters for An. nuneztovari from Colombia is proposed.

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.

Genetic structure of natural populations of the sand fly Lutzomyia longipalpis (Diptera: Psychodidae) from the Brazilian northeastern region

In Latin America, Lutzomyia longipalpis is the principal vector of Leishmania chagasi, and is associated with the majority of active foci of visceral leishmaniasis. In spite of the fact that this sand fly is spread practically throughout the entire Neotropical Region, its distribution is not uniform due to geographic and environmental barriers. Geographic isolation coupled with reduced flight abilities may contribute to the appearance of cryptic species of Lutzomyia longipalpis, which may differ in their capacity to transmit L. chagasi. In this work, we describe the genetic structuring patterns based on polymorphism analysis of 24 RAPD-PCR loci of 7 natural populations of Lutzomyia longipalpis obtained from Brazil's northeastern region. The estimated degree of genetic differentiation between populations, based on the population subdivision index theta(ST) (0.136), suggests a moderate degree of genetic structuring as a result of geographical isolation and restricted gene flow. Genetic distances were found to be compatible with those found between members of a single species, suggesting a taxonomic uniformity of Lutzomyia longipalpis in the region studied.

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.