SciELO - Scientific Electronic Library Online

 
vol.62 issue4Diversity and potential distribution of culicids of medical importance of the Yucatan Peninsula, MexicoDetermination of insecticides’ lethal concentrations and metabolic enzyme levels in Triatoma dimidiata author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Salud Pública de México

Print version ISSN 0036-3634

Salud pública Méx vol.62 n.4 Cuernavaca Jul./Aug. 2020  Epub May 16, 2022

https://doi.org/10.21149/10111 

Artículos originales

Composition and abundance of anopheline species according to habitat diversity in Mexico

Composición y abundancia de especies de anofelinos según la diversidad de hábitats en México

Cuauhtémoc Villarreal-Treviño(1) 

Jana C Ríos-Delgado(1) 

R Patricia Penilla-Navarro(1) 

Américo D Rodríguez(1) 

José Higinio López(2) 

José A Nettel-Cruz(1) 

David A Moo-Llanes(1  3) 

Gabriel Fuentes-Maldonado(1) 

(1) Centro Regional de Investigación en Salud Pública, Instituto Nacional de Salud Pública. Tapachula, Chiapas, Mexico.

(2) Laboratorio de Análisis de Información Geográfica y Estadística, El Colegio de la Frontera Sur. Tapachula, Chiapas, Mexico.

(3) Facultad de Medicina, Universidad Nacional Autónoma de México. Mexico City, Mexico.


Abstract:

Objective:

To determine the abundance and geographic distribution of the main malaria vectors, which are influenced by habitat characteristics and ecological factors that directly impact adult density and the dynamics of malaria transmission in Mexico.

Materials and methods:

Samples of larvae were collected from 19 states in Mexico. Each larval habitat was characterized in situ determining the following parameters: water depth, turbidity, percentage of vegetation cover, amount of detritus, presence of algae, light intensity, type of vegetation, amount of predators, habitat stability, altitude, and hydrologic type.

Results:

A total of 21 687 larvae corresponding to 13 anopheline species were obtained from 149 aquatic habitats. The most abundant species were Anopheles pseudopunctipennis (52.91%), An. albimanus (39.14%) and An. franciscanus (5.29%). The multiple logistic regression analysis showed a negative association between An. pseudopunctipennis and water turbidity (ß=-1.342; Wald=6.122; p=0.013) and the amount of detritus (ß=-2.206; Wald=3.642; p=0.050). While in An. albimanus, there was a significant positive association with water turbidity (ß=1.344; Wald=4.256; p=0.039), a negative correlation was found with the altitude (ß=-3.445; Wald=5.407; p=0.020). The highest mosquito species diversity index was found in Chiapas (Fisher’s α=1.20) and the lowest diversity in Chihuahua (Fisher’s α=0.26). The greatest richness was found in streams (n=11).

Conclusions:

The two most abundant species were: An. albimanus and An. pseudopunctipennis. Detailed knowledge of the distribution and characteristics of their larval habitats will be useful for the effective implementation of control strategies in Mexico.

Keywords: Anopheles; malaria; larval habitat; Mexico

Resumen

Objetivo:

Determinar la abundancia y la distribución geográfica de los principales vectores de la malaria, las cuales están influenciadas por las características del hábitat y los factores ecológicos que afectan directamente la densidad de los adultos y la dinámica de la transmisión de la malaria en México.

Material y métodos:

Se obtuvieron muestras de larvas de 19 estados de México. Cada hábitat larvario se caracterizó in situ determinando los siguientes parámetros: profundidad del agua, turbidez, porcentaje de cobertura vegetal, cantidad de detritus, presencia de algas, intensidad de luz, tipo de vegetación, cantidad de depredadores, estabilidad del hábitat, altitud y tipo hidrológico.

Resultados:

Se identificaron un total de 21 687 larvas pertenecientes a 13 especies de anofelinos, de 149 hábitats acuáticos. Las tres especies más abundantes fueron Anopheles pseudopunctipennis (52.91%), An. albimanus (39.14%) y An. franciscanus (5.29%). El análisis de regresión logística múltiple mostró una asociación negativa para An. pseudopunctipennis y la turbidez del agua (ß=-1.342; Wald= 6.122; p=0.013) y la cantidad de detritus (ß=-2.206; Wald= 3.642; p=0.050). Para An. albimanus se encontró una asociación positiva significativa con la turbidez del agua (ß=1.344; Wald= 4.256; p=0.039) y una correlación negativa con la altitud (ß=-3.445; Wald=5.407; p=0.020). El índice de diversidad más alto se encontró en Chiapas (α de Fisher=1.20) y la diversidad más baja en Chihuahua (α de Fisher=0.26). La mayor riqueza se encontró en los arroyos (n=11).

Conclusiones:

Las dos especies más abundantes fueron An. albimanus y An. pseudopunctipennis. El conocimiento detallado de la distribución y características de sus hábitats larvales será útil para la implementación efectiva de las estrategias de control en México.

Palabras clave: Anopheles; malaria; hábitat larvario; México

Malaria is the most important anopheles-borne disease and causes great mortality worldwide. This disease is endemic to 106 countries, primarily in the tropical and subtropical regions. In 2015, there were approximately 212 million new cases and 429 000 deaths, mostly among children. Currently, half of the world population, or about 3.2 billion people, are at risk of contracting malaria.1 In Mexico, the number of cases have decreased significantly, from 7 259 in 2000 to 518 in 2016.2,3 Malaria cases are concentrated in four foci: a) Sinaloa and Nayarit, b) Durango and Chihuahua, c) Campeche and Quintana Roo, and d) Chiapas (Lacandon Forest). The latter is the main focus of transmission in the country, having registered 63.70% of the cases in 2016.3

There are 475 species of the genus Anopheles in the world, 41 of which have been incriminated for transmitting the Plasmodium parasite.4,5 In Mexico, there are 26 species of Anopheles6,7 distributed along the country, and three were incriminated as the main vectors of malaria: Anopheles pseudopunctipennis, An. albimanus and An. vestitipennis.8,9,10,11 Studies of the larval ecology of anophelines in Mexico date back to the fifties.6 However, the geographic distribution and abundance of the species have changed due to human activities such as deforestation, increase of croplands and cattle grazing areas, road building, creation of new human settlements, and water control systems, among others.12,13

Recent studies on the anopheline larvae ecology were carried out mainly in southeastern Mexico (Chiapas), pointing to An. pseudopunctipennis as the main vector in the mountain regions. Their larvae inbreed in puddles along river margins, with a positive association with the presence of filamentous green algae, and the presence of Heteranthera.14,15 While An. albimanus is the main malaria vector in the coastal plain of the Pacific and Atlantic. 6,8,16 In Chiapas, larvae of this species are positively associated to planktonic algae both in the dry and wet seasons, and negatively associated with altitude.15 During the wet season the larvae can be found in the margins of permanent water bodies, positively associated with floating plants, and during the rainy season they are positively correlated with the presence of emergent plants, particularly Cyperaceae and flood grasses.15,17 The main vector in the Lacandon Forest (Chiapas) is An. vestitipennis and its larval habitat has not been described in Mexico.6,9 Recently, the larval habitats of An. darlingi, the secondary vector in the Lacandon Forest, were studied and characterized.18 Part of the knowledge obtained in those studies has been the basis for the implementation of successful antilarval measures in Chiapas and Oaxaca, applying the technique of “Elimination of Anopheline Breeding Sites”, which consists mainly in removing green algae from rivers through community participation.19 This strategy reduced the density of An. pseudopunctipennis larvae and impacted the adult population in such a way that the numbers of malaria cases were significantly reduced.20 However, this knowledge cannot be applied equally throughout the Mexican territory because of the high environmental variability caused by its biogeography, because the northern part of the country is located in the Nearctic region, and the southern part, in the Neotropical region.21 Each biogeographical region has specific climatic, hydrologic and orographic characteristics that are unique to the larval habitats of each anopheline species.22 It is essential to know the ecology and type of larval habitats of anophelines in order to design control measures of the immature stadia using larvicides, environmental management or the removal of filamentous algae.

Materials and methods

Study área

A descriptive cross-sectional research was carried out in 19 states in Mexico from April to February 2016. The study area was selected based on previous records on Mexican anophelines.6 The area covers two main biotic regions: Nearctic and Neotropical. Ten states are located in the Nearctic region and nine in the Neotropical region, which include a wide variety of climates, topographic diversity, different land uses and vegetation. The climatic conditions in the Nearctic zone are very varied, but the predominant climates are arid and semi-arid (Koppen BW, BS), with precipitation ranging between 100 and 600 mm, and an annual median temperature ranging between 22°C and 26°C in some regions, and 18°C and 22°C in the cold and temperate areas. The Neotropical region comprises the states of southern Mexico. The climate is hot and humid with dry and rainy seasons (Koppen Aw) or tropical with year-round rains (Koppen Af), a median annual temperature of 18°C to 22°C, and precipitation ranges between 2 000 and 4 500 mm.23,24 The search for anopheline larvae was carried out at any available body of water found along the trip both inside and outside urban or rural settlements, by the side of the road, in natural and man-made habitats.

Larvae collection

Each site was visited every two to three months (from February 2012 to April 2016). The larval habitats were sampled along their margin using a 500 mL dipper, according to a standardized method.17 The water samples were poured in a white 30×30×5 cm tray in order to collect larvae of all stadia (I to IV). The larvae were counted and placed in 1 L plastic containers to be transported to the insectary (Centro Regional de Investigación en Salud Pública [CRISP]) for rearing to adulthood. The samples were collected in each habitat for at least 1 hr., between 8:00 and 14:00 hours. Google Earth was utilized in the search for water body images, and the positive sites were georeferenced using a GPS (Garmin eTrex Garmin Ltd., Kansas USA).

Larvae identification and growth

The larvae were raised in the CRISP insectary in a climate-controlled room at 25±1 °C temperature, 70±10% relative humidity and a photoperiod of 12:12 (L: D) h, using a standardized technique.25 In some cases where the distance did not allow the transportation of biological material to CRISP, an insectary was improvised in the study area under natural conditions of temperature, humidity and photoperiod. All the larvae were raised using water from the same collection site in order to minimize mortality, but taking care of using clean water free of predatory insects.26 The larvae were fed with finely ground and sterilized commercially available mouse food (Laboratory Rodent Diet 5001, LabDiet, St. Louis, MO). All the adults obtained were mounted and accounted for in a database built for later analysis. Larval densities (LD) were quantified as total number of larvae/the number of dips.18 All specimens were identified using the Wilkerson and Strickman (1990) keys, based on morphological characteristics.7

Characterization of larval hábitats

Several water bodies were visited, selecting those that were positive and negative for anopheline larvae. Each larval habitat was characterized in situ determining the following parameters: water depth, turbidity, plant cover (%), amount of detritus, presence of algae, light intensity, vegetation type, amount of predators, water movement, habitat stability, altitude, and hydrologic type. Water depth was measured using a graduated string with a lead weight, and was categorized in three strata: 0-200 cm, 201-400 cm, and 401-600 cm. Turbidity was determined according to color as “clear or turbid”. Plant cover was classified in the following groups: 0-30%, 31-60%, and 61-100%. Detritus, composed of twigs, dry leaves and dead insects, was determined according to presence in the habitat as “low, moderate or high”. The algae were characterized as “present or absent”. Light intensity was measured according to the amount of vegetation surrounding the habitat and categorized as “sunny, partially sunny, or shaded”. Vegetation was categorized as “emergent, floating/submerged, or none”. Emerging plants include aquatic vegetation and submerged terrestrial vegetation. Predators were categorized by larval habitat according to abundance: low (≤5 predators), moderate (6 to 10 predators) or abundant (≥11 predators). Water movement was determined as “stagnant, slow, or fast”. Habitat stability was considered permanent if the water body persisted throughout the year and temporary if it was only present during part of the year. Altitude was determined by portable GPS as follows: 0-759 meters above the sea level (masl), 760-1 519 masl and 1 520-2 279 masl. The hydrologic type was classified according to the nature of the water body.18

Data analysis

A one-way analysis of variance was used (Anova) to compare the absolute larval index (ALI) between hydrologic types, and for species abundance, the results from the identification of the obtained adults were used. A multiple logistic regression analysis was used to determine the association between environmental variables and the occurrence of An. pseudopunctipennis and An. albimanus larvae. The presence of larvae was categorized as one (1) and the absence of larvae was categorized as zero (0). Data were analyzed using statistical software SPSS 20.0 (SPSS, Inc., Chicago, IL). We determined the species richness with Fisher’s α test (number of species found in each state and type of habitat), using the software “Species Diversity & Richness version 4”.27,28 Geographic distribution maps of the species were created utilizing the QGIS software Version 2.18.1.

Results

Composition and abundance of anopheline species

In this study, 19 states of the Mexican Republic were examined, selecting 169 larval habitats: 91.94% were natural, and the rest were artificial; 149 tested positive for anopheline larvae, and 20 larval habitats tested negative (figure 1). All water bodies, positive and negative, were sampled throughout the study. Lagoons and fast flowing rivers were the main negative habitats. Of the total positive larval habitats, 21 were located in the Nearctic region and 128 in the Neotropical region. In the Nearctic region of Mexico, five anopheline species were found: An. pseudopunctipennis, An. albimanus, An. franciscanus, An. punctipennis and An. crucians (figure 2). Of these species, An. franciscanus and An. punctipennis were specific to this biogeographic region; the rest of the species were found across Mexico. In central Mexico, the transition between the Nearctic and Neotropical regions (the states of Morelos and Puebla), only An. pseudopunctipennis was detected, in river margins. Eleven species were found in southern Mexico: An. pseudopunctipennis, An. albimanus, An. vestitipennis, An. darlingi, An. punctimacula, An. crucians, An. hectoris, An. apicimacula, An. gabaldoni, An. argyritarsis, and An. eiseni. The southern state of Chiapas presented 112 aquatic habitats and 11 mosquito species, in the Lacandon Forest area and the Soconusco area in the coastal plain (figure 2).

A total 21 687 larvae of 13 Anopheles species were captured and reared until adulthood with 84.48% of breeding success, i.e. obtaining 18 322 adults. The most abundant species throughout the study was An. pseudopunctipennis (52.92%), followed by An. albimanus (39.14%) and An. franciscanus (5.29%), which comprise 97.35% of the total specimens. The remaining 2.54 % is constituted by the least abundant species: An. punctipennis (1.0%), An. vestitipennis (0.70%), An. darlingi (0.49%), An. Punctimacula (0.13%), An. crucians (0.11%), An. hectoris (0.11%), An. apicimacula (0.08%), An. gabaldoni (0.03%), An. argyritarsis (0.01%), and An. eiseni (0.01%) (figure 3). An. pseudopunctipennis had a significantly higher ALI (±StD) (14.05±21.11), compared to other species (F=1.897; df=11; p< 0.05) (table I). As for the spatial distribution of the main malaria vectors in Mexico, An. albimanus was found in 13 states, An. pseudopunctipennis in 10 states and An. vestitipennis was found only in two states in southern Mexico (figure 4).

Figure 1 Study area indicating the positive and negative sampling sites of anopheline larvae. February 2010 to April 2015. Mexico 

Figure 2 Spatial distribution of the occurrence points of Anopheles species in the Nearctic (A) and Neotropical (B) areas of Mexico. February 2010 to April 2015 

Figure 3 Percentage of anopheline species in Mexico. February 2010 to April 2015 

Figure 4 Spatial distribution of the occurrence points of the main malaria vectors in Mexico. a) An. pseudopunctipennis b) An. albimanus, c) An. vestitipennis. February 2010 to April 2015 

Habitat diversity and larval abundance

Twelve hydrologic types were found positive for anopheline larvae in Mexico. The most common hydrologic types were river margins (29.53%), rain puddles (16.77%), ponds (12.75%) and river pools (10.73%). The least frequent were swamps (6.71%), streams (6.71%), ditches (4.02%), pools (4.02%), lagoons (3.35%), gravel pits (2.68%), irrigation canals (2.01%), and drinking troughs (0.67%). The average ALI (±StD) was 14.55 larvae per dip (l/d) (±25.28) with a range of 0 to 153.0 l/d. The greatest ALI averages (±StD) were found in river pools 19.40 l/d (±19.46), ponds 18.99 l/d (±35.75), river margins 18.97 l/d (±31.19), gravel pits 17.65 l/d (±9.13), streams 17.20 l/d (±27.41) and ditches 16.25 l/d (±28.89). No significant differences were found in the absolute larval index (ALI) between hydrologic types (F=0.748; df=11,137; p=0.690) (table I).

Larval presence and associated parameters in An. Pseudopunctipennis and An. Albimanus

The multiple logistic regression analysis showed a significant negative association between the presence of An. pseudopunctipennis and water turbidity (ß=-1.342; 232 Wald=6.122; p=0.013) and the amount of detritus (ß=-2.206; Wald=3.642; p=0.050) (table II). In contrast An. albimanus had a significant positive association with water turbidity (ß=1.344; Wald=4.256; p=0.039) and a negative correlation with altitude (ß= −235 3.445; Wald=5.407; p=0.020) (table III).

Richness and diversity of anopheline species

The diversity of Anopheles species in 19 states exhibited a general Fisher’s α value of 1.77. The state with the most species richness was Chiapas (n=11), followed by Coahuila (n=5). In other states, only 1-3 species were detected. The greatest diversity index was found in Chiapas (Fisher’s α=1.20), and the lowest diversity index was found in Chihuahua (Fisher’s α=0.26) (table IV). The greatest richness was found in streams (n=11), followed by river margins (n=9) and ponds (n=6).

Table I: Average adult density (±StD) of anophelines obtained by hydrologic type, in 19 States of Mexico. February/2010 to April/2015 

Hydrologic type n Ap Aa Af Apu Av Ad Apc Ac Ah Api Ag Aar Ae
River margin 56 18.35±27.43 9.65±16.80 4.83±0.28 5.00 0 0.01±0.008 0.55±0.07 0.55±0.63 2.00 0.30±0.28 0.50 0.10 0.10
Rain puddle 25 1.20±0.90 6.87±17.65 0 0 2.10 0.07±0.03 0.10 0 0 0 0 0 0
Pond 19 23.62±28.53 11.91±16.60 0 0 0.65±0.77 0.10±0.06 0.50 0.10 0 0 0 0 0
River pool 16 16.93±14.81 4.80±5.18 2.00 5.00 0 0 0 0 0 0 0 0 0
Swamp 10 0.53±0.58 5.76±2.51 0 0 0 0 0 0.40 0 0.80 0 0 0
Stream 10 10.46±11.65 3.93±3.61 2.30±1.85 2.76±0.49 1.00 0.18±0.15 0 0.40 0 0 0 0 0
Ditch 6 0.87±1.03 3.25±1.06 70.00 0 2.33±1.75 0.25±0.17 0 0 0 0 0 0 0
Ground pool 6 1.70±0.70 5.44±3.88 0 0 0 0.21±0.08 0 0 0 0 0 0 0
Lagoon 14 0 2.85±3.04 0 0 0 0 0 0 0 0 0 0 0
Gravel pit 4 14.25±6.35 2.95±2.41 0 0 0 0 0 0 0 0 0 0 0
Irrigation canal 3 0 11.6±16.60 0 0 0.75±0.35 0 0.60 0 0 0 0 0 0
Drinking trough 1 0 0 1.20 0 0 0 0 0 0 0 0 0 0
Total 149 14.05±21.11 7.54±13.45 9.69±21.26 3.66±1.27 1.43±1.23 0.47±0.36 0.46±0.20 0.40±0.36 2.00 0.46±0.35 0.50 0.10 0.10

Ap=An. pseudopunctipennis; Aa=An. albimanus, Af=An. franciscanus, Apu=An. punctipennis, Av=An. vestitipennis, Ad=An. darlingi, Apc=An. punctimacula, Ac= An. crucians, Ah=An. hectoris, Api=An. apicimacula, Ag=An. gabaldoni, Aar=An. argyritarsis, Ae=An. eiseni.

Table II: Multiple logistic regression analysis of the larval abundance and environmental parameters for Anopheles pseudopunctipennis in Mexico. February 2010 to April 2015 

Variables Coefficient(ß) S.E. Wald df Sig Exp(ß) CI95% Exp(ß)
Lower limit Upper limit
Interception 1.126 2.557 0.194 1 0.660
Water depth
0-200 0 0
201-400 -1.156 0.000 1 0.315 0.315 0.315
401-800 0 0
Turbidity
Clear -1.342 0.542 6.122 1 0.013 0.261 0.090 0.757
Turbid 0 0
Coverage %
0-30 -0.564 0.808 0.487 1 0.485 0.569 0.117 2.771
31-60 -0.592 0.945 0.393 1 0.531 0.553 0.087 3.525
61-100 0 0
Detritus
Low -1.148 0.692 2.756 1 0.097 0.317 0.082 1.230
Moderate 0 0
High -2.206 1.156 3.642 1 0.050 0.110 0.011 1.061
Algae
Absent 0.656 0.675 0.943 1 0.331 1.927
Present 0 0
Light intensity
Sunny -0.694 0.892 0.605 1 0.437 0.500
Sun/shade -1.320 0.976 1.829 1 0.176 0.267
Shade 0 0
Vegetation
Emergent -0.007 1.362 0.000 1 0.996 0.993 0.069 14.325
Floating/submerged 0.385 1.396 0.076 1 0.783 1.469 0.095 22.675
None 0.119 1.471 0.007 1 0.935 1.127 0.063 20.137
Predators
Low 0.452 0.539 0.705 1 0.401 1.572 0.547 4.517
Moderate 0 0
Abundant -1.121 1.315 0.727 1 0.394 0.326 0.025 4.291
Water movement
Stagnant 0 0
Slow 0 0
Fast 0.052 1.165 0.002 1 0.964 1.053 0.107 10.344
Habitat stability
Permanent 0.417 0.706 0.348 1 0.555 1.517 0.380 6.054
Temporary 0 0
Altitude
0 -759 1.264 1.018 1.541 1 0.214 3.540 0.481 26.056
760 -1 519 0
1 520 -2 279 0.926 1.558 0.353 1 0.552 2.524 0.119 53.487
Hydrologic type
Ditch -1.439 1.488 0.935 1 0.333 0.237 0.013 4.381
Drinking trough 20.230 0.000 1 6.104E-5 6.104E-5 6.104E-5
Gravel pit -19.127 7 132.062 0.000 1 0.998 4.934E-9 0.000
Pools -1.223 1.488 0.675 1 0.411 0.294 0.016 5.439
Irrigation canal 18.154 8 502.676 0.000 1 0.998 76 627 495.02 0.000
Lagoon 18.715 6 306.399 0.000 1 0.998 134 229 308.46 0.000
Pond 0.663 1.232 0.290 1 0.590 1.941 0.174 21.703
Rain puddle 1.471 1.299 1.282 1 0.257 4.355 0.341 55.594
River margins -0.752 1.031 0.532 1 0.466 0.472 0.063 3.555
River pools -1.395 1.242 1.262 1 0.261 0.248 0.022 2.826
Creeks -0.408 1.304 0.098 1 0.755 0.665 0.052 8.572
Swamps 0 0

Sig: significant; df: degrees of freedom

Table III: Multiple logistic regression analysis of larval abundance and environmental parameters for Anopheles albimanus in Mexico. February 2010 to April 2015 

Variables Coefficientß S.E. Wald df Sig Exp(ß) CI95% Exp(ß)
Lower limit Upper limit
Interception 1.317 2.580 0.261 1 0.610
Water depth
0-200 0.054 0.797 0.005 1 0.946 1.056 0.221 5.040
201-400 0 0
401-600 0 0
Turbidity
Clear 1.344 0.652 4.256 1 0.039 3.836 1.069 13.760
Turbid 0 0
Coverage %
0-30 -0.005 0.922 0.000 1 0.995 0.995 0.163 6.058
31-60 -0.703 1.229 0.327 1 0.568 0.495 0.045 5.509
61-100 0 0
Detritus
Low 1.528 0.829 3.394 1 0.065 4.609 0.907 23.421
Moderate 0 0
High 1.421 1.287 1.219 1 0.270 4.140 0.333 51.534
Algae
Absent -0.544 0.835 0.424 1 0.515 0.581 0.113 2.985
Present 0 0
Light intensity
Sunny 0.225 0.937 0.058 1 0.810 1.252 0.199 7.865
Sun/shade 0.843 0.966 0.762 1 0.383 2.324 0.350 15.422
Shade 0 0
Vegetation
Emergent 0.234 0.694 0.114 1 0.736 1.264 0.324 4.923
Floating/submerged 0 0
None 0.651 0.944 0.475 1 0.491 1.917 0.301 12.204
Predators
Low -0.157 0.637 0.061 1 0.805 0.854 0.245 2.980
Moderate 0 0
Abundant -16.759 7 018.994 .0000 1 0.998 5.27E8 0.000
Water movement
Stagnant 0 0
Slow 0 0
Fast 1.886 1.335 1.996 1 0.158 6.594 0.482 90.252
Habitat stability
Permanent -1.456 0.828 3.092 1 0.079 0.233 0.046 1.182
Temporary 0 0
Altitude
0 - 759 -3.445 1.486 5.407 1 0.020 0.032 0.002 0.581
760 - 1 519 0 0
1 520 - 2 279 15.470 0.000 1 5 230 658.306 5 230 658.306 5 230 658.306
Hydrologic type
Ditch 1.531 1.993 0.590 1 0.443 4.621 0.093 229.885
Drinking trough 34.453 0.000 1 9.180E-14 9.180E-14 9.180E-14
Gravel pit -19.824 7 540.604 0.000 1 0.998 2.457E-9 0.000
Pools 0.969 1.935 0.251 1 0.617 2.635 0.059 116.950
Irrigation canal -17.419 9 242.478 0.000 1 0.998 2.722E-8 0.000
Lagoon 2.077 1.674 1.540 1 0.215 7.984 0.300 212.396
Pond -1.091 1.623 0.452 1 0.501 0.336 0.014 8.079
Rain puddle -3.157 1.936 2.660 1 0.103 0.043 0.001 1.891
River margins 0.808 1.380 0.343 1 0.558 2.244 0.150 33.555
River pools 0.791 1.505 0.276 1 0.599 2.205 0.115 42.124
Creeks -0.282 1.711 0.027 1 0.869 0.754 0.026 21.568
Swamps 0 0

Sig: significant; df: degrees of freedom

Table IV: Ecologic diversity of anophelines in 19 states in Mexico. February 2010 to April 2015 

Geographic region/State Number of habitats Number of larvae Species richness Fisher’s α
Nearctic region
B. California 1 10 1 0.71
Sinaloa 2 1 160 2 0.38
Nayarit 1 25 1 0.50
Chihuahua 1 750 1 0.26
Coahuila 3 959 5 1.00
Durango 1 12 1 0.65
Nuevo León 6 4 245 3 0.42
Tamaulipas 1 60 2 0.65
Puebla 4 285 1 0.33
Morelos 1 65 1 0.39
Neotropical region
Jalisco 1 55 1 0.50
Colima 1 32 1 0.51
Guerrero 1 56 1 0.41
Oaxaca 3 987 2 0.38
Chiapas 112 11 667 11 1.20
Veracruz 6 816 3 0.59
Tabasco 1 88 1 0.50
Campeche 1 65 1 0.41
Quintana Roo 2 350 2 0.45
Total 149 21 687 13 1.77

Discussion

There were 13 Anopheles species identified throughout Mexico in this study. An. pseudopunctipennis and An. albimanus were the most abundant species and occupied the greatest diversity of habitats, and both are considered the main malaria vectors in Mexico.8,11 Of the two species mentioned, An. pseudopunctipennis exhibited the greatest larval density and was widely distributed throughout the country, in accordance with the published geographic distribution.29 This species is reported to be the main vector between 200 and 500 masl, and to be associated with river margins with filamentous algae14 and positively associated with altitude in southern Mexico.15 In this national study, this species was found in an altitude range that goes from 3 to 2 279 masl, and no association was found with the altitude in the multiple logistic regression analysis. The reason for this lack of association is that this species did not exhibit at the national level a similar pattern to that reported in Tapachula, Chiapas, where populations are abundant at the mountain foothills, at an altitude above 200 masl and up to 500 masl, and is not commonly found in the coastal plain.14,15 According to the results, we infer that, nationally, the altitude range where this vector may transmit malaria could be greater than the previously reported one of 200 to 500 masl.3,10 These results agree with the findings of a study in the Yucatán Peninsula, which reports An. pseudopunctipennis in a cenote (its natural, permanent habitat) at 20 masl.30On the other hand, there is a negative association with water turbidity, which indicates that females prefer to lay their eggs in clear rather than turbid water. This coincides with a previous report from 10 countries, spanning most of the geographic range of the species, whose larvae were found in clear-water breeding sites in 90% of the cases (i.e. in 54 out of 60 larval habitats), but is not consistent with the presence of An. pseudopunctipennis in turbid water habitats contaminated with cattle feces in the city of Salinas, Ecuador, and other turbid breeding sites caused by flooding in Monterrey, Mexico.31 This study is also consistent with various studies about the presence of larval habitats with filamentous algae and their disappearance during the rainy season caused by the strong flow of water.14,31

An. albimanus was found in the greatest diversity of hydrologic types (11 out of 12), indicating that this species is generalist when it comes to its preference for larval breeding sites.32 It is the species with the broadest geographic distribution, from Chiapas to Sinaloa on the Pacific side, and from Yucatán to Tamaulipas on the Gulf side, which is consistent with previous reports.29,32 Larval habitats for An. albimanus were characterized as sunny and were positively associated with turbidity, but negatively associated with altitude, which is consistent with reports for southern Chiapas.15,17 On the Pacific side, there are breeding sites for this species all along the coast up to Sinaloa, but no positive larval habitats were found further north. This habitat limitation may be due to the high temperatures and arid climate persistent in Sonora, which is consistent with previous studies.32

Another abundant species was An. franciscanus, which was found only in northern Mexico, the Nearctic region. The geographic distribution of An. franciscanus in this study in northern Mexico apparently is not consistent with a previous report by Vargas and Martínez-Palacios (1956),6 who reported this species only in northern Baja California. At the time of the study, in 1956, the status of An. franciscanus as a species was unclear, a problem that was solved when genetic incompatibility was found to exist between An. pseudopunctipennis and An. franciscanus, and it was concluded that these were two different species.33 The literature also reported An. pseudopunctipennis var Willardi in Sonora, Chihuahua, Coahuila and San Luis Potosí, but this variety disappeared and was integrated to the species An. franciscanus.34,35 Finally the results of the present study are consistent with those studies regarding the geographic distribution reported for northern México.6In this study, this species was not found in Tamaulipas, as was previously reported,36 but was found for the first time in Nuevo León and Durango. The lack of reports for this species in Nuevo León could be derived from the lack of taxonomic studies in this geographic area, or possibly because of a range expansion derived from anthropogenic changes or climate change.12,13 Apparently An. franciscanus has little importance as a malaria vector today or in the past;34 however, it is necessary to deepen taxonomical and binomial analyses in the states of Sinaloa, Nayarit, Durango and Chihuahua, where malaria cases have been reported.3 The vector An. vestitipennis was found only in southern Mexico, in the Neotropical region, where it is one of the main malaria vectors in the Lacandon Forest.9 Larval habitats for this species in the Lacandon Forest were found in ponds, ditches and rain puddles. The most abundant larval habitat were ditches, which are caused by man-made changes in orography, because road building avoids natural river paths, allowing the formation of water bodies that are attractive for this species. The same pattern coincides in a report from Colombia, where road building favored the development of several species of Anopheles.32 Further studies are required about the incrimination of vectors in the regions of Soconusco, Chiapas, and Catemaco, Veracruz, where this species was found and its role as a malaria vector is unknown. This study describes for the first time the characteristics of the larval habitats of this species in the Lacandon Forest area in Chiapas.

The greatest diversity of anopheline species was found in the Lacandon Forest area in Chiapas, with 11 species and the highest Fisher’s α diversity value: 1.20. These results are consistent with studies in other regions of the world, where it has been confirmed that tropical zones have the greatest biodiversity of insect species.37 It must be noted that one of the main historical references of abundance and distribution of Mexican anophelines (published in 1956),6 does not mention any anopheline species for the Lacandon Forest in Chiapas. At that time in the 1950s, despite there being populated areas in the Lacandon Forest, there were few taxonomic works,38,39,40 mainly due to the difficult access to these areas.41 At present, despite the Lacandon Forest being the most important malaria focus in Mexico,3 very little has been written about vectors.9 This is the first report that exposes the high anopheline biodiversity for this area. The high incidence of malaria cases in this region could be accounted for by the high diversity of species, as was reported in Colombia, where there was a positive correlation between richness of anopheline species and risk of transmission.42,43 Another phenomenon that may cause the high incidence of malaria could be the movement of people between communities for work or business, or the emigration of people from endemic areas who are traveling north or to the United States, which results in elevated transmission risk for P. vivax and P. falciparum in this area.3,44

The second region with the greatest species diversity, with a Fisher’s α diversity index of 1.0, is the state of Coahuila, in northeastern Mexico, where we found five species: An. crucians, An. franciscanus, An. punctipennis, An. pseudopunctipennis and An. albimanus. In this area, the climate is extreme and semi-desertic (Koppen BS). This area in northeastern Mexico is free of malaria despite the high species diversity. The explanation is complex and may be attributable to environmental factors such as extreme climate with very hot summers that reach up to 45 °C and cold winters, as low as -8 °C, most human dwellings are made of concrete and use screens in the windows, compared to Oaxaca, Chiapas and Tabasco in southeastern Mexico, where most houses in the rural areas are made of plant matter or wood with palm roofs.23 On the other hand, vector characteristics such as abundance, seasonal distribution, life expectancy, degree of antropophilia and vector competition may influence the malaria transmission cycle.45,46

The success of focalized control is based on community involvement or the search and treatment of sick people and environmental management through elimination of anopheline breeding sites (the ECA, “Eliminación de Criaderos de Anofelinos”).3,19 This control strategy had a resounding success in Oaxaca in the nineties and may be used as a model to be implemented in other areas with active transmission or during new disease outbreaks. In order to ensure the elimination of malaria in Mexico, it is necessary to eliminate or manage those larval habitats at risk in or near communities, applying antilarval measures and based on the knowledge about larval ecology in every specific region and for every anopheline species.46

Acknowledgements

We thank Dr. Lilia González Cerón, PhD, for her comments and suggestions in relation to the manuscript. We also thank José Antonio Zavala and Óscar Reyes Gálvez for their technical assistance.

Funding

This work would not have been possible without the financial support of Consejo Nacional de Ciencia y Tecnología (National Council for Science and Technology, Conacyt) Mexico City, grant 105806.Funding 140. This work would not have been possible without the financial support of Consejo Nacional de Ciencia y Tecnología (National Council for Science and Technology, Conacyt) Mexico City, grant 105806.

References

World Health Organization. World Malaria Report 2016. Geneva: WHO, 2016 [cited April 1, 2017]. Available form:Available form:https://www.who.int/malaria/publications/world-malaria-report-2016/report/en/Links ]

Secretaría de Salud. Reducción de los casos de paludismo en México. México: SS, 2016 [cited May 10, 2017]. Available from:Available from:https://www.gob.mx/salud/prensa/reduccion-del-93-por-ciento-los-casos-de-paludismo-en-mexicoLinks ]

Betanzos-Reyes A. La malaria en México. Progresos y desafíos hacia su eliminación. Bol Med Hosp Infant Mex. 2011;68(2):159-68. [ Links ]

Harbach R. Genus Anopheles Meigen, 1818 [Internet]. Mosquito Taxonomic Inventory. 2008 [cited April 1, 2017]. Available from:Available from:http://mosquito-taxonomic-inventory.info/simpletaxonomy/term/6047Links ]

Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, et al. Developing global maps of the dominant Anopheles vectors of Human Malaria. PLoS Med. 2010;7(2):e1000209. https://doi.org/10.1371/journal.pmed.1000209 [ Links ]

Vargas L, Martínez-Palacios A. Anofelinos mexicanos, taxonomía y distribución. México City: Comisión Nacional para la Erradicación del Paludismo/Secretaría de Salubridad y Asistencia, 1956. [ Links ]

Wilkerson RC, Strickman D, Litwak TR. Illustrated key to the female anopheline mosquitoes of Central America and Mexico. J Am Mosq Control Assoc. 1990;6:7-34. [ Links ]

Rodríguez MH, Loyola EG. Situación epidemiológica actual y perspectivas de la investigación entomológica en México. In: Sociedad Mexicana de Entomología. Memorias de IV Simposio Nacional de Entomología Médica y Veterinaria. Oaxtepec, Morelos, Mexico: SME, 1989:15-40. [ Links ]

Loyola EG, Arredondo-Jiménez JI, Rodríguez MH, Bown DN, Vaca-Marín MA.Anopheles vestitipenis, the probable vector of Plasmodium vivax in the Lacandon forest of Chiapas, Mexico. Trans R Soc Trop Med Hyg. 1991;85:171-4. https://doi.org/10.1016/0035-9203(91)90010-V [ Links ]

Fernández-Salas I, Rodríguez MH, Roberts DR, Rodríguez MC, Wirtz RA. Bionomics of adult Anopheles pseudopunctipennis (Diptera: Culicidae) in the Tapachula foothill area of southern Mexico. J Med Entomol. 1994;31:663-70. https://doi.org/10.1093/jmedent/31.5.663 [ Links ]

Villarreal-Treviño C, Arredondo-Jiménez JI, Rodríguez MH. Bionomía de los principales vectores del paludismo en México. In: Kumate J, Martínez-Palomo A (eds.). A cien años del descubrimiento de Ross. El Paludismo en México. Mexico City: El Colegio Nacional, 1998;149-65. [ Links ]

Patz JA, Graczyk TK, Geller N, Vittor AY. Effects of environmental change on emerging parasitic diseases. Int J Parasitol. 2000;30(12-13):1395-405. https://doi.org/10.1016/S0020-7519(00)00141-7 [ Links ]

Vora N. Impact of anthropogenic environmental alterations on vector-borne diseases. Medscape J Med. 2008;10(10):238. [ Links ]

Fernández-Salas I, Roberts DR, Rodríguez MH, Marina-Fernández CF. Bionomics of larval populations ofAnopheles pseudopunctipennisin the Tapachula foothills area, southern Mexico. J Am Mosq Control Assoc. 1994;4:477-86. [ Links ]

Savage HM, Rejmankova E, Arredondo-Jiménez JI, Roberts DR, Rodríguez MH. Limnological and botanical characterization of larval habitats for two primary malarial vectors, Anopheles albimanus and Anopheles pseudopunctipennis, in coastal areas of Chiapas State, Mexico. J Am Mosq Control Assoc. 1990;6:612-20. [ Links ]

Faran ME. Mosquito studies (Diptera, Culicidae) XXXIV. A revision of the Albimanus section of the subgenus Nyssorhynchus of Anopheles. Contrib Am Entomol Inst. 1980;15(7):1-215. [ Links ]

Rodríguez AD, Rodríguez MH, Meza RA, Hernández JE, Rejmankova E, Savage HM, et al. Dynamics of population densities and vegetation associations of Anopheles albimanus larvae in a coastal area of southern Chiapas, Mexico. J Am Mosq Control Assoc. 1993;9(1):46-58. [ Links ]

Villarreal-Treviño C, Penilla-Navarro RP, Vázquez-Martínez MG, Moo-Llanes DA, Ríos-Delgado JC, Fernández-Salas I, Rodríguez AD. Larval habitat characterization of Anopheles darlingi from its northernmost geographical distribution in Chiapas, México. Malar J. 2015;14:517. https://doi.org/10.1186/s12936-015-1037-0 [ Links ]

Chanon KE, Méndez-Galván JF, Galindo-Jaramillo JM, Olguín-Bernal H, Borja-Aburto VH. Cooperative actions to achieve malaria control without the use of DDT. Int J Hyg Environ Health. 2003;206:387-94. https://doi.org/10.1078/1438-4639-00235 [ Links ]

Bond JG, Rojas JC, Arredondo-Jiménez JI, Quiroz-Martínez H, Valle J, Williams T. Population control of the malaria vectorAnopheles pseudopunctipennisby habitat manipulation. Proc Biol Sci. 2004;271:2161-9. https://doi.org/10.1098/rspb.2004.2826 [ Links ]

Morrone JJ, Márquez J. Biodiversity of Mexican terrestrial arthropods (Arachnida and Hexapoda): A biogeographic puzzle. Acta Zoológica Mexicana. 2008;24:15-41. https://doi.org/10.21829/azm.2008.241613 [ Links ]

Roiz D, Ruiz S, Soriguer R, Figuerola J. Climatic effects on mosquito abundance in Mediterranean wetlands. Parasit Vectors. 2014;7:333. https://doi.org/10.1186/1756-3305-7-333 [ Links ]

Instituto Nacional de Estadística y Geografía. Climatología. Mapa digital. México: INEGI [cited January 5, 2017]. Available from: Available from: https://www.inegi.org.mx/temas/climatologia/Links ]

Weatherspark. El clima típico de cualquier lugar del mundo. Datos ambientales [cited May 16, 2017]. Available from:Available from:https://es.weatherspark.com/Links ]

Villarreal C, Arredondo-Jiménez JI, Rodríguez MH, Ulloa A. Colonization of Anopheles pseudopunctipennis from Mexico. J Am Mosq Control Assoc. 1998;14:369-72. [ Links ]

Becker N, Zgomba M, Petric D, Dahl C, Boase C, Lane J, et al. Mosquitoes and their control. New York: Springer Science & Business Media, 2013. [ Links ]

Magurran AE. Measuring ecological diversity. Oxford: Blackwell Publishing, 2004. [ Links ]

Seaby RM, Henderson PA. Species diversity and Richness IV. Lymington, Hampshire, England: Pisces Conservation, 2007. [ Links ]

Malaria Atlas Project. The Malaria Atlas Project website [cited July 1, 2017]. Available from:Available from:http://www.map.ox.ac.ukLinks ]

Baak-Baak CM, Cigarroa-Toledo N, Arana-Guardia R, Chi-Chim WA, Chan Orilla JA, Machain-Williams C, et al. Mosquito fauna associated with Aedes aegypti (Diptera: Culicidae) in Yucatán State of southeastern México, and checklist with new records. Fla Entomol. 2016;99(4):703-9. https://doi.org/10.1653/024.099.0420 [ Links ]

Manguin S, Roberts DR, Peyton EL, Rejmankova E, Pecor J. Characterization ofAnopheles pseudopunctipennislarval habitats. J Am Mosq Control Assoc . 1996 12(4):619-26. [ Links ]

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. 2012;30:11:13. https://doi.org/10.1186/1476-072X-11-13 [ Links ]

Smithson TW. Species rank for Anopheles franciscanus based on failure of hybridization with Anopheles pseudopunctipennis pseudopunctipennis. J Med Entomol. 1972;9(6):501-5. https://doi.org/10.1093/jmedent/9.6.501 [ Links ]

Russell PE, Rozeboom LE, Stone A. Keys the anopheline mosquitoes of the world with notes on their identification, distribution, biology and relation to malaria. Philadelphia, PA: The American Entomological Society/The Academy of Natural Science,1943. [ Links ]

Darsie RF Jr. A survey and bibliography of the mosquito fauna of Mexico (Diptera: Culicidae). J Am Mosq Control Assoc . 1996;2:298-306. [ Links ]

Ortega-Morales AI, Zavortink TJ, Huerta-Jiménez H, Sánchez-Rámos FJ, Valdés-Perezgasga MT, Reyes-Villanueva F, et al. Mosquito records from México: The Mosquitoes (Diptera: Culicidae) of Tamaulipas State. J Med Entomol. 2015;52(2):171-84. https://doi.org/10.1093/jme/tju008 [ Links ]

Gadagkar R, Chendrashekara K, Nair P. Insect species diversity in the tropics: Sampling methods and a case study. J Bombay Nat.1990;87(3):337-53. [ Links ]

Instituto Nacional de Ecología. Programa de manejo Reserva de la biosfera Montes Azules, México. Mexico City: Instituto Nacional de Ecología Secretaría del Medio Ambiente, Recursos Naturales y Pesca (National Institute of Ecology/Ministry of the Environment, Natural Resources and Fisheries), 2000. [ Links ]

Kumm HW, Bustamante ME, Herrera JR. Report concerning certain anophelines found near the Mexican-Guatemalan frontier. Am J Trop Med. 1943;32:373-6. https://doi.org/10.4269/ajtmh.1943.s1-23.373 [ Links ]

Vargas L. El Anopheles darlingi Root, 1926 en México. Revista del Instituto de Salubridad y Enfermedades Tropicales. 1946;7(4): 221-6. [ Links ]

Knight KL, Stone A. A catalog of the Mosquitoes of the World (Diptera: Culicidae). 2nd ed. College Park, Maryland: Entomological Society of America, 1977. [ Links ]

Fuller DO, Alimi T, Herrera S, Beier JC, Quiñones ML. Spatial association between malaria vector species richness and malaria in Colombia. Acta Trop. 2016;158:197-200. https://doi.org/10.1016/j.actatropica.2016.03.008 [ Links ]

Herrera-Varela M, Orjuela LI, Peñalver C, Conn JE, Quiñones ML. Anopheles species composition explains differences in Plasmodium transmission in La Guajira, northern Colombia. Mem Inst Oswaldo Cruz. 2014;109(7):952-6. https://doi.org/10.1590/0074-0276140126 [ Links ]

Secretaría de Salud. Programa de Acción Específico. Prevención y Control del Paludismo 2013-2018. Mexico City: Secretaría de Salud, 2014 [cited May 1, 2018]. Available from: Available from: http://www.cenaprece.salud.gob.mx/descargas/pdf/PAE_PrevencionControlPaludismo2013_2018.pdfLinks ]

Day JF. Host-seeking strategies of mosquito disease vectors. J Am Mosq Control Assoc . 2005;21(4 Suppl):17-22. https://doi.org/10.2987/8756-971X(2005)21[17:HSOMDV]2.0.CO;2 [ Links ]

Beerntsen BT, James AA, Christensen BM. Genetics of mosquito vector competence. Microbiol Mol Biol Rev. 2000;64(1):115-37. https://doi.org/10.1128/MMBR.64.1.115-137.2000 [ Links ]

Received: October 08, 2018; Accepted: June 27, 2019

Corresponding author: Cuauhtémoc Villarreal Treviño. Centro Regional de Investigación en Salud Pública, Instituto Nacional de Salud Pública. 4 av. Norte y 19 calle poniente. Tapachula, Chiapas, México. Email: cvilla@insp.mx

Declaration of conflict of interests. The authors declare that they have no conflict of interests.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License