La raza Limousin, originaria de Francia, como raza pura o en esquemas de cruzamiento^1,2 presenta cualidades productivas, reproductivas y de adaptación que han permitido su distribución en gran cantidad de países y de sistemas de producción^3,4,5; además, se ha utilizado en el desarrollo de razas sintéticas⁶. A México llegó por importaciones de Canadá y Estados Unidos que datan de la década de 1970; la Asociación Mexicana de Criadores de Ganado Limousin (AMCGL) se constituyó en 1989^7,8. Actualmente se distribuye en 17 estados como raza pura especialmente; aunque también es utilizada en esquemas de cruzamiento abierto y como base para la conformación de razas sintéticas, tales como Limousan (5/8 Limousin y 3/8 Angus)⁹ y Brahmousin (5/8 Limousin y 3/8 Brahman)¹⁰.

La AMCGL coordina el registro genealógico de pureza de raza y grados de pureza, así como los registros productivos que definen los criterios y objetivos de selección de la raza⁷. Los datos productivos asociados al crecimiento incluyen el peso vivo al nacer, a 120, 210 y 365 días de edad, con mediciones en el intervalo de más o menos 45 días a la edad especificada. Las mediciones de peso vivo generan una distribución de observaciones a través de la vida del animal, que en su conjunto se pueden utilizar para caracterizar y evaluar la curva de crecimiento. Los modelos no lineales (MNL) caracterizan y analizan la curva de crecimiento animal con base en la interpretación biológica y aplicaciones de los coeficientes de regresión, así como parámetros de crecimiento derivados de los coeficientes de regresión^11,12,13. Los coeficientes de regresión y parámetros de crecimiento son de trascendencia en la toma de decisiones para los programas de manejo, nutrición, reproducción y mejora genética^14,15,16,17. Con base en lo anterior, el objetivo del presente estudio fue la selección y ajuste de un MNL para describir y evaluar la curva de crecimiento en bovinos Limousin de México.

La base de datos estuvo conformada por mediciones de peso vivo, en el intervalo de peso al nacer a 500 días de edad en bovinos Limousin (PRZ; pureza de raza). Para definir la curva de crecimiento se evaluaron cuatro MNL: Brody (BRO), von Bertalanffy (BER), Gompertz (GOM), y logístico (LOG); todos conformados por tres coeficientes (β₁, β₂ y β₃) de regresión^12,13,18. En las ecuaciones de los MNL (Cuadro 1), y _i representa el peso vivo (kg) medido al tiempo t ; β₁, es el valor asintótico cuando t tiende a infinito, interpretado como el parámetro de peso adulto (PAD); β₂, es un parámetro de ajuste cuando y ≠ 0 y t ≠ 0; y β₃, es la tasa de crecimiento (TAC), expresando la ganancia de peso como proporción del peso total. Los modelos BER, GOM y LOG se caracterizan por describir el crecimiento con base en una curva sigmoide, para los cuales se calculó la edad (EPI; meses) y el peso (PPI; kg) al punto de inflexión. El modelo BRO presenta una curva de crecimiento con TAC constante sin punto de inflexión. Con los coeficientes de regresión se estimaron la edad para alcanzar el 50 % de madurez (E50M), la madurez obtenida a los 15 meses (GM15) de edad^19,20, así como la correlación (r_ac) entre TAC y PAD.

Los análisis se realizaron para cada sexo, con el método de Gauss-Newton del procedimiento NLIN del programa para análisis estadístico SAS²¹. La selección del modelo con mejor ajuste se realizó en función de^18,19: criterio de información Akaike [AIC= n*ln(sce/n) + 2k]; criterio de información Bayesiano [BIC= n*ln(sce/n) + k*ln(n)]; coeficiente de determinación [R²= (1 - (sce/sct))]; y, error estándar general o del modelo (EEG=scen-p-1. Donde: n = número total de datos; sce= suma de cuadrados del error; sct = suma de cuadrados total; k= número de parámetros en el modelo; ln = logaritmo natural. Para AIC y BIC, el modelo con el menor valor se consideró como de mejor ajuste.

La AMCGL administró un registro genealógico con diferentes grados de pureza (GPZ) con el objetivo de incrementar la población de bovinos Limousin a través de cruzamiento absorbente, con base en vacas cruzadas y sementales de PRZ. Con el modelo seleccionado como de mejor ajuste en la población de PRZ se caracterizó la curva de crecimiento en poblaciones definidas por cinco GPZ o generaciones: primera (PG) con ½ de Limousin; segunda (SG) con – de Limousin; tercera (TG) con 7/8 de Limousin; cuarta (CG) con 15/16 de Limousin; y, quinta (QG) con 31/32 de Limousin. En el Cuadro 2 se describe la base de datos analizada en PRZ y GPZ.

En la selección de modelos, dentro de sexo con AIC y a través de sexo con BIC, el modelo de mejor ajuste fue BER seguido de BRO y GOM; en todos los modelos el R² fue superior al 95 % (Cuadro 3). En el Cuadro 4 se presentan los resultados para los coeficientes de regresión y parámetros de crecimiento producto de los MNL evaluados. La estimación del PAD fue mayor para PRZ vs GPZ, en contraste, la TAC fue superior en los resultados de GPZ. El esquema de mejora genética del ganado Limousin en México incluye el peso al destete ajustado a 205 días⁷, con posible trascendencia en las curvas de crecimiento, dado que el punto de inflexión está ubicado en el periodo predestete. El modelo de BRO fue el segundo en la clasificación de modelos; sin embargo, presentó resultados atípicos para PAD, E50M y GM15. Todos los modelos presentaron r_ac debajo de -0.75 (Cuadro 4), lo cual señala que altos PAD no derivan de altas TAC. En la Figura 1 para machos y Figura 2 para hembras, se describe el crecimiento con base en el modelo BER para todos los genotipos evaluados.

Figura 1 Curvas de crecimiento para machos Limousin. Pureza, animales de raza pura; PG, 1/2 de Limousin; SG, 3/4 de Limousin; TG, 7/8 de Limousin; CG, 15/16 de Limousin; QG, 31/32 de Limousin 

Figura 2 Curvas de crecimiento para hembras Limousin. Pureza, animales de raza pura; PG, 1/2 de Limousin; SG, 3/4 de Limousin; TG, 7/8 de Limousin; CG, 15/16 de Limousin; QG, 31/32 de Limousin 

En bovinos Limousin de raza pura, Igarzabal et al³ en tres sistemas de producción reportaron como modelo de mejor ajuste GOM. En esquemas de cruzamiento de Limousin con Angus, Hereford y MARC III, Zimmermann et al¹⁷ utilizaron el modelo de BRO para caracterizar la curva de crecimiento y evaluar el peso vivo a la madurez. En bovinos Limousin x Friesian, representaron el crecimiento con base en el modelo GOM²². En la raza Madrasin, producto del cruzamiento de Limousin con Madura, el crecimiento presentó una curva tipo sigmoide, caracterizada con el modelo LOG²³. Curvas de crecimiento evaluadas con el modelo BER, fueron reportadas en bovinos Holstein²⁴, Pirenaica y Blonda³.

En México, diversos estudios han discutido contrastes en el tipo de curva de crecimiento a través de razas. Para curvas de crecimiento sin punto de inflexión, Domínguez-Viveros et al²⁵ en cinco razas de cebú en ganadería tropical, reportaron que los MNL de mejor ajuste fueron Brody, Meloum III y Mitscherlich; caso particular del modelo BRO, fue seleccionado con mejor ajuste en vacas Romosinuano²⁰, en bovinos Tropicarne¹⁹ y Salers²⁶. Para curvas de crecimiento tipo sigmoide, Contreras et al.²⁷ en vacas Jersey, Holstein y cruzas Jersey con Holstein los MNL seleccionados fueron GOM, LOG y BER, respectivamente; el modelo BER se ha reportado en bovinos Hereford²⁶.

La incorporación de vaquillas de reemplazo a la fase reproductiva es de trascendencia para el progreso genético y rentabilidad del hato. Este procedimiento transcurre en tres etapas²⁸: inicia maduración de la hipófisis detonada a cierta edad y peso; seguida con el desarrollo de los ovarios y crecimiento corporal; maduración del útero como consecuencia del desarrollo pituitario y su influencia hormonal sobre el crecimiento corporal y actividad ovárica, permite a la vaquilla el empadre y desarrollar la gestación. En diversos estudios se ha analizado la influencia de los parámetros de crecimiento sobre variables reproductivas^20,24,29; de esta forma, el punto de inflexión se ha asociado^13,30,31 con el inicio de la fase reproductiva. La edad al primer parto es un indicador del tiempo que tarda un animal en alcanzar su madurez sexual y reproducirse por primera vez, el empadre alrededor de los 15 meses, con edad al primer parto en torno a los 24 meses, tiene efectos positivos en la longevidad y productividad de la vaca^32,33. Con base en el modelo BER, se observan diferencias en las hembras de GPS vs PRZ para los componentes de la curva de crecimiento (Cuadro 4), las cuales se pueden atribuir a las diferencias genéticas a través de razas y los efectos de heterosis, producto del esquema de cruzamiento. En vaquillas de PRZ, la EPI se estimó a 3.7 meses con 131.2 kg para PPI; en GPZ, la EPI y el PPI estuvieron en los intervalos de 2.9 a 3.7 meses y 110.9 a 124.4 kg, con valores promedio de 3.3 y 116.4, respectivamente. Para E50M en hembras, en PRZ se estimó a 10.6 meses y para GPZ estuvo en el intervalo de 8.9 a 10.5 meses, con un valor promedio de 9.5. Para GM15 en hembras, en GPZ el valor promedio fue de 90.5 % y la estimación en PRZ fue de 87.9 %. En contraste, para hembras de otras poblaciones y con base en el modelo BER: Contreras et al²⁷ para Holstein, Jersey y cruzas estimaron EPI (meses) y PPI (kg) en los intervalos de 7.4 a 9.8 y 115.0 a 151.7, respectivamente; Domínguez-Viveros et al²⁵ en cinco razas de cebú, reportaron estimaciones de EPI y PPI en los intervalos de 3.9 a 11.7 y 107.2 a 230.9, respectivamente; en Romosinuano²⁰, Tropicarne¹⁹ y Siboney²⁹ los resultados de EPI - PPI fueron de 15.5 - 132.5, 7.7 - 180.5 y 5.9 - 152.4, respectivamente.

Con respecto a los machos, la selección de sementales se realiza en PRZ y se incorporan a la reproducción a partir del año de edad; no obstante, reducir la edad de incorporación a la reproducción reduce el intervalo generacional y trasciende en el progreso genético³⁴. La curva de crecimiento puede incidir en el desarrollo de la fase reproductiva; en razas Bos taurus, los eventos fisiológicos asociados a la reproducción inician a los seis u ocho meses; la madurez y capacidad reproductiva se precisan por la calidad del semen, con variaciones por efectos del peso vivo, tasa de crecimiento, circunferencias escrotal, entre otros factores³⁵. Los resultados señalan que los machos en PRZ alcanzan el 50 % de madurez a partir de los 13 meses, con valores superiores al 80 % a partir de los 15 meses (Cuadro 4). En contraste²⁶, en bovinos Hereford y Salers se reportaron grados de madurez de 68.2 % y 76.6 %, al año, respectivamente. Por otro lado, los indicadores de la curva de crecimiento están asociados a la rentabilidad en la producción; la TAC tiene efecto en la edad y peso al sacrificio; grado de madurez trascendencia en la eficiencia y composición de la canal^16,22. Diferencias y derivaciones de la curva de crecimiento con relación a la producción, para machos de raza pura y diversos cruzamientos, han sido valoradas por diversos autores^15,36,37.

El modelo de mejor ajuste fue von Bertalanffy, el cual describió una curva de crecimiento tipo sigmoide, con diferencias en los parámetros de crecimiento a través de los genotipos evaluados. Las estimaciones del punto de inflexión están en el contexto del crecimiento predestete.