RUBÉN BUENO MAR͹ y RICARDO JIMÉNEZ PEYDRÓ¹

¹ Ph.D. Laboratorio de Entomología y Control de Plagas Instituto Cavanilles de Biodiversidad y Biología Evolutiva Universitat de València (Estudi General). Apartado oficial 22085.46071 Valencia (España). ruben.bueno@uv.es autor para correspondencia.

Recibido: 28-jun-2010 - Aceptado: 30-abr-2011

Profundizar en el conocimiento de los factores ambientales que afectan a la ocurrencia de las diferentes especies de mosquitos (Diptera: Culicidae) es sumamente importante para el desarrollo de exitosos programas de control poblacional (Axtell 1979; Bueno Marí y Jiménez Peydró 2011). En este sentido, existen diversos estudios que se basan en el empleo de regresiones logísticas para predecir la incidencia de los culícidos a partir de diversas variables ambientales (Udevitz et al. 1987), sobre todo aquéllas relacionadas con las características físicas y químicas de los biotopos larvarios (Vniska y Pappas 1984; Blaustein y Margalit 1995; Sunish y Reuben 2001). Con la finalidad de profundizar en el conocimiento de estas cuestiones y aportar información interesante acerca de las variables que mejoren la predictibilidad de la ocurrencia de los culícidos mediterráneos, se investigaron las poblaciones larvarias presentes en el Levante español de estos dípteros nematóceros. El objetivo fue establecer qué variables cuantitativas y cualitativas de los biotopos larvarios se asocian mejor con la presencia de las distintas especies.

Área de estudio y metodología de análisis de campo. Entre los meses de marzo a octubre de los años comprendidos entre 2005 y 2008 se llevó a cabo un muestreo larvario de culícidos en la Comunidad Valenciana (Fig. 1); región de clima típicamente Mediterráneo situada en el este de la Península Ibérica, que ocupa una extensión aproximada de 23.260 km2 y que presenta una importante heterogeneidad de ambientes hídricos. La captura de larvas se realizó mediante la técnica 'dipping" (Service 1993) y para la determinación específica de los ejemplares recolectados se emplearon los criterios taxonómicos de Encinas Grandes (1982), Darsie y Saminadou Voyadjoglou (1997) y Schaffner et al. (2001). Paralelamente a la colecta de juveniles, se procedió a la determinación del pH, temperatura, conductividad y salinidad in situ empleando un equipo portátil multiparamétrico (pIONneer65) con una sonda combinada para pH y temperatura (pHC5977), y otra para la conductividad y temperatura (CDC30T). La determinación de la salinidad se calculó automáticamente con el instrumental basándose en el 'Extended Practical Salinity Scale" de 1978 (Lewis 1980). Para la realización de estos análisis se sumergieron las respectivas sondas en el cuerpo hídrico muestreado hasta la estabilización de la medida. Del mismo modo, los puntos de muestreo fueron georreferenciados a partir de coordenadas GMS (Grados: Minutos: Segundos) establecidas mediante un geoposicionador (modelo GPS 12 de Garmin, datum WGS 84), a través del cual también se determinó la altitud. Además, se tipificó cada punto de muestreo recogiendo cuáles de las siguientes variables cualitativas se ajustaba mejor a su condición:

− Tipo de ambiente acuático

* Léntico

* Lótico

* Río/arroyo/riachuelo/escorrentía

* Laguna temporal o permanente

* Embalse

* Encharcamiento temporal/arrozal/zona de marismas

* Canal o acequia

* Gran receptáculo: albergan cantidades de agua superiores a 10 litros (balsa, bidón, fuente, abrevadero, pozo, piscina,…)

* Pequeño recipiente: albergan cantidades de agua inferiores a 10 litros (cubo, macetero, neumático,…)

− Tipo de sustrato (se considera el mayoritario en el biotopo en caso de coexistencia)

* Rocoso

* Arenoso

* Cenagoso

* Hormigón

* Plástico/caucho/metal

− Profundidad del foco (medición en el punto más profundo muestreado)

* Aguas someras (menos de 40 cm)

* Aguas profundas (más de 40 cm)

− Grado de antropización del ambiente

* Bajo (medios naturales o salvajes, sin actividad humana aparente en áreas circundantes)

* Alto (medios urbanos o periurbanos)

Análisis estadístico. Mediante el empleo del programa estadístico SPSS 18 (SPSS 2009) se realizaron diversas regresiones logísticas binarias por pasos hacia delante, según la razón de verosimilitud, con el objetivo de analizar la predictibilidad de la incidencia (presencia/ausencia) de las diferentes especies de culícidos capturadas, a partir de los datos cuantitativos y cualitativos registrados de cada biotopo larvario. Cada regresión logística construyó una ecuación de regresión para predecir la incidencia de cada especie, a partir de una combinación lineal de las variables predictoras cuantitativas y cualitativas que se fueron recopilando durante el presente estudio.

Dados los beneficios reportados a partir de la dicotomización (Aguayo Canela y Lora Monge 2007) de las variables continuas (físicas y químicas del agua y altitud del biotopo) y de las categóricas policotómicas (variables cualitativas de foco acuático), se decidió transformar a sistema binario (0/1) todas las variables observadas (Tabla 1).

Puesto que la decisión acerca de la inclusión definitiva de una variable en la ecuación depende únicamente de los criterios establecidos del investigador, debido a la ausencia de pruebas estadísticas al respecto, se emplearon las pautas propuestas por Aguayo y Lora (2007) para decidir qué variables quedaban incluidas en el modelo de predicción. Por tanto fueron incorporadas aquellas variables que al ser introducidas:

– Modificaron el valor de la 'estimación de la medida de asociación" (Exp(B)) de la variable principal de forma considerable (la mayoría de autores afirman que este cambio debe ser de, al menos, un 10% del valor original).

– No aumentaron sustancialmente los intervalos de confianza de la Exp(B) de la variable principal, puesto que de lo contrario habría aumentado la imprecisión de la estimación.

– No cambiaron la significación estadística del contraste de Wald para el coeficiente de regresión de la variable principal.

– Presentaron un intervalo de confianza al 95% que no contenía el valor nulo (= 1).

Se identificaron un total de 679 criaderos larvarios y se capturaron 11.279 ejemplares de culícidos pertenecientes a 29 especies así: Anopheles algeriensis Theobald, 1903, Anopheles atroparvus Van Thiel, 1927, Anopheles claviger (Meigen, 1804), Anopheles maculipennis Meigen, 1818, Anopheles marteri Senevet & Prunnelle, 1927, Anopheles petragnani Del Vecchio, 1939, Anopheles plumbeus Stephens, 1828, Aedes vexans (Meigen, 1830), Aedes vittatus (Bigot, 1861), Ochlerotatus echinus (Edwards, 1920), Ochlerotatus geniculatus (Olivier, 1791), Ochlerotatus gilcolladoi (Sánchez- Covisa, Rodriguez & Guillén, 1985), Ochlerotatus berlandi (Séguy, 1921), Ochlerotatus caspius (Pallas, 1771), Ochlerotatus detritus (Haliday, 1833), Ochlerotatus pulcritarsis (Rondani, 1872), Culex modestus Ficalbi, 1889, Culex laticinctus Edwards, 1913, Culex mimeticus Noe, 1899, Culex pipiens Linnaeus, 1758, Culex theileri Theobald, 1903, Culex hortensis Ficalbi, 1889, Culex impudicus Ficalbi, 1890, Culex territans Walker, 1856, Culiseta longiareolata (Macquart, 1838), Culiseta annulata (Schrank, 1776), Culiseta subochrea (Edwards, 1921), Orthopodomyia pulcripalpis (Rondani, 1872) y Uranotaenia unguiculata Edwards, 1913. Sin embargo, únicamente se obtuvieron óptimos porcentajes acerca de las siguientes cuatro especies:

- Oc. caspius: especie paleártica muy abundante en terrenos pantanosos costeros (Edwards 1921; Rioux 1958) gracias a su destacada resistencia frente a la salinidad del agua, pudiendo incluso desarrollarse en ambientes hipersalinos, influenciados por la evaporación, con concentraciones de hasta 106 g/l (López Sánchez 1989). Es multivoltina, exófila, exofágica, eurígama y sus huevos resisten la desecación y las temperaturas extremas, siendo, consecuentemente habituales, los estados hibernante y estivante. La abundancia y distribución de estos huevos en medios concretos parece guardar cierta relación con el tipo de vegetación predominante, la cual se utilizaría como un elemento indicativo y fiable de las condiciones ambientales (López Sánchez 1989). Las hembras se alimentan preferentemente en horas diurnas, picando con gran agresividad y frecuencia al ser humano, y pueden desplazarse hasta 40 kilómetros en busca de hospedadores adecuados (Schaffner et al. 2001). Su interés vectorial es extenso, englobando a la transmisión de filarias (Aranda et al. 1998) y arbovirus como West Nile o Tahyna (Schaffner et al. 2001).

La inclusión de las variables 'TDS alto" y 'Encharcamientos" en la ecuación de regresión permitió la correcta predicción de la presencia de Oc. capius en un 67,2 % de los casos (Tablas 2 y 3). Ambas variables están íntimamente relacionadas con las características bioecológicas descritas para la especie ya que Oc. caspius tiende a anidar en ambientes de intermitencia hídrica de elevada salinidad (Bueno Marí 2011).

- Cx. pipiens: Especie distribuida por toda la región Holártica, el este y sudeste de África y Sudamérica (Stone et al. 1959). Es multivoltina, pudiendo completar hasta seis ciclos gonotróficos y, dada su ubicuidad, también es conocida como el 'mosquito común". Su adaptabilidad a los ambientes urbanos, unida a su elevada endofilia, le permite ser una especie muy frecuente en el interior de las viviendas, picando ávidamente al ser humano durante la noche. Es capaz de transmitir numerosos arbovirus, tales como West Nile, Sindbis, Tahyna, Batai y Usutu, además de filarias como Dirofilaria immitis (Leidy, 1856) y plasmodios de afección aviar (Aranda et al. 1998; Schaffner et al. 2001; Bueno Marí y Jiménez Peydró 2010).

Las regresiones logísticas revelan que los 'canales y acequias" 'oligosalobres" caracterizados por una 'antropización alta", permiten el correcto pronóstico de la presencia de Cx. pipiens en el 60% de los casos estudiados (Tablas 4 y 5).

- Cx. territans: Especie propia de la región Holártica ampliamente representada en el continente europeo (Horsfall 1972). Es multivoltina, eurígama y típicamente asociada a ambientes montañosos (Bueno Marí et al. 2009). Las hembras hibernan en cavidades tanto naturales como antrópicas, pican a reptiles y anfibios, y transmiten filarias de afección en batracios (Encinas Grandes 1982; Schäefer 2004).

Los ambientes fluviales 'ríos o arroyos" de baja salinidad 'dulces o hiposalinos" posibilitaron pronosticar con éxito la presencia de Cx. territans en el 60% de los casos (Tablas 6 y 7). Al respecto cabe señalar que Petersen y Chapman (1969) y Vniska y Pappas (1984), también evidenciaron una elevada asociación de Cx. territans con cuerpos hídricos de baja salinidad.

- Cs. longiareolata: Especie de distribución paleártica, oriental y afro-tropical (Schaffner et al. 2001). En Europa su presencia es más común en los países mediterráneos, siendo excepcionales sus hallazgos en regiones más norteñas del continente. Es multivoltina, estenógama, autógena, sus hábitats larvarios son muy variados y, mientras en regiones templadas hiberna en estado larvario; en áreas más frías son las hembras quienes se encargan de la supervivencia de la especie durante el periodo desfavorable (Encinas Grandes 1982; Schaffner et al. 2001). Se han descrito casos de canibalismo en larvas L4 de Cs. longiareolata, alimentándose de estados larvarios menos desarrollados de la misma u otras especies o incluso también de adultos en los momentos previos a su emergencia (Al-Saadi y Mohsen 1988). Sus preferencias hematofágicas parecen dirigirse hacia las aves, aunque ocasionalmente piquen también al hombre, y pueden transmitir el paludismo aviar y el Virus West Nile (Schaffner et al. 2001). Las variables cualitativas 'Gran receptáculo" y 'Sustrato plástico" permitieron predecir correctamente el 67% de las presencias de Cs. longiareolata (Tablas 8 y 9). Por tanto, se ratifica la predilección de la especie por desarrollar su etapa preimaginal en recipientes que habitualmente se encuentran en entornos domésticos o peridomésticos (Bueno Marí et al. 2010). La no inclusión de variables químicas del agua en el modelo predictivo, está en concordancia con investigaciones realizadas por otros investigadores que sugieren que la especie es euriónica y eurihalina (Dimentman y Margalit 1981; Margalit et al. 1988).

Existen diferentes variables cuantificables en los biotopos larvarios de los culícidos que pueden ser utilizadas para predecir la incidencia de estos dípteros nematóceros cuyo interés sanitario es de primer orden mundial. En este sentido es necesario, no sólo profundizar en estos estudios, sino también regionalizarlos, ya que es de sobra conocido que las mismas especies pueden exhibir comportamientos claramente dispares en función de la hidrografía y clima local. En conclusión, los sistemas de predicción de mosquitos culícidos nos permiten disponer de una herramienta útil desde el punto de vista del control, ya que posibilita la identificación de criaderos de especies concretas con un elevado grado de fiabilidad.

Los autores agradecemos a Mario Sendra, Profesor Titular del Departamento de Estadística e Investigación Operativa de la Universitat de València, por su asesoramiento en el tratamiento estadístico de los datos. Asimismo, queremos indicar nuestra gratitud a la Conselleria de Medi Ambient, Aigua, Urbanisme i Habitatge de la Generalitat Valenciana, por la concesión de los permisos pertinentes para la captura de insectos. Por último, señalar que este trabajo ha sido financiado, en parte, por el proyecto CGL 2009-11364 (BOS) otorgado por el Ministerio de Ciencia e Innovación del Gobierno de España.

AGUAYO CANELA, M.; LORA MONGE, E. 2007. Cómo hacer una Regresión Logística binaria "paso a paso" (II): análisis multivariante. DocuWeb-fabis. Huelva: Fundación Andaluza Beturia para la Investigación en Salud. Disponible en: http://www.fabis.org/html/?page_name=DocuWeb [Fecharevisión: 15 junio 2010]        [ Links ]

AL-SAADI, M.; MOHSEN, Z. H. 1998. Predatory and cannibalistic behavior of Culiseta longiareolata (Macquart) (Diptera: Culicidae) in Iraq. Journal of Biological Sciences Research 19 (2): 339-351.        [ Links ]

AXTELL, R. C. 1979. Principles of integrated pest management (IPM) in relation to mosquito control. Mosquito News 39 (4): 709-718.        [ Links ]

ARANDA, C.; PANYELLA, O.; ERITJA, R; CASTELLA, J. 1998. Canine filariasis Importance and transmission in the Baix Llobregat area, Barcelona (Spain). Veterinary Parasitology 77 (4): 267-275.        [ Links ]

BLAUSTEIN, L.; MARGALIT, J. 1995. Spatial distributions of Culiseta longiareolata (Culicidae: Diptera) and Bufo viridis (Amphibia: Bufonidae) among and within desert pools. Journal of Arid Environments 29 (2): 199-211.        [ Links ]

BUENO MARÍ, R. 2011. Bioecología, diversidad e interés epidemiológico de los culícidos mediterráneos (Diptera, Culicidae) [tesis doctoral]. Servei de Publicacions de la Universitat de València. Valencia (España). 425 pp.        [ Links ]

BUENO MARÍ, R.; CHORDÁ OLMOS, F. A.; BERNUÉS BAÑERES, A.; JIMÉNEZ PEYDRÓ, R. 2009. Aportaciones al conocimiento de los mosquitos (Diptera, Culicidae) de alta montaña presentes en la Península Ibérica. Pirineos 164: 49-68.        [ Links ]

BUENO MARÍ, R.; JIMÉNEZ PEYDRÓ, R. 2010. Situación actual y eco-epidemiología de las arbovirosis transmitidas por mosquitos culícidos (Diptera: Culicidae) en España. Revista Española de Salud Pública 84 (3): 255-269.        [ Links ]

BUENO MARÍ, R.; CORELLA LÓPEZ, E.; JIMÉNEZ PEYDRÓ, R. 2010. Culícidofauna (Diptera: Culicidae) presente en los distintos enclaves hídricos de la ciudad de Valencia (España). Revista Colombiana de Entomología 36 (2): 235-241.        [ Links ]

BUENO MARÍ, R.; JIMÉNEZ PEYDRÓ, R. 2011. Classification of Spanish mosquitoes in functional groups. Journal of the American Mosquito Control Association 27 (1): 1-7.        [ Links ]

EDWARDS, E. W. 1921. A revision of the mosquitoes of the Paleartic Region. Bulletin of Entomological Research 12 (3): 263-351.        [ Links ]

ENCINAS GRANDES, A. 1982. Taxonomía y biología de los mosquitos del área salmantina (Diptera, Culicidae) [tesis doctoral]. CSIC, Centro de edafología y Biología aplicada. Ed. Universidad de Salamanca. Salamanca (España). 437 pp.        [ Links ]

DARSIE, R. F.; SAMINADOU VOYADJOGLOU, A. 1997. Keys for the identification of the mosquitoes of Greece. Journal of the American Mosquito Control Association 13 (3): 247-254.        [ Links ]

DIMENTMAN, C.; MARGALIT, J. 1981. Rainpools as breeding and dispersal sites of mosquitoes and other aquatic insects in the Central Negev. Journal of Arid Environments 4: 123-129.        [ Links ]

HORSFALL, W. R. 1972. Mosquitoes: their bionomics and relation to disease. Ed. Hafner Publications. New York (United States of America). 723 pp.        [ Links ]

LEWIS, E.L. 1980. The practical salinity scale 1978 and its antecedents. IEEEJ. Oceanic Engineering OE-5:3.        [ Links ]

LÓPEZ SÁNCHEZ, S. 1989. Control integral de mosquitos en Huelva. Ed. Junta de Andalucía, Consejería de Salud y Servicios Sociales. Sevilla (España). 340 pp.         [ Links ]

MARGALIT, J., DIMENTMAN, C., TAHORI, A.S. 1988. Geographical, seasonal and ecological distribution of mosquito larvae (Diptera: Culicidae) in Southern Israel. Archives für Hydrobiologie 112 (2): 233-249.        [ Links ]

RIOUX, J. A. 1958. Les culicides du midi Mediterreneen. Etude systematique et ecologique. Le Chevalier. Paris (France). 304 pp.        [ Links ]

PETERSEN, J. J.; CHAPMAN, H. C. 1969. Chemical factors of water in tree holes and related breeding of mosquitoes. Mosquito News 29: 29-36.        [ Links ]

SCHÄEFER, M. 2004. Mosquitoes as a Part of Wetland Biodivesity [doctoral thesis]. Faculty of Science and Technology, Ed. University of Uppsala. Uppsala (Sweden). 63 pp.        [ Links ]

SCHAFFNER, F.; ANGEL, G.; GEOFFROY, B.; HERVY, J. O.; RHAEIM, A. 2001. The mosquitoes of Europe / Les moustiques d' Europe [programa de ordenador]. Montpellier, France: IRD Éditions and EID Méditerranée.        [ Links ]

SERVICE, M. W. 1993. Mosquito Ecology. Field Sampling Methods. Ed. Elsevier Science Publishers. Liverpool (United Kingdom). 988 pp.        [ Links ]

SPSS (STATISTICAL PACKAGE FOR THE SOCIAL SCIENCES). 2009. Base 18.0 User's Guide. Ed. SPSS Inc., Chicago (United States of America).        [ Links ]

STONE, A.; KNIGHT, K.; STARCKE, H. 1959. A synoptic catalog of the Mosquitoes of the world (Diptera, Culicidae). Ed. Entomological Society of America, Thomas Say Foundation, 6. Washington (United States of America). 358 pp.        [ Links ]

SUNISH I. P.; REUBEN, R. 2001. Factors influencing the abundance of Japanese encephalitis vectors in rice fields in India - I. Abiotic. Medical and Veterinary Entomology 15 (4): 381-392.        [ Links ]

UDEVITZ M. S.; BLOOMFIELD P.; APPERSON, C. S. 1987. Prediction of the occurrence of four species of mosquito larvae with logistic regression on water-chemistry variables. 1987. Environmental Entomology 16 (1): 281-285.        [ Links ]

VNISKA, L. A.; PAPPAS, L. G. 1984. Chemical analysis of mosquito larval habitats southeastern Nebraska. Mosquito News 44: 506-509.        [ Links ] ]]> 2007 1998 19 2 2 339-351 1979 39 4 4 709-718 1998 77 4 4 267-275 1995 29 2 2 199-211 2011 425 2009 164 49-68 2010 84 3 3 255-269 2010 36 2 2 235-241 2011 27 1 1 1-7 1921 3 3 263-351 1982 437 1997 13 3 3 247-254 1981 4 4 123-129 1972 723 1980 5 1989 340 1988 . 112 2 2 233-249 1958 304 1969 29 29 29-36 2004 63 2001 1993 988 STATISTICAL PACKAGE FOR THE SOCIAL SCIENCES 2009 1959 358 2001 15 4 4 381-392 1987 16 1 1 281-285 1984 44 44 506-509