693 Revista MVZ Córdoba Rev. MVZ Córdoba 0122-0268 1909-0544 Universidad de Córdoba Colombia revistamvz@gmail.com 69360025017 https://doi.org/10.21897/rmvz.1662 Artículos Originales Evaluación genética para características de peso en ganado Brahman comercial Genetic evaluation for weight traits in commercial Brahman cattle http://orcid.org/0000-0002-2981-6287 Bedoya O Christian cfbo09@gmail.com http://orcid.org/0000-0003-3186-0778 Alzate R Jovany cfbo09@gmail.com http://orcid.org/0000-0001-5339-5913 Ángel C Juan cangelca@gmail.com http://orcid.org/0000-0002-3955-484X Escobar R Carlos carlos99992@gmail.com http://orcid.org/0000-0003-3400-5208 Calvo C Samir samirmejoramiento@gmail.com Universidad Católica de Oriente, Facultad de Ciencias Agropecuarias, Grupo de investigación en agronomía y zootecnia GIAZ, Rionegro–Antioquia. Universidad Católica de Oriente Colombia Universidad Católica de Oriente, Facultad de Ciencias Agropecuarias, Grupo de investigación en agronomía y zootecnia GIAZ, Rionegro–Antioquia. Universidad Católica de Oriente Colombia Universidad Católica de Oriente, Facultad de Ciencias Agropecuarias, Grupo de investigación en agronomía y zootecnia GIAZ, Rionegro–Antioquia. Universidad Católica de Oriente Colombia Universidad de Antioquia, Facultad de Ciencias Agrarias, Grupo de Investigación en Ciencias Agrarias-GRICA, Medellín-Antioquia. Universidad de Antioquia Colombia Universidad Católica de Oriente, Facultad de Ciencias Agropecuarias, Grupo de investigación en agronomía y zootecnia GIAZ, Rionegro–Antioquia. Universidad Católica de Oriente Colombia Mayo-Agosto 2019 24 2 7225 7230 01 05 2018 01 11 2018 01 05 2019 Resumen

Objetivo. Evaluar genéticamente las características de peso al nacimiento y peso al destete en una población de ganado Brahman comercial. Materiales y métodos. Se analizaron 1.015 registros de peso al nacimiento (PN) y peso al destete ajustado a los 270 días (PA270) de crías entre los años 2002 y 2014 pertenecientes a la hacienda La Envidia, ubicada en el municipio de Planeta Rica, Córdoba. Se utilizó un modelo animal bi-carácter que incluyó efectos fijos de sexo, mes, año y edad al destete, efecto genético aditivo y el efecto materno, estimando los parámetros genéticos con el software MTDF-REML. Resultados. Se encontró un promedio para PN de 29.62±3.13kg y coeficiente de variación (CV) del 10.6%. La media para PA270 fue de 216.713±26.97 kg y CV 12.5%. Los estimados de heredabilidad directa (h2d) fueron de 0.44±0.14 y 0.39±0.12 para PN y PA270, respectivamente. La heredabilidad materna (h2m) fue de 0.09±0.092 y 0.42±0.087 para PN y PA270, respectivamente. La correlación genética y fenotípica estimada entre las características evaluadas fue de 0.11 y 0.13, respectivamente; con progreso genético por año para PN de -2.02 kg y 7.81 kg para PA270. Conclusiones. Las estimaciones de heredabilidad indican la existencia de variabilidad genética aditiva directa, evidenciando la oportunidad de ganancia genética por medio de selección. Sin embargo, las características de peso al nacer y peso del destete presentan gran influencia del ambiente, por lo cual se debe poner especial atención al manejo y la nutrición de las hembras durante la gestación y de los terneros desde el nacimiento hasta el destete.

Abstract

Objective. Evaluate genetically the traits of birth weight and weaning weight in commercial Brahman cattle. Materials and methods. A total of 1.015 records of birth weight (BWT) and weaning weight adjusted to 270 days (WW270), were evaluated. Data correspond to the offspring born between 2002 and 2014 belonging to “La Envidia” farm, located in the municipality of Planeta Rica, Córdoba. A bi-character animal model that included, additive genetic effect, maternal effect and sex, month, year and age at weaning as fixed effects, was used. Genetic parameters were estimated using the MTDF-REML software. Results. The average and coefficient of variation (CV) of BWT was 29.62±3.13 kg and 10.59%, respectively. The mean and CV for WW270 was 216.71±26.97 kg and 12.45%, respectively. The estimates of direct heritability were 0.44±0.143 and 0.39±0.12 for BWT and WW270, respectively. Maternal heritability was 0.09±0.092 and 0.42±0.087 for BWT and WW270, respectively. The estimated genetic and phenotypic correlation between the evaluated characteristics was 0.11 and 0.13, respectively; with genetic progress per year of -2.02 kg for BWT and 7.81 kg for WW270. Conclusion. The estimations of heritability indicate the existence of direct additive genetic variability, evidencing the opportunity of genetic gain through selection. However, environment have a great influence on birth weight and weaning weight traits, which is why special attention should be paid to the handling and nutrition of females during pregnancy and of calves from birth to weaning.

Palabras clave Correlación genética ganado de carne heredabilidad modelo animal progreso genético Keywords Animal model beef cattle genetic correlation genetic progress heritability ¿Cómo citar? Bedoya Ospina, C., Alzate Ramirez, J., Ángel Cardona, J., Escobar R, C., & Calvo Cardona, S. (2019). Evaluación genética para características de peso en ganado Brahman comercial. Revista MVZ Córdoba, 24(2), 7225-7230. https://doi.org/10.21897/rmvz.1662 <bold>INTRODUCCIÓN</bold>

Según las perspectivas mundiales de la FAO, se prevé que la producción mundial de carne en 2018 aumente a 336 millones de toneladas en equivalente de peso en canal; es decir, un 1,7 por ciento (6 millones de toneladas) más que en 2017 (1). En Colombia la producción de carne bovina fue de aproximadamente 905 mil toneladas en canal en el año 2017; 5 mil toneladas menos que en el año anterior (2). Teniendo en cuenta esto, es importante mencionar que Colombia es un país que presenta condiciones agroecológicas tropicales favorables para la producción de forraje durante todo el año; lo que se considera una condición positiva en la producción de carne bovina (3). Específicamente el Valle del río Sinú es una zona ganadera por excelencia y se caracteriza por tener una temperatura promedio de 27°C, humedad relativa de 83% y zonificación agroecológica de bosque seco tropical (4). La mayoría de centros productivos ubicados en estas zonas, se dedican a la cría, levante y ceba de animales comerciales para carne, en su mayoría poblaciones de ganado Brahman producidos a partir de apareamientos no direccionados, ni con criterios de selección (3).

Si se quiere hacer de la ganadería de carne un negocio más rentable, es necesaria la toma de registros de pesajes al nacimiento y al destete como estrategia para mejorar estos parámetros de crecimiento de generación en generación. Estas características de peso están controladas por factores genéticos aditivos del padre y de la madre, factores genéticos aditivos maternos (producción de leche) y efecto ambiental directo, como manejo nutricional etc. (5, 6, 7). Un excelente programa de evaluación genética, selección, apareamiento y manejo de estas características de peso, hacen que las poblaciones seleccionadas generen progenies mejoradas, más precoces y más rentables para las ganaderías de carne comerciales (3).

El objetivo de esta investigación fue evaluar genéticamente las características de peso al nacimiento y peso al destete ajustado a los 270 días en un hato de Brahman comercial en el valle del Sinú-Colombia.

 <bold>MATERIALES Y MÉTODOS</bold>

Animales y condiciones agroecológicas. Se realizó la colecta de información fenotípica y genealógica de madres, padres y crías de una vacada no seleccionada de Brahman comercial durante el año 2002 hasta el año 2014. Los registros fueron tomados de la base de datos de la hacienda La Envidia ubicada en el valle del Sinú, departamento de Córdoba, Colombia. La hacienda se ubica a una altura de 18 msnm, temperatura promedio de 35°C y precipitación media anual de 1.156 mm. Se encuentra en zona de vida bosque seco tropical (bs-T), con suelos clasificados en el orden molisol, ricos en materia orgánica y con pendientes inferiores al 5% (8).

Alimentación y manejo. La finca es una ganadería de cría con hembras bovinas de la raza Brahman comercial, manejada bajo un sistema de apareamiento por monta natural e inseminación artificial a partir del año 2010. Los animales se manejaban bajo un sistema de pastoreo rotacional con potreros entre 4 y 6 hectáreas, establecidos con Braquipará (Brachiaria plantaginea), Pangola (Digitaria decumbens) y Angleton (Dichantium aristatumbenth), con 45-90 días de descanso dependiendo de la época del año; adicionalmente se les suministraba sal proteinada, grasa de sobrepaso y agua a voluntad.

Datos genealógicos y fenotípicos. Para la evaluación genética se recolectó la siguiente información: Identificación del animal, identificación del padre y de la madre, fecha de nacimiento de la cría, peso al nacimiento y sexo de la cría, peso y edad al destete. El número de registros para las características de la investigación fue de 1,015 datos en total. Los datos genealógicos fueron depurados usando el software especializado pedigree viewer (9) y los datos fenotípicos usando el programa R-project versión 3.5.1 (10).

A partir de los datos recolectados en finca, fue calculado el peso al destete ajustado a los 270 días usando la siguiente fórmula:

Dónde: PA270 es peso ajustado a los 270 días y PN es el peso al nacimiento.

Modelo estadístico. La evaluación genética se realizó usando el siguiente modelo animal bi-carácter:

Y 1 Y 2 = X 1 0 0 X 2 β 1 β 2 + Z 1 0 0 Z 2 a 1 a 2 + W 1 0 0 W 2 m 1 m 2 + e 1 e 2

Donde;

Y1, Y2: Vector de observaciones (Y1= PN yY2: PA270).

X1, X2: Matrices de incidencia que relaciona las observaciones a los efectos fijos (sexo, año de nacimiento, mes de nacimiento, año de destete, mes de destete).

Β1,β2: Soluciones para los efectos fijos de peso al nacimiento y peso al destete.

Z1, Z2: Matriz de incidencia que relaciona las observaciones de efectos genéticos aditivos directos.

a1,a2: Soluciones para los efectos aleatorios directos (valor de cría esperado).

W1, W2: Matriz de incidencia que relaciona el efecto aleatorio materno.

m1, m2: Soluciones para los efectos maternos.

e1, e2: Error estadístico.

Análisis de correlación. Para estimar la correlación fenotípica entre las características de PN y PA270, se utilizó la ecuación planteada por Falconer y Mackay (11), descrita a continuación:

r p 1 , 2 = ( r g 1 , 2 * h 1 2 * h 2 2 ) + ( r a 1 , 2 * 1 - h 1 2 * 1 - h 2 2 )

Donde:

rp1,2: Correlación fenotípica entre las características 1 y 2

rg1,2: Correlación genética entre las características 1 y 2

ra1,2: Correlación ambiental entre las características 1 y 2

h2 1: Heredabilidad de la característica 1

h2 2: Heredabilidad de la característica 2

La estimación de componentes de varianza y los parámetros genéticos se realizó por medio del método de máxima verosimilitud restricta (REML) utilizando el software MTDF-REML (12).

El progreso genético fue estimado por medio de la siguiente fórmula (13):

∆ G / a ñ o = h 2 * S I G

Donde:

ΔG: Progreso genético

h2: Heredabilidad de la característica

S: Diferencial de selección

IG: Intervalo generacional, visto como el promedio ponderado de machos y hembras de la población evaluada.

 <bold>RESULTADOS</bold>

En la matriz de parentesco se encontraron 1501 animales distribuidos en tres generaciones (83 padres, 314 madres y 795 individuos en la progenie); 66 familias paternas y 10 familias maternas. El coeficiente de endogamia promedio fue de 0.0 con 0 animales consanguíneos en la genealogía.

El promedio para PN fue de 29.62±3.13 Kg con un CV de 10.6% y valores mínimos de 20 y máximo de 40 kg. La media para PA270 fue 216,71±26.97 Kg y CV de 12.5%; con un valor mínimo y máximo de peso de 109.6 y 260.64 kg, respectivamente. En la tabla 1 se puede observar la estadística descriptiva para las características de peso al nacimiento y peso al destete en la población de Brahman evaluada.

<bold>Tabla 1</bold>. Estadística descriptiva para las características del peso al nacimiento y peso al destete en ganado Brahman comercial. Tabla 1 Tabla 1. Estadística descriptiva para las características del peso al nacimiento y peso al destete en ganado Brahman comercial.
Características N X ± D.E C.V Min Max
PN (kg) 1015 29.62 ±3.13 10.6 20 40
PA270(kg) 1015 216.71 ±26.97 12.5 109.6 260.6
N= Número de muestras, X= Media, D.E= Desviación estándar, C.V= Coeficiente de variación, Max= Peso máximo, Min= Peso mínimo.

Los componentes de varianza, parámetros genéticos y fenotípicos para las características PN y PD encontrados en ganado Brahman comercial se pueden observar en la tabla 2 y figura 1. La heredabilidad directa y materna para peso al nacimiento fue de 0.23±0.086 y 0.05±0.05, respectivamente y para PA270 fue 0.39±0.127 y 0.39±0.091, respectivamente. La correlación genética (0.11) y fenotípica (0.13) entre las características de PN y PA270 en la población de Brahman comercial evaluada, fueron positivas y muy bajas.

El progreso genético por año para PN y PA270 fue de -2.02 kg y 7.81 kg, respectivamente. En la tabla 3 se pueden observar los 10 mejores animales y las 10 mejores madres para PN y PA270 con respecto a la evaluación genética realizada por el presente estudio para Brahman comercial. El mejor animal para peso al nacimiento fue la hembra “1022” con una DEP de -3.12 kg y para peso al destete la hembra “48” con una DEP de 38.34 kg. La mejor madre para peso al nacimiento fue la “1020” con una DEP materna de -1.48 kg y para peso al destete la mejor madre fue el animal “326” con una DEP materna de 23.35 kg.

<bold>Tabla 2</bold>. Componentes de varianza, parámetros genéticos y fenotípicos para las características del peso al nacimiento y peso al destete en ganado Brahman comercial. Tabla 2. Tabla 2. Componentes de varianza, parámetros genéticos y fenotípicos para las características del peso al nacimiento y peso al destete en ganado Brahman comercial.
Características σ²a σ²m σ²e σ²p h²d±EE h²m ±EE ΔG/año
PN 4.04 0.81 5.4 9.21 0.44 ±0.14 0.09 ±0.092 -2.02
PA270 348.9 376.72 224.27 657.79 0.39 ±0.12 0.42 ±0.087 7.81
σ²a: Varianza aditiva. σ²m: Varianza materna. σ²e: Varianza ambiental. σ²p: Varianza fenotípica. h²d ± EE: Heredabilidad directa ± error estándar. h²m ± EE: Heredabilidad materna ± error estándar. ΔG/año: Progreso genético por año.

<bold>Figura 1.</bold> Comportamiento de los valores genéticos para peso al nacimiento y peso al destete por año en una población de ganado Brahman comercial no seleccionada. Figura 1 Figura 1. Comportamiento de los valores genéticos para peso al nacimiento y peso al destete por año en una población de ganado Brahman comercial no seleccionada.

<bold>Tabla 3</bold>. Ranking de animales con los mejores valores genéticos directos y maternos para las características de peso al nacimiento (PN) y peso al destete (PD) en una población de Brahman comercial. Tabla 3 Tabla 3. Ranking de animales con los mejores valores genéticos directos y maternos para las características de peso al nacimiento (PN) y peso al destete (PD) en una población de Brahman comercial.
Peso al Nacimiento Peso al Destete
ID Sexo DEPdir ID Madres DEPmat ID Sexo DEPdir ID Madres DEPmat
1022 H -3.12 1020 -1.48 48 H 38.34 326 23.35
399 H -2.68 399 -1.27 1184 H 32.19 554 18.12
14 M -2.68 91 -1.27 1194 M 31.04 510 17.77
91 H -2.59 931 -1.23 1401 H 29.08 168 16.45
933 H -2.35 1332 -1.11 69 H 28.62 356 15.91
1335 H -2.14 407 -1.02 1424 H 27.51 655 15.75
342 M -2.14 1417 -1.01 1055 M 27.31 129 15.71
407 H -2.06 184 -0.98 73 H 27.30 213 15.70
883 M -1.95 1054 -0.92 1144 H 26.18 365 15.58
1420 H -1.94 375 -0.92 1091 H 25.40 368 15.19
ID: Identificación del animal. H: Hembra. M: Macho. DEPdir: Diferencia esperada de progenie directa. DEPmat: Diferencia esperada de progenie materna.

<bold>DISCUSIÓN</bold>

En Colombia con ganado Brahman, otros autores han encontrado pesos al nacimiento bajo condiciones de pastoreo promedios de 35±3 kg y 34.88±3.9 kg, respectivamente (3,14). Esto demuestra que hay diferencias principalmente en el manejo de los hatos comerciales en el trópico colombiano. Es posible que el peso promedio de la población evaluada en el presente estudio haya sido menor, dado que las condiciones de manejo y alimentación en cada finca comercial son diferentes y este centro de producción específicamente tiene un manejo extensivo, sin suplementación y cargas de animales muy altas.

Los pesos al nacimiento obtenidos en este trabajo son similares a los obtenidos por otros autores en Brasil y México, donde evaluaron ganado Nelore y Brahman comercial y encontraron medias de 32.2±1.8 kg y 32.30±3.80 kg, respectivamente (15, 16). En Colombia también evaluaron poblaciones de ganado criollo comercial y encontraron promedios de 28.92±3.2 Kg en el departamento de Arauca (17) y en ganado criollo Costeño con cuernos de 29±3.71 kg (18). Esto resultados están dentro de los rangos encontrados en este estudio ya que para la característica peso al nacimiento existe una menor influencia del ambiente agroecológico, nutricional y de manejo de las diferentes fincas.

La media hallada para PA270 fue de 216.71±26.97 kg, fue un promedio normal para cebuínos comerciales pesados al destete y es superada por las reportadas por autores colombianos cuyos valores obtenidos fueron de 241±30 kg y 237.08±35.60 kg, respectivamente (3,14). Valores similares de 210 kg para PA270 fueron encontrados para la raza Brahman (19); demostrando así, una variabilidad en el peso de la misma raza por selección y tipos de manejo en cada hato, pudiendo alcanzar pesos entre los 210 y 250 kg bajo condiciones tropicales colombianas. También existen trabajos que reportan pesos al destete en ganado criollo colombiano en el trópico que oscilan entre 169±24.7 kg y 137.65±28.04 kg, respectivamente (18,20). En Brasil, se obtuvieron pesos para la raza Nelore de 185±29.40 kg (15) y en Brahman de México promedios de 180.7±32.6 kg (16); estos resultados son diferentes a los del presente estudio, seguramente por las condiciones zootécnicas y ambientales diferentes bajo las cuales fueron evaluadas estas poblaciones de Brasil y México. Al comparar estos resultados se puede inferir que hay gran variabilidad genética y fenotípica entre poblaciones de ganado cebú en diferentes países y que existe un rango de peso que se puede trabajar por selección para mejorar los pesos al destete de la población evaluada en el presente estudio.

La heredabilidad directa y materna, encontrada en este estudio para peso al nacimiento fue baja (0.23) y casi inexistente (0.05), respectivamente. En estudios similares se reportan valores bajos en Brahman comercial de 0.33 y 0.08; 0.32 y 0.16, respectivamente (14, 21) y en ganado Romosinuano de 0.25 y 0.06, respectivamente (20). Otros autores en ganado Romosinuano criado en condiciones similares a las del presente estudio, encontraron valores por debajo de 0.17 y 0.12 (4,22). En Brasil en raza Nelore reportan valores de 0.25 y 0.25, respectivamente (15,23). Las heredabilidades maternas para peso al nacimiento en la mayoría los trabajos fueron muy bajas, lo que significa que la influencia de la genética de la madre es muy poca durante la gestación y la genética aditiva directa y el ambiente uterino tienen un efecto mayor sobre esta característica. La heredabilidad directa fue baja por lo que se debe poner especial atención en mejorar la alimentación y el manejo en el centro de producción evaluado; además, se debe aumentar la presión de selección en las hembras para de esta forma aumentar el desempeño productivo de los terneros y de las hembras madres de las próximas generaciones.

Para PA270 se obtuvo una h²d de 0.39±0.13 y h²m de 0.39±0.091, estos valores son medios. Los valores encontrados de h²d para PA270 en el presente trabajo están por encima de los reportados en Brahman, Nelore y Romosinuano de 0.22, 0.28 y 0.34, respectivamente (4,21,23). De esto se puede inferir que la raza Brahman puede obtener individuos con mejor desempeño en su progenie si se realiza un buen proceso de selección genética y apareamiento direccionado. Otros autores en Brahman comercial encentraron valores similares (0.37) de h²d (14) y superiores (0.45 y 0.63) en ganado Brangus y Blanco Orejinegro (7,24). Los resultados obtenidos son importantes ya que una heredabilidad media tanto directa como materna, provee la oportunidad para realizar selección y mejoramiento de las características productivas tanto en los terneros como en las madres y de esta forma aumentar la eficiencia productiva y reproductiva de los hatos bovinos comerciales. La raza Brahman posee gran potencial por su adaptación al trópico, pero los centros de producción de bovinos de carne comercial no realizan evaluación genética de toros y vacas y por lo tanto no tienen criterio alguno para la selección de los futuros reproductores. Los sistemas productivos que vienen haciendo selección sobre animales con gran potencial genético para la producción de carne tienden a tener mayor productividad de kilogramos/hectárea.

Las correlaciones genéticas (ϒg) y fenotípicas (ϒp) para PN y PA270 fueron bajos posiblemente porque en la población de ganado evaluada no se tenían criterios de selección; además, la característica peso al nacimiento tiende a ser menos influenciada por el ambiente agroecológico, nutricional y de manejo y el peso al destete tiene una mayor influencia del ambiente materno y nutricional; esto hace que las dos características así sean de peso no se correlacionen en gran proporción. Sin embargo, en Colombia otros autores reportan en la misma raza valores de correlación genética de 0.64, 0.51 y -0.2 (14,19,21) y en Brasil en la raza Nelore valores de 0.36 y 0.27 (15,23). La diferencia entre los valores reportados con los encontrados, se debe al trabajo de selección de animales mejorantes para las dos características que se hace en otras ganaderías especializadas. Las correlaciones fenotípicas encontradas son similares en razas criollas y a las reportadas en otros trabajos realizados en Brahman comercial, con valores de 0.13, 0.14 y 0.16, respectivamente (14,18,22).

El progreso genético por año para PN y PA270 fue de -2.02 kg y 7.81 kg, respectivamente. En la figura 1 se puede observar la ganancia genética del PA270 durante los años 2002 a 2014. En esta figura se nota un aumento en el progreso genético entre los años 2010 a 2014, debido principalmente a la inclusión de programas de inseminación artificial a tiempo fijo, lo que obliga al productor a usar toros evaluados, probados y que poseen gran potencial genético para las características de peso.

En conclusión, las estimaciones de heredabilidad indican la existencia de variabilidad genética aditiva directa, evidenciando la oportunidad de ganancia genética por medio de selección. Sin embargo, las características de PN y PA270 presentan también gran influencia del ambiente, por lo cual se debe poner especial atención al manejo zootécnico y la nutrición de las hembras durante la gestación y de los terneros del nacimiento hasta el destete.

La inseminación artificial es un óptimo medio de introducción de material genéticamente superior en los hatos ganaderos, contribuyendo así con la ganancia genética en las futuras generaciones.

Es importante usar los resultados de la evaluación genética y de la valoración de los animales como estrategia para los apareamientos futuros que permitan un mejoramiento genético de las características de peso en la población de Brahman comercial.

La implementación de programas de mejora genética y de evaluaciones genéticas en ganaderías comerciales puede hacer más eficiente la selección de animales con gran potencial genético, los cuales deben ser incluidos en programas de apareamientos dirigidos; mejorando así, los parámetros productivos de generación en generación y haciendo más rentable la ganadería de carne en Colombia.

Conflicto de intereses

Los Autores declaramos expresamente a la revista MVZ de la Universidad de Córdoba Colombia, que el manuscrito remitido no ha sido publicado previamente; y que no tenemos conflicto de intereses de ninguna índole.

 Agradecimientos

Agradecemos a la hacienda La Envidia por la disponibilidad de los registros y al grupo de investigación en agronomía y Zootecnia-GIAZ de la Universidad Católica de Oriente, por el acompañamiento durante la investigación.

 <bold>REFERENCIAS</bold> 1. Organización de las Naciones Unidas para la Alimentación y la Agricultura-FAO. Perspectivas alimentarias (Resumen de mercado). FAO [en línea] 2018 [fecha de acceso 26 de octubre de 2018]; 12. URL disponible en: http://www.fao.org/giews/reports/food-outlook/es/ Perspectivas alimentarias (Resumen de mercado). Organización de las Naciones Unidas para la Alimentación y la Agricultura-FAO 2018 http://www.fao.org/giews/reports/food-outlook/es/ 2. DANE. La estructura de la producción de carne bovina en Colombia. Sipsa [en línea]. 2012 [fecha de acceso 26 de octubre de 2018]; 4:75. URL disponible en: https://www.dane.gov.co/files/investigaciones/agropecuario/sipsa/insumos_factores_de_produccion_octubre_2012.pdf La estructura de la producción de carne bovina en Colombia. Sipsa DANE. 2012 https://www.dane.gov.co/files/investigaciones/agropecuario/sipsa/insumos_factores_de_produccion_octubre_2012.pdf 3. Montes D, Vergara O, Prieto E, Barragán W. Estimación de la repetibilidad y factores que afectan el peso al nacer y al destete en ganado bovino cebú brahmán. Rev Colombiana Cienc Anim. 2009; 1(1):19-36. https://doi.org/10.24188/recia.v1.n1.2009.405 Montes D Vergara O Prieto E Barragán W Estimación de la repetibilidad y factores que afectan el peso al nacer y al destete en ganado bovino cebú brahmán. Rev Colombiana Cienc Anim. 2009 1 1 19 36 10.24188/recia.v1.n1.2009.405 4. Martínez R y Pérez E. Parámetros y tendencias genéticas para características de crecimiento en el ganado criollo colombiano Romosinuano. Corpoica Cienc Tecnol Agropecu. 2006; 7(1):25-32. https://doi.org/10.21930/rcta.vol7_num1_art:56 Martínez R Pérez E Parámetros y tendencias genéticas para características de crecimiento en el ganado criollo colombiano Romosinuano. Corpoica Cienc Tecnol Agropecu. 2006 7 1 25 32 10.21930/rcta.vol7_num1_art:56 5. Osorio M, Segura C. Factores que afectan el peso al nacer y al destete de becerros de doble propósito en el trópico. Livest Res Rural Dev. 2008; 20(1):Article15. URL Disponible en: http://www.lrrd.org/lrrd20/1/osor20015.htm Osorio M Segura C Factores que afectan el peso al nacer y al destete de becerros de doble propósito en el trópico. Livest Res Rural Dev. 2008 20 1 http://www.lrrd.org/lrrd20/1/osor20015.htm 6. Quintero J, Triana J, Quijano J, Arboleda E. Influencia de la inclusión del efecto materno en la estimación de parámetros genéticos del peso al destete en un hato de ganado de carne. Rev Col Cienc Pec. 2007; 20(2):117-123. https://aprendeenlinea.udea.edu.co/revistas/index.php/rccp/article/view/324128 Quintero J Triana J Quijano J Arboleda E Influencia de la inclusión del efecto materno en la estimación de parámetros genéticos del peso al destete en un hato de ganado de carne Rev Col Cienc Pec. 2007 20 2 117 123 https://aprendeenlinea.udea.edu.co/revistas/index.php/rccp/article/view/324128 7. Cañas J, Ramírez J, Arboleda O, Ochoa J, Vergara O, Cerón M. Estimación de parámetros genéticos para peso al destete en ganado Blanco Orejinegro (BON) en el noroccidente colombiano. Rev MVZ Córdoba. 2008; 13(1):1138-1145. https://doi.org/10.21897/rmvz.405 Cañas J Ramírez J Arboleda O Ochoa J Vergara O Cerón M Estimación de parámetros genéticos para peso al destete en ganado Blanco Orejinegro (BON) en el noroccidente colombiano. Rev MVZ Córdoba 2008 13 1 1138 1145 10.21897/rmvz.405 8. Castro D. Ensayo sobre tipología de suelos colombianos-Énfasis en génesis y aspectos ambientales. Rev Acad Colomb Cienc. 2003; 27(104):319-341. http://www.accefyn.com/revista/Vol_27/104/319-341.pdf Castro D Ensayo sobre tipología de suelos colombianos-Énfasis en génesis y aspectos ambientales. Rev Acad Colomb Cienc. 2003 27 104 319 341 http://www.accefyn.com/revista/Vol_27/104/319-341.pdf 9. Kinghorn BP. An algorithm for efficient constrained mate selection. Genet Sel Evol. 2011; 43(1):4. DOI: https://doi.org/10.1186/1297-9686-43-4 Kinghorn BP An algorithm for efficient constrained mate selection Genet Sel Evol. 2011 43 1 4 10.1186/1297-9686-43-4 10. R Development Core Team. R: a language and environment for statistical computing. [En linea] Version 3.5.1. Vienna: R Foundation for Statistical Computing. 2018. [fecha de acceso 26 de octubre de 2018]; URL disponible en: http://www.Rproject.org/ R: a language and environment for statistical computing Vienna: R Foundation for Statistical Computing 2018 http://www.Rproject.org/ 11. Falconer D, Mackay T. Introduction to Quantitative Genetics. Harlow, UK: Longman Group Ltd; 1996. Falconer D Mackay T Longman Group Ltd 1996 Harlow, UK Introduction to Quantitative Genetics. 12. Boldman K, Kriese L, Van Vleck L, Van Tessell C, Kachman S. A manual for use of MTDFREML. A set of programs to obtain estimates of variances and covariances (Draft). USA; USDA-Agric Res Serv: 1995. Boldman K Kriese L Van Vleck L Van Tessell C Kachman S A manual for use of MTDFREML. A set of programs to obtain estimates of variances and covariances (Draft) USDA-Agric Res Serv 1995 13. Daniel Gianola KH. Advances in Statistical Methods for Genetic Improvement of Livestock. vol 18. Lexington, KY: Springer; 1990. Daniel Gianola KH Advances in Statistical Methods for Genetic Improvement of Livestock Lexington, KY: Springer 1990 vol 18 14. Montes V, Vergara G, Prieto M, & Rodríguez P. Estimación de los parámetros genéticos para el peso al nacer y al destete en ganado bovino de la raza Brahman. Rev MVZ Córdoba. 2008; 13(1):1184-1191. https://doi.org/10.21897/rmvz.409 Montes V Vergara G Prieto M Rodríguez P Estimación de los parámetros genéticos para el peso al nacer y al destete en ganado bovino de la raza Brahman. Rev MVZ Córdoba 2008 13 1 1184 1191 10.21897/rmvz.409 15. Chud T, Caetano S, Buzanskas M, Grossi D, Guidolin, D, Nascimento G, Munari D. Genetic analysis for gestation length, birth weight, weaning weight, and accumulated productivity in Nellore beef cattle. Livest Sci. 2014; 170:16-21. https://doi.org/10.1016/j.livsci.2014.09.024 Chud T Caetano S Buzanskas M Grossi D Nascimento G Munari D Genetic analysis for gestation length, birth weight, weaning weight, and accumulated productivity in Nellore beef cattle. Livest Sci. 2014 170 16 21 10.1016/j.livsci.2014.09.024 16. Parra G, Martínez J, Cienfuegos E, García F, Ortega E. Parámetros genéticos de variables de crecimiento de ganado Brahman de registro en México. Veterinaria México. 2007; 38(2):217-229. https://doi.org/10.4067/s0301-732x2009000200004 Parra G Martínez J Cienfuegos E García F Ortega E Parámetros genéticos de variables de crecimiento de ganado Brahman de registro en México. Veterinaria México. 2007 38 2 217 229 10.4067/s0301-732x2009000200004 17. Salamanca C, Quintero V, Benitez M. Características de crecimiento predestete en becerros del sistema doble propósito en el municipio de Arauca. Zoot Trop. 2011; 29(4):455-465. Disponible en http://www.scielo.org.ve/scielo.php?script=sci_arttext&pid=S0798-72692011000400007&lng=es&nrm=iso. Salamanca C Quintero V Benitez M Características de crecimiento predestete en becerros del sistema doble propósito en el municipio de Arauca. Zoot Trop. 2011 29 4 455 465 http://www.scielo.org.ve/scielo.php?script=sci_arttext&pid=S0798-72692011000400007&lng=es&nrm=iso 18. García E, Mitat A, Saraz A, Simanca C, Garay D. Componentes de covarianza y parámetros genéticos para características de crecimiento en una población de ganado Costeño con Cuernos en Colombia. Livest Res Rural Dev. 2015; 27(1):Article6. URL Disponible en: http://www.lrrd.org/lrrd27/1/pere27006.html García E Mitat A Saraz A Simanca C Garay D Componentes de covarianza y parámetros genéticos para características de crecimiento en una población de ganado Costeño con Cuernos en Colombia. Livest Res Rural Dev. 2015 27 1 http://www.lrrd.org/lrrd27/1/pere27006.html 19. Plasse D, Verde O, Fossi H, Romero R, Hoogesteijn R, Bastidas P, Bastardo J. (Co) variance components, genetic parameters and annual trends for calf weights in a pedigree Brahman herd under selection for three decades. J Anim Breed Genet. 2002; 119(3):141-153. https://doi.org/10.1046/j.1439-0388.2002.00321.x Plasse D Verde O Fossi H Romero R Hoogesteijn R Bastidas P Bastardo J (Co) variance components, genetic parameters and annual trends for calf weights in a pedigree Brahman herd under selection for three decades. J Anim Breed Genet. 2002 119 3 141 153 10.1046/j.1439-0388.2002.00321.x 20. Ossa G, Pérez J, Suárez M. Valores genéticos de caracteres productivos y reproductivos en bovinos Romosinuano. Corpoica Cienc Tecnol Agropecu. 2008; 9(1):93-101. https://doi.org/10.21930/rcta.vol14_num2_art:503 Ossa G Pérez J Suárez M Valores genéticos de caracteres productivos y reproductivos en bovinos Romosinuano. Corpoica Cienc Tecnol Agropecu. 2008 9 1 93 101 10.21930/rcta.vol14_num2_art:503 21. Jiménez A, Manrique C, Martínez C. Parámetros y valores genéticos para características de composición corporal, área de ojo de lomo y grasa dorsal medidos mediante ultrasonido en la raza Brahman. Rev Med Vet Zoot 2010; 57(3):159-167. https://revistas.unal.edu.co/index.php/remevez/article/view/18236 Jiménez A Manrique C Martínez C Parámetros y valores genéticos para características de composición corporal, área de ojo de lomo y grasa dorsal medidos mediante ultrasonido en la raza Brahman. Rev Med Vet Zoot 2010 57 3 159 167 https://revistas.unal.edu.co/index.php/remevez/article/view/18236 22. Ossa G, Suárez M, Pérez E. Efecto del medio y la herencia sobre el peso al nacimiento, destete y a los 16 meses de edad en terneros de la raza Romosinuano. Corpoica Cienc Tecnol Agropecu. 2007; 8(2):81-92. https://doi.org/10.21930/rcta.vol8_num2_art:98 Ossa G Suárez M Pérez E Efecto del medio y la herencia sobre el peso al nacimiento, destete y a los 16 meses de edad en terneros de la raza Romosinuano. Corpoica Cienc Tecnol Agropecu. 2007 8 2 81 92 10.21930/rcta.vol8_num2_art:98 23. Boligon A, Albuquerque G, Mercadante E, Lôbo B. Herdabilidades e correlações entre pesos do nascimento à idade adulta emrebanhos da raça Nelore. R Bras Zootec. 2009; 38(12):2320-2326. https://doi.org/10.1590/s1516-35982009001200005 Boligon A Albuquerque G Mercadante E Lôbo B Herdabilidades e correlações entre pesos do nascimento à idade adulta emrebanhos da raça Nelore. R Bras Zootec. 2009 38 12 2320 2326 10.1590/s1516-35982009001200005 24. Vergara O, Martínez N, Almanza R, Patiño R, Calderón A. Parámetros y tendencias genéticas para característica de crecimiento predestete en una población bovina multirracial en Colombia. Rev Facult Cienc Vet, UCV. 2014; 55(2):68-77. Disponible en http://www.scielo.org.ve/scielo.php?script=sci_arttext&pid=S0258-65762014000200002&lng=es&nrm=iso. Vergara O Martínez N Almanza R Patiño R Calderón A Parámetros y tendencias genéticas para característica de crecimiento predestete en una población bovina multirracial en Colombia. Rev Facult Cienc Vet, UCV. 2014 55 2 68 77 http://www.scielo.org.ve/scielo.php?script=sci_arttext&pid=S0258-65762014000200002&lng=es&nrm=iso.