Genética poblacional de camarón blanco Litopenaeus vannamei en Sinaloa, México
Ricardo Perez-Enriquez a, *, Ana María Millán-Márquez b, Pedro Cruz-Hernández a, y César Julio Saucedo-Barrón c
a Centro de Investigaciones Biológicas del Noroeste, Calle IPN 195, 23096 La Paz, Baja California Sur, México
b Universidad de Valle, # 100-00, 760032 Cali, Colombia
c Instituto Sinaloense de Acuacultura y Pesca, Fray Servando Teresa de Mier 1870, 80129 Culiacán, Sinaloa, México
*Autor para correspondencia: rperez@cibnor.mx (R. Perez-Enriquez)
Resumen
El camarón blanco Litopenaeus vannamei es una especie de importancia pesquera en Sinaloa, México cuya producción ha declinado en años recientes. Para su recuperación, la repoblación mediante la liberación de larva cultivada se ha considerado como alternativa. Esta estrategia requiere incluir información sobre las características genéticas de las poblaciones silvestres, para garantizar que tales programas no impacten negativamente sobre su eficacia biológica en el largo plazo. En este trabajo se estudió la estructura genética del camarón blanco en Sinaloa, México, mediante los perfiles genéticos de 6 microsatélites en muestras de 10 localidades ubicadas en 5 lagunas costeras. Se registraron altos niveles de diversidad genética (15.6 alelos por locus y 0.87 heterocigosidad esperada), así como desviaciones del equilibrio de Hardy-Weinberg (EHW) en la mayor parte de los sitios. A pesar de observarse diferencias genéticas significativas entre localidades (Fst = 0.012; p < 0.001), la falta de un patrón geográfico específico y una baja diferenciación entre estuarios (Fst = 0.005; p = 0.062), indicó homogeneidad genética regional. Las desviaciones del EWH y las diferencias entre localidades se explican por el efecto Wahlund. Asimismo, se discuten recomendaciones para el manejo de las actividades de repoblación sobre la base de criterios genéticos.
Palabras clave:
Diversidad genética; Estructura genética; Microsatélites; Repoblación; Flujo genético
© 2018 Universidad Nacional Autónoma de México, Instituto de Biología. Este es un artículo Open Access bajo la licencia CC BY-NC-ND
Population genetics of whiteleg shrimp Litopenaeus vannamei in Sinaloa, Mexico
Abstract
Whiteleg shrimp Litopenaeus vannamei is an important fishery species in Sinaloa, Mexico, which production has declined in recent years. For its recovery, the release of cultivated larvae has been suggested as an alternative. This strategy requires to include information about the genetic characteristics of the wild populations to warrant that these programs do not negatively impact their long-term biological fitness. In this work, the genetic structure of Whiteleg shrimp in Sinaloa, Mexico was studied through the genetic profile of 6 microsatellites of samples from 10 sites located in 5 coastal lagoons. High levels of genetic diversity (15.6 number of alleles per locus and 0.87 expected heterozygosity), and deviation from the Hardy-Weinberg equilibrium (HWE) in most sites were registered. In spite of the significant genetic differences observed among locations (Fst = 0.012; p < 0.001), the lack of a specific geographical pattern and a low differentiation among estuaries (Fst = 0.005; p = 0.062) indicated regional genetic homogeneity. The departures from HWE and the differences among locations are explained by the Wahlund effect. Also, recommendations for the management of restocking activities under genetic criteria are discussed.
Keywords:
Genetic diversity; Genetic structure; Microsatellites; Restocking; Gene flow
La liberación de larva o juveniles producidos en cautiverio para repoblación, ya sea para la recuperación del tamaño de la población desovante o para el incremento del tamaño poblacional al nivel permisible por la capacidad de carga del ecosistema para recuperar los niveles de captura, es una de las medidas con potencial (Bell et al., 2006). Sin embargo, persiste una preocupación inherente al potencial impacto negativo sobre las poblaciones silvestres, particularmente en términos genéticos (Perez-Enriquez et al., 2013). Miller y Kapuscinsky (2003) señalan como potenciales impactos genéticos la reducción de la variabilidad genética en la población silvestre a nivel intrapoblacional y diferenciación interpoblacional, entre otros.
La pérdida de variabilidad genética que se da cuando poblaciones silvestres son suplementadas con poblaciones de cultivo, ocurre debido al efecto “Ryman-Laikre” (ver Miller y Kapuscinsky, 2003) en el cual hay una reducción del tamaño efectivo de la población total (Ne), debido a una proporción desigual entre la cantidad de individuos reclutas respecto al número efectivo de reproductores de cada grupo (silvestre y de cultivo). Este fenómeno sucede cuando la larva producida proviene de un número reducido de padres, situación que ocurre comúnmente en especies en cultivo (Taniguchi, 2004). La reducción en Ne puede generar un incremento de la deriva génica y/o de la endogamia, y por consecuencia, una mayor probabilidad de depresión por endogamia; esto a su vez, tiene el potencial de afectar negativamente la eficacia biológica de la población, así como su viabilidad y su potencial evolutivo (Hallerman, 2003).
La introducción de larvas de un mismo origen, en subpoblaciones genéticamente diferenciadas, tiene el potencial de causar disrupciones en la conformación de complejos de genes coadaptados que se expresan ante ambientes en particular, dando como resultado una menor eficacia biológica de dichas subpoblaciones (Miller y Kapuscinsky, 2003). Bajo estas consideraciones, la instrumentación de un programa de repoblación requiere, como punto de partida, la caracterización genética de las poblaciones objetivo, en términos de su diversidad genética, nivel de estructuración en subpoblaciones y flujo genético, entre otros parámetros.
En el caso particular de L. vannamei, la caracterización genética de las poblaciones silvestres en un amplio rango geográfico (miles de kilómetros) desde México hasta Ecuador ha mostrado resultados contrastantes (Sunden y Davis, 1991; Valles-Jiménez et al., 2004), probablemente debido a la naturaleza de los marcadores genéticos. Sunden y Davis (1991) encontraron bajos niveles de diversidad y baja diferenciación genética interpoblacional con aloenzimas, mientras que Valles-Jiménez et al. (2004) utilizando microsatélites reportaron una mayor diversidad genética y diferencias significativas entre las regiones de estudio.
A nivel regional (decenas a cientos de kilómetros) no existen estudios que determinen el grado de diferenciación genética en L. vannamei entre localidades del golfo de California. De acuerdo con un período larvario planctónico de aproximadamente 3 semanas en mar abierto (Hendrickx, 2001) y corrientes estacionales intensas en la parte sureste del golfo capaces de transportar partículas por largas distancias en períodos cortos de tiempo (Marinone et al., 2011), la hipótesis es que las poblaciones de Sinaloa no presentan diferenciación genética significativa. Sin embargo, en otras especies de peneidos (e.g., Penaeus (Litopenaeus) stylirostris, F. californiensis) varios autores indican evidencias de diferenciación genética dentro del propio golfo (Aubert y Lightner, 2000; De la Rosa-Vélez et al., 2000; Ramos-Paredes y Grijalva-Chon, 2003).
El objetivo del presente estudio fue caracterizar la composición genética de la población del camarón blanco (L. vannamei) de Sinaloa, México, determinar si existe diferenciación genética entre las localidades de estudio y definir recomendaciones de manejo, incluyendo la información genética, en función de potenciales actividades de repoblación de la especie.
Materiales y métodos
Se realizaron extracciones de ADN mediante el kit comercial de extracción y purificación de ADN (Wizard® Genomic DNA Purification-Promega). Se analizaron 6 microsatélites (Pvan1758, Pvan1815, TumxLv8.256, TumxLv10.312, TumxLv9.103 y Lv5) marcados con flourescencia (6-FAM, VIC, NED o PET), siguiendo el protocolo de amplificación reportado por Perez-Enriquez et al. (2009). La genotipificación se realizó mediante el analizador genético ABI3130 (Applied Biosystems), los alelos fueron leídos y los tamaños clasificados con el programa GeneMapper® version 4.0, homogenizando las lecturas entre localidades con el estándar de tamaño GeneScan 500 LIZ (Applied Biosystems). Las lecturas de alelos se confirmó al menos una vez en aquellas muestras que presentaron patrones de genotipos no esperados. Se calcularon las frecuencias alélicas y se estimó la diversidad genética [número de alelos por locus (na), número efectivo de alelos por locus (ne), heterocigosidad observada (Ho) y heterocigosidad esperada (He)] con el programa Arlequin 3.5 (Excoffier et al., 2005). Con el mismo programa se estimó la desviación del equilibrio de Hardy-Weinberg (EHW) para cada locus y cada población, y se verificó la presencia de alelos nulos mediante el programa Microchecker (Van Oosterhout et al., 2004).
La estructura genética de la población de Sinaloa se probó mediante el valor de Fst obtenido mediante un análisis de varianza molecular (Amova) utilizando Arlequin 3.5. Se realizó un análisis de Fst pareado por localidades usando el paquete Adegenet de R (Jombart, 2008). Los valores de Fst se visualizaron como distancias genéticas en un árbol “vecino más cercano” de Saitou y Nei (1987) estimado con el paquete Ape de R (Paradis et al., 2004). Para visualizar la distribución de genotipos por individuo/localidad en 2 dimensiones se corrió un análisis de componentes principales (PCoA) utilizando el paquete
Resultados
Se obtuvo un valor significativo de diferenciación genética global entre las 10 localidades (Fst = 0.012, p < 0.001). El análisis de Amova dio diferencias significativas entre las lagunas costeras (Fst = 0.004; p < 0.05) y no arrojó diferencias significativas entre los sitios dentro de cada laguna costera (bahía de Jitzamuri, bahía del Colorado, bahía de Navachiste, bahía Santa María y ensenada del Pabellón-Altata; fig. 1), (Fst = 0.005; p = 0.196). El análisis pareado de Fst indicó que 11 de las 45 comparaciones simultáneas fueron estadísticamente significativas (fig. 3). Sin embargo, ninguno de los análisis de diferenciación [i.e. distancias genéticas (fig. 4a), PCoA (fig. 4b), método de grupos (fig. 4c)] mostraron una agrupación acorde a la ubicación geográfica de las localidades y las lagunas costeras. Un segundo análisis de Amova agrupando los datos por laguna costera (bahía de Jitzamuri, bahía del Colorado, Bahía de Navachiste, bahía Santa María y ensenada del Pabellón-Altata; fig. 1), reveló que no hubo diferencias genéticas significativas entre estos cuerpos de agua (Fst = 0.005; p = 0.062).
Discusión
El análisis de estructura genética de la población silvestre mostró un patrón poco consistente, ya que aun cuando el índice Fst indicó diferencias genéticas significativas entre localidades, dichas diferencias no fueron acordes con su ubicación geográfica. Sin embargo, dado que el análisis en donde se agruparon muestras de las 2 localidades de cada sistema estuarino como una sola localidad no indicó valores pareados de Fst significativos, se concluye que el camarón que habita las 5 lagunas costeras del estudio constituyen una sola población panmíctica, con flujo genético entre ellas.
Si bien, para el camarón blanco Litopenaeus vannamei no hay estudios de diferenciación genética de corta escala geográfica, los análisis a amplia escala indican niveles de diferenciación significativos, derivados de barreras geográficas u oceanográficas (Valles-Jiménez et al., 2004). Supungul et al. (2000) con Penaeus monodon en Tailandia obtienen resultados similares con lo inferido para L. vannamei, con una población genéticamente diferenciada respecto a su área de distribución, pero homogeneidad local.
La homogeneidad genética de L. vannamei a nivel regional pudiera ser explicada por la dispersión larvaria de organismos entre lagunas costeras debido a las corrientes y a la carencia de barreras al flujo. La combinación de un periodo larvario de aproximadamente 25 días (FAO, 2006) con los patrones de corrientes de la parte sur del golfo de California (la cual contiene al área de estudio), favorecen el transporte de partículas en una escala de cientos de kilómetros (Marinone et al., 2011), y soportarían la hipótesis de la falta de diferenciación genética regional.
Sin embargo, los resultados contrastan con lo reportado para otras especies de peneidos dentro de la región de estudio con periodos larvarios de duración similar (e.g., L. stylirostris y F. californiensis), los cuales indican flujo genético limitado y, por ende, diferenciación genética en
De cualquier manera, no se puede descartar que exista un fenómeno de corta escala asociado a la diferenciación genética, tal como el aislamiento y la presencia de remolinos en el Alto golfo (Marinone et al., 2011), que sirvan como barreras oceanográficas y/o geográficas. Esto puede ser similar a lo reportado en L. schmitii por Borell et al. (2007), peneido del Atlántico con un periodo larvario de duración similar a L. vannamei (Boschi, 1977), quienes reportaron diferenciación genética en un rango geográfico
Por otra parte, una posible alternativa para explicar la desviación del EHW aun habiendo panmixia entre localidades, es el efecto Wahlund, el cual establece que la mezcla de poblaciones con diferentes frecuencias genotípicas tiende a presentar un déficit de heterocigotos (ver Hedrick, 2000). En este sentido, considerando que el cultivo de camarón en el estado de Sinaloa es una actividad de alta relevancia (Conapesca, 2013), se propone la hipótesis de la mezcla de poblaciones de cultivo con silvestres. Esta hipótesis podría ser soportada por el hecho de que la composición genética de las poblaciones de cultivo es diferente a las silvestres (Mendoza-Cano et al., 2013; Perez-Enriquez et al. 2009), y a que se ha observado la presencia de ejemplares con perfiles genéticos de camarón cultivado en muestras de poblaciones silvestres del estado de Sinaloa (Perez-Enriquez et al., 2018).
Las implicaciones de los resultados del presente estudio para efectos de programas potenciales de repoblación de la especie van en 3 vertientes: 1) dado que no hay elementos que indiquen que cada localidad deba ser tratada, desde el punto de vista genético, de manera independiente, la constitución de un lote de reproductores de origen silvestre puede tener como origen cualquiera de los cuerpos de agua comprendidos en el estudio; 2) dado el potencial de disrupción de la composición genética derivada de las actividades acuícolas, se sugiere que se evite la liberación de larvas proveniente de laboratorios dedicados a la producción comercial de larva; 3) el manejo por repoblación requiere de manera adicional la consideración de factores ecológicos (e.g., razones que explican la disminución de las poblaciones, sitios adecuados para la liberación, calidad del hábitat, etc.) y pesqueros (distribución del esfuerzo, cuotas de captura, etc.), entre otros (ver Miller y Kapuscinsky, 2003; Perez-Enriquez et al., 2013) para incrementar la probabilidad de éxito.
Agradecimientos
Referencias
Bell, J. D., Bartley, D. M., Lorenzen, K. y Loneragan, N. R. (2006). Restocking and stock enhancement of coastal fisheries: Potential, problems and progress. Fisheries Research, 80, 1–8.
Borrell, Y. J., Arenal, F., Mirelles-Mbemba, Z., Santana, O., Díaz-Férnandez, R., Vázquez, E. et al. (2007). Spatial and temporal genetic analysis of the Cuban white shrimp Penaeus (Litopenaeus) schmitti. Aquaculture, 272 (Suplem.), S125–S138.
Boschi, E. (1977) Biología de los crustáceos cultivables en América Latina. En (FAO) La acuicultura en América Latina. Actas del Simposio sobre Acuicultura en América Latina, Volumen 2. Montevideo: FAO.
Conapesca (Comisión Nacional de Acuacultura y Pesca). (2013). Anuario estadístico de acuacultura y pesca, Edición 2013. Mazatlán, Sinaloa: Comisión Nacional de Acuacultura y Pesca.
De la Rosa-Vélez, J., Escobar-Fernández, R., Correa, F., Maqueda-Cornejo, M. y de la Torre-Cueto J. (2000). Genetic structure of two commercial penaeids (Penaeus californiensis and P. stylirostris) from the Gulf of California, as revealed by allozyme variation. Fishery Bulletin, 98, 674–683.
Earl, D. A. y von-Holdt, B. M. (2012). Structure harvester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources, 4, 359–361.
Evanno, G., Regnaut, S. y Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular Ecology, 14, 2611–2620.
Excoffier, L., Laval, G. y Schneider, S. (2005). Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online, 1, 47–50.
FAO. (2006). Cultured Aquatic Species Information Programme. Penaeus vannamei. Text by Briggs, M. En FAO Fisheries and Aquaculture Department [online]. Rome. Updated 7 April 2006. [Cited 21 October 2016]. http://www.fao.org/fishery/culturedspecies/Penaeus_vannamei/en
Freeland, J. R. (2005). Molecular Ecology. West Sussex, England: John Wiley & Sons.
Hallerman, E. M. (2003). Population viability analysis. En E. M. Hallerman (Ed.), Population genetics: principles and applications for fisheries scientists (pp. 403–417). Bethesda, Maryland: American Fisheries Society.
Hedrick, P. W. (2000). Genetics of populations. 2nd edition. Sudbury, Massachusetts: Jones and Bartlett Publishers.
Hendrickx, M. E. (2001). Taxonomía, biología y zoogeografía de los peneidos de importancia comercial del Pacífico mexicano. En F. Páez Osuna (Eds.), Camaronicultura y medio ambiente (pp. 24–39). Mazatlán, Sinaloa: Instituto de Ciencias del Mar y Limnología.
Holthuis, L. B. (1980). Shrimps and prawns of the World. An annotated catalogue of species of interest to fisheries. Catalogue Vol.1. Rome: FAO Fisheries Synopsis No.125.
Jombart, T. (2008). Adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics, 24, 1403–1405.
Marinone, S. G., Lavín, M. F. y Parés-Sierra, A. (2011). A quantitative characterization of the seasonal Lagrangian circulation of the Gulf of California from a three-dimensional numerical model. Continental Shelf Research, 31, 1420–1426.
Mendoza-Cano, F., Grijalva-Chon, J. M., Perez-Enriquez, R., Ramos-Paredes, J. y Varela-Romero, A. (2013). Genetic diversity of mitochondrial DNA from Litopenaeus vannamei broodstock used in northwestern Mexico. Ciencias Marinas, 39, 401–412.
Miller, M. R. y Kapuscinsky, A. R. (2003). Genetic guidelines for hatchery supplementation programs. En E. M. Hallerman (Eds.), Population genetics: principles and applications for fisheries scientists (pp. 329–355). Bethesda, Maryland: American Fisheries Society.
Novelo, N. D., Gomelsky, B. y Pomper, K. W. (2010). Inheritance and reliability of random amplified polymorphic DNA-markers in two consecutive generations of common carp (Cyprinus carpio L.). Aquaculture Research, 41, 220–226.
Paradis E., Claude J. y Strimmer K. (2004). APE: analyses of phylogenetics and evolution in R language. Bioinformatics, 20, 289–290.
Pérez-Enríquez, R., Cruz-Hernández, P. y Saucedo-Barrón, C. J. (2013). Liberación de juveniles para la repoblación de cuerpos de agua, consideraciones y recomendaciones de manejo para minimizar los riesgos de impacto genético negativo. Industria Acuícola, 9, 18–23.
Pérez-Enríquez, R., Hernández-Martínez, F. y Cruz, P. (2009). Genetic diversity status of white shrimp Penaeus (Litopenaeus) vannamei broodstock in Mexico. Aquaculture, 297, 44–50.
Perez-Enriquez, R., Medina-Espinoza J. A., Max-Aguilar A. y Saucedo-Barrón, C. J. (2018) Genetic tracing of farmed shrimp (Decapoda: Penaeidae) in wild populations from a main aquaculture region in Mexico. Revista de Biología Tropical, 66, 381–393.
Pritchard, J. K., Stephens, M. y Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics, 155, 945–959.
Ramos-Paredes, J. y Grijalva-Chon, J. M. (2003). Allozyme genetic analysis in hatchery strains and wild blue shrimp, Penaeus (Litopenaeus) stylirostris (Stimpson), from the Gulf of California. Aquaculture Research, 34, 221–234.
Sagarpa. (2012). Acuerdo por el que se da a conocer la Actualización de la Carta Nacional Pesquera. Diario Oficial de la Federación. 24 de agosto de 2012, Segunda Sección. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación. México D.F.
Saitou, N. y Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution, 4, 406–425.
Saucedo-Barrón, C. J. (2015). Metas y compromiso de gobierno en pesca y acuacultura. Instituto Sinaloense de Acuacultura y Pesca. Encuentro regional para el intercambio de innovaciones y tecnologías aplicadas en el sector pesquero y acuícola 2015. Recuperado el 05 octubre, 2016 de http://ager-net.org/EncuentrosINAPESCA2015/
Sunden, S. L. F. y Davis, S. K. (1991). Evaluation of genetic variation in domestic population of Penaeus vannamei (Boone): a comparison with three natural populations. Aquaculture, 97, 131–142.
Supungul, P., Sootanan, P., Klinbunga, S., Kamonrat, W., Jarayabhand, P. y Tassanakajon, A. (2000). Microsatellite polymorphism and the population structure of the black tiger shrimp (Penaeus monodon) in Thailand. Marine Biotechnology, 2, 339–347.
Taniguchi, N. (2004). Broodstock management for stock enhancement programs of marine fish with assistance of DNA marker (a review). En K. M. Leber, S. Kitada, H. L. Blankenship y T. Svåsand (Eds.), Stock enhancement and sea ranching. Developments, pitfalls and opportunities (pp. 329–338). Oxford, UK: Blackwell Publishing.
Valles-Jiménez, R., Cruz, P. y Pérez-Enríquez, R. (2004). Population genetic structure of Pacific White shrimp (Litopenaeus vannamei) from Mexico to Panama: microsatellite DNA variation. Marine Biotechnology, 6, 475–484.
Van-Oosterhout, C., Hutchinson, W. F., Wills, D. P. M. y Shipley, P. (2004). MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes, 4, 535–538.