Modelado de curvas de lactancia estándar en vacas holstein con diferente paridad y escenarios climáticos

Introducción

La producción de leche es un proceso complejo destinado a proporcionar nutrición y protección inmunológica a la cría. Durante la gestación hay una fase de activación celular, reemplazando las células adiposas de la glándula mamaria por tejido glandular especializado y una fase de regresión celular después del parto terminando con el secado de la vaca (Strucken et al., 2015). Estas dos fases fisiológicas durante una lactancia completa generan un patrón en forma de curva. La función gamma incompleta ha sido el modelo más utilizado para describir la curva de lactancia en mamíferos tales como bovinos, búfalos y pequeños rumiantes (Macciotta et al., 2004). El rendimiento máximo y persistencia de una curva de lactancia son características que definen el rTot (Jeretina et al., 2013). Se ha documentado que la producción de leche por vaca es dependiente de la trayectoria de la curva de lactancia y está influenciada por factores ambientales y genéticos (Torshizi, 2016). En el noroeste de México, anualmente se presenta el efecto del estrés por calor debido a las condiciones de humedad y temperatura. En esta zona, la programación de partos en la población de ganado lechero está altamente modelada por la estación climática del año, debido a que la zona de confort térmico es superada durante siete meses aproximadamente. La curva de lactancia representa una importante herramienta para investigaciones destinadas a aumentar la capacidad de producir leche por vaca y adaptación a distintos escenarios de producción por vaca. La descripción de la forma de la curva de lactancia, según la paridad y estrés por calor, puede proporcionar una estimación del patrón de rTot. Por lo tanto, el objetivo del presente estudio fue analizar los parámetros y componentes de la curva de lactancia y su relación con el rendimiento total de leche en vacas lecheras raza Holstein con diferente paridad y manejadas bajo diferentes escenarios climáticos.

Materiales y métodos

La toma de datos inicio en el año 2015 y finalizó en el 2018. Participaron cinco establos lecheros localizados en el Valle del Yaqui, Sonora. Se colectaron inicialmente 472,712 registros diarios de leche (rLeche), después de una depuración exhaustiva fueron seleccionados 278,317 rLeche de 1,352 lactancias completas y 1,086 vacas Holstein. Las lactancias se agruparon según paridad (Par1, 2, 3, 4 y 5 incluyendo vacas con > 5 lactancias dentro de Par5). Además, utilizando el Índice de Temperatura y Humedad (ITH; Hahn, 1999) según el rLeche diario, se formaron tres escenarios climáticos-productivos (ClasITH): 1= Confort (Conf; ITH< 68), 2= Umbral de estrés (Umb; ITH 68-71) y 3= Estrés medio-moderado (Mod; ITH 72- 79). Para calcular los parámetros y componentes de la curva de lactancia, se utilizó la función gamma incompleta siguiendo la metodología descrita por Macciotta et al. (2004). El rTot fue calculado a partir de los rLeche. Se examinaron los efectos de Par y ClasITH mediante PROC GLM de SAS. Para la comparación de cuadrados medios se utilizó la prueba de Duncan.

Resultados y discusión

En la tabla 1, se observa la comparación de medias de acuerdo al número de paridad para los parámetros y componentes de la curva de lactancia. Solamente los parámetros a y c fueron estadísticamente diferentes (P< 0.001), mientras que todos los componentes de la curva de lactancia fueron diferentes (P< 0.001) entre cada Par. El mayor tPmax y Pst se encontró en Par1. El mayor rTot se observó en las vacas de Par3 y 5. En la tabla 2, según la ClasITH los parámetros y componentes de la curva fueron estadísticamente diferentes (P< 0.0001), excepto rMax. Las vacas con estrés Mod fueron las que presentaron valores inferiores en los componentes de la curva de lactancia en comparación con los grupos Conf y Umb. El rTot para vacas en estrés moderado fue 9% inferior. El comportamiento de los parámetros y componentes de la curva fue similar a lo encontrado en otras investigaciones (Tekerli et al., 2000; Jeretina et al., 2013; Strucken et al., 2015; Torshizi, 2016; Duque et al., 2018). Una explicación fisiológica para el rTot de Par1 es que las vacas primerizas están pasando por un proceso de desarrollo de la glándula mamaria y sus células adiposas son remplazadas por tejido glandular contrarrestando la disminución del rendimiento total de leche a medida que avanza su número de lactancia (Stanton et al., 1991). Lo anterior implica que el tPmax y Pst se prolonge. Además, su capacidad de contener leche es reducida por lo que rMax y rTot fueron menores en comparación con las vacas de más de un parto (Stanton et al., 1991; Strucken et al., 2015). A partir de un ITH de 68 la vaca experimenta estrés por calor (Zimbelman et al., 2009), y en este punto la producción de leche tiende a reducirse en la mayoría de las vacas. Sin embargo, Berman (2011) menciona que el ganado puede sufrir un proceso de adaptación a las altas temperaturas, lo que al parecer ocurrio con un gran número de vacas dentro de el grupo Umb y Mod.

La paridad y el estrés por calor influyeron en los parámetros y componentes de la curva de lactancia, así como en el rendimiento total de leche. Se recomienda el modelaje de los componentes de la curva de la lactancia para incrementar la precisión en los esquemas de selección animal enfocada al rendimiento total de leche.

Tabla 1. Comparación de medias entre número de paridad sobre los parámetros y componentes de la curva de lactancia estándar en vacas Holstein en producción.

Tabla 2. Comparación de medias según clasificación de ITH sobre los parámetros y componentes de la curva de lactancia estándar en vacas Holstein.

Modelado de curvas de lactancia estándar en vacas holstein con diferente paridad y escenarios climáticos

Biodisponibilidad relativa de Mepron® respecto a Smartamine® M medida según el método de selenometionina en leche

ADM