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

Molecular diversity of genes related to biological rhythms (period and timeless) and insecticide resistance (Na V and ace-1) in Anopheles darlingi

BACKGROUND

Malaria is a public health concern in the Amazonian Region, where Anopheles darlingi is the main vector of Plasmodium spp. Several studies hypothesised the existence of cryptic species in An. darlingi, considering variations in behaviour, morphological and genetic aspects. Determining their overall genetic background for vector competence, insecticide resistance, and other elements is essential to better guide strategies for malaria control.

OBJECTIVES

This study aimed to evaluate the molecular diversity in genes related to behaviour and insecticide resistance, estimating genetic differentiation in An. darlingi populations from Amazonian localities in Brazil and Pacific Colombian region.

METHODS

We amplified, cloned and sequenced fragments of genes related to behaviour: timeless (tim) and period (per), and to insecticide resistance: voltage-gated sodium channel (Na V ) and acetylcholinesterase (ace-1) from 516 An. darlingi DNA samples from Manaus, Unini River, Jaú River and Porto Velho - Brazil, and Chocó - Colombia. We discriminated single nucleotide polymorphisms (SNPs), determined haplotypes and evaluate the phylogenetic relationship among the populations.

FINDINGS

The genes per, tim and ace-1 were more polymorphic than Na V . The classical kdr and ace-1 R mutations were not observed. Phylogenetic analyses suggested a significant differentiation between An. darlingi populations from Brazil and Colombia, except for the Na V gene. There was a geographic differentiation within Brazilian populations considering per and ace-1.

CONCLUSIONS

Our results add genetic data to the discussion about polymorphisms at population levels in An. darlingi. The search for insecticide resistance-related mechanisms should be extended to more populations, especially from localities with a vector control failure scenario.

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.

Genetic structure and phenotypic variation of Anopheles darlingi in northwest Colombia

This study evaluated the influence of environmental heterogeneity on Anopheles darlingi genetic and morphometric traits at a microgeographic level. Specimens of An. darlingi collected from multiple municipalities in the Colombian malaria endemic region Urabá-Bajo Cauca and Alto Sinú (UCS) were analyzed using 13 microsatellite loci. Spatial genetic structure, population variation and wing geometric morphometric analyses were performed. Microsatellite results showed low genetic differentiation and high gene flow among populations; four highly admixed subpopulations were detected with no particular association to the municipalities. Wing geometric morphometrics analysis showed a subtle but significant difference in wing shape for El Bagre vs. Mutatá populations, possibly influenced by geographical distance. Discrimination among populations in the morphospace showed a slight separation of the Tierralta population. There was no significant correlation between the genetic and geographic or genetic and environmental distances. We hypothesize that environmental heterogeneity in the UCS region does not reach a threshold to affect population structure of An. darlingi. Another possibility is that microsatellites are not sensitive enough to detect existing structure. It remains to be determined which local factors govern phenotypic variation among these populations and how, or whether these may affect mosquito biology and transmission capacity.

Genetic Differentiation of Colombian Populations of Anopheles darlingi Root (Diptera: Culicidae)

Anopheles darlingi Root is a primary vector of malaria in the neotropic region, a species not just highly anthropophilic but very efficient in transmitting Plasmodium species and considered the most important vector in the Amazon region. The main goal of this study was to determine the genetic structure of the A. darlingi populations using microsatellites (STR) in western and eastern regions of Colombia. DNA extraction was done with the cited protocol of band using the Genomic Prep™ cell and tissue isolation commercial kits. We used the STR reported by Conn et al (Mol Ecol Notes 1: 223-225, 2001). The analysis with STR proved there was a high genetic diversity and significant alterations of the Hardy-Weinberg equilibrium. The greatest genetic diversity was recorded in Mitu (Vaupes) (Na = 14, Ho = 0.520). The lowest was in Pueblo Nuevo (Cordoba) (Na = 12, Ho = 0.457). The eastern region and the Mitu (Vaupes) populations presented the highest number of primer alleles (Ap = 30; Ap = 13; Ap = 9), with variations between 0.010 and 0.097. The AMOVA revealed that the whole population underwent moderate genetic differentiation (F _ST = 0.063, p < 0.05). The same differentiation was noticed (0.06 < F _ST > 0.06, p < 0.05) with five of the six populations included in this job, and there was a low differentiation in the Las Margaritas (Santander) area (F _ST = 0.02s3, p < 0.05). Our results suggest a slight positive correlation, which does not show a statistical significance between the geographic and genetic distances, probably suggesting that the moderate genetic differentiation found between pairs of populations does not need to be explained for the hypothesis of separation by distance.

Behavior and population structure of Anopheles darlingi in Colombia

Anopheles darlingi is a widely distributed and important malaria vector in Colombia. Biogeographical and ecological heterogeneity across the Colombian distribution led to the hypothesis of behavioral and genetic differentiation among A. darlingi populations. A total of 2017 A. darlingi specimens were collected during 222 h of sampling. This vector was the most abundant anopheline species in most of the localities sampled. Subdivision between samples collected west and east of the Andes was indicated by 1) mitochondrial COI and nuclear CAD sequences from NW-W and CE-S populations (COI ΦST=0.48761-0.81974, CAD FST=0.11319-0.21321), 2) a COI haplotype network, and 3) SAMOVA. Endo- and exophagy were detected in populations west of the Andes, whereas exophagy was evident in PTG, a locality east of the Andes. Isolation by resistance was significant for COI and explained 26% of the genetic differentiation. We suggest that at a macrogeographic scale, the Andes influence the differentiation of A. darlingi in Colombia and may drive divergence, and, at a microgeographic scale, ecological differences have a significant impact on structure. These data could constitute a baseline for the design of effective vector interventions, locality-specific for the east and similar for panmictic populations west of the Andes.

Inter- and intraspecific identification of the screwworm, Cochliomyia hominivorax, using random amplified polymorphic DNA-polymerase chain reaction

The screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), is one of the most devastating arthropod pests of livestock in the Western Hemisphere. Early instars are very difficult to distinguish morphologically from several closely related blow fly species. Random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) markers were developed for identifying C. hominivorax from other wound inhabiting species. Forty decameric primers were screened; nine showed clear reproducible RAPD profiles suitable for distinguishing all life stages of C. hominivorax from 7 other species, including C. macellaria (Fabricius). The results from RAPD-PCR with field-collected samples of unknown first instars agreed with morphological identification that the samples were not C. hominivorax. Three different primers showed DNA polymorphisms (intraspecific) for samples originating from Mexico, Costa Rica, Panama, Jamaica, and Brazil. Therefore, RAPD-PCR may be useful for determining the geographic origin of C. hominivorax samples. Comparing products from these primers, used with known and unknown screwworm samples from an outbreak in Mexico, clearly showed that the outbreak did not originate from the mass rearing facility. Accurate identification of suspected C. hominivorax samples is possible using RAPD-PCR. Further development to identify the geographic origin of samples would benefit the ongoing surveillance programs against C. hominivorax and the decision process during suspected outbreaks of this important pest.

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.

Genetic differences based on AFLP markers in the mosquito species Anopheles darlingi collected in versus near houses in the region of Porto Velho, RO, Brazil

Anopheles darlingi is the most important malaria vector in Central and South America. After a dramatic reduction of malaria cases in the whole Brazilian territory, with the lowest abundance being attained by 1970, the disease resurged in the Amazon region, where it is now a great public health concern. Consequently, better knowledge about vector species became urgent. We examined the genetic diversity and population structure of A. darlingi, sampled in four localities in the State of Rondônia, Brazil, using 139 amplified fragment length polymorphism marker loci. In each locality, samples were collected in two environments: a peri-domicile one (in the balconies of houses) and an extra-domicile environment (about 15 m from the house). Estimates of expected heterozygosity, Shannon diversity index and percentage of polymorphic loci showed medium to high values, with the samples from the areas closer to Porto Velho exhibiting the smallest values. There was evidence of small population differences, evaluated by F(st), genetic distance and analysis of molecular variance. Comparison between peri- and extra-domicile samples showed greater values of F(st) and genetic distance than between pairs of localities, indicating influence of environmental conditions on the genetics of populations.

Microgeographic genetic variation of the malaria vector Anopheles darlingi root (Diptera: Culicidae) from Cordoba and Antioquia, Colombia

Anopheles darlingi is an important vector of Plasmodium spp. in several malaria-endemic regions of Colombia. This study was conducted to test genetic variation of An. darlingi at a microgeographic scale (approximately 100 km) from localities in Córdoba and Antioquia states, in western Colombia, to better understand the potential contribution of population genetics to local malaria control programs. Microsatellite loci: nuclear white and cytochrome oxidase subunit I (COI) gene sequences were analyzed. The northern white gene lineage was exclusively distributed in Córdoba and Antioquia and shared COI haplotypes were highly represented in mosquitoes from both states. COI analyses showed these An. darlingi are genetically closer to Central American populations than southern South American populations. Overall microsatellites and COI analysis showed low to moderate genetic differentiation among populations in northwestern Colombia. Given the existence of high gene flow between An. darlingi populations of Córdoba and Antioquia, integrated vector control strategies could be developed in this region of Colombia.

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.

Population genetic structure of the major malaria vector Anopheles darlingi (Diptera: Culicidae) from the Brazilian Amazon, using microsatellite markers

The population genetic structure of Anopheles darlingi, the major human malaria vector in the Neotropics, was examined using seven microsatellite loci from nine localities in central and western Amazonian Brazil. High levels of genetic variability were detected (5-25 alleles per locus; H E = 0.519-0.949). There was deviation from Hardy-Weinberg Equilibrium for 59.79% of the tests due to heterozygote deficits, while the analysis of linkage disequilibrium was significant for only two of 189 (1.05%) tests, most likely caused by null alleles. Genetic differentiation (F ST = 0.001-0.095; Nm = 4.7-363.8) indicates that gene flow is extensive among locations < 152 km apart (with two exceptions) and reduced, but not absent, at a larger geographic scale. Genetic and geographic distances were significantly correlated (R(2) = 0.893, P < 0.0002), supporting the isolation by distance (IBD) model. The overall estimate of Ne was 202.4 individuals under the linkage disequilibrium model, and 8 under the heterozygote excess model. Analysis of molecular variance showed that nearly all variation (approximately 94%) was within sample locations. The UPGMA phenogram clustered the samples geographically, with one branch including 5/6 of the state of Amazonas localities and the other branch the Acre, Rondônia, and remaining Amazonas localities. Taken together, these data suggest little genetic structure for An. darlingi from central and western Amazonian Brazil. These findings also imply that the IBD model explains nearly all of the differentiation detected. In practical terms, populations of An. darlingi at distances < 152 km should respond similarly to vector control measures, because of high gene flow.