INTRODUCCIÓN

En Colombia existe una diversidad de ovinos, denominados ovinos criollos, resultado del cruce de diferentes razas, desde la época de la colonización y que se adaptaron a diferentes regiones del país (^{Vergara Garay et al. 2017a}). Dentro de estos ovinos, se encuentran los Sudán y los Etíope, considerados como dos poblaciones con características y potencial morfoestructural diferente (^{Flórez et al. 2018}; ^{Vivas A. et al. 2020}).

En Colombia, el sector ovino ha tomado importancia dentro de los diferentes sistemas de producción y mercadeo. Este sector, se caracteriza por la producción de cárnicos, la crianza de animales como reproductores y vientres, para abastecer las granjas que quieran multiplicar sus rebaños o iniciar nuevos proyectos (^{Acero-Plazas, 2014}). Además, la producción ovina del país, se caracteriza por la utilización de razas introducidas y criollas y de cruzamientos entre ellas, con el fin de fortalecer el vigor híbrido, para la producción de carne. Estos sistemas productivos son extensivos y semi-intensivos, con baja a mediana incorporación de tecnología y, por lo general, pertenecientes a pequeños productores (^{Vergara-Garay et al. 2017b}). Estas explotaciones son el eje central de la agricultura familiar y el abastecimiento de pequeños nichos de mercado, cumpliendo con una función económica importante (^{Acero-Plazas, 2014}).

La conservación de los recursos zoogenéticos tiene importancia, a nivel mundial, ya que, para la implementación de programas de mejoramiento, en los cuales, se seleccionen aquellos animales superiores para características productivas y que tenga adaptabilidad a un entorno determinado, facilita este proceso y se puede lograr mayor progreso genético, en un tiempo determinado (^{Noriega M. et al. 2018}). Se debe tener en cuenta que, bajo condiciones de altas temperatura, se activan mecanismos de termorregulación, en los cuales, animales que no estén adaptados a estas altas temperaturas, como las de los bosque seco tropical y bosque muy seco tropical, aumentan su frecuencia cardiaca y consumo de agua y disminuyen el consumo de alimento (^{Pérez et al. 2020}). En este sentido, los ovinos de pelo son un recurso genético, que debe ser valorado.

A pesar de los esfuerzos, acompañamientos y campañas realizadas para incentivar la crianza de ovinos en el país, el inventario nacional es bajo. Según el Instituto Colombiano Agropecuario (^{ICA, 2020}), la explotación de la especie en el país cuenta con solo 1´682,767 de animales. La mayor parte de la población, se registra en los Departamentos de La Guajira, con 42,1 %; Magdalena, con 11,8 %; Cesar, con 8,8 %; Boyacá, con 7,8 % y Córdoba, con 5,4 %, que agrupan el 75,9 % del total. De acuerdo con estas cifras, las poblaciones de ovino de pelo criollo colombiano OPC habitan, principalmente, en zonas de bosque muy seco tropical, que se presenta una época del año en la que la disponibilidad de forraje es escaza (^{Acero-Plazas, 2014}).

Aunque se debe considerar que en los Departamentos de Córdoba y Cesar la ovinocultura se encuentra en un periodo de crecimiento, en donde se trata de producir estos animales como una fuente primaria de ingresos y no como una especie secundaria a otros sistemas de producción y donde los animales relevantes son los bovinos. Cabe resaltar que el sacrifico de ovinos aumentó, a nivel nacional, entre el segundo trimestre del 2016 y el segundo trimestre del 2020, en un 46,2 %, mientras que las exportaciones se han reducido, ya que entre 2017 y del 2018, las exportaciones fluctuaron entre 56 y 60 toneladas de carne año, respectivamente, pero para el segundo trimestre del 2020, no se reportan exportaciones (^{DANE, 2022}).

Por lo anterior, y teniendo en cuenta las características de adaptabilidad del OPC a condiciones del Trópico bajo colombiano y la información que aportan los diferentes modelos no lineales para describir el crecimiento de las poblaciones animales, se deben llevar a cabo investigaciones, que permitan describir el crecimiento de los ovinos y demuestren las ventajas que brindarían estos animales a los sistemas productivos, por lo que el objetivo de esta investigación fue determinar la curva de crecimiento en ovinos de pelo colombiano, en zonas de bosque muy seco tropical, del Departamento de Córdoba y bosque muy seco tropical, del Departamento del Cesar, Colombia.

MATERIALES Y MÉTODOS^{Gompertz, 1825 Nelder, 1961 Von Bertalanffy, 1938 Brody, 1945}

Localización. Se utilizaron las bases de datos provenientes de cinco granjas, ubicadas en la zona de bosque seco tropical, del Departamento de Córdoba y seis bases de datos de granjas pertenecientes al bosque seco muy tropical, del Departamento del Cesar. Para la utilización de las bases de datos, se contó con los permisos respectivos de los propietarios de las granjas. Dichas bases de datos tenían información sobre la identificación del padre, de la madre, fecha de nacimiento y de destete de los animales, peso al nacer y los pesajes periódicos cada 30 días.

Toma de datos y manejo de animales. Se estudiaron 99 animales puros de los biotipos Sudán (59 animales) y Etíope (40 animales), de los cuales, 42 (19 machos y 23 hembras), se encontraban en la zona de bosque seco tropical y 57 (46 machos y 11 hembras), en el bosque muy seco tropical. Los nacimientos ocurrieron entre diciembre de 2015 y julio de 2016. A partir del nacimiento, cada 30 días, los animales fueron pesados, hasta una edad aproximada de 12 meses, tiempo, en el cual, alcanzan el peso al sacrificio. Los pesajes, se realizaron en una báscula digital AWS DINA 50k, con capacidad para 50 kg, considerando para análisis los animales con un mínimo de 10 pesajes y máximo 13. El destete se realizó a los 90 días.

Los animales, se manejaron bajo condiciones de pastoreo en praderas de pasto natural angleton (Dichantium aristatum) y guínea (Megathyrsus maximus), Bothriochloa pertusa, Brachiaria humidicola y Brachiaria brizanta, de acuerdo con las pasturas de cada sistema de producción; además, contaban con disponibilidad de agua y sal mineralizada a voluntad. Asimismo, los animales fueron suplementados en la época de escases de lluvia con ensilado de maíz Zea mays, debido a la poca disponibilidad de forraje en los sistemas de producción. Esta suplementación se hizo, con el fin de que los animales no tuvieran pérdidas de peso en este periodo.

Curvas de Crecimiento. A partir de la información generada por los pesajes, se construyó una base de datos, que se utilizó para evaluar las curvas de crecimiento de los siguientes modelos no lineales:

Dónde: yi representa el peso del animal en el tiempo t; β ₀ representa el peso asintótico del animal cuando “t” tiende a infinito; β ₁ es un parámetro de ajuste cuando y ≠ 0 o t ≠ 0; β ₂ es un índice de madurez o medida de precocidad y su valor, a medida que aumenta, sugiere una tasa de crecimiento mayor.

El porcentaje de madurez a una edad de cuatro (M4) y seis meses (M6) y edad al alcanzar el 75 y 95 % de madurez (EM75 y EM95, respectivamente), se estimaron a partir del modelo que mejor ajuste tuvo para describir el crecimiento de los animales en estudio, de acuerdo con las siguientes expresiones:

Análisis estadístico. Para la estimación de los parámetros de las curvas de crecimiento evaluadas, se utilizó el procedimiento NLIN del SAS (^{SAS, 2007}). Como criterios para seleccionar el modelo que mejor describiera el crecimiento de los ovinos estudiados, se utilizó el coeficiente de determinación (R²), el criterio de información Bayesiano (BIC), el criterio de información de Akaike (AIC) y el porcentaje de curvas convergentes (PCC), dado por el número de curvas que convergieron respecto al total de curvas. A partir del modelo que mejor ajuste tuvo, se determinó, para los parámetros del modelo (β ₀ , β ₁ , y β ₂ ), los porcentajes de madurez y las edades a la madurez; sí hubo diferencia estadística entre biotipo racial, sexo, zona y época de nacimiento (época de lluvias y época de sequía), a través de un análisis de varianza, mediante el procedimiento GML del SAS (^{SAS, 2007}).

RESULTADOS Y DISCUSIÓN

En la tabla 1, se presentan los estimadores de los parámetros β ₀ , β ₁ , β ₂ , y los criterios para la selección del modelo que mejor se ajustó a los datos. De acuerdo con los resultados, se puede observar que, en los modelos de Gompertz, Bertalaffy y Logístico, las diferencias entre los criterios de selección son muy pequeñas, lo que supondría que cualquiera de ellas ajustaría, satisfactoriamente, para describir el crecimiento de la población ovina estudiada; sin embargo, teniendo en cuenta los valores AIC y BIC estimados para el modelo de Brody, se infiere que es el modelo que mejor se ajusta, además el valor de R² no dista mucho de los otros modelos, similar a lo reportado por ^{Nimase et al. (2018)}. Asimismo, ^{Ullah et al. (2013)}, ^{Sieklicki et al. (2016)}, ^{Balan et al. (2017)}, ^{Gautam et al. (2018)} y ^{Mohammadiam et al. (2019)} encontraron que el modelo Brody fue el que mejor ajuste tuvo al describir el crecimiento de ovinos machos comerciales Texel, Mecheri, Sonadi y Kordi, respectivamente.

A diferencia de este estudio, se reportó que el modelo de Brody no convergió bien, al describir el crecimiento en ovinos de diferentes grupos genéticos y sexo (^{López et al. 2018}). En este mismo sentido, en corderos Norduz fueron evaluados los modelos de Gompertz, Von Bertalanffy, Brody, Logístico y Exponencial Negativo y se demostró que el modelo Logístico fue el que mejor describió el crecimiento en los ovinos (^{Daskiran et al. 2010}). ^{Da Silva et al. (2012)} estudiaron los modelos de Brody, Richards, Von Bertalanffy, Gompertz y Logístico para describir el crecimiento en ovinos Santa Inés y sugirieron el modelo Logístico, como el de mejor ajuste. A la misma conclusión llegaron ^{Lupi et al. (2015)}, al estudiar en ovinos de la raza Segureña utilizando los modelos de Brody Verhulst, Gompertz, Von Bertalanffy y Logístico. Por otro lado, ^{Paz et al. (2018)} utilizaron los modelos de Brody, Richards, Von Bertalanffy, Gompertz y Logístico para describir el crecimiento en ovinos de la raza Morada Nova, siendo el modelo de Gompertz el que presentó el mejor ajuste. En ovinos Rahmani fueron evaluados los modelos Brody, Logístico, Gompertz y Von Bertalanffy, siendo, para este estudio, el modelo de Gompertz, el que presentó mejor ajuste para describir su crecimiento (^{Rashad et al. 2017}). Por su parte, otros investigadores, al describir las curvas de crecimiento en ovejas Mehraban iraníes, utilizando los modelos de Brody, Exponencial Negativo, Logístico, Gompertz y Von Bertalanffy, encontraron que el modelo que mejor describió la curva de crecimiento en la población de estudio fue el Von Bertalanffy (^{Hojjati & Hossein-Zadeh, 2018}).

De acuerdo con lo expresado, se puede evidenciar que es específico de cada población tener un modelo particular para describir en forma más apropiada su crecimiento, lo cual, está directamente relacionado con el tipo racial, sistema de alimentación, condiciones ambientales, en las que se desarrollan los individuos. Al ser el modelo de Brody el de mejor ajuste en la población de estudio, puede servir para hacer ajustes en el manejo nutricional de los animales y establecer metas, de acuerdo con la edad y el peso, adecuados para el sacrifico.

Al comparar los valores estimados en esta investigación con otros estudios, se encontraron valores superiores para β₀ y β₁ en ovinos Santa Inés por Criollo (β ₀= 46,8; β ₁=0,95) (^{Simanca et al. 2016}) y en ovinos de Pelo Criollos (β₀= 58,5; β₁=0,98) (^{Vergara-Garay et al. 2016}), utilizando el modelo Brody. Los valores estimados para β ₂, por ^{Sarmento et al. (2006)} y ^{Vergara-Garay et al. (2016)}, fueron superiores a los encontrados en el presente estudio (0,0094 y 0,0054, respectivamente). En ovinos criollos, en Colombia, se encontraron valores inferiores para los diferentes parámetros (β₀=33,8; β₁= 0,91 y β₂=0,0046) (^{Simanca S. et al. 2017}). Por otro lado, en Ovejas Madgyal fueron encontrados valores inferiores para β₀=25,2, β₁=0,85 y superior β₂=0,21 (^{Nimase et al. 2018}). En hembras F1 Santa Inés por Criollo y Criollo, se estimó β₀ (43,5) superior al presente estudio, β₁ similar (0,94) y β₂ (0,0046) inferior (^{Vergara-Garay et al. 2017b}). En este estudio, se debe considerar que en los sistemas de producción en la zona de estudio se presentan limitaciones en la toma de información, por ello, solo se utilizaron datos hasta alrededor del año de edad, por lo que podría haber una sobrestimación del parámetro β ₀.

Con relación al porcentaje de madurez alcanzado a los seis meses, el cual, es importante en la especie ovina, ya que en las granjas y mercados regionales es el referente, como la edad adecuada, para sacrificar a los corderos y para utilizar las hembras como reproductoras de reemplazo. De acuerdo con los resultados obtenidos (Tabla 2), se sugiere la necesidad de mejorar el manejo alimentico de los animales y hacer selección de animales superiores, con el fin de mejorar los porcentajes de madurez y las edades a la madurez, para obtener animales más pesados al sacrificio en menor tiempo (^{Vergara Garay et al. 2019}). En la figura 1, se muestra la curva de crecimiento promedio de los ovinos de pelo colombiano, de acuerdo con el modelo Brody, realizada a partir de los parámetros estimados para este modelo, descritos en la tabla 1. Para la realización del gráfico, se hizo una proyección hasta los 720 d de edad.

Por otro lado, como el mejor modelo para describir el crecimiento en ovinos de pelo colombiano fue el Brody, al considerar los efectos biotipo racial, sexo, zona y época para este modelo (Tabla 3), se encontró que para el efecto sexo no se detectaron diferencias (P >0,05), entre las medias de los parámetros estimados. Para el efecto biotipo racial fue significativo (P ≤0,05) para todos los parámetros estimados, con excepción de β ₁ , mientras que el efecto zona fue significativo para β ₀ , β ₂ y M4. Con relación al efecto época, solo fue significativo (P ≤0,05) para el parámetro β ₁ . Se puede observar en la figura 2, la representación de las curvas de crecimiento por biotipo racial, sexo, zona y época de los ovinos de pelo colombiano, de acuerdo con el modelo Brody.

Similar a lo estimado en este trabajo, se ha reportado diferencia significativa del sexo sobre los estimados de β ₀ , demostrando superioridad de los pesos de los machos respecto a las hembras. Esto se explica, por el efecto anabólico producido por la hormona testosterona en los machos, que actúa como promotor de crecimiento, reflejando un mayor peso a mayor madurez (^{Daskiran et al. 2010}; ^{Lupi et al. 2015}).

De acuerdo con los resultados obtenidos, los animales de la raza Sudán son más precoces, debido a que tienen mayor valor de β ₂ , ya que este representa la tasa de madurez y la velocidad de crecimiento del animal al llegar a su tamaño adulto. Animales con mayor valor de β ₂ representan madurez precoz, en comparación con individuos de valores inferiores de β ₂ y de peso inicial semejante, por lo que, el valor de β ₂ representa con buena precisión, la variación en la velocidad relativa de crecimiento del animal (^{Simanca S. et al. 2017}). Esto se confirma, por los valores obtenidos en los porcentajes de madurez y edades de los biotipos estudiados (Tabla 3), en la cual, el biotipo Sudán fue superior en un 34,76 %, para el porcentaje de madurez, a los seis meses y alcanzando una edad al 95 % de madurez, 824,3 días antes, respecto al biotipo Etíope. Estos resultados coinciden biológicamente con la alta variabilidad genética encontrada en OPC, por ^{Noriega M. et al. (2018)}, lo cual, hace interesante el biotipo Sudan para futuros planes de selección asistida por genes, en pro de aumentar la productividad en los sistemas ovinos, manejándolo como base genética para su explotación.

Se debe considerar que, al evaluar el crecimiento en ovejas, se pueden encontrar divergencias en la convergencia de los modelos, por las particularidades de cada población, en cuanto a tipo racial, alimentación, sexo, entre otros factores (^{Lupi et al. 2015}). Además, al evaluar diferentes sistemas de producción, se requiere determinar todos los parámetros de forma colectiva para poder decidir cuál es modelo de curva de crecimiento que mejor se ajusta a determinado sistema de producción (^{Teixeira et al. 2012}). Por ello, para cada sistema de producción, debe existir un modelo que mejor ajusta para describir el crecimiento de los animales.

Como conclusión, el modelo de Brody se puede utilizar en la descripción del crecimiento de ovinos, de acuerdo con los parámetros de la curva de crecimiento, para instaurar proyecciones en el comportamiento productivo. Asimismo, los estimados, los porcentajes y las edades, a la madurez en las poblaciones en estudio, muestran que la raza, con mejoras en factores no genéticos, como la alimentación, permiten que sean candidatos para implementar programas de selección, con el fin de incrementar los pesos de los animales al llegar a la edad adulta y al sacrificio, mejorando sus índices productivos. También, se puede afirmar que el biotipo Sudán tuvo una mayor tasa de crecimiento, aunque un peso adulto inferior al biotipo Etíope, por lo que puede ser un biotipo adecuado para su explotación en las zonas en estudio.