Introducción

La variación morfológica intraespecífica puede ser el resultado del desarrollo ontogénico (^{Frédérich, Adriaens, & Vandewalle, 2008}; ^{Díaz & Campos, 2014}), la variación genética (^{Frédérich, Liu, & Dai, 2012}) o la influencia de la heterogeneidad ambiental (^{Cardini, Jansson, & Elton, 2007}; ^{Cardini & Elton, 2009}). En el último caso, los organismos están expuestos a diversas condiciones ambientales que influyen en sus procesos biológicos, lo cual, se puede reflejar en el comportamiento y las adaptaciones morfológicas de las especies (C^{ardini, Filho, Polly, & Elton, 2010}).

En el ambiente marino se presentan múltiples variables ambientales que podrían explicar la variación morfológica de las especies marinas a lo largo de gradientes ambientales y latitudinales (^{Márquez & Van Der Molen, 2010}; ^{Nava-Ferrer & Severeyn-Valbuena, 2011}; ^{Sánchez, Sepúlveda, Brante, & Cárdenas, 2011}). Particularmente, en los arrecifes de coral del océano Atlántico mexicano, se han reportado diferencias morfológicas intraespecíficas del pez Abudefduf saxatilis (^{Piñeros, Ríos-Cárdenas, Gutiérrez-Rodríguez, & Mendoza-Cuenca, 2015}) en respuesta a que los individuos se distribuyen en parches geográficamente separados, a lo largo de gradientes latitudinales, soportando diferentes condiciones ambientales (^{Arias-González, Legendre, & Rodríguez-Zaragoza, 2008}).

El conjunto de arrecifes que se distribuyen a lo largo del litoral veracruzano denominado como el Corredor Arrecifal del Suroeste del Golfo de México (CASGM) por ^{Ortiz-Lozano et al. (2013}) tiene variables ambientales que cambian a lo largo del gradiente latitudinal, como es el caso de la salinidad, que disminuye en dirección norte-sur o la temperatura superficial del mar y la concentración de clorofila-a que aumentan en la misma dirección (^{Zavala-Hidalgo, Morey, & O’Brien, 2003}; Salas-Pérez, Ocaña-Valencia, & ^{González-Gándara, 2015}). Por su lado, los patrones de corrientes marinas dentro del CASGM también cambian en el tiempo y el espacio (^{Salas-Monreal et al., 2017}).

Se ha sugerido que dichas variaciones ambientales influyen en la composición, riqueza y distribución de algunas especies arrecifales. Tal es el caso de los equinodermos, que muestran variabilidad espacial de composición y riqueza entre el norte y sur de Veracruz, probablemente como resultado de la presión ambiental que ejerce la condición oligotrófica de las aguas en el norte y las eutróficas del sur (^{González-Gándara et al., 2015}). Asimismo, se ha reportado que el ensamblaje coralino está estructurado en función de los gradientes de turbidez y productividad oceánica que hay a lo largo del CASGM (^{Jordán-Garza, González-Gándara, Salas-Pérez, & Morales-Barragán, 2017}). De este modo, el ambiente heterogéneo del corredor puede influir en la biología de los organismos que ahí se distribuyen. A pesar de ello, aún no se han desarrollado estudios para determinar el efecto geográfico y ambiental en la biología de las especies, específicamente en su morfología.

Se ha documentado que las características de algunas estructuras en erizos de mar tales como el grosor de la testa o caparazón, así como número y tamaño de las espinas se modifican por la intensidad de oleaje, la temperatura y por la presencia o ausencia de depredadores (^{Ebert, 1968}; ^{Dix, 1970}; ^{Marcus, 1980}, 1983; ^{Byrne & Przeslawski, 2013}; ^{Byrne et al., 2014}). Las estructuras funcionales de los erizos de mar como la linterna de Aristóteles (aparato masticador) también pueden presentar variación morfológica como respuesta al cambio en la disponibilidad del alimento ya que la fuerza de forrajeo aumenta cuando el recurso es limitado o cuando sus principales presas tienen estructuras calcáreas duras (Ebert, 1980; ^{Hagen, 2008}; ^{Gianguzza & Bonaviri, 2013}; ^{Epherra et al., 2015}). De igual manera, la fuerza de forrajeo cambia dependiendo del tipo de sustrato en el que se encuentren. Si el sustrato es duro o plano, la linterna crecerá más respecto al diámetro de la testa, lo contrario ocurre si se desarrolla en sustratos blandos o con presencia de cavidades (^{Hernández & Russell, 2009}). También, las variaciones de salinidad generan cambios en el tamaño y peso de la linterna (^{Lau, Lau, Qian, & Qiu, 2009}), motivo por el cual puede ser una estructura que proporcione información sobre la plasticidad morfológica y al mismo tiempo sea un indicador de heterogeneidad ambiental.

Por lo anterior, los erizos de mar pueden ser modelos adecuados para estudiar la variación del tamaño y la forma como respuesta a las distintas presiones ambientales y gradientes latitudinales que se presentan en el litoral Veracruzano. Por ello, el presente estudio analiza la variación morfológica de la linterna de Aristóteles de Eucidaris tribuloides a lo largo del CASGM.

Materiales y métodos

Área de estudio: El Corredor Arrecifal del Suroeste del Golfo de México (CASGM) cubre una distancia de 680 km lineales dentro de la plataforma continental del oeste del Golfo de México (^{Geissert-Kientz, 1999}). Está conformado por sistemas de arrecifes coralinos de diferentes dimensiones que se distribuyen de norte a sur frente a la costa del estado de Veracruz, desde Tamiahua hasta Coatzacoalcos (^{Ortiz-Lozano et al., 2013}). No obstante, se pueden distinguir tres grandes sistemas: el Sistema Arrecifal Lobos-Tuxpan (DOF, 2009), el Parque Nacional Sistema Arrecifal Veracruzano (DOF, 1992) y los Arrecifes de los Tuxtlas.

En el norte, el Sistema Arrecifal Lobos-Tuxpan (SALT) se caracteriza por presentar una salinidad promedio de 36.5 ups (^{Zavala-Hidalgo et al., 2003}), una temperatura superficial promedio de 23.5 °C, variando de 20 a 30 °C y una concentración de clorofila-a de 0.76 mg/m³, con mínimos de 0.2 mg/m³ y máximos de 1.6 mg/m³. Asimismo, está influenciado principalmente por el Río Tuxpan, cuyo promedio de descarga es de 138 m³/s (^{Salas-Pérez et al., 2015}). Los sitios de recolección en el norte fueron los arrecifes; Blanquilla (21°32'24.80" N & 97°16'44.67" W), Medio (21°30'32.75" N & 97°15'02.64" W), Tuxpan (21°1'39.77" N & 97°11'38.45" W) y Blake (20°45'30" N & 96°59'35" W).

En el centro, el Sistema Arrecifal Veracruzano (SAV) está definido por una salinidad promedio de 36 ups con un rango entre 34 y 36.5 ups, una temperatura superficial del mar de 24 °C en promedio, con un mínimo de 22 °C y un máximo de 31 °C, mientras que la concentración de clorofila-a promedio es de 2.5 mg/m³ y oscila entre 0.9 mg/m³ y 4 mg/m³. Este sistema está influenciado principalmente por el río Jamapa, del cual se ha estimado una descarga promedio de 45 m³/s con máximos de 380 m³/s y mínimos de 8 m³/s (^{Perales-Valdivia, Sanay-González, & Marín-Hernández, 2015}; ^{Salas-Pérez et al., 2015}) y también hay influencias de los ríos La Antigua y Papaloapan (^{Krutak, 1997}). Los sitios de recolección fueron los arrecifes; La Segunda Laja (19°28'44" N & 96°18'24" W) en las inmediaciones de la localidad de Chachalacas y en la Anegada de adentro (19°13'54.8" N & 96°3'27.5" W).

En el sur, los Arrecifes de los Tuxtlas (AT) están caracterizados por presentar salinidad promedio de 35.5 ups (^{Zavala-Hidalgo et al., 2003}), temperaturas superficiales promedio de 25.64 °C con mínimos de 22 y 30.5 °C, una concentración de clorofila-a promedio de 4.9 mg/m³ y un rango de 0.82 y 10.91 mg/m³. Este sistema está influenciado principalmente por el río Coatzacoalcos con una descarga promedio de 353 m³/s (^{Salas-Pérez et al., 2015}). Los sitios de recolección fueron los arrecifes; Zapotitlán (18°30'56.9" N & 94°47'12.8" W), Balzapote 18°37'20" N & 95°3'58.3" W), La Perla (18°32'33.8" N & 94°49'4.1" W), Tripié (18°10’33.6" N & 94°22'03" W) y Palo Seco (18°10'33.1" N & 94°31'32" W).

Las corrientes superficiales varían en el tiempo y el espacio dentro del corredor. La dirección de las corrientes costeras de Veracruz es de norte-sur de septiembre a marzo, mientras que de abril a agosto la dirección es de sur-norte (^{Zavala-Hidalgo et al., 2003}). Por otro lado, la información sobre las corrientes superficiales del mar Caribe y Golfo de México proporcionadas por la Administración Nacional de la Atmósfera y el Océano, 2015 (NOAA) a través del Nodo Regional del Caribe y Golfo de México Coast-Watch (http://cwcaribbean.aoml.noaa.gov/data.html) evidencian giros ciclónicos y anticiclónicos temporales frente a las costas de Veracruz, lo cual, puede influir en la dispersión larval de las especies a lo largo del CASGM (^{Salas-Monreal et al., 2018}).

Manejo de ejemplares: Se tomaron 26 ejemplares de Eucidaris tribuloides de la Colección de Invertebrados de la Universidad Veracruzana (número de catálogo: DF-CC-279-13; muestras tomadas de marzo a octubre 2013). Además, se recolectaron 77 individuos en campo para completar el tamaño de la muestra previamente establecida (marzo a junio 2016). Considerando la máxima talla de madurez sexual registrada por ^{McPherson (1968}), se procuró capturar ejemplares mayores de 25 mm de diámetro de testa para evitar variación relacionada con la madurez sexual. En total, se analizaron 104 individuos de los cuales 30 corresponden al SALT (17 en campo y 13 de colección), 35 al SAV y 39 (26 en campo y 13 de colección) al SAT.

Una vez capturados los ejemplares, se narcotizaron dejándolos reposar en un recipiente oscuro con suficiente agua marina y una cucharada de cloruro de magnesio durante 1 h. Luego se fijaron en alcohol etílico al 96 % durante 24 h, se etiquetaron con los siguientes datos: recolector, localidad, profundidad, fecha de colecta y posición geográfica del sitio de recolecta y finalmente se conservaron en seco. Las espinas de cada erizo de mar se removieron con una inmersión de hipoclorito de sodio durante 10 a 20 min, esto para obtener las medidas lineales más precisas del diámetro de la testa. Se extrajo la linterna de Aristóteles de cada individuo y se registró la altura de la estructura, todas las medidas se obtuvieron con un vernier digital con precisión de 0.01 mm. Posteriormente, se removió el tejido adherido a la linterna de Aristóteles con una inmersión en hipoclorito de sodio al 10 % durante 24 h para desarticular los osículos de la linterna.

Análisis estadístico de morfometría clásica: Con el diámetro de la testa y la altura de la linterna se realizó una regresión para cada grupo (norte, centro y sur) con el fin de calcular la alometría entre el diámetro de la testa versus la altura de la linterna. Un análisis de covarianza (ANCOVA) fue utilizado para comparar las pendientes de esas regresiones a un nivel de significancia α < 0.05, donde la altura de la linterna de Aristóteles fue utilizada como variable dependiente, las zonas; norte, centro y sur como variables agrupadas y el diámetro de la testa como covariable. Lo anterior se realizó probando previamente la homogeneidad de las pendientes y con el análisis de Levene se probó homogeneidad de varianzas de los datos mediante el Software estadístico SPSS.

Morfometría geométrica: Se seleccionó una rótula por cada linterna de Aristóteles para obtener imágenes fotográficas de su parte aboral (Fig. 1A). Las imágenes se tomaron en el Laboratorio de Microscopía y Fotografía de la Biodiversidad en el Instituto de Biología de la UNAM utilizando una cámara LEICADF-C490 de 8 megapíxeles que viene integrada a un microscopio estereoscópico multifocal LEICA-Z16APOA, manipulado desde un ordenador mediante el software “Leica Application Suite” (LAS por sus siglas en inglés). Asimismo, se obtuvieron imágenes fotográficas de una semipirámide por cada linterna orientado en su plano lateral interno con una cámara de 16 megapíxeles de la marca Panasonic DMH-ZS10 (Fig. 1C).

Posteriormente, se colocaron landmarks y semilandmarks con el software TpsDig 2.26 (^{Rohlf, 2016}). Estos son puntos homólogos que se colocan sobre las imágenes de las estructuras a estudiar considerando los criterios de ^{Bookstein (1991}). En este estudio, se asumió simetría bilateral para el caso de la rótula. Por lo tanto, se analizó el lado derecho de cada estructura para evitar redundancia de información. La configuración consistió de seis landmarks y seis semilandmarks para caracterizar la forma del lado derecho de la rótula (Fig. 1B). El criterio que se utilizó para colocar los semilandmarks consistió en poner seis líneas equidistantes entre sí y transversales al contorno de la rótula con el software Makefan8 (IMP) (^{Sheets, 2014}a), posterior a ello, se colocaron los semilandmarks en el punto de intercepción entre las líneas y el contorno (Fig. 1B). Para el caso de la semipirámide se establecieron seis landmarks figura (Fig. 1D).

Fig. 1.  Fig. 1 (A) Linterna de Aristóteles vista en el plano superior, (B) Rótula parte aboral y sus landmarks y semilandmarks, (C) Linterna de Aristóteles vista desde el plano vertical, (D) Semipirámide vista lateral interno y sus landmarks. Imagen modificada de ^{Märkel (1979}). (A) Aristotle's lantern view from above, (B) Rotula view from above and its landmarks and semilandmarks, (C) Aristotle's lantern view from the side, (D) Left demipyramid inner side and its landmarks. From ^{Märkel (1979}). 

Se realizó un ajuste de Procrustes a la configuración de coordenadas de cada estructura analizada, método de superposición que consiste en eliminar todos los efectos por tamaño, posición y rotación de los objetos, descartando así la variación de la forma no biológica. Para el caso particular de la rótula, los semilandmarks se ajustaron con el método de “minimum bending energy” usando el Software TpsRelw 1.65 (Rholf, 2016). También, se extrajo el valor del tamaño centroide (información completa del tamaño) de cada estructura y se aplicó la función de logaritmo natural para normalizar los datos. Con ello, se realizaron ANOVAs para comparar el tamaño centroide de las rótulas y semipirámides entre los grupos norte, centro y sur del CASGM. Consecutivamente, se aplicó una prueba de Tukey para comparar las medias de las muestras en caso de presentar diferencias significativas (α < 0.05).

Se efectuó una regresión multivariada de las coordenadas Procrustes (variable dependiente) sobre el tamaño centroide de la rótula y semipirámide (variable independiente), con la finalidad de observar si el cambio de la forma de la rótula y semipirámide se debía al tamaño. También se realizó un Análisis de Componentes Principales Entre Grupos (BGPCA por sus siglas en inglés) utilizando la configuración promedio de las coordenadas Procrustes de cada grupo (norte, centro y sur), esto con base en las sugerencias de ^{Mitteroecker y Bookstein (2011}). Además, se utilizaron gradillas de deformación para observar los patrones e intensidad de la variación morfológica a lo largo de los ejes del BGPCA en el software PcaGeen8 (IMP) (^{Sheets, 2014}b), Dichas gradillas permiten visualizar con mejor precisión la ubicación anatómica con mayor variación morfológica de las estructuras analizadas. Finalmente, se extrajo los “partial warp scores” o variables de la forma con ayuda del Software CoordGen8 (^{Sheets, 2014}c) de la serie IMP8 para realizar un Análisis Multivariado de la Varianza (MANOVA por sus siglas en ingles) para determinar si existen diferencias morfológicas entre los individuos del norte, centro y sur del CASGM.

Resultados

Morfometría clásica: El diámetro de la testa de E. tribuloides osciló entre 20 y 45 mm, con un promedio de 30.44 mm y una desviación estándar de 4.89 mm, mientras que la altura de la linterna de Aristóteles osciló entre 9 y 18 mm con un promedio de 14.5 mm y una desviación estándar de 2.42 mm. La relación entre esas dos estructuras fue lineal en todos los sitios. Los individuos del norte presentaron las linternas más grandes, mientras que los del sur las más pequeñas (Fig. 2). La prueba de Levene no mostró diferencias significativas entre grupos (F = 2.18, P = 0.119) y la interacción entre las zonas y la altura de la linterna indicó que las pendientes fueron homogéneas (F=0.964, P = 0.386). Por su parte, el ANCOVA manifestó diferencias significativas entre grupos (F= 10.74, P < 0.001) con un coeficiente de regresión de 86 %, mientras que la comparación por pares con ajuste de Bonferroni señaló que el norte es diferente al centro (SE = 0.249, P = 0.002) y el sur (SE = 0.253, P < 0.001). Sin embargo, no hubo diferencias entre centro y sur (SE = 0.226, P = 0.771).

Fig. 2.  Fig. 2 Relación del diámetro de la testa versus altura de la linterna de Aristóteles de Eucidaris tribuloides para cada sitio estudiado (norte, centro y sur de Corredor Arrecifal del Suroeste del Golfo de México).  Allometric relationship between test diameter and Aristotle’s lantern of Eucidaris tribuloides length to each site (north, center and south from Corridor of the Southwestern Gulf of Mexico). 

Morfometría geométrica: El análisis de varianza (ANOVA) confirmó la existencia de diferencias significativas entre los grupos norte, centro y sur para el logaritmo natural del tamaño centroide de la rótula (F = 8.02; P = 0.001) y semipirámide (F = 5.96; P = 0.004). La prueba post hoc de Tukey indica que el tamaño de la rótula de los erizos del norte es significativamente mayor que los del centro y sur. Por otro lado, el tamaño de la semipirámide de los organismos del norte es diferente a la de los individuos del sur. Sin embargo, no hubo diferencias significativas entre individuos del centro y norte, tampoco entre los individuos del centro y el sur.

El análisis de regresión múltiple explicó el 4.6 % de la variación de la forma para la rótula y el 11.26 % para semipirámide, lo que sugiere una posible ausencia de correlación entre el tamaño y la forma de la linterna de Aristóteles para Eucidaris tribuloides. Los MANOVAs confirmaron la existencia de diferencias significativas entre los grupos norte, centro y sur para la rótula (Lambda de Wilks = 0.384, P = 0.001) y semipirámide (Lambda de Wilks = 0.6432, P < 0.0001). La comparación por pares mediante la T² de Hotelling mostró que el norte fue significativamente diferente al sur para la rótula y semipirámide (Tabla 1).

Por otro lado, el Análisis de Componentes Principales Entre Grupos para la rótula mostró la distribución de los individuos en dos ejes, el primero contiene 75.36 % de la variación de la forma. La proyección de los individuos sobre los dos ejes canónicos reveló solapamiento entre los sistemas arrecifales. Pese a ello, se observa que los individuos del norte presentan mayor variación morfológica en comparación con los individuos del sur (Fig. 3A). Para el caso de la semipirámide, la distribución de los individuos fue sobre dos ejes canónicos de los cuales el primer eje contuvo el 87.38 % de la variación. La proyección de los individuos sobre los dos ejes canónicos presenta gran solapamiento entre los grupos del norte, centro y sur. Sin embargo, se observa que los individuos del norte presentaron mayor variación que los individuos del sur (Fig. 3B).

Fig. 3.  Fig. 3 Análisis de Componentes Principales entre grupos de la rótula (A) y Semipirámide (B) de Eucidaris tribuloides. Between-Groups Principal Component Analyses of rotula (A) and demi-pyramid (B) of Eucidaris tribuloides. 

Las rejillas de deformación revelaron algunas diferencias en la forma de la rótula, principalmente en el cóndilo, que es más robusto en los individuos del centro. También, la fosa lateral superior de los individuos del norte es angosta (Fig. 4). Para el caso de la semipirámide, se revelaron algunas diferencias morfológicas donde los individuos que se distribuyeron sobre los valores negativos de la primera función discriminante de la figura 3B son más robustas. Por el contrario, los individuos que se distribuyeron sobre los valores positivos de la primera función discriminante de la figura 3B poseen semipirámides más estrechas. Por otro lado, existe una notoria variación de la forma en la parte superior de la semipirámide (Fig. 4).

Fig. 4.  Fig. 4 Variación de la forma de la rótula y semipirámide a lo largo de la primera función discriminante de la figura 3, cada rejilla mostró la predicción de la forma para A) -0.1 y B) 0.1, las figuras de color negro muestran la forma promedio y las figuras de color gris muestran la variación de la forma.  Variation of the shape from rotula and semipyramid along the first discriminant function of the figure 3, B. Each grid showed the prediction of shape to A) -0.1 and B) 0.1, the black figures show the average shape and the gray figures show the variation of the shape. 

Discusión

La morfometría geométrica es una herramienta que sirve para detectar cambios sutiles de la forma (^{Zelditch, Swiderski, & Sheets, 2012}). En biología, estos cambios pueden ser atribuidos a los procesos de especiación (^{Binks, Evans, Kennington, & Price, 2011}), a la baja o nula interacción genética entre poblaciones (^{Frédérich et al., 2012}) o a la plasticidad morfológica (^{Vega-Sánchez, Camacho-Morales, Chassin-Nora, & Mendoza-Cuenca, 2012}). Debido a esto, los estudios de morfometría geométrica a lo largo de gradientes latitudinales o ambientales han servido en la práctica para determinar stocks pesqueros (^{Cadrin & Friedland. 1999}; ^{Márquez, Robledo, Peñaloza, & Van der Molen., 2009}) o establecer estructuras espaciales finas de poblaciones locales o regionales (^{Schachter-broide, Dujardin, Kitron, & Gürtler, 2004}). En el presente estudio, fue una herramienta que, aplicada a la linterna de Aristóteles de Eucidaris tribuloides, sirvió como indicador de heterogeneidad ambiental.

La relación entre la longitud de la linterna de Aristóteles y el diámetro de la testa de E. tribuloides fue diferente entre los grupos analizados (norte, centro y sur del CASGM), probablemente como resultado del desarrollo ontogénico, o bien por la influencia de los factores genéticos o ambientales que pueden limitar o estimular el crecimiento de alguna de las estructuras analizadas, influyendo así en la variación alométrica entre grupos poblacionales (^{Klingenberg, 2016}). Precisamente, para el caso de los erizos de mar se ha demostrado que la presión ambiental puede influir en el crecimiento de la testa y la linterna de Aristóteles de los erizos de mar (^{Lawrence & Jangoux, 2013}).

Mediante la morfometría geométrica se observó la ausencia de una relación alométrica entre el tamaño centroide y la forma de la rótula y semipirámide de la linterna de Aristóteles de E. tribuloides. Esto sugiere que el tamaño de las estructuras analizadas no afectó la variación morfológica observada en esta investigación. Asimismo, ^{Epherra et al. (2015}) obtuvieron resultados similares para el caso de la semipirámide de Arbacia dufresnii. Sin embargo, ocurrió lo contrario para la rótula de dicha especie. Esto puede ser debido a las diferencias en los diseños de muestreo, ya que en A. dufresnii se analizaron individuos con diámetros de testa que van de 10 a 50 mm, mientras que en el presente estudio se recolectaron individuos mayores a 20 mm de diámetro de testa para evitar la variación morfológica debido al estado de madurez sexual. Con base en lo anterior, es posible que la variación morfológica en erizos de mar también se manifieste por efecto del desarrollo ontogénico (^{Frédérich et al., 2008}; ^{Díaz & Campos, 2014}).

Por otro lado, la variación morfológica de la linterna de Aristóteles para erizos herbívoros está influenciada por la disponibilidad de algas, es decir, cuando la cobertura o disponibilidad de algas es baja, la fuerza de forrajeo aumenta. Por lo tanto, la linterna de Aristóteles crece más en relación al diámetro de la testa (^{Ebert 1980}; ^{Black, Johnson, & Trendall, 1982}; ^{Levitan, 1988}). Esto podría explicar la variación morfológica observada en el presente estudio, ya que, posiblemente la disponibilidad de alimento representado por la cobertura de algas difiere entre los sitios estudiados, lo que estimularía la variación morfológica de la linterna de Aristóteles E. tribuloides.

Sin embargo, E. tribuloides es una especie omnívora (^{Lawrence & Jangoux, 2013}) y por ello, puede alimentarse de organismos blandos como: tapetes de algas, pastos marinos, esponjas y organismos intersticiales de sedimentos limosos, así como de organismos con estructuras duras, tales como: algas calcáreas, foraminíferos, corales, gasterópodos y bivalvos (^{McPherson, 1968}; ^{Kier, 1975}; Lawrence & Glynn, 1984; ^{Cutress, 1965}; ^{Santos, Coutinho, & Hajdu, 2002}). Dicho comportamiento pudo estimular el crecimiento de la linterna de Aristóteles en erizos que se distribuyeron sobre sustratos duros o compuestos principalmente por organismos con estructuras duras y limitar el crecimiento de aquellas linternas de erizos que se encontraron sobre sustrato blando o compuesto principalmente por organismos blandos (^{Hagen, 2008}; ^{Gianguzza & Bonaviri, 2013}).

Por otra parte, el CASGM tiene un gradiente de salinidad que disminuye con dirección norte-sur, así como gradientes de temperatura superficial de mar y concentración de clorofila-creciente en la misma dirección (^{Zavala-Hidalgo et al., 2003}: ^{Salas-Pérez et al., 2015}). Lo anterior puede estar relacionado con la longitud relativa de la linterna de Aristóteles y el tamaño centroide de la rótula y semipirámide de E. tribuloides que gradualmente disminuyen hacia el sur. Esto es posible, por un lado, debido las fluctuaciones de salinidad que pueden limitar el crecimiento de la testa y semipirámide de algunos erizos de mar cuando la salinidad es baja (^{Lau et al., 2009}). Por otro lado, la temperatura superficial del mar y la concentración de clorofila-a también pueden estar asociados a los cambios morfológicos de la linterna de Aristóteles de E. tribuloides. Sin embargo, no hay suficiente información que demuestre el efecto de las fluctuaciones de temperatura y concentración de clorofila-a sobre la morfología y crecimiento de la linterna de Aristóteles de los erizos de mar. No obstante, estos factores junto con los atributos propios de cada arrecife tales como el tipo de sustrato, profundidad, complejidad y la distancia a la costa, representan un papel importante en la definición del sustrato y la estructura comunitaria de los organismos bentónicos de cada arrecife (^{Arias-González et al., 2006}, 2008; ^{Jordán-Garza et al., 2017}). Esto a su vez puede influir en la dieta de E. tribuloides y consecutivamente podría tener un efecto indirecto en la variación morfología de la linterna de Aristóteles debido a la capacidad de adaptación alimenticia que presenta esta especie (^{Lawrence & Jangoux, 2013}).

Dicho lo anterior, es importante considerar que algunas muestras se capturaron en el año 2013 y otras en el año 2016. Ello, podría influir en los resultados de esta investigación, dado que las comunidades bentónicas de los arrecifes coralinos son complejos y dinámicos, por lo que la estructura de las comunidades bentónicas de los arrecifes puede variar a lo largo de un corto periodo de tiempo (Ceccarelli, Richars, Pratchett, & Cvitanovic, 2011). Esto, también podría influir en la dieta y por lo tanto en la variación morfología de la linterna de Aristóteles de E. tribuloides. Sin embargo, los resultados de esta investigación muestran diferencias entre los grupos del norte y sur, sitios que tienen proporción similar de ejemplares de colección capturados en 2013. Además, no se encontraron diferencias entre el sur y el centro, en este último, todos los ejemplares se capturaron en el año 2016. Por lo tanto, se infiere que este factor no es relevante en los resultados encontrados.

Por otro lado, la disminución del flujo genético debido a la distancia geográfica entre localidades puede dar lugar a la divergencia fenotípica (^{Slatkin, 1987}; ^{Hendry, Day, & Taylor, 2001}) detectada en la morfología de la rótula y la semipirámide de los erizos del norte (SALT) versus centro (SAV) y sur (SAT). Por otra parte, el patrón de las corrientes marinas en el corredor arrecifal veracruzano (^{Salas-Monreal et al., 2018}) conecta a las poblaciones de erizos pero en algunos casos pueden actuar como barreras para el intercambio genético entre poblaciones geográficamente separadas (^{Cowen, Gawarkiewicz, Pineda, Thorrold, & Werner, 2007}), especialmente cuando hay barreras locales, tal como se ha reportado para caso del pez Stegastes partitus que mostró un limitado intercambio genético entre subpoblaciones de los sistemas arrecifales del SAV (^{Villegas-Sánchez, Pérez-España, Rivera-Madrid, Salas-Monreal, & Arias-González, 2014}).

En contraste, las subpoblaciones del pez Abudefduf saxatilis que se distribuyen en el norte y centro de Veracruz presentan homogeneidad genética. Ello, a pesar de las diferencias morfológicas detectadas entre dichas subpoblaciones lo cual, indica que las poblaciones de esta especie están interconectadas genéticamente entre los arrecifes veracruzanos, pero el ambiente ejerce presión que afecta su morfología (^{Piñeros et al., 2015}). Esto mismo puede estar ocurriendo con E. tribuloides, debido a que se reproduce sexualmente por fecundación externa, presenta larvas planctónicas y son sedentarios en su etapa adulta (^{Hendler, Miller, Pawson, & Kier, 1995}; Lawarence & Jangoux, 2013). Dicho estilo de vida podría permitir que esta especie mantenga un buen intercambio genético entre subpoblaciones y se disperse a lo largo del CASGM, además de soportar la presión ambiental local durante su etapa juvenil y adulta, sugiriendo que el origen de la variación morfológica detectada en esta investigación se debe a la presión ambiental local y no a la baja interacción genética entre los grupos analizados.

En esta investigación, la mayoría de los individuos que representaron al centro (con 31 ejemplares recolectados) provienen del arrecife La Segunda Laja, el cual se encuentra 10 km al norte del SAV. Pese a ello, dicho arrecife tiene características similares a los arrecifes del sur, ya que es de tipo marginal, somero, con baja cobertura coralina y un sustrato compuesto principalmente por sedimento y tapetes de algas (observaciones personales). Lo anterior también puede explicar la ausencia de diferencias morfológicas entre los individuos del centro y el sur del CASGM.

En conclusión, los individuos de E. tribuloides a lo largo del Corredor Arrecifal del Suroeste del Golfo de México mostraron variación morfológica en sus estructuras analizadas. Estas diferencias morfológicas pueden estar relacionadas al diseño de muestreo, factores genéticos o bien a las múltiples variables ambientales, incluyendo la salinidad, temperatura, concentración de clorofila-a y las variables propias de cada arrecife, tales como la cobertura bentónica y el tipo de sustrato en el que se encuentran. Nosotros sugerimos que estas variables ambientales pueden tener un efecto directo o indirecto sobre la variación morfología de la linterna de Aristóteles de E. tribuloides al adaptar sus hábitos alimenticios y ajustar el rendimiento de ramoneo de acuerdo a la disponibilidad del alimento y las características del entorno particular de cada arrecife.

Finalmente, es necesario estudiar con mayor profundidad el tema del efecto geográfico, genético y ambiental sobre la variación morfológica de los organismos en CASGM, con la finalidad de tener información más precisa sobre las causas que estimulan la variación morfológica en el corredor, esto puede ser útil para identificar heterogeneidad ambiental, proponer unidades o subunidades de manejo, y (en caso de requerirse) establecer stocks pesqueros, información que es necesaria para el manejo y administración de los recursos del Corredor Arrecifal Suroeste del Golfo de México.

Declaración de ética: los autores declaran que todos están de acuerdo con esta publicación y que han hecho aportes que justifican su autoría; que no hay conflicto de interés de ningún tipo; y que han cumplido con todos los requisitos y procedimientos éticos y legales pertinentes. Todas las fuentes de financiamiento se detallan plena y claramente en la sección de agradecimientos. El respectivo documento legal firmado se encuentra en los archivos de la revista.