ISSN L 2644-3856

Revista Investigaciones Agropecuarias

Volumen 2, N°2. pp. 49-68

Junio 2020 - Noviembre 2020

Panamá

Recepción: 16 de abril de 2020   Aceptación: 4 de mayo de 2020

PERFIL DE LA CURVA DE LACTACIÓN Y PRODUCCIÓN LÁCTEA COMERCIAL EN VACAS CEBUÍNAS Y CRUZADAS EN FINCAS LECHERAS CON BAJA TECNOLOGÍA EN EL TRÓPICO

Edil E. Araúz¹, Eduardo E. Araúz² y Ángela Montezuma³

^1* Universidad de Panamá, Facultad de Ciencias Agropecuarias, Departamento de Zootecnia. Panamá

² Ingeniero Agronomo Zootecnista, Servicios profesionales y asesorías en sistemas e industria de la Producción Animal, Bioseguridad, Procesamiento y Conservacion de los Productos de Origen Animal. Panamá

³ Zootecnista y Extensionista en Ganadería, Ministerio de Desarrollo Agropecuario, Región 1 en la Provincia de Chiriquí, Panamá.

*edil.arauz@up.ac.pa, eduarauz_20@hotmail.com, angelamontezuma25@gmail.com

RESUMEN

La curva de lactación y la producción de leche comercial fueron evaluadas en 250 vacas multíparas cebuínas y cebuínas x europeo tipo leche en catorce fincas lecheras Grado C; bajo pastoreo, estrés calórico moderado y un ordeño manual diario con ternero en la época lluviosa en el trópico húmedo de bajura; donde el estrés calórico diurno máximo ocurrió a la 01 PM (HR 76.71%, TAbs 32.7°C, RSD 676.3 kcal/m2hr, ITH 86.59 e ITHRSDVVmax 92.35). La producción comercial acumulada y la longitud lactacional en el grupo cebuíno resultó en 692.07 kg en 177 días; ½ Bos taurus tipo leche x ½Cebú 1267.35 kg en 213 días y en las ¾ Bos taurus tipo leche x ¼ Cebú 2012.16 kg en 262 días (P<.01). El perfil típico de la curva de lactación mostró una producción diaria inicial, máxima y final en las cebuínas de 3.42, 4.62 a los 14 días y 2.23 kg/vaca (P<.05); 5.63, 8.63 a los 26 días y 3.27 kg/vaca (P<.01) en las ½ Bos taurus tipo leche x ½Cebú y 7.50, 10.95 a los 34 días y 4.40 kg/vaca (P<.01 en las 3/4 Bos taurus x ¼ Cebú). La curva de lactación en las cebuínas fue Y(kg/vaca día)=13.09572121X^{- 0.0733187}e^{-0.005064242X} (R²=0.38); en las ½Cebú x ½ Bos taurus Y(kg/vaca-día)=8.5935X^0.13048e^{-0.005115067X} (R²=0.68) y en las ¾ Bos taurus lechero x ¼Cebú Y(kg/vaca-día)=12.7632X^0.24051e^{- 0.007117078X} (R²=0.72). La producción lechera fue potenciada con el mayor aporte genético tipo leche en los cruces; sin embargo, el entorno ambiental y el manejo afectaron la curva de lactación y redujeron la habilidad lechera fenotípica en todos los grupos genéticos en el trópico húmedo.

PALABRAS CLAVES: cebuínas, cruzadas, estrés calórico, lactación, lecherías, trópico, potencial lechero

PROFILE OF THE LACTATION CURVE AND COMERTIAL MILK PRODUCTION IN ZEBU AND CROSSBRED COWS IN LOW TECHNOLOGYCAL SMALL DAIRY FARMS IN THE TROPIC

ABSTRACT

The lactation curve and commertial milk production were studied in 250 lactating multiparous zebu and crossbred cows (zebu x Bos taurus milk type) in 14 small dairy farms Grade C under grayzing, mild heat stress, calf stimulation and one milking by day across raining time in lowland of humid tropic. The maximum heat stress occurred at 01 PM (RH 76.71%, ATdb 32.7°C, DSR 676.3 kcal/m2hr, THI 86.59 and THIDSRWS 92.35). Total commertial milk production and lactation length resulted for zebu 692.07 kg in 177 days, ½ Bos taurus milk type x ½ Zebu 1267.35 in 213 days and for ¾ Bos taurus type milk x ¼ Zebu 2012.16 kg in 262 days (P<.01). The tipic profile of the lactation curve in the zebu group showed initial, peak and final daily milk yield of 3.42, 4.62 at 14 days and 2.23 kg/cow (P<.05); in ½ Bos taurus milk type x ½ Zebu 5.63, 8.63 at 26 days (P<.01) and in 3/4Bos taurus milk type x ¼ Zebu 7.50, 10.95 at 34 days and 4.40 kg/cow (P<.01). The lactation curve resulted in Zebu cows Y(kg/cow-day)=13.09572121X^-0.0733187e^{-0.005064242X} (R²=0.38); in ½ Zebu x ½ Bos taurus Y(kg/cow-day)=8.5935X^0.13048e^{-0.005115067X} (R²=0.68) and in ¾ Bos taurus lechero x ¼ Zebu Y(kg/cow-day)=12.7632X^0.24051e^{-0.007117078X} (R²=0.72). Milk production was potentiated by the proportion of Bos taurus milk type in the crossbred; however, the environment and management altered the lactation curve and reduced phenotype milk ability in all genetic groups in the humid tropical climate.

KEYWORDS: crossbred, dairy farms, heat stress, lactation, tropic, milk potential, Zebu

INTRODUCCIÓN

La producción lechera tropical presenta varios modelos tecnológicos de fincas; entre ellas, que poseen las menores condiciones en infraestructuras, facilidades físicas, inversiones, tipo de alimentación basada en el forraje verde a través del pastoreo, sistema de registros simplificados o ausentes, animales con bajo potencial genético lechero, ausencia de equipos para la preservación y el cuidado de la leche (Araúz y Guevara, 2015). Estas fincas han sido clasificadas de acuerdo con la legislación sanitaria y de calidad láctea en Panamá en Grado C; siendo el mayor grupo de productores que alcanzaron unas 5197 (89.73% de las fincas lecheras en Panamá entre el 2010 y 2011) según el MIDA (2012); las cuales representan en la actualidad unas 6453 fincas (MIDA, 2018).

Las fincas lecheras Grado C producen el 47% de la producción nacional en Panamá; al mismo tiempo que contribuyen como fuente de trabajo y sustento económico para las familias rurales y constituyen un apoyo para el funcionamiento y la operación de las pequeñas y medianas industrias de lácteos nacionales (MIDA, 2012, 2018). Sin embargo, con la evolución de la situación lechera en Panamá, las fincas con la menor tecnología han mostrado fortalezas y ventajas operativas ya que tienen menores inversiones fijas, animales más resistentes al trópico, menor dependencia en la mano de obra, menores costos de producción y un menor compromiso financiero (Araúz, 2019). A pesar de la prevalencia de las fincas lecheras Grado C en Panamá y otros países latinoamericanos; aún no se cuenta con estudios sectoriales nacionales que describan el potencial lechero y la curva de lactación; destacando aspectos como: producción inicial, máxima y al cierre, tiempo al pico lactacional, persistencia y tiempo en producción. El objetivo de esta investigación fue establecer cuáles son los descriptores críticos de la curva de lactación en vacas multíparas cebuínas y cruzadas cebú x europeo tipo leche (Holstein, Pardo Suizo); así como su rendimiento lechero potencial con fines comerciales con baja tecnología lechera, alimentación a base de forraje verde, un ordeño manual diario con ternero y periodo de acceso al amamantamiento libre de unas cinco horas.

Características e indicadores de la Curva de lactación y el potencial lechero en la vaca

El ciclo de vida biológico de la hembra bovina incluye los sectores del crecimiento, desarrollo, pubertad, establecimiento reproductivo efectivo, la gestación, el parto, la lactación, la nueva gestación, el periodo del secado y el nuevo parto y que luego de iniciar la fase reproductiva se describe como el intervalo entre partos (Araúz, 2019). Este ciclo es importante para las acciones y directrices del manejo apropiado en nutrición, reproducción, salud, genética, ambiente, bienestar y bioseguridad (Crowe et al., 2018).

La trayectoria de la curva de lactación es esencial para establecer el manejo optimo en el manejo nutricional, reproductivo y sanitario de la vaca lechera moderna para hacer los modelos de producción más eficientes fisiológicamente y alcanzar una mayor rentabilidad económica (Cabrera, 2018).

Uno de los modelos más utilizados para describir la curva de lactación en la vaca ha sido la función gamma incompleta descrita por la expresión: Y_t = a t ^b e ^{–c t} donde Yt es la producción de leche diaria en kilogramos, a es la producción al inicio de la lactación, t es el tiempo lactación en días, b es la pendiente hacia la máxima producción y c es la pendiente del máximo al cierre lactacional y e es el valor neperiano 2.718282 como modelo integral (Wood, 1967). El modelo linearizado ln (yt)= ln (a) + b ln (t) – ct permite estimar el tiempo para la máxima producción láctea (b/c), la máxima producción (a(b/c) b_e−b) y la persistencia láctea [−(b + 1) ln (c)] según Atashi et al., (2012). La producción de leche es afectada por varios factores, los cuales modifican la curva lactacional; sobresaliendo el plano nutricional (Walker et al., 2004) y el estrés calórico que afecta la producción (Gupta et al., 2016; Das et al., 2016) y la reproducción (West, 2003; Hansen, 2009; Crowe y Williams, 2018) debido a las alteraciones en la fisiología general, el metabolismo, la carga calórica corporal y en la cinética de la termorregulación (Araúz, 2017).

Los modelos regresivos lactacionales predicen la trayectoria de la producción de leche y contribuyen con la adecuación del manejo y el aprovechamiento de la capacidad productiva en las fincas lecheras (Cabrera, 2018). Por ende, los registros biolactacionales son esenciales para evaluar la capacidad funcional de las vacas y definir el manejo integral; en especial en el medio trópico, donde hay limitantes nutricionales y el estrés térmico que reducen simultáneamente la productividad lechera (Araúz, 2010).

Requerimientos energéticos y alimentación de la vaca lechera en el trópico

El patrón de los requerimientos nutricionales energéticos para el ganado lechero ha sido establecido por la NRC (2001) y la ARC (1999); a partir de las cuales se han derivado las normativas de alimentación y los ajustes del manejo para potenciar la capacidad lechera de los bovinos según su valor genético. El perfil energético para la producción de leche aun para la ganadería con baja tecnología lechera es un factor de valor productivo (Araúz, 2015).

Existe un balance favorable para la producción lechera tropical a través de las lecherías grado C ya que la alimentación depende del forraje verde, los animales presentan un menor tamaño corporal (McDowell, 1981), la habilidad para utilizar la energía para producir leche es baja (NRC, 2001) y el perfil genético favorece el almacenamiento de reservas corporales y limita la movilización de los nutrientes dietéticos y endógenos para la producción de leche (Stothard et al., 2004). La habilidad del tejido mamario, el patrón endocrino y el perfil metabólico se encuentra bajo el control de los genes que regulan el metabolismo del tejido adiposo y el uso de otras reservas durante la lactación en las hembras cebuínas; limitando el uso de nutrientes para la biosíntesis láctea (Nayeri y Stothard, 2016).

El gasto energético para la locomoción en el ganado lechero es 0.00045 Mcal/kg por cada kilómetro recorrido y para el pastoreo 0.0012 Mcal por cada kg de peso corporal (NRC, 2001). En el trópico se gasta energía extra para realizar la termólisis obligatoria que demanda la homeotermia para evitar el distanciamiento calórico extremo corporal y la alteración normotérmica (McNab, 2001). Las condiciones del trópico incluyen el estrés calórico efectivo durante la fase diurna (Araúz, 2010); lo que altera el balance calórico y el bienestar biotérmico (Polsky et al., 2017). En consecuencia, el pastoreo, la locomoción y la termólisis de la vaca en producción demandan energía que reduce la habilidad lechera (Araúz, 2006).

El periodo de lactación representa la fase metabólica en la que se utilizan las reservas orgánicas (grasa corporal, músculo y huesos); así como las reservas hepáticas (glucógeno, vitaminas liposolubles) y minerales óseos que son utilizadas para la biogénesis de la leche; dependiendo del potencial lechero, del estado lactacional y de la composición de la leche (Knegsel et al., 2007). Las vacas cebuínas y cruzadas Holstein x cebuínas bajo diferentes perfiles de alimentación y nutrición energética restringida mostraron que utilizan la energía a medida que aumenta la producción lechera (Araujo et al., 2018).

Los forrajes tropicales son reconocidos por su bajo valor energético para favorecer la producción lechera al contener entre 0.12 y 0.35 Mcal/kg de EN leche; incluyendo los pastos de corte y las gramíneas pasto señal (Brachiaria decumbens) y el pasto estrella (Cynnodon dactylon) según varios estudios (Rosas et al., 1985). La presión energética sobre las reservas corporales durante la lactación en el trópico está acompañada de un rápido deterioro de la condición corporal de la vaca en producción aun en las fincas lecheras grado C (Araúz, 2018).

PARTE EXPERIMENTAL

Ubicación del estudio y área geográfica

El estudio fue desarrollado en el Distrito de Alanje (El Tejar, Santo Thomas, Canta Gallo, La Pita, Mostrenco, Querébalos, La Martina y Guarumal) en la Provincia de Chiriquí; utilizando catorce fincas lecheras Grado C en la época lluviosa del 2008 al 2012.

Animales experimentales para la medición lactacional y los registros de producción

Se utilizó 250 vacas multíparas entre dos y cuatro partos de un total de 355 animales candidatos que fueron distribuidas en cebuínas 58, ½ Cebú x ½ Europeo tipo leche 106 y ¾ Europeo tipo leche x ¼ Cebú 86. El tipo lechero incluyó los aportes de la raza Holstein y Pardo Suizo como parte de los cruzamientos. Se estableció el registro de la producción de leche semanal con la identificación individual según la fecha postpartal de producción.

Peso corporal, peso metabólico, superficie corporal, alimentación y manejo en el ordeño

El peso corporal fue determinado por perimetría torácica utilizando la cinta OVNY de tela y plástico resistente confeccionada en Colombia (Inalmet, 2007) calibrada para cebuínos y bovinos cruzados de doble propósito.

El peso metabólico determinó según la NRC (2001), la superficie corporal para cebuínos según Curtis (1981) para los cebuínos SC (M²) = 0.12 (PCkg) ^0.60; mientras que para los ½ Cebuínos x ½ Europeo Lechero (SC (M²) = 0.135(PC_kg) ^0.58) y en los ¾ Bos taurus lechero x ¼ Cebuíno SC (M²) = 0.143(PC _kg)^0.57 según Araúz (2018).

El manejo animal incluyó acceso al pastoreo libre en el día y en la noche; con acceso a las crías entre las 8 AM y 01 PM); separando los terneros a la 01 PM hasta el siguiente día en la mañana. El manejo incluyó el amamantamiento restringido, higiene de los pezones con agua clorinada, el ordeño manual y el pesaje de la producción cada siete días. La alimentación dependió del forraje verde con pastoreo libre y se utilizó el pasto suazi (Digitaria swazilandensis), señal (Brachiaria decumbens) y en algunas fincas (3) se utilizó la asociación suazi y pangola (Digitaria decumbens Stent). Se utilizó sal mineralizada (Sal Cruda 65%) y una fuente de Minerales en un 35%) a base de calcio, fósforo, magnesio, manganeso, zinc, hierro, cobre, cobalto y zinc complementada con vitaminas (Laboratorios Bayer, 2015).

Entorno microambiental, sensibilidad y presión calórica ambiental

El medio microambiental de las fincas fue evaluado en el campo abierto; tomando una medición psicrométrica en la mañana (06 am) y en la tarde a la 01:00 PM cada 15 días en la época lluviosa. Se utilizó el termohigrómetro digital para registrar la Temperatura Ambiental de bulbo seco (oC) y la Humedad Relativa (%) en cada finca codificada para conformar la base de datos microclimáticos. La evaluación psicrométrica (Martínez, 2012), el ITH según Armstrong, (1994) y Dikmen y Hansen (2006) y el ITHRSDVV según Mader et al. (2006). La radiación solar directa y la velocidad del viento se obtuvieron en Unión Fenosa para Alanje. Los índices de sensibilidad calórica normal y ajustada se determinaron según: ITH = (TAbs x 0.80) + [(HR/100 x (TAbs – 14.4)] e ITHRSDVV = 4.51 + ITH – (1.992 x VV m seg-1) + (0.0068 x W m^-2 hr^-1).

Medición de la Producción de Leche y el contenido de grasa láctea

Se peso la producción de leche obtenida por ordeño manual diario, registrada según la finca, la vaca y la fecha a través del formulario de campo; generando 8168 observaciones para la producción de leche (1508 en cebuínos, 3392 en ½ Bos taurus tipo leche x ½ Cebú y 3268 en ¾ Bos taurus tipo leche x ¼ Cebú). Se registró el peso del recipiente vacío, el recipiente más la leche obtenida y el peso de la leche (kg); utilizando una romana digital marca Constant (HQ), seriada 14192-115E con capacidad de 50 kg. La fecha de parto de cada animal y la producción de leche fue anotada a partir del 8^vo al 10^mo día postparto hasta culminar su ordeño comercial. Una muestra de leche fue tomada cada 14 días de cada vaca para el análisis de la grasa láctea, analizada en el laboratorio de Fisiología Animal y Ciencia Lechera del Departamento de Zootecnia.

Diseño Experimental y Análisis Estadístico

El modelo utilizado incluyo: años (A_i), grupos genéticos (B_k), estado lactacional semanal (C_m) y las covariables pesos corporal en la primera semana de producción y el cambio de temperatura ambiental diurna; conformando el arreglo factorial más covarianza; descrito por como:

Y_ijk = u + A_i + B_k + C_m + (BC)_ik + b₁ (X_ij – û₁) + b2 (X_ij – û₂) + e_ikj donde

Y_ijk = es la producción de leche diaria del animal jmo en el año imo, en el grupo genético kmo en función del estado lactacional mmo intervenido por el peso corporal (kg) en la primera semana de producción y el cambio térmico diurno (oC) entre las 7 AM y la 1 PM.

u = Media general

A_i: Efecto del imo año del estudio (imo = 1 (2008), 2(2009), 3(2010), 4(2011), 5(2012).

B_k: grupo genético (kmo = cebuínas (1), 1/2Cebu x ½ taurina (2) y 3/4BTTL x 1/4Cebu (3) Cm: Estado lactacional semanal (mmo = 1, …..,37º semanas de producción láctea diaria) (CB)km: Interacción del Grupo genético x el Estado lactacional semanal

b₁(Xij – û1): Peso corporal de la vaca en la primera semana de producción comercial en kg. b2 (Xij – û2): Cambio térmico diurno en oC cada 15 días.

e_ijkm = Residuo o Error Experimental; incluyendo el efecto de fincas en el residuo o error.

Las medias cuadradas fueron ajustadas por covarianza y se utilizó el error tipo III para ajustar el análisis de varianza. La curva de lactación por grupos genético fue generada según la Función Logística Gamma Incompleta (Wood, 967): Y_t (kg/día) = at^be^{–c t}; siendo Y la producción de leche diaria (kg) en el momento lactacional t (día), a es la producción de leche inicial (comercial), b es la pendiente desde el inicio al máximo de producción, c es la pendiente del pico al cierre lactacional y e es el número neperiano; integrados según la programación en SAS (2004) descrita por Herrera y Barreras (2001).

RESULTADOS Y DISCUSIÓN

Microambiente, presión calórica y estrés térmico efectivo

El entorno microambiente diurno en la época lluviosa para el periodo en estudio mostró una temperatura ambiental promedio antes de la salida del sol de 23.5°C, Humedad Relativa 88.3%, ITH 73.33 e ITH ajustado 73.74; indicando la menor presión calórica y la menor influencia térmica para los bovinos (Tabla I). Para la 01 Pm se alcanzó un comprometimiento térmico microambiental (TAbs 32.7 °C, HR 76.7%, ITH 86.59 e ITH ajustado 92.35) que represento un estrés calórico severo; aunque es evidente que la tensión calorica trasciende a partir de las 10 en adelante; siendo perjudicial para la salud térmica y productiva de la vaca lechera según Dikmen y Hansen (2009). La modificación térmico ambiental (AM vs PM) fue de +25.23%; incluyendo la contribución de la temperatura ambiental, la humedad relativa y la radiación solar directa (Tabla 1). La mayor presión y estrés calórico se alcanzó después del medio día, cuyo nivel fue severo con las implicaciones negativas en la fisiología sistémica, la termorregulación, la conducta y la producción lechera (West, 2003; Araúz et al., 2015; Polsky et al., 2017). Por ende, se requieren ajustes en el manejo y medio físico en las horas más calurosas para cuidar el bienestar térmico endógeno, reducir los efectos sistémicos negativos (Araúz, 2006; Arias et al., 2008), evitar o reducir las alteraciones metabólicas (Slimen et al., 2016; Gupta et al., 2016), reproductivas (Jordan et al., 2003; Hansen, 2009; Das et al., 2016) y las pérdidas económicas implícitas (St-pierre et al., 2003).

Tabla 1. Indicadores del perfil microclimático diurno promedio en el Distrito de Alanje enla época lluviosa durante el estudio lactacional.

☼ Araúz, E. E. (2010, 2018) ╪NOAA (2014) *Knigth, J. (2006, 2020) Ϯ (Mader et al., 2006)

Peso corporal y metabólico y estado hematimétrico de las vacas experimentales

La categorización del peso corporal al parto en las vacas multíparas fue diferente entre los tres grupos genéticos (P<.01); evidenciando que a medida que aumentó la contribución genética del Bos taurus tipo leche en los cruces también aumentó el peso corporal, el peso metabólico y la superficie corporal (tabla II). El perfil hematimétrico no fue diferente para el hematocrito entre los grupos genéticos (P>.05); sin embargo, el hematocrito, hemoglobina y el conteo de glóbulos rojos mostraron una ligera disminución al aumentar el aporte genético del ganado taurino lechero en los cruzamientos; siendo la hemoglobina y el conteo de glóbulos rojos menores ante los cebuínos (P<.05). McDowell (1985) indica que las vacas cebuínas pesan 350 kg a diferencia del presente estudio con 402 kg; siendo superior en un 14.85%. Las vacas media sangre promediaron 437 kgy el grupo con más contribución lechera pesó 478 kg (Tabla 2); siendo inferior en 16.8%. Estas diferencias obedecen a la influencia del ambiente; formando parte de la variación fenotípica que desarrollan los bovinos cebuínos y cruzados en las condiciones tropicales; donde el forraje verde es la principal fuente de alimentación y el estrés calórico se combinan; moldeando la biológico de la vaca cebuíno y cruzada bajo el perfil subnutricional y la adaptación somática para minimizar la influencia negativa del estrés térmico.

Tabla 2. Indicadores somáticos y hematimétricos al parto según el grupo genético.

Medias a b y b c difieren al 5% (P<.05); a c difieren al 1% (P<.01); letras iguales no difieren al 5% (P>.05)

Producción de leche comercializada y prevalencia biolactacional en las fincas

La producción de leche comercializada fue heterogénea a través de fincas evaluadas debido a las contribuciones de la disponibilidad del forraje verde, tipo de pastura, planificación forrajera y manejo de las pasturas, composición genética de los animales y composición genética de las vacas en producción al momento de la evolución. El potencial lechero de cada finca lechera Grado C fue considerado en un ordenamiento secuencial en base al promedio diario potencia por vaca; pudiendo observarse que un total de cuatro fincas (28.35%) mostraron una producción diaria de 3.12 a 3.97 kg/vaca; constituyendo las fincas de menor potencial lechero. Las fincas intermedias oscilaron entre 4.04 y 4.94 kg/vaca e incluyeron un total de seis fincas que representaron el 42.86%. El tercer grupo de fincas lecheras grado C presentó un rango diario promedio entre 5.4 y 9.81 kg/vaca y representaron el 28.57% de las unidades ganaderas evaluadas. Las fincas estuvieron en el mismo entorno microclimático, pero la variación de la producción lechera diaria comercial fue 39.18% con un promedio de 5.23 ± 2.05 kg/vaca (ver Figura 1).

Figura 1. Ubicación de las fincas lecheras Grado C evaluadas en el distrito de Alanje.

La ponderación e importancia de la producción lechera diaria en las fincas Grado C representa un aporte esencial en los sistemas ganaderos tropicales con tecnología limitada ya que los ingresos dependen de la comercialización de la leche y de la venta de los terneros luego del destete (Guerra, 2015). Sin embargo, estos resultados indican que los ingresos y la productividad ganadera se encuentra afectada por la baja habilidad lechera de orden genético; además de confrontar las limitaciones nutricionales de los pastos tropicales (Rosas et al.,1988, 1994), los efectos directos del estrés calórico en la fase diurna (Araúz, 2006, 2017), la influencia negativa de una baja frecuencia del ordeño (Gómez, 2002) y la práctica de permitir el a amamantamiento por cuatro a cinco horas en las vacas en producción con el propósito de alimentar parcialmente las crías (Araúz, 2019). Estos factores se conjugan alrededor de la ganadería lechera de doble propósito y conducen a la generación de la baja productividad ganadera en Panamá (Corella, 2016).

Indicadores de la curva de lactación y la producción de leche total comercializada

El empleo de la función gamma incompleta sobre la producción de leche y el estado lactacional permitió determinar los descriptores de la curva de lactación; los cuales fueron diferentes entre los grupos genéticos (P<.001); así como la duración lactacional y la producción láctea comercial (P<.001). La producción de leche inicial, máxima y al cierre en las cebuínas fue 3.42, 5.58 y 2.33 kg/día; en las vacas ½ Cebú x ½ Bos taurus tipo leche fueron 5.63, 8.63 y 3.27 kg/día; y en el grupo ¾ Bos taurus tipo leche x ¼ Cebú resultó en 7.5, 10.95 y 4.4 kg/día; evidenciando la mayor diferencia al aumentar el potencial lechero (ver Tabla 3).

Las inflexiones diferenciales al inicio de la lactación de los cruces de media y tres cuartos ganado europeo tipo leche x cebuínos superaron la producción inicial cebuína en 64.5 y 128%; sin embargo, al considerar los máximos de producción absoluta, se observó un incremento de 54.66 en el cruce ½ Bos taurus tipo leche x ½ Cebú y de 96.24% en las vacas producto del cruce ¾ Bos taurus tipo leche x ¼ Cebú. El cierre lactacional en las cebuínas (2.33 kg/vaca día) fue superado por las vacas ½ Bos taurus tipo leche x ½ Cebú en 40% y en las ¾ Bos taurus tipo leche x ¼ Cebú en un 89% sobre las cebuínas como referentes (Figura 2).

Tabla 3. Descriptores lactacionales según el grupo genético en lecherías Grado C.

Medias con las letras a b y b c entre los grupos genéticos difieren al 1% (P<.01).

Figura 2. Valor crítico referente en la producción lechera inicial, máxima y al final de la lactación en fincas lecheras Grado C según el grupo genético y fenotípico.

La curva de lactación fue diferente entre los grupos genéticos; destacando que el periodo de producción fue más corto ocurrió en las cebuínas (177 ± 12 días); mientras que en las vacas ½ Bos taurus tipo leche x ½ Cebú fue 213 ± 16.75 días; siendo ambos superados por las vacas ¾ Europeo x ¼ Cebú con 262 ± 21.14 días (Tabla 4). El término comercial se refiere al periodo que las vacas se mantuvieron en producción para la venta de leche; aunque las mismas continuaron con sus crías como parte del perfil de manejo en las lecherías grado C.

La producción de leche y la longitud lactacional fueron diferentes entre los grupos genéticos (P<.001; P<.04) y según el estado lactacional (P<.001). La producción láctea total de índole comercial por vaca resultó en las cebuínas 692.07 kg en 177 días, ½Cebú x ½ Bos taurus lechero 1267.35 kg en 213 días y en las ¾ Bos taurus tipo leche x ¼ cebú 2012.16 kg en 262 días con un ordeño manual diario y acceso al amamantamiento por la cría (P<.001).

La producción de leche comercial aumentó entre las vacas cebuínas y las ½ Cebú x ½ Bos taurus tipo leche en +575.28 kg (+1.8312X) con un diferencial lactacional de +40 días (+1.2259X); mientras que el grupo genético ¾ Bos taurus tipo leche x ¼ Cebú aumentó con respecto a las cebuínas en +1320.00 kg (+2.908X) en 262 días de producción. La grasa láctea disminuyó desde el grupo cebuíno hasta el grupo con la mayor influencia genética taurina tipo leche. La producción láctea ajustada al 4% de grasa resultó en 798.99 kg en las cebuínas; en las ½ Bos taurus lechero x ½ Cebú 1339.59 kg y en las ¾ Bos taurus tipo leche x ¼ Cebú fue 1933.69 kg. La producción de leche acumulada para las vacas ½ Bos taurus tipo leche x ½ Cebú según McDowell (1985) fue 2830 kg contra 1267.35 kg en 213 días (44.78%); mientras que el grupo con más genética lechera alcanzó 2012.16 kg (39.53%) de 5090 kg.

Tabla 4. Componentesdecurva de lactación según la funcióngammaincompletaparalas vacas cebuínas, ½ cebú x ½ Bos taurus y ½ Bos taurus lechero x ½ Cebuínas.

Figura 3. Trayectoria de la curva de lactación según la función logística gamma incompleta en vacas cebuínas, ½ Bos taurus tipo leche x ½ Cebú y ¾ Bos taurus x ¼ Cebú.

La forma de la curva de lactación en el grupo de las vacas multíparas ¾ Bos taurus x ¼ Bos indicus siguió la trayectoria con una alta aproximación al patrón lactacional en las razas lecheras Holstein, Pardo Suizo y Jersey (Visser y Wilson, 2006); aunque otros modelos sectoriales (Grossman y Koops, 1988) y polinomiales (Araúz, 2008) pueden adaptarse. La curva de lactación en las vacas cruzadas ½ Bos taurus lechero x ½ Cebú se redujo en su expresión inicial, máxima, final y en la producción total según el grupo genético.

En las vacas cebuínas multíparas se observó que la curva de lactación fue reducida en duración (177 días), producción inicial (3.42 ± 0.49 kg), máximo de producción de 5.58 ± 0.61 kg y en el cierre lactacional 2.33 ± 0.53 kg. Estos indicadores estuvieron asociados con el tiempo de producción con fines comerciales y el ordeño manual diario en la mañana en las fincas lecherías Grado C evaluadas en el trópico de bajura durante la época lluviosa.

CONCLUSIONES

El modelo de producción lechera tropical con baja tecnología presentó marcadas diferenciasentre los grupos genéticos, cuyo potencial lechero fue limitado por la subnutrición energética según la alimentación a base de forraje verde, la presencia el estrés calórico diurno moderado, el modelo de amamantamiento y el ordeño con ternero al pie; los cuales de manera conjunta disminuyeron la habilidad lechera fenotípica, comercial y productiva de las vacas y fincas.

Los puntos críticos y la trayectoria de la lactación en las vacas cebuínas y cruzadas fue heterogénea; pero a su vez, evidencio la importancia del aporte proporcional de la genética tipo leche con las razas Holstein y Pardo Suizo en los cruces lecheros comerciales para el medio tropical, aun bajo condiciones de tecnología limitada y en presencia de otros factores ambientales negativos como la alimentación y la presión calórica ambiental propias del medio tropical húmedo. El aumento del aporte genético del Bos taurus tipo leche en los cruzamientos incrementó los indicadores lactacionales y el rendimiento lechero comercial; incrementando la producción de leche inicial, máxima, final y la duración del periodo postpartal con valor comercial y económico en el marco operativo de las fincas lecheras a pesar de los factores limitantes.

El bajo potencial energético del forraje verde, el estrés calórico tropical, el manejo del ternero y el ordeño redujeron el potencial biológico lactacional y comercial en las vacas cebuínas y cruzadas; evidenciando que el mayor aprovechamiento productivo en las fincas Grado C demanda ajustes en alimentación, nutrición, microambiente y manejo para favorecer la habilidad lechera e incrementar la productividad en las fincas lecheras con bajo nivel tecnológico ubicadas en las zonas bajas del trópico húmedo, aun en la época lluviosa.

AGRADECIMIENTO

En primera instancia, se les agradece a los productores de leche Grado C del Distrito de Alanje en los corregimientos de Santo Tomas, La Pita, El Tejar, Guarumal, Mostrenco, Querébalos, la Martina y Canta Gallo por haber colaborado en esta investigacion al permitir utilizar sus facilidades físicas, los animales y el personal que labora en las fincas lecheras Grado C para el levantamiento de la información de campo; a ustedes nuestro eterno agradecimiento por el apoyo brindado.

DECLARACION DE INTERESES

Los autores declaran que no existe conflicto entre los autores, con los productores, dentro de la Universidad de Panamá u otras instancias por haber elaborado este articulo con valor científico y técnico a partir de la investigación desarrollada originalmente por el autor principal. Por tal motivo declaramos reiteradamente que no hay diferencias sobre el contenido, conformación y orientación del documento científico presentado aquí para la Revista Científica de la Facultad de Ciencias Agropecuarias de la Universidad de Panamá.

REFERENCIAS BIBLIOGRÁFICAS

Araujo Carvalho Pedro Henrique de, Ana Luiza da Costa Cruz Borges, Ricardo Reis e Silva, Helena Ferreira Lage, Paolo Antonio Dutra Vivenza, Jose´ Reinaldo Mendes Ruas, Elias Jorge Facury Filho, Rodrigo Liberio Arau´jo Palhano, Lucio Carlos Gonc¸alves, Iran Borges, Eloısa de Oliveira Simões Saliba, Diogo Gonzaga Jayme, Antonio Ultimo de Carvalho. Energy metabolism and partition of lactating Zebu and crossbred Zebu cows in different planes of nutrition. 2018. Citado: 15/03/2020, Disponible: PLOS ONE | https://doi.org/10.1371/journal.pone.0202088. doi: 10.1371/journal.pone.0202088; disponible: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6097685/

Araúz, E. E. (2019). Sectores críticos en el ciclo de vida de la hembra bovina tipo leche; cuidados biotécnicos y manejo para el buen desarrollo, salud, bienestar y productividad. Citado: 24/03/2020. Disponible: www.engormix.com

Araúz, E. E. (2018). Valor térmico de las superficies para los animales en las pérdidas calóricas pasivas. En: Fisiología de la Adaptación y Producción Animal. Carrera de Zootecnia, Departamento de Zootecnia, Facultad de Ciencias Agropecuarias, Universidad de Panamá.

Araúz, E. E. (2017). Influencia del color del pelaje sobre el comportamiento térmico corporal, cinética de la sobrecarga calórica y alteración cardiorespiratoria circadiana en vacas lecheras cruzadas (6/8 Bos taurus x 2/8 Bos indicus) bajo estrés calórico en el trópico húmedo. Citado: 05/05/2018. Disponible en http://www.veterinaria.org/revistas/redvet 2017 Volumen 18 Nº 7, Revista electrónica de Veterinaria - ISSN 1695-7504.

Araúz, E. E. (2010). Principales registros biológicos para evaluar la capacidad funcional de la vaca lechera y su importancia para mejorar el manejo y la eficiencia en la producción lechera. Citado: 153/2020. Disponible: www.engormix.com

Araúz. E. E. (2006). El estrés calórico y sus efectos negativos sobre la fisiología, metabolismo, reproducción y eficiencia de la producción en el ganado bovino de leche. Departamento de Zootecnia, Facultad de Ciencias Agropecuarias, Universidad de Panamá.

Araúz, E. E., Araúz, E. E. y Caballero, S. (2015). Situación actual de la producción bovina de leche en Panamá. En: Generalidades y situación actual de la ganadería en Panamá. La Ganadería en América Latina y el Caribe. Editores: Rafael Núñez, Rodolfo Ramírez, Salvador Fernández, Omar Araujo, Miguel García y Tito Efraín Diaz. Pag 422 – 426. ALPA, FAO, IICA, BBA, México.

Arias, R. A., Mader, T. L. y Escoba, P. C. (2008). Factores climáticos que afectan el desempeño productivo del ganado bovino de carne y leche. Arch Med Vet 40, 7-22.

Atashi, H., Abdolmohammadi, A. R., Asaadi , A. Akhlaghi , A., Dadpasand, M. y Jafari Ahangari, Y. (2011). Using an incomplete gamma function to quantify the effect of dystocia on the lactation performance of Holstein dairy cows in Iran. Citado: 03/04/2020. https://doi.org/10.3168/jds.2011- 4954

Cabrera, V. E. (2018). Invited review: Helping dairy farmers to improve economic performance utilizing data-driving decision support tools. Department of Dairy Science, University of Wisconsin- Madison, WI 53706, USA. Animal (2018), 12:1, pp 134–144 © The Animal Consortium 2017 doi: 10.1017/S1751731117001665.

Corella, J. (2016). Situación económica de la ganadería de leche y carne en Panamá. Conferencia en las Jornadas Agropecuarias del Banco Nacional, Bugaba, Chiriquí, República de Panamá.

Crowe M. A. y Williams E. J. (2018). Triennial lactation symposium: effects of stress on postpartum reproduction in dairy cows. Journal of Animal Science, 90,1722–7.

Curtis, E. (1981). Psychometric approach to determine heat stress in animal production. Environmental management in agriculture. USA.

Das R., L. Sailo, N. Verma, P. Bharti, J. Saikia, I. y Kumar, R. (2016). Impact of heat stress on health and performance of dairy animals: A review. Veterinary World, EISSN: 2231-0916 Citado en internet: 02 de Abril, 2020. Disponible: www.veterinaryworld.org/Vol.9/March-2016/7.

Dikmen, S. y Hansen, P. J (2009). Is the temperature – humidity in the best indicator of heat stress in lactating dairy cows in a subtropical environment? Journal of Dairy Science, 92,109 – 116.

Grossman, M. y Koops W. J. (1988). Multiphase analysis of lactation curves in dairy cattle. Journal of Dairy Science, 71, 1598.

Guerra, P. (2015). Situación actual de la ganadería de carne y las perspectivas futuras para el mejoramiento integral. En: La Ganadería en América Latina y el Caribe. Editores: Rafael Núñez, Rodolfo Ramírez, Salvador Fernández, Omar Araujo, Miguel García y Tito Efraín Diaz. Pag 4266 – 431. ALPA, FAO, IICA, BBA, México.

Gómez, J. (2002). La frecuencia del ordeño y sus ventajas. Boletín informativo, IDIAP. Estación experimental de Gualáca.

Gupta, S. K., Shinde, K. P., One, S. A., Thakur, A. y. Kumar, N. (2016). The potential impact of heat stress on production and reproduction of dairy animals: consequences and possible solutions: E-mail: sgshailesh786@gmail.com International Journal of Science, Environment ISSN 2278-3687 (O) and Technology, 5, (3), 903 – 911

Hansen, P. (2009). Effects of heat stress on mammalian reproduction. Philos. Trans. Royal Society London Biologycal Science, 364(1534): 3341–3350.

Herrera, H. J. y Barreras S. A. (2000). Estadística de la Curva de lactación en la vaca de leche y el Modelo de Wood. En; Manual de Procedimientos: Análisis estadísticos de experimentos pecuarios. Montecillo, Texcoco, Estado de México, México. Citado 09/03/2020. 

http://137.222.110.150/Calnet/vetrep5/image/int%20morpholgy%20of%20scrotum%20cut%20caudally2-dog.jpg

Inalment. (2007). Cinta para la estimación del peso corporal en bovinos cebuínos, cruzados y puros por perimetría torácica. Colombia.

Jordan, E. R., Mader, T. L., Davis, M. S. y Brown-Brand, T. (2003). Effects of Heat Stress on Reproduction. Journal of Dairy Science, 86, E104 E114.

Knegsel, A. T. M., Van den Brand, H., Graat, E. A. M., Dijkstra, J., Jorritsma, R., Decuypere, E., Tamminga, S. y Kemp, B. (2007). Dietary Energy Source in Dairy Cows in Early Lactation: Metabolites and Metabolic Hormones. Journal of Dairy Science, 90, 1477–1485.

Knight, J. (2006). Psychrometrics calculator and evaluation of the microenvironment. Based on ASHRAE, Microenvironment and Fundamentals. Citado: 26/03/2020. Disponible:https://www.kwangu.com/work/psychrometric.htm https://es.weatherspark.com/y/16739/Clima-promedio-en-Alanje-Panam%C3%A1-durante-todo-el-a%C3%B1o.

Laboratorios Bayer. (2015). Boletín informativo sobre el pecutrin vitaminado como fuente de minerales. Bayer, Alemania.

Mcnab, B. (2001). Patterns of Physical exchange with the environment. En: The Physiological Ecology of Vertebrates: A View from Energetics, Publisher: Cornell University Press, NY, USA. Pag. 8 – 23.

Mader, T. L., Davis, M. S. y Brown-Brandl, B. (2006). Environmental factors influencing heat stress in feedlot cattle1,2*University of Nebraska, Concord 68728; †Koers-Turgeon Consulting Service, Inc., Salina, KS 67401; and ‡USDA-ARS U.S. Meat Animal Research Center, Clay Center, NE 68933.

McDowell, R. E. (1985). Crossbreeding in Tropical Areas with Emphasis on Milk, Health, and Fitness. 1985 Journal of Dairy Science, 68, 2418-2435

National Oceanic and Atmosphere Administration. 2014. Environment and Psicromteric Evaluation and heat stress. Internet, Google.

National Research Council (NRC), 1989, 2001. Nutrient Requirements of Dairy Cattle. 7th ed., National Academy Press. Washington, DC. USA.

Nayeri, S. y Stothard, P. (2016). Tissues, Metabolic Pathways and Genes of Key Importance in Lactating Dairy Cattle. Citado: 23/03/2020. Disponible: Springer Science Reviews (2016) 4:49–77 DOI 10.1007/s40362-016-0040-3

Polsky, L., Marina A. y. Von Keyserlingk1, G. (2017). Invited review: Effects of heat stress on dairy cattle welfare. Animal Welfare Program, 2357 Main Mall, Faculty of Land and Food Systems, University of British Columbia, Vancouver, BC, V6T 1Z4 Canada.

Rosas, H., Gómez, J., Quintero, S. y Pimentel, N. (1985). Composición bromatológica y valor nutricional de los pastos en Panamá. En: Composición de los alimentos en Panamá y Nutricion Animal. Editorial Universitaria, Universidad de Panamá.

Rosas, H., Gómez, J. Quintero, S. y Pimentel, N. (1994). Aporte energético de las gramíneas y valor nutricional de los pastos en Panamá. En: Nutrición Animal y Composición de los alimentos en Panamá. Editorial Universitaria, Universidad de Panamá.

SAS, (1999, 2004). The analysis of variance and covariance. General Lineal Models (GLM), Statistical Analyses System, NC, USA.

Slimen, I., Belhadj, T., Najar A., Ghram y Abdrrabba, M. (2015). Heat stress effects on livestock: molecular, cellular and metabolic aspects, a review. Journal of Animal Physiology and Animal Nutrition.

Slimen, B., Najar, T., Ghram A y Abdrrabba, M. (2016). Heat stress effects on livestock: molecular, cellular and metabolic aspects. Journal of Animal Physiology and Animal Nutrition,100, 401–412.

Statistical Analysis System (SAS). (2004). General Lineal Models and the Analysis of Variance – Covariance. North Caroline State University, Raleigh, NC, USA.

St-Pierre, N. R., Cobanov, B. y Schnitkey, G. (2003). Economic Losses from Heat Stress by US Livestock Industries1 *Department of Animal Sciences Journal of Dairy Science, 86, E52–E77.

Stothard. S., Paul, A., Mandal, B., Mandal, P. Kannan, A. y. Pathak. N. N. (2004). Deriving Nutrient Requirements of Lactating Indian Cattle under Tropical Condition Using Performance and Intake Data Emanated from Feeding Trials Conducted in Different Research Institutes. Citado: 10/02/2020. Disponible: Asian-Aust. Journal of Animal Science, 2004 (6), 769-776)

Visser, R. y Wilson, R. (2006). Potencial de la producción lechera según los grupos raciales tipo leche. Horizons, CRI. Revista distribuida por Cooleche, Concepción, Bugaba, Provincia de Chiriquí.

Walker G. P., Dunshea, C. F. R. y Doyle, A. (2004). Effects of nutrition and management on the production and composition of milk fat and protein: a review. CSIRO PUBLISHING www.publish.csiro.au/journals/ajar Australian Journal of Agricultural Research, 2004, 55, 1009–1028. Australian Journal of Agricultural Research, 55, 1009–1028

West, J. W. (2003). Effects of Heat Stress on Production in Dairy Cattle. Journal of Dairy Science, 86:2131- 2144.

Wood P. D. P. (1967). Algebraic models of the lactation curve in cattle. Nature, 216, 164–165.