Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

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

Abstract

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.

Introduction

Malaria is an infectious disease caused by protozoans of the genus Plasmodium; it is transmitted to humans by bites of female mosquitoes of the genus Anopheles [1]. In Latin America, Colombia occupies the third place in the number of malaria cases after Brazil and Venezuela, with 10% of the total number of reported cases [2]. In particular, the Urabá-Bajo Cauca and Alto Sinú (UCS, Fig 1) regions have often reported the highest numbers of malaria cases in Colombia [3]; only recently, in 2015, this area dropped to second in number, with 16.6% of total cases in the country [4]. For decades, Plasmodium vivax has been the predominant malaria parasite in UCS, representing an annual parasite index (API) of 28.7/1,000 inhabitants in 2015; while the API for Plasmodium falciparum was 10.2/1,000 [5]. The three main Colombian vectors Anopheles darlingi, Anopheles nuneztovari and Anopheles albimanus have an important role in malaria transmission in the UCS region [3,6,7].

thumbnail
Fig 1. Study area and collection sites for Anopheles nuneztovari, Anopheles darlingi, and Anopheles albimanus across the Urabá-Bajo Cauca and Alto Sinú endemic region (UCS), Colombia.

https://doi.org/10.1371/journal.pone.0179093.g001

Ecological niche models (ENMs) are used to understand ecological requirements of species, aspects of their biogeography, predict geographic distributions, identify sites for translocations and reintroductions, select areas for conservation and forecast effects of environmental change [8]. Studies of geographic distributions of vector species using ENM techniques relate occurrence records and environmental characteristics across species’ distributions [912]. For instance, ENM analyses of anopheline species (subgenus Nyssorhynchus) in Amazonian Brazil revealed diversification in habitat-use: An. triannulatus is a generalist, whereas An. oryzalimentes and An. janconnae are specialists [13]. ENMs were also used to predict distributions of An. bellator, An. cruzii, and An. marajoara of the Riviera Valley in southern Brazil, which revealed specific associations with land cover types [14]. Finally, low tolerance to dry environments was documented for An. darlingi; projected climate change would significantly reduce its suitable habitat mainly in Amazonian biomes, influencing both its distribution and abundance, in contrast to species of the Albitarsis complex [15].

Implementation of effective methods for ecological and epidemiological data analysis is essential in a country like Colombia, considering the rapid reduction of native habitats and their conversion to agriculture, livestock, and mining uses [16]. These changes make it a priority to identify variables that determine and constrain distributions of disease vector species. The limited resources available for regular monitoring of vector species and vector-borne diseases and control of infectious diseases would be improved by development of an early warning system regarding transmission and outbreaks of vector-borne diseases. Hence, we used ENM to evaluate potential vector distributions and to assess the hypothesis that environmental heterogeneity is a driver for allopatric distributions of possible competing niche-related species.

Materials and methods

Occurrence data

The endemic malaria region Urabá Bajo-Cauca Alto-Sinú (UCS) is located at northwest Colombian (Fig 1), has an estimated population at-risk of malaria equivalent to 2,500,000 people, distributed in the departments of Antioquia and Córdoba [17]. The predominant economic activities are mining, the banana agro-industry, cattle production and the timber industry; some of these activities have been linked to malaria transmission [3]. Twenty municipalities of UCS were sampled (Fig 1) between December 2012 and March 2015, using human landing-catches, under an informed consent agreement and collection protocol reviewed and approved by the Universidad de Antioquia Institutional Review Board (Bioethics Committee, Facultad Nacional de Salud Pública-Universidad de Antioquia, Acta 063). In light of the risk that these vectors represent for the transmission of the disease to humans, collections were performed in or near human residences. Two localities were selected per municipality (Table 1), based on criteria such as high numbers of reported malaria cases in 2012 and accessibility. Collections were conducted during three nights at each locality, by a team of four people, two persons sampling one site each night, for a total of six sites by locality, from 18:00 to 24:00 h, both indoors and outdoors. Adult mosquitoes were identified using a morphological key [18]; species assignation for a percentage of the specimens and all those presenting problems in morphological identification were confirmed by PCR-RFLP-ITS2 [1921] and COI barcoding [2224].

thumbnail
Table 1. Summary of mosquito collection data at 40 localities of Urabá-Bajo Cauca and Alto Sinú region, Colombia.

https://doi.org/10.1371/journal.pone.0179093.t001

Environmental variables

To characterize environmental variation across the study region, the Normalized Difference Vegetation Index (NDVI) was used. NDVI is a measure of photosynthetic activity; its variation reflects the spatial and temporal dynamics of vegetation that indirectly influence mosquito reproduction and development [25]. NDVI data were obtained in the form of imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra satellite, at 250 m spatial resolution and 16-day temporal resolution (https://reverb.echo.nasa.gov/reverb/). In all, 69 images from 2012 to 2014 were used. Images were reprojected to the MAGNA-SIRGAS coordinate system using the MODIS Reprojection Tool [26]. To reduce inter-correlation among data layers, a principal components analysis (PCA) was performed using all 69 NDVI images as variables [27]. Climatic information was not included, owing to the coarse spatial resolution of such data sets, which would not permit fine-grained predictions. Regardless, rainfall is indeed reflected in the NDVI data via vegetational responses and consequently increased or decreased photosynthetic mass [28].

Ecological niche modeling

To evaluate potential distributions of the vector species, ENMs were developed using a maximum entropy approach implemented in Maxent, version 3.3.3k [29]. Maxent assesses suitability for species by integrating occurrence records with relevant environmental predictors across a defined geographic space [29]. Maxent attempts to estimate a target probability distribution by finding the probability distribution of maximum entropy (i.e., that which is most spread out, or closest to uniform), subject to a set of constraints that represent the incomplete nature of information about the target distribution [30]. A sampling bias layer was designed to improve niche and range estimates with Maxent and to point process models by integrating spatially explicit information [30], with our sampling locations defined as 1; the remaining pixels that were not sampled had “no data” in the grid. Given the broad potential distribution previously reported for An. nuneztovari, An. albimanus, and An. darlingi in South America [15,31,32] and the known distribution of the species across the country [3,7,18,21], a very broad hypothesis of accessible area (M) was used, considering the entire endemic region as accessible [33].

To assess how robust and predictive the models were, occurrence data for each species were split into training (50%) and testing (50%) subsets; this random splitting was done five times. No clamping or extrapolation was permitted in MaxEnt; remaining parameters were left as default. Logistic output formats were used. ENMs were calibrated for each species, with 10 bootstrapped replicates; the median across replicates was used as a basis for further analysis. All maps were converted to binary via a conservative least presence thresholding approach, consisting of the lowest predicted value corresponding to any occurrence record of the species in the calibration data set [34,35].

The partial receiver operating characteristic (ROC) approach was used to assess models performance [36]; the evaluation dataset was bootstrapped, and probabilities obtained by direct count of area under the curve (AUC) ratios falling ≤ 1 via a Visual Basic script (N. Barve, pers.comm.; https://kuscholarworks.ku.edu), with 100 iterations. Model performance was tested for different combinations of principal components (the first 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, and 55 components). The omission rate and partial ROC AUC ratios were used as criteria to select optimal environmental data sets for each species. Final models were produced for each species using all available data [37]. We inspected the loading values of each raw variable (16-day composite NDVI) on each of the first two principal components and how they related to monthly rainfall averages in the study area to have a better understanding of how vegetation dynamics reflected in NDVI related to suitability for each species. We also compared available (i.e., within M) and occupied environments in terms of NDVI values across the endemic area. The occupied niche was represented in two ways by: (1) considering NDVI values over the study period for each occurrence coordinate, and (2) taking only NDVI values matching the date of the mosquito’s capture.

Niche overlap

To evaluate a hypothesis of niche similarity or difference among the three species within this endemic area, we used a background similarity test implemented in ENMTools v.1.4.3 [38]. The test generates a similarity measure by overlaying predictions of the two species compared [38]; niche overlap values were calculated using the Schoener’s D metric, with values ranging from 0 (no overlap) to 1 (complete overlap). The background similarity test uses the original predictions for each species and compares each of them with overlaps generated from models based on random occurrences selected from across the accesible area for each species. Observed similarity was compared with the distribution of similarities between the focal species and the suite of null models. Probabilities were determined by direct count, and two comparisons were made for each pair of species, with each species serving in turn as the focal species. In all, 100 pseudoreplicates were generated for each pair [38,39]. Finally NicheA was used to visualize overall overlap based on minimum volume ellipsoids for the species in three PCA dimensions [40]; the Jaccard index (IJ) was used as a numerical estimation of environmental overlap among species [41].

Results

Presence records totaled 100 for Anopheles nuneztovari, 69 for Anopheles darlingi, and 88 for Anopheles albimanus. All models developed performed statistically significantly better than random expectations (all p<0.01; S1 Table). According to omission rate and ROC partial values, optimal principal component combinations were, An. nuneztovari first 25 PCAs, An. darlingi first 35 PCAs, and An. albimanus first 25 PCAs.

Potential distributions estimated for the three species indicate that the area with highest environmental suitability for An. darlingi was Bajo Cauca (39% suitable; Fig 2A). In contrast, the broadest potential distribution for An. nuneztovari was in parts of Bajo Cauca and Alto Sinú (48.3% suitable; Fig 2B). The most suitable areas for An. albimanus were in coastal areas of Urabá and Alto Sinú and some areas of Bajo Cauca (31% suitable; Fig 2C).

thumbnail
Fig 2. Ecological niche models.

Showing environmental suitability for the three Colombian main malaria vectors across the Urabá-Bajo Cauca and Alto Sinú endemic region. Red circles indicate location of the records used for model calibration and green squares indicate absence records.

https://doi.org/10.1371/journal.pone.0179093.g002

The relationship between the original variables (NDVI values) and the principal components showed that the first principal component (PC1) was positively associated with NDVI for all dates throughout the three years, while PC2 exhibited positive and negative associations with NDVI in each of the three years. Both PC1 and PC2 exhibited temporal trends associated with NDVI variation, reflecting seasonality and corresponding leaf phenology, partially in response to fluctuation of rain (Fig 3), which is evident in two areas that are unoccupied and that present contrasting seasonality, was well as overall levels of NDVI (Fig 4). Areas with highest PC1 values had high NDVI values and relatively little seasonal fluctuation, whereas areas with the highest PC2 values had clear fluctuations in both years. PC1 and PC2 explained 39.7% and 11.0% of the total variance, respectively. Anopheles albimanus and An. darlingi occupied areas with lower NDVI values throughout the year (Fig 5) as compered than An. nuneztovari. All species were apparently absent from areas with high NDVI values throughout the year (evergreen forests) and areas with very pronounced NDVI fluctuations.

thumbnail
Fig 3. Loading plot from a principal component analysis using NDVI values for the Urabá-Bajo Cauca and Alto Sinú region of Colombia.

Black bars indicate the relationship (positive or negative) of the principal components with the dates of the NDVI values. Black line represents the average monthly rainfall for the study area (data from Instituto de Hidrología, Meteorología y Estudios Ambientales de Colombia IDEAM, Octubre 2016).

https://doi.org/10.1371/journal.pone.0179093.g003

thumbnail
Fig 4. Results of the principal component analysis.

Based on NDVI values for the Urabá-Bajo Cauca and Alto Sinú region of Colombia. Blue circles indicate NDVI values for sites across the study area. Other colored circles represent NDVI values for occurrence data for the main malaria vectors An. nuneztovari, An. albimanus, and An. darlingi. Red arrows indicate regions of available, but unoccupied environmental space. Box plots represent NDVI values in those unoccupied areas of environmental space (1 and 2).

https://doi.org/10.1371/journal.pone.0179093.g004

thumbnail
Fig 5. NDVI values at all occurrence sites for the three malaria vectors in northern Colombia.

The black dots represent specific NDVI values corresponding to mosquito collection dates.

https://doi.org/10.1371/journal.pone.0179093.g005

Assessment of both geographic and environmental dimensions showed high overlap among the three species (D > 0.35). Nonetheless, observed overlap between An. nuneztovari and An. albimanus was significantly lower than expected from a null distribution (p = 0.032 in the background similarity tests, S1 Fig). NicheA visualization revealed broad overlap among the niches of the three species, and showed that niche breadth was greater for An. nuneztovari than for the other two species (Fig 6). The Jaccard index indicated that environmental overlap between An. nuneztovari and An. darlingi was 0.63, that for An. nuneztovari versus An. albimanus was 0.36, and that for An. albimanus versus An. darlingi was 0.13. These results are consistent with Schoener’s D metric.

thumbnail
Fig 6. Environmental space modeled in NicheA.

Visualization in three dimensions in terms of principal components 1, 2 and 3; ecological niche models were displayed as minimum-volume ellipsoids used to illustrate limits of environmental distributions. The background cloud (in gray), was derived from a random sample of 3000 random points from across the study area.

https://doi.org/10.1371/journal.pone.0179093.g006

Discussion

The three main Colombian Anopheles vector species, An. nuneztovari, An. albimanus, and An. darlingi, occur in sympatry in the Urabá-Bajo Cauca and Alto Sinú region [3,7,42]. Using remotely sensed environmental data with high spatial and temporal resolution, we were able to identify small but meaningful differences in environments occupied by these three species. These details can improve the utility of ENMs with the aim of improving the knowledge of the biology of these important species, and guiding integrated vector management efforts [43]. Deforestation and changes in land use are occurring at a higher speed than changes in climate; as a result, use of environmental variables indicative of changes in vegetation cover will be a useful input in the prediction and prevention of vector-borne disease transmission.

These results agree with previous efforts to map Anopheles species’ distributions in the Americas, in which models indicated their presence across our study region [12,15,31,32]. At difference from those studies, we used environmental layers with higher spatial and temporal resolution and thus greater information content at finer scales, allowing documentation of more specific and detailed distributional patterns. The three main Colombian vectors An. darlingi, An. nuneztovari and An. albimanus have an important role in malaria transmission in the UCS region [3,6,7]. Anopheles nuneztovari and An. darlingi, have been reported naturally infected by Plasmodium spp. in UCS [8], and reports have indicated that these two species are able to maintain transmission even at low abundances [4]. While An. albimanus has a predominant coastal distribution and in the Colombian Pacific Coast was detected infected with Plasmodium vivax and Plasmodium falciparum [6]. Our final results showed that An. nuneztovari has the widest geographic distribution, occupying almost the entire study region; this breadth is probably related to the ability of this species to colonize both natural and disturbed environments [3,44]. Anopheles nuneztovari is the most frequently detected species in some localities of the Bajo Cauca and Alto Sinú region [3,7]. Anopheles darlingi is also an important malaria vector in other regions of Colombia [45,46], it has been associated with breeding sites in riverine or gallery forest [4750]. Our model for An. darlingi showed that most suitable localities for this species are in Bajo Cauca; specifically, along the Cauca and Nechí rivers, that conform an important aquifer system in this region [51]. The riparian ecosystems in this subregion and the potential for rapid adaptation by this species to novel environments [49,52,53] may favor the establishment of populations of this species. Particularly, An. darlingi is associated with riparian ecosystems in the upper Orinoco River in southern Venezuela, where overflow of the river creates lagoons that constitute suitable larval habitats [54]. Therefore, the results of the present study suggest the importance of authorities enforcing regulations on deforestation for the Cauca and Nechí rivers basins, where illegal activities such as gold mining may be increasing human contact exposition to An. darlingi, as has been previously suggested [55].

Finally, in Colombia, An. albimanus is a species with a mostly coastal distribution [7,56,57]. The model identified suitable areas for An. albimanus in Urabá and Bajo Cauca. Absence of this species from some areas modeled as suitable may relate to factors such as competitive exclusion, existing vector control measures, or rapid changes in land use as a result of agricultural activities and mining, as previously suggested for this species in Mesoamerica and the Caribbean region [32].

Previous studies have shown that distributions of mosquito species are partly related to land use factors such as presence or absence of wetlands, type of surrounding vegetation, and agricultural practices [58]. Various studies on species modeling are using NDVI instead of land use factors, given that this index is a measure of photosynthetic activity; thus, allows an indirect approximation to the suitable conditions for mosquito development [25]. Our result concerning the relationship between precipitation and NDVI patterns in different time periods showed that PC1 was more indicative of overall NDVI values (ie; forest versus non-forest), whereas PC2 was more indicative of the magnitude of seasonality in NDVI. Apparently, the three main Colombian Anopheles vector species in Urabá-Bajo Cauca and Alto Sinú do not occupy environments either with high seasonality, or with low seasonality and high NDVI. Because NDVI has been used to estimate vegetation characteristics [59], the relation between NDVI and species’ occurrence suggests similar affinities of the three species for particular vegetation types; An. albimanus tended to occupy places with relatively low NDVI values, representing ecosystems with bare soils and low forest cover [60]. Previous reports indicate that An. albimanus is associated with larval habitats exposed to sun [48], which as in the present study, are related to scarce forest cover, and also its larvae can tolerate some salinity [61]. Places with low seasonality and high NDVI values probably represent residual forests unoccupied by people.

The background similarity tests and NicheA analysis showed high overlaps in the geographic and environmental conditions occupied by An. darlingi, An. nuneztovari and An. albimanus. The dispersal capacity of these species and their ability to colonize different ecosystems are well documented [3,53,62,63]. These aspects may facilitate their sympatry across areas presenting diverse environmental conditions.

The models developed in this study have important applications, since they could be projected or replicated for different eco-epidemiological malaria zones of Colombia, as has been done elsewhere in the world [64]. In other applications of these approaches, models of the potential geographic distribution of Ebola and Marburg viruses in Africa [65] were able to anticipate the potential for Marburg outbreaks in Angola [66]. Variation in NDVI in relation to the potential distribution of the vector species may be reflected in patterns of malaria transmission in Urabá-Bajo Cauca and Alto Sinú region. Historically and until recently, this region reported the highest numbers of malaria cases in Colombia [

3]. The main vectors An. nuneztovari and An. darlingi have been detected naturally infected by Plasmodium spp. in UCS [7], and recent reports indicate that these species are the most important malaria vectors in some localities of the region [3]. Vector distribution and malaria transmission might be highly dynamic because of antrophic or natural causes; thus a good characterization of environmental variation, through the use of NDVI at high spatio-temporal resolution [58,67,68] to predict mosquito spreading, will provide baseline information that can be used to reduce malaria risk in this region.

Supporting information

S1 Table. Partial AUC ratios of Anopheles vector species.

https://doi.org/10.1371/journal.pone.0179093.s001

(DOCX)

Acknowledgments

We thank JD Sánchez Rodríguez, Y Espinosa, and J Rodríguez Zabala at Laboratorio de Microbiologia Molecular-UdeA for mosquito collection support. The Biodiversity and Macroecology Laboratory at Kansas University for data analysis support. The Instituto de Hidrología, Meteorología y Estudios Ambientales de Colombia- IDEAM for providing crucial meteorological information.

Author Contributions

  1. Conceptualization: MA MMC.
  2. Formal analysis: MA MMC SA JLP ATP.
  3. Funding acquisition: MMC.
  4. Investigation: MA MMC.
  5. Methodology: MA MMC SA JLP ATP.
  6. Project administration: MMC.
  7. Supervision: MMC.
  8. Writing – original draft: MA MMC.
  9. Writing – review & editing: MA MMC SA JLP ATP.

References

  1. 1. Cohuet A, Harris C, Robert V, Fontenille D. Evolutionary forces on Anopheles: what makes a malaria vector? Trends Parasitol. 2010;26: 130–136. pmid:20056485
  2. 2. World Health Organization. World Malaria Report 2015. Geneva: World Health Organization. 312 p.
  3. 3. Naranjo- Díaz N, Rosero DA, Rua-Uribe G, Luckhart S, Correa MM. Abundance, behavior and entomological inoculation rates of anthropophilic anophelines from a primary Colombian malaria endemic area. Parasit Vectors. 2013;6: 1–11.
  4. 4. Ospina M, Mancel E, Pacheco O, Quijada H. Boletin Epidemiológico Semanal Sivigila. Bogotá;2015: 1–72.
  5. 5. INS. Instituto Nacional de Salud, Boletín epidemiológico Semanal. Estadísticas del sistema de vigilancia en salud pública- SIVIGILA, Casos totales en la Semana Epidemiológica 52 y acumulados del año, Subdirección de Vigilancia y Control en Salud Pública. 2016.
  6. 6. Gutiérrez LA, Naranjo N, Jaramillo LM, Muskus C, Luckhart S, Conn JE, et al. Natural infectivity of Anopheles species from the Pacific and Atlantic Regions of Colombia. Acta Trop. 2008;107: 99–105. pmid:18554564
  7. 7. Gutiérrez L, González JJ, Gómez GF, Castro MI, Rosero DA, Luckhart S, et al. Species composition and natural infectivity of anthropophilic Anopheles (Diptera: Culicidae) in the states of Córdoba and Antioquia, Northwestern Colombia. Mem Inst Oswaldo Cruz. 2009;104: 1117–24. pmid:20140372
  8. 8. Peterson AT, Soberón J, Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, et al. Ecological Niches and Geographic Distributions. Princeton University Press; 2011. 330p.
  9. 9. Foley DH, Wilkerson RC, Birney I, Harrison S, Christensen J, Rueda LM. MosquitoMap and the Mal-area calculator: new web tools to relate mosquito species distribution with vector borne disease. Int J Health Geogr. 2010; 1–8.
  10. 10. Fuller DO, Parenti MS, Hassan AN, Beier JC. Linking land cover and species distribution models to project potential ranges of malaria vectors: an example using Anopheles arabiensis in Sudan and Upper Egypt. Malar J. Malaria Journal; 2012;11: 264. pmid:22866895
  11. 11. Overgaard HJ, Ekbom B, Suwonkerd W, Takagi M. Effect of landscape structure on anopheline mosquito density and diversity in northern Thailand: implications for malaria transmission and control. Landsc Ecol; 2003;18: 605–619.
  12. 12. Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, et al. The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010;3: 1–26.
  13. 13. McKeon SN, Schlichting CD, Povoa MM, Conn JE. Ecological suitability and spatial distribution of five Anopheles species in Amazonian Brazil. Am J Trop Med Hyg. The American Society of Tropical Medicine and Hygiene; 2013;88: 1079–86. pmid:23546804
  14. 14. Laporta GZ, Ramos DG, Ribeiro MC, Sallum MAM. Habitat suitability of Anopheles vector species and association with human malaria in the Atlantic Forest in south-eastern Brazil. Mem Inst Oswaldo Cruz; 2011;106: 239–245. pmid:21881779
  15. 15. Laporta GZ, Linton Y-M, Wilkerson RC, Bergo ES, Nagaki SS, Santana DC, et al. Malaria vectors in South America: current and future scenarios. Parasit Vectors. Parasites & Vectors; 2015;8: 426. pmid:26283539
  16. 16. Cabrera E, Vargas DM, Galindo G. García MC, Ordoñez MF, Vergara, LK, Pacheco, AM, Rubiano, JC. y Giraldo, P. 2011. Memoria técnica de la cuantifi cación de la deforestación histórica nacional–escalas gruesa y fi na. Instituto de Hidrología, Meteorología, y Estudios Ambientales-IDEAM-. Bogotá D.C., Colombia. 106 p.
  17. 17. Montoya R, Chaparro PE, Araújo RM, Narváez PC, Valencia SH. Epidemiology and control of malaria in Epidemiology and control of malaria in Colombia. 2011; 1:114–22.
  18. 18. González R, Carrejo N. Introducción al estudio taxonómico de Anopheles de Colombia: claves y notas de distribución. 2nd ed. Cali: Universidad del Valle; 2009. 260p.
  19. 19. Cienfuegos A, Gómez G, Córdoba L, Luckhart S, Conn J, Correa M. Diseño y evaluación de metodologías basadas en PCR–RFLP de ITS2 para la identificación molecular de mosquitos Anopheles spp. (Diptera:Culicidae) de la Costa Pacífica de Colombia. Rev Biomed. 2008;19: 35–44.
  20. 20. Cienfuegos A V, Rosero DA, Naranjo N, Luckhart S, Conn JE, Correa MM. Acta Tropica Evaluation of a PCR–RFLP–ITS2 assay for discrimination of Anopheles species in northern and western Colombia. Acta Trop.; 2004;118: 128–135.
  21. 21. Zapata MA, Cienfuegos A V, Quirós OI, Quiñones ML, Luckhart S, Correa MM. Discrimination of seven Anopheles species from San Pedro de Uraba, Antioquia, Colombia, by polymerase chain reaction-restriction fragment length polymorphism analysis of its sequences. Am J Trop Med Hyg. 2007;77: 67–72. pmid:17620632
  22. 22. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3: 294–9. pmid:7881515
  23. 23. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. Proc Biol Sci. 2003;270: 313–21. pmid:12614582
  24. 24. Gómez G, Jaramillo L, Correa MM. Wing geometric morphometrics and molecular assessment of members in the Albitarsis Complex from Colombia. Mol Ecol Resour. 2013;13: 1082–92. pmid:23702155
  25. 25. Lourenço PM, Sousa CA, Seixas J, Lopes P, Novo MT, Almeida APG. Anopheles atroparvus density modeling using MODIS NDVI in a former malarious area in Portugal. J Vector Ecol. 2011;36: 279–91. pmid:22129399
  26. 26. Dwyer J, Schmidt G. The MODIS Reprojection Tool. Earth science satellite remote sensing. Berlin, Heidelberg: Springer Berlin Heidelberg; 2006. pp. 162–177. https://doi.org/10.1007/978-3-540-37294-3_9
  27. 27. Feilhauer H, He KS, Rocchini D. Modeling species distribution using niche-based proxies derived from composite bioclimatic variables and MODIS NDVI. Remote Sens; 2012;4: 2057–2075.
  28. 28. Dantur Juri MJ, Estallo E, Almirón W, Santana M, Sartor P, Lamfri M, et al. Satellite-derived NDVI, LST, and climatic factors driving the distribution and abundance of Anopheles mosquitoes in a former malarious area in northwest Argentina. J Vector Ecol. 2015;40: 36–45. pmid:26047182
  29. 29. Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Modell. 2006;190: 231–259.
  30. 30. Merow C, Allen JM, Aiello-Lammens M, Silander JA. Improving niche and range estimates with Maxent and point process models by integrating spatially explicit information. Glob Ecol Biogeogr. 2016;25: 1022–1036.
  31. 31. Alimi TO, Fuller DO, Qualls WA, Herrera S V., Arévalo-Herrera M, Quiñones ML, et al. Predicting potential ranges of primary malaria vectors and malaria in northern South America based on projected changes in climate, land cover and human population. Parasit Vectors. BioMed Central; 2015;8: 431. pmid:26289677
  32. 32. Fuller DO, Ahumada ML, Quiñones ML, Herrera S, Beier JC. Near-present and future distribution of Anopheles albimanus in Mesoamerica and the Caribbean Basin modeled with climate and topographic data. Int J Health Geogr. BioMed Central; 2012;11: 13. pmid:22545756
  33. 33. Barve N, Barve V, Jiménez-Valverde A, Lira-Noriega A, Maher SP, Peterson AT, et al. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol Modell; 2011;222: 1810–1819.
  34. 34. Pearson RG, Raxworthy CJ, Nakamura M, Peterson AT. Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J Biogeogr; 2006;34: 102–117.
  35. 35. Gómez-Palacio A, Arboleda S, Dumonteil E, Peterson AT. Ecological niche and geographic distribution of the Chagas disease vector, Triatoma dimidiata (Reduviidae: Triatominae): Evidence for niche differentiation among cryptic species. Infect Genet Evol. 2015;36: 15–22. pmid:26321302
  36. 36. Peterson AT, Papes M, Sober J. Rethinking receiver operating characteristic analysis applications in ecological niche modeling. Ecol Modell. 2007;3: 63–72.
  37. 37. Idohou R, Peterson AT, Assogbadjo AE, Vihotogbe RL, Padonou E, Glèlè Kakaï R. Identification of potential areas for wild palm cultivation in the Republic of Benin through remote sensing and ecological niche modeling. Genet Resour Crop Evol. 2016; 1–11.
  38. 38. Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution; 2008;62: 2868–2883. pmid:18752605
  39. 39. Mota-Vargas C, Rojas-Soto OR. Taxonomy and ecological niche modeling: Implications for the conservation of wood partridges (genus Dendrortyx). J Nat Conserv. 2016;29: 1–13.
  40. 40. Qiao H, Peterson AT, Campbell LP, Soberón J, Ji L, Escobar LE. NicheA: Creating virtual species and ecological niches in multivariate environmental scenarios. Ecography. 2016;
  41. 41. Qiao Huijie, Escobar Luis E., Peterson AT. Accessible areas in ecological niche comparisons of invasive species: Recognized but still overlooked. Sci Rep. 2017;2: 1–24.
  42. 42. Parra-Henao G, Alarcón Pineda EP. Observaciones sobre la bionomía de Anopheles spp. (Diptera: Culicidae) en el municipio Valencia, departamento Córdoba, Colombia. Boletín Malariol y Salud Ambient. 2008;48: 95–98.
  43. 43. Alimi TO, Fuller DO, Quinones ML, Xue R-DD, Herrera S V, Arevalo-Herrera M, et al. Prospects and recommendations for risk mapping to improve strategies for effective malaria vector control interventions in Latin America. Malar J. BioMed Central; 2015;14: 519. pmid:26694047
  44. 44. Tadei WP, Thatcher BD, Santos JM, Scarpassa VM, Rodrigues IB, Rafael MS. Ecologic observations on anopheline vectors of malaria in the Brazilian Amazon. Am J Trop Med Hyg. 1998;59: 325–35. pmid:9715956
  45. 45. Ahumada ML, Pareja PX, Buitrago LS, Quiñones ML. Biting behavior of Anopheles darlingi Root, 1926 (Diptera: Culicidae) and its association with malaria transmission in Villavicencio (Meta, Colombia). Biomedica. 2013;33: 241–50. pmid:24652134
  46. 46. Brochero HL, Rey G, Buitrago LS, Olano VA. Biting activity and breeding sites of Anopheles species in the municipality Villavicencio, Meta, Colombia. J Am Mosq Control Assoc. 2005;21: 182–6. pmid:16033120
  47. 47. Manguin S, Roberts DR, Andre RG, Rejmankova E, Hakre S. Characterization of Anopheles darlingi (Diptera: Culicidae) larval habitats in Belize, Central America. J Med Entomol. 1996;33: 180–194.
  48. 48. Villarreal-Treviño C, Penilla-Navarro RP, Vázquez-Martínez MG, Moo-Llanes DA, Ríos-Delgado JC, Fernández-Salas I, et al. Larval habitat characterization of Anopheles darlingi from its northernmost geographical distribution in Chiapas, Mexico. Malar J. BioMed Central; 2015;14: 517. pmid:26693816
  49. 49. Vittor AY, Pan W, Gilman RH, Tielsch J, Glass G, Shields T, et al. Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009;81: 5–12. pmid:19556558
  50. 50. Roberts DR, Manguin S, Rejmankova E, Andre R, Harbach RE, Vanzie E, et al. Spatial distribution of adult Anopheles darlingi and Anopheles albimanus in relation to riparian habitats in Belize, Central America. J Vector Ecol. 2002;27: 21–30. pmid:12125869
  51. 51. Betancur T, Mejia O, Palacio C. Conceptual hydrogeology model to Bajo Cauca antioqueño: a tropical aquifer system. Rev Fac Ing. 2009; 107–118.
  52. 52. Vittor AY, Gilman RH, Tielsch J, Glass G, Shields T, Lozano WS, et al. The effect of deforestation on the human-biting rate of Anopheles darlingi, the primary vector of falciparum malaria in the Peruvian Amazon. Am J Trop Med Hyg. 2006;74: 3–11. pmid:16407338
  53. 53. Lainhart W, Bickersmith SA, Nadler KJ, Moreno M, Saavedra MP, Chu VM, et al. Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru. Malar J. BioMed Central; 2015;14: 375. pmid:26415942
  54. 54. Magris M, Rubio-Palis Y, Menares C, Villegas L. Vector bionomics and malaria transmission in the Upper Orinoco River, southern Venezuela. Mem Inst Oswaldo Cruz; 2007;102: 303–312. pmid:17568935
  55. 55. de Castro MC, Monte-Mór RL, Sawyer DO, Singer BH. Malaria risk on the Amazon frontier. Proc Natl Acad Sci U S A. National Academy of Sciences; 2006;103: 2452–7. pmid:16461902
  56. 56. Bejarano EE. Occurrence of the malaria Vector Anopheles albimanus Wiedemman (Diptera: Culicidae) in Isla Fuerte, Colombia. Neotrop Entomol; 2003;32: 517–518.
  57. 57. Gómez GF, Márquez EJ, Gutiérrez LA, Conn JE, Correa MM. Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation. Acta Trop; 2014;135: 75–85. pmid:24704285
  58. 58. Stefani A, Dusfour I, Corrêa APS, Cruz MC, Dessay N, Galardo AK, et al. Land cover, land use and malaria in the Amazon: a systematic literature review of studies using remotely sensed data. Malar J. 2013;12: 192. pmid:23758827
  59. 59. Turvey CG, Mclaurin MK, Turvey CG, Mclaurin MK. Applicability of the normalized difference vegetation index (NDVI) in index-based crop insurance design. Weather Clim Soc. 2012;4: 271–284.
  60. 60. Pinault LL, Hunter FF. Characterization of larval habitats of Anopheles albimanus, Anopheles pseudopunctipennis, Anopheles punctimacula, and Anopheles oswaldoi s.l. populations in lowland and highland Ecuador. J Vector Ecol. 2012;37: 124–136. pmid:22548546
  61. 61. Faran ME. Mosquito studies Diptera: Culicidae. A revision of the albimanus section of the subgenus Nyssorhynchus of Anopheles. Contrib Am Entomol Inst. 1980;15: 1–215.
  62. 62. Hiwat H, Bretas G. Ecology of Anopheles darlingi Root with respect to vector importance: a review. Parasit Vectors; 2011;4: 177. pmid:21923902
  63. 63. Loaiza JR, Scott ME, Bermingham E, Sanjur OI, Wilkerson R, Rovira J, et al. Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics. Mol Phylogenet Evol. 2010;57: 1341–6. pmid:20888924
  64. 64. Sinka ME, Bangs MJ, Manguin S, Rubio-Palis Y, Chareonviriyaphap T, Coetzee M, et al. A global map of dominant malaria vectors. Parasit Vectors. 2012;5: 69. pmid:22475528
  65. 65. Peterson AT, Bauer JT, Mills JN. Ecologic and geographic distribution of filovirus disease. Emerg Infect Dis; 2004;10: 40–7. pmid:15078595
  66. 66. Lash RR, Brunsell N, Peterson AT. Spatiotemporal environmental triggers of Ebola and Marburg virus transmission. Geocarto Int. 2008;23: 451–466.
  67. 67. Afrane YA, Githeko AK, Yan G. The ecology of Anopheles mosquitoes under climate change: case studies from the effects of deforestation in East African highlands. Ann Acad Sci. 2012;1249: 204–210. pmid:22320421
  68. 68. Fuller DO, Troyo A, Alimi TO, Beier JC. Participatory risk mapping of malaria vector exposure in northern South America using environmental and population data. Appl Geogr. 2014;48: 1–7. pmid:24976656