INTRODUCCIÓN

La ontogenia de la mayoría de las larvas de peces se caracteriza por cambios drásticos que determinan su capacidad fisiológica y de comportamiento durante el desarrollo (^{Peña, y Dumas, 2009}). Los cambios funcionales y corporales de las larvas en los primeros días después de la eclosión, están relacionados con patrones característicos de crecimiento alométrico y su cronología se asocia con importantes eventos en su historia de vida (^{Gisbert, 1999}). El análisis del crecimiento de las larvas, en etapas iniciales de desarrollo, proporciona información sobre su comportamiento y ecología y puede utilizarse para evaluar relaciones filogenéticas (^{Gisbert, y Doroshov, 2006}). Se conoce que existen variaciones interespecíficas en los tiempos de desarrollo y funcionabilidad, por lo que es importante orientar estudios de crecimiento para optimizar las técnicas de cría de larvas y de condiciones de alimentación (^{Padrós, Villalta, Gisbert, y Estévez, 2011}); así, el crecimiento de una especie está determinado por la respuesta celular a varios factores internos y externos, como se ha observado en larvas de peces y anfibios (^{Mansano, De Stéfani, Pereira, Nascimento, y Macente, 2014}; ^{Mansano et al, 2017}). Al respecto, la explicación del crecimiento en larvas de peces y anfibios usa modelos matemáticos que pretenden hacer una descripción, control y predicción de las variables asociadas (peso, talla, entre otras) (^{Mansano et al, 2017}). Los patrones de crecimiento alométrico son los más comunes en larvas de peces (^{Fuiman, 1983}) y representan prioridades en sus funciones vitales (^{Osse y van den Boogaart, 1995}; ^{Osse, van den Boogaart, van Snik, y van der Sluys, 1997}). Para larvas de peces de agua dulce los reportes de ^{Bettinelli-Nogueira, Lima-Godinho, y Pereira-Godinho (2014)}, ^{Oliveira, Bialetzki, Gomes, Santin, y Taguti (2012)}; ^{Souza, Melo, Caramaschi, Andrade, y Monteiro (2016)} y ^{Souza, Caramaschi, y Monteiro (2017)}, sustentan lo descrito anteriomente.

De manera general, la historia de vida de los peces puede ser malinterpretada si únicamente se consideran, para su estudio, sus formas adultas (^{Balon, 1981}, 1986). Además, pequeñas fallas en la identificación de larvas de peces pueden generar grandes errores en la interpretación de fenómenos ecológicos y taxonómicos, en vista de que el estado larval puede tener diferentes funciones ecológicas en diferentes animales (^{Orton, 1953}). Por lo tanto, considerar el desarrollo temprano en las larvas puede anexar detalles que ayudan a validar su identificación (^{Powles y Markle, 1983}). Se sabe que la anatomía larval es dinámica y que usualmente posee pocas tipologías asociadas a los adultos, haciendo su identificación más compleja; para ello, estrategias como las de reproducir en ambientes controlados u obtener huevos y larvas de padres conocidos es una aproximación válida para la identificación de estadios tempranos de desarrollo (^{Leiby, 1984}).

Todos estos cambios, regulados por la relación ambiente-genotipo (^{Gilbert, y Bolker, 2003}), determinan fenotipos que varían en sus coeficientes de crecimiento, cambiando sus proporciones alométricas (^{Gisbert, y Doroshov, 2006}) y la forma cómo las larvas aprovechan el ambiente (^{Devlin, Vandersteen, Uh, y Stevens, 2012}). De esta manera, la característica más importante para la identificación del desarrollo temprano es la merística, ya que puede ser usada a diferentes niveles taxonómicos (^{Kendall, Ahlstrom, y Moser, 1983}).

La dorada Brycon moorei, es una especie en estado de vulnerabilidad (categoría A2c,d -Nacional en Peligro Crítico (A2c) - Regional) por la reducción de su tamaño poblacional en los últimos años (^{Mojica, Usma, Álvarez-León, y Lasso, 2012}); además, no se encuentran reportes sobre su morfología, patrones de crecimiento y relaciones corporales. Por lo tanto, el objetivo del presente trabajo fue realizar el seguimiento morfométrico completo y las relaciones merísticas en su desarrollo temprano para sentar las bases de su gestión y conservación en el medio natural, así como también para su manejo en acuicultura.

MATERIALES Y MÉTODOS

Localización. La fase experimental, con peces vivos, se desarrolló en la Estación Piscícola Doradal localizada a 165 km de Medellín, en el corregimiento de Doradal, Municipio de Sonsón, Antioquia, Colombia; a una altitud de 150 m s. n. m.

Plantel de reproductores, mantenimiento en cautiverio, eventos reproductivos y desarrollo embrionario. Reproductores de dorada de un año de edad del río Magdalena, completamente adaptados al confinamiento, fueron alimentados con una dieta comercial del 38% de proteína cruda (PC) para tilapia (Oreochromis sp.), cuya composición nutricional fue 38% de PC mínimo; grasa mínima 4%; cenizas máximas 12%; humedad máxima 13% y fibra máxima 4%. El registro y la evaluación de los eventos reproductivos se realizaron según ^{David-Ruales and Castañeda-Álvarez, (2014)}. El seguimiento del desarrollo embrionario hasta el momento de la eclosión de la larva se hizo, con algunas modificaciones, según lo reportado por ^{Kimmel, Ballard, Kimmel, Ullmann, y Schilling (1995)} y ^{Tsai, Chang, Liu, Abe, y Ota (2013)}. De los huevos fecundados se obtuvieron las variables: diámetro total del huevo (ED), diámetro del vitelo (YD) y espacio perivitelino (PS), el cual fue caracterizado de acuerdo con su participación en el volumen total del huevo (^{Nakatani, Agostinho, Baumgartner, Bialetzki, Sanches and Cavicchioli, 2001}). (Figura 1A y 1B).

Figura 1 Diagrama de las medidas tomadas: A. huevo. B. juveniles de Brycon moorei. ED: diámetro total del huevo, YD: diámetro del saco de vitelo, PS: espacio perivitelino, TL: longitud total, SL: longitud estándar, Sn.L: longitud hocico, ED: diámetro ojo, HL: longitud cabeza, HH: altura cabeza, BH: altura del cuerpo, SnDF: distancia hocico a aleta dorsal, SnPF: distancia hocico a aleta pectoral, SnPeF: distancia hocico a aleta pélvica, SnAF: distancia hocico a aleta anal.  

Larvicultura, alevinaje, protocolo alimenticio y recolección de muestras. Se realizaron en un sistema de recirculación, con 25 tanques de 60 L de volumen efectivo, unidos a un sistema de filtración con desgasificador, filtro biomecánico y lámpara UV, el cual mantuvo los parámetros ambientales dentro de rangos de confort de la especie: OD (5,5 ± 0,5 mg L^-1), pH (6,5 ± 0,8), temperatura en °C (26 ± 0,4), alcalinidad y dureza (20,8 ± 1,7 mg L^-1 CaCO₃) y NH₃ (0,002 ± 0,00001 mg L^-1). El registro de los parámetros anteriores de calidad del agua se realizó una vez al día a las 7:00 am, con un equipo medidor de calidad de agua (Kit Hach FF2) y con una sonda multiparamétrica YSI (Profesional plus). El caudal de entrada de cada tanque fue controlado por medio de flujómetros (Parker MR) y se mantuvo en 3,5 L min^-1; las demás características ambientales y de manejo del sistema se realizaron según ^{David-Ruales, y Castañeda-Álvarez (2014)}. Para ello, a las larvas de 26 h post eclosión (HPE) (17 Larvas L^-1), se les realizaron biometrías cada 3 h hasta las 47 HPE; cada 4 h hasta las 71 HPE; cada 6 h hasta las 95 HPE; cada 8 h hasta las 111 HPE; cada 12 h hasta las 135 HPE y una muestra diaria hasta el día 15 (^{Faustino, 2010}; ^{Neumann, 2008}). A todas las larvas se les aplicó un protocolo alimenticio ad libitum que en los tres primeros días incluyó larvas forrajeras de bocachico (Prochilodus magdalenae), luego una combinación con Artemia salina y, a partir del sexto día, una ración con alimento balanceado del 45% de PC, que incrementó hasta el día 15; la Tabla 1 ilustra la técnica utilizada para la alimentación de las larvas.

Las larvas recolectadas (10) se clasificaron con base en lo descrito por ^{Kendall et al (1983)} y modificado por ^{Nakatani et al (2001)}. El grado de flexión de la notocorda y el desarrollo de estructuras como la aleta caudal y las escamas, entre otras, determinaron siete etapas: 1. Larval Vitelino (LV), etapa entre la eclosión y el inicio de la alimentación exógena; 2. Pre-Flexión (PF), etapa entre el inicio de la alimentación exógena, hasta el inicio de la flexión de la región terminal de la notocorda; 3. Inicio Flexión (iF), etapa intermedia del inicio de la flexión de la notocorda; 4. Flexión (F), etapa donde aparece la aleta pélvica y radios de la aleta caudal; se observa claramente la flexión de la notocorda; 5. Inicio Post-Flexión (iPtF) - 6. Post-Flexión (PtF), marcados por la flexión total de la notocorda, formación completa de huesos hipurales y 7. Juvenil (J), etapa donde se observa la formación total de los radios de las aletas, aparición de escamas y absorción total de la aleta primigenia. Todas las larvas se fijaron en formol tamponado al 10% y fueron observadas en un estereomicroscopio SteREO Discovery.V12 (Carl Zeiss Microscopy, LLC) con analizador de imagen.

Análisis morfométrico. Las variables morfométricas para larvas y juveniles fueron: longitud total (TL), longitud estándar (SL), longitud del hocico (Sn.L), diámetro del ojo (ED), longitud de la cabeza (HL), altura de la cabeza (HH), longitud del hocico a la aleta dorsal (Sn.D.F), longitud del hocico a la aleta pectoral (Sn.P.F), longitud del hocico a la aleta pélvica (Sn.Pe.F), longitud del hocico a la aleta anal (Sn.A.F) y altura del cuerpo (HB); la Figura 1-b, muestra específicamente estas medidas para juveniles. Las variables merísticas también abarcaron el número total, el número pre-anal y post-anal de miómeros.

Las medidas corporales fueron transformadas a logaritmo en base 10, ajustando la ecuación alométrica Y=aX^b (^{Fuiman, 1983}), la cual pasó a ser, después de la transformación, a la ecuación Log ₁₀ Y=Log ₁₀ a+bLog ₁₀ X (^{Stewart, 2008}), permitiendo mejorar el coeficiente de determinación (R²), arrojando un modelo más confiable para pronósticos en nuevas observaciones y análisis de máximos y mínimos de las variables relacionadas; siendo Y la variable dependiente, X la variable independiente, a el valor del intercepto en logaritmo y b , el valor de la pendiente o coeficiente alométrico (Fuiman, 1983). Los conceptos de isometría y alometría basados en el valor de la pendiente dan como resultado, si b=1 , Y y X, crecen a la misma velocidad (isometría); si b>1 , Y crece más rápido que X (alometría positiva) y si b<1 , Y crece más despacio que X (alometría negativa) (^{Peña, y Dumas, 2009}). Para evaluar el crecimiento en todo el ciclo experimental, los valores de longitud total y estándar y su relación con respecto al tiempo en HPE se ajustaron al modelo alométrico transformado a logaritmo como fue descrito previamente (Log₁₀Y=Log₁₀a+bLog₁₀X). Para obtener los intervalos de confianza al 95%, se usó el estadístico “t-student”, con n-2 grados de libertad, para la pendiente (coeficiente de alometría) del modelo transformado (^{Walpole, Myers, Myers, y Ye, 2012}). Las relaciones establecidas fueron: longitud total (LT) (X) y longitud estándar (LE) (Y), LE (X) y altura de la cabeza (HCBz)(Y), LE (X) y longitud cabeza (LCBz)(Y), LE(X) y altura del cuerpo (H-Cuerpo) (Y) y LCBz(X) y diámetro del ojo (Diam.OJO)(Y). Cada relación se representó en un diagrama de dispersión, ajustando los datos a la línea recta con base en la ecuación alométrica descrita anteriormente (^{Ricker, 1979}). El desempeño en crecimiento para la especie fue evaluado con base en la Tasa de Crecimiento Instantáneo (TCI), de acuerdo con la fórmula TCI=(LnLT ₂ -LnL ₁ )/t ₂ -t ₁ ; siendo LnLT₂ y LnLT₁ los logaritmos naturales de los promedios de las longitudes totales en sus respectivos tiempos (t₂ y t₁) (Ricker, 1979). Todos los datos se analizaron con el paquete estadístico Stat-Graphics Centurion XV, con licencia de la Corporación Universitaria Lasallista.

RESULTADOS

Periodo embrionario. Presenta una duración de 13 h que van desde la fertilización hasta la eclosión. En 190 huevos se midió diámetro total del huevo (ED), diámetro del vitelo (YD) y espacio perivitelino (PS), desde los 90 s hasta la 4ª hora después de la fertilización. ED obtuvo un promedio de 2,72 ± 0,12 mm (± desviación estándar), con un mínimo de 2,077 mm y un máximo de 3,21 mm; YD, presentó una longitud promedio de 1,44 ± 0,05 mm, con un valor mínimo de 0,76 mm y un máximo de 1,80 mm; PS, con relación al volumen total del huevo, fue “restricto” (3,73%) y “moderado” con respecto a ED (15,4%), con un promedio de 0,42 ± 0,05 mm, un valor mínimo de 0,12 mm y un valor máximo de 0,72 mm. El embrión midió 2,24 ± 0,21 mm y su saco vitelino 0,51± 0,1 mm³. La Tabla 2 resume los valores encontrados en esta etapa.

Periodo larval. La larva eclosiona a las 12 h post-fertilización; este periodo va desde las 0 h post-eclosión (HPE) hasta la 254 HPE, que corresponde desde larval vitelino (LV) hasta Post-Flexión (PtF). La Figura 2, describe las etapas de desarrollo larval, comprendidas entre las 0 HPE hasta los 15 días post-eclosión (DPE); la Tabla 2 indica todas las medidas y sus relaciones desde LV hasta juvenil (J).

Estadio larval vitelino (LV) (Figuras 2A y 2B-Tabla 2). Periodo de tiempo comprendido entre la hora 13 post-fertilización (HPF), hasta las 24 HPE. Las larvas alcanzaron valores promedios en TL de 4,75 mm; SL de 4,21 mm; ED de 0,37 mm; BH de 1,26, HL de 1,02 mm y HH de 0,88 mm. La relación ED/HL del 34,89 ± 0,63%, indica un ojo grande; HL/SL del 24,06 ± 0,37% y BH/SL del 30,91 ± 0,26, indican respectivamente una longitud de cabeza y una longitud del cuerpo moderados. Al eclosionar la larva ya presenta un ojo en su mayoría pigmentado y algunos melanocitos dendríticos en la parte anterior del saco vitelino, patrón de pigmentación que incrementa durante toda la fase larval. Hacia la 8 HPE, el brote de la aleta pectoral comienza a girar hacia la parte posterior. El peso promedio en esta etapa fue 1,8 ± 0,2 mg.

Tabla 2 Mínimo (MIN) y máximo (MAX), promedios (X) y desviación estándar (SD) para las variables morfométricas y merísticas en el periodo embrionario, larvas y juveniles de Brycon moorei (adaptada de ^{Nakatani et al, 2001}) 

n: número de individuos, ED: diámetro total, YD: diámetro del vitelo, PS: espacio perivitelino, LV: larval vitelino, PF: pre-flexión, iF: inicio flexión, F: flexión, iPtF: inicio post-flexión, PtF: post-flexión, J: juvenil, TL: longitud total, SL: longitud estándar, SnL: longitud hocico, ED: diámetro ojo, HL: longitud cabeza, HH: altura cabeza, BH: altura cuerpo, CF: aleta caudaL, DF: aleta dorsal, PF: aleta pectoral, PeF: aleta pélvica, AF: aleta anal, SnPF: hocico-aleta pectoral, SnPeF: hocico-aleta pélvica, SnDF: hocico-aleta dorsal, SnAF: hocico-aleta anal, YL: longitud del vitelo, YH: altura del vitelo, AAU: aleta ausente, NV: no visible y S.D.: sin dato.

Figura 2 Etapas de desarrollo larval de Brycon moorei. A. 0-HPE-LV; B. 24-HPE-LV; C. 25-53 HPE-PF; D. 56-71 HPE-iF; E. 82-110 HPE-F; F. 132-182 HPE-iPtF; G. 206-254 HPE-PtF; H. 278-326 HPE-J. 

Estadio pre-flexión (PF) (Figura 2C-Tabla 2). Periodo de tiempo comprendido entre las 25 hasta las 53 HPE. Las larvas alcanzaron valores promedios en TL de 6,85 mm, incremento que corresponde al 44% con relación a la etapa anterior; SL de 6,32 mm; ED de 0,39 mm; BH de 1,19, el cual redujo su crecimiento en un 5,6% debido a la reducción del vitelo; HL de 1,29 mm y HH de 1,18 mm. La relación ED/HL del 30,41 ± 0,18%, indica un ojo moderado; HL/SL del 20,3 ± 0,74% y BH/SL del 18,95 ± 0,79, indican respectivamente una longitud de cabeza moderada y una longitud del cuerpo largo. Aún se observan vestigios de la aleta primigenia; la vejiga gaseosa a las 44 HPE se ve inflada, transparente y externamente pigmentada en su región dorsal; los patrones de pigmentación en general incrementan en la región ventral y cefálica; el opérculo cubre toda la abertura branquial y llega hasta el brote de la aleta pectoral.

Estadio inicio de flexión (iF) (Figura 2D-Tabla 2). Periodo de tiempo comprendido entre las 56 hasta las 71HPE. Las larvas alcanzaron valores promedios en TL de 7,89 mm, incremento que corresponde al 15% con relación a la etapa anterior; SL de 7,35 mm; ED de 0,46 mm; BH de 1,45 mm, el cual incrementó en un 21,8% debido al consumo de alimento exógeno; HL de 1,32 mm y HH de 1,37 mm. La relación ED/HL del 35,67 ± 0,36%, indica un ojo grande; HL/SL del 17,85 ± 0,3% y BH/SL del 19,74 ± 0,6, evidencian, respectivamente, una longitud de cabeza pequeña y una longitud del cuerpo largo. Se conservan las mismas características morfológicas de la etapa anterior, con ligera pigmentación puntiforme a nivel caudal.

Estadio de flexión (F) (Figura 2E-Tabla 2). Periodo de tiempo comprendido entre las 82 hasta las 110 HPE. Las larvas alcanzaron valores promedios en TL de 9,85 mm, incremento que corresponde al 24,8% con relación a la etapa anterior; SL de 9,26 mm; ED de 0,66 mm; BH de 2,12 mm, el cual aumentó en un 60,6% debido al consumo de alimento exógeno; HL de 1,96 mm y HH de 1,82 mm. La relación ED/HL del 34,2 ± 0,24%, indica un ojo grande; HL/SL del 21,14 ± 0,36% y BH/SL del 22,94 ± 0,31, indican respectivamente una longitud de cabeza y una longitud del cuerpo moderados. La notocorda alcanza la flexión total a las 95 HPE, siendo visible en el pedúnculo caudal con incremento notorio de la pigmentación dendrítica y puntiforme rodeando la notocorda y los radios de la aleta caudal primigenia; además, las aletas dorsal y anal comienzan su desarrollo y se observan también las aberturas nasales de forma alargada.

Estadio de inicio de post-flexión (iPtF) (Figura 2F-Tabla 2). Periodo de tiempo comprendido entre las 123 hasta las 182 HPE. Las larvas alcanzaron valores promedios en TL de 13,81 mm, incremento que corresponde al 40,3% con relación a la etapa anterior; SL de 12,44 mm; ED de 0,96 mm; BH de 3,13 mm, el cual incrementó en un 47,64%, aunque menor comparado al incremento del 60,6% de la etapa de flexión, probablemente debido a la transición a la dieta seca; HL de 3,24 mm, el cual incrementó un 66,8% con respecto a la etapa anterior y HH de 2,65 mm. La relación ED/HL del 30,55 ± 0,21%, indica un ojo moderado; HL/SL del 25,76 ± 0,38%, y BH/SL del 24,91 ± 0,4, indican respectivamente una longitud de cabeza y una longitud del cuerpo moderados. La aleta primigenia ha sido absorbida en gran mayoría, solo se observa una tira delgada en el segmento preanal. Se observa un ensanchamiento en la base del pedúnculo caudal que formará la aleta adiposa. Se observan también pliegues y estrías a nivel del sistema digestivo y aletas anal, dorsal y caudal con radios. Los miómeros ya no son visibles a partir de esta etapa.

Estadio de post-flexión (PtF) (Figura 2G-Tabla 2). Periodo de tiempo comprendido entre las 206 hasta las 254 HPE. Las larvas alcanzaron valores promedios en TL de 15,82 mm, incremento que corresponde al 14,55% con relación a la etapa anterior; SL de 13,68 mm; ED de 1,1 mm; BH de 3,6 mm, el cual incremento en un 15,01% con respecto a la etapa anterior, menor probablemente debido a la transición a dieta seca; HL de 3,47 mm y HH de 3,29 mm. La aleta caudal, tienen entre 24 a 26 radios; la aleta dorsal entre 9 a 10 radios, la aleta anal entre 20 a 22 radios; los radios de las aletas pectorales y pélvicas son difíciles de observar; con respecto a la pigmentación de las aletas, la única que no presentó patrones de pigmentación fue la aleta pectoral, independiente de su grado de desarrollo. La relación ED/HL del 29,91 ± 0,087%, indica un ojo moderado; HL/SL del 27,28 ± 0,43%, y BH/SL del 28,62 ± 0,57, indican respectivamente una longitud de cabeza y una longitud del cuerpo moderados. Se mantienen la mayoría de las características morfológicas descritas en la etapa anterior.

Estadio de juvenil (J) (Figura 2H-Tabla 2). Periodo de tiempo comprendido entre las 278 hasta las 326 HPE. Los juveniles alcanzaron valores promedios en TL de 25,66 mm, incremento que corresponde al 62,2% con relación a la etapa anterior; SL de 21,49 mm; ED de 1,71 mm; BH de 6,13 mm, el cual incremento en un 70,2% con respecto a la etapa anterior, mayor probablemente debido a la adaptación a la dieta seca; HL de 6,33 mm y HH de 4,89 mm. La aleta caudal, tienen entre 25 a 27 radios; la aleta dorsal entre 11 a 13 radios, la aleta pectoral entre 15 a 16 radios, la aleta anal entre 29 a 32 radios y la pélvica entre 9 a 10. La relación ED/HL del 27,17 ± 0,071%, indica un ojo moderado; HL/SL del 29,45 ± 0,095%, y BH/SL del 28,56 ± 0,082, indican respectivamente una longitud de cabeza y una longitud del cuerpo moderados. Se observan escamas y todas las estructuras natatorias bien desarrolladas, incluyendo la aleta adiposa. Se observa además el punto negro característico en la base del pedúnculo caudal, con proyecciones rojizas sobre la aleta caudal e iridiscencias en todo el cuerpo, típicas en esta especie.

Figura 3 Relación lineal en logaritmo para LT y LE, con respecto a las HPE, durante el ciclo de desarrollo larval a juvenil en B. moorei (n=470). Larval Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación. * indica el inicio de la dieta seca. 

Relaciones corporales lineales. Los modelos de crecimiento de las diferentes variables corporales relacionadas para B. moorei, se presentan desde la Figura 3 a la 8; en ellas se pueden observar las relaciones lineales entre longitud total (LT) y longitud estándar (LE) vs horas post-eclosión (HPE); LT vs LE; LE vs longitud de la cabeza (LCBz); LE vs altura de la cabeza (HCBz); LE vs altura del cuerpo (H-Cuerpo) y LCBz vs diámetro del ojo (DIAM.Ojo). Además, se incluyen los intervalos de confianza (95%), que determinan las relaciones de isometría o alometría positiva o negativa y si existen diferencias entre los rangos para todas las etapas de desarrollo temprano que incluye desde larva hasta juvenil. La Figura 3 muestra la curva de crecimiento en LT y LE, con respecto al tiempo desde larva vitelina hasta juvenil, explicado a través del modelo logarítmico, con un coeficiente de determinación del 97% para LT y del 98% para LE, con una alta correlación entre las variables (98%). Se indica también un periodo en el que el crecimiento fue menor (144 a 250 HPE), el cual coincide con el inicio del protocolo de alimentación con dieta seca (* 144 HPE).

La Figura 4 indica la relación entre la longitud total y la longitud estándar, con una correlación superior al 90% y un coeficiente de determinación elevado en todas las etapas; existe isometría (b=1) en LV, iF y F, es decir que LT y LE crecen en la misma proporción; para PF, iPtF, PtF y J la alometría es menor que 1, por lo tanto la LE no crece en la misma proporción que la LT; los intervalos de confianza permitieron establecer que la pendiente ( b ) de LV tiene diferencias significativas con PF, iPtF y J, a su vez la b de PF es diferente significativamente de F y b de F es significativamente diferente de J.

Figura 4 Relaciones lineares en logaritmo e intervalos de confianza (95%) para LT y LE, durante las etapas Larval Vitelino (n=180), Pre-Flexión (n=90), inicio Flexión (n=50), Flexión (n=50), inicio Post-Flexión (n=40), Post-Flexión (n=30) y Juvenil (n=30). Larva Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación y r: Coeficiente de correlación. 

La Figura 5 indica la relación entre LE vs LCBz, en la cual para las etapas LV, PF, iF e iPtF, la longitud de la cabeza crece de manera más rápida (b>1) que la longitud estándar, a su vez los coeficientes de correlación y determinación son elevados, mostrando alta correlación entre las variables y un modelo confiable. Para F y PtF existe isometría (b=1) entre las variables, con elevada correlación y un modelo ajustado; para J, aunque el crecimiento fue alométrico negativo, no se encontró relación entre las variables (r=24%) y el modelo no logra explicar dicha relación (R²=6%) por la elevada dispersión de los datos. Los intervalos de confianza permitieron establecer que b de LV tiene diferencias significativas con F; PtF y J y que b de J tiene diferencias estadísticas significativas con todas las pendientes de las otras etapas.

La Figura 6 muestra la relación entre la longitud estándar y la altura de la cabeza, donde LV e IPtF, presentan un crecimiento alométrico positivo (b>1), con un modelo confiable y una alta correlación entre las variables; para PF, iF y F, existe isometría, con un R² y r elevados a excepción de PF, donde existió una mayor dispersión de los datos. Para PtF, el modelo no se ajustó por la elevada dispersión de los datos y en J, existe alometría negativa, con elevada correlación y confiabilidad del modelo.

Figura 5 Relaciones lineares en logaritmo e intervalos de confianza (95%) para LE y LCBz, durante las etapas Larval Vitelino (n=180), Pre-Flexión (n=90), inicio Flexión (n=50), Flexión (n=50), inicio Post-Flexión (n=40), Post-Flexión (n=30) y Juvenil (n=30). Larva Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación y r: Coeficiente de correlación. 

Figura 6 Relaciones lineares en logaritmo e intervalos de confianza (95%) para LE y HCBz, durante las etapas Larval Vitelino (n=180), Pre-Flexión (n=90), inicio Flexión (n=50), Flexión (n=50), inicio Post-Flexión (n=40), Post-Flexión (n=30) y Juvenil (n=30). Larva Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación y r : Coeficiente de correlación. 

La Figura 7, representa la relación entre LE y la altura del cuerpo, la cual en la primera etapa (LV) muestra de manera clara el efecto del gasto de las reservas endógenas (vitelo), por lo cual la correlación es negativa, es decir entre la larva más crece más delgada se vuelve; luego se presenta isometría para F, con elevada correlación; para PF y J también se presenta isometría, solo que la dispersión de los datos es alta, en este caso la correlación solo llega al 56%; iF e iPtF, con elevada correlación y modelo ajustado la alometría es positiva.

Figura 7 Relaciones lineares en logaritmo e intervalos de confianza (95%) para LE y H-Cuerpo, durante las etapas Larval Vitelino (n=180), Pre-Flexión (n=90), inicio Flexión (n=50), Flexión (n=50), inicio Post-Flexión (n=40), Post-Flexión (n=30) y Juvenil (n=30). Larva Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación y r : Coeficiente de correlación. 

Por último, la Figura 8, indica la relación entre la longitud de la cabeza y el diámetro del ojo, donde LV, iF, iPtF, PtF y J, presentan elevada correlación y un modelo ajustado, mostrando alometría negativa; F, también presentó alometría negativa pero el modelo no es confiable, además presenta elevada dispersión de los datos, por lo tanto, la correlación entre las variables es baja; se presentó isometría en PF, pero con baja correlación y elevada dispersión de los datos. Los intervalos de confianza permitieron establecer que b de LV tiene diferencias significativas con todas las etapas, en las relaciones indicadas por las Figuras 5, 6 y 7; de manera general los rangos de los intervalos de confianza que más diferencias estadísticas tuvieron para todas las relaciones fueron los de la etapa juvenil.

Figura 8 Relaciones lineares en logaritmo e intervalos de confianza (95%) para LCBz y Diam.OJO, durante las etapas Larval Vitelino (n=180), Pre-Flexión (n=90), inicio Flexión (n=50), Flexión (n=50), inicio Post-Flexión (n=40), Post-Flexión (n=30) y Juvenil (n=30). Larva Vitelino (LV), Pre-Flexión (PF), Inicio Flexión (iF), Flexión (F), Inicio de Post-Flexión (iPtF), Post-Flexión (PtF) y Juvenil (J). (b<1); (b>1) y (b=1): alometría negativa, positiva o isometría. R²: Coeficiente de determinación y r : Coeficiente de correlación. 

Los valores de TCI para LT examinados en el transcurso del desarrollo de B. moorei, son mostrados en la Figura 9.

Figura 9 Tasa de crecimiento instantáneo (TCI) para LT durante el desarrollo larval de Brycon moorei. 

En las primeras 77 HPE, se observaron las mayores variaciones en la TCI para B. moorei, incluyendo los valores más altos entre la eclosión hasta la 3 HPE, con una caída hasta las 12 HPE y luego con una serie de altos y bajos; el valor más bajo se encontró a las 63 HPE, siguiendo con una ligera estabilidad, para luego, entre las 133 hasta las 181 HPE tener otra caída.

DISCUSIÓN

Los cambios en el desarrollo temprano, son inherentes a cada especie (^{Balon, 1981}) y a su vez, estos cambios están muy relacionados con las condiciones ambientales prevalentes. La temperatura, por ejemplo, puede acelerar o retrasar el desarrollo embrionario en los peces (^{Nakaghi, Neumann, Faustino, Mendes, y De Braga, 2014}; ^{Osse, y Van den Boogaart 1995}). Para B. moorei la temperatura, entre otros, se mantuvo dentro de los rangos de confort para la especie (^{David-Ruales, y Castañeda-Álvarez, 2014}). Muchas estructuras se convierten en claves para determinar parámetros de clasificación, como el tamaño del huevo (diámetro), el espacio perivitelino y el brote de la aleta pectoral, entre otras (^{Nakatani et al, 2001}). Brycon moorei presenta un huevo pequeño (2,72 ± 0,12 mm); el diámetro del huevo se puede relacionar con la estrategia reproductiva, por lo general especies que no presentan cuidado parental presentan este tipo de huevo (^{Kolm, y Ahnesjö, 2005}), otros autores afirman que muchos de los peces teleósteos tienen esta característica, donde las reservas de vitelo se consumen rápidamente y las larvas producto de estas estructuras son capaces de alimentarse de presas pequeñas y adaptarse al medio que los rodea (^{Osse, 1990}; ^{Osse et al, 1997}).

Dentro del género Brycon al eclosionar, ciertas características anatómicas se comparten, por ejemplo, ausencia de pigmentación, ano, intestino y vejiga gaseosa, además, la vesícula óptica está indiferenciada y la aleta embrionaria es hialina y rodea todo el cuerpo (^{Faustino, Makino, Neumann, y Nakaghi, 2018}; ^{Marques et al, 2017}). Estos rasgos se observaron en el presente estudio y también en larvas de: B. insignis (piabanha) (^{Andrade-Talmelli, Kavamoto, Romagosa, y Fenerich-Verani, 2001}). B. orbignyanus (piracanjuba) (^{Bettinelli-Nogueira, Lima-Godinho, y Pereira-Godinho, 2014}; ^{Raposo-Maciel et al, 2010}; ^{Reynalte-Tataje, Zaniboni-Filho, y Esquivel, 2004}); B. amazonicus (yamú) (^{Da Silva et al, 2017}; ^{Mira-López, Medina-Robles, Velasco-Santamaría, y Cruz-Casallas, 2007}; ^{Neumann, Paes, Mendes, Braga, y Nakaghi, 2018}); B. cephalus (sinonimia=B. amazonicus) (matrinxã) (^{Romagosa, Narahara, y Fenerich-Verani, 2001}; ^{Sversut de Alexandre et al, 2010}); B. gouldingi (matrinchã) (^{Faustino, Nakaghi, y Neumann, 2011}; ^{Faustino, Makino, Neumann, y Nakaghi, 2015}) y B. nattereri (peripetinga) (^{Maria, Ninhaus-Silveira, Orfão, y Viveiros, 2017}). A diferencia de otros bricónidos (^{Andrade-Talmelli et al, 2001}; ^{Mira-López et al, 2007}; ^{Nakatani et al, 2001}; ^{Reynalte-Tataje et al, 2004}), la pigmentación del ojo en B. moorei inició más rápido (11 h después de la fertilización), no después de la eclosión. Se conoce que la visión es fundamental en la percepción de la distancia para la obtención de alimento y evitar la depredación (^{Hubbs, y Blaxter, 1986}), además el género Brycon se caracteriza por tener ojos grandes y pigmentados, que indican mayor facilidad para la caza y evitar la depredación (^{Ceccarelli, 1997}). En B. moorei, los patrones de pigmentación corporal, incrementan con el cambio de fase de crecimiento, intensificándose desde el inicio de la flexión, las máculas típicas se presentan a partir de la etapa de post-flexión localizadas en el pedúnculo caudal y en la parte posterior del opérculo; similares características son descritas para B. amazonicus (^{Gomes, y Urbinati, 2005}).

En teleósteos existen patrones secuenciales de crecimiento (alometría) (^{Osse, y van den Boogaart, 1995}), estos patrones se presentaron en B. moorei con diferencias dentro de cada etapa de desarrollo. Tales diferencias pueden estar relacionadas con la competencia por el alimento (^{Osse et al, 1997}; ^{van Snik, Van den Boogaart, y Osse, 1997}), por el canibalismo exacerbado que presenta esta especie (^{Baras, 2005}; ^{Baras, y Jobling, 2002}; ^{Baras, y Lucas, 2010}; ^{Baras, Maxi, Ndao, y Mélard, 2000a}; ^{Baras et al, 2000b};) y por factores medioambientales (^{Fuiman, y Higgs, 1997}; ^{Fuiman, Poling, y Higgs, 1998}).

El crecimiento larval por lo general sigue el modelo exponencial (^{Fuiman, 1983}). En el presente estudio, los patrones de crecimiento alométrico fueron explicados a través de la transformación de este modelo, el valor de LT para B. moorei al final del periodo experimental fue de 25,66 ± 0,64 mm; usando el mismo modelo, se reportó un valor similar en B. orbignyanus (24,9 mm) (^{Bettinelli-Nogueira et al, 2014}); para la etapa de post-flexión en el presente estudio se encontró un valor de 15,82 ± 0,41 mm, se reportan valores similares en la misma etapa de desarrollo para larvas de B. insignis con 16,48 ± 1,38 mm, larvas de Salminus brasiliensis (14,84 ± 1,70 mm) (^{Souza et al, 2017}) y P. vimboides con 15,68 ± 1,0 mm (^{Souza et al, 2016}) y, sin reportar el modelo, menor valor para B. amazonicus, con una LT de 18,89 ± 3,75 (^{Neumann et al, 2018}); menor valor para B. orthotaenia al día 7 con 9,58 ± 0,67 mm (^{Zeferino-Gomes, Sato, Rizzo, y Bazzoli, 2011}), siendo que en B. moorei en el mismo periodo alcanzaba los 13,66 ± 1,27 mm.

Con respecto al crecimiento en el periodo experimental para B. moorei reflejado en la relación de la LT y LE versus HPE, se observa entre las 150 y las 250 HPE una queda en el crecimiento, comparado con las etapas anteriores y la posterior a las 250 HPE, este patrón en el crecimiento puede estar influenciado por el cambio en el protocolo de alimentación, el cual coincide con el inicio de la dieta; se conoce que el crecimiento y la composición química corporal de los peces pueden ser influenciados por la manipulación de los nutrientes de las dietas (^{Dabrowski, 1984}), además en larvas de desarrollo indirecto la transición a dieta seca, después del alimento vivo, puede presentar este tipo de comportamientos, por lo que se recomienda el diseño de protocolos de alimentación (^{Portella, Leitão, Takata, y Lopes, 2012}), asociado al seguimiento de la calidad del agua, el tipo de contendor y aspectos etológicos de la especie (^{David-Ruales, y Castañeda-Álvarez, 2010})

Respecto a la relación LT vs LE, en B. moorei, se encontraron variaciones entre isometría y alometría negativa con rangos muy parecidos, este resultado puede ser el reflejo del crecimiento de la aleta caudal; al respecto, ^{Osse et al (1997)} indican que el desarrollo de esta estructura está relacionado con las respuestas a la captura de presas y la huida de depredadores. La correspondencia entre la LCBz y HCBz vs LE, tubo en principio alometría positiva, con acelerado crecimiento hasta PtF, indicando priorización y diferenciación de estructuras bucales, sensoriales y otras que tienen que ver con el consumo de alimento y el cambio de respiración cutánea a branquial (funciones vitales) (^{Fuiman, 1983}; 2002; ^{Fuiman et al, 1998}; ^{Osse, 1990}; ^{Osse, y van den Boogaart, 1995}; ^{Osse et al, 1997}; ^{Van Snik et al, 1997}; ^{Vandewalle, Germeau, Besancenet, Parmentier, y Baras, 2005}). Este tipo de desarrollo también se presentó en B. amazonicus (^{Neumann et al, 2018}; ^{Sampaio-Nakauth, Villacorta-Correa, Figueiredo, Bernardino, y França, 2016}), B. cephalus (^{Romagosa et al, 2001}), B. orbignyanus (^{Bettinelli-Nogueira et al, 2014}) y en otras especies como Prochilodus vimboides (^{Souza et al, 2016}), Pseudoplatystoma reticulatum (^{Andrade et al, 2016}), Paralichthys californicus (^{Gisbert et al, 2002}), Acipenser medirostris (^{Gisbert, y Doroshov, 2006}), Ciprinus carpio, Clarias gariepinus (^{Van Snik et al, 1997}) y para Nannacara anomala (^{Kupren, Prusi, Arski, Krejszeff, y Kucharczyk, 2014}). En general, para la relación LCBz vs DIAM.OJO hubo una marcada alometría negativa, coincidiendo con la especie P. vimboides (^{Souza et al, 2016}) y también para la especie P. reticulatum cuando comparado con la LT. Al respecto, se conoce que los cambios morfológicos se asocian con variaciones alométricas en órganos y otras partes del cuerpo con relación a las presiones medioambientales; sin embargo, órganos como el ojo y el cerebro crecen a tasas asociadas más con la edad que con el tamaño del cuerpo. Esto implica provisión de habilidades sensoriales y cognitivas necesarias para la sobrevivencia (^{Gonzalez-Voyer, Winberg, y Kolm, 2009}; ^{Pankhurst, y Montgomery, 1994}). La interacción entre HCuerpo vs LE, pasó de alometría negativa en la etapa LV, alcanzando la isometría en la etapa de PF, para luego ser predominante la alometría positiva. El resultado en LV es directamente relacionado con el alto consumo de vitelo, el cual disminuye el tamaño del saco vitelino, influenciando directamente sobre la altura del cuerpo, resultados iguales se observaron en las especies: P. vimboides (^{Souza et al, 2016}); C. carpio y C. gariepinus (^{Van Snik et al, 1997}).

Con respecto a la TCI, el mayor valor al inicio de la eclosión se deba al hecho de presentar todas las reservas endógenas, con la correspondiente caída, dada al consumo de las mismas debido a la elevado metabolismo por la formación de estructuras relacionadas con la locomoción y la alimentación (^{Osse et al, 1997}). Además del estrés generado por el ambiente, el protocolo de alimentación y hacia las 22 h post-eclosión el incremento exacerbado del canibalismo, que se extiende durante todo el proceso. Las tasas de crecimiento no vuelven a su estado inicial debido al gasto energético por la actividad de caza y/o depredación. La otra caída que se observa hacia las 180 HPE, puede estar asociada probablemente con la transición alimenticia a dieta seca, en adelante se observa una ligera estabilidad hasta el final del experimento.

Los valores morfológicos reportados para B. moorei (Tabla 2), son similares a otros especímenes de la familia Characidae y también del género Brycon por ejemplo: B. orbignyanus (^{Nakatani et al, 2001}), S. brasiliensis (dorado), Leporinus steindachneri (piau), P. lineatus (curimba) y B. insignis (piabanha) (^{Souza et al, 2017}), P. vimboides (curimbatá) (^{Souza et al, 2016}) y mayores al ser comparados con las especies B. hilarii (piraputanga) (^{Oliveira et al, 2012}), Pseudoplatystoma reticulatum (bagre) (^{Andrade et al, 2016}) y Pimelodus britskii (bagre) (^{Almeida et al, 2018}).

El modelo de crecimiento y las relaciones morfométricas encontradas para B. moorei, sugiere la priorización de estructuras que le dan a esta especie la habilidad para la locomoción y la captura de alimento en respuesta a ser considerada la especie de agua dulce con mayor canibalismo reportado para larvas de peces de agua dulce (^{Baras, y Jobling, 2002}; ^{Baras y Lucas, 2010}; ^{Baras et al, 2000b}). Esta priorización es determinante para la sobrevivencia de los estadios iniciales en peces (^{Osse, 1990}; ^{Osse et al, 1997}). Se conoce de la dificultad que se presenta para determinar el desarrollo temprano de las especies del neotrópico, la poca información, escasez de claves taxonómicas y el gran parecido entre especies en estos estados, constituye uno de los principales cuellos de botella a la hora de realizar este tipo de estudios (^{Lopes, Oliveira, Bialetzki, y Agostinho, 2015}). Con los resultados obtenidos, se puede indicar que el manejo en confinamiento de la especie B. moorei, necesita de un protocolo de alimentación que incluya alimento vivo de diferentes tamaños y, que la transición a dieta seca, debe hacerse de manera paulatina, sin dejar en principio la co-alimentación con alimento vivo, de acuerdo con los patrones de crecimiento presentados; este procedimiento permite el mejor manejo para una especie de alto canibalismo, asegurando mejores tasa de sobrevivencia, que permiten certificar su adaptación para la acuicultura. Además, las relaciones morfométricas presentadas, se constituyen en valiosa información que servirá como claves merísticas para los estudios de ictioplancton, toda vez que la especie está declarada como vulnerable.