Introducción

Amphipterygium adstringens (Schltdl.) Schiede ex Standl. o cuachalalate, es una de las especies de mayor importancia en la herbolaria mexicana, actualmente su corteza es utilizada para aliviar más de 30 enfermedades, entre ellas la malaria, cáncer de estómago, ulcera gástrica y enfermedades renales (^{Casas, Blancas et al., 2016}; ^{Cuevas, 2005}; ^{Osorio-García et al., 2018}; ^{Solares et al., 2012}). Amphipterygium adstringens ha sido catalogada como una especie altamente vulnerable de acuerdo con la amplitud de su área de distribución, el tipo de vegetación en que se desarrolla, las características de su ciclo de vida y la gran cantidad de usos que se le da (^{Beltrán-Rodríguez et al., 2017}; ^{Casas, Lira et al., 2016}).

^{Solares et al. (2012)} indicaban que en la región centro sur de México se recolecta un aproximado de 57.5 toneladas anuales de corteza de A. adstringens. El descortezamiento tradicional puede ser letal, se ha descrito que 60% de los árboles llega a morir, ya que se afecta el cambium vascular y el floema del individuo ocasionando su muerte (^{Solares y Gálvez, 2002}; ^{Solares et al., 2006}). En un estudio con descortezamiento experimental en A. adstringens, se encontró una mortalidad de 13.9% en 2 años, la cual está asociada a la época en que se realiza el descortezamiento y la cantidad de corteza extraída (^{Beltrán-Rodríguez et al., 2021}). Por otro lado, el cuachalalate presenta un alto porcentaje de frutos sin semilla, y el número de semillas por fruto es muy variable (^{Guzmán-Pozos y Cruz-Cruz, 2014}), sugiriendo que la reproducción sexual es una limitante en el mantenimiento de las poblaciones.

En el municipio de Chietla, Puebla, la baja abundancia relativa de individuos adultos (19 a 56 individuos/ha), indica que las poblaciones desaparecerán en el corto plazo si no se realizan acciones de conservación (^{Martínez, 2011}). Además de la demanda regional y nacional, la explotación del cuachalalate ha aumentado en los últimos años debido a la demanda de empresas transnacionales, a partir de que se descubrió un compuesto con alto potencial en el tratamiento de adenocarcinoma (^{Solares et al., 2012}). Tal explotación, sumada al cambio de uso de suelo y el cambio climático, ha provocado la desaparición de poblaciones en diversos sitios y/o la disminución del número de individuos y de su diversidad genética, además, se ha proyectado un decremento de 10% para las siguientes 3 generaciones (^{Guzmán-Pozos y Cruz-Cruz, 2014}; ^{Martínez et al., 2020}; ^{Solares et al., 2006}).

A pesar de su importancia, A. adstringens no se encuentra protegida por las leyes nacionales (^{DOF, 2019}), solo se ha catalogado como especie vulnerable internacional- mente por la IUCN (^{Martínez et al., 2020}), aunque no se conoce el estatus de sus poblaciones. Particularmente en aquellas especies silvestres que se encuentran bajo explotación, el conocimiento sobre su dinámica poblacional y ecología es imperativo (^{Beltrán-Rodríguez et al., 2017}). El análisis de los parámetros poblacionales ayuda a establecer estrategias de conservación, uso y manejo para especies en riesgo (^{Krebs y Elwood, 2008}; ^{Mollers, 2008}). Una población se caracteriza por la distribución de sus individuos y la abundancia de estos. El tamaño poblacional, el patrón de distribución espacial, la estructura de tamaños/edades, la proporción de sexos, las tasas de natalidad y mortalidad, y la tasa finita de crecimiento, son indicadores del estado de la población (^{Caswell, 1986}; ^{Mollers, 2008}).

El patrón de distribución espacial es particularmente importante cuando se pretende dar un manejo a la población o recuperar las poblaciones naturales. Cada especie presenta un patrón particular en el espacio, que está determinado por factores bióticos y abióticos, tales como la dispersión de semillas, distribución de las plantas parentales, polinización, depredación, presencia y ubicación de especies facilitadoras, condiciones climáticas y microclimáticas (^{Sousa y Cunha, 2018}). Un patrón al azar ocurre en ambientes en donde las condiciones son homogéneas y hay la misma probabilidad de encontrar a los organismos en cualquier lugar, un patrón uniforme es indicativo de que los organismos se separan de manera consistente a distancias mayores de lo esperado en una distribución azarosa, es decir, la distancia entre vecinos se maximiza, mientras que un patrón agregado muestra que los organismos se encuentran en un ambiente heterogéneo y se establecen principalmente bajo ciertas condiciones (^{Fortin y Dale, 2005}).

La estructura de edades permite conocer el estado de una población, puede indicar periodos reproductivos exitosos, periodos de baja o alta supervivencia, y puede predecir si una vez que los individuos adultos mueran, habrá o no un reemplazamiento. En una población estable la cantidad de organismos pre-reproductivos será muy similar a la de aquellos reproductivos, asegurando el reemplazamiento de individuos en el tiempo; mientras que en una población en crecimiento los organismos prereproductivos serán más abundantes que los reproductivos, incrementando el número poblacional al paso del tiempo; finalmente, una población se considera decreciente cuando no hay producción de nuevos individuos o ésta es muy baja, de manera que a largo plazo, se pone en peligro la existencia de la población (^{Begon et al., 2006}; ^{Elzinga et al., 2009}; ^{Mollers, 2008}).

El monitoreo de la edad (en años) para especies caducifolias arbóreas es complicado por el crecimiento lento y reducido que ocurre particularmente durante la época de lluvias (^{Brienen et al., 2016}; ^{Verdú y García-Fayos, 2002}). Sin embargo, diversos estudios han demostrado que el diámetro a la altura del pecho (DAP) de varias especies se correlaciona positiva y significativamente con la cantidad de anillos de crecimiento, permitiendo establecer estrategias de manejo, y teniendo al DAP como una medida indirecta de la edad (^{Brienen et al., 2016}; ^{Paredes-Villanueva et al., 2013}; ^{Tolera et al., 2013}, ^{Worbes et al., 2003}).

Las tasas de natalidad, mortalidad y crecimiento son más difíciles de determinar, ya que requieren de estudios a mediano y largo plazo para evitar interpretaciones erróneas, por ejemplo, en una especie arbórea y caducifolia como A. adstringens, la producción de semillas, germinación, establecimiento y crecimiento estará determinado por factores bióticos y abióticos, cuya variación anual es difícil de predecir.

El objetivo de este trabajo fue determinar el estado actual de las poblaciones de A. adstringens en 10 localidades a lo largo de su distribución en México. En cada población se evaluó, la densidad de individuos, el patrón de distribución espacial, la estructura de edades (reproductiva y pre-reproductiva), densidad de individuos por edad y su distribución espacial, así como la posible asociación con especies nodriza y la evidencia de extracción de corteza. De acuerdo con los antecedentes de explotación y estatus de vulnerabilidad de A. adstringens, se hipotetiza que las poblaciones de cuachalalate tendrán una baja densidad de individuos y una estructura de edades decreciente, esto es, una mayor proporción de adultos en comparación con los estadios jóvenes y, por tanto, un bajo establecimiento de nuevos individuos, lo que podría poner en riesgo a las poblaciones.

Materiales y métodos

Amphipterygium adstringens (Anacardiaceae) es un árbol endémico de México, su distribución se restringe al bosque tropical caducifolio (BTC) de Sinaloa, Nayarit, Jalisco, Colima, Michoacán, Morelos, Estado de México, Puebla, Guerrero, Oaxaca y Chiapas (^{Cuevas, 2005}; ^{Guzmán-Pozos y Cruz-Cruz, 2014}), a una altitud entre 5 y 1,500 m (^{Cuevas, 2005}). En Morelos, la especie es considerada un elemento abundante del bosque, tanto por su densidad como por su tamaño; se han reportado individuos entre 2 y 16 m de altura con un DAP de hasta 35 cm (circunferencia a la altura del pecho o CAP aproximada de 110 cm) y hasta 12.85 m² de cobertura (^{Beltrán-Rodríguez et al., 2018}).

Amphipterygium adstringens es dioica, la floración ocurre durante la época de lluvias (usualmente de mayo a junio) y dura menos de 2 meses, mientras que la fructificación ocurre en la época seca y tiene una duración mayor a 5 meses, el pico de dispersión de semillas se presenta al inicio de la época seca (octubre) (^{Luna-Nieves et al., 2017}). El fruto es alado sugiriendo la dispersión por viento, presenta de 1 a 4 semillas, los frutos más grandes tienen mayor probabilidad de desarrollar semillas viables, las semillas abortadas son comunes (^{Guzmán-Pozos y Cruz-Cruz, 2014}).

El trabajo se realizó en los estados de Guerrero, Morelos, Michoacán, Puebla y Oaxaca durante la tempo- rada seca 2016-2017. Se estudiaron 2 localidades por estado (Tabla 1). Se obtuvo el apoyo de los pobladores quienes compartieron (de manera informal) sus conocimientos sobre la especie, tanto sobre su presencia en las localidades como del aprovechamiento en cada región. En la mayoría de los sitios se realizan actividades agropecuarias y ganaderas (cultivos de diferentes especies y pastoreo de ganado bovino, vacuno y caprino). Los terrenos son de propiedad ejidal, excepto en Michoacán y Oaxaca, en donde son de propiedad comunal.

Las localidades de Morelos (Tlayca y Tlayecac) se eligieron debido a que son sitios con remanentes de BTC fuera de la Reserva de la Biósfera Sierra de Huautla, en donde la especie se encuentra ampliamente estudiada. Particularmente, en Ostula y Pómaro (Michoacán), los sitios fueron propuestos por el Regidor de Ecología del Ayuntamiento del municipio de Aquila, con la finalidad de contribuir al poco conocimiento de la zona en materia de recursos naturales y que el equipo de trabajo se encontraría protegido ante la inseguridad social de zonas aledañas. En las localidades de Puebla, Oaxaca y Guerrero el equipo de trabajo ha realizado diferentes estudios por varios años lo que facilitó la obtención de permisos y permitió la observación previa. El uso en los sitios de estudio que se da al cuachalalate en general es local, aunque en Tlayca y Tlayecac, las personas indicaron que los árboles de esta especie y de Bursera han disminuido en los últimos años, ya que la gente de otros pueblos realiza el saqueo de corteza y resina de manera constante y eso provoca la muerte de los árboles, al grado de que ellos ya no pueden aprovechar nada.

En cada localidad se establecieron 4 unidades de muestreo de 30 × 40 m, únicamente se consideró que las unidades de muestreo no fueran parte de un camino o vereda y que los pobladores reconocieran el sitio como hábitat de A. adstringens. En cada unidad de muestreo se realizó una búsqueda intensiva de las plántulas, juveniles y adultos. No se registraron plántulas durante los censos, para juveniles se midió la altura, circunferencia a la altura de la base (CAB). Para los adultos se midió circunferencia a la altura del pecho (CAP_1.30m), altura y tamaño de la copa; el CAP se obtuvo de acuerdo con ^{Eshete et al. (2011)} y ^{Tolera et al. (2013)}.

Las edades de los individuos se asignaron apegándose en la medida de lo posible, a la clasificación de ^{Komarov et al. (2003)}, de acuerdo con la edad pre-reproductiva (plántulas y juveniles) y reproductiva (adultos, con evidencia de estructuras reproductivas). Particularmente, en organismos reproductivos de A. adstringens no hay reportes sobre los cambios que sufre el crecimiento vegetativo o la actividad reproductiva durante la vida de los individuos. Por lo tanto, para la estructura de edades también se consideró la altura y diámetro, de manera que se obtuvieron las siguientes edades: juvenil 1 o individuos pre-reproductivos de reciente establecimiento, sin cotiledones y sin un sistema de ramas laterales desarrollado, con altura < 50 cm y CAB < 10 cm; juvenil 2, individuos pre-reproductivos con altura de 50.1-130 cm y CAB < 10cm; adulto 1, individuos con CAP = 10.1-50 cm; adulto 2, individuos con CAP = 50.1-100 cm y adulto 3, individuos con CAP > 100.1 cm (modificado de [^{Ponce [2013}]).

La cantidad de individuos registrados en las localidades se comparó con una prueba de χ ² , la densidad absoluta de individuos se obtuvo por unidad de muestreo y se expresó como densidad promedio por hectárea para cada localidad.

Para determinar el patrón de distribución espacial se utilizó el índice estandarizado de Morisita (IEM). Este índice se ha utilizado ampliamente para diferentes organismos y no es sensible al tamaño de la unidad de muestreo, ni al cambio en la densidad de individuos (^{Dalmagro y Vieira, 2005}; ^{García-González et al., 2016}; ^{García-Ibáñez et al., 2007}; ^{Marín et al., 2009}; ^{Wattanachaiyingcharoen et al., 2008}); el análisis se hizo con el programa Past 4.03 (^{Hammer et al., 2001}). Los datos utilizados consisten en una serie de recuentos de puntos en cuadrantes o unidades muestreo del mismo tamaño. El índice estandarizado de Morisita (IEM) fue sugerido por ^{Smith-Gill (1975)}, como una modificación al índice de Morisita I _δ (1959), y es calculado por las siguientes fórmulas:

En donde n es el número de cuadrantes y x son los conteos. La prueba de significancia que sigue al índice de Morisita Iδ dada por la razón de F consiste en:

Los grados de libertad son n-1 y ∞, adicionalmente se incluye una prueba de Monte Carlo con 9999 repeticiones, cada una con una distribución al azar de puntos en los cuadrantes. Posteriormente se calculan los límites de confianza superior e inferior (a 95%) con respecto al índice de Morisita I _δ . El límite inferior también es llamado índice uniforme (M _u ), mientras que el límite superior es conocido como índice agrupado (M _c ).

En donde es el percentil 97.5 de la distribución chicuadrada con n-1 grados de libertad. El índice estandarizado de Morisita (IEM) se calcula como sigue:

Cuando el índice de Morisita es mayor que el límite agrupado (I _δ >M _u ) el patrón de distribución espacial es agregado. En cambio, cuando el índice de Morisita es menor que el límite uniforme (I _δ < M _c ), el patrón es uniforme. El índice estandarizado de Morisita (IEM) indica una distribución al azar cuando el valor obtenido se encuentra entre -0.5 y 0.5, una distribución uniforme cuando es menor que -0.5 y un patrón agregado cuando es mayor que 0.5 (^{Hammer et al., 2001}; ^{Krebs, 1999}; ^{Sousa y Cunha, 2018}).

La distribución espacial se analizó primero considerando el número de individuos en cada sitio. Posteriormente, también se analizó la distribución de adultos y de individuos pre-reproductivos por separado, en el análisis se omitieron las localidades con menos de 3 individuos pre-reproductivos.

La cantidad de individuos registrados por categoría de edad se comparó con una prueba de χ ² . La estructura de edades se analizó para cada localidad cuantitativa y cualitativamente a través de la densidad y mediante histogramas horizontales considerando que en la base se coloca el estadio más joven y por encima de éste los estadios sucesivos. Una estructura con forma triangular es considerada una población en crecimiento, ya que cuenta con una gran proporción de individuos jóvenes e indica alta tasa de natalidad; una población con forma de campana es una población estable, ya que la proporción de juveniles y adultos es similar e indica tasas de natalidad y mortalidad similares. Finalmente, una población decreciente tendrá forma de urna o bulbo, con estadios jóvenes reducidos en comparación con la proporción de adultos, este tipo de estructura indica baja tasa de natalidad o establecimiento de nuevos individuos y pronostica la desaparición de la población si los estadios adultos desaparecen (^{Krebs, 2009}; ^{Smith y Smith, 2007}).

La facilitación en zonas áridas y en el BTC es un proceso clave en el reclutamiento de nuevos individuos, por lo que consideramos importante evaluar estos datos (^{Contreras, 2015}; ^{González, 2012}, ²⁰¹³; ^{Khurana y Singh, 2004}; ^{Schleicher et al., 2011}; ^{Valiente-Banuet y Verdú, 2007}). Cuando se encontraron individuos pre-reproductivos bajo la copa de plantas perennes adultas, se registró la orientación del individuo con respecto al centro de la copa considerando 4 puntos cardinales, las plantas adultas fueron identificadas a nivel de especie, se midió la altura, CAP_1.30m, tamaño de la copa y se registró su forma de crecimiento según ^{Judd et al. (2002)}. Para estimar el tamaño de la copa se midieron 2 diámetros perpendiculares y el tamaño se calculó por medio del área de una elipse.

Se utilizó una prueba de t para determinar si había diferencias entre la cantidad de individuos pre-reproductivos creciendo a cielo abierto y aquellos que crecían bajo la copa de plantas adultas. Se analizó la riqueza de especies de plantas bajo las que se encontraron los individuos pre-reproductivos. Para determinar cuántos organismos asociados había por cada individuo adulto, se calculó la razón de asociación por especie nodriza dividiendo el número de individuos pre-reproductivos de A. adstringens en la unidad de muestreo, entre el total de individuos de la especie nodriza con registros asociados. Se utilizaron correlaciones de Spearman para determinar si el número de individuos asociados tenía alguna relación con el tamaño de la copa y con el número de plantas adultas, considerando la totalidad de éstas, así como las especies arbóreas y arbustivas por separado. Se utilizó la correlación de Pearson para determinar si había una relación entre el número de individuos pre-reproductivos creciendo bajo conespecíficos y el número de A. adstringens adultos.

Finalmente, se analizó si existían diferencias significativas entre el número de individuos de A. adstringens establecidos bajo las copas considerando las diferentes orientaciones mediante pruebas de χ ² , bajo el supuesto de que existe la misma probabilidad de establecimiento hacia cualquier punto cardinal. Para definir qué orientación era significativa, se hizo un análisis de residuales estandarizados cuya fórmula es R=(O-E)/√E, en donde O son los datos observados y E son los datos esperados; una orientación se consideró significativa para un valor residual absoluto R ≥ 1.96 (^{Sheskin, 2011}). Los análisis se realizaron con GraphPad Prism 7, previa prueba de normalidad, con error alfa de 0.05.

En todos los sitios se buscaron los árboles con evidencia de extracción de corteza considerando que al retirar ésta, se observa el xilema secundario desnudo (cuando la extracción es reciente), o el árbol produce nueva corteza y su desarrollo es evidente por diferencias en el grosor y coloración de la peridermis. Para determinar si la cantidad de árboles presentes en los sitios está correlacionada con la cantidad de árboles sometidos a extracción de corteza, se utilizó una correlación de Pearson previa prueba de normalidad con error alfa de 0.05 en el programa GraphPad Prism 7. Los porcentajes de extracción se analizaron con respecto a la edad de los árboles utilizados.

Resultados

Se registraron 839 individuos en total en las unidades de muestreo. El número de individuos varió entre localidades (χ ² = 444.34, g.l. = 27, p < 0.0001). La mayor densidad de individuos se registró en Puebla (Tabla 2), en las localidades de Coxcatlán (279 ± 163 individuos/ha) y Tilapa (265 ± 28 individuos/ha), seguido de Cascalote en Oaxaca (250 ± 84 individuos/ha) y Pómaro en Michoacán (242 ± 94 individuos/ha). Las localidades con menor densidad fueron Tlayca y Tlayecac en el estado de Morelos (50 ± 11 individuos/ha y 77 ± 40 individuos/ha, respectivamente). En algunas unidades se registró 1 solo individuo y en el caso de Tlayecac, una de las unidades no presentó individuos.

De los 839 individuos registrados, 27.2% eran juveniles de reciente establecimiento (juvenil 1), 4.9% eran individuos pre-reproductivos juvenil 2; 44.6% adultos 1, 19.6% adultos 2 y 3.7% adultos 3. Esta diferencia entre la cantidad de individuos registrados en cada categoría de edad fue significativa (χ ² = 422.9, g.l. = 36, p < 0.0001), de manera que 2/3 de los organismos fueron individuos reproductivos (67.9%). La alometría promedio de los individuos se presenta en la Tabla 3. En ninguno de los sitios se registraron plántulas.

En Michoacán no se encontraron organismos de la categoría juvenil 1 (Tabla 4), la menor densidad promedio se registró en las 2 localidades de Morelos (4 ± 2 individuos/ ha en Tlayca y 6 ± 4 individuos/ha en Tlayecac), mientras que la mayor densidad fue en Cascalote, Oaxaca (146 ± 76 individuos/ha) y Tilapa, Puebla (115 ± 47 individuos/ ha). En cuanto a los individuos de la categoría juvenil 2, en 5 localidades (Ostula, Pómaro, Cuambio, Tlayca y Tlayecac) los individuos se encontraron en una sola unidad de muestreo, mientras que en Papalutla y Cascalote no hubo registros. La densidad promedio más alta se registró en Coxcatlán en Puebla (31 ± 20 individuos/ha), seguida del cañón el Sabino en Oaxaca y Tilapa en Puebla (21 ± 12 y 15 ± 10 individuos/ha, respectivamente).

A su vez, los organismos de la categoría adulto 1 estuvieron bien representados en todas las localidades, con una densidad mínima de 31 ± 14 individuos/ha en Cuambio en Guerrero, y la máxima en Tilapa en Puebla (113 ± 35 individuos/ha), otras localidades con alta densidad fueron Papalutla, Coxcatlán y el cañón el Sabino, sin embargo, también presentaron una fuerte variación entre los resultados por unidad de muestreo. Los registros de la categoría adulto 2 en general presentaron bajas densidades (0 a 5 individuos/ha), excepto en Pómaro en Michoacán (19 ± 7 individuos/ha). Finalmente, la densidad de adultos 3 fue nula en 5 localidades, mientras que en Coxcatlán en Puebla y Ostula en Michoacán solo se encontraron

En Ostula se localizaron 2 individuos pequeños (1.5 m y 1.54 m de altura, ambos con 9 cm de CAB) con estructuras reproductivas, no presentaban cobertura, el ápice estaba representado por una rama delgada sin ramificación, el primero presentaba 3 botones florales y el segundo, 2 flores en antesis, sugiriendo que estas características de tamaño ya son un primer indicio de la etapa adulta para la especie, al menos en esa región.

El patrón de distribución espacial de los individuos (pre-reproductivos y reproductivos) en las poblaciones fue al azar en Cuambio (Guerrero), Tilapa (Puebla) y Tlayca (Morelos), en el resto de las localidades se encontró un patrón agregado (Tabla 2). Sin embargo, al analizar por separado la distribución de adultos e individuos prereproductivos asociados y no asociados, únicamente en Cuambio y Tlayca se conservó el patrón, en el resto de los sitios se observaron cambios.

Número de individuos pre-reproductivos y reproductivos, densidad promedio por hectárea (± error estándar), índice estandarizado de Morisita (IEM) y respectivo patrón de distribución espacial de A. adstringens registrado por localidad.

El conjunto de individuos adultos por localidad mostró una distribución al azar únicamente en Cuambio (IEM = -0.43, n = 34 individuos) y Tlayca (IEM = -0.03, n = 21 individuos), en el resto de las localidades se observó una distribución agregada (IEM = 0.51 a 0.82, n = 33 en Tlayecac, n = 50 en El Cascalote, n = 53 en Ostula, n = 54 en el cañón el Sabino, n = 65 en Tilapa, n =72 en Papalutla, n =73 en Coxcatlán, n =115 en Pómaro). En cuanto a los individuos pre-reproductivos, 56% se registró creciendo a cielo abierto y 44% bajo la sombra de alguna planta adulta (n = 250).

Los organismos no asociados a plantas (n = 140 individuos totales, 98.2% juvenil 1 y 1.8% juvenil 2), mostraron una distribución agregada en Tilapa, Coxcatlán, cañón el Sabino y El Cascalote (IEM = 0.54 a 0.62), mientras que en Ostula, Cuambio, Tlayca y Tlayecac se encontró una distribución al azar (IEM = -0.35 a 0.47); en Papalutla y Pómaro solo se registró 1 individuo. Los organismos asociados a plantas adultas (n =110, 75.7% juvenil 1 y 24.3% juvenil 2) presentaron una distribución agregada en Tilapa, Coxcatlán y El Cascalote (IEM = 0.53 a 0.69), en Cuambio y en el cañón el Sabino la distribución fue al azar (IEM = 0 a 0.08), Papalutla solo tuvo un registro y finalmente, Ostula, Pómaro, Tlayca y Tlayecac no tuvieron registros.

No se encontraron poblaciones en crecimiento o estables, en cada localidad la estructura poblacional es decreciente ya que la proporción de adultos es mucho mayor que la de individuos pre-reproductivos (Fig. 1). La proporción de juvenil 1 fue de 5 a 8% en Papalutla y las 2 localidades de Morelos, y de 30 a 34% en Características alométricas de los individuos de A. adstringens. Altura, circunferencia a la altura de la base (CAB) o el pecho (CAP) y cobertura de la copa. Datos tomados en 10 localidades a lo largo de la distribución de la especie; se presenta el valor promedio ± error estándar.

Figura 1 Estructuras de edad de A. adstringens en los sitios de estudio. Los números representan el porcentaje de individuos en cada categoría de edad. 

En cuanto a los organismos adultos, el patrón general indica que dominan los individuos adulto 1 y adulto 2, mientras que los adulto 3 se encuentran en baja proporción o no existen. En la mitad de los sitios (Ostula, Papalutla, Tlayca, Tlayecac y el cañón el Sabino), la proporción de adulto 1 estuvo por arriba de 52%; en Pómaro la mayor proporción fue de adulto 2 (64%); Cuambio presentó proporciones decrecientes a medida que se observan árboles más grandes (28% en adulto 1, 21% en adulto 2 y 15% en adulto 3), en Ostula se observó un patrón similar, sin embargo, la diferencia entre las categorías de adulto fue más grande que en Cuambio (64% en adulto 1, 25% en adulto 2 y 5% en adulto 3). Salvo por la baja proporción de juvenil 2, la estructura poblacional de Cuambio, Guerrero podría considerarse estable.

La proporción de individuos asociados y no asociados mostró un comportamiento distinto en cada localidad, en Cascalote se encontró la misma proporción de individuos asociados que no asociados. En Coxcatlán y el cañón el Sabino la mayor proporción fue de individuos no asociados (62.3% y 72.4%, respectivamente), mientras que en Cuambio y Tilapa la mayor proporción fue de individuos asociados (53.3% y 56.5%, respectivamente). Sin embargo, estas diferencias no fueron suficientes para establecer un patrón general de asociación (t = 0.92, df = 4, p = 0.4).

Los individuos asociados se encontraron principalmente bajo especies arbóreas (77.1%) y en menor proporción bajo especies arbustivas (22.9%); solo se encontró 1 individuo asociado a una roca. La mayor proporción de organismos asociados se encontró en Cascalote y Tilapa (31.8% en cada localidad), seguido de Coxcatlán (21%). Se registró un total de 54 individuos de 20 especies de plantas asociadas, pertenecientes a 8 familias, la más común fue Fabaceae con 6 especies, seguida de Cactaceae con 5, y Euphorbiaceae y Burseraceae con 2 especies cada una (Tabla 5).

Los sitios con mayor riqueza de plantas asociadas fueron El Cascalote (10 especies), Coxcatlán (8 especies) y Tilapa (6 especies). En 31.2% de los casos, la razón de asociación fue igual a 1, es decir, 1 organismo asociado por cada individuo adulto, en el resto la razón de asociación varió entre 1.3 y 6.5 individuos por planta adulta. En varios sitios destacan algunas especies, tal es el caso de Tilapa en donde Ceiba parvifolia, Fouquieria formosa y Parkinsonia praecox (razón de asociación igual a 6.5, 5.5 y 5, respectivamente), son las especies con mayor número de individuos asociados. Se encontró que los individuos no reproductivos de A. adstringens se establecen bajo la cobertura de conespecíficos adultos en 5 de 6 localidades (en Papalutla solo hubo un registro bajo Lysiloma sp.), además, la razón de asociación fue de 1.5 a 3 asociados por adulto, sugiriendo que no es una relación azarosa. Especies como Mimosa luisana y Parkinsonia praecox, también resultaron con asociaciones importantes en las localidades de Puebla y Oaxaca.

Se observó también que el número de individuos asociados no tiene relación con el tamaño de la copa de las plantas adultas (r_s = 0.07, p = 0.77, n = 20), sino con el número de plantas adultas presentes (r_s = 0.86, p < 0.0001, n = 20), incluso dicha asociación es significativa cuando se considera solo la cantidad de especies arbóreas (r_s = 0.89, p < 0.0001, n = 14), pero se pierde con las especies arbustivas (r_s = 0.58, p = 0.24, n = 6). Excepto en Tilapa, donde se encontraron las copas más grandes, se observó que la amplitud de la copa no es relevante en el número de asociaciones, ya que la mayor cantidad de asociados se encuentra en donde hay mayor cantidad de adultos.

Por otro lado, también se observó que en todas las localidades hay al menos 1 especie arbustiva importante en la interacción (Fig. 2). En todas las localidades se encontró que la asociación de A. adstringens y sus conespecíficos adultos aumenta cuando hay más adultos, aunque la relación no fue significativa (r = 0.86, p = 0.057, n = 5). Estos resultados sugieren que la asociación está determinada por patrones de dispersión, condiciones de facilitación y características del desarrollo de A. adstringens.

Figura 2 Relación entre el número de individuos asociados de A. adstringens (triángulos), el tamaño de la copa (área azul) y el número de plantas adultas (barras) en los sitios de estudio. Las barras vacías indican especies arbóreas y las barras llenas especies arbustivas. El tamaño de la copa representa el área promedio de las plantas adultas. 

La orientación de individuos establecidos bajo la copa mostró diferencias significativas (χ ² = 65.86, df = 3, p < 0.0001), indicando que hay una mayor cantidad de individuos en el oeste, el norte y el este (R > |2| para cada orientación) con respecto al sur. El 78% del total de los datos se concentró en 5 especies (A. adstringens, C. parvifolia, M. luisana, F. formosa y P. praecox) en las localidades de Puebla y Oaxaca; para este porcentaje se observó el mismo patrón de orientación con mayor cantidad de individuos en el oeste, el norte y el este (R > |2| para cada orientación), mientras que los individuos en el sur no fueron significativos (χ ² = 45.67, df = 3, p < 0.0001). Los individuos asociados a A. adstringens (28% del total de los datos), mostraron una tendencia a establecerse al oeste (R > |2|) de la planta adulta (χ ² = 15.71, df = 3, p = 0.0013).

Las marcas de extracción de corteza se localizaron en 37 individuos adultos en las localidades de Guerrero y Puebla y en el cañón el Sabino (Oaxaca). Se observó una relación negativa pero no significativa entre el total de individuos adultos y el número de organismos con marcas (r = -0.63, p = 0.25, n = 5). Se destaca que en Cuambio 35.3% se presentó alguna marca de extracción, esta es la población con menor densidad de adultos y la proporción de adulto 1 es el doble de las 2 categorías restantes; sin embargo, las marcas se encontraron en árboles adulto 2 (14.7%) y adulto 3 (20.6%). En Coxcatlán solo se encontraron marcas en 12.5% de los adultos, 2/3 partes eran adulto 2 y 1/3 era adulto 1, sin embargo, ésta es una de las poblaciones con mayor densidad de adultos y solo 1 individuo en categoría adulto 3. En el cañón el Sabino se registró la extracción en 18.9% de los adultos, en este sitio la densidad de adulto 2 es muy baja aun cuando los organismos registraron marcas (3.8%), el restante 15.1% fue localizado en adulto 1; en este sitio no había adulto 3. Estos resultados sugieren que los pobladores tienden a utilizar la corteza de los árboles con mayor tamaño.

Discusión

Amphipterygium adstringens es un árbol dominante en algunas localidades de BTC de México (^{Beltrán-Rodríguez, 2018}; ^{Cortés-Flores et al., 2017}; ^{Luis, 2017}; ^{Maldonado et al., 2013}; ^{Méndez-Toribio et al., 2014}; ^{Sánchez et al., 2018}). En el presente estudio encontramos variaciones en la densidad de individuos entre los sitios. Los valores más altos de densidad se encontraron en Tilapa, Coxcatlán, Pómaro y El Cascalote, seguido del cañón el Sabino, tales valores pueden deberse a la baja perturbación y las pocas actividades antrópicas que existen en estos sitios, ya que 4 de ellos se encuentran en un área destinada a la conservación de especies dentro de la Reserva de la Biósfera Tehuacán-Cuicatlán (^{Conanp, 2021}).

Las densidades más bajas encontradas fueron en los sitios Tlayca y Tlayecac, similares o mayores a las reportadas para Chietla, Puebla (19 a 56 individuos/ha; ^{Martínez, 2011}) y para las descritas en 2 localidades de San Juan Bautista Jayacatlán, Oaxaca (54.9 y 29.6 individuos/ ha; ^{Zamora, 2003}). Las bajas densidades encontradas en el presente trabajo en Tlayca y Tlayecac pueden deberse a que estos sitios son los más fragmentados de las localidades estudiadas; por una parte, en Tlayecac solo quedan remanentes de vegetación natural debido a los cultivos, el crecimiento de la mancha urbana y el pastoreo de ganado vacuno (^{Messenguer, 2013}). Por otro lado, en Tlayca también se ha observado ganado vacuno y el sitio de estudio se encuentra rodeado por cultivos, actividades piscícolas y extracción de roca (^{Flores, 2018b}). Así mismo, el estado de Morelos se ha descrito como uno de los sitios con mayor extracción de corteza para uso medicinal, además de uso para postes y combustible (^{Anzures et al., 2004}; ^{Solares et al., 2012}). Por otra parte, en el estado de Morelos se ha observado una disminución en el patrón de precipitación en los últimos 50 años, con un aumento en la frecuencia de fuertes precipitaciones concentradas en pocos días durante la temporada de lluvias (^{CICESE, 2017}), lo cual tiene efectos en las poblaciones.

Las variaciones en el número de individuos de A. adstringens también se deben a que cada sitio de estudio presenta condiciones ambientales particulares, en donde los procesos ecológicos y factores tanto abióticos como bióticos, tales como el clima, la topografía, la precipitación, la germinación, el establecimiento, la competencia, la fenología, la reproducción, la producción de semillas, la dispersión de éstas, y el disturbio presente, influyen de diferente manera en la densidad, distribución espacial de los individuos y en la estructura de las poblaciones (^{CICESE, 2017}; ^{Ceccon et al., 2006}; ^{Fortin y Dale, 2005}; ^{Hubbell, 1979}).

El patrón de distribución al azar encontrado en Cuambio, Tilapa y Tlayca pueden ser resultado de la homogeneidad ambiental, donde los individuos de la misma especie pueden encontrarse en cualquier lugar (^{Begon et al., 2006}; ^{Fortin y Dale, 2005}), que posiblemente es causada por la fragmentación del hábitat presente en estos sitios, y por la capacidad de A. adstringens de establecerse a cielo abierto o bajo el dosel de plantas adultas. Por otro lado, el patrón de distribución agregado encontrado en Ostula, Pómaro, Papalutla, Coxcatlán, El Cascalote y el cañón el Sabino, es similar a lo descrito para la especie en San Juan Bautista Jayacatlán, Oaxaca (^{Zamora, 2003}). Este patrón de distribución de individuos se debe a la disposición heterogénea de los recursos, en dónde los individuos de la misma especie se encuentran cercanos en áreas con condiciones favorables y suficientes recursos (^{Begon et al., 2006}; ^{Fortin y Dale, 2005}). Sin embargo, este patrón se contrapone con la disposición espacial azarosa, por lo que es importante indagar que factores o procesos intervienen en la supervivencia de A. adstringens a lo largo del ciclo de vida y determinan un patrón u otro, ya que a nivel de paisaje solo se observó diferente grado de fragmentación.

En el BTC, los recursos como los nutrientes y el agua varían estacionalmente en el ambiente, son fundamentales para el establecimiento, supervivencia y desarrollo de los individuos, lo cual tiene efectos en la distribución espacial (^{Balvanera et al., 2011}; ^{Begon et al., 2006}; ^{Mooney et al., 2009}; ^{Rahman et al., 2009}). Así como el alto grado de deforestación presente en Cuambio, Tlayca y Tlayecac (^{Flores-Alta et al., 2019}; ^{Flores, 2018b}), y la presencia de eventos esporádicos de lluvias torrenciales que se han presentado en los últimos 2 sitios (^{CICESE, 2017}), pueden ocasionar la pérdida de recursos por deterioro o arrastre por agua, lo que afecta la productividad del ambiente (^{Bertol et al., 2003}; ^{Borrelli y Oliva, 2001}; ^{Jie et al., 2002}), y por el contrario, en los sitios donde existe muy baja perturbación, los recursos tienen mayor disponibilidad. Así mismo, el tipo de dispersión de semillas que presenta la especie de estudio (anemocoria) puede influir en el tipo de distribución en algunos sitios, además del ambiente y del espacio (^{Balvanera et al., 2011}; ^{Condit et al., 2000}), por lo que las diferencias en los patrones de distribución pueden deberse también a la variedad de interacciones entre organismos y el medio abiótico (^{Fortin y Dale, 2005}; ^{Sapkota et al., 2009}).

Por otra parte, encontramos que el mayor porcentaje de individuos de A. adstringens corresponde a la categoría adulto 1, es decir adultos jóvenes, lo que nos indica que existe poca regeneración reciente de la especie, lo cual es característico de poblaciones que se encuentran en declive (^{Mollers, 2008}; ^{Rahman et al., 2009}). Este patrón también fue observado para A. adstringens en Chamela, Jalisco (^{Fernández, 2014}), y en el Ejido el Limón, Morelos (^{Beltrán, 2018}), así como para otras especies caducifolias como Bursera glabrifolia en Jayacatlán, Oaxaca (^{Hernández-Apolinar et al., 2006}).

En las poblaciones de A. adstringens del presente estudio, el bajo número de individuos pre-reproductivos juega un papel fundamental en la estructura, la composición y en la diversidad de los bosques (^{Johnson et al., 2017}). La ausencia de plántulas puede ser resultado de la temporada en la que se realizaron los muestreos, ya que la mayoría de éstos se realizaron en temporada de secas, y la germinación de semillas del BTC ocurre durante la temporada de lluvias (^{Ceccon et al., 2006}). Así mismo, puede deberse a un rápido desarrollo de las plántulas (plantas que ya han perdido los cotiledones), que no permitió observarlas durante los muestreos y lo que se registró fueron juveniles. Además, se ha descrito que especies del BTC con semillas de tamaño pequeño, como es el caso de A. adstringens, tienen menor sobrevivencia que las especies que presentan semillas de tamaño grande (^{Khurana y Singh, 2004}; ^{Maass y Burgos, 2011}).

Por otro lado, los cambios en el patrón de precipitación de los últimos años y la disminución de ésta en Coxcatlán, Tlayca y Aquila (^{CICESE, 2017}; ^{Serrano, 2018}; ^{Villalobos, 2021}), pueden tener consecuencias negativas sobre las plántulas y juveniles, ya que se ha descrito que los eventos de lluvia torrencial, pueden afectar las semillas y plántulas de especies filogenéticamente cercanas (Beiselia y Bursera, Burseraceae), debido a que éstas se saturan de agua y mueren “ahogadas” o son arrastradas por la corriente debido al incipiente desarrollo del sistema radicular (^{Villalobos, 2021}), así mismo, estas especies requieren un periodo de lluvia constante durante al menos 4 meses para que se diferencie eficientemente el sistema vascular y puedan resistir el periodo de sequía, un proceso que puede estar limitado por los cambios en el régimen de lluvias o las sequías prolongadas (^{Azcón-Bieto y Talón, 2008}; ^{Ortega y Collazo, 2013}; ^{Serrano, 2018}).

El bajo número de individuos de la categoría juvenil 1 se puede explicar por una limitación en la producción de semillas, la germinación, la sobrevivencia y establecimiento de plántulas, ya que los individuos de primeras etapas son más susceptibles a la muerte que los individuos adultos (^{Ceccon et al., 2006}; ^{Fortin y Dale, 2005}; ^{Hubbell, 1979}; ^{Lieberman y Li, 1992}). La producción de semillas es afectada negativamente por una reducción de la fecundidad, limitación por polen o errores en el desarrollo de la semilla, entre otros factores (^{Callaway, 2007}; ^{Jordano et al., 2008}). Mientras que la sobrevivencia de plántulas y juveniles depende de factores bióticos y abióticos como las especies de plantas bajo las cuales crecen, las características de los micrositios en donde se lleva a cabo la germinación y establecimiento de plántulas, además de la maduración de tejidos del organismo, la depredación, y de las actividades antrópicas (^{Callaway, 2007}; ^{Beltrán, 2018}).

Particularmente en el municipio de Aquila (al que pertenecen Ostula y Pómaro), se sabe que A. adstringens es una especie agroforestal, la hojarasca producida por ésta y otras especies arbóreas es usada como forraje para el ganado, lo que contribuye a su conservación (^{Herrera et al., 2018}). Sin embargo, el pisoteo del ganado durante el forrajeo puede ocasionar la muerte de plántulas y/o juveniles o incluso puede haber ingesta de pequeños individuos, ya que es poco probable que hayan desarrollado metabolitos secundarios y sean palatables para el ganado (^{Serrano, 2018}). Otro factor observado en Aquila es la quema de hierbas al inicio de la época seca para permitir que la hojarasca esté “limpia” para el ganado; aunque se cuida que los árboles adultos permanezcan sin daños, es probable que se influya negativamente en el establecimiento de nuevos individuos.

Durante las primeras etapas de vida, las plantas asignan la mayor cantidad de recursos al crecimiento y expansión de tejidos, quedando en segundo término la construcción de defensas químicas y físicas, lo que las hace más susceptibles a herbívoros y patógenos (^{García y Collazo, 2013}). Particularmente, en las especies de Bursera de Coxcatlán y Tilapa (en donde habita A. adstringens) se ha observado que las plántulas no desarrollan defensas químicas o físicas, es alrededor de los 4 meses de edad, cuando ya son juveniles, que comienza la producción de defensas contra la herbivoría (^{Serrano, 2018}), es factible que A. adstringens tenga un desarrollo similar al crecer bajo las mismas condiciones ambientales y estar relacionado filogenéticamente, sin embargo, es un tema que debe estudiarse.

Adicionalmente, en A. adstringens se han reportado escarabajos xilófagos (Cerambycidae), considerados parásitos para varias especies de plantas de importancia comercial, ya que las larvas perforan los tejidos conductores provocando la caída de las ramas y facilitan la entrada de otra fauna nociva (^{Calderón-Cortés et al., 2016}; ^{López-Martínez et al., 2015}; ^{Vargas-Cardoso et al., 2018}). Si bien esta afectación en los árboles adultos puede no ser importante, en el caso de los individuos pre-reproductivos que comienzan a producir ramas, podría significar una reducción en el crecimiento del individuo o incluso la muerte, ya que la cantidad de meristemos para el crecimiento lateral y/o apical aún es reducido y se afecta directamente el sistema vascular.

En el presente estudio no encontramos poblaciones en crecimiento o estables debido a que en todas las localidades la estructura poblacional es decreciente, ya que la proporción de adultos es mayor que la de individuos pre-reproductivos, como es el caso de la población de A. adstringens en el Ejido el Limón, Morelos (^{Beltrán, 2018}). La categoría juvenil 2 es la menos representada en 8 sitios, contrario a lo descrito para A. adstringens en La Barranca de Horcones y la laguna en San Juan Bautista Jayacatlán, Oaxaca, en donde la categoría mejor representada es la que presenta una altura mayor de 100 cm, pero menor de 10 cm de diámetro a la altura del pecho (^{Zamora, 2003}). Por otro lado, se ha observado un patrón similar al encontrado para A. adstringens en especies filogenéticamente relacionadas como; Bursera fagaroides, B. arida, B. submoniliformis y B. schlechtendalii en Coxcatlán, Puebla, así como para B. submoniliformis, en Guadalupe Victoria, Puebla (^{Gonzalez, 2019}), para Bursera cuneata en la Ciudad de México (^{Hernández, 2018}), y para el género Bursera en Pómaro (^{Hernández y Ramos-Ordoñez, 2018}).

La baja representatividad de individuos juveniles (1 y 2) en el presente estudio indica que en algún momento la especie presentó un establecimiento irregular, y que su producción de semillas y el establecimiento no es constante en el tiempo, debido a las condiciones a las que se encontraron expuestos, por lo que posiblemente A. adstringens es una especie sensible a los cambios que afectan su establecimiento (^{Zamora, 2003}). Así mismo, esta disminución en el número de individuos también se puede atribuir al pastoreo que existe en los sitios estudiados, principalmente en Tilapa, Coxcatlán, Tlayca, Tlayecac y Papalutla. Tal disminución en el número de individuos podría traer consecuencias negativas a las poblaciones generando cuellos de botella (^{Groenendijk et al., 2012}; ^{Tolera et al., 2013}), ya que los individuos de estas categorías en los siguientes años remplazaran a las subsecuentes, resultando en una reducción en el número de individuos adultos.

Uno de los factores que determina la sobrevivencia de individuos juveniles es el microclima, debido a que éste afecta su crecimiento. Se ha descrito que especies de plantas con las que se asocian proporcionan mayor humedad, menor compactación del suelo, menor radiación solar y bajo éstas existe mayor cantidad de materia orgánica (^{Callaway, 2007}; ^{Verdú y García-Fayos, 2002}). Esta asociación también se ha descrito en otras especies en el BTC, la cual permite el establecimiento de diferentes plantas en este ambiente altamente estacional (^{Arriaga et al., 1993}; ^{Contreras, 2015}; ^{Méndez et al., 2006}; ^{Ruiz, 2007}).

No obstante, en el presente estudio se encontró que los individuos de A. adstringens no necesariamente requieren de una planta perenne para establecerse, como es el caso de especies pertenecientes al género Bursera (^{Gonzalez, 2019}). Sin embargo, los individuos juveniles de A. adstringens que se encuentran asociados se localizan principalmente bajo el dosel de árboles de diferentes géneros, debido a que estos son los elementos dominantes en este tipo de vegetación (^{Flores, 2018b}; ^{Rivera-Ortíz et al., 2013}; ^{Trejo, 1998}; ^{Vázquez, 2007}).

Encontramos que un gran número de juveniles de A. adstringens asociados se encontraron bajo el dosel de la misma especie en El Cascalote, Coxcatlán, Cuambio y el cañón el Sabino, a pesar de que se ha descrito que existe mayor mortalidad de plántulas bajo individuos conespecíficos, debido a patógenos y herbivoría (^{Howe y Smallwood, 1982}). Además, se sabe que existe mayor competencia entre especies filogenéticamente cercanas debido a que tienen nichos similares (^{Valiente-Banuet y Verdú, 2008}). Dicha asociación puede deberse a la dispersión de semillas anemócora propia de la especie (dispersión por viento), que influye en la distribución limitada de las semillas (^{Balvanera et al., 2011}), las cuales pueden caer cerca de la misma especie, incluso del progenitor. Por lo que sería importante evaluar en estudios futuros el sexo de los individuos adultos de A. adstringens con los cuales se asocian los juveniles, lo cual podría ayudar a explicar un poco más el proceso de facilitación con individuos de la misma especie.

Esta asociación se puede explicar ya que las especies vecinas a los individuos asociados de A. adstringens, son taxa filogenéticamente alejados y neutralizan los efectos negativos como la competencia, dado que puede existir mayor complementariedad, y el efecto de los depredadores es menor; además, la presencia de una red de micorrizas que une a las especies vegetales permite la transferencia de nutrientes y protección contra patógenos (^{Cerón, 2015}; ^{Montesinos-Navarro et al., 2019}; ^{Valiente-Banuet y Verdú, 2013}).

Debido a que A. adstringens puede establecerse en espacios desprovistos de vegetación y otras plantas se establecen bajo su copa en sitios como el Clarín (Guanajuato), Talyca (Morelos), la cueva del Maíz (Coxcatlán, Puebla), y en el cañón el Sabino en Oaxaca (Godinez-Manriquez y Contreras-González, 2019; ^{Gonzalez, 2019}; ^{Segundo-Hernández y Ramos-Ordoñez, 2016}; ^{Serrano-Rosas et al., 2018}), se puede considerar a A. adstringens como especie facilitadora de otras especies, como es el caso de Mimosa luisana en el BTC y en el matorral xerófilo (^{Contreras, 2015}; ^{Valiente-Banuet y Ezcurra, 1991}; ^{Valiente-Banuet y Verdú, 2007}).

Sin embargo, en algunos puntos de su distribución, A. adstringens se extrae en gran medida para uso medicinal (^{Solares et al., 2012}), presenta diferencias fisiológicas entre sexos (que podría afectar más a un sexo en ambientes estresantes [^{Beltrán, 2018}]), produce un bajo número de semillas por fruto (^{Guzmán-Pozos y Cruz-Cruz, 2014}), los juveniles tienen una tasa relativa de crecimiento baja (^{Rincón y Huante, 1993}), así como una limitada transición entre etapas de desarrollo. Todos estos son factores que pueden llegar a poner en riesgo la dinámica de las comunidades en donde se encuentra presente la especie, ya que es facilitadora de otras especies de plantas del BTC.

En este trabajo también resultaron asociaciones importantes entre los juveniles de A. adstringens y los adultos de M. luisana y Parkinsonia praecox en El Cascalote y el cañón el Sabino, C. parvifolia y M. luisana en Tilapa. Parkinsonia praecox, también ha sido descrita como nodriza de otras especies de plantas ya que puede favorecer su crecimiento (^{Contreras, 2015}; ^{Valiente-Banuet y Escurra, 1991}; ^{Valiente-Banuet y Verdú, 2007}). Se sabe que bajo el dosel de M. luisana y P. praecox, se presentan suelos con mayor cantidad de nutrientes que los sitios desprovistos de vegetación, y hay mayor sobrevivencia de plantas asociadas a éstas (^{García-Chávez et al., 2014}; ^{Valiente-Banuet et al., 1991}), tales interacciones presentes entre plantas pueden influir positivamente en la estructura de las comunidades (^{Valiente-Banuet y Verdú, 2013}; ^{Valiente-Banuet et al., 2014}).

Los individuos juveniles de A. adstringens asociados al dosel de una planta adulta principalmente se encontraron orientados al norte, este y oeste, lo que nos indica que evitan establecerse al sur de la planta con la cual se asocian, ya que al norte la radiación solar anual es menor, mientras que al oeste reciben mayor radiación en la tarde, sin embargo, es una zona menos cálida (^{Valiente-Banuet et al., 1991}). Además, en especies suculentas se ha descrito mayor sobrevivencia de juveniles de reciente establecimiento al este, y existe una mayor distribución de suculentas al norte con respecto a las plantas perennes con las que se asocian en regiones áridas de México (^{Valiente-Banuet y Ezcurra, 1991}). No obstante, en la Reserva de Chamela-Cuixmala se encontró que el crecimiento de los individuos juveniles de A. adstringens no es diferente en una orientación sur y norte (^{Fernández, 2014}).

Los individuos de A. adstringens que mostraron marcas de extracción de corteza fueron pocos, y éstos se encontraron en Guerrero, Puebla y en el cañón el Sabino en Oaxaca. En general, encontramos que los individuos que presentaron marcas de extracción principalmente fueron los individuos de mayores tamaños, lo que nos indica que al parecer existe una extracción selectiva de la corteza. El sitio que presentó mayor porcentaje de individuos con marcas de extracción de corteza fue Cuambio, y es el sitio que presenta menor número de individuos adultos con respecto a los otros sitios estudiados.

Durante la toma de datos en Cuambio se observó que los sitios con individuos de A. adstringens son parches de vegetación rodeados de parcelas dedicadas a la agricultura y ganadería, y Guerrero es uno de los estados que presenta mayor velocidad de degradación de la vegetación natural primaria (2.46 % anual; ^{Semarnat, 2016}). Por el contrario, en los sitios de Morelos no se encontraron árboles con marcas de extracción, pero es el sitio con menor cantidad de individuos, Tlayecac se encuentra rodeado de parcelas dedicadas a la agricultura y a la ganadería y Tlayca se encuentra sometido al pastoreo y al turismo; además, Morelos también es uno de los estados con mayor velocidad de degradación de vegetación primaria (1.26 % anual; ^{Semarnat, 2016}).

A pesar de que se ha descrito una elevada demanda de la corteza de esta especie para venta (^{Solares et al., 2012}), la extracción de corteza en Cuambio, en el cañón el Sabino y en Coxcatlán se ha indicado de manera informal por los pobladores de estos sitios, que la especie solamente es utilizada por personas de las comunidades cercanas, es decir, es aprovechada localmente y la extracción no es con fines comerciales, al menos por parte de los pobladores. Esta información coincide con un trabajo previo en la comunidad de San Rafael, dentro del municipio de Coxcatlán (^{Dávila y Lira-Saade, 2002}).

Finalmente, los resultados de este trabajo complementan estudios anteriores desde una perspectiva poblacional y confirman que A. adstringens es una especie en alto riesgo, tanto por sus características ecológicas y socioeconómicas (^{Beltrán-Rodríguez et al., 2017}, ²⁰¹⁸), como por las características antrópicas de las poblaciones naturales. El tamaño poblacional, particularmente en sitios con poblaciones pequeñas, se podría incrementar con ayuda de planes de manejo y reforestación. Las estructuras poblacionales en declive coincidieron en una baja densidad de individuos pre-reproductivos como se esperaba, por lo que es importante analizar los factores y procesos que causan esto en diferentes poblaciones a lo largo de la distribución de la especie, e intentar cambiar dicho patrón hasta establecer poblaciones en crecimiento.

Los resultados obtenidos sobre el sitio de estable- cimiento de A. astringens contribuyen al conocimiento de la especie, pero también pueden dar la pauta para la investigación sobre el crecimiento de los individuos. El patrón de distribución espacial agregado y azaroso y la falta de diferencias entre los individuos establecidos a cielo abierto y los que se establecen bajo otras plantas, sugieren que es posible realizar plantaciones en sitios perturbados aprovechando tanto la vegetación remanente como los sitios abiertos. Por otro lado, la evidencia de extracción de corteza en árboles de talla grande podría indicar la aplicación del conocimiento local por parte de los pobladores en los sitios de estudio, sin embargo, es necesario un estudio etnobotánico para evaluar la sostenibilidad de esta práctica. No obstante, en algunos puntos de su distribución, la extracción elevada con fines comerciales (Anzures et al., 2006; ^{Beltrán-Rodríguez et al., 2017}; ^{Casas, Lira et al., 2016}; ^{Solares et al., 2012}) puede poner en riesgo a la especie, como ha ocurrido con algunas otras del género Bursera que han sido sobre explotadas para la obtención de combustibles, aceites esenciales, medicinas locales, y en la elaboración de artesanías, con una evidente disminución en sus poblaciones (^{dos Santos, 2017}; ^{Hersch, 2009}; ^{Weeks y Tye, 2009}).