SciELO - Scientific Electronic Library Online

 
vol.18 issue1Report of new alleles at BG loci in Camperos chickensComparación de ordenaciones de muestras a nivel de secuencias de aminoácidos, nucleótidos y marcadores RFLP author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

  • Have no cited articlesCited by SciELO

Related links

Share


BAG. Journal of basic and applied genetics

On-line version ISSN 1852-6233

BAG, J. basic appl. genet. vol.18 no.1 Ciudad Autónoma de Buenos Aires Jan./June 2007

 

ARTÍCULOS ORIGINALES

Análisis multicaracter para estimar parámetros genéticos para crecimiento en bovinos de raza Braford en Corrientes

López, María V.1, Arias Mañotti, Adolfo A.2, Slobodzian, Ana3

1.- Universidad Nacional del Nordeste. Facultad de Cs. Exactas y Naturales y Agrimensura. ). 9 de Julio 1449. CP: 3400. Corrientes. Argentina. TE: (03783) 423126.
2.- Instituto Nacional de Tecnología Agropecuaria. Estación Experimental Agropecuaria Corrientes (INTA). C.C. 57. CP: 3400. Corrientes. Argentina. TE: (08783) 421786.
3.- Instituto Nacional de Tecnología Agropecuaria. Estación Experimental Agropecuaria Las Breñas (INTA). C.C. 38. C.P.: 3722. Las Breñas. Chaco. Argentina. TE: (03731) 460033.

Dirección de envío de correspondencia: Dr. Adolfo A. Arias Mañotti. San Juan Nº 1443. CP: 3400. Corrientes. Argentina. mvlopez@exa.unne.edu.ar; ariasm@corrientes.inta.gov.ar

RESUMEN

Se ha efectuado un estudio en La Estación Experimental Agropecuaria Corrientes del INTA sobre un rodeo Braford estabilizado 3/8 Brahman x 5/8 Hereford, que tiene como propósito estimar los parámetros genéticos y ambientales para peso al nacer (PNAC), peso ajustado a 210 días (P_AJ_210) y peso ajustado a 570 días (P_AJ_570), mediante el Modelo Animal Multicaracter, utilizando el software MTDFREML. Se procesó una base de datos de 658 registros de producción y genealogía de bovinos de raza Braford, correspondientes al período 1992-2004. Los componentes de varianza y covarianza, y parámetros genéticos y ambientales estimados son: a) Para PNAC: s2a = 6,57, s2m = 2,21, sam = -0,32, s2pe = 0,72, s2e = 7,31, s2p = 16,50, h2a = 0,40, h2m = 0,13, ram = -0,09, c2 = 0,04, e2 = 0,44; b) Para P_AJ_210: s2a = 153,89, s2m = 62,18, sam = -56,47, s2pe = 10,35, s2e = 223,05, s2p = 393,02, h2a = 0,39, h2m = 0,16, ram = -0,58, c2 = 0,02, e2 = 0,57; c) Para P_AJ_570: s2a = 292,95, s2e = 419,02, s2p = 711,97, h2a = 0,41, e2 = 0,59. Además se calcularon las correlaciones genéticas entre los caracteres analizados: PNAC-P_AJ_210: 0,83, PNAC-P_AJ_570: 0,31, P_AJ_210-P_AJ_570: 0,74. Finalmente se obtuvieron los valores genéticos de los reproductores (DEPs), cuya combinación apropiada permitirá tomar decisiones de selección para el mejoramiento de esta raza.

Palabras clave: Modelo animal; Parámetros genéticos y ambientales; Bovinos para carne; Braford; Peso al nacer; Peso ajustado a 210 días; Peso ajustado a 570 días.

SUMMARY

This study was carried out at INTA Corrientes Experimental Station with data from a stabilized Braford 3/8 Brahman x 5/8 Hereford herd. The objective was to estimate the genetic and environmental parameters for birth weight (PNAC), weight adjusted to 210 days of age (P_AJ_210), and weight adjusted to 570 days of age (P_AJ_570). We used a Multivariate Animal Model and MTDFREML software. A database of growth and genealogy data from 658 samples of Braford males and females corresponding to the period 1992-2004 was processed. Estimates of variance and covariance components and the genetic and environmental parameters were: a) For PNAC: s2a = 6,57, s2 m = 2,21, sam = -0,32, s2pe = 0,72, s2e = 7,31, s2p = 16,50, h2a = 0,40, h2m = 0,13, ram = -0,09, c2 = 0,04, e2 = 0,44; b) For P_AJ_210: s2a = 153,89, s2m = 62,18, sam = -56,47, s2pe = 10,35, s2e = 223,05, s2p = 393,02, h2a = 0,39, h2m = 0,16, ram = -0,58, c2 = 0,02, e2 = 0,57; c) For P_AJ_570: s2a = 292,95, s2e = 419,02, s2p = 711,97, h2a = 0,41, e2 = 0,59. Genetic correlations among the characters analyzed were also calculated: PNAC-P_AJ_210: 0,83, PNAC-P_AJ_570: 0,31, P_AJ_210-P_AJ_570: 0,74. Finally, the genetic value of each individual was obtained (Expected Progeny Differences - EPDs). The combination of these values will allow to improve the characteristics of this breed.

Introducción

La raza Braford se forma a partir del cruzamiento entre el Brahman principalmente (Bos indicus sp.) y el Hereford (Bos taurus sp.), obteniéndose una raza sintética que combina la rusticidad del Cebú y su adaptación al medio, con la calidad de carnes, precocidad, eficiencia, mansedumbre y fertilidad del Hereford. De esta forma se logra un biotipo que se caracteriza por su alta producción aún en condiciones ambientales adversas.
La razón de ser de la raza Braford, desde su formación en Argentina, hace ya más de cincuenta años, fue entonces la de cubrir una necesidad existente en la región subtropical donde las razas europeas tenían menor productividad como consecuencia de las rigurosas condiciones climáticas y nutricionales en que debían desenvolverse. Esta raza se ha posicionado en las regiones Litoral y Norte del país satisfaciendo una demanda cada vez más exigente de su ganadería en expansión y orientada a cubrir sus requerimientos en áreas marginales y ambientes realmente difíciles para dicha ganadería.
Tanto los terneros como los novillos son muy valorados por los compradores. Se reconocen por su buena sanidad y por su adaptación a condiciones adversas así como por producir un animal terminado que se adapta a los distintos mercados. En cuanto a la vaca, se destacan su eficiencia con altos porcentajes de preñez y parición, su facilidad de parto, su habilidad materna y su larga vida útil. La identifican con un vientre productivo aún en condiciones desfavorables. Con su población creciendo año tras año, hoy se presenta como la raza más demandada en la vasta zona de cría del Norte argentino, distribuyéndose geográficamente a lo largo y a lo ancho del país (Sitio Web de la Asociación Braford Argentina).
Este estudio se ha efectuado en el INTA EEA Corrientes sobre un rodeo Braford estabilizado 3/8 Brahman x 5/8 Hereford, y tiene como propósito estimar los parámetros genéticos y ambientales para peso al nacer (PNAC), peso ajustado a 210 días (P_AJ_210) y peso ajustado a 570 días (P_AJ_570), mediante la implementación del Modelo Animal Multicaracter, utilizando el software MTDFREML.

Materiales y Métodos

Preparación del archivo de datos para analizar los caracteres peso al nacer, peso ajustado a 210 días y peso ajustado a los 570 días.
Se procesó una base de datos de 658 registros de producción y genealogía de machos y hembras de la raza Braford, correspondientes al período 1992-2004, disponibles en el INTA EEA Corrientes, para analizar las variables peso al nacer (PNAC), peso ajustado a 210 días (P_AJ_210) y peso ajustado a los 570 días de edad (P_AJ_570). El número de toros es de 4 y el de vacas madre de 193.
Esta base de datos fue sistematizada por el Área de Producción Animal del INTA EEA Corrientes. El archivo final posee la siguiente estructura:

• I_NOM: Código del individuo.
• T_NOM: Código del padre del individuo.
• V_NOM: Código de la madre del individuo.
• I_MESNAC: Mes de nacimiento del individuo. Se consideran los meses de junio a octubre.

• GC: Código de grupo contemporáneo. Estos grupos fueron definidos como animales nacidos en el mismo año (carimbo) y del mismo sexo (GC = I_CAR + I_SEXO), obteniéndose un total de 26 grupos para los caracteres peso al nacer y peso ajustado a 210 días y 24 grupos para el carácter peso ajustado a 570 días.
I_CAR: Año de nacimiento (carimbo) del individuo.
Tiene 13 niveles.
I_SEXO: Código de sexo del individuo. Tiene dos niveles (1=macho, 2=hembra).
• V_EDAD: Edad de la madre. Está expresada en años y está tomada al momento del parto. Tiene 8 niveles, se inicia con la edad de 3 años y finaliza con la edad de 10 años.
• V_LAC_PREV: Código de lactancia previa o anterior. Indica si la madre viene de amamantar su cría el año anterior y tiene 2 niveles (1=Sí, 2=No).
• TIPO_DES: Código de tipo de destete del individuo (1=Normal, 2=Precoz).
• PNAC: Peso al nacer del individuo.
• I_EDAD_210: Edad del animal en días al momento del destete normal cercano a los seis o siete meses. Se han aplicado dos tipos de destete: normal en marzo y precoz en noviembre-diciembre. Por tanto, sólo en algunos animales el peso a los 210 días coincide con el destete. Se usa como covariable para analizar el peso ajustado a los 210 días.
• P_AJ_210: Peso del individuo en kg ajustado a 210 días. Es el peso vivo del animal tomado lo más cercano a los 6-7 meses y ajustado luego a los 210 días de edad. La fórmula de cálculo es:

donde: P210 = peso vivo del animal tomado a los 6-7 meses aproximadamente; PNAC = Peso al nacer; I_EDAD_210 = ya descripto.

• I_EDAD_2A: Edad del animal en días al momento del pesaje final cercano a los diecinueve meses. Se usa como covariable para analizar el peso ajustado a los 570 días.
• P_AJ_570: Peso del individuo ajustado a 570 días (kg). A partir de los registros de campo, se tomó el peso final cercano a los diecinueve meses y se ajustó luego a los 570 días de edad, tal que:

donde: PESO_2A = Peso final; I_EDAD_2A = ya descripto.

Modelo para el análisis multivariado de los caracteres peso al nacer, peso ajustado a 210 días y peso ajustado a los 570 días
Las estimaciones de los componentes de varianza, parámetros genéticos y ambientales, y las diferencias esperadas de la progenie (DEPs) fueron obtenidas por medio del Modelo Animal Multivariado, utilizando el software MTDFREML (Multiple Trait Derivative Free Restricted Maximum Likelihood) (Boldman y Van Vleck, 1991; Boldman et al., 1993), un conjunto de programas para estimar los componentes de (co)variancia usando el modelo animal y el algoritmo no derivativo REML (Smith y Graser, 1986 y Graser et al., 1987).
En primer lugar, fueron ajustados modelos univariados de los tres caracteres analizados, haciendo las covarianzas igual a cero, y valores supuestos para las variancias tomados de Snelling et al. Luego, se corrió el modelo multivariado usando como valores iniciales las variancias estimadas en el paso anterior y valores supuestos para las covariancias (obtenidos de Crews, D.H. Jr), manteniendo las variancias constantes. Esto permitió obtener las estimaciones de las covariancias entre los tres caracteres considerados. Finalmente, las estimaciones obtenidas en los pasos anteriores fueron utilizadas como valores iniciales de una nueva corrida del modelo multivariado, comenzando desde cero, y permitiendo que todas las variancias y covariancias varíen libremente durante las iteraciones. Este modelo fue ejecutado varias veces, continuando cada vez el procesamiento anterior, hasta que se alcanzó el mínimo global de -2 log L (cuando éste no se modificó en la tercera posición decimal luego de sucesivas corridas). La convergencia se alcanzó después de 6075 iteraciones.

Se estimaron los efectos aleatorios de los valores genéticos aditivos directos y maternos (sólo para los caracteres peso al nacer y peso ajustado a 210 días), efectos aleatorios del ambiente permanente de las madres (sólo para los caracteres peso al nacer y peso ajustado a 210 días), efectos aleatorios residuales y los siguientes efectos fijos:

Peso al nacer: mes de nacimiento, grupo contemporáneo, edad de la madre, lactancia previa de la madre.
Peso ajustado a 210 días: mes de nacimiento, grupo contemporáneo, edad de la madre, lactancia previa de la madre, tipo de destete, edad en días al momento del pesaje de marzo que es cercano a los 210 días (covariable).
Peso ajustado a 570 días: mes de nacimiento, grupo contemporáneo), edad de la madre, tipo de destete, edad en días al momento del pesaje final en mayo (covariable).

El modelo matemático propuesto es:

donde:
Yi = vector de observaciones para el i-ésimo carácter.
bi = vector de efectos fijos para el i-ésimo carácter.
ui = vector de efectos genéticos aditivos directos de los animales,
mi = vector efectos genéticos aditivos maternos de los animales,
pei = vector de efectos ambientales permanentes de las madres,
ei = vector de efectos aleatorios residuales,
Xi, Zi, Wi y Si son matrices de incidencia que relacionan registros del i-ésimo carácter a los efectos fijos, genéticos directos, genéticos maternos y ambientales permanentes respectivamente.

La estructura de (co)variancias supuesta es:

donde s2ai con i=1,..,3 son las variancias genéticas aditivas directas para los caracteres 1 (peso al nacer), 2 (peso ajustado a 210 días) y 3 (peso a 570 días), s2mi con i=1,..,2 son las variancias genéticas aditivas maternas para los caracteres 1 y 2, s2pei con i=1,..,2 son las variancias de ambiente permanente para los caracteres 1 y 2, con i=1,..,3 son las variancias residuales para los caracteres 1, 2 y 3, saiaj con i=1,.., 3; j=1,..,3 son las covariancias entre los efectos genéticos aditivos directos de los caracteres 1, 2 y 3, saimj y sajdi con i=1,.., 3; j=1,..,2 son las covariancias entre los efectos genéticos aditivos directos de los caracteres 1, 2 y 3 y los efectos genéticos aditivos maternos de los caracteres 1 y 2, speipej con i=1,.., 2; j=1,..,2 son las covariancias entre los efectos de ambiente permanente de los caracteres 1 y 2, seiej con i=1,.., 3; j=1,..,3 son las covariancias entre los efectos residuales de los caracteres 1, 2 y 3, A es la matriz de parentesco entre los individuos, Ipe y Ie son matrices identidad.

Resultados y Discusión

Principales estadísticas descriptivas
El Cuadro 1 presenta las estadísticas descriptivas generales para los tres caracteres analizados, informando el número de observaciones (N), valor medio, desvío standard (SD), coeficiente de variación (CV), valor mínimo y valor máximo.

Efectos aleatorios
Las estimaciones de los componentes de varianza y covarianza, parámetros genéticos y ambientales se presentan en el Cuadro 2, donde: σ2a = varianza genética aditiva directa, σ2 m= varianza genética aditiva materna, σam = covarianza entre efectos genéticos directo y materno, σ2pe = varianza de ambiente permanente, σ2e = varianza residual, σ2p = varianza fenotípica, h2a = heredabilidad para efectos genéticos directos, h2m = heredabilidad para efectos genéticos maternos, ram = correlación genética entre los efectos directos y maternos, e2 = proporción de la varianza fenotípica total debida a la varianza residual y c2 = proporción de la varianza fenotípica total debida a la varianza ambiental permanente.


Figura 1: Proporciones de los componentes de varianza estimados para peso al nacer
Figure 1: Proportions of estimated variance components for birth weight


Figura 2: Proporciones de los componentes de varianza estimados para peso ajustado a 210 días
Figure 2: Proportions of estimated variance components for weight adjusted to 210 days


Figura 3: Proporciones de los componentes de varianza estimados para peso ajustado a 570 días
Figure 3: Proportions of estimated variance components for weight adjusted to 570 days

Las correlaciones genéticas y ambientales entre los tres caracteres analizados se presentan en el Cuadro 3.

Se observa que el peso al nacer y el peso ajustado a los 210 días presentan un alto grado de correlación, al igual que el peso ajustado a los 210 días con el peso ajustado a los 570 días. En ambos casos esto indica una importante respuesta correlacionada. En el segundo caso, la alta correlación genética implica que al seleccionar por el peso ajustado a los 570 días habrá una respuesta positiva en el peso ajustado a los 210 días, sin incrementar la habilidad materna. En cuanto a los caracteres peso al nacer y peso ajustado a los 570 días, se observa que, si bien se correlacionan positivamente, el valor es bajo, indicando que al seleccionar para peso ajustado a los 570 días, los incrementos en el peso al nacer no serán importantes. No obstante, como la correlación es positiva, se debe esperar que los animales de mayor peso a 570 días tiendan a producir crías más pesadas al nacer.

Valores genéticos predichos
Las DEPs calculados a partir de los valores genéticos predichos (BVs) fueron graficadas en función del número de observaciones, dando una idea general de su distribución en la población analizada.
La observación de los gráficos muestra claramente que la distribución de los DEPs se aproxima a la normal. En las Figuras 4 y 5 se muestra la representación gráfica de la distribución de los DEPs para los efectos directo y materno en el carácter peso al nacer, en las Figuras 6 y 7 se puede ver la distribución de los DEPs para los efectos directo y materno en el carácter peso ajustado a 210 días, y en la Figura 8 se observan los DEPs para el efecto directo en el peso ajustado a 570 días.


Figura 4: Distribución de los DEPs para el efecto directo en el peso al nacer
Figure 4: Distribution of DEPs for direct effect for birth weight


Figura 5: Distribución de los DEPs para el efecto materno en el peso al nacer
Figure 5: Distribution of DEPs for maternal effect for birth weight


Figura 6: Distribución de los DEPs para el efecto directo en el peso ajustado a 210 días
Figure 6: Distribution of DEPs for direct effect for weight adjusted to 210 days


Figura 7: Distribución de los DEPs para el efecto materno en el peso ajustado a 210 días
Figure 7: Distribution of DEPs for maternal effect for weight adjusted to 210 days


Figura 8: Distribución de los DEPs para el efecto directo en el peso ajustado a 570 días
Figure 8: Distribution of DEPs for direct effect for weight adjusted to 570 days

Finalmente, se ha elaborado un Resumen de Padres, en el cual se han incluido los DEPs calculados para los caracteres analizados, obteniéndose un resumen de padres, que contiene los siguientes campos:

• Código de padre: Es el código que recibe cada toro padre de plantel.
• Nro. de hijos: Número de hijos del padre evaluado.

Para cada carácter se incluye, tanto para el efecto directo como para el efecto materno:

• DEP: Diferencia esperada en la progenie de cada toro. Se expresa en las unidades del carácter.
• Precisión: Indica la confiabilidad del valor del DEP.

Conclusiones

Se ha procesado una base de datos con información de producción y genealogía de bovinos de raza Braford disponible en el INTA EEA Corrientes, mediante el Modelo Animal Multivariado, obteniéndose estimaciones de los factores ambientales y genéticos que influencian los caracteres peso al nacer, peso ajustado a 210 días y peso ajustado a los 570 días. Se estimaron los componentes de varianza, parámetros genéticos y ambientales, y además se obtuvo información de predicción del valor genético de los reproductores (DEPs).
Mediante la combinación apropiada de la información proporcionada por los DEPs para los tres caracteres analizados, será posible identificar reproductores mejoradores en cada una de esas características de crecimiento lo cual permitirá tomar decisiones de selección para el mejoramiento de esta raza y de los rodeos individuales que la componen.

Bibliografía

1. Boldman, K.G. y Van Vleck, L.D. (1991). Derivative- Free Restricted Maximum Likelihood Estimation in Animal Model with a Sparse Matrix Solver. J. Dairy Sci, v.74.        [ Links ]

2. Boldman, K.G., Kriese, L.A., Van Vleck, L.D. y Kachman, S.D. (1993). A Manual for Use of MTDFREML: A Set of Programs to Obtain Estimates of Variances and Covariances. USDAARS.        [ Links ]

3. Crews, D.H Jr. Multiple Trait Selection for maternal Productivity: The Hereford Maternal Productivity Index. Agriculture and Agri-Food Canada Research Centre. Lethbridge. Alberta. Canadá.        [ Links ]

4. Graser, H.U., Smith, S.P. y Tier, B. (1987). A derivative-free approach for estimating variance components en animal models by restricted maximum likehood. J. Anim. Sci., v.64.        [ Links ]

5. Sitio Web de la Asociación Braford Argentina. Pasaje Rodolfo Rivarola 111, piso 5, oficina 15. CP:1015. Capital Federal. http://www.braford.org.ar/frame.htm.        [ Links ]

6. Smith, S.P. y Graser, H.U. (1986). Estimating variance components in a class of mixed models by restricted maximum likelihood. J. Dairy Sci., v.69.        [ Links ]

7. Snelling, W.M., Macneil, M.D., Kress, D.D., Anderson, D.C. y Tess, M.W. Factors Influencing Genetic Evaluations of Linebred Hereford Cattle in Diverse Environments. Ft. Keogh LARRL, ARS, USDA, Miles City, MT 59301. Department of Animal Science, Montana State University, Bozeman 59717. Northern Agric. Research Center, Havre, MT 59501        [ Links ]

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License