Introducción

El dimorfismo sexual, es decir la presencia de variaciones morfológicas debido a una presión selectiva actuando de diferente manera en cada sexo, es una de las características más conspicuas en el reino animal que aún plantea interrogantes en diferentes organismos (^{Andersson, 1994}; ^{Delph, 2005}; ^{Lande, 1980}). Para dar cuenta de este fenómeno, ^{Darwin (1871)} propuso el concepto de selección sexual, indicando que existen caracteres sexuales secundarios en los seres vivos, que se desarrollan y seleccionan sobre otros para alcanzar mayor éxito reproductivo. Los machos pueden usar los caracteres sexuales secundarios ya sea para llamar la atención de las hembras y así ser elegidos por éstas, o bien para los enfrentamientos con otros machos, previos a la cópula (^{Álvarez, 2017}). El dimorfismo sexual en el tamaño corporal es un carácter que puede encontrarse influenciado por la selección sexual, la selección natural y por un crecimiento alométrico relativo al tamaño corporal (regla de Rensch) (^{Clutton-Brock, 2017}). Según esta regla, el dimorfismo sexual en especies filogenéticamente relacionadas se acentúa con el aumento del tamaño corporal cuando el macho es el sexo de mayor tamaño, y disminuye con el aumento corporal cuando la hembra es mayor (^{Rensch, 1950}, ¹⁹⁶⁰).

La familia Gonyleptidae Sundevall, 1833 (Arachnida, Opiliones) tiene como característica más notable su marcado dimorfismo sexual, el cual se manifiesta principalmente en la pata IV, desde la coxa hasta el fémur, a veces también en patela y tibia. En los machos, este apéndice está provisto de una armadura (tubérculos, apófisis) más compleja y desarrollada, así como con mayor grosor, curvatura o longitud de cada artejo que las hembras, cuya morfología es mucho más simple (^{Kury y Pinto-da Rocha, 2007}). Dicha armadura estaría directamente relacionada con estrategias reproductivas, ya que estos opiliones utilizan la pata IV para enfrentarse con otros machos por el recurso hembras, y con una armadura más fuerte, mayores son sus posibilidades de imponerse (^{Buzatto et al., 2014}). Un fenómeno que ha llamado la atención en machos de Gonyleptidae es que la armadura muestra alometría positiva respecto del tamaño corporal (^{Acosta, 2002}; ^{Acosta et al., 2007}; ^{Buzatto et al., 2014}). El término “alometría positiva” refiere a que en los individuos de mayor tamaño (dentro de una población), determinada estructura está proporcinalmente más desarrollada que en los individuos menores (si se mantuviera la misma proporción en el rango de tamaños hablaríamos de isometría). En estos machos, se detectó que tal alometría no tiene una distribución perfectamente continua, sino que es factible reconocer al menos 2 grupos discretos, es decir, diferentes morfotipos en los machos adultos de una misma especie (^{Buzatto et al., 2014}). Este dimorfismo intrasexual observado en los machos adultos ha recibido el nombre de pecilandria (^{Regalin, 1997}; ^{Ringuelet, 1959}). En cualquier caso, cabe cuestionarse si este fenómeno refleja realmente la existencia de morfotipos discretos, o si en cambio se trata de los 2 fenotipos más abundantes o más fácilmente identificables, en una distribución poblacional continua en dichos rasgos morfológicos. Los estudios realizados hasta el momento apoyarían la primera interpretación, es decir, la de morfotipos discretos (^{Buzzato et al., 2014}; ^{Painting et al., 2015}). Este dimorfismo de machos ha sido largamente descrito e investigado en especies de coleópteros cuyos machos utilizan cuernos, mandíbulas y otras estructuras para entablar luchas por el acceso a las hembras (^{Emlen et al., 2007}; ^{Gotoh et al., 2014}; ^{Lavine y Miura, 2014}; ^{Lavine et al., 2015}; ^{Marlowe et al., 2015}; ^{McCullough et al., 2015}; ^{Vera-Cano et al., 2017}).

Se estima que los costos de producir los caracteres sexuales secundarios en machos, como ornamentaciones o armadura, pueden ser muy altos, llevando a generar variaciones discretas, con machos más grandes y ornamentados y otros más pequeños y menos llamativos (^{Clutton-Brock, 2017}; ^{Gadgil, 1972}; ^{Petrie, 1992}). En poblaciones de Gonyletpidae generalmente se identifican machos adultos grandes, con pata IV más desarrollada y/o armadura proporcionalmente más fuerte, con frecuencia llamados “machos α” o “majors”, y otros más pequeños, con armadura menos desarrollada, a veces incluso con apariencia decididamente feminiforme, los “machos β” o “minors” (^{Buzatto, 2012}; ^{Buzatto et al., 2014}). Estas diferencias son a veces tan acentuadas que pueden llevar a confusiones taxonómicas (como creer erróneamente que los machos pequeños y grandes son diferentes especies) o incluso, dar la falsa idea de que los machos menores representarían una suerte de adultos de primer estadio capaces de mudar (^{Gnaspini et al., 2004}). Si fuera cierta la supuesta existencia de 2 estadios adultos en Gonyleptidae, sería una excepción inédita en toda la Clase Arachnida.

En Gonyleptidae, cada morfotipo ha sido asociado a una estrategia reproductiva diferente: los machos α defienden un harén y/o la puesta de huevos, mientras los machos β son solitarios y aprovechan alguna distracción de los majors para inseminar hembras (^{Buzatto et al., 2014}). En 2 opiliones neozelandeses del suborden Eupnoi, Pantopsalis cheliferoides (Colenso, 1882) y Forsteropsalis pureora Taylor, 2013 (Neopilionidae) se ha descrito incluso la existencia de trimorfismo (3 morfos discretos) en los machos, en este caso afectando las proporciones de largo y grosor de artejos en los quelíceros, que son los apéndices que usan para enfrentar a otros machos (^{Painting et al., 2015}; ^{Powell et al., 2020}). Solo en este último trabajo, las mediciones fueron analizadas con morfometría geométrica, mientras en los restantes de dimorfismo en machos de opiliones, se usó morfometría lineal o tradicional.

Por la complejidad y variaciones de la armadura del macho, el goniléptido Discocyrtus prospicuus (Holmberg, 1876) y sus formas afines representan un modelo muy adecuado para abordar estas cuestiones. Según el conocimiento actual, D. prospicuus (sensu stricto) se distingue por caracteres muy sutiles de al menos 2 especies nominales, D. bucki Mello-Leitão, 1935 y D. vegetus Canals, 1939, formando en conjunto lo que denominamos el complejo D. prospicuus (nombre basado en su especie más característica). Mientras D. vegetus fue considerada una especie válida por ^{Ringuelet (1959)}, los ejemplares de Misiones que este autor refirió como D. prospicuus, fueron luego asignados a D. bucki por ^{Acosta y Guerrero (2011)}. Observaciones inéditas muestran que la distribución de D. bucki se prolonga desde Misiones, a lo largo del río Uruguay, hasta los 29°43’ S en la provincia de Corrientes (fig. 1), con una aparente y sutil variación clinal N-S de algunos caracteres. A partir de las poblaciones de Entre Ríos hacia el sur, se detecta la morfología típica de D. prospicuus s.s.

Figura 1 Ubicación de las 11 localidades (círculos rojos) de donde provienen los ejemplares estudiados. Círculos pequeños: distribución conocida del complejo Discocyrtus prospicuus: D. prospicuus s.s., círculos celestes; D. bucki, círculos blancos; D. vegetus, círculos amarillos. Números en recuadros: las 4 zonas de la distribución del complejo D. prospicuus referidas en el texto. 

Junto con D. dilatatus Sørensen, 1884 y D. testudineus (Holmberg, 1876), D. prospicuus s.s. es una de las 3 especies del género Discocyrtus Holmberg, 1878 que define los límites occidental y austral del extenso género, con mayoría de especies en el este de Brasil (^{Kury, 2003}). Las 3 muestran una amplia distribución en el área opiliológica mesopotámica “sensu stricto” de Argentina, definida por ^{Acosta (2002)}. Dicha área comprende las extensas planicies, húmedas y semihúmedas, situadas entre los ríos Paraná y Uruguay, así como una franja al oeste del primero; contiene, en conjunto, un mosaico de arbustales, bañados, pastizales y bosques en galería (^{Acosta, 2002}). En el extremo nordeste de la mesopotamia (provincia de Misiones), la fisonomía es de selva subtropical y posee componentes opiliofaunísticos distintivos, siendo reconocida como un área diferente, el área misionera o paranense (^{Acosta, 2002}). Según los registros conocidos hasta el momento (^{Acosta y Guerrero, 2011}; obs. inéd.), el complejo D. proscipuus, tomado en conjunto, se concentra en Argentina en al menos 4 zonas, indicadas en figura 1. Las 3 primeras corresponden a lo que llamamos D. prospicuus s.s., en coincidencia con el concepto taxonómico adoptado por ^{Acosta y Guerrero (2011)}. Estas zonas son: 1) un área principal rodeando el estuario del río de la Plata, el delta del río Paraná y la parte baja del río Uruguay, tanto del lado argentino (hasta Entre Ríos) como uruguayo, 2) un sector asociado a las Yungas en el NOA (= noroeste argentino, provincias de Tucumán, Salta y Jujuy), y 3) una zona central que abarca las planicies de la provincia de Córdoba, hasta las primeras estribaciones serranas. La zona 4) de este complejo de especies, se ubica en la provincia de Misiones, con la mencionada prolongación sobre el borde oriental de la vecina provincia de Corrientes, sector que contiene a D. vegetus y D. bucki (fig. 1). Sin abrir juicio sobre la identidad de las 3 especies nominales (sin dudas muy próximas), en este trabajo se estudian poblaciones asignadas a 2 entidades del mencionado complejo D. prospicuus.

A pesar de la compleja estructura que ofrecen los machos de Gonyleptidae, y su aparente variabilidad, no se ha aplicado hasta el momento la técnica de morfometría geométrica (MG) para abordar diversas problemáticas, como probar la existencia de esta variación y de posibles morfotipos discretos con diferentes tamaños centroide y conformación. Ésta es una herramienta matemática con fundamentos biológicos que ha demostrado ser muy efectiva para descomponer la variación resultante de la fisiología del individuo, de aquella más estable, propia de una población y producto principalmente de un componente genético (^{Dujardin, 2008}). Debido a ésto, la MG ha logrado resolver exitosamente problemas taxonómicos, ha sido utilizada para identificar marcadores de hábitats, en estudios de estructuración poblacional, de dispersión y diferenciación sexual, en arácnidos e insectos entre otros grupos (^{Bechara y Liria, 2012}; ^{Carvalho et al., 2010}; ^{Crews y Hedin, 2006}; ^{Dujardin, 2008}, ²⁰¹¹; ^{Hernández et al., 2011}, ²⁰¹³, ²⁰¹⁵, ²⁰¹⁸, ²⁰²⁰; ^{Pretorius, 2005}; ^{Sasakawa, 2016}).

Si bien existen varias hipótesis que intentan explicar estas diferencias morfométricas, las investigaciones en Gonyleptidae, hasta el momento, han abordado el problema a través de medidas lineales de una o pocas dimensiones corporales (como largo del escudo, del fémur IV, quelíceros, etc.), o más frecuentemente de la covariación de tales caracteres (e.g., ^{Buzzatto et al., 2014}; ^{Painting et al., 2015}; ^{Zatz et al., 2011}). El único antecedente de uso de MG en Opiliones corresponde a otra familia (^{Powell et al., 2020}). La medición de caracteres lineales deja de lado la geometría de la estructura corporal medida y no tiene en cuenta que dicha estructura puede sufrir cambios de manera diferente en cada punto anatómico de referencia. Como se dijo, no existen, de momento, estudios en Gonyleptidae que utilicen MG basada en puntos de referencia para intentar describir estas variaciones, aplicarlas en la delimitación de especies y explicar las causas de la variación.

El presente trabajo busca evaluar la potencialidad de la MG para determinar el dimorfismo sexual en el tamaño y la conformación de 2 caracteres: el escudo dorsal y la coxa IV de machos y hembras dentro del complejo D. prospicuus. También se procura evaluar la existencia de pecilandria a través de la variabilidad morfométrica del escudo dorsal y coxa IV de los machos, así como la influencia geográfica sobre la variabilidad morfométrica y el dimorfismo sexual. El desarrollo de la coxa IV es un carácter muy conspicuo, previamente usado en la identificación con medidas lineales de los 2 tipos de dimorfismo en machos de Opiliones (^{Buzatto et al., 2014}). En cuanto al escudo, solamente se ha utilizado como referencia para estandarizar medidas lineales de la pata IV en el análisis de cambios alométricos (^{Willemart et al., 2009}; ^{Zatz et al., 2011}). El contorno del escudo ha sido utilizado mayormente con enfoque taxonómico pero solo como un carácter cualitativo, identificándose diversas tipologías que contemplan las dimensiones relativas de carapax y mesotergo, curvaturas laterales de este último y el aspecto general (^{Kury y Medrano, 2016}). Este es, por tanto, el primer intento de abordar esta estructura desde la MG.

Materiales y métodos

La identificación taxonómica de los ejemplares utilizados en los análisis fue realizada siguiendo los criterios de ^{Acosta (1999)}, ^{Acosta y Guerrero (2011)}, ^{Capocasale (1966)} y ^{Ringuelet (1959)}. Este material, conservado en etanol 70%, forma parte de la colección aracnológica de la Cátedra de Diversidad Biológica II, Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba (CDA). Los ejemplares corresponden a 3 provincias argentinas (Entre Ríos, Córdoba y Corrientes), representativas, respectivamente, de los sectores 1, 3 y 4 según se definen más arriba y se muestran en la figura 1. Se trabajó con 113 individuos (58 machos y 55 hembras) de 11 localidades (4 localidades en Entre Ríos: n = 18♂♂ 21♀♀, 4 en Córdoba: n = 13♂♂ 24♀♀, 3 en Corrientes: n = 27♂♂ 10♀♀). El detalle de localidades (con coordenadas), número y sexo de los ejemplares examinados se presenta en el apéndice 1. Para los análisis de dimorfismo sexual se trabajó con la totalidad de machos y hembras como un único grupo, es decir sin identificar su procedencia. En los análisis para evaluar posible dimorfismo en machos se tuvo en cuenta el sector (Corrientes, Entre Ríos o Córdoba) a fin de poner a prueba si esta característica sufre modificaciones según el área geográfica.

El mapa de localidades fue diseñado con el software QGIS 2.8.1 - Wien (Sistema de Información Geográfica de código abierto, https://qgis.org/), utilizando capas espaciales libremente disponibles en http://www.diva-gis.org/Data y https://www.worldwildlife.org/publications/terrestrial-ecoregions-of-the-world. Para los análisis morfométricos (de escudo y coxa IV), los ejemplares fueron colocados en una cápsula de Petri con etanol 70% cubriendo completamente su cuerpo. Cada ejemplar fue fotografiado por sus caras dorsal y ventral con ayuda de un estereomicroscopio Leica M5 (10X) equipado con cámara digital Leica EC3. Las fotografías fueron capturadas, procesadas y guardadas en formato .jpg utilizando el software Leica Application Suite (versión 3.4, http://www. leica-microsystems.com/education). Se seleccionaron 9 puntos de referencia tipo 2 y 3 (definidos según el nivel de precisión; ^{Bookstein, 1991}) para el escudo dorsal (fig. 2A) y 8 puntos de referencia tipo 2 para la coxa de la cuarta pata, en vista ventral (fig. 2B). Las variables de conformación fueron obtenidas por medio del análisis generalizado de Procrustes y la subsiguiente proyección de los residuos en el espacio euclidiano (^{Rohlf, 1999}). Los componentes no uniformes y uniformes fueron usados como variables de conformación. Los componentes uniformes describen la variación global como estiramientos y compresión, y los componentes no uniformes corresponden a la variación local (^{Zelditch et al., 2004}). Estos 2 componentes describen las diferencias en la conformación como desviaciones de la configuración de puntos de referencia promedio. Para comparaciones del tamaño del escudo y de las coxas se utilizó el estimador denominado tamaño centroide (CS, por sus siglas en inglés), derivado de las coordenadas. Éste es definido como la raíz cuadrada de la suma de las distancias al cuadrado entre el centro de la configuración de puntos de referencia y cada punto individual (^{Bookstein, 1991}). La varianza en la conformación de escudo dorsal y coxas IV se determinó utilizando el índice de disparidad métrica (MD). Este índice proporciona un valor único que mide la variabilidad en la conformación en una muestra y se basa en las distancias euclidianas de cada forma con respecto al consenso (^{Zelditch et al., 2004}).

Figura 2 Puntos de referencia utilizados para el escudo dorsal (A, vista dorsal) y para la coxa IV (B, vista ventral) de los ejemplares del complejo Discocyrtus prospicuus. En blanco: puntos de referencia tipo 2; en rojo: puntos de referencia tipo 3 (res. 2,048×1,536 pixeles, 24 bits por pixel). 

Dos diferentes paquetes de software fueron usados para la digitalización, análisis y obtención de resultados: el paquete CLIC version 97 (^{Dujardin y Slice, 2007}) (https://xyom-clic.eu) y el reciente paquete morfométrico online, XYOM (https://xyom.io). La digitalización de puntos de referencia, el análisis generalizado de Procrustes y los análisis de tamaño centroide y conformación se realizaron con los apartados “digitization”, “characterization” y “miscellaneous” del paquete XYOM. Los análisis de alometría se realizaron con el módulo CLIC49. Se utilizó el análisis de función discriminante (DFA, por sus siglas en inglés) para comparaciones entre grupos de diferentes provincias. La significancia estadística de los análisis de tamaño y conformación fue evaluada usando un test no paramétrico basado en permutaciones (5,000 ciclos).

Resultados

Los primeros análisis permitieron la caracterización del dimorfismo sexual. El tamaño centroide (CS) de la coxa IV resultó estadísticamente diferente entre sexos. Las hembras mostraron menor tamaño en sus coxas IV (mediana: 7.29, desviación estándar DS: 0.40) en relación a los machos (mediana: 12.94, desviación estándar: 1.57) p < 0.0001 (fig. 3A). Además, en el grupo de los machos se observaron 2 subgrupos con diferentes tamaños, mientras que en las hembras se observa muy poca dispersión en los valores de CS de la coxa IV y ausencia de subgrupos como los observados en machos.

Figura 3 A) Valores medios del tamaño centroide (CS) de coxas IV para machos y hembras del complejo D. prospicuus; los puntos representan los valores de CS de cada individuo. Las letras distintas (a, b) indican diferencias estadísticamente significativas para hembras y machos (p < 0.05). B) Variación de la conformación de la coxa IV; la figura muestra los objetos promedio luego del análisis generalizado de Procrustes y el desplazamiento de los puntos de referencia para cada sexo. El valor cero en los ejes indica el centroide de la estructura luego del AGP. H: Hembras, M: machos. 

El análisis discriminante con las variables de conformación de la coxa IV mostró diferenciación significativa de los grupos de machos y hembras (p < 0.001, 1,000 permutaciones), con una distancia de Mahalanobis de 7.59. Los puntos de referencia de esta estructura muestran importantes movimientos entre machos y hembras, principalmente en los relacionados con el ancho de la coxa y con la presencia de apófisis en machos y que son prácticamente inexistentes en hembras (fig. 3B).

Con respecto al tamaño centroide del escudo dorsal, este parámetro resultó mayor en hembras (mediana: 12.95, desvío: 1.76) en relación a machos (mediana: 11.88, desvío: 1.53), p < 0.001 después de 1,000 permutaciones. El gráfico de cajas muestra que tanto en la distribución de machos como de hembras, se observan 2 grupos discretos bien diferenciados, con mayor y menor tamaño centroide (fig. 4A).

Figura 4 A) Valores medios del tamaño centroide (CS) del escudo dorsal para machos y hembras del complejo D. prospicuus; los puntos representan los valores de CS de cada individuo. Las letras distintas (a, b) indican diferencias estadísticamente significativas para hembras y machos (p < 0.05). B) Variación de la conformación del escudo dorsal; la figura muestra los objetos promedio luego del análisis generalizado de Procrustes y el desplazamiento de los puntos de referencia para cada sexo. El valor 0 en los ejes indica el centroide de la estructura luego del AGP. H: Hembras, M: machos. 

En cuanto a la conformación del escudo dorsal, el gráfico de objetos promedio muestra que los puntos de referencia con mayores cambios son los relacionados con el ancho máximo y posterior del escudo (a la altura de la apófisis coxal): 5 y 6, 7 y 8. En la región delimitada por esos puntos de referencia (fig. 2A), las hembras evidenciaron escudos mas angostos que los machos si se considera el ancho máximo, pero con mayor ancho posterior (fig. 4B). El análisis multivariado discriminante muestra diferencias significativas (p < 0.001, usando la corrección de Bonferroni) en la conformación del escudo de machos y hembras. La distancia de Mahalanobis entre grupos fue de 2.115.

En cuanto a los análisis de dimorfismo en machos y variación geográfica, los resultados muestran que el tamaño centroide de coxa IV resultó significativamente mayor en individuos de Corrientes (12.36) en relación con los de una área contigua, Entre Ríos (10.68, p < 0.001). El grupo de Córdoba, aunque levemente menor que Corrientes, no se diferencia significativamente de ninguno de los otros grupos (mediana: 11.57) (fig. 5A). En los grupos de Córdoba y Corrientes, se observan 2 subgrupos marcados en cuanto al tamaño centroide, evidenciando machos con coxas de mayor tamaño y otro con coxas menores dentro de estas provincias. El grupo de Entre Ríos no mostró una gran dispersión ni subgrupos marcados en cuanto a esta variable (fig. 5A).

Figura 5 Valores medios del tamaño centroide (CS) para machos del complejo D. prospicuus por zona de procedencia. A) Coxa IV, B) escudo dorsal. Las cajas muestran valores de medianas, mínimos y máximos. Los puntos representan los valores de CS de cada individuo. Las letras diferentes (a, b) indican diferencias estadísticamente significativas para los grupos (p < 0.05). 

Con respecto a la conformación de la coxa IV en machos, los primeros ejes del análisis DFA explican el 92% de la variación (51% y 41% para el primer y segundo factor, respectivamente). El mapa factorial muestra un solapamiento de los 3 grupos, y una ausencia de diferenciación significativa entre machos de los grupos de Corrientes, Entre Ríos y Córdoba. Las distancias de Mahalanobis, aunque no muestran diferencias significativas entre grupos, resultaron mayores entre los grupos de Córdoba y Entre Ríos (1.92) y menores entre los grupos de Corrientes y Entre Ríos (1.49) (fig. 6A). El efecto alométrico fue de 5% para las variables canónicas CV1 y 13% CV2. Eliminando el efecto alométrico las diferencias entre los machos de las 3 provincias resultaron significativamente diferentes (p < 0.05). La disparidad métrica (MD) resultó menor en el grupo de Corrientes (0.0029) diferenciándose significativamente de Córdoba (0.0046), quien muestra la mayor DM. Entre Ríos muestra un valor intermedio de DM, diferenciándose también de Córdoba y Corrientes con valores significativos (0.0036), (p < 0.005, 5,000 bootstraps).

Figura 6 Análisis discriminante con variables de conformación de coxas (A) y variables de conformación del escudo dorsal de machos del complejo D. prospicuus por provincia de procedencia (B). Cada punto representa un individuo en el análisis. CV1 y CV2 corresponden a los ejes discriminantes. 

Al igual que sucede con la variable coxa IV, el grupo de Entre Ríos (mediana: 11.36) es el que muestra menor dispersión al analizar el tamaño centroide del escudo dorsal, sin diferenciación de subgrupos (fig. 5B). Córdoba presenta menor valor de mediana (11.24), y al igual que en el grupo de Corrientes (mediana: 13.19), se observan subgrupos discretos con mayores y menores valores de tamaño centroide del escudo dorsal. Las diferencias de CS resultaron estadísticamente significativas entre los grupos de Córdoba y Corrientes, y entre los de Corrientes y Entre Ríos (p < 0.05). Por su parte, Córdoba y Entre Ríos no difieren en el tamaño centroide del escudo dorsal (fig. 5B).

En cuanto a la conformación del escudo dorsal, los primeros factores discriminantes explican 94 % de la variación (64% y 30% para CV1 y CV2). El análisis discriminante no encuentra diferencias estadísticamente significativas entre los grupos de las 3 provincias para el escudo dorsal. (p < 0.05, 1000 ciclos). La menor distancia de Mahalanobis se da entre los grupos de Córdoba y Corrientes (0.884) mientras que Corrientes y Entre Ríos muestran la mayor distancia (1.65) (fig. 6B). El efecto alométrico es del 6%. Extrayendo el efecto alométrico, las diferencias entre grupos no resultaron estadísticamente significativas. La disparidad métrica no resultó significativamente diferente entre grupos (Córdoba 0.0109, Corrientes 0.0092, Entre Ríos 0.0108).

Discusión

Las diferencias entre machos y hembras son el resultado evolutivo de presiones selectivas actuando diferencialmente entre sexos. Las presiones selectivas pueden ser intra o intersexuales, aunque a veces presentan una función ecológica diferente, que es reducir el solapamiento de nichos y la competencia entre sexos (^{Pianka, 1982}). En nuestros análisis se pudo observar un marcado dimorfismo sexual en el tamaño y la conformación de la coxa IV y del escudo dorsal, aunque en cada estructura este dimorfismo se mostró en distinto sentido. Las coxas IV revelaron menor tamaño centroide en hembras en relación a machos, mientras para el escudo dorsal las hembras fueron las que mostraron mayor CS. Estos resultados evidencian que las fuerzas de canalización (es decir, la capacidad de un genotipo de producir el mismo fenotipo a pesar de la variabilidad de su entorno) actúan de manera diferente en las 2 estructuras medidas en este trabajo. Observamos que la canalización es más fuerte en el escudo en comparación con las coxas IV. Estas últimas resultaron las estructuras con mayor variabilidad morfométrica, sobre todo en los machos; ellas son un elemento usado en los enfrentamientos con otros machos por los apareamientos con hembras. Por otro lado, el escudo es una estructura donde la evolución priorizaría un desarrollo más armónico a lo largo de su ciclo de vida, en comparación con las coxas IV, en las cuales la importancia de su desarrollo se relaciona al estadio adulto (para competencia de machos por luchas intrasexuales). Las diferencias en la dirección de desarrollo en coxas y escudo apoyan la hipótesis de canalización genética, como se ha sugerido también para mediciones morfométricas de cabeza y alas en insectos (^{Caro-Riaño, 2009}; ^{Dujardin, 2011}; ^{Hernández et al., 2011}), y en experimentos de plasticidad fenotípica (^{Hallson y Bjӧrklund, 2012}). ^{Fairbairn (2005)}, mostró la existencia de una plasticidad diferencial en machos y hembras de Aquarius remigis (Say, 1832) (Hemiptera, Gerridae). Este tipo de plasticidad plantea una canalización adaptativa sobre los rasgos más estrechamente asociados con la aptitud reproductiva (por ejemplo, ancho del abdomen en hembras, longitud genital en machos, etc.) como una hipótesis alternativa a la regla de Rensch.

Por otro lado, los resultados aquí expuestos demuestran el dimorfismo en machos (pecilandria) para las 2 estructuras medidas, en al menos 2 de las 3 zonas geográficas estudiadas. En nuestros datos, los machos evidenciaron 2 subgrupos bien marcados en cuanto al tamaño centroide, tanto de las coxas IV como del escudo dorsal. ^{Buzatto et al. (2014)} encontraron que la selección sexual puede llevar a evitar el gasto en rasgos sexuales secundarios en algunos machos. De esta manera, la existencia de machos alfa y beta, con coxas de mayor y menor tamaño, o las diferencias en el tamaño de los escudos, pueden estar indicando habilidades diferentes para conseguir llegar a la hembra y finalmente aparearse (^{Gross, 1996}; ^{Shuster y Wade, 2003}; ^{Taborsky et al., 2008}). El comportamiento de cópulas furtivas en machos pequeños fue registrado en Opiliones para Longiperna concolor (Mello-Leitão, 1923) (Gonyleptidae: Mitobatinae), donde parece ser parte de una táctica reproductiva alternativa a la defensa de los recursos que se da en los machos a (^{Zatz et al., 2011}). De modo similar, ^{Powell et al. (2020)} lograron dilucidar un caso de trimorfismo en machos de Neopilionidae, utilizando morfometría geométrica de quelíceros, que en esta familia son los apéndices utilizados en los enfrentamientos entre machos. Según sus observaciones, cada morfotipo tendría una estrategia diferente en el uso de sus quelíceros durante las luchas entre machos (^{Powell et al., 2020}).

Aunque el dimorfismo en machos se detectó en las poblaciones de Corrientes y Córdoba, no pudo observarse en los ejemplares de Entre Ríos, en donde la variabilidad fenotípica en los rasgos medidos resultó muy baja. Esta característica puede estar indicando diferencias geográficas dadas por el ambiente. Es posible que debido a una estrategia reproductiva o a una limitante ambiental (competencia, predación, variación fenotípica local dada por el ambiente) en los ejemplares de Entre Ríos, el dimorfismo de machos no sea evidente, o no sea una estrategia adaptativa preferida para esta región geográfica.

Para finalizar, los resultados mostraron que al analizar el escudo dorsal, las hembras también exhibieron subgrupos discretos con diferentes tamaños, algo que no se ha descripto previamente en hembras de Opiliones. El escudo dorsal es un indicador del tamaño corporal, por lo que los resultados estarían indicando diferencias de tamaño dentro del grupo de las hembras que podrían deberse a diferencias ambientales regionales. Un examen detallado evidencia que las hembras más pequeñas que se observan en los gráficos de tamaño centroide (fig. 4A), corresponden a ejemplares recolectados en Entre Ríos. Este resultado puede relacionarse con lo observado en machos de Entre Ríos (baja variabilidad en CS de coxa y escudo) y ausencia de dimorfismo. Esto sugiere diferencias regionales en Entre Ríos que impactan directamente en el tamaño de los ejemplares, tanto machos como hembras. Se sabe que el tamaño centroide es una variable altamente influenciada por el estado nutricional y por la temperatura (regla de Bergmann). Es posible que el tamaño en D. prospicuus de Entre Ríos tenga una influencia ambiental local, se trate de poblaciones que habrían sufrido un efecto fundador o presenten muy bajo movimiento dispersivo con poblaciones vecinas. En balance, los ejemplares de Corrientes (“D. bucki”) tienden a diferenciarse más acentuadamente de los de otras zonas (por lo general hacia el mayor tamaño), si bien no estamos en condiciones de abrir juicio sobre su posible conespecificidad con los otros grupos del complejo.

Por su parte, se analizan por primera vez los caracteres morfológicos y morfométricos en hembras de Opiliones, caracterizando su variación. En cuanto al uso de la técnica, los resultados permiten concluir que la aplicación de MG logró la identificación de morfotipos discretos en machos del complejo D. prospicuus. Se identificó por primera vez cómo ocurre la variación de las coxas IV y escudo dorsal, proponiendo una variabilidad diferencial en cada una de las estructuras medidas. Este trabajo demuestra que es posible dilucidar la complejidad morfológica presente en D. prospicuus (sensu lato) a través de MG, lo cual constituye un avance en la comprensión de la ocurrencia de dimorfismos y las variaciones geográficas asociadas con esta variación.