Introducción
En los últimos años el Grupo Intergubernamental de Expertos sobre el Cambio Climático (IPCC)1, La Organización de las Naciones Unidas para la Alimentación y Agricultura (FAO)2, han señalado que el sector agropecuario es la principal fuente de gases de efecto invernadero (GEI) de vida corta, metano (CH4) y óxido nitroso (N2O).
Algunas estrategias para la mitigación de metano proveniente del ganado lechero incluyen la reducción del hato bovino, el cambio de alimentación, el uso de suplementos, la inmunización contra las arqueas metanogénicas y la selección de animales con menor producción de CH43.
Para crear estrategias de selección de animales con menor producción de metano es necesario el conocimiento de las correlaciones genéticas entre la producción de CH4 y otras características de importancia productiva y económica4.
El índice de selección es una metodología para maximizar la mejora genética en un objetivo específico5. Los índices de selección se han aplicado extensivamente en la estimación del valor de reproducción del ganado lechero para rasgos individuales, así como para combinaciones de rasgos para propósitos de selección6.
En bovinos y ovinos, se han demostrado la existencia variación de las emisiones de CH4 entre los individuos alimentados con la misma dieta7. De Haas et al8 mencionan que existe la posibilidad de seleccionar vacas con baja emisión de CH4, ya que la variación genética sugiere que las reducciones serian de 11 a 26 % en 10 años, y podrían ser aún mayores en un programa de selección genómica. Sin embargo, en la actualidad, la información disponible sobre las oportunidades para mitigar el CH4 entérico a través del mejoramiento genético es escasa, pero se podría pensar que la selección genética en contra de las emisiones de metano podría tener repercusiones en características productivas de interés económico. El objetivo de este trabajo fue estimar la respuesta a la selección a través de diferentes índices de selección entre producción de metano y la producción de leche y sus componentes en tres sistemas de producción de leche en México.
Material y métodos
El presente estudio se realizó en tres unidades de producción (UP) del sistema de doble propósito (DP), dos de lechería tropical especializada (LTE) y cuatro de sistema familiar (SF) (Cuadro 1). Se midió los componentes de la leche y emisiones de metano a un total de 274 vacas (98, 74 y 102 en los sistemas de DP, LTE y SF respectivamente).
Rancho | Sistema | n | Localización | Razas |
---|---|---|---|---|
La Posta | DP | 33 | Veracruz | HOC y SPC |
El Zapato | DP | 16 | Veracruz | HOC |
La Doña | DP | 49 | Puebla | HOC, SPC y SMC |
Santa Elena | LTE | 37 | Puebla | HO, SP y HOSP |
Aguacatal | LTE | 37 | Puebla | HO, SP y HOSP |
Rancho 5 | SF | 16 | Jalisco | HO |
Rancho 6 | SF | 32 | Jalisco | HO |
Rancho 7 | SF | 24 | Jalisco | HO |
Rancho 8 | SF | 30 | Jalisco | HO |
DP= doble propósito; LTE= lechería tropical especializada; SF= familiar.
HOC= Holstein x Cebú, SPC= Suizo Pardo x Cebú, SMC= Simmental x Cebú, HO= Holstein, SP= Suizo Pardo, HOSP= Holstein x Suizo Pardo.
Dos de las tres UP de DP se encuentran ubicadas en el municipio Medellín de Bravo, Veracruz, con un clima tropical subhúmedo cálido Aw(o) y una altitud de 12 msnm9. La temperatura media anual es de 25 °C, así como la precipitación pluvial media anual es de 1,460 mm9. La tercera UP de DP, al igual que las UP de LTE, se encuentra en el municipio de Hueytamalco, Puebla; con una altitud de 240 msnm. El clima es subtropical húmedo Af(c) con temperatura media anual de 23 ºC y precipitación pluvial media que varía de 2,200 a 2,500 mm9.
Las cuatro UP del SF se encuentran en el municipio de Tepatitlán, Jalisco con un clima templado subhúmedo cálido (A)C(w1) (e)g a una altitud de 1,927 msnm., con temperatura media anual de 18 ºC y precipitación pluvial media anual de 715 mm9.
En el sistema de DP se utilizan principalmente ganado cruzado Bos taurus taurus y Bos taurus indicus. En el caso de las razas Bos taurus indicus las más utilizadas son Brahman, Gyr y Sardo Negro; mientras que para las razas Bos taurus taurus se encuentran principalmente Holstein, Suizo Pardo y Simmental10.
Una de las variantes de los sistemas tropicales de producción de leche bovina es la LTE. Este sistema se caracteriza por la utilización de razas puras tales como Holstein y Suizo Pardo. El manejo general de la LTE es muy similar al manejo del sistema de DP; salvo en la crianza de becerros, que se realiza de manera artificial, y en la ordeña que se realiza sin el apoyo del becerro10.
Los SF se caracterizan por pequeñas unidades de producción que fluctúan entre 3 y 30 vacas. Las unidades de producción están condicionadas a pequeñas superficies de terreno, aledañas a las viviendas, por lo que se le llama también de "traspatio". Pueden ser de tipo estabulado o semiestabulado, de acuerdo con las condiciones del campo de cultivo. Los animales son en su mayoría de la raza Holstein. El nivel tecnológico se puede considerar escaso; ya que los productores no realizan prácticas adecuadas de alimentación, manejo reproductivo, medicina preventiva o mejoramiento genético. En este sistema se carece de registros de producción y las instalaciones son rudimentarias, principalmente se realiza el ordeño manual. La alimentación es basada en el pastoreo o en el suministro de forrajes y esquilmos provenientes de los cultivos del productor11.
Muestreo de metano
Para el muestreo de CH4 se utilizó la metodología desarrollada por Garnsworthy et al12, mediante un análisis continuo de concentraciones de gas, a través de un sistema de doble longitud de onda no disperso, el cual recupera la información segundo a segundo de la concentración de gases en el ambiente, esto por medio de los Guardian NG- Infrared Gas Monitor (Edinburgh Instruments, Escocia, Reino Unido).
Los equipos se instalaron en los comederos donde se les ofrecía alimento a las vacas durante la ordeña. Se realizaron adaptaciones a los diferentes tipos de comederos, de tal manera que se creara una atmosfera cerrada, para evitar que las corrientes de aire sesgaran las concentraciones de CH4. Lo anterior con el objetivo de generar el menor disturbio en la rutina de ordeño. Lo que permitiría tomar una muestra de la atmósfera del comedero mientras el animal se alimenta normalmente.
Se realizó un periodo de adaptación de una semana a la presencia de los nuevos comederos. El periodo de medición consistió en dos semanas de muestreo de CH4 durante la ordeña; con el objetivo de contar con un mínimo de 10 días efectivos de medición en cada UP.
Muestreo de leche
Se obtuvo una muestra de leche de manera individual durante el ordeño; cada muestra consistió en un mínimo de 50 ml de leche obtenida directamente de los pesadores al inicio de las mediciones. Las muestras se conservaron con bronopol y se identificaron mediante el número de la UP y el número de identificación del animal. El análisis de la leche se realizó en el laboratorio de calidad de leche de la Asociación Holstein de México A.C. mediante la técnica de infrarrojo medio para medir el porcentaje de grasa y porcentaje de proteína.
Muestreo y extracción de ADN
Se obtuvieron muestras de pelos con folículos pilosos directamente de la cola de todos los animales involucrados en el experimento. Las muestras de pelo fueron identificadas y enviadas al laboratorio GENESEEK (Lincoln, Nebraska). En este laboratorio se extrajo el ADN y se obtuvieron los genotipos mediante microarreglos de alta densidad. Para los animales provenientes del SF, se utilizó el arreglo GGP BOVINE LD V4, mismo que permite obtener 30,125 SNP. En el caso de los animales cruzados del LTE y DP se optó por el arreglo de GGDP BOVINE 150K que hace posible la identificación de 138,962 SNP por animal; lo anterior obedeció a que los animales cruzados requieren un mayor número de marcadores para que la información sea útil. En este estudio solo se incluyeron los SNP ubicados en los 29 cromosomas autosómicos. El control de calidad de los genotipos se realizó a través del software PLINK 1.713, y consistió en: 1) remover los individuos con menos del 90 % de la información genotípica, 2) remover los animales con un valor menor a 5 % en la frecuencia del alelo menor, y 3) remover los animales con una proporción menor de 90 % de marcadores útiles. Al final del control de calidad y manteniendo los marcadores en común en ambas plataformas el número de marcadores disponibles fue de 20,776 SNP para cada animal.
Análisis estadístico
Estimación de relaciones genéticas genómicas
La estimación de componentes de varianza se realizó mediante la metodología de modelos mixtos. Debido a que no se contó con información genealógica completa para la construcción de la matriz de relaciones aditivas (Matriz A), se utilizaron los marcadores moleculares para construir la matriz de relaciones genómicas entre todos los animales (G). La matriz G se basó en el método propuesto por VanRaden14, que consiste en construir la matriz M a través de las dimensiones número de individuos (n) por el número de marcadores (m). Los elementos de la matriz se codificaron como -1 (homocigótico para un alelo), 0 (heterocigotos) y 1 para (homocigótico para el otro alelo). A la matriz M se le substrae la matriz P(nxm) que contiene columnas con todos los elementos 2(pi-0.5), donde pi es la frecuencia del segundo alelo en el locus i y se genera la matriz Z (Z = M - P). Finalmente, la matriz G se calculó como:
Estimación de componentes de varianza
La estimación de componentes de varianza para las emisiones de CH4 durante el ordeño y los componentes de la leche (porcentaje de grasa y proteína), se implementaron mediante el programa ASReml-R15.
La selección del modelo se realizó con los efectos de producción diaria de leche durante las mediciones, días en lactancia, periodo de lactancia, número de lactancia, sistema de producción, número de hato y raza dentro del sistema, para las variables respuesta producción de metano durante la ordeña, porcentaje de grasa, porcentaje de proteína. Se probaron todas las combinaciones lógicas dentro de los efectos fijos y aleatorios que lograran converger con las variables respuesta, de donde se obtuvo el presente modelo. El modelo univariado se representa de la siguiente manera:
Donde,
Estimación de componentes de covarianza
La estimación de componentes de covarianza entre las emisiones de CH4 durante el ordeño y los componentes de la leche, se calcularon mediante el programa ASReml-R15. Se utilizaron las varianzas estimadas mediante los modelos univariados como valores iniciales para la estimación de las covarianzas.
El modelo bivariado se representa en términos matriciales de la siguiente manera:
Donde,
Estimación de parámetros genéticos
Las h2 fueron obtenidas a partir de los componentes de varianza
estimados con los modelos univariados. Las correlaciones genéticas (
Las
Índice de selección
La estimación de la respuesta a la selección se llevó a cabo mediante diferentes índices de selección entre producción de metano y la producción de leche y sus componentes. Se realizó un análisis de sensibilidad, en el cual se identificaron diferentes escenarios en los que se podía seleccionar emisiones de CH4, porcentaje de grasa y porcentaje de proteína, y de esta manera observar la dinámica entre la exactitud de los índices y su ganancia genética (Cuadro 2). Los índices de selección se realizaron a cinco generaciones. A las características que se incluyeron en los índices se les asigno un valor dependiendo la importancia e intensidad de selección, la suma de los valores en cantidades absolutas debe ser igual a 100. A las emisiones de CH4 se le dieron valores que iban desde el 0 al -100 y al porcentaje de grasa y porcentaje de proteína valores que iban del 0 al 100.
Índice | CH4 | Grasa | Proteína | Índice | CH4 | Grasa | Proteína |
---|---|---|---|---|---|---|---|
INDEX1 | -100 | 0 | 0 | INDEX34 | 0 | 10 | 90 |
INDEX2 | -90 | 0 | 10 | INDEX35 | -40 | 40 | 20 |
INDEX3 | -90 | 10 | 0 | INDEX36 | -20 | 20 | 60 |
INDEX4 | -80 | 0 | 20 | INDEX37 | -30 | 30 | 40 |
INDEX5 | -80 | 10 | 10 | INDEX38 | -40 | 50 | 10 |
INDEX6 | -80 | 20 | 0 | INDEX39 | -10 | 20 | 70 |
INDEX7 | -70 | 0 | 30 | INDEX40 | 0 | 20 | 80 |
INDEX8 | -70 | 10 | 20 | INDEX41 | -40 | 60 | 0 |
INDEX9 | -60 | 0 | 40 | INDEX42 | -30 | 40 | 30 |
INDEX10 | -70 | 20 | 10 | INDEX43 | -20 | 30 | 50 |
INDEX11 | -70 | 30 | 0 | INDEX44 | -10 | 30 | 60 |
INDEX12 | -60 | 10 | 30 | INDEX45 | -30 | 50 | 20 |
INDEX13 | -50 | 0 | 50 | INDEX46 | -20 | 40 | 40 |
INDEX14 | -60 | 20 | 20 | INDEX47 | 0 | 30 | 70 |
INDEX15 | -50 | 10 | 40 | INDEX48 | -30 | 60 | 10 |
INDEX16 | -40 | 0 | 60 | INDEX49 | -10 | 40 | 50 |
INDEX17 | -60 | 30 | 10 | INDEX50 | -30 | 70 | 0 |
INDEX18 | -60 | 40 | 0 | INDEX51 | -20 | 50 | 30 |
INDEX19 | -50 | 20 | 30 | INDEX52 | -20 | 60 | 20 |
INDEX20 | -30 | 0 | 70 | INDEX53 | -10 | 50 | 40 |
INDEX21 | -40 | 10 | 50 | INDEX54 | -20 | 70 | 10 |
INDEX22 | -50 | 30 | 20 | INDEX55 | 0 | 50 | 50 |
INDEX23 | -20 | 0 | 80 | INDEX56 | 0 | 40 | 60 |
INDEX24 | -40 | 20 | 40 | INDEX57 | -20 | 80 | 0 |
INDEX25 | -10 | 0 | 90 | INDEX58 | -10 | 60 | 30 |
INDEX26 | -30 | 10 | 60 | INDEX59 | 0 | 60 | 40 |
INDEX27 | -50 | 40 | 10 | INDEX60 | -10 | 70 | 20 |
INDEX28 | 0 | 0 | 100 | INDEX61 | -10 | 80 | 10 |
INDEX29 | -50 | 50 | 0 | INDEX62 | 0 | 70 | 30 |
INDEX30 | -20 | 10 | 70 | INDEX63 | -10 | 90 | 0 |
INDEX31 | -40 | 30 | 30 | INDEX64 | 0 | 80 | 20 |
INDEX32 | -30 | 20 | 50 | INDEX65 | 0 | 90 | 10 |
INDEX33 | -10 | 10 | 80 | INDEX66 | 0 | 100 | 0 |
Los componentes de varianza y covarianza utilizados fueron los obtenidos mediante los modelos previamente descritos para los componentes de la leche (porcentaje de grasa y porcentaje de proteína) y producción de CH4 durante la ordeña.
La especificación original del índice de selección prevé el uso de una variable correlacionada (I) basada en el rendimiento fenotípico de cada animal para varias características5. Por lo tanto, se define como:
I = b p
Donde p es un vector de valores fenotípicos para el conjunto de criterios de selección y b son los factores de ponderación utilizados en la toma de decisiones de selección. Con el fin de maximizar la correlación de I con la contribución de cualquier candidato para la selección como posible padre, la información se combina como:
Ga = Pb
Donde G es una matriz nxm de variancias genéticas y covarianzas entre todos los rasgos m, a es un vector mx1 de valores relativos a la intensidad de selección para todos los rasgos. P es una matriz nxn de variancias fenotípicas y covarianzas entre los n rasgos medidos y disponibles como criterios de selección y b es un vector nx1 de factores de ponderación que se aplicará a los rasgos utilizados en la toma de decisiones de selección. La ecuación precedente se resuelve entonces como:
P -1 Ga = b
Para obtener los factores de ponderación contenidos en b y los candidatos para la selección se clasifican en base al índice (I).
La precisión del índice (
Donde, b es un vector de factores de ponderación que se aplicará a las características utilizadas en la toma de decisiones de selección; P es una matriz de variancias fenotípicas y covarianzas entre las características medidas y disponibles como criterios de selección; a es un vector de valores relativos para todas las características; Q es una matriz de varianzas genéticas y covarianzas entre todas las características consideradas parte del sistema.
Se estimó la ganancia genética
Donde, i= es la intensidad de la selección; G= es la matriz de
la varianza-covarianza genética de las características; b= es un vector de
factores de ponderación que se aplicará a las características utilizadas en la
toma de decisiones de selección.;
La desviación estándar del índice se calculó de la siguiente manera:
Donde, b= es un vector de factores de ponderación que se aplicará a las características utilizadas en la toma de decisiones de selección; P= es una matriz de variancias fenotípicas y covarianzas entre las características medidas y disponibles como criterios de selección.
Resultados
Las emisiones de CH4 en el sistema de LTE fueron de 0.08 mg de CH4/ L, y en SF y DP en ambos casos fue de 0.06 mg de CH4/ L (Cuadro 3). El promedio general para los tres sistemas fue de 0.065 mg de CH4/ L. En el caso de los componentes de la leche, el promedio del porcentaje de grasa en los tres sistemas fue de 4.82 %, y para cada uno de los sistemas fue 3.69 % (LTE), 3.72 % (SF) y 6.84 % (DP). En el caso del porcentaje de proteína en el sistema de LTE, SF y DP fue de 3.20 %, 3.29 % y 3.21 %, respectivamente. Siendo 3.23 % el promedio combinado para los tres sistemas.
Sistema | Metano (mg/L) | % de grasa | % de proteína | |||
---|---|---|---|---|---|---|
Media | DE | Media | DE | Media | DE | |
DP (n=98) | 0.06 | 0.039 | 6.84 | 4.938 | 3.21 | 0.405 |
SF (n=102) | 0.06 | 0.014 | 3.72 | 0.632 | 3.29 | 0.315 |
LTE (n=74) | 0.08 | 0.016 | 3.69 | 0.492 | 3.20 | 0.417 |
Promedio | 0.065 | 0.028 | 4.828 | 3.344 | 3.234 | 0.380 |
DP= sistema de doble propósito, SF= lechería Familiar, LTE= lechería tropical especializada, n= número de observaciones, DE= desviación estándar.
La h2 estimada para emisiones de CH4 durante la ordeña con el modelo univariado fue de 0.19, de la misma manera porcentaje de grasa presento una h2 de 0.39 y porcentaje de proteína de 0.18.
En contraste, las h2 estimadas mediante los modelos bivariados para las emisiones de CH4 durante la ordeña fue de 0.32±0.245 y el porcentaje de grasa fue de 0.46±0.278, respectivamente. En cambio, la h2 para porcentaje de proteína no mostro cambios con respecto al modelo univariado, sin embargo, la h2 de CH4 es similar a la encontrada en el análisis bivariado con porcentaje de grasa (0.35). Las correlaciones genéticas entre el porcentaje de grasa y proteína con las emisiones de CH4 durante la ordeña fueron de -0.090 ± 0.080 y -0.18 ± 0.575, respectivamente.
La exactitud de los índices de selección (
Índice | rHI | CH4 (mg/L) | Grasa
(%) |
Proteína (%) |
Índice | rHI | CH4 (mg/L) | Grasa
(%) |
Proteína (%) |
---|---|---|---|---|---|---|---|---|---|
INDEX1 | 19.58 | -0.021 | 0.030 | 0.073 | INDEX34 | 10.72 | -0.013 | 0.015 | 0.117 |
INDEX2 | 18.38 | -0.021 | 0.029 | 0.078 | INDEX35 | 10.63 | -0.021 | 0.030 | 0.086 |
INDEX3 | 17.92 | -0.021 | 0.030 | 0.073 | INDEX36 | 10.41 | -0.018 | 0.023 | 0.110 |
INDEX4 | 17.23 | -0.021 | 0.029 | 0.082 | INDEX37 | 10.32 | -0.020 | 0.027 | 0.100 |
INDEX5 | 16.72 | -0.021 | 0.030 | 0.078 | INDEX38 | 10.12 | -0.021 | 0.031 | 0.078 |
INDEX6 | 16.26 | -0.021 | 0.030 | 0.072 | INDEX39 | 9.95 | -0.016 | 0.020 | 0.115 |
INDEX7 | 16.14 | -0.021 | 0.028 | 0.088 | INDEX40 | 9.70 | -0.014 | 0.016 | 0.117 |
INDEX8 | 15.57 | -0.021 | 0.029 | 0.083 | INDEX41 | 9.70 | -0.021 | 0.032 | 0.069 |
INDEX9 | 15.13 | -0.020 | 0.027 | 0.093 | INDEX42 | 9.62 | -0.020 | 0.029 | 0.094 |
INDEX10 | 15.06 | -0.021 | 0.030 | 0.078 | INDEX43 | 9.56 | -0.018 | 0.025 | 0.107 |
INDEX11 | 14.61 | -0.021 | 0.031 | 0.072 | INDEX44 | 9.01 | -0.017 | 0.022 | 0.113 |
INDEX12 | 14.49 | -0.020 | 0.028 | 0.089 | INDEX45 | 9.01 | -0.020 | 0.030 | 0.087 |
INDEX13 | 14.21 | -0.020 | 0.026 | 0.099 | INDEX46 | 8.77 | -0.019 | 0.027 | 0.103 |
INDEX14 | 13.92 | -0.021 | 0.029 | 0.084 | INDEX47 | 8.70 | -0.014 | 0.018 | 0.117 |
INDEX15 | 13.50 | -0.020 | 0.027 | 0.095 | INDEX48 | 8.49 | -0.021 | 0.032 | 0.078 |
INDEX16 | 13.41 | -0.019 | 0.025 | 0.104 | INDEX49 | 8.12 | -0.017 | 0.024 | 0.111 |
INDEX17 | 13.40 | -0.021 | 0.030 | 0.078 | INDEX50 | 8.09 | -0.020 | 0.033 | 0.067 |
INDEX18 | 12.96 | -0.021 | 0.031 | 0.071 | INDEX51 | 8.05 | -0.020 | 0.029 | 0.097 |
INDEX19 | 12.86 | -0.020 | 0.028 | 0.090 | INDEX52 | 7.41 | -0.020 | 0.031 | 0.089 |
INDEX20 | 12.74 | -0.018 | 0.023 | 0.110 | INDEX53 | 7.28 | -0.018 | 0.026 | 0.107 |
INDEX21 | 12.63 | -0.019 | 0.026 | 0.101 | INDEX54 | 6.90 | -0.020 | 0.033 | 0.078 |
INDEX22 | 12.27 | -0.021 | 0.030 | 0.085 | INDEX55 | 6.78 | -0.016 | 0.022 | 0.114 |
INDEX23 | 12.23 | -0.017 | 0.020 | 0.114 | INDEX56 | 6.78 | -0.015 | 0.020 | 0.116 |
INDEX24 | 11.90 | -0.020 | 0.027 | 0.097 | INDEX57 | 6.52 | -0.020 | 0.034 | 0.064 |
INDEX25 | 11.89 | -0.015 | 0.017 | 0.117 | INDEX58 | 6.52 | -0.019 | 0.029 | 0.100 |
INDEX26 | 11.89 | -0.018 | 0.024 | 0.107 | INDEX59 | 5.89 | -0.017 | 0.025 | 0.110 |
INDEX27 | 11.76 | -0.021 | 0.031 | 0.078 | INDEX60 | 5.86 | -0.020 | 0.032 | 0.091 |
INDEX28 | 11.75 | -0.013 | 0.014 | 0.118 | INDEX61 | 5.35 | -0.020 | 0.034 | 0.077 |
INDEX29 | 11.32 | -0.021 | 0.032 | 0.070 | INDEX62 | 5.08 | -0.017 | 0.028 | 0.104 |
INDEX30 | 11.30 | -0.017 | 0.021 | 0.112 | INDEX63 | 5.02 | -0.019 | 0.036 | 0.057 |
INDEX31 | 11.23 | -0.020 | 0.029 | 0.092 | INDEX64 | 4.40 | -0.018 | 0.032 | 0.092 |
INDEX32 | 11.08 | -0.019 | 0.025 | 0.104 | INDEX65 | 3.90 | -0.018 | 0.035 | 0.073 |
INDEX33 | 10.91 | -0.015 | 0.018 | 0.116 | INDEX66 | 3.67 | -0.017 | 0.037 | 0.044 |
En los índices con mayor
Discusión
Las emisiones de CH4 durante la ordeña estimadas en el presente estudio son menores a los publicados por Bell et al17., posiblemente debido a la hetereogeneidad de las dietas en el presente estudio. Las explotaciones medidas en este trabajo tiene los animales en pastoreo y en los sistemas especializados tienen dietas a base de concentrados. Las h2 aquí estimadas para la producción de CH4 durante la ordeña son similares a los publicados por otros autores, los cuales midieron en cámaras metabólicas18,19 o incluso con las obtenidas mediante ecuaciones de predicción8. Estos resultados concuerdan con los estimadores obtenidos mediante metodologías similares para la medición de CH4 publicado por otros autores20. Lo anterior, sugiere que la medición de CH4 durante la ordeña es una característica que puede ser utilizada en programas de mejoramiento genético, debido a que puede ser incorporada fácilmente a programas de control de producción, y su medición utiliza equipos poco sofisticados en comparación con las cámaras respiratorias; además de que puede ser movilizado a lugares de difícil acceso. Las h2 estimadas para el porcentaje de grasa en este trabajo son mayores a la observada en otros estudios21-22. La h2 estimada para el porcentaje de proteína es menor a la publicada por Othmane et al23 con una h2 de 0.23, lo cual puede ser reflejo de la heterogeneidad propia de las unidades de producción de leche en función de la alimentación y la composición racial de los hatos.
Las correlaciones genéticas obtenidas entre emisiones de CH4 y los componentes de la leche sugieren un antagonismo genético; sin embargo, los estimadores de este trabajo no fueron diferentes de cero. En comparación a los presentados por Pszczola et al4) cuyo valor para la correlación entre CH4 y gramos de grasa en leche fue de mayor magnitud (0.21), y Lassen et al24 para los resultados de correlación de CH4 y gramos de proteína en leche fue de la misma magnitud (0.39); en ambos casos el signo contrario; cabe mencionar que esta diferencia se debe a las unidades de medida, ya que la producción de leche y los porcentajes de los componentes de la leche tienen una correlación negativa, mientras que la producción de leche con el contenido de los componentes de la leche tiene una correlación positiva25.
En la actualidad, se ha propuesto la factibilidad de reducir las emisiones de CH4 a través de la selección genética, misma que puede disminuir en 10 años las emisiones de CH4 entre un 11 y 26 % en el ganado lechero8 y un 5 % para el ganado de carne26.
Los índices de selección elaborados por Kandel et al27 incluyen producción de grasa y proteína, los cuales obtuvieron correlaciones positivas con la producción de CH4, si se toma en cuenta las unidades de medida, el resultado es similar al obtenido en este trabajo, ya que como se mencionó la producción de leche y los porcentajes de los componentes de la leche tienen una correlación negativa, mientras que la producción de leche con el contenido de los componentes de la leche tiene una correlación positiva27.
Conclusiones e implicaciones
Las correlaciones genéticas aquí estimadas entre las emisiones de CH4 con los componentes de la leche (porcentaje de grasa y porcentaje de proteína) sugieren la factibilidad de proponer un programa de mejoramiento genético para disminuir las emisiones de CH4, y simultáneamente aumentar los porcentajes de los componentes de la leche. En otras palabras, los resultados aquí presentados hacen posible el seleccionar genéticamente para disminuir las emisiones de CH4 y sin afectar de manera negativa la composición de la leche. Lo anterior se confirma mediante las ganancias genéticas predichas por generación a partir de los índices de selección. Lo anterior, con reducciones de entre 0.013 a 0.021 mg/L de emisiones de CH4 durante la ordeña en cinco generaciones, sin decremento de los porcentajes de grasa y proteína en la leche.