Introducción

Los cambios en el uso de suelo derivados de las actividades humanas son una de las principales causas de que más del 58% de la superficie terrestre presente un empobrecimiento importante en su biodiversidad, medida como la proporción promedio de biodiversidad remanente en ecosistemas locales (^{Newbold et al., 2016}). Esta pérdida de la biodiversidad está asociada directa o indirectamente con una disminución en la provisión de servicios ambientales que repercuten en el bienestar humano (Millennium ^{Ecosystem Assessment, 2005}). Por ello, es necesario proteger la biodiversidad subsistente en el contexto global de una creciente población humana, que requiere continuar incrementado su producción de alimentos de forma cada vez más eficiente (^{Crist et al., 2017}). Ante este ineludible reto, una de las líneas a seguir es evaluar la capacidad de los paisajes agrícolas para conservar y promover la biodiversidad (^{Johnson et al., 2017}). Los agroecosistemas pueden tener una composición de especies cada vez más diferente, o más similar, a la de los ecosistemas nativos, lo que genera heterogeneidad u homogeneidad biótica en la región (^{Solar et al., 2015}). La homogeneidad biótica entre los agroecosistemas y los ecosistemas nativos indica que el manejo humano permite la conservación de la diversidad regional, por el contrario, la heterogeneidad biótica entre estos ambientes ocurre por la pérdida de especies, o por la presencia de especies invasoras. Por ello, el conocimiento sobre los cambios en la riqueza, diversidad y composición de especies entre áreas de uso agrícola y ecosistemas nativos es esencial para la protección de la biodiversidad regional y el manejo ambiental sostenible de los sistemas agrícolas (^{Bush et al., 2016}; ^{Socolar et al., 2016}).

En el territorio mexicano, cerca del 34% de la superficie terrestre está dedicada a actividades agrícolas y ganaderas (^{Velázquez et al., 2005}), pero poca atención se ha dado al potencial de los sistemas agrícolas como sitios de conservación de la biodiversidad. Por ejemplo, se ha estudiado el papel de algunos agroecosistemas tropicales como sistemas de milpa (e.g., ^{Mateos-Maces et al., 2016}), sistemas agroforestales (e.g., ^{Campos-Salas et al., 2016}) y plantaciones de café bajo sombra (e.g., ^{González-Zamora et al., 2016}) como reservorios de la biodiversidad regional. En las zonas áridas y semiáridas de matorral xerófilo en el Altiplano Central de México uno de los principales cultivos es el de nopal (Opuntia spp., Cactaceae). Las plantaciones comerciales de nopal son agroecosistemas perennes que utilizan distintas variedades nativas domesticadas del género Opuntia, principalmente para la producción de fruta (tuna), nopal verdura y como forraje para el ganado (^{Financiera Rural, 2011}). Sin embargo, a pesar de su importancia económica, son aún pocos los trabajos acerca del potencial de los cultivos de nopal para la conservación de la biodiversidad (^{Riojas-López y Mellink, 2005}). En el estado de Jalisco se ha visto que los agroecosistemas de nopal pueden formar estructuras complejas de hábitats por los estratos de la vegetación, que contribuyen a mantener la biodiversidad y los servicios ecosistémicos (^{Riojas-López, 2006}, ²⁰¹²; ^{Riojas-López y Mellink, 2005}).

Las hormigas (Hymenoptera: Formicidae) se han considerado como un buen grupo indicador para evaluar los cambios ambientales derivados de las actividades humanas, tales como cambios en el uso del suelo por agricultura (e.g., ^{Andersen et al., 2002}). Son un grupo diverso, tanto taxonómica como funcionalmente, y son abundantes en casi todos los ecosistemas y hábitats. En México la fauna de hormigas incluye 973 especies, y las subfamilias más diversas son Myrmicinae y Formicinae, que agrupan el 51 y 20% de la diversidad de especies, respectivamente (^{Ríos-Casanova, 2014}). Entre las ventajas del grupo como indicador del impacto de prácticas agrícolas están su elevada biomasa, su sensibilidad a distintas alteraciones ambientales y la variedad de funciones que realizan en los ecosistemas (^{Andersen et al., 2002}). Estas características hacen que las comunidades de hormigas sean sensibles a los cambios ambientales que reducen la complejidad del hábitat (^{Perfecto y Snelling, 1995}), pues a mayor complejidad estructural en la vegetación, hay una mayor variación microambiental que favorece la riqueza de especies de hormigas (^{Hernández-Ruiz y Castaño-Meneses, 2006}).

En la zona árida del valle del Mezquital, en el estado de Hidalgo, los sistemas agrícolas con bajo nivel de disturbio y buena calidad de agua de riego albergan mayor riqueza y abundancia de hormigas que los sistemas agrícolas irrigados con aguas residuales, y se han identificado los géneros Tapinoma y Odontomachus como bioindicadores de poco disturbio en estos agroecosistemas (^{Hernández-Ruiz y Castaño-Meneses, 2006}). Sin embargo, cuando las prácticas agrícolas se vuelven intensivas, pueden depauperar la diversidad y composición de las comunidades de hormigas (^{Asfiya et al., 2015}; ^{de Castro-Solar et al., 2016}; ^{Hernández-Ruiz y Castaño-Meneses, 2006}; ^{Perfecto y Vandermeer, 2002}). Por ello, se deben mantener remanentes de vegetación natural, aunque sean pequeños, pues pueden tener un potencial importante para el mantenimiento de la diversidad de hormigas en paisajes agrícolas (^{de Castro-Solar et al., 2016}; ^{Gray et al., 2015}; ^{Pacheco et al., 2013}).

Este trabajo tiene como objetivos: 1) comparar la riqueza, abundancia, diversidad y composición de especies epigeas de hormigas entre zonas de matorral xerófilo dominado por Opuntia spp. (de aquí en adelante, nopaleras naturales) y cultivos de nopal tunero (Opuntia albicarpa Scheinvar); 2) detectar posibles especies indicadoras en cada tipo de vegetación y 3) evaluar la influencia de algunas variables del hábitat (volumen de las plantas de nopal, complejidad vertical de la vegetación, número de flores y frutos de nopal), en la riqueza, abundancia y diversidad de las comunidades de hormigas. En el género Opuntia se ha registrado una asociación mutualista con las hormigas que visitan los nectarios extraflorales y brindan protección antiherbívoros y depredadores de semillas, lo que repercute en un incremento del 50% en el éxito reproductivo de la planta, representado por el número de frutos producidos (^{Oliveira et al., 1999}). Sin embargo, no todas las especies de Opuntia tienen nectarios extraflorales y el mutualismo con hormigas es facultativo (^{Miller, 2007}). El estado de Hidalgo es uno de los centros de producción de tuna más importantes en el país (^{Jolalpa-Barrera et al., 2011}), y es el tercer estado con el mayor número de géneros y especies de hormigas (57 y 229 respectivamente), solo por detrás de Veracruz y Chiapas (^{Ríos-Casanova, 2014}), aunque estos números pueden ser menores debido a las sinonimias que se han presentado en los últimos años (^{Vásquez-Bolaños, 2015}). Por ello, se espera encontrar una alta riqueza y diversidad de hormigas tanto en las nopaleras naturales como en los cultivos, y considerando que los cultivos de nopal pueden proveer un hábitat similar al de las nopaleras naturales, predecimos que su composición de especies será similar , lo que contribuiría a mantener la homogeneidad biótica en la región.

Materiales y métodos

El área de estudio se ubica en la región del Altiplano en el estado de Hidalgo, la cual se constituye de amplios valles seccionados por cadenas montañosas volcánicas (^{Morrone, 2001}). La vegetación dominante es matorral xerófilo y el clima es semiseco templado, con precipitación anual promedio de 350.80 mm y temperatura promedio de 17.15 °C (^{Inegi, 2010}). La zona de estudio incluyó los municipios de Actopan, Singuilucan, Pachuca, Zempoala, Epazoyucan y Mineral de la Reforma, todos ellos con suelos de origen volcánico (^{Inegi, 2010}). Los cultivos de nopal estudiados pertenecen a la variedad de tuna blanca de la especie Opuntia albicarpa Scheinvar cv. “Reyna”, registrada como NOP-020-221104 en el Catálogo Nacional de Variedades Vegetales de la Sagarpa. En los sitios de nopaleras naturales se identificaron 4 especies: Opuntia streptacantha, O. heliabravoana, O. spinulifera y O. robusta.

El muestreo se realizó en 13 sitios: 8 cultivos de nopal tunero y 5 nopaleras naturales. Los sitios se ubicaron a una distancia mínima de 1 km entre ellos. Los muestreos se realizaron una vez al mes en cada sitio, durante 4 meses (entre mayo y agosto del 2014). En cada muestreo se utilizaron 9 trampas de caída por sitio, separadas 100 metros, las cuales consistieron en recipientes de 20 cm de altura y 10 cm de diámetro, llenos a un 20 % de su capacidad con una mezcla de alcohol al 70%, etilenglicol para evitar la evaporación y jabón para romper la tensión superficial, sin ningún cebo como atrayente. Las trampas fueron enterradas a nivel del suelo y permanecieron abiertas por 48 horas. Cada muestra fue almacenada en etanol al 70% para su posterior montaje en seco. La recolecta de ejemplares se realizó con el permiso Semarnat-FAUT 0254. La identificación se llevó a cabo con las claves de ^{Cuezzo y Guerrero (2012)}, ^{Mackay y Mackay (2005)}, ^{MacGown et al. (2014)} y ^{Snelling y Snelling (2007)}, y con la base de datos AntWeb (^{The California Academy of Sciences, 2012}). Ejemplares de referencia se depositaron en la colección entomológica del Instituto de Ecología, A.C. (IEXA; Reg. SEMARNAT: Ver. IN.048.0198).

Se midieron 4 características del hábitat: volumen de las plantas, complejidad vertical, número de flores y número de frutos, tanto en los cultivos como en las nopaleras naturales, para explorar su relación con los atributos de las comunidades de hormigas. A continuación se describe brevemente la metodología realizada para estas 4 características del hábitat, pero una descripción más amplia de los 2 tipos de hábitat se incluye en el Material Suplementario. El volumen (ancho × largo × altura) se midió en 9 plantas de nopal de cada sitio, separadas por una distancia de al menos 100 m. La complejidad vertical se refiere a la densidad vertical de la vegetación a diferentes alturas, la cual se estimó utilizando una barra de 2 m x 10 cm con rayas blancas y negras de 10 cm de ancho. La barra se colocó a 50 m del centro de cada sitio de muestreo, perpendicularmente a las plantas de nopal. Desde el punto central se estimó el porcentaje de rayas que eran visibles para un observador. Se utilizó el inverso de este porcentaje como el índice de la cubierta vertical (^{Griffith y Youtie, 1988}) y se usó el promedio obtenido de 10 medidas hechas a lo largo de 2 transectos de 100 m de longitud en cada sitio. Los números de flores y de frutos (tunas) se contaron en 2 trayectos lineales en banda (100 ( 1 m), desde el centro de cada uno de los sitios (^{Gardner, 1998}; ^{Rich et al., 2005}).

Para comparar la riqueza de especies acumulada entre las nopaleras naturales y los cultivos, se realizaron curvas de acumulación de especies y se calcularon intervalos de confianza al 95%. Debido a que en las nopaleras naturales el número de individuos fue bajo, se realizó una extrapolación de la riqueza de especies para predecir la curva de acumulación al doble del número de individuos encontrados. Además, se calculó el estimador de riqueza Chao1 para predecir la riqueza máxima y evaluar qué tan completos fueron los inventarios. Las curvas de acumulación y el estimador Chao1 se obtuvieron con el programa EstimateS 9.1.0. (^{Colwell, 2013}). Considerando la amplia diferencia en el número de individuos capturados entre ambientes, se comparó la riqueza total de especies estandarizando a un mismo esfuerzo de muestreo (medido como cobertura de la muestra, i.e., la proporción que representan los individuos de cada especie en la muestra, con respecto al número total de individuos), en lugar de un mismo número de individuos como se hace en la rarefacción tradicional (^{Chao y Jost, 2012}; ^{López-Mejía et al., 2017}). Esta estandarización por cobertura de la muestra, y sus intervalos de confianza al 95%, se obtuvieron en el programa iNEXT (^{Hsieh et al., 2013}).

Se comparó el promedio por sitio de la abundancia (número de individuos), la riqueza y la diversidad de especies entre las nopaleras naturales y los cultivos. La diversidad se midió con el exponencial del índice de Shannon, que corresponde a la serie de números de Hill cuando el parámetro q = 1, y cuyas unidades son números de especies efectivas (^{Jost, 2006}). Para las comparaciones se utilizaron pruebas de t, verificando previamente la normalidad de los datos con pruebas de Kolmogorov-Smirnov (p > 0.05) en el programa SigmaStat 3.5 (^{Systat Software, 2006}).

Para comparar la composición de especies entre sitios de nopaleras naturales y de cultivos de nopal, se realizó un análisis de ordenación no paramétrico multidimensional (NMDS) con el índice de Bray-Curtis. Este análisis permite ubicar la posición relativa de sitios de muestreo en un espacio de 2 dimensiones, de forma que la distancia entre los sitios es un reflejo de la disimilitud en su composición de especies. Posteriormente, para evaluar las diferencias estadísticas en la composición entre las nopaleras naturales y nopaleras cultivadas se realizó un análisis de similitud (Anosim). El Anosim es análogo a un Anova y se basa en el parámetro R, que mide la variación en la disimilitud de especies (1-índice de Bray-Curtis) entre grupos de sitios y la compara con la variación dentro de los grupos. En este trabajo, la hipótesis nula (un valor de R cercano a cero) es que no hay diferencias estadísticas en la composición de especies entre nopaleras naturales y cultivos de nopal. Tanto el NMDS como el Anosim se realizaron en el programa PC-ORD (^{McCune y Mefford, 1999}).

Para la identificación de especies indicadoras de cada tipo de hábitat se calculó el índice de valor indicador (IndVal) propuesto por ^{Dufrêne y Legendre (1997)}. Este índice se basa en la abundancia y la frecuencia relativa de cada especie en cada tipo de hábitat, y su valor (que va de 0 a 100) se compara con el valor esperado de una prueba de aleatorización de Monte Carlo para probar su significación estadística (α = 0.05). Para cada especie, el IndVal y su valor P, se calcularon en el programa PC-ORD (^{McCune y Mefford, 1999}).

Finalmente, para evaluar la influencia de las variables del hábitat en las comunidades de hormigas se usaron modelos lineares generalizados (GLM) simples tomando el volumen de las plantas, la complejidad vertical, el número de flores y el número de frutos como variables explicativas; la riqueza, abundancia, y diversidad de hormigas se consideraron como variables de respuesta. Se usó la distribución Poisson para la riqueza por ser una variable discreta (conteos), la distribución binomial para la abundancia con el fin controlar la sobredispersión de los datos en una variable discreta, y la distribución gaussiana para la diversidad porque los datos son normales (^{Crawley, 2007}). Los GLM se realizaron en el programa R 3.3.2 (^{R Development Core Team, 2016}).

Resultados

Se recolectaron en total 12,839 hormigas (2,643 en nopaleras naturales y 10,196 en cultivos de nopal), pertenecientes a 25 especies de 17 géneros y 6 subfamilias (Tabla 1). Pogonomyrmex barbatus (Myrmicinae) fue la especie más abundante en ambos hábitats. En las nopaleras naturales se registraron 19 especies, lo cual representa el 97.56 de la riqueza estimada, mientras que en los cultivos de nopal se observaron 23 especies y el estimador Chao1 proporcionó el mismo número, por lo que el inventario está 100% completo de acuerdo a este estimador. Al extrapolar la curva de acumulación de especies de las nopaleras naturales, se observa que la riqueza no es distinta de la de cultivos (Fig. 1A). Además, al estandarizar ambos ambientes a un mismo esfuerzo de muestreo (cobertura de la muestra 1.00) la riqueza no difiere estadísticamente pues la extrapolación de las nopaleras naturales genera intervalos de confianza muy amplios (Fig. 1B). Tampoco hubo diferencias significativas en la riqueza (t < 0.01, g.l. 11, p = 1.00), abundancia (t = 0.91, g.l. = 11, p = 0.38) y la diversidad (t = -0.17, g.l. = 11, p = 0.86) promedio de los sitios entre nopaleras naturales y cultivos de nopal (Fig. 2). Cabe resaltar que la abundancia promedio es mucho menor en las nopaleras naturales, pero entre los sitios de cultivos presenta mucha mayor dispersión.

Figura 1 A) Riqueza de especies total acumulada y B) riqueza estandarizada a un mismo nivel de completitud del inventario (cobertura de muestra = 1.0) de hormigas epigeas en sitios de nopalera natural (matorral crasicaule dominado por Opuntia spp.) y cultivos de nopal tunero (Opuntia albicarpa). Las curvas de acumulación (A) son la interpolación (línea continua) y extrapolación (línea punteada) de la riqueza acumulada con intervalos de confianza al 95% (área sombreada). Las barras representan la riqueza estandarizada promedio (B) y las barras de error son intervalos de confianza del 95% obtenidos con 100 aleatorizaciones. 

Figura 2 A) Riqueza de especies, B) abundancia y C) diversidad promedio de hormigas epigeas en sitios de nopalera natural (n = 5) y cultivos de nopal tunero (n = 8). La diversidad es el exponencial del índice de Shannon (especies efectivas). Las barras de error son la desviación estándar. 

Diecisiete especies de hormigas (68%) se comparten en ambos hábitats, pero Camponotus pellarius, Formica sp., Myrmecocystus melanoticus, M. mexicanus y Leptothorax sp. se encontraron únicamente en las nopaleras cultivadas, mientras que Pseudomyrmex mexicanus y Temnothorax sp. sólo en las nopaleras naturales (Tabla 1). La disimilitud promedio entre pares de sitios fue de 0.73, con la mayor disimilitud (0.99) entre 2 sitios de cultivos y la menor (0.33) entre un sitio de cultivo y una nopalera natural. En el NMDS no se separan los sitios de diferentes ambientes por su composición de especies (Fig. 3), y el análisis de disimilitud mostró que no hay diferencias significativas en la composición de especies de hormigas entre las nopaleras naturales y los cultivos de nopal (Anosim R = -0.003, p = 0.53).

Figura 3 Ordenación no paramétrica multidimensional (NMDS) de los sitios en nopaleras naturales (puntos verdes) y cultivos de nopal (puntos anaranjados) con base en su disimilitud de especies de hormigas, con el índice de Bray-Curtis. 

Se encontraron únicamente 3 especies con valor indicador significativo (p < 0.05) para cultivos de nopal tunero: Dorymyrmex sp. (IndVal = 58.9), Liometopum apiculatum (IndVal = 55.9) y Tapinoma sessile (IndVal = 53.3). Por otro lado, en las nopaleras naturales Camponotus pudorosus y Pheidole obtusospinosa tuvieron valores de IndVal altos, aunque dichos valores no fueron significativos (Tabla 2).

En las nopaleras naturales ninguno de los parámetros de las comunidades de hormigas (riqueza, abundancia y diversidad) estuvo relacionado con las variables de estructura de la vegetación (volumen de las plantas y complejidad vertical) ni con la producción de flores y frutos (Tabla 3). Sin embargo, en los cultivos de nopal, la riqueza y la diversidad de especies sí se relacionan significativamente con el número de flores y el número de frutos, aunque no con la estructura de la vegetación (Tabla 3, Fig. 4).

Figura 4 Relaciones entre la diversidad de hormigas (medida con el exponencial del índice de Shannon, en número de especies efectivas) y variables de la estructura de la vegetación y producción de flores y frutos en cultivos de nopal tunero. Las líneas y los valores de p corresponden a modelos lineales generalizados. n.s.: p > 0.05. 

Discusión

Como se esperaba, la riqueza de hormigas en los 2 tipos de hábitat estudiados es alta (25 especies), en comparación con las 21 especies encontradas en sistemas agrícolas del valle del Mezquital (^{Hernández-Ruiz y Castaño-Meneses, 2006}) y las 26 especies encontradas en la zona árida del valle de Tehuacán (^{Ríos-Casanova et al., 2006}). Esta riqueza de especies es relevante tomando en cuenta que, a diferencia de nuestro estudio, en ambos trabajos referidos el muestreo abarcó toda la variación climática anual. En este trabajo, el muestreo realizado en los 2 tipos de hábitat alcanzó un buen nivel de completitud, por lo que las listas de especies son representativas de estos ambientes.

La riqueza de especies encontrada en los cultivos de nopal (23 especies) es mayor a la encontrada en otros estudios realizados en agroecosistemas (^{Hernández-Ruiz y Castaño-Meneses, 2006}; ^{Perfecto y Vandermeer, 1996}; ^{Risch y Carroll, 1982}) y es similar a la encontrada en las nopaleras naturales. No siempre se encuentra una relación directa entre la intensidad de manejo humano y la riqueza de especies de hormigas, y en consecuencia un ambiente modificado puede representar una alta o baja diversidad en comparación con el ecosistema nativo (^{Chanatásig-Vaca et al., 2011}). Nuestros resultados muestran una clara semejanza en la riqueza, abundancia, diversidad y composición de especies de hormigas entre los cultivos de nopal y las nopaleras naturales. Esta semejanza corrobora lo reportado en otros estudios que indican que algunos agroecosistemas pueden ser reservorios de la biodiversidad de hormigas (e.g., ^{Chanatásig-Vaca et al., 2011}; ^{Valenzuela-González et al., 2008}). En nuestro caso, es posible que las condiciones ambientales en los cultivos de nopal contribuyan a mantener la homogeneidad biótica entre estos agroecosistemas y los remanentes del ecosistema nativo.

La mayor abundancia y riqueza de hormigas en nuestro estudio corresponde a las subfamilias Myrmicinae y Formicinae, que son las mejor representadas a nivel nacional (^{Guzmán-Mendoza et al., 2010}; ^{Ríos-Casanova, 2014}; ^{Ríos-Casanova et al., 2004}; ^{Rojas-Fernández, 2001}; ^{Vásquez-Bolaños, 2015}). Los géneros Pogonomyrmex, Dorymyrmex y Pheidole, al igual que en los trabajos de ^{Hernández-Ruíz y Castaño-Meneses (2006)} y ^{Guzmán-Mendoza et al. (2010)}, fueron los más abundantes. La especie Liometopum apiculatum es una de las indicadoras de los cultivos de nopal y tiene una alta importancia por el valor gastronómico de sus larvas (escamoles) en la región (^{Cruz-Labana et al., 2014}). Además, esta especie se ha reportado que visita frecuentemente a las plantas de Opuntia sp. dando protección a los nopales contra herbívoros y depredadores (^{Miller, 2007}). También Tapinoma sessile resultó ser una especie indicadora de los cultivos de nopal, y este género ya había sido propuesto como indicador de agroecosistemas con poco disturbio del valle del Mezquital (^{Hernández-Ruiz y Castaño-Meneses, 2006}). Por lo tanto, es una especie con alta frecuencia de aparición y abundancia relativa en sitios con las condiciones ambientales de los cultivos de nopal. Será muy interesante evaluar cuáles variables del ambiente y del paisaje influyen en la permanencia de las especies en estos agroecosistemas.

La alta riqueza de especies de hormigas en los cultivos de nopal puede estar relacionada con la variación micro-ambiental generada por la complejidad estructural de la vegetación (^{Hernández-Ruiz y Castaño-Meneses, 2006}). Por ello, en ambientes modificados el cambio en la composición de la vegetación puede disminuir, o no afectar, la estructura y composición de las comunidades de hormigas (^{Perfecto y Snelling, 1995}). Nuestros resultados indican que la complejidad estructural de la vegetación (volumen de las plantas de Opuntia y complejidad vertical) no afecta a las comunidades de hormigas de los cultivos ni de las nopaleras naturales. Sin embargo, en los cultivos de nopal el número de flores y frutos de Opuntia sí se relaciona con la riqueza y la diversidad de hormigas. Estas relaciones podrían deberse a la protección que dan algunas especies de hormigas a las plantas de Opuntia contra hervíboros (^{Miller, 2007}). Sin embargo, esta hipótesis requerirá un diseño experimental apropiado para su estudio. Cabe resal-tar que nuestros muestreos se realizaron durante la época de floración y principio de la fructificación de los nopales, pero no hay datos sobre las comunidades de hormigas en los cultivos de nopal durante el resto del año.

En la zona árida de Tehuacán, en el centro de México, la distribución y abundancia de las especies dominantes se relaciona con características particulares de los suelos, como el porcentaje de arena (^{Ríos-Casanova et al., 2006}). Por lo tanto, para la protección de las hormigas epigeas y otros elementos de la biodiversidad del suelo, así como de los procesos del ecosistema, será necesario evaluar el impacto de los agroquímicos que se utilizan en los cultivos de nopal tunero. Un esquema de manejo sustentable y orgánico, libre de plaguicidas, herbicidas y fertilizantes químicos favorecerá la capacidad de estos cultivos perennes para conservar la biodiversidad. Además, será necesario evaluar los posibles costos y beneficios que la alta diversidad de hormigas podría tener para la producción de tunas. Por ejemplo, la herbivoría podría causar daños severos en las raíces y tallos de Opuntia, pero la presencia de las hormigas también podría beneficiar la productividad al dar protección contra otros herbívoros, frugívoros o robadores de néctar en estos cultivos.