|
ISSN: 1659-455X • e-ISSN: 1659-407X REVMAR REVISTA CIENCIAS MARINAS Y COSTERAS
Recibido: 8 enero 2021 • Corregido: 18 marzo 2021 Identificación de cuatro especies de Clupeiformes (Actinopterygii) mediante análisis de secuencias de ADN mitocondrial en zonas de explotación pesquera del Golfo de Nicoya, Costa Rica Identification of four species of Clupeiformes (Actinopterygii) through analysis of mitochondrial DNA sequences in fishing exploitation areas of the Gulf of Nicoya, Costa Rica Carolina Sancho-Blanco1, Luis Hernández-Noguera2, Luis Vega-Alpízar2, Rosa Soto-Rojas2 y Rodolfo Umaña-Castro1* RESUMEN Las sardinas junto con las anchovetas constituyen unos de los grupos de mayor importancia en la pesquería del golfo de Nicoya. Son alimento para peces de mayor tamaño y aves marinas, por tanto, no solo sufren explotación pesquera, sino también presentan altas tasas de mortalidad natural. Debido a la importancia comercial y ecológica que representan estas especies, es necesario contar con una correcta identificación taxonómica, principalmente para la elaboración de estudios pesqueros. Antes, la determinación de especies de peces se fundamentaba únicamente en características morfológicas externas, sin embargo, no siempre es posible, debido a su similitud morfológica, por lo que el uso combinado de datos morfológicos junto con los genéticos basados en genes mitocondriales puede contribuir en su reconocimiento. Históricamente, estudios taxonómicos han agrupado a las especies O. libertate, O. medirastre y O. bulleri como integrantes del complejo Opisthonema spp., dada la gran similitud de forma. Mientras, para la familia Engraulidae, las indagaciones taxonómicas son escasas. En el presente trabajo, se identificaron molecularmente especies del complejo Opisthonema spp. y una especie del género Cetengraulis sp., colectadas entre junio y diciembre del 2017, mediante secuenciación de genes mitocondriales: COI, Ctyb y 16S ARNr. Los hallazgos confirman, con datos moleculares, la identidad y presencia de tres especies del complejo Opisthonema spp. y una especie de anchoveta, localizadas en zonas de explotación pesquera en el golfo de Nicoya, Costa Rica, información que facilita el uso y manejo de estos recursos pesqueros en un área marina de importancia comercial en el país. Palabras clave: Opisthonema spp., Cetengraulis sp., ADNmt, ADNr, posicionamiento taxonómico ABSTRACT Sardines, together with anchovies, constitute one of the most important groups in the fishery of the Gulf of Nicoya. They are consumed by large fish and seabirds; therefore, they not only suffer from fishing exploitation but also have high rates of natural mortality. Due to their commercial and ecological importance, it is necessary to have a reliable taxonomic identification of these species, mainly for the preparation of fisheries studies. Previously, the identification of fish species was based only on external morphological characteristics; however, this is not always possible due to their morphological similarity, but the combined use of morphological and genetic data based on mitochondrial genes can contribute to their accurate recognition. Historically, taxonomic studies have grouped the species O. libertate, O. medirastre and O. bulleri as members of the Opisthonema spp. complex due to their great morphological similarity, while, there are few taxonomic studies of the family Engraulidae. In the present work, species of the Opisthonema spp. complex and a species of the genus Cetengraulis sp., collected between June-December 2017, were molecularly identified by sequencing mitochondrial genes: COI, Ctyb, and 16S rRNA. The findings of the molecular data analysis confirm the identity and presence of three species of the Opisthonema spp. complex and a species of anchovy, in areas of fishing exploitation in the gulf of Nicoya, Costa Rica, information that facilitates the use and management of these fishing resources in a marine area of commercial importance for the country. Keywords: Opisthonema spp., Cetengraulis sp., mtDNA, rDNA, taxonomic placement El orden Clupeiformes incluye especies conocidas como arenques, sardinas y anchovetas, todas de gran importancia comercial para las pesquerías globales. Asimismo, abarca familias como Chirocentridae, Clupeidae, Denticipitidae, Dussumieriidae, Engraulidae Pristigasteridae y Sundasalangidae; sin embargo, la taxonomía de especies de las familias Clupeidae y Engraulidae aún sigue siendo controversial en términos morfológicos y genéticos (Li & Ortí, 2007; Bloom & Lovejoy, 2014; Vicente et al. 2020). Tanto los Clupeidos como los Engraulidos son individuos de tallas pequeñas generalmente, se agrupan en grandes cardúmenes, son filtradores y conforman la base de la cadena alimenticia de los ecosistemas (Fischer et al. 1995); por tanto, no solo sufren de la explotación pesquera, sino que también presentan altas tasas de mortalidad natural, al ser el alimento para otros organismos como peces de mayor tamaño y aves marinas (Bussing & López, 1994). En las costas del Pacífico, el género Opisthonema Gill, 1861 está representado por cuatro especies: Opisthonema libertate Günther, 1867, Opisthonema medirastre Berry & Barrett, 1963, Opisthonema bulleri Regan, 1904 y Opisthonema berlangai Berry & Barrett, 1963 (Froese & Pauly, 2019). En Costa Rica, debido a su difícil diferenciación morfológica, históricamente las especies O. libertate, medirastre y bulleri se han agrupado como integrantes del complejo Opisthonema spp. y han sido utilizadas para la industria de enlatados con un mismo valor comercial (Soto & Rodríguez, 1999; Vega-Corrales, 2010). Por su parte, la familia Engraulidae se encuentra representada, principalmente, por Cetengraulis mysticetus, el cual es un recurso pesquero de importancia, al emplearse como carnada viva por la flota pesquera artesanal en la zona interna del golfo de Nicoya, Costa Rica (Rodríguez & Gómez, 1998; Soto & Rodríguez, 1999; Murase et al. 2014). Estudios enfocados en la identificación y taxonomía del género Opisthonema señalan si las diferencias encontradas en el nivel morfológico son el resultado de variaciones intraespecíficas o de características diagnósticas (Pérez-Quiñónez et al. 2017). Por otro lado, para la familia Engraulidae, los hallazgos en cuanto a taxonomía molecular son escasos, a pesar de su importancia pesquera y del impacto ecológico potencial que las pesquerías pueden tener en el ecosistema pelágico (Wang et al. 2018). A partir del desarrollo de técnicas moleculares, se inició la incorporación de datos genéticos, basados en el ADN mitocondrial y nuclear, que permitieran esclarecer interrogantes taxonómicas no resueltas mediante características morfológicas y merísticas (Ward et al. 2005; Teletchea, 2009; Ardura et al. 2010; Carvalho et al. 2011; Pereira et al. 2013; Pérez-Quiñónez, 2014; Durand et al. 2017; Pérez-Quiñónez et al. 2017; Goodbody-Gringley et al. 2019; Pérez-Quiñónez et al. 2019). En peces, los marcadores genéticos más utilizados para la identificación de especies se basan en regiones del ADN mitocondrial (Jahan et al. 2017), los cuales, al mostrar alta tasa de mutación y gran variabilidad interespecífica, hacen posible una adecuada distinción entre especies (Rocha et al. 2005; Paine et al. 2007). Entre dichos marcadores mitocondriales, se encuentran el gen del 16S del ARN ribosómico (16S ARNr), el gen del citocromo oxidasa I (COI), el gen citocromo oxidasa b (Cytb), entre otros. El 16S ARNr es altamente conservado; variaciones de unos pocos nucleótidos en su estructura permiten inferir relaciones taxonómicas en diferentes niveles, inclusive entre taxa estrechamente relacionados (Chakraborty & Iwatsuki, 2006). Por su parte, el gen COI, a pesar de ser muy conservado, posee regiones polimórficas que facilitan la identificación de especies (Hebert et al. 2003). Mientras, la secuencia de nucleótidos del gen Cytb contiene información específica de la especie; posibilita reconstruir tanto relaciones inter- o intraespecíficas como de nivel superior, por lo cual se considera un gen con alto rendimiento en filogenia molecular (Peng et al. 2004). La identificación taxonómica de estas especies en las capturas, a través de los diferentes artes de pesca en el golfo de Nicoya, permite evaluar el estado actual de sus poblaciones; fijar la cuota de pesca, el establecimiento correcto de las vedas; recomendar medidas de ordenamiento y pautas de captura, así como para garantizar los niveles de biomasa del stock desovante (Kochzius et al. 2010; Pérez-Quiñónez et al. 2018). De esta forma, se mantiene, además, la sostenibilidad de otras especies de interés pesquero -corvinas, pargos, jureles- y ecológico -las aves- (FAO, 2018). En el presente estudio, se identificó, molecularmente, a partir de secuencias del ADN mitocondrial (ADNmt) junto con caracteres morfológicos, algunas especies que integran las familias Engraulidae y Clupeidae. Con ello, se proporcionan elementos suficientes que instan a diferenciar entidades taxonómicas en una zona de explotación pesquera en el golfo de Nicoya, Costa Rica. La información obtenida facilitará el uso y manejo de estos recursos pesqueros en una zona marina de importancia comercial en el país, al detectar posibles sesgos en la identificación tradicional basada en caracteres morfológicos diagnósticos. Se recolectaron especímenes del complejo Opisthonema spp. y Cetengraulis sp. distribuidos específicamente en el golfo de Nicoya, Pacífico Central, Costa Rica, de junio a setiembre de 2017. Se identificaron, en la zona externa del golfo, individuos de Opisthonema bulleri (sardina azul), O. medirastre y O. libertate (sardina gallera), en un rango de longitudes estándar de 137 mm a 184 mm. Por otra parte, en la zona interna del golfo, se detectaron individuos de C. mysticetus (anchoveta) (Fig. 1). La identificación preliminar se realizó a partir de caracteres morfológicos diagnósticos, utilizando la metodología previamente descrita por Berry & Barret, 1963; Bussing & López, 1994, Robertson & Allen, 2002; Pérez Quiñones, 2014. (Cuadro 1).
Fig. 1. Áreas de captura de las especies analizadas. Fuente de datos de figura: Atlas 2014 (TEC) y Laboratorio de Análisis Biológicos Pesqueros (UNA) (QGIS, 2018) Fig. 1. Areas in which the species analyzed were captured. Data source: Atlas 2014 (TEC) and Fisheries Biological Analysis Laboratory (UNA) (QGIS, 2018) Cuadro 1. Características morfológicas diagnósticas evaluadas para las especies de Opisthonema spp. analizadas en este estudio para la localidad del Golfo de Nicoya, Costa Rica Table 1. Diagnostic morphological characteristics evaluated for Opisthonema spp., analyzed in this study for the location of the Gulf of Nicoya, Costa Rica
bc: branquiespinas ceratobranquiales, std: longitud estándar, ph: pecas humerales (ausentes = /, presentes = +) Extracción de ADN, amplificaciones y secuenciación de genes mitocondriales Se realizó una extracción de ADN total de especies del complejo Opisthonema spp., O. medirastre (n = 11), O. libertate (n = 4), O. bulleri (n = 7) y de la familia Engraulidae C. mysticetus (n = 5), a partir de 30 mg de tejido, aproximadamente, mediante el kit comercial PureLink Genomic DNA Mini Kit (Thermo Fisher Scientific), bajo las recomendaciones del fabricante. La cantidad de muestras seleccionadas por especie, en este estudio, se respalda en un trabajo previamente realizado por Bingpeng et al. 2018. El ADN obtenido fue amplificado mediante una PCR punto final (volumen final: 25 µL) compuesta por 1X de PCR Máster Mix, 0.9 µM de cada cebador y alrededor de 100 ng de ADN. Se analizaron los genes mitocondriales citocromo oxidasa subunidad I (COI) y citocromo b (Cytb), así como el gen del ARN ribosomal 16S (ARNr 16S) (Cuadro 2). La PCR fue llevada a cabo en las siguientes condiciones: 94°C por 3 min (desnaturalización inicial); 40 ciclos de 94°C por 30 s (desnaturalización), 48/55°C por 30 s (alineamiento), 72°C por 45 s (extensión) y, finalmente, 72°C por 10 min como extensión final. Los productos de PCR se resolvieron por electroforesis en gel de agarosa al 1.5%. La corrida se ejecutó a 80 V por 60 min en una solución amortiguadora TBE 1X (pH 8.0). Los amplicones obtenidos de los tres genes fueron secuenciados en ambas direcciones, utilizando la secuenciación de Sanger en un analizador genético ABI 3130, bajo una química de secuenciación BigDye Terminator v3.1. Cuadro 2. Cebadores utilizados en este estudio y sus respectivas secuencias de ADN. Se indica la región mitocondrial (mt) objetivo, las secuencias del cebador en dirección 5ʹ-3ʹ, la referencia (Ref.) y el tamaño esperado del amplicón (Te) en pares de bases (pb) Table 2. Primers used in this study and their respective DNA sequences. Target mitochondrial region (mt), 5ʹ-3ʹ primer sequences direction, reference (Ref.), and expected amplicon size (Te) in base pairs (bp) are indicated
Análisis molecular y de posicionamiento taxonómico Utilizando el programa Geneious versión R9 (Biomatters), se editaron las secuencias nucleotídicas obtenidas, mediante una edición manual de los bordes no alineados. Luego, para corroborar, preliminarmente, su identidad taxonómica en el nivel de género, esas secuencias fueron evaluadas en el programa en línea BLASTn (Altschul et al. 1990), disponible en el National Center for Biotechnology Information (NCBI) (Benson et al. 2012); al final, fueron depositadas en el GenBank (Cuadro 3). Por otra parte, se realizó un alineamiento múltiple de secuencias, usando el programa MAFFT 7.0 (Katoh et al. 2009) y aplicando el método iterativo de refinamiento G-INS-i (1PAM/k = 2). Para determinar el mejor modelo de sustitución nucleotídica, se analizó el alineamiento múltiple mediante el JModelTest v2.1.10 (Darriba et al. 2012), empleando parámetros por defecto, bajo el criterio de información de Akaike con 95% de intervalo de confianza y tomando en cuenta la corrección AICc. El programa arrojó los seis mejores modelos para las matrices correspondientes a los genes 16S ARN ribosomal [GTR + G, GTR + I, GTR + I + G, TVM + I, TVM + I + G, TIM3 + G], citocromo oxidasa subunidad 1 [TVM + G, TPM2uf + G, TVM + I + G, TPM2uf + I + G, GTR + G, TIM2 + G], para el cual TVM + G = GTR + G, y citocromo oxidasa b [TrN + G, TIM1 + G, TrN + I + G, TIM3 + I + G, TIM2 + I + G, GTR + G], para el cual TrN + G = GTR + G. La elaboración de las topologías multigénicas de posicionamiento taxonómico se realizó mediante dos métodos: el primero, utilizando el programa raxmlGUI v.7.4.2 (Stamatakis et al. 2005), bajo un algoritmo de máxima verosimilitud (ML), con un modelo reversible en tiempo general (GTR-GAMMA) y aplicando 5000 permutaciones como valor bootstrap. El segundo método se efectuó a través de inferencia bayesiana (BI), con el programa MrBayes v.3.2.6 (Huelsenbeck & Ronquist, 2001), un modelo de sustitución GTR-GAMMA y los siguientes parámetros: nst = 6, rates = gamma, cuatro cadenas simultáneas de Markov Monte Carlo (mcmc) y ngen = 1000000. La visualización y edición de los árboles obtenidos se confeccionó con el programa FigTree v1.4 (Rambaut, 2009); después, se concatenaron manualmente, según la similitud de las topologías y utilizando secuencias nucleotídicas ya reportadas en la base de datos del NCBI (Benson et al. 2012). Las características morfológicas diagnósticas evaluadas para las tres especies del género Opisthonema spp. muestran que su identificación, tras observar las pecas humerales, permite distinguir O. libertate (presencia de pecas humerales, ph) de O. medirastre y O. bulleri (ausencia de ph); sin embargo, este análisis se debe realizar con el espécimen fresco o recién capturado, de lo contrario, se dificulta su diferenciación. Por otro lado, medir la longitud estándar (std) y el conteo de las unidades de branquiespinas ceratobranquiales (bc) muestra una media de 86.3 para O. libertate (n = 3); 68.0, O. medirastre (n = 9), y 36.6, O. bulleri (n = 6). Estos resultados se encuentran entre los intervalos de unidades bc y longitud estándar señalados en la literatura (Cuadro 1). A partir de las reacciones de PCR, fue posible obtener amplificaciones positivas (tamaño esperado del amplicón en pares de bases), para los tres cebadores utilizados (16S, COI y Cytb), en las especies O. medirastre, O. libertate y O. bulleri del complejo Opisthonema spp., así como para una especie de la familia Engraulidae. Los amplicones recuperados se purificaron y las secuencias parciales conseguidas se depositaron en el Genbank (Cuadro 3). Cuadro 3. Secuencias depositadas en el Genbank (https://www.ncbi.nlm.nih.gov/genbank) de las tres especies del complejo Opisthonema spp., derivadas de tres marcadores moleculares y una especie del género Cetengraulis sp. obtenidas de este estudio Table 3. Sequences deposited on Genbank (https://www.ncbi.nlm.nih.gov/genbank) of three species of Opisthonema spp. complex, derived from three molecular markers and one species from the genus Cetengraulis sp. obtained from this study
Las topologías de los árboles filogenéticos construidos con base en regiones parciales de los genes Cytb, COI y 16S (Figuras 2A, 2B, 2C) confirman la existencia, en el golfo de Nicoya, de un grupo monofilético conformado por tres clados genéticamente diferentes, según la alta divergencia de las agrupaciones formadas. Tales clados se definen como O. medirastre, O. libertate y O. bulleri; así, se demuestra la presencia de las tres especies que componen el complejo complejo Opisthonema spp. en el golfo de Nicoya, Puntarenas. Además, se observa que para los genes Cytb y 16S las especies O. libertate y O. medirastre son similares entre sí a nivel genético. Estos resultados moleculares son congruentes con los datos morfométricos y merísticos obtenidos, acorde con las características diagnósticas (pecas humerales y branquiespinas ceratobranquiales) indicadas por Bussing y Lopez (1994). Por otro lado, los individuos de C. mysticetus, analizados mediante marcadores mitocondriales, se posicionan taxonómicamente en el clado marino Cetengraulis: Cetegraulis edentulus y C. mysticetus, separando en ramas bien soportadas a los especímenes originarios del Pacífico y del Atlántico, demostrando la especificidad en el método de identificación taxonómica con herramientas moleculares.
Fig. 2. Topologías de posicionamiento taxonómico mediante Máxima Verosimilitud/Inferencia Bayesiana (ML/BI), entre secuencias genéticas de Opisthonema spp., Cetengraulis sp. del Pacífico costarricense y secuencias parciales de los genes A) citocromo b (Cytb), B) citocromo oxidasa subunidad 1 (COI) y C) 16S ARN ribosomal (16S), obtenidas de la base de datos Genbank. El número en las ramas indica el valor de soporte boostrap (porcentaje derivado de 5 000 permutaciones para ML y 1 000 000 generaciones MCMC para BI). Sardinella aurita es definido como grupo externo. La región y el código de colecta de las muestras de este estudio se presentan entre paréntesis. Para las secuencias obtenidas del Genbank, los números de accesión se muestran entre corchetes. Valores de bootstrap o de probabilidad posterior menor a 50% o no obtenidos no se muestran en la topología o se indican con - Fig. 2. Taxonomic placement topologies using Maximum Likelihood / Bayesian Inference (ML / BI), between genetic sequences of Opisthonema spp., Cetengraulis sp. from the Costa Rican Pacific and partial sequences of the genes A) cytochrome b (Cytb), B) cytochrome oxidase subunit 1 (COI) and C) 16S ribosomal RNA (16S), obtained from the Genbank database. The number in the branches indicates the bootstrap support value (percentage derived from 5 000 permutations for ML and 1 000 000 MCMC generations for BI). Sardinella aurita is defined as an external group. The region and collection code of the samples in this study are shown in parentheses. For sequences obtained from the Genbank, the accession numbers are shown in square brackets. Boostrap or posterior probability values less than 50% or not obtained are not shown in the topology or indicated with a - Según Berry y Barret (1963), las diferenciaciones taxonómicas de estas especies se basan, principalmente, en dos aspectos: distribución geográfica y número o cantidad de branquiespinas de la región ceratobranquial (bc), asociadas a la longitud del individuo. Se reporta que O. libertate presenta 71-107 bc; O. medirastre, 50-68 bc, y O. bulleri, 28-35 bc. Posteriormente, en un estudio realizado por Rodríguez-Sánchez (1985), se considera que la densidad de branquiespinas disminuye en relación con el tamaño del espacio interbranquial, hecho que discrepa de la veracidad de los caracteres merísticos descritos por Berry y Barret (1963). Fue, entonces, menor la cantidad de branquinespineas en O. bulleri, seguido de O.medirastre, y mayor en O.libertate, por lo que se concluye que el conteo de las branquiespinas permite identificar únicamente O. libertate y O. bulleri. En este estudio, el conteo de las branquiespinas de la región ceratobranquial concuerda con los rangos de las unidades bc para cada especie del complejo Opisthonema spp., señaladas en la literatura. De la misma forma, los valores medidos de la longitud estándar se encuentran entre los límites descritos para las especies en estudio, lo que facilita la identificación en especímenes recién capturados (Cuadro 1). Por otro lado, la presencia de las pecas humerales permitió distinguir, con mayor facilidad, entre los especímenes de O. libertate y los de O. medirastre y O. bulleri. Una investigación realizada por Lagúnez y Rodríguez (1992), en poblaciones de estas mismas especies en el pacífico mexicano, a partir de análisis de proteinogramas, arrojó la existencia de seis patrones en O. libertate muy similares a O. medirastre; mientras, estos difieren de los de O. bulleri, por lo tanto, los autores consideran que las dos primeras especies son entidades taxonómicamente muy cercanas y que más que especies distintas podrían ser subpoblaciones de una misma especie. Por su parte, Pérez-Quiñónez (2014), a través de patrones morfométricos de las tres especies y un análisis de componentes principales mediante variables canónicas, obtuvo una mayor similitud entre las formas de O. medirastre y O. libertate, mientras que O. bulleri presentaba más diferenciación con respecto a la forma corporal. Nuestros hallazgos muestran una mayor cercanía filogenética entre O. libertate y O. medirastre, de acuerdo con las topologías taxonómicas obtenidas para los genes 16S y Cytb, lo cual concuerda con lo reportado por los autores citados anteriormente. La adición de datos moleculares en estudios morfológicos mejora significativamente el soporte y la resolución de los análisis filogenéticos (Wortley & Scotland, 2006); con ello, proporciona datos relevantes cuando se evalúan especies de interés comercial. En estas últimas, no es posible (o resulta complejo) determinar las capturas por especie mediante características morfológicas (Pérez-Enríquez et al. 2016), como es el caso de las sardinas de la familia Clupeidae, razón por la cual dichas capturas son registradas, generalmente, como sardina crinuda (Opisthonema spp.) tanto en Costa Rica como en otros países de la región -México entre ellos- (Pérez-Quiñonez et al. 2020). Lo anterior conduce a una posible subestimación o sobreestimación de la riqueza de especies e incluso a la sobreexplotación de las poblaciones de peces. Li y Ortí (2007), a partir de análisis de máxima parsimonia, máxima verosimilitud y análisis bayesianos de datos con ADN mitocondrial (12S y 16S) y ADN nuclear (genes activadores de recombinación, RAG 1 y RAG 2), encontraron monofilia en subfamilias Engraulidae (Engraulis y Anchoa), pero no en la familia Clupeidae. Pérez-Quiñónez et al. (2017), completaron análisis moleculares del gen citocromo oxidasa I (COI) del ADNmt para especies reconocidas del complejo Opisthonema spp., en la región sur del golfo de California, México, y con fundamento en sus resultados morfológicos y genéticos se soporta la existencia de tres entidades evolutivas para la región: O. bulleri, O. medirastre y O. libertate. Estos datos coinciden con nuestros hallazgos, de acuerdo con la topología observada en el árbol de posicionamiento taxonómico para los tres genes (Cytb, COI y 16S ARNr), lo que apoya la identidad y existencia de las tres especies para la región y, específicamente, para el golfo de Nicoya, Costa Rica. La correcta identificación de estas especies de interés pesquero, mediante características morfológicas y genéticas, facilita el uso y manejo de estos recursos marinos dentro de una zona de relevancia comercial en el país como lo es el golfo de Nicoya. Ello con el objetivo de minimizar la alteración del porcentaje de la biomasa virginal que debe contener el ecosistema y así mantener el resto de los eslabones de la cadena alimentaria. Con base en los resultados obtenidos en este estudio sobre los caracteres morfométricos y merísticos diagnósticos, así como en herramientas moleculares, se confirma la identificación de tres entidades taxonómicas diferentes en el complejo Opisthonema spp. descrito en el golfo de Nicoya: O. libertate, O. medirastre y O. bulleri, confirmando su identidad y presencia en sitios de pesca en el exterior de dicho golfo. Por otro lado, se corrobora el posicionamiento taxónomico de los individuos de C. mysticetus en el clado marino Cetengraulis sp., C. edentulus y C. mysticetus, lo cual comprueba su identidad y presencia en los sitios de pesca dentro del golfo de Nicoya. Con la información obtenida en este trabajo y la disponible en la base mundial de datos genéticos, es posible realizar, a futuro, investigaciones comparativas de diversidad y estructura genética (por ejemplo, redes de haplotipos y RAD-seq), en el golfo de Nicoya y otras zonas de captura de Clupeidos de importancia comercial. Al proyecto 0079-14: Aportes a la gestión pesquera artesanal y semi-industrial del golfo de Nicoya, Costa Rica, financiado por la Ley de Pesca, Costa Rica. También, a Karolina Ramos Jiménez y Karen Oviedo Bolaños, por el soporte técnico en los ensayos moleculares. Por último, a los evaluadores, por sus recomendaciones y aportes para la publicación de este artículo. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol., 215(3), 403-410. https://doi.org/10.1016/S0022-2836(05)80360-2 Ardura, A., Linde, A. R., Moreira, J. C. & García-Vázquez, E. (2010). DNA bar-coding for conservation and management of Amazonian commercial fish. Biol. Conserv., 143(6), 1438-43. https://doi.org/10.1016/j.biocon.2010.03.019 Baldwin, C. C., Mounts, J. H., Smith, D. G. & Weigt, L. A. (2009). Genetic identification and color descriptions of early life/history stages of Belizean Phaeoptyx and Astrapogon (Teleostei: Apogonidae) with comments on identification of adult Phaeoptyx. Zootaxa, 2008, 1-22. Benson, D. A., Cavanaugh, M., Clark, K., Karsch-Mizrachi, I., Lipman, D. J., … & Sayers, E. W. (2012). GenBank. Nucleic Acids Res., 41(D1), D36-D42. https://doi.org/10.1093/nar/gks1195 Berry, F. H. & Barrett, I. (1963). Análisis de las branquiespinas y denominación del arenque de hebra Opisthonema. Inter-A. Tuna Trop. Comm. Bull., 7(2), 110-190. Bingpeng, X., Heshan, L., Zhilan, Z., Chunguang, W., Yanguo, W. & Jianjun, W. (2018). DNA barcoding for identification of fish species in the Taiwan Strait. PLoS ONE, 13(6), e0198109. https://doi.org/10.1371/journal.pone.0198109 Bloom, D. D. & Lovejoy, N. R. (2014). The evolutionary origins of diadromy inferred from a time/calibrated phylogeny for Clupeiformes (herring and allies). Proc. R. Soc. Lond. B. Biol. Sci., 281(1778), 2013-2081. Bussing, W. & López, M. (1994). Peces demersales y pelágicos costeros del Pacífico de Centroamérica Meridional. Guía Ilustrada. Rev. Biol. Trop., 47 (Supplemento), 47-164. Carvalho, D. C., Oliveira, D. A., Pompeu, P. S., Leal, C. G., Oliveira, C. & Hanner, R. (2011). Deep barcode divergence in Brazilian freshwater fishes: the case of the São Francisco River basin. Mitochondrial DNA, 22(1), 80-86. Chakraborty, A. & Iwatsuki. Y. (2006). Genetic variation at the mitochondrial 16S rRNA gene among Trichiurus lepturus (Teleostei:Trichiuridae) from various localities: Preliminary evidence of a new species from West Coast of Africa. Hydrobiologia, 563, 501-513. https://doi.org/10.1007/s10750-006-0105-4 Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. (2012). jModelTest 2: more models, new heuristics, and parallel computing. Nat. Methods, 9(8), 772. https://doi.org/10.1038/nmeth.2109 Durand, J. D., Hubert, N., Shen, K. N. & Borsa, P. (2017). DNA barcoding grey mullets. Rev. Fish. Biol. Fish., 27, 233-43. https://doi.org/10.1007/s11160-016-9457-7 FAO. Organización de las Naciones Unidas para la Alimentación y la Agricultura. (2018). El estado mundial de la pesca y la acuicultura 2018. Cumplir los objetivos de desarrollo sostenible. Italia: FAO. Fischer, W., Krupp, F., Schneider, W., Sommer, C., Carpenter, K. E. & Miem, V. H. (1995). Guía FAO para la identificación de especies para los fines de la pesca. Pacífico Centro Oriental. Italia: FAO. Froese, R. & Pauly, D. (2019). FishBase. World Wide Web electronic publication. Goodbody-Gringley, G., Strand, E. & Pitt, J. M. (2019). Molecular characterization of nearshore baitfish populations in Bermuda to inform management. PeerJ, 7, e7244. Hebert, P. D., Ratnasingham, S. & de Waard, J. R. (2003). Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lond. B. Biol. Sci., 270 Suppl 1(Suppl 1), S96-S99. https://doi.org/10.1098/rsbl.2003.0025 Huelsenbeck, J. P. & Ronquist, F. (2001). MRBAYES: Bayesian inference of phylogeny. Bioinf., 17(8), 754-755. https://doi.org/10.1093/bioinformatics/17.8.754 Ivanova, N. V., deWaard, J. R. & Hebert, P. D. N. (2006). An inexpensive automatization friendly protocol for recovering high quality DNA. Mol. Ecol. Notes, 6(4), 998-1002. https://doi.org/10.1111/j.1471-8286.2006.01428.x Jahan, H., Akter, M., Begum, R. A. & R. Shahijahan. (2017). Identification and comparison of three carp fishes based on mitochondrial 16S rRNA gene. J. Biol. Sci., 26 (2): 167-174. Katoh, K., Asimenos, G. & Toh, H. (2009). Multiple alignment of DNA sequences with MAFFT. In D. Posada (Ed.), Bioinformatics for DNA Sequence Analysis. Methods in Molecular Biology (Methods and Protocols) (pp. 39-64). EE. UU.: Humana Press. https://doi.org/10.1007/978-1-59745-251-9_3 Kocher, T. D., Thomas, W. K., Meyer, A., Edwards, S. V., Paabo, S., … & Wilson, A. C. (1989). Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc. Nat. Acad. Sci. U. S. A., 86(16), 6196-6200. https://doi.org/10.1073/pnas.86.16.6196 Kochzius, M., Seidel, C., Antoniou, A., Kumar, S., Campo, D., Cariani, A., ... & Blohm, D. (2010). Identifying fishes through DNA barcodes and microarrays. Plos One, 5(9), e12620. https://doi.org/10.1371/journal.pone.0012620 Lagúnez, L. & Rodríguez, F. (1992). Contribución al conocimiento genético de Opisthonema spp. del noreste de México. Inv. Mar. CICIMAR, 7(1),15-24. Li, C. & Ortí, G. (2007). Molecular phylogeny of Clupeiformes (Actinopterygii) inferred from nuclear and mitochondrial DNA sequences. Mol. Phylogenet. Evol., 44, 386-398. https://doi.org/10.1016/j.ympev.2006.10.030 Meyer, A., Kocher, T. D., BasasibwakI, P. & Wilson, A. C. (1990). Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial DNA sequences. Nature, 347(6293), 550-553. https://doi.org/10.1038/347550a0. Murase, A., Angulo, A., Miyazaki, Y., Bussing, W. & López, M. (2014). Marine and estuarine fish diversity in the inner Gulf of Nicoya, Pacific coast of Costa Rica, Central America. Check List, 10(6), 1401-1413. https://www.biotaxa.org/cl/article/view/10.6.1401 Doi: https://doi.org/10.15560/10.6.1401 Paine, M., McDowell, J. & Graves, J. (2007). Specific identification of Western Atlantic Ocean scombrids using mitochondrial DNA cytochrome c oxidase subunit I (COI) gene region sequences. Bull. Mar. Sci., 80(2), 353-367. Palumbi, S. R. (1996). Nucleic acids II: the polymerase chain reaction. In D. M. Hillis, C. Moritz & B. K. Mable (Eds.), Molecular systematics (pp. 205-247). EE. UU.: Sinauer & Associates, Inc. Peng, Z., He, S. & Zhang, Y. (2004). Phylogenetic relationships of glyptosternoid fishes (Siluriformes: Sisoridae) inferred from mitochondrial cytochrome b gene sequences. Mol. Phylogenet. Evol., 31(3), 979-987. https://doi.org/10.1016/j.ympev.2003.10.023 Pereira, L. H. G., Hanner, R., Foresti, F. & Oliveira, C. (2013). Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna? BMC Genet., 14, 20. https://doi.org/10.1186/1471-2156-14-20 Pérez-Enríquez, R., Díaz-Viloria, N., Cruz-Hernández, P., Aranceta-Garza, F., Gutiérrez-González, J. L., … & Max-Aguilar, A. (2016). Estudios de genética en poblaciones de abulón y sus aplicaciones en ordenamiento pesquero. Rec. Nat. Soc., 2(2), 24-39. https://doi.org/10.18846/renaysoc.2016.02.02.02.0002 Pérez-Quiñonez, C. I. (2014). Discriminación de las especies del género Opisthonema Gill, 1861 en el Sur del Golfo de California usando análisis morfométricos y genéticos. (Tesis de Maestría no publicada). Centro Interdisciplinario de Ciencias Marinas (CICIMAR), Instituto Politécnico Nacional, México. Pérez-Quiñonez, C. I., Quiñonez-Velázquez, C. & García-Rodríguez, F. J. (2020). A simple method for the genetic identification of commercially important species in the Opisthonema genus Gill, 1861 in the southern Gulf of California. Ciencias Marinas, 46(3), 145-154. https://doi.org/10.7773/cm.v46i3.3059 Pérez-Quiñonez, C. I., Quiñonez-Velázquez, C. & García-Rodríguez, F. J. (2018). Detecting Opisthonema libertate (Günther, 1867) phenotypic stocks in northwestern coast of Mexico using geometric morphometrics based on body and otolith shape. Lat. Am. J. Aquat. Res., 46(4), 779-790 Pérez-Quiñonez, C. I., Quiñonez-Velázquez, C. & García-Rodríguez, F. J. (2019). Genetic homogeneity of the Pacific thread herring (Opisthonema libertate) (Günther, 1867) in the Eastern Pacific, inferred from mtDNA sequences. Mitochondrial DNA Part A, 30(3), 517-524. https://doi.org/10.1080/24701394.2019.1570173 Pérez-Quiñónez, C. I., Quiñónez-Velázquez, C., Ramírez-Pérez, J. S, Vergara-Solana, F. J. & García-Rodríguez, F. J. (2017). Combining geometric morphometrics and genetic analysis to identify species of Opisthonema Gill, 1861 in the eastern Mexican Pacific. J. Appl. Ichthyol., 33(1), 84-92. QGIS. Development Team. (2018). QGIS Geographic Information System. EE. UU: Open Source Geospatial Foundation Project. http://qgis.osgeo.org Rambaut, A. (2009). FigTree v1.4 2012-2014: Tree Figure Drawing Tool. Edinburgh. Scotland: Institute of Evolutionary Biology, University of Edinburgh. http://tree.bio.ed.ac.uk/software/figtree Robertson, D. R. & Allen, G. R. (2002). Shore fishes of the Tropical Eastern Pacific: an Information System. Smithsonian Tropical Research Institute. Rocha, A., Garber, N., Garber, A. & Stuck, K. (2005). Structure of the mitochondrial control region and flanking tRNA genes of Mugil cephalus. Hidrobiológica, 15(2), 139-149. Rodríguez, J. A. & Gómez, K. R. (1998). Aspectos relevantes en la biología de Cetengraulis mysticetus (Günther) (Pisces: Engraulidae) en el Golfo de Nicoya, Costa Rica. Uniciencia, 15(1), 57-60. Rodríguez-Sánchez, R. (octubre, 1985). Aspectos de dinámica poblacional en apoyo de la separación de las especies del género Opisthonema por el método de índice de densidad de branquiespinas. CalCOFI abstracts, annual, conference. EE. UU. Soto, R. & Rodríguez, J. (1999). Dinámica poblacional de Opisthonema medirastre (Pisces: Clupeidae) en la costa Pacífica de Costa Rica. Uniciencia, 15(16), 61-64. Stamatakis, A., Ludwig, T. & Meier, H. (2005). RAxML/III: a fast program for maximum likelihood/ based inference of large phylogenetic trees. Bioinf., 21(4), 456-463. Teletchea, F. (2009). Molecular identification methods of fish species: reassessment and possible applications. Rev. Fish. Biol. Fish., 19(3), 265-293. https://doi.org/10.1007/s11160-009-9107-4 Vega-Corrales, L. A. (2010). Evaluación poblacional del stock explotable del complejo Opisthonema (Pisces: Cupleidae) en el Golfo de Nicoya, Costa Rica. Rev. Cien. Mar. Cos., 2, 83-94. Vicente, F., Loeb, M. V., Paiva, A. C. G. D., Sampaio, C. L., Argolo, L. A. & Jacobina, U. P. (2020). Integrative systematics unveils the controversial identity of Engraulidae fishing stocks in a Neotropical estuary, northeast Brazil. Neotrop. Ichthyol., 18(4), 1-17. https://doi.org/10.1590/1982-0224-2020-0037 Wang, P., Zhao, C., Fan, S., Yan, L. & Qiu, L. (2018). The complete mitochondrial genome of Thryssa hamiltonii and phylogenetic analysis of Engraulidae (Clupeiformes; Clupeoidei). Mitochondrial DNA Part B., 3(2), 538-540. https://doi.org/10.1080/23802359.2018.1467232 Ward, R. D., Zemlak, T. S., Innes, B. H., Last, P. R. & Hebert, P. D. N. (2005). DNA barcoding Australia’s fish species. Phil. Trans. R. Soc., 360(1462), 1847-1858. https://doi.org/10.1098/rstb.2005.1716 Wortley, H. & Scotland, R. (2006). The effect of combining molecular and morphological data in Published Phylogenetic Analyses. Syst. Biol., 55(4), 677-685. https://doi.org/10.1080/10635150600899798 1 Universidad Nacional. Escuela de Ciencias Biológicas, Laboratorio de Análisis Genómico (LAGEN). 2 Universidad Nacional. Estación de Biología Marina Juan Bertoglia Richards, Escuela de Ciencias Biológicas, Laboratorio de Análisis Biológicos-Pesqueros.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
/revmar13-1.jpg)
/Art3-Fig-1.png)
/Art3-Fig-2a.png)
/Art3-Fig-2b.png)
/Art3-Fig-2c.png)
