Introducción
Mugil incilis es un pez eurihalino conocido comúnmente como lisa en Colombia. Habita ecosistemas acuáticos costeros con conexión al mar en el Caribe, desde Haití y Panamá bordeando la costa atlántica occidental de Suramérica hasta el sureste de Brasil (Bustos-Montes et al., 2012; Harrison, 2002). En Colombia está distribuida desde Bahía de Cispatá hasta Bahía Portete (Rey & Acero, 2002). Tiene una importancia ecológica por su rol trófico en los estuarios y lagunas costeras dado que es un pez detritívoro que al ser depredado ofrece transferencia de materia y energía a niveles tróficos superiores; de igual manera las transfiere al mar adyacente durante el periodo de migración reproductiva (Blanco-Racedo, 1983; Osorio, 1998; Cogua, Jiménez-Reyes, & Duque, 2013). También tiene una importancia pesquera por ser aprovechada comercialmente al registrarse capturas que aportan cerca del 25 % de la producción pesquera de la Ciénaga Grande de Santa Marta (INVEMAR, 2014) y el 3 % del Caribe de Colombia (De La Hoz-Maestre et al., 2013). A pesar de su importancia ecológica y comercial, estudios ecológicos han documentado que los mayores registros de captura se presentan cuando sale de las lagunas costeras hacia el mar adyacente para reproducirse. Esto limita la renovación poblacional puesto que los ejemplares capturados generalmente están maduros; dentro de las lagunas y estuarios es capturada por debajo de la talla media de madurez sexual con efectos importantes en su abundancia y supervivencia (Blanco-Racedo, 1983; Sánchez-Ramírez, Rueda, & Santos-Martínez, 1998; Narváez-Barandica, Herrera-Pertuz, & Blanco-Racedo, 2008; INVEMAR, 2014). A lo anterior se suman otros factores que afectan su biología poblacional como la calidad ambiental de la mayoría de las aguas de las lagunas y estuarios donde habita debido a la presencia de contaminantes de origen industrial, urbano y agrícola, principalmente en la Ciénaga Grande de Santa Marta, Bahía de Cispatá y la Ciénaga de la Virgen (INVEMAR, 2017). Por ejemplo, estudios de contenidos de organoclorados y metales pesados de M. incilis documentaron la presencia de estos contaminantes en peces de esas dos lagunas costeras (Plata, Campo & Ramírez, 1993; Alonso, Pineda, Olivero, González, & Campos, 2000; Burgos-Núñez, Navarro-Frómeta, Marrugo-Negrete, Enamorado-Montes, & Urango-Cárdenas, 2017). Como consecuencia de lo anterior, el estrés ambiental está incidiendo en altos niveles de prevalencia de parásitos en M. incilis como una indicación del estado vulnerable de sus poblaciones (Bustos-Montes et al., 2012; Jaramillo-Colorado, Arroyo-Salgado, & Ruiz-Garcés, 2015).
Con referencia a todo lo anterior, M. incilis es un modelo biológico interesante para conocer la influencia del comportamiento reproductivo en la dispersión y conectividad genética de las poblaciones de peces en el Caribe de Colombia. Así mismo, por tratarse de una especie sobreexplotada que habita ecosistemas marino-costeros con altos niveles de contaminación puede revelar cómo los peces están enfrentando una reducción de los niveles de variabilidad genética y así permitir dilucidar cuáles son las capacidades para superponerse a los cambios del medio ambiente. Estudios previos han demostrado que mientras mayor sea la variación genética de las poblaciones, mejor será la respuesta de una especie a los cambios ambientales (Bickham, Sandhu, Hebert, Chikhi, & Athwal, 2000; Belfiore & Anderson, 2001; Ross, Cooper, & Bidwell, 2002; McMillan, Bagley, Jackson, & Nacci, 2006). Aunque mucho se conoce sobre la ecología y el impacto en el sector pesquero de esta especie en el Caribe de Colombia, no existen estudios genéticos que describan el estado de conservación de las poblaciones en ambientes naturales. Por lo anterior, en este estudio se seleccionaron siete loci microsatélites con el objetivo de determinar los patrones de distribución de la variabilidad genética y el nivel de estructuración genética poblacional de M. incilis en áreas de media y alta presión pesquera del Caribe de Colombia. Esto con el fin de proponer nuevos criterios científicos que permitan desarrollar futuros programas de manejo pesquero para este pez.
Materiales y métodos
Muestreo: Las muestras de M. incilis fueron recolectadas en seis sectores del Caribe colombiano donde se concentra la mayor actividad pesquera que aprovecha esta especie; tres sistemas lagunares costeros de acceso abierto a la pesca (Ciénaga La Virgen = CVIR, Ciénaga La Caimanera = CCAI y Bahía de Cispatá = BCIS) y tres del sistema de Parques Nacionales Naturales de Colombia (Santuarios de Fauna y Flora Los Flamencos = SFFF, Ciénaga Grande de Santa Marta = CGSM y El Corchal = SFFC) (Tabla 1 y Fig. 1). En total 249 muestras fueron recolectadas entre marzo y julio de 2014. La muestra de tejido fue de la aleta caudal y fueron depositados en alcohol al 96 % y congelados a -20 °C para su posterior procesamiento.
Abreviación | Localidad | N | Coordenadas | Departamento |
SFFF | Santuario de Fauna y Flora Los Flamencos | 50 | 11°22'48.5" N 73°07'23.2" W | La Guajira |
CGSM | Ciénaga Grande de Santa Marta | 40 | 10°52'27.3" N 74°24'56.3" W | Magdalena |
CVIR | Ciénaga La Virgen | 39 | 10°26'14.8" N 75°29'55.0" W | Bolívar |
SFFC | Complejo cenagoso del Santuario de Fauna y Flora El Corchal | 40 | 10° 02´00" N 75°32´00" W | Sucre |
CCAI | Ciénaga La Caimanera | 40 | 9°25'58.3" N 75°37'53.8" W | Sucre |
BCIS | Bahía de Cispatá | 40 | 9°24'08.7" N 75°49'29.8" W | Córdoba |
N= número de muestras recolectadas por sitio.
N= number of samples collected per site.
Extracción ADN y amplificación de microsatélites: El ADN genómico fue extraído utilizando el protocolo de sales de Lopera-Barrero et al. (2008). Siete microsatélites loci previamente descritos para Mugil cephalus fueron seleccionados en la genotipificación de las muestras de M. incilis: Mcs-15 AM (Miggiano et al., 2005), Mce-2, Mce-3, Mce-7, Mce-11 (Shen et al., 2010), Muce-55 y Muce-80 (Xu, Sun, Shi, & Wang, 2010). Es importante resaltar que la amplificación de estos microsatélites fue exitosa en especies hermanas como Mugil liza (Mai et al., 2014) y Mugil curema (Pacheco-Almanzar, Ramírez-Saad, Velázquez-Aragón, Serrato, & Ibáñez, 2017). Por este motivo, para aceptar y trabajar con los siete microsatélites citados se realizó una amplificación con la reacción en cadena de la polimerasa (PCR) con las siguientes condiciones: volumen final de 10µl que contenía una concentración de buffer de reacción en 1 X, MgCl2 2 mM, dNTPs 0.2 mM, 0.2 μM de cada primer y taq-polimerasa 0.22 U (Bioline Meridian Life Science).
La reacción de PCR incluyó un primer ciclo de desnaturalización a 95 °C por 5 minutos; luego 30 ciclos a 94 °C por 30 segundos, una temperatura de annealing de 62 °C por 30 segundos y una temperatura de extensión a 72 °C por 90 segundos. La reacción finalizó con una temperatura de extensión final a 72 °C por 30 minutos. Los productos de la PCR fueron verificados en geles de agarosa al 2 % y corridos a 80 V por 1.5 horas en una cámara de electroforesis horizontal. Luego fueron analizados a través de electroforesis capilar con el equipo QIAxcel Advanced (QIAGEN). El tamaño de cada locus amplificado se determinó con el programa computacional QIAxcel ScreenGel v1.0, utilizando un marcador de alineamiento entre 15 y 1 000 pb.
Análisis de datos: Se detectó la presencia de alelos nulos y la identificación de posibles errores de genotipificación (amplificación como bandas “stutter” y dominancia de alelos pequeños) utilizando el programa MICRO-CHECKER v2.2.3 (Van Oosterhout, Hutchinson, Wills, & Shipley, 2004). El programa GenAlEx v.6.5 (Peakall & Smouse, 2012) se utilizó para determinar la variabilidad genética intrapoblacional con los estimativos de: número de alelos por locus, frecuencias de alelos y las heterocigosidad esperada (He) y observada (Ho). Con GENEPOP v4.0 (Rousset, 2008) se calculó el coeficiente de endogamia (FIS) (Weir & Cockerham, 1984), el desequilibrio de ligamiento (LD) y se probaron desvíos de proporciones genotípicas en el equilibrio de Hardy-Weinberg (E-HW). Se ajustó el nivel de alfa < 0.05 para comparaciones múltiples con el procedimiento de Bonferroni (Rice, 1989).
La reducción reciente en el tamaño de la población se verificó usando BOTTLENECK v1.2.02 (Piry, Luikart, & Cornuet, 1999) con 1 000 repeticiones. Para probar si existían evidencias de un cuello de botella genético, primero se consideraron dos modelos mutacionales para los microsatélites, un modelo de mutación gradual (SMM) y uno de dos fases (TPM) con una proporción de 70 % del modelo SMM y 30 % de varianza. La significancia del resultado fue evaluada con una prueba de Wilcoxon (Wilcoxon, 1945) que evaluó el exceso o déficits de individuos heterocigotos en caso de un cuello de botella reciente. Estas pruebas se realizaron utilizando todos los loci y todas las localidades.
El análisis de la estructura poblacional se analizó con el estadístico FST de Wright entre pares de localidades. Para esto se realizó un análisis molecular de varianza (AMOVA) con el programa computacional el programa Arlequin v3.5 (Excoffier & Lischer, 2010). Este análisis permitió evaluar el nivel de diferenciación en diferentes niveles geográficos de organización de las muestras. Los niveles de agrupamiento que se emplearon fueron: un primer análisis asumiendo a todas las localidades como un sólo grupo jerárquico; en el segundo se asumió un grupo para las muestras recolectadas al suroeste de la desembocadura del río Magdalena (CVIR, SFFC, CCAI, BCIS) contrastadas con las del noreste del río (CGSM y SFFF); un tercer análisis asumiendo las muestras recolectadas al noreste del Golfo de Morrosquillo (CVIR, SFFC, CGSM y SFFF) versus las que pertenecen al golfo (CCAI y BCIS). Como complemento a lo anterior, se realizó un análisis de comparación múltiple entre las localidades con el test exacto de diferenciación genotípica (Raymond & Rousset, 1995). Una vez obtenidos los resultados las comparaciones múltiples fueron ajustadas con una prueba de Bonferroni.
Por otro lado, se realizó una prueba de Mantel en IBDWS v3.23 (Jensen, Bohonak, & Kelley, 2005) para examinar la correlación entre la distancia genética y la distancia geográfica (medida en kilómetros a lo largo del borde de la costa) y evaluar si la población de M. incilis se ajusta al modelo de aislamiento por distancia. Para analizar las relaciones genéticas entre las localidades se construyó un dendrograma en Mega v6.0 (Tamura, Stecher, Peterson, Filipski, & Kumar, 2013), utilizando el método de agrupamiento Vecinos más cercano con la matriz de distancia genética (δµ2) de (Goldstein, Ruiz Linares, Cavalli-Sforza, & Feldman, 1995). Esta última se determinó en el programa Arlequin V3.5 (Excoffier & Lischer, 2010).
Se evaluó la presencia de grupos genéticos diferentes utilizando el método de agrupamiento bayesiano implementado en STRUCTURE v2.3.3 (Hubisz, Falush, Stephens, & Pritchard, 2009). Para estimar el número de grupos de poblaciones genéticas (K) dentro del conjunto completo de datos, seis simulaciones independientes para K = 1-6 se realizaron con un burn-in de 200 000, 800 000 iteraciones y diez réplicas. El análisis se realizó utilizando un modelo de mezcla (Admixture Model) de estructura de poblacional, y se correlacionó la frecuencia de alelos entre las poblaciones. Se estimó el número más probable de grupos de población mediante el estadístico de Evanno, Regnaut y Goudet (2005) utilizando la aplicación en línea STRUCTURE HARVESTER (Earl & VonHoldt, 2012).
Se aplicó el análisis discriminante exploratorio de componentes principales (DAPC), utilizando el paquete Adegenet (Jombart, 2008; Jombart, Devillard, & Balloux, 2010), el cual cuenta con la función DAPC en el software R (R Core Team, 2016). Este análisis se realizó sin información previa sobre las poblaciones individuales. El número de grupos se evaluaron mediante la función find.clusters, que corre sucesivos promedios del K con el aumento de número de agrupaciones (K). Para seleccionar el número óptimo de las agrupaciones, se aplicó el criterio de información bayesiano (BIC) para evaluar el modelo con mayor soporte (Jombart, 2008).
Se estimó el flujo de genes contemporáneo (en las últimas generaciones) entre las poblaciones de M. incilis muestreadas en diferentes puntos del Caribe de Colombia, utilizando el método bayesiano implementado en BAYESASS v3 (Wilson & Rannala, 2003). Inicialmente ajustamos los valores de delta para las frecuencias alélicas, las tasas de migración y los coeficientes de endogamia para asegurar que las tasas de aceptación de los cambios en estos parámetros estuvieran entre 40 y 60 % (Wilson & Rannala, 2003). Realizamos 10 ejecuciones repetidas con 1 x 107 iteraciones MCMC y una fase de burn-in de 106 iteraciones. Cada ejecución inició con una semilla de generador de números aleatorios diferente para verificar la convergencia según la concordancia entre las ejecuciones en las estimaciones de los parámetros medios posteriores. La convergencia del modelo se evaluó mediante la comparación de las densidades de probabilidad posterior de los coeficientes de endogamia y las frecuencias de los alelos en 10 series repetidas (Faubet, Waples, & Gaggiotti, 2007). La medida de la de varianza bayesiana se utilizó para determinar la estimación que mostraba el mejor ajuste del modelo (Spiegelhalter, Best, Carlin, & van der Linde, 2002) de cada una de las réplicas de cadenas de MCMC utilizando las recomendación y script de R de Meirmans (2014). Para esta ejecución de mejor ajuste se volvió a ejecutar el análisis utilizando la semilla de la ejecución que mejor se ajusta, pero aumentamos la duración de la ejecución a 5 × 107 iteraciones. Los resultados presentados corresponden a la ejecución final.
Resultados
Diversidad genética: Un total de 94 alelos con promedio de 11.6 alelos por locus en los sietes loci microsatélites examinados, los cuales presentaron polimorfismo en cada locus. El promedio de alelos por locus en cada localidad no tuvo mucha variación (Tabla 2). Los alelos más frecuentes en el locus Mcs-15 AM fueron el alelo 266 para las localidades Santuario de Fauna y Flora Los Flamencos (SFFF), Santuario de Fauna y Flora Ciénaga Grande de Santa Marta (CGSM), Ciénaga La Virgen (CVIR) y Santuario de Fauna y Flora El Corchal (SFFC) y el alelo 214 en el locus Muce-80 para las localidades de SFFF, CGSM, SFFC y Ciénaga La Caimanera (CCAI). Se encontraron 23 alelos privados, seis de ellos pertenecen a SFFC, los demás alelos están distribuidos en cada una de las localidades, siendo la CGSM el sitio con menor número de alelos privados (un alelo).
Todas las localidades mostraron desviación en el equilibrio de Hardy-Weinberg, con una probabilidad significativa en cada localidad (P < 0.001) (Tabla 2). Todos los loci microsatélites mostraron evidencia de alelos nulos y alelos “stutter” pero no de perdida alélica; lo que sugiere que puede existir un exceso de genotipos homocigotos. Los valores promedio de Ho mostraron una baja variabilidad genética en todas las localidades, con datos que oscilaron entre 0.108 para la localidad de CCAI y 0.402 en SFFF. Con respecto a la heterocigosidad esperada, los valores oscilaron entre 0.772 en CVIR y 0.868 en Bahía de Cispatá (BCIS). La localidad más crítica en términos de variabilidad genética fue CCAI debido a los valores de HO (0.108) y FIS (0.873).
Se determinaron ocho asociaciones no aleatorias (desequilibrio de ligamiento) en las 123 comparaciones realizadas que involucraron pares diferentes de loci para cada caso, por lo que se asumió la independencia de los datos para la realización de cada uno de los análisis.
La prueba de rangos de Wilcoxon mostró dos tallos para excesos (P < 0.05) o deficiencia de individuos heterocigotos (P < 0.001), siendo significativo en cuatro de las localidades estudiadas bajo el modelo de mutación IAM y demostrando que han pasado por un reciente cuello de botella debido a una disminución en el tamaño de la población. Bajo el modelo TPM solamente BCIS ha pasado por cuello de botella y según el modelo SMM en tres localidades se evidenció cuello de botella en un pasado reciente.
Locus | Mcs15AM | Mce2 | Mce7 | Mce11 | Mce3 | Muce55 | Muce80 | Promedio |
SFFF | ||||||||
N | 50 | 50 | 50 | 50 | 50 | 50 | 49 | 50 |
Na | 14 | 21 | 12 | 5 | 13 | 7 | 11 | 12 |
RA | 22.67 | 22.1 | 21.19 | 5.52 | 18.24 | 10.33 | 19.22 | 17.04 |
Ho | 0.58 | 0.22 | 0.52 | 0.00 | 0.54 | 0.3 | 0.653 | 0.402 |
He | 0.907 | 0.942 | 0.897 | 0.663 | 0.814 | 0.685 | 0.854 | 0.823 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.362 | 0.768 | 0.422 | 1 | 0.3385 | 0.564 | 0.236 | 0.528 |
CGSM | ||||||||
N | 40 | 40 | 40 | 40 | 40 | 40 | 40 | 40 |
Na | 14 | 17 | 13 | 5 | 1 | 9 | 11 | 12 |
RA | 20.41 | 20.77 | 16.45 | 5.54 | 16.21 | 8.25 | 19.87 | 15.36 |
Ho | 0.5 | 0.3 | 0.55 | 0.00 | 0.325 | 0.125 | 0.55 | 0.33 |
He | 0.901 | 0.921 | 0.814 | 0.546 | 0.863 | 0.625 | 0.858 | 0.790 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.448 | 0.676 | 0.327 | 1 | 0.626 | 0.801 | 0.361 | 0.606 |
CVIR | ||||||||
N | 39 | 39 | 39 | 39 | 39 | 39 | 39 | 39 |
Na | 14 | 20 | 11 | 5 | 10 | 6 | 17 | 12 |
RA | 19.59 | 20.97 | 15.05 | 5.14 | 14.79 | 8.63 | 22.5 | 15.24 |
Ho | 0.538 | 0.231 | 0.333 | 0.00 | 0.41 | 0.103 | 0.692 | 0.330 |
He | 0.891 | 0.927 | 0.802 | 0.393 | 0.737 | 0.743 | 0.909 | 0.772 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.398 | 0.753 | 0.58 | 1 | 0.446 | 0.863 | 0.240 | 0.613 |
SFFC | ||||||||
N | 40 | 40 | 40 | 40 | 40 | 40 | 40 | 40 |
Na | 13 | 18 | 13 | 5 | 11 | 6 | 12 | 11 |
RA | 17.61 | 19.68 | 18.41 | 5.44 | 17.59 | 8.84 | 17.81 | 15.05 |
Ho | 0.55 | 0.15 | 0.45 | 0.025 | 0.5 | 0.1 | 0.6 | 0.339 |
He | 0.859 | 0.910 | 0.863 | 0.438 | 0.828 | 0.753 | 0.859 | 0.787 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.3628 | 0.836 | 0.481 | 0.943 | 0.399 | 0.868 | 0.303 | 0.600 |
CCAI | ||||||||
N | 40 | 37 | 40 | 39 | 40 | 37 | 26 | 40 |
Na | 13 | 13 | 10 | 9 | 6 | 7 | 14 | 10 |
RA | 13.96 | 15.5 | 12.95 | 9.93 | 6.76 | 8.06 | 15.02 | 11.74 |
Ho | 0.125 | 0.108 | 0.125 | 0.128 | 0.025 | 0.054 | 0.192 | 0.108 |
He | 0.851 | 0.843 | 0.847 | 0.808 | 0.677 | 0.75 | 0.931 | 0.815 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.8547 | 0.8732 | 0.854 | 0.843 | 0.963 | 0.928 | 0.796 | 0.873 |
BCIS | ||||||||
N | 40 | 40 | 40 | 40 | 40 | 40 | 39 | 40 |
Na | 14 | 20 | 14 | 11 | 10 | 9 | 14 | 13 |
RA | 20.31 | 19.9 | 19.06 | 12.03 | 14.42 | 10.77 | 20.12 | 16.66 |
Ho | 0.375 | 0.25 | 0.425 | 0.125 | 0.35 | 0.1 | 0.462 | 0.298 |
He | 0.914 | 0.933 | 0.87 | 0.816 | 0.819 | 0.822 | 0.902 | 0.868 |
H-WE | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
Fis | 0.592 | 0.734 | 0.514 | 0.848 | 0.575 | 0.879 | 0.491 | 0.593 |
Número de muestras amplificadas (N); número de alelos (Na); Riqueza Alelica (RA); heterocigosidad esperada (He); heterocigosidad observada (Ho); desviación del equilibrio Hardy-Weinberg (H-WE) y coeficiente de endogamia (Fis). El nivel de significancia fue de P < 0.001.
Number of amplified samples (N); number of alleles (Na); allelic richness (AR); expected heterozygosity (He); observed heterozygosity (Ho) deviation from the Hardy-Weinberg equilibrium (H-WE) and inbreeding coefficient (Fis). The level of significance was P < 0.001.
Divergencia genética: Los resultados de la AMOVA mostraron diferenciación genética entre todas las localidades (FST(5, 248) = 0.106, P = 0.001), siendo el mayor nivel de diferenciación genética cuando las localidades del Golfo de Morrosquillo (CCAI y BCIS) se contrastaron con las demás (SFFF, CGSM, CVIR y SFFC; Tabla 3). El análisis de comparación múltiple permitió observar que los valores de FST entre las seis localidades variaron de 0.061 entre CGSM y SFFC hasta 0.140 entre CVIR y CCAI (Tabla 4), presentando una moderada diferenciación genética de acuerdo con la escala propuesta por Wright (Hartl & Clark, 1997).
El análisis de la prueba de Mantel no mostró evidencia significativa de una relación directa entre la distancia genética (Fst pareadas) y la distancia geográfica (Km) entre las seis localidades muestreadas (Z = 932.4847, r = 0.0101, P = 0.2920).
Variación | Agrupamiento | ||
1. Todas las localidades | 2. Noreste del río Magdalena vs. Suroeste | 3. Golfo de Morrosquillo vs. Noreste del golfo | |
Entre grupos (FCT) | N/A | 0.004 P= 0.3 | 0.018 P= 0.06 |
Entre poblaciones dentro de grupos (FSC) | N/A | 0.104 P= 0.001 | 0.097 P= 0.001 |
Dentro de poblaciones (FST) | 0.106 P= 0.001 | 0.004 P= 0.001 | 0.018 P= 0.001 |
Ver sección métodos los agrupamientos. See methods for grouping.
Localidades | SFFF | CGSM | CVIR | SFFC | CCAI | BCIS |
SFFF | 0 | 0.049 | 0.078 | 0.001 | 0.000 | 0.002 |
CGSM | 0.084* | 0 | 0.034 | 0.076 | 0.001 | 0.029 |
CVIR | 0.129* | 0.159* | 0 | 0.006 | 0.000 | 0.008 |
SFFC | 0.094* | 0.061* | 0.126* | 0 | 0.000 | 0.000 |
CCAI | 0.108* | 0.115* | 0.140* | 0.123* | 0 | 0.014 |
BCIS | 0.082* | 0.097* | 0.107* | 0.115* | 0.063* | 0 |
(*) indica comparaciones significantes (P < 0.05).
(*) indicates significant comparisons (P < 0.05).
El dendrograma construido a partir de la matriz de distancia genética (δµ2) permitió distinguir dos subpoblaciones para M. inclis en el Caribe de Colombia. Una gran población conformada por la relación genética entre las localidades del SFFF, CGSM, CVIR y SFFC; una segunda por BCIS y CCAI (Fig. 2A). Esta agrupación fue similar a lo observado en el procedimiento realizado con STRUCTURE (Evanno et al., 2005), ya que se identificaron dos grupos poblacionales (K = 2, Fig. 2B). El primero está conformado por los individuos de las localidades de SFFF, CGSM, CVIR y SFFC (rojo) y un segundo por CCAI y BCIS (verde), siendo esta última localidad con la mayor mezcla de genotipos del primer grupo (Fig. 2C).
El análisis de DAPC fue inicialmente realizado sin algún grupo asignado a priori. Para obtener el número óptimo de cluster con la función find.cluster, 75 ejes representaron el 80 % del total de varianza retenida. El programa cubrió un ámbito de posibles grupos de 1-25. El valor bajo de BIC (96.6117) correspondió a K = 6. Para el análisis de DAPC, 75 ejes de PCA y cinco funciones discriminantes fueron retenidos por contener el 95 % de la varianza. Se formaron tres agrupaciones, una gran agrupación conformada por los sitios SFFF, CGSM, CVIR y SFFC, mientras que CCAI y BCIS se diferenciaron como dos poblaciones independientes (Fig. 3).
Flujo genético: Las estimaciones bayesianas de flujo genético contemporáneo entre las seis poblaciones (expresado como proporciones de la población que migra a otras regiones en las últimas generaciones) indicó que hay una fracción de individuos de la localidad BSCI que son migrantes de la población CCAI (0.136) (Tabla 5). Para el caso de las localidades del norte se observó que las tasas de flujo de genes fueron muy bajas (< 0.01).
SFFF | CGSM | CVIR | SFFC | CCAI | BCIS | |
SFFF | 0.960 | 0.007 | 0.010 | 0.006 | 0.006 | 0.011 |
CGSM | 0.009 | 0.948 | 0.014 | 0.008 | 0.008 | 0.013 |
CVIR | 0.009 | 0.013 | 0.943 | 0.019 | 0.008 | 0.008 |
SFFC | 0.015 | 0.013 | 0.008 | 0.949 | 0.007 | 0.008 |
CCAI | 0.008 | 0.015 | 0.008 | 0.008 | 0.945 | 0.017 |
BCIS | 0.018 | 0.008 | 0.008 | 0.008 | 0.136 | 0.822 |
Poblaciones de individuos que migran (filas) y población de origen de los migrantes (columnas). Los valores a lo largo de la diagonal son las proporciones de individuos derivados de las poblaciones originales en cada generación. Las tasas de migración ≥ 0.10 están en cursiva. Las desviaciones estándar para todas las distribuciones fueron < 0.05.
Populations of individuals that migrates (rows) and population origin of migrants (columns). The values along the diagonal are the proportions of individuals derived from the original populations in each generation. Migration rates ≥ 0.10 are in italics. The standard deviations for all distributions were < 0.05.
Discusión
Se presenta el primer estudio genético de Mugil incilis en el Caribe colombiano, empleando partidores no específicos para la especie. La amplificación cruzada usando partidores diseñados para M. cephalus permitió tener información relevante sobre la variabilidad y estructura genética de esta especie. Todos los loci fueron polimórficos y los resultados obtenidos son comparables con estudios realizados para otras especies de mugílidos (Miggiano et al., 2005; Pacheco-Almanzar et al., 2017; Shen et al., 2010; Whitfield, Panfili, & Durand, 2012; Xu et al., 2010). La riqueza alélica observada en M. incilis difirió de la documentada para poblaciones de M. cephalus y M. liza (Mai et al., 2014; Miggiano et al., 2005; Whitfield et al., 2012). Por ejemplo, algunos microsatélites presentaron menos alelos (Mcs15AM, Mce3 y Mce7) que otros (Mce2 y Mce11), en comparación con las poblaciones de M. cephalus (Miggiano et al., 2005; Shen et al., 2010). Sin embargo los bajos valores de alelos se han reportado en especies hermanas como Mugil liza (Mai et al., 2014) y Mugil curema (Pacheco-Almanzar et al., 2017) que han utilizado los microsatélites diseñados para M. cephalus.
Con respecto a los valores promedios de heterocigosidad observada y esperada para M. incilis, se obtuvieron valores menores a los documentados para poblaciones de M. cephalus (HO = 0.2593-0.8966, He = 0.3047-0.8454) (Miggiano et al., 2005; Shen et al., 2010; Xu et al., 2010). Por lo anterior, fue evidente un déficit de heterocigotos asociado a la alta endogamia poblacional (Fis > 0.5) debido probablemente a la reducción drástica de los tamaños efectivos poblacionales ocasionada por la sobrepesca y la contaminación ambiental. Los análisis de cuello de botella indicaron que todas las localidades experimentaron este fenómeno en el pasado reciente y coincide con la disminución de tallas y de las capturas desembarcadas que experimentan los peces de M. incilis en la mayoría de las lagunas costeras del Caribe de Colombia (Blanco-Racedo, 1983; Sánchez Ramírez et al., 1998; Narváez-Barandica et al., 2008; INVEMAR, 2014; Martínez et al., 2014). Así mismo, la desviación del equilibrio Hardy-Weinberg en todos los loci puede también explicar el déficit de heterocigotos por la acción de algún tipo de selección y de la deriva genética (Hartl & Clark, 1997). Por otro lado, se detectó la presencia de alelos nulos en todos los loci y esto pudo afectar las estimaciones de diversidad genética debido a los comunes errores de amplificación/genotipificación de los loci microsatélites y no por alguna señal de diferenciación genética poblacional (Hoffman & Amos, 2005; Landínez-García, Ospina-Guerrero, Rodríguez-Castro, Arango-Isaza, & Márquez-Fernández, 2009; Orozco-Berdugo & Narváez-Barandica, 2014; Meirmans, 2015). A pesar de lo anterior, siempre será difícil probar las reales causas que afectan la diversidad genética de las poblaciones; sin embargo, este estudio genera una primera aproximación del estado genético de las poblaciones muestreadas de M. incilis en el Caribe de Colombia.
Para M. incilis se determinó una diferenciación genética poblacional entre todas las localidades estudiadas. La mayor diferenciación genética la presentaron BCIS y CCAI con el resto de las localidades. Es así que esta diferenciación genética coincidió en el análisis de dendrograma y en el bayesiano de STRUCTURE, donde los peces muestreados en SFFF, CGSM, CVIR y SFFC conforman la subpoblación 1; mientras que los de CCAI y BCIS conforman la subpoblación 2. Con respecto a esta última subpoblación, las dos localidades que la conforman presentan la mayor relación genética revelado a partir del análisis de las tasas de migración (Tabla 5). Sin embargo, fue evidente que algunos genotipos de las localidades de la subpoblación 1 se mezclaran con los de la subpoblación 2, principalmente en la localidad BCIS (Fig. 2A, Fig. 2C). A diferencia de lo anterior, el análisis de DAPC determinó que los peces muestreados en SFFF, CGSM, CVIR y SFFC también conforman la subpoblación 1; pero esta vez aparece BCIS conformando una subpoblación 2 y CCAI una subpoblación 3. La principal explicación se debe a que los peces recolectados en BCIS tienen similar proporción de genotipos de la subpoblación 1 y 3 (Fig. 2C). Por este motivo, y con principio conservacionista se decide por optar la sub-estructuración genética de tres subpoblaciones para M. incilis en el Caribe de Colombia. De igual manera es importante resaltar que las actividades pesqueras y las condiciones ambientales son diferentes en estas lagunas (INVEMAR, 2014). Por otro lado, el análisis de autocorrelación espacial no evidenció la existencia del modelo poblacional de aislamiento por distancia en la población de M. incilis asociada al Caribe de Colombia.
Con respecto a otros estudios, la subestructuración genética de M. incilis propuesta aquí es similar a lo registrado en América para otras especies del mismo género (Mai et al., 2014; Miggiano et al., 2005; Siccha-Ramirez, Menezes, Nirchio, Foresti, & Oliveira, 2014; Whitfield et al., 2012). Es así como Mai et al. (2014) documentaron para M. liza la existencia de dos poblaciones diferenciadas genéticamente a lo largo de las costas de Brasil y Argentina, los cuales son separados por la dinámica oceanográfica y ambiental de las corrientes de Brasil y Malvinas que confluyen en el área. Otro ejemplo lo presentaron Pacheco-Almanzar et al. (2017), que describieron tres poblaciones entre el Golfo de México y una localidad ubicada en el Noroeste del Atlántico, las cuales son definidas por la variación espacial de las condiciones oceanográficas del Golfo, diferentes tiempos de desoves y diferencias en las condiciones ecológicas y geomorfológicas de las lagunas.
Lo que concierne a M. incilis en el Caribe de Colombia, en las localidades estudiadas en esta primera aproximación genética el patrón de diferenciación genética observado probablemente está influido por el comportamiento de las corrientes que predominan durante la época de migración reproductiva desde las lagunas hacia el mar (septiembre-diciembre) (Bustos-Montes et al., 2012; Mármol, Viloria, & Blanco, 2010). Estas corrientes ocasionan un punto de divergencia entre ambas poblaciones puesto que una porción de la corriente Panamá-Colombia, denominada contracorriente del Darién, viaja adyacente a la costa en dirección noreste y puede estar limitando el flujo genético desde la subpoblación 1 (noreste) hacia la 2 (interior del Golfo de Morrosquillo) (Centurioni & Niiler, 2003; Sheng & Tang, 2003; Gyory, Mariano, & Ryan, 2008; Miloslavich et al., 2010). Sin embargo, durante los máximos valores de los vientos Alisios del NE (enero-marzo), la contracorriente del Darién se debilita en las capas de aguas superficiales y permite una mayor influencia de la corriente Caribe para arrastrar larvas de la subpoblación 1 hacia el sur del Caribe de Colombia (Centurioni & Niiler, 2003; Gyory et al., 2008). Esto podría explicar por qué los peces de SFFF, CGSM, CVIR y SFFC no se diferencian genéticamente en el análisis bayesiano; y por qué una considerable proporción de genotipos de la subpoblación 1 están presentes en BCIS y en una mínima fracción en CCAI (Fig. 2C). Por otro lado, la diferenciación genética de los peces de CCAI con el resto de localidades se debe a que las condiciones locales de vientos, olas y corrientes dentro del Golfo de Morrosquillo (Ricaurte Villota, Rangel Buitrago, Galeano, & Coca Domínguez, 2012) probablemente influye en el autoreclutamiento biológico local y también limite la dispersión de las larvas hacia la subpoblación 1.
En conclusión, esta investigación representa el primer estudio genético poblacional de M. incilis en el Caribe colombiano. Los resultados de este estudio mostraron indicadores de variabilidad genética bajos en comparación con las especies del mismo género y especialmente en las localidades de los tres santuarios de conservación para M. incilis. Esta condición demuestra que las áreas no son suficientes para la conservación de M. incilis, puesto que no abarcan todo el espacio geográfico que necesita la especie para cumplir su ciclo de vida (Blanco-Racedo, 1983; Sánchez Ramírez et al., 1998). Se sabe que M. incilis es un recurso muy vulnerable a la actividad pesquera, principalmente cuando los peces salen desde las áreas de alimentación y desarrollo hacia el mar para realizar sus migraciones reproductivas (Narváez-Barandica et al., 2008).
La estructura genética de M. incilis presentó tres subpoblaciones en el Caribe de Colombia. La primera tiene mayor ocupación en un ámbito geográfico que va desde el departamento de La Guajira (SFFF) hasta el de Sucre (SFFC). Las otras dos están ubicadas en las lagunas costeras dentro del Golfo de Morrosquillo; una segunda subpoblación en BCIS y una tercera en CCAI. Esto se pudo explicar por el comportamiento migratorio del pez y la dinámica oceanográfica de la corriente Caribe, de la contracorriente del Darién y a las condiciones oceanográficas del interior del Golfo de Morrosquillo. Se recomienda en futuros estudios generar iniciadores específicos para la especie que permitan comprobar los resultados obtenidos en este y de igual manera ampliar el número de sitios de muestreo.
La variabilidad genética baja y la subdivisión poblacional de M. incilis debe motivar a las autoridades ambientales y a las comunidades de pescadores adelantar procesos para su conservación y manejo pesquero en todo el Caribe de Colombia. Se deben hacer esfuerzos interinstitucionales para evitar la pérdida de la variabilidad genética que experimenta la especie debido a la sobreexplotación pesquera y a la mala calidad ambiental de las aguas donde habita. Para iniciar con estos procesos se proponen las siguientes recomendaciones: i). Establecer una veda temporal que coincida con el periodo de migración reproductiva de la especie. ii). Acoger tres unidades de manejo (UM) con el objetivo de establecer las estrategias de manejo y conservación para M. incilis de acuerdo a la historia de vida y a la forma como se aprovecha con la pesca en cada unidad (Cabral, Mamauag, & Aliño, 2015; Moritz, 1994; Roberts et al., 2003). La primera UM está ubicada en el área que cubre la subpoblación 1 (SFFF, CGSM, CVIR y SFFC); la segunda está en BCIS y la tercera está en CCAI, donde se deben desarrollar estudios biológicos-pesqueros que permitan establecer puntos de referencias para futuros programas de manejo. iii). Finalmente, generar espacios de concertación entre las comunidades de pescadores, las instituciones administradoras de los recursos pesqueros y naturales para acordar medidas de manejo pesquero más acertadas para el aprovechamiento sostenible de M. incilis.