INTRODUCCIÓN
Actualmente el programa de mejoramiento genético de búfalos en Cuba se basa en las evaluaciones genéticas de las madres para producción de leche y en los resultados de las pruebas de comportamiento de los machos en pastoreo que provienen de las mejores madres para producción de leche. Las pruebas de comportamiento que se realizan son evaluaciones del crecimiento de machos durante aproximadamente un año en condiciones de pastoreo de bucerros de una misma empresa a partir del destete a 8 meses, pero de sus diferentes rebaños bubalinos. Al finalizar la prueba se dictamina entre estos animales los mejores para futuros sementales atendiendo a su crecimiento y, además, la producción de leche de la lactancia correspondiente de sus madres en términos porcentuales, ambos ponderados respectivamente (70:30) según Mitat. et al, 2010.
No se ha realizado un estudio estadístico de la evolución de los indicadores de las pruebas que se iniciaron en 2011 a fin de conocer su evolución y tomar medidas oportunas con datos ajustados por los efectos ambientales más sobresalientes. Este artículo tiene el objetivo de realizar un estudio de la evolución del peso al destete, peso vivo a 20 meses, ganancia en peso de 8-20 meses y peso por edad a 20 meses en las pruebas de comportamiento de búfalos en la Empresa Genética Pecuaria Los Naranjos en el período del 2011 - 2021 a partir de las CMC que representan medias ajustadas para los efectos de año de entrada, valores genéticos y error. En área del trópico no se señalan programas de selección para crecimiento a pesar de destacar su importancia económica (Martínez-Reina et al., 2020) y estimaciones de parámetros genéticos, por lo que no existen estudios de tendencias del crecimiento. No obstante, sí se ha informado sobre indicadores reproductivos en Egipto (Amin et al., 2021).
MATERIALES Y MÉTODOS
Se tomó la información de la Empresa Genética Pecuaria “Los Naranjos” Artemisa en Cuba donde se realizaron pruebas de comportamiento (PC) de búfalos en condiciones de pastoreo sobre pastos mejorados (Cynodon nlemfuensis) a razón de 0.6 a 1,0 animales /ha, durante un periodo de 12,72 meses como promedio, a la que sólo se enviaron los animales que en 3 ocasiones durante la campaña de nacimiento de cada año, tuvieron buenos pesos al destete y procedían de madres que evaluaban de buenas productoras de leche en los años de nacimiento comprendidos entre 2011 y 2019 con 435 individuos totales. Se conformó una base de datos a partir de las informaciones recogidas en las pruebas de comportamiento que incluyeron el peso al destete a 8 meses (PD), peso final ajustado a 20 meses de la prueba de comportamiento (PF), ganancia en peso vivo desde el inicio hasta el final de la PC (Gan), peso por edad al final de la PC (PPE), edad a inicio de la prueba (EI), edad al final de la PC (EF) y duración de la PC (Dur).
Se utilizó un modelo mixto en el SAS (2013) con el Proc GLIMMIX (Wolfinger y O’Connell 1993; Schabenberger 2006) para procesar los indicadores de PD, PF, Gan y PPE después de identificar con el Proc Severity la distribución más ajustada a los datos. Se estudió como efecto fijo al AE a la PC que tuvo 9 niveles correspondientes a los años 2011/2012 (unidos debido a insuficiente información del primer año) hasta 2020 y como efecto aleatorio al animal anidado al AE a la PC. Se efectuaron además como efectos fijos, ajustes por covarianza del PD a la edad de entrada a PC, el peso final y PPE a la edad final de la PC mientras que la Gan se ajustó a la Dur. Se plotearon las CMC con los EE para conocer sus cambios a través de los AE. Se empleó el modelo mixto siguiente:
Yijklm= μ + αi + aj(αi) + Βk (Xl -μl) + eijklm
μ = Media general o constante común a todas las observaciones
αi = Efecto fijo del i-ésimo año de entrada a la PC (i=2011/2012…, 2020).
Βk (Xl -μl)= Efecto del ajuste lineal por covarianza k de la edad inicial a la PC para PD, EF para PF y PPE o Dur para Gan. (l=1,2,3,4 y k=1,2,3,4).
aj(αi) = Efecto aleatorio del j-ésimo animal o identificación individual IDT (j=1,2…, 435) anidado en el i-ésimo AE a la PC.
eijklm= Error aleatorio asociado a las observaciones N (o, δ2 e)
Con los resultados obtenidos se realizaron gráficas de líneas, considerando en el eje de las X los años (2011/2012-2020 o sea 9 años) y en las Y las medias de comportamiento PD o PF en kg y GAN o PPE en g/día que provenían de las 435 observaciones totales.
RESULTADOS Y DISCUSIÓN
El AE afectó todos los indicadores evaluados, mientras que la covarianza no alcanzó la significación solamente con el PF; lo que pudo deberse a que individualmente se ajustó para esta edad (Tabla 1). Esto significa que cuando no se ajusta el PF por edad final se precisa de la covarianza por EF, pero que en el caso de la Gan y PPE que depende de dos mediciones, si es preciso su consideración en la modelación. El ajuste a la normalidad no fue necesario para el PD contrariamente a lo esperado pues a PC solo ingresaron los individuos con buenos PD y provenientes de madres con las mejores producciones de leche en la lactancia correspondiente. Las otras medidas sí precisaron de las transformaciones Gamma con link log
De acuerdo con la gráfica 1 se aprecia una cierta variabilidad del peso al destete a través de los AE, lo cual se expresa porque a pesar del ajuste por edad de entrada estos difirieron entre sí con los valores más bajos próximos a los 130 kg y más altos en 150 kg. O sea, 20 kg de diferencia, aunque estable si se atiende a que a través de los años no se mantuvo una caída ni pérdida sostenida. No obstante, debe velarse por los factores que generaron los bajos valores observados de individuos que entran a PC en algunos años(Figura 1).
Indicadores | Efectos | Ajuste a la normalidad/función de enlace | |||
---|---|---|---|---|---|
Año de entrada | Covarianza | IDT(Año) | Error | ||
Peso destete | *** | EI(***) | - | - | Gauss/Identidad |
Peso final | *** | EF(NS) | - | - | Gauss/Log |
Ganancia | *** | Dur(***) | - | - | Gauss/Log |
Peso por edad | *** | EF(***) | - | - | Gauss/Log |
Grados de libertad Numerador/denominador | 8/398 | 1/27 | 407 | 398 | - |
***: (P<0,0001) NS: no significativo (P<0,3277)
Lo mismo aconteció para los restantes indicadores evaluadas, donde el efecto de AE a la PC fue siempre importante, pero en estos casos se observa una caída a partir de los AE 2016-2017. Esto significa que a partir de entonces se han producido diferentes efectos de carácter ambiental como genéticos que están conspirando contra la mejor expresión del potencial genético. Entre estos pudiera pensarse en:
De carácter genético:
Incremento de la consanguinidad en la población bajo estudio.
Deficiencia en la selección para crecimiento de los individuos de la población particularmente de los progenitores.
De carácter ambiental:
Inadecuada alimentación expresada en la calidad y cantidad de los pastos y forrajes recibidos y otros suplementos alimenticios.
Inadecuado manejo expresado por la carga o encierre nocturno con insuficiente alimento.
Estudios realizados con empleo del programa ENDOG en esta misma población (Fraga, 2022) sugieren que no existen elementos que puedan atribuirse a que la consanguinidad sea una posible causa del deterioro ya que la consanguinidad esperada a través de la contribución desbalanceada, la consanguinidad media computada, el coeficiente de parentesco promedio y el número de antecesores explicando 50% de los genes tienen valores que evidenciaron la no existencia de peligros de consanguinidad.
Lo anterior significa que existe necesidad de optar por métodos más precisos de mejora genética y no por efectos derivados de la consanguinidad.
Los resultados del peso al destete pueden considerarse aceptables pero no así los del peso vivo a 20 meses (Figura 2) que son bajos si se comparan con los 400 kg o más referenciados por Borghese et al. (2010) en Italia con búfalos a 15 meses confinados o de Brasil en pastos con 300- 351 kg por Malhado et al. (2008) y Jorge et al. (2005) y los de 336 kg a 18 meses por Crudeli et al. (2007) con búfalos Mediterráneos aunque similares a los de Colombia con 277 kg informados por Agudelo et al. (2009) en Colombia (citado por Bolivar et al., 2012). Similar efecto ocurre con la ganancia y el peso por edad (Figuras 3 y4).
En lo referente a los efectos ambientales que están gravitando en esta empresa u otras de Cuba, se debe garantizar la mejor expresión del potencial genético de los animales en estos momentos donde se aprecia en general una reducción de indicadores y en la que los búfalos, necesitan de ingentes esfuerzos también para su mejor desempeño. Datos colaterales con esta base de datos evidencian que en los últimos años de entrada (2018, 2019 y 2020) las CMC de la PLM de estos individuos tuvieron también descensos ya que mostraron valores de 905, 877 y 851 kg superiores al señalado de 742 kg por Mitat (2022) en similar período en la misma empresa en toda la masa en producción de la empresa.
CONCLUSIONES
El peso vivo final ajustado individualmente a 20 meses no precisa de ajustes por covarianza a edad final lo cual sí es necesario en el caso del PD, PPE o la Gan con respecto a la edad de entrada, final o duración de la PC. Estos elementos deben ser tomados en cuenta en la selección que se realice para garantizar la uniformidad de las mediciones.
Se evidenció descenso del peso vivo final a partir de 2016, particularmente en último AE donde alcanzó valores similares a los del período 2011-2012. Se mantiene, no obstante, estabilidad a través de los AE del peso al destete lo cual indica que deben atenderse aspectos de manejo que están afectando el período de PC aunque no se descarta posibilidad de aplicar una selección más eficiente y precisa