Objetivo. Caracterizar las curvas de crecimiento y el comportamiento reproductivo de vacas Romosinuano en México. Materiales y métodos. Se analizaron 928 registros de peso y edad de 127 hembras para características de crecimiento, y 459 registros de partos de 113 vacas para características reproductivas. Se ajustaron cuatro modelos no lineales mixtos para describir la curva de crecimiento, y se estimaron los indicadores de edad al 50% de madurez (E50M) y porcentaje de madurez a los 18 meses de edad (PM18). Para evaluar el comportamiento reproductivo se estimó la edad al primer parto (EPP), intervalo entre partos (IEP), y tasa de pariciones (TP). Resultados. El modelo de Brody fue el que mejor ajustó la curva de crecimiento. El parámetro A, que estima el peso maduro, osciló entre 430.2 y 466.2 kg; el . entre 0.94 y 5.45; y el . entre 0.001 y 0.002; dependiendo del modelo utilizado. La E50M varió entre 20 y 21 meses, y el PM18 entre 44.5 y 45.5%. Las medias de EPP, IEP, y TP a los 550 d fueron 1059±282 d, 553±68 d, y 74±19%, y mejorando (p<0.05) cuando las hembras presentan mayor porcentaje de su peso maduro, dependiendo de la época de parto. Conclusiones. El peso maduro estimado indica que el ganado Romosinuano puede considerarse como una raza de tamaño pequeño o mediano, con aceptable aptitud reproductiva. Las hembras con menor tamaño maduro estimado presentaron edades al primer parto más temprano que vacas más pesadas.

La producción del ganado bovino en las regiones tropicales cálidas se caracteriza por tener un bajo nivel técnico de manejo y un sistema de producción tradicional basado en una baja calidad de praderas, escaso uso de tecnología y limitados recursos genéticos adaptados al trópico (1). La producción bovina en zonas tropicales está afectada por varios factores, siendo los climáticos los principales que afectan el crecimiento animal (2). Los bovinos criollos son animales que surgen del aislamiento y adaptación de ganado desde su llegada al continente americano, originando diferentes razas que se desarrollaron en ambientes hostiles de América (3). Los bovinos Romosinuano son animales criollos adaptados a condiciones de clima adverso, lo cual sugiere que su crecimiento y comportamiento reproductivo es aceptable considerando las contingencias que conllevan los ambientes tropicales.

La evaluación del crecimiento de los animales es una práctica importante para evaluar la rentabilidad y la eficiencia en la producción de carne. El crecimiento corporal se ha representado mediante modelos no lineales, que describen el crecimiento a lo largo del tiempo en función de un número reducido de parámetros (4). Por otro lado, un factor que restringe la productividad de la ganadería tropical es la baja eficiencia reproductiva. Ossa et al (5) mencionan que una característica sobresaliente del ganado Romosinuano es su adaptación a condiciones tropicales, lo cual se refleja en excelentes índices de fertilidad.

Desde hace dos décadas, productores mexicanos establecieron hatos de bovinos Romosinuano (3), cuya población ha ido en aumento (6), por lo que es conveniente la evaluación productiva y reproductiva de este recurso genético. El objetivo de este estudio fue caracterizar el crecimiento y la actividad reproductiva de hembras Romosinuano en México.

Origen de la información. Los datos provinieron de la Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical, particularmente de una unidad de producción localizada en el municipio de Tihuatlán, Veracruz. El rancho está ubicado en un área tropical de México (Noreste del país), en donde los animales eran alimentados con base en pastoreo, principalmente Brachiaria brizantha y Cynodon plectostachyus. El sitio presenta clima cálido subhúmedo con lluvias en verano, temperatura media anual entre 24 y 26°C, y un rango de precipitación de 1100 a 1300 mm (7).

Se estimaron las curvas de crecimiento de hembras Romosinuano, mediante información que provino de 127 animales nacidos entre 2001 y 2018, considerando un total de 928 registros de peso vivo y edad. Los registros de peso vivo y edad se realizaban a lo largo de la vida de los animales: al momento del nacimiento y destete de las becerras (aproximadamente a los ocho meses de edad) (8), y durante las palpaciones para verificar gestación posterior al empadre y previo a cada parto. Las hembras consideradas tuvieron en promedio 7.3 registros por animal y al menos tres registros.

La información para realizar los análisis reproductivos provino de 113 vacas con registros de parto entre 2004 y 2019, considerando un total de 459 registros. Las variables reproductivas consideradas fueron la edad al primer parto (EPP), el intervalo entre partos (IEP) y la tasa de pariciones (TP), y en los análisis finales sólo se consideraron registros de animales dentro de ±3σ de la media. Para el análisis de la EPP se utilizaron registros de 113 vacas, agrupadas en 31 grupos contemporáneos (año y estación de parto). Se definieron dos estaciones de parto, lluvias (junio a noviembre) y secas (diciembre a mayo). Para los análisis de IEP y TP, se utilizaron 346 registros de partos en 81 vacas, agrupados en 28 grupos contemporáneos (año-estación de parto), con lo que se tuvo un promedio de 4.3 partos por vaca. Los números de registros para calcular los intervalos entre partos fueron 74 del primero al segundo, 70 del segundo al tercero, 53 del tercero al cuarto, 48 del cuarto al quinto, y 101 intervalos con vacas de quinto parto o más.

Análisis de crecimiento. Los análisis estadísticos para ajustar las curvas de crecimiento fueron realizados mediante los modelos no lineales mixtos siguientes:

Brody: yij=f(tij;φ.) +εij=A.(1-be-kit )+εij

Logístico: yij=f(tij;φ.)+εij=A./(1+be-kit)+ εij

Von Bertalanffy: yij=f(tij;φ.)+εij=A.(1-be-kit).+εij

Gompertz: yij=f(tij;φ.)+εij=A.(e-be-kit )+εij

Las suposiciones en los modelos fueron: φ.=φ.+η., η.~N(0,Ω), εij~N(0,σ..), donde . es el peso corporal del animal en kilogramos, . es la edad del animal en días, . es la base del logaritmo natural, . es el peso maduro o asintótico en kilogramos, . es una constante de integración, . es la tasa de madurez del animal, φ. y η. son los vectores de efectos fijos y aleatorios respectivos en los parámetros . y k, Ω es la matriz de (co)varianzas de los efectos aleatorios, εij es el residual, y σ.. es la varianza del residual.

Para definir qué modelo de crecimiento utilizar en el cálculo de indicadores del crecimiento, se realizaron pruebas de bondad de ajuste utilizando el criterio de información de Akaike (AIC), el criterio de información Bayesiano (BIC), y el logaritmo de máxima verosimilitud (logL), escogiendo el mejor modelo como aquél con menores valores de AIC y BIC, y con mayor valor de logL. En los cuatro modelos de crecimiento comparados en este estudio, los parámetros relacionados con el peso asintótico (.) y la tasa de madurez (.) se ajustaron como efectos aleatorios, dada la consideración de la parte fija más la aleatoria. Las estimaciones de parámetros y componentes de (co)varianza se estimaron mediante máxima verosimilitud, utilizando el procedimiento nlme (9) en el software R (10). Con los parámetros estimados del mejor modelo, se calcularon los indicadores de crecimiento de edad (EPI) y peso (PPI) al punto de inflexión (con excepción del modelo de Brody), edad al 50% de madurez (E50M) y porcentaje de madurez a los 18 meses de edad (PM18) (Tabla 1).

Análisis reproductivo. Las variables de edad al primer parto (EPP), intervalo entre partos (IEP) y tasa de pariciones (TP) fueron evaluadas como características reproductivas de las vacas. Para el análisis de EPP se consideraron los efectos de grupo contemporáneo (clases formadas con año y estación de parto) y la covariable del peso estimado al parto (generado a partir del modelo de Brody). Para el análisis de IEP se consideraron los efectos incluidos en EPP, además del número de parto como covariable, y los efectos aleatorios de semental y vaca. Para ambas variables de respuesta se evaluó la significancia de los efectos fijos considerados en los modelos (p<0.05), con el propósito principal de ajustarlas por esos efectos ambientales. Los modelos se ajustaron por máxima verosimilitud restringida mediante el paquete de R lme4 (11).

Para estimar la TP, se generó una variable con distribución binaria, codificando como uno cuando la hembra parió, dado que había estado en empadre, y como cero en caso contrario; con el fin de estimar la probabilidad que una vaca esté parida en determinado día después del nacimiento. Para lo anterior, se utilizó un modelo de regresión logística múltiple, usando la función siguiente:

E(y)=1/(1+exp{β0 - x1β1 - x2β2 - x3β3})

Donde β0 es el intercepto de la regresión, β1 es el parámetro relacionado con la época de parto (lluvias y secas), β2 es el parámetro relacionado con el número de parto (primero, segundo, tercero, cuarto y cinco+ partos), y β3 es el parámetro relacionado con la covariable del día después del parto. La estimación de los parámetros se realizó mediante la metodología de máxima verosimilitud, utilizando el paquete glm2 (12) del software R (10).

Para tener puntos de comparación entre diferentes épocas del año y números de partos de las vacas, y para comparar con lo obtenido en estudios similares, se generaron indicadores reproductivos complementarios, estableciendo puntos de corte fijados a momentos específicos. Considerando que el periodo más largo de un animal en este hato para alcanzar el primer parto fue aproximadamente dos años y la expresión de indicadores reproductivos usados para otros estudios, se establecieron los indicadores (usando las probabilidades de vacas paridas a los 550 d y del número de días transcurridos hasta que en el hato ocurría el 50% de los partos): tasa de pariciones a los 550 días (TP550), y días al 50% de partos (D50P).

Análisis de crecimiento. En la tabla 2 se presentan las estimaciones de los parámetros y componentes de (co)varianza para el análisis de crecimiento en los cuatro modelos no lineales mixtos estudiados. Con base en las estimaciones de los parámetros de crecimiento y los componentes de (co)varianza aleatorios, se ajustaron los modelos para predecir y describir el peso de los animales a través del tiempo, para cada uno de los modelos evaluados (Figura 1). En general, los cuatro modelos describieron de forma apropiada las curvas de comportamiento del crecimiento a través del tiempo.

En la tabla 3 se describen los resultados para los indicadores de bondad de ajuste utilizados en la selección del mejor modelo. Con base en los valores de los criterios de información de Akaike, Bayesiano (el menor es mejor) y el logaritmo de máxima verosimilitud (el mayor es mejor), el modelo que presentó mejor ajuste fue el modelo de Brody, seguido de los modelos Logístico y Gompertz, y por último Von Bertalanffy.

Los parámetros estimados sirvieron para calcular indicadores como EPI y PPI de la función, los cuales fueron de 628 d y 215 kg, 309 d y 132 kg, y 418 d y 162 kg para los modelos Logístico, Von Bertalanffy y Gompertz, respectivamente. La E50M y el PM18 estimados para las hembras Romosinuano fueron muy similares en los modelos estudiados, con valores de 634 d y 45.42%, 628 d y 44.54%, 619 d y 45.54%, y 620 d y 45.31%, para los modelos de Brody, Logístico, von Bertalanffy y Gompertz, respectivamente.

Análisis reproductivo. Las medias de EPP e IEP fueron 1059 (±282) y 553 (±68) días. Como era de esperarse, los grupos contemporáneos afectaron a ambas variables (p<0.05). El peso al parto de la vaca afectó la EPP e IEP (p<0.05). Las vacas con menor peso maduro fueron más precoces (1.8 d/kg) y presentaron EPP más temprana. Para IEP, un mayor número de partos significó menores valores para esta variable (-11.23 d por cada parto); mientras que un menor peso estimado al momento del parto significó un mayor IEP (-0.10 d/kg).

Las estimaciones de los parámetros para la regresión logística múltiple en tasa de pariciones resultaron significativas (p<0.01) y los valores se presentan en la tabla 4. Los indicadores reproductivos TP550 y D50P se dividieron para los grupos de época de lluvias y secas, así como para los grupos de número de parto (tabla 5). En general, cuando los partos se presentaron en temporada de lluvias, las hembras mostraron de 5.5 a 12.3% mejor desempeño reproductivo en la variable TP550, dependiendo el número de parto que presentaba la hembra. La variable D50P se redujo 34 días por cada aumento en el número de parto de la hembra, y 32 días si se presentaba el parto anterior en temporada de lluvias.

La TP promedio a los 550 d fue 74±19%. La TP fue afectada (p<0.05) por la época de parto anterior (Figura 2), siendo mayores TP cuando el parto anterior fue en lluvias. Las diferencias de la TP entre diferentes números de parto se muestran en la Figura 3, siendo las TP mayores conforme aumentó el número de parto (p<0.05).

La caracterización de la curva de crecimiento en las hembras Romosinuano usando modelos no lineales, presentó un mejor ajuste cuando se incluyó en su estructura los efectos aleatorios de peso maduro y de la tasa de madurez. Otros estudios relacionados con el ajuste de modelos no lineales en la producción animal, también han presentado mejor ajuste cuando se incluyen efectos aleatorios al modelo, tanto en modelos de curvas de crecimiento (13) como de lactancia (14). Esto indica la importancia del efecto individual de los animales en la expresión de su producción.

Con base en los resultados de las pruebas de bondad de ajuste, el mejor modelo para describir la curva de crecimiento en hembras Romosinuano fue el modelo de Brody. En otros estudios con ganado bovino, el modelo de Brody también ha sido el más apropiado para ajustar su patrón de crecimiento (15,16,17). El modelo de Brody ha sido el más frecuentemente utilizado para caracterizar el crecimiento de ganado bovino, debido, en parte, a que posee una interpretación biológica que facilita la descripción del crecimiento de los animales.

La comparación directa de los valores en los parámetros de crecimiento obtenidos con los diferentes modelos no es posible, dado que miden aspectos ligeramente diferentes; sin embargo, el parámetro para estimar el peso maduro considera una interpretación biológica similar. A pesar de la gran variabilidad entre las estimaciones de EPI y PPI en los modelos estudiados, las estimaciones de peso maduro variaron poco entre los modelos (430 a 466 kg), lo que indica cierta constancia en cercanía del valor estimado con el real, independientemente del modelo utilizado. El valor estimado de ., relativamente pequeño, sugiere que la madurez de las hembras de esta raza podría ubicarse en edades tempranas de la vida de los animales. Los valores estimados en el presente estudio para peso maduro fueron similares a los reportados por Ramírez et al. (18) para seis grupos de hembras bovinas cruzadas (incluyendo el F1 Romosinuano-Cebú) en condiciones tropicales. Los valores de peso maduro estimados concuerdan con la ubicación de Romosinuano como una raza de talla pequeña a mediana. Esta característica en razas criollas y nativas indica que son animales de menores requerimientos nutricionales comparado con razas de tallas mayores, lo que podría representar también ventajas en características reproductivas para sistemas de producción con restricciones alimenticias (19,20). En este estudio, los componentes de varianza de peso maduro oscilaron entre 61.2 y 62.2% de la varianza total, siendo mayor con el modelo de Brody. La varianza de la tasa de madurez no sobrepasó el 0.06% de la varianza total en todos los modelos.

El modelo de Brody, que produjo el mejor ajuste entre los modelos estudiados, no permite estimar la EPI y PPI, debido a que no tiene punto de inflexión. La alta variabilidad en las estimaciones de estos indicadores dependiendo el modelo utilizado, evidencia la necesidad de utilizar la información del mejor modelo. En los tres modelos donde sí se puede calcular la EPI se obtuvieron edades mayores a las de la edad al destete común de la raza en México (240 d) (8), a diferencia de lo observado en la raza Tropicarne, y en ganado Cebú (16,21). Esto sugiere que el destete en los animales Romosinuano no cesa significativamente la aceleración del crecimiento alrededor de esta edad. La E50M estimada entre 20 y 21 meses y el PM18 cerca del 45%, ubican a este recurso genético como tardío en alcanzar la madurez, comparado con otros recursos genéticos tropicales. Domínguez-Viveros et al. (16) estimaron la E50M y el PM18 de bovinos Tropicarne en México aproximadamente al año de edad y 65%, respectivamente. En general, estos valores relativamente bajos indican que los bovinos en el trópico tienen bajas tasas de madurez.

La media estimada para EPP en este estudio (35.3 meses) fue menor en aproximadamente tres meses a las estimaciones reportadas (22,23) para Romosinuano en Colombia, y mayor en tres meses a la reportada para esta raza en México (24). Sin embargo, contrastando con otras razas en México, fue similar al reportado para la raza Tropicarne (25), y de uno a tres meses menor que los reportados por otros autores con ganadería tropical (26,27). Vera (28) menciona que se debe procurar que las novillas mantengan bajas ganancias de peso vivo para reducir la edad al primer parto. Por otra parte, Núñez-Domínguez et al. (25) encontraron que vacas de tamaño mediano en la raza Tropicarne eran más precoces que las de tamaño pequeño o grande.

El IEP estimado en este estudio (18.4 meses) fue entre cuatro y cinco meses mayores que los publicados para ganado Romosinuano de Colombia y México (23,24,29), y que los publicados para otras razas bovinas tropicales en México (26,30,31). Para esta variable, es posible que las vacas Romosinuano de este estudio no pudieron satisfacer sus requerimientos nutricionales en el postparto temprano, provocando retrasos en el reinicio de la actividad ovárica postparto; y consecuentemente, generando aumentos en el IEP. Para IEP y TP, el comportamiento de las vacas Romosinuano fue afectado por el número de partos y la época de partos. Los mayores valores de IEP en los primeros partos podría deberse a que se presentan en edades donde las vacas no habían alcanzado su madurez en peso (sólo 67.3 % del peso maduro de la hembra al primer parto). Ossa et al (22,29) estimaron una tendencia a disminuir el IEP cuando se incrementó el número de partos en vacas Romosinuano de Colombia. Asimismo, Mejía-Baustista et al (26) reportaron menores IEP a medida que aumentó el número de partos en animales Bos taurus, Bos indicus y sus cruzas, en condiciones tropicales de México.

En conclusión, los modelos no lineales mixtos utilizados permiten describir adecuadamente el crecimiento de las hembras Romosinuano en el trópico con alimentación basada en pastoreo, aunque el modelo de Brody es más apropiado que los otros. En las condiciones ambientales de la unidad de producción estudiada, esta población de hembras presenta una tasa de crecimiento lenta, por lo que puede considerarse como raza de tamaño pequeño o mediano; con una aceptable eficiencia reproductiva, considerando las restricciones ambientales y genotipos disponibles en zonas tropicales. La actividad reproductiva es mejor cuando las hembras tienen un mayor porcentaje de peso maduro y se tienen partos en épocas de lluvias. Las hembras Romosinuano con menor tamaño maduro presentan edades al primer parto más temprano que vacas más pesadas.

En general, la caracterización del crecimiento y la aptitud reproductiva de hembras Romosinuano en las condiciones tropicales del hato en estudio, dan una idea de la situación que ocurre en los otros hatos similares que utilizan esta raza en el país, y es complementaria a las estimaciones realizadas con Romosinuano en otros países. Sin embargo, una mayor recolección de información permitiría una caracterización más precisa y representativa del potencial productivo de los hatos Romosinuano en México.

Los autores declaran no tener conflictos de interés.