Livestock Research for Rural Development 18 (4) 2006 Guidelines to authors LRRD News

Citation of this paper

Estimativas de factores de corrección para duración de la lactancia, edad y época de parto en búfalas de la Costa Atlántica Colombiana

A Gutiérrez-Valencia, N Hurtado-Lugo y M Cerón-Muñoz

GRICA, Grupo de Investigación en Ciencias Agrarias, Facultad de Ciencias Agrarias, Universidad de Antioquia, Colombia
mceronm@agronica.udea.edu.co


Resumen

El propósito de este trabajo fue estimar factores multiplicativos de ajuste de la producción de leche a 240 y 270 días de lactación, edad adulta y época de parto en búfalas de la Costa Atlántica colombiana. Fueron empleados para la producción de leche a 240 días 2548 lactancias de 1059 búfalas, para producción de leche a 270 se utilizaron 1736 lactancias de 959 búfalas y para la estimación de los factores a edad adulta se utilizaron 3153 lactancias de 1150 búfalas. Se emplearon los registros de producción de la Hacienda Altamar, ubicada en el Municipio de Buenavista, en el Departamento de Córdoba. Los análisis estadísticos fueron realizados mediante la utilización de un modelo animal y el método de máxima verosimilitud restricta libre de derivadas. Los efectos fijos considerados en el modelo para el ajuste a los 240 y 270 días en leche fueron: grupo contemporáneo (año- época), clases de duración de la lactancia (meses de edad) y la covariable edad de la búfala al parto (efecto lineal y cuadrático). Para el ajuste de producción de leche a edad adulta fueron considerados como efectos fijos: época de parto y clases de edad de la búfala. Los efectos aleatorios para todos los análisis fueron: genético aditivo, de ambiente permanente y residual. Los factores de ajuste para 240 y 270 fueron eficientes con correlaciones entre producciones reales y ajustadas de 0.92. Los factores de ajuste para edad adulta fueron eficientes, ya que removieron la variación causada por el efecto de la edad y la época de parto.

Palabras Claves: Búfalas, factores de ajuste, mejoramiento genético, modelo animal, producción de leche



Estimates of correction factors for lactation length, age and calving season in buffaloes of Colombian Atlantic Coast

Abstract

The purpose of this research was to estimate multiplicative adjustment factors for the milk yield at 240 and 270 days in milk (DIM), age of maturity and season of calving,  in buffaloes in the Atlantic Colombian Coast. The number of lactations used for milk yield at 240 DIM was 2548 from 1059 buffaloes and 270 DIM was for 1736 lactations of 959 buffaloes. For effect of age of maturity, 3153 lactations of 1150 buffaloes of the Altamar farm from municipality of Buenavista in the department of Córdoba were used. The statistical analyses were done by the application of the animal model and the derivate-free restricted maximum likelihood method. For the adjustment of milk yield at 240 and 270 DIM, the fixed effects considered in the model were: contemporary group (year-season), and classes of lactation length (months of age) and the variable age of the buffaloes at calving (linear and quadratic effect). For milk yield at mature age the fixed effects considered in the model were: parity (year of calving), age and season of calving and lactation length (months of age). The random effects for all analyses were: genetic additive, permanent environment and residual effects. The adjustment factors for 240 and 270 DIM were significant with correlations between real and adjustment productions of 0.92. The adjustment factors of age of buffalo and parity season were significant.

Key words: Adjustment factors, animal model, buffalos, genetic improvement, milk yield


Introducción

La selección aplicada sobre una o varias características productivas y reproductivas implica la escogencia de los individuos genéticamente superiores. Sin embargo, existen varios factores ambientales que influyen en el desempeño productivo y reproductivo de los animales (Schmidt y Van Vleck 1974. En virtud de esto, cuando es necesario comparar individuos por medio de registros de producción se requiere corregir los efectos ambientales que alteran el potencial productivo.

Algunos factores ambientales que influyen directamente sobre la producción de leche, pueden ser removidos utilizando el ajuste previo de los datos o pueden ser incluidos en los modelos estadísticos de evaluación genética. Los principales efectos ambientales controlados con factores de ajustes para producción de leche por lactancia, envuelven otras características de desempeño de la vaca como: duración del periodo seco, duración del periodo parto-concepción; efectos causados por el manejo ó nivel de producción de las haciendas como: sistemas de alimentación, sistemas y número de ordeños, efectos causados por el ciclo de vida del animal como: edad y número de partos de la hembra y circunstancias ambientales, como época de parto (Searle 1962; Miller 1973; Schmidt y Van Vleck 1974; Warwick y Legates 1980)

En Colombia no se tienen criterios definidos para la selección de búfalos lecheros y no se cuenta con estudios de curvas de lactancia y factores ambientales que influyen en la producción de leche de los animales, por lo tanto, el propósito de este trabajo fue estimar factores de corrección para la producción de leche de búfalas en la región de la costa atlántica colombiana a los 240 y 270 días de lactancia, edad adulta y época de parto.


Materiales y Métodos

Se utilizaron los registros de producción de leche de búfalas de la hacienda Altamar desde el año de 1999 hasta el 2004. El hato que se encuentra ubicado en el Municipio de Buenavista (Córdoba), y esta compuesto por búfalos mestizos con alto linaje de la raza Murrah, criados en pastoreo, y ordeñados manualmente en la mañana y en la tarde. La suplementación de la hacienda se lleva a cabo en la época de verano (diciembre a marzo), es realizada con heno de pasto Climacuna (Dichanthium annulatum).

Para el cálculo de las producciones hasta los 240 y 270 días en leche (PL240 y PL270), se utilizó el procedimiento recomendado por el Comité Internacional de registro animal (ICAR, 2002), según la expresión [1]

 

En que:

PL = producción acumulada de leche hasta 240 ó 270 días (PL240 ó PL270 días)
Ci = Producción de leche en el i-ésimo control, i Î [1, n];
D1 = Días al primer control;
Di = Días en leche, i Î [2, n];
Ds = Dia de secado.
C1 = Producción de leche en el primer control
Dn = Días al último control
Cn = Producción de leche en el último control

Las lactancias con mas de 240 días se agruparon al azar en 15 clases de días en leche (de 100 a 249 días con intervalos de 10 días) y se recalculó la producción de leche acumulada hasta el día de corte de cada clase (PD240), utilizando la expresión [2].

En que:

Pd =Producción de leche total al corte al azar;
C1 =Producción de leche en el primer control;
Ci = Producción de leche en el i-ésimo control, [1-u];
Cu =Producción de leche del último control antes del corte al azar;
D1 = Intervalo entre el dato del parto y del primer control;
Di = Dias en leche, i Î [2≤ u];
Dd= Días de leche donde se hizo el corte al azar;
Du= Días en leche del último control, donde u ≤ d

De igual forma, se recalculó la producción de leche hasta el día de corte (PD270) de 18 clases, para lactancias con más de 270 días en lactancia (Tabla 1). Para la PD240 y PD270 el número de lactancias asignadas en cada clase varió entre 148 y 213 y 74 y 118 para un total de 2548 y 1736 lactancias de 1059 y 959 búfalas, respectivamente. (Tabla 1).

Tabla 1. Distribución de las lactaciones de búfalas en la hacienda Altamar en la Costa Atlántica Colombiana, en clases de duración de la lactación entre 100-240 días y 100-270

240 días en leche

270 días en leche

Clases

Días en lactancia

No

Días en lactancia

No

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

100-109

110-119

120-129

130-139

140-149

150-159

160-169

170-179

180-189

190-199

200-209

210-219

220-229

230-239

240-249

 

179

184

157

161

152

176

171

161

213

148

162

161

172

170

181

460±128

508±122

571±147

587±126

642±143

676±152

716±160

758±182

799±165

846±182

899±206

926±186

957±219

1030±214

1034±213

100-109

110-119

120-129

130-139

140-149

150-159

160-169

170-179

180-189

190-199

200-209

210-219

220-229

230-239

240-249

250-259

260-269

270-279

118

107

74

85

95

77

93

107

92

98

112

86

98

93

85

100

109

107

452±104

505±112

544±150

595±148

644±156

734±188

716±176

784±179

796±205

866±190

899±190

945±201

985±182

985±206

1051±241

1107±227

1098±229

1158±242

No= número de lactaciones, = promedio y desviación estándar de la producción de leche

Para la estimativa de las soluciones de cada clase para D240 y D270 se utilizó un modelo animal, con medidas repetidas. Fue utilizado el programa computacional MTDFREML (Boldman et al 1995), se consideró una matriz de parentesco con 5498 animales de ocho generaciones. Las hembras con registros fueron hijas de 149 machos reproductores. Se empleo el siguiente modelo:

yijklm = bi + dj + ak + pl + eijklm,

En donde:

yijklm = variables dependientes (producciones de leche cortadas al azar).
bi = efecto fijo del i-ésimo grupo contemporáneo (año, época de parto);
d j = efecto fijo de la j-ésima clase de duración de la lactación;
ak = efecto aleatorio del k-ésimo animal, con media 0 y varianza σ2 a ;
pl = efecto aleatorio del medio ambiente permanente, con media 0 y varianza σ2 c ;
eijklm = efecto aleatorio residual, con media 0 y varianza σ2 e.

Se consideraron las épocas de parto así: verano (de Diciembre a Febrero), época de transición de verano a invierno (Marzo a Mayo), época de transición de verano a invierno (Julio a Septiembre) y época de lluvia (Septiembre a Noviembre).

Las estimativas para las clases de la duración de la lactación fueron obtenidas sumando la media general de cada clase a las soluciones de las clases de duración de las lactaciones, considerando las ecuaciones del modelo anteriormente mencionado. Estas estimativas fueron suavizadas usando un polinomio de segundo grado.

Los factores de corrección fueron estimados mediante la expresión:

Los factores de corrección estimados fueron aplicados a las producciones de leche para diferentes días en leche, obteniendo así, las producciones corregidas. Estas producciones fueron comparadas con las producciones observadas hasta los 240 y 270 días, con el objetivo de determinar la eficiencia de los factores de ajuste estimados mediante un análisis de correlación.

Para estimar los factores de ajuste a edad adulta a los 240 y 270 días en leche fueron utilizadas 3153 lactancias de 1150 búfalas respectivamente, estas producciones fueron distribuidas en clases de edad y época del parto de las búfalas.

Los factores de corrección fueron estimados mediante la utilización del siguiente modelo:

Yijklm= Hi+Ij +ak+pl+eijklm

Donde:

Yijklm= Producción de leche hasta 240 y 270 días.
Hi = Efecto fijo del i-ésimo grupo contemporáneo (año del parto)
Ij = Efecto fijo de la j-ésima clase de edad - época de parto.
ak = Efecto aleatorio del k -ésimo animal, con media 0 e varianza σ2a ;
pl= Efecto aleatorio de ambiente permanente de la búfala, con media 0 de varianza σ2c, y
eijklm = Efecto aleatorio residual, con media 0 e varianza σ2e

Los factores de ajuste fueron calculados mediante la expresión utilizada por Iloeje (1980), definida como:

Fj = Mjb / [Mjb -(Ajb -Aj)],

Donde:

Fj es el factor de ajuste multiplicativo para la clase de edad - estación de parto j.
Mjb es la media de producción de leche de la clase considerada como base búfalas a edad adulta.
Ajb es la solución obtenida en la clase de duración de la lactación escogida como base y
Aj es la solución en la clase de edad - y estación de parto j.

Los factores multiplicativos Fj fueron suavizados dentro de cada estación del año y utilizando el método descrito por Wiggans (1981).

Las producciones de leche ajustadas y no ajustadas a edad adulta fueron comparadas por la reducción de la variación asociada a los efectos de la edad y la estación de parto utilizando análisis de varianza, incluyendo efectos época - edad (grupo época y edad de la búfala), año de parto y búfala con el propósito de determinar la reducción de la variación asociada a la edad del animal (Searle y Henderson 1959; Miller et al 1968; Searle 1962).


Resultados

Los factores de ajuste estimados para los 240 y 270 días en leche se encuentran en la tabla 2, se observó que los factores de ajuste a los 100 días fueron de 2.404 y de 2.465 respectivamente.

Tabla 2.  Factores de corrección multiplicativos para el ajuste de producción de leche a 240 y 270 días en lactancia en búfalas de la costa atlántica colombiana

Días en lactancia

Factor de ajuste

Para 240 días

Para 270 días

100-109

110-119

120-129

130-139

140-149

150-159

160-169

170-179

180-189

190-199

200-209

210-219

220-229

230-239

240-249

250-259

260-269

270-279

2.404

2.163

1.969

1.809

1.675

1.561

1.463

1.378

1.304

1.238

1.180

1.128

1.081

1.039

1.000

2.465

2.237

2.051

1.896

1.766

1.654

1.558

1.474

1.400

1.335

1.277

1.225

1.178

1.135

1.097

1.062

1.029

1.000

Las medias de producción de leche hasta y ajustada a los 240 días fueron de 1025 ± 213 y 1063 ± 243 kg, respectivamente, con una correlación entre las dos variables de 0.92. Las producciones de leche hasta y ajustada a los 270 días fueron de 1152 ± 229 y 1143 ± 259 kg, respectivamente, con una correlación de 0.92.

Lo anterior indica que los factores de ajuste fueron eficientes y confiables y pueden ser utilizados para ajustar lactancias incompletas. Sin embargo, la eficiencia del ajuste es menor en hembras con pocos días en leche, con correlaciones entre las producciones de leche hasta y ajustadas de 0.85 para 240 días y 0.79 para 270 días en la clase 1, pero esta eficiencia va mejorando a medida que la hembra tiene mas días en leche, con correlaciones entre las producciones de leche hasta y ajustadas de 0.99 para 240 días y 0.99 para 270 días en la clase 15 y 18, respectivamente.

En las tablas 3 y 4 se presentan los factores de ajuste para edad adulta. En hembras de 20 a 31 meses de edad, el factor de ajuste a edad adulta varió de 1.129 a 1.215 en la producción de leche hasta los 240 días, indicando que hembras de esta edad producen entre el 79 y el 83% de una hembra adulta, dependiendo de la época de parto (Tabla 3).

Tabla 3. Factores de corrección multiplicativos1 para el ajuste a edad adulta de la producción de leche hasta 240 días de búfalas de la costa atlántica colombiana

Edad de la búfala

(meses)

Factores de ajuste

Diciembre-Febrero

Julio-Septiembre

Octubre-Noviembre

20-31

32-38

39-44

45-51

52-59

60-67

68-75

76-85

86-95

96-108

109-125

126-299

1.129

1.093

1.065

1.044

1.029

1.020

1.016

1.018

1.025

1.038

1.056

1.082

1.103

1.103

1.071

1.047

1.028

1.014

1.005

1.000

1.000

1.004

1.012

1.026

1.215

1.169

1.133

1.107

1.088

1.077

1.072

1.074

1.082

1.098

1.121

1.152

1Clase Base: búfalas entre 86 y 95 meses de edad y ocurridos en el mes de Julio-Septiembre

En el caso de la producción de leche hasta 270 días, las hembras de 20 a 31 meses producen entre un 83 y 91% de una hembra adulta (Tabla 4).

Tabla 4.  Factores de corrección multiplicativos1 para el ajuste a edad adulta de la producción de leche hasta 270 días de búfalas de la costa atlántica colombiana

Edad de la búfala

(meses)

Factores de ajuste

Diciembre-Febrero

Julio-Septiembre

Octubre-Noviembre

20-31

32-38

39-44

45-51

52-59

60-67

68-75

76-85

86-95

96-108

109-125

126-299

1.094

1.066

1.045

1.030

1.019

1.013

1.011

1.013

1.021

1.033

1.049

1.072

1.118

1.084

1.057

1.035

1.019

1.008

1.001

0.998

1.000

1.006

1.016

1.031

1.167

1.131

1.102

1.081

1.067

1.058

1.056

1.058

1.067

1.082

1.103

1.132

1Clase Base: búfalas entre 86 y 95 meses de edad y ocurridos en el mes de Julio-Septiembre

Estos factores fueron aplicados a las de producciones de leche hasta los 240 días, donde las producciones de leche sin ajustar y ajustada a edad adulta fueron de 954.18±214.93 y 1009.306 ±224.76, respectivamente (Tabla 5).

Tabla 5. Análisis de varianza para lactaciones hasta 240 días ajustadas y no ajustadas a edad adulta de búfalas de la costa atlántica colombiana

Fuente de Variación

GL

Producción de leche

Producción de leche ajustada a edad adulta

CM

Prob F

CM

Prob F

Edad de la búfala y época de parto

35

104086

0.0036*

62026

0.24n.s.

Año de parto

7

743188

0.0001*

811478

0.0001*

Búfala

1062

110714

0.0001*

123751

0.0001*

Error

1713

46197

 

50517

 

Total

2817

 

 

 

 

Coeficiente de determinación

64.8

64.3

Promedio de leche (kg)

954

1009

Coeficiente de Variación

22.5

22.7

GL= grados de libertad, CM=cuadrado medio, n.s.=diferencias estadísticas no significativas y *=diferencias estadísticas altamente significativas (P≤0.01)

En los análisis de varianza para las producciones de leche sin ajustar se encontró que el efecto edad-época de parto fue altamente significativo (P≤0.01), y en la producción de leche ajusta fue no significativo (P≥0.05), indicando que los factores de ajuste removieron el efecto del grupo de edad de la búfala y época de parto en la producción de leche hasta los 240 días (Tabla 5).

También se aplicaron los factores de ajuste a las producciones de leche hasta los 270 días, donde las producciones de leche sin ajustar y ajustada a edad adulta fueron de 1027.874±237.01 y 1069.59 ±244.86, respectivamente (Tabla 6).

Tabla 6.  Análisis de varianza para lactaciones hasta 270 días ajustadas y no ajustadas a edad adulta de búfalas de la costa atlántica colombiana

Fuente de Variación

GL

Producción de leche

Producción de leche ajustada a edad adulta

CM

Prob F

CM

Prob F

Edad de la búfala y época de parto

35

99733

0.0001*

69433

0.17n.s.

Año de parto

7

997945

0.0001*

1060714

0.0001*

Búfala

1062

133129

0.0001*

144649

0.0001*

Error

1713

46197

 

59957

 

Total

2817

 

 

 

 

Coeficiente de determinación

64.7

64.6

Promedio de leche (kg)

1028

1070

Coeficiente de Variación

23.1

22.9

En los análisis de varianza para las producciones de leche sin ajustar se encontró que el efecto edad-época de parto fue altamente significativo (P≤0.01), y en la producción de leche ajusta fue no significativo (P≥0.05), indicando que los factores de ajuste removieron el efecto del grupo de edad de la búfala y época de parto (Tabla 6), en la producción de leche hasta los 270 días.


Discusión

Las medias de producción de leche fueron bajas, comparadas con las encontradas en Brasil, donde la producción de leche corregida para los 240 y 270 días fueron de 1535 ± 566 y 1712 ± 347, respectivamente (Tonhati et al 2004), las encontradas en Italia, donde la producción de leche corregida para los 270 días fue de 2286± 492 kg (Rosati y Van Vleck 2002), en Pakistán, donde la producción de leche ajustada a 308 días de búfalos Nilli-Ravi fue de 2139 kg (Khan y Chaudry 2001). La mayor producción de leche en Italia, Brasil y otros países europeos se debe probablemente a las diferencias en la alimentación animal, duración de la lactancia y la constitución genética.

Se encontró que al aplicar los factores de ajuste para edad adulta se pudo remover el efecto del grupo de edad de la búfala. La mayor producción de leche fue alcanzada por búfalas que parieron en la estación que comprendía los meses de julio a septiembre, con una edad de 86 a 95 meses (Esta fue considerada la edad adulta). No se tuvo en cuenta la estación de Marzo a Junio en las tablas de factores de corrección para edad de la búfala (Tablas 3 y 4) por no existir datos suficientes para un análisis confiable, debido a que las búfalas tienden a parir en determinadas épocas del año.

No se encontraron investigaciones realizadas en los reportes de la literatura para la comparación y discusión de factores de ajuste a edad adulta en búfalos. Sin embargo, en bovinos se han encontrado edades adultas en ganado Holstein en Colombia entre los 76 a 86 meses de edad (Stanton 1990; Martinez 1993; Cerón-Muñoz et al 2003).


Conclusiónes


Agradecimientos

Este Proyecto fue financiado por el "Fondo para apoyar trabajos de Grado en los programas de Pregrado" (Vicerrectoría de Investigación de la Universidad de Antioquia). Los autores agradecen a la hacienda Altamar S.A, por facilitar la información para el desarrollo de este trabajo y al centro de investigaciones de la Facultad de Ciencias Agrarias (CIAG) por el apoyo.


Bibliografía

Boldman K, Kriese L, Van Vleck L, Van Tassell C and Kachman S 1995 A manual for use of MTDFREML. A set of programs to obtain estimates of variances and covariances. Lincon: Department of Agriculture Research Service. 120 p.

Cerón Muñoz M, Tonhati H, Costa C, Solarte C y Benavides O 2003 Factores de ajuste para producción de leche en bovinos Holstein colombiano. Revista Colombiana de Ciencias Pecuarias. Volume 16, pp. 26-32. http://kogi.udea.edu.co/revista/16/16-1-4.pdf

ICAR 2002 Section 2.4 Guidelines for buffalo milk recording for low to medium and medium to high input production systems. In: International agreement of recording practices. International Commitee in Animal Recording www.icar.org. Accesado el 20 de mayo de 2005.
http://216.239.59.104/search?q=cache:sZwwoIKV6CoJ:www.icar.org/docs/technical_series/tec_series_4_bled_buffalo.pdf+Guidelines+for+buffalo+milk+recording+for+low+to+medium+and+medium+to+high+input+production+systems&hl=fr&client=firefox-a

Iloeje M 1980 Age-adjustment factors for Alpine, LaMancha, Nubian ,Saanen, and Toggenburg dairy goats, Journal of Dairy Science. Volume 63, pp 1309-1316.

Khan M and Chaudry H 2001 Prediction of lactaction yield from last- record day and average daily yield in Nili- Ravi Buffaloes. Bulletin of the FAO Inter-regional Cooperative Research Network on Buffaloes, Volume 15, pp 7-10.

Martinez M 1993 Determinación de factores de ajuste para ganado Holstein en la Sabana de Bogotá. Tesis de Maestría, Facultad de Ciencias Pecuárias, Universidad Nacional de Colombia, Bogotá. 86 p.

Miller R 1973 A recent study of age adjustment. Journal of Dairy Science, Volume 56, pp 952-959.

Miller R, Mc Daniel B and Myowman R 1968 Effects of errors in the edge adjustment of first lactactions. Journal of Dairy Sience, Volume 51, pp 378-384,

Rosati A and Van Vleck L 2002 Estimation of genetic parameters for milk, fat, protein and mozzarella cheese yield in the Italian river buffaloes population. Livestock Production Science, Volume 74, pp185-190.

Schmidt G and Van Vleck L 1974 Bases científicas de la producción lechera. Edición Acribia; España, 583 p.

Searle S 1962 Age and herd effects in New Zealand dairy cow records. Journal of Dairy Science Volume 45, pp 82-85.

Searle S and Henderson C 1959 Establishing age-correction factors related to the level of herd yield. Journal of Dairy Science, Volume 42, pp 824-835.

Stanton T 1990 Investigation of genotype by environment interaction for Holstein milk yield in Colombia, Mexico y Puerto Rico. Tesis de Doctorado, Department of Animal Science, Cornell University., Ithaca, 623p.

Tonhati H, Muñoz-Cerón M, Duarte J, Reichert R, Oliveira J and Lima a 2004 Estimates of correction factors for lactation length and genetic parameters for milk yield in buffaloes. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, Volume 56, pp 251-257 http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-09352004000200016&lng=en&nrm=iso

Warwick E and Legates J 1980 Cria y mejora del ganado. 3ra edición México: McGraw_Hill,. 623 p.

Wiggans G 1981 Smoothed Age- Season Adjustments Factors for Dairy Goat Lactation Milk and fat records. Journal of Dairy Science, Volume 64, pp 350-352.


Received 29 June 2005; Accepted 26 January 2006; Published 11 April 2006

Go to top