Unraveling the genetic basis of methane emission in dairy cattle: a comprehensive exploration and breeding approach to lower methane emissions

Ruminant animals, such as dairy cattle, produce CH4, which contributes to global warming emissions and reduces dietary energy for the cows. While the carbon foot print of milk production varies based on production systems, milk yield and farm management practices, enteric fermentation, and manure management are major contributors togreenhouse gas emissions from dairy cattle. Recent emerging evidence has revealed the existence of genetic variation for CH4 emission traits among dairy cattle, suggests their potential inclusion in breeding goals and genetic selection programs. Advancements in high-throughput sequencing technologies and analytical techniques have enabled the identification of potential metabolic biomarkers, candidate genes, and SNPs linked to methane emissions. Indeed, this review critically examines our current understanding of carbon foot print in milk production, major emission sources, rumen microbial community and enteric fermentation, and the genetic architecture of methane emission traits in dairy cattle. It also emphasizes important implications for breeding strategies aimed at halting methane emissions through selective breeding, microbiome driven breeding, breeding for feed efficiency, and breeding by gene editing.

Genetic parameter estimates for methane emission from breath during lactation and potential inaccuracies in reliabilities assuming a repeatability versus random regression model

Methane emissions will be added to many national ruminant breeding programs in the coming years. Little is known about the covariance structure of CH4 traits over a lactation, which is important for optimizing recording strategies and establishing optimal genetic evaluation models. Our aim was to study CH4 over a lactation using random regression (RR) models, and to compare the accuracy to a fixed regression repeatability model under different phenotyping strategies. Data were available from repeated measurements of CH4 concentrations (ppm) recorded in the feed bins of milking robots on 52 commercial dairy farms in the Netherlands. In total, 36,370 averaged weekly records were available from 4,664 cows. Genetic parameters were estimated using a fixed regression model, and a RR model with first- to fifth-order Legendre polynomials for the additive genetic and within-lactation permanent environmental effect. The mean heritability (± SE) was 0.17 ± 0.04, and the mean within-lactation repeatability was 0.56 ± 0.03. The genetic correlations between DIM were high and ranged from 0.34 ± 0.36 to 1.00 ± <0.01. Permanent environmental correlations showed large deviations and ranged from -0.73 ± 0.08 to 1.00 ± <0.01. With a large number of full lactation daughter CH4 records per bull, the reliability was not sensitive to using the fixed versus the RR model. However, when shorter periods were recorded at the start and end of the lactation, the fixed regression model resulted in a loss of reliability up to 28% for bulls. Assuming the fixed model when the true (co)variance structure is reflected by the RR model, more than twice as long of a recording from the start of lactation was required to achieve maximum reliability for a bull. Thus, a too simplistic model could result in implementing too little recording, and in lower genetic gains than predicted from the reliability.

Feed intake in housed dairy cows: validation of a three-dimensional camera-based feed intake measurement system

Measuring feed intake accurately is crucial to determine feed efficiency and for genetic selection. A system using three-dimensional (3D) cameras and deep learning algorithms can measure the volume of feed intake in dairy cows, but for now, the system has not been validated for feed intake expressed as weight of feed. The aim of this study was to validate the weight of feed intake predicted from the 3D cameras with the actual measured weight. It was hypothesised that diet-specific coefficients are necessary for predicting changes in weight, that the relationship between weight and volume is curvilinear throughout the day, and that manually pushing the feed affects this relationship. Twenty-four lactating Danish Holstein cows were used in a cross-over design with four dietary treatments, 2 × 2 factorial arranged with either grass-clover silage or maize silage as silage factor, and barley or dried beet pulp as concentrate factor. Cows were adapted to the diets for 11 d, and for 3 d to tie-stall housing before camera measurements. Six cameras were used for recording, each mounted over an individual feeding platform equipped with a weight scale. When building the predictive models, four cameras were used for training, and the remaining two for testing the prediction of the models. The most accurate predictions were found for the average feed intake over a period when using the starting density of the feed pile, which resulted in the lowest errors, 6% when expressed as RMSE and 5% expressed as mean absolute error. A model including curvilinear effects of feed volume and the impact of manual feed pushing was used on a dataset including daily time points. When cross-validating, the inclusion of a curvilinear effect and a feed push effect did not improve the accuracy of the model for neither the feed pile nor the feed removed by the cow between consecutive time points. In conclusion, measuring daily feed intake from this 3D camera system in the present experimental setup could be accomplished with an acceptable error (below 8%), but the system should be improved for individual meal intake measurements if these measures were to be implemented.

Short communication: Correlation of methane production, intensity, and yield with residual feed intake throughout lactation in Holstein cows

The environmental impact of dairy production can be reduced in several ways, including increasing feed efficiency and reducing methane (CH4) emissions. There is no consensus on their relationship. This study aimed at estimating the correlations between residual feed intake (RFI) and CH4 emissions expressed in g/d methane production (MeP), g/kg of fat- and protein-corrected milk methane intensity (MeI), or g/kg of DM intake methane yield (MeY) throughout lactation. We collected CH4 data using GreenFeed devices from 107 Holstein cows, as well as production and intake phenotypes. RFI was predicted from DM intake, fat- and protein-corrected milk, BW, and body condition score. Five-trait random regression models were used to estimate the individual variance components of the CH4 and production traits, which were used to calculate the correlations between RFI and CH4 traits throughout lactation. We found positive correlations of RFI with MeP and MeI ranging from 0.05 to 0.47 throughout the lactation. Correlations between RFI and MeY are low and vary from positive to negative, ranging from -0.18 to 0.17. Both MeP and MeI are favorably correlated with RFI, as is MeY during the first half of lactation. These correlations are mostly favorable for genetic selection, but the confirmation of these results is needed with genetic correlations over a larger dataset.

The Resilient Dairy Genome Project-A general overview of methods and objectives related to feed efficiency and methane emissions

The Resilient Dairy Genome Project (RDGP) is an international large-scale applied research project that aims to generate genomic tools to breed more resilient dairy cows. In this context, improving feed efficiency and reducing greenhouse gases from dairy is a high priority. The inclusion of traits related to feed efficiency (e.g., dry matter intake [DMI]) or greenhouse gases (e.g., methane emissions [CH4]) relies on available genotypes as well as high quality phenotypes. Currently, 7 countries (i.e., Australia, Canada, Denmark, Germany, Spain, Switzerland, and United States) contribute with genotypes and phenotypes including DMI and CH4. However, combining data are challenging due to differences in recording protocols, measurement technology, genotyping, and animal management across sources. In this study, we provide an overview of how the RDGP partners address these issues to advance international collaboration to generate genomic tools for resilient dairy. Specifically, we describe the current state of the RDGP database, data collection protocols in each country, and the strategies used for managing the shared data. As of February 2022, the database contains 1,289,593 DMI records from 12,687 cows and 17,403 CH4 records from 3,093 cows and continues to grow as countries upload new data over the coming years. No strong genomic differentiation between the populations was identified in this study, which may be beneficial for eventual across-country genomic predictions. Moreover, our results reinforce the need to account for the heterogeneity in the DMI and CH4 phenotypes in genomic analysis.

Genetic parameters for feed intake and body weight in dairy cattle using high-throughput 3-dimensional cameras in Danish commercial farms

Recording complex phenotypes on a large scale is becoming possible with the incorporation of recently developed new technologies. One of these new technologies is the use of 3-dimensional (3D) cameras on commercial farms to measure feed intake and body weight (BW) daily. Residual feed intake (RFI) has been proposed as a proxy for feed efficiency in several species, including cattle, pigs, and poultry. Dry matter intake (DMI) and BW records are required to calculate RFI, and the use of this new technology will help increase the number of individual records more efficiently. The aim of this study was to estimate genetic parameters (including genetic correlations) for DMI and BW obtained by 3D cameras from 6,000 cows in commercial farms from the breeds Danish Holstein, Jersey, and Nordic Red. Additionally, heritabilities per parity and genetic correlations among parities were estimated for DMI and BW in the 3 breeds. Data included 158,000 weekly records of DMI and BW obtained between 2019 and 2022 on 17 commercial farms. Estimated heritability for DMI ranged from 0.17 to 0.25, whereas for BW they ranged from 0.44 to 0.58. The genetic correlations between DMI and BW were moderately positive (0.58-0.65). Genetic correlations among parities in both traits were highly correlated in the 3 breeds, except for DMI between first parity and late parities in Holstein where they were down to 0.62. Based on these results, we conclude that DMI and BW phenotypes measured by 3D cameras are heritable for the 3 dairy breeds and their heritabilities are comparable to those obtained by traditional methods (scales and feed bins). The high heritabilities and correlations of 3D measurements with the true trait in previous studies demonstrate the potential of this new technology for measuring feed intake and BW in real time. In conclusion, 3D camera technology has the potential to become a valuable tool for automatic and continuous recording of feed intake and BW on commercial farms.

Novel genetic parameters for genetic residual feed intake in dairy cattle using time series data from multiple parities and countries in North America and Europe

Residual feed intake is viewed as an important trait in breeding programs that could be used to enhance genetic progress in feed efficiency. In particular, improving feed efficiency could improve both economic and environmental sustainability in the dairy cattle industry. However, data remain sparse, limiting the development of reliable genomic evaluations across lactation and parity for residual feed intake. Here, we estimated novel genetic parameters for genetic residual feed intake (gRFI) across the first, second, and third parity, using a random regression model. Research data on the measured feed intake, milk production, and body weight of 7,379 cows (271,080 records) from 6 countries in 2 continents were shared through the Horizon 2020 project Genomic Management Tools to Optimise Resilience and Efficiency, and the Resilient Dairy Genome Project. The countries included Canada (1,053 cows with 47,130 weekly records), Denmark (1,045 cows with 72,760 weekly records), France (329 cows with 16,888 weekly records), Germany (938 cows with 32,614 weekly records), the Netherlands (2,051 cows with 57,830 weekly records), and United States (1,963 cows with 43,858 weekly records). Each trait had variance components estimated from first to third parity, using a random regression model across countries. Genetic residual feed intake was found to be heritable in all 3 parities, with first parity being predominant (range: 22-34%). Genetic residual feed intake was highly correlated across parities for mid- to late lactation; however, genetic correlation across parities was lower during early lactation, especially when comparing first and third parity. We estimated a genetic correlation of 0.77 ± 0.37 between North America and Europe for dry matter intake at first parity. Published literature on genetic correlations between high input countries/continents for dry matter intake support a high genetic correlation for dry matter intake. In conclusion, our results demonstrate the feasibility of estimating variance components for gRFI across parities, and the value of sharing data on scarce phenotypes across countries. These results can potentially be implemented in genetic evaluations for gRFI in dairy cattle.

The relationship between feed efficiency and behaviour differs between lactating Holstein and Jersey cows

In dairy production, high feed efficiency (FE) is important to reduce feed costs and negative impacts of milk production on the climate and environment, yet little is known about the relationship between FE, eating behaviour and activity. This research communication describes how cows differing in FE, expressed as daily energy corrected milk production per unit of feed intake, differed in eating behaviour and activity. We used data from a study of 253 lactations obtained from 97 Holstein and 91 Jersey cows milked in an automatic milking system. Automated feed troughs recorded feed intake behaviour and cows wore a sensor that recorded activity from 5 to 200 d in milk (DIM). We used a mixed linear model to estimate random solutions for individual cows for traits of steps, lying and eating behaviour and calculated their correlation with FE during four periods (5-35, 36-75, 76-120 and 121-200 DIM). Separate analyses were performed for each breed and period. We found that individual level correlations between FE and behaviour traits were stronger in Jersey than in Holstein cows. Eating rate correlated weakly negatively to FE in Holstein cows and more strongly so in Jersey cows, such that efficient Jerseys were slower eaters. The physical activity of Jersey cows was weakly and negatively correlated to FE, but this was not the case in Holstein cows. We conclude that eating rate was consistently negatively associated with FE throughout lactation for Jersey cows, but not for Holstein cows.

Heritability and genetic correlations between enteric methane production and concentration recorded by GreenFeed and sniffers on dairy cows

To reduce methane (CH₄) emissions of dairy cows by animal breeding, CH₄ measurements have to be recorded on thousands of individual cows. Currently, several techniques are used to phenotype cows for CH₄, differing in costs and applicability. However, there is uncertainty about the agreement between techniques. To judge the similarity and repeatability between measurements of different recording techniques, the repeatability, heritability, and genetic correlation are useful metrics. Therefore, our objective was to estimate (1) the repeatability and heritability for CH₄ and carbon dioxide production recorded by GreenFeed (GF) and for CH₄ and carbon dioxide concentration measured by cost-effective but less accurate sniffers, and (2) the genetic correlation between CH₄ recorded with these 2 different on farm and high throughput techniques. Data were available from repeated measurements of CH₄ production (grams/day) by GF units and of CH₄ concentration (ppm) by sniffers, recorded on commercial dairy farms in the Netherlands. The final data comprised 24,284 GF daily means from 822 cows, 170,826 sniffer daily means from 1,800 cows, and 1,786 daily means from 75 cows by both GF and sniffer (in the same period). Additionally, CH₄ records were averaged per week. For daily and weekly mean GF CH₄ the heritabilities were 0.19 ± 0.02 and 0.33 ± 0.04, and for daily and weekly mean sniffer CH₄ the heritabilities were similar and were 0.18 ± 0.01 and 0.32 ± 0.02, respectively. Phenotypic correlations between GF CH₄ production and sniffer CH₄ concentration were moderate (0.39 ± 0.03 for daily means and 0.37 ± 0.05 for weekly means). However, genetic correlations were high; 0.71 ± 0.13 for daily means and 0.76 ± 0.15 for weekly means. The high genetic correlation indicates that selection on low CH₄ concentrations (ppm) recorded by the cost-effective sniffer method, will result in reduced CH₄ production (grams/day) as recorded with GF.

Implementation of Feed Efficiency in Iranian Holstein Breeding Program

This study aimed to evaluate the economic impact of improving feed efficiency on breeding objectives for Iranian Holsteins. Production and economic data from seven dairy herds were used to estimate the economic values of different traits, and a meta-analysis was conducted to analyze the genetic relationships between feed efficiency and other traits. Economic weights were calculated for various traits, with mean values per cow and per year across herds estimated at USD 0.34/kg for milk yield, USD 6.93/kg for fat yield, USD 5.53/kg for protein yield, USD −1.68/kg for dry matter intake, USD −1.70/kg for residual feed intake, USD 0.47/month for productive life, and USD −2.71/day for days open. The Iranian selection index was revised to improve feed efficiency, and the feed efficiency sub-index (FE$) introduced by the Holstein Association of the United States of America was adopted to reflect Iran’s economic and production systems. However, there were discrepancies between Iranian and US genetic coefficients in the sub-index, which could be attributed to differences in genetic and phenotypic parameters, as well as the economic value of each trait. More accurate estimates of economic values for each trait in FE$ could be obtained by collecting dry matter intake from Iranian herds and conducting genetic evaluations for residual feed intake.

Detection of potential functional variants based on systems-biology: the case of feed efficiency in beef cattle

BACKGROUND

Potential functional variants (PFVs) can be defined as genetic variants responsible for a given phenotype. Ultimately, these are the best DNA markers for animal breeding and selection, especially for polygenic and complex phenotypes. Herein, we described the identification of PFVs for complex phenotypes (in this case, Feed Efficiency in beef cattle) using a systems-biology driven approach based on RNA-seq data from physiologically relevant organs.

RESULTS

The systems-biology coupled with deep molecular phenotyping by RNA-seq of liver, muscle, hypothalamus, pituitary, and adrenal glands of animals with high and low feed efficiency (FE) measured by residual feed intake (RFI) identified 2,000,936 uniquely variants. Among them, 9986 variants were significantly associated with FE and only 78 had a high impact on protein expression and were considered as PFVs. A set of 169 significant uniquely variants were expressed in all five organs, however, only 27 variants had a moderate impact and none of them a had high impact on protein expression. These results provide evidence of tissue-specific effects of high-impact PFVs. The PFVs were enriched (FDR < 0.05) for processing and presentation of MHC Class I and II mediated antigens, which are an important part of the adaptive immune response. The experimental validation of these PFVs was demonstrated by the increased prediction accuracy for RFI using the weighted G matrix (ssGBLUP+wG; Acc = 0.10 and b = 0.48) obtained in the ssGWAS in comparison to the unweighted G matrix (ssGBLUP; Acc = 0.29 and b = 1.10).

CONCLUSION

Here we identified PFVs for FE in beef cattle using a strategy based on systems-biology and deep molecular phenotyping. This approach has great potential to be used in genetic prediction programs, especially for polygenic phenotypes.

Selecting for Feed Efficient Cows Will Help to Reduce Methane Gas Emissions

In the last decade, several countries have included feed efficiency (as residual feed intake; RFI) in their breeding goal. Recent studies showed that RFI is favorably correlated with methane emissions. Thus, selecting for lower emitting animals indirectly through RFI could be a short-term strategy in order to achieve the intended reduction set by the EU Commission (-55% for 2030). The objectives were to 1) estimate genetic parameters for six methane traits, including genetic correlations between methane traits, production, and feed efficiency traits, 2) evaluate the expected correlated response of methane traits when selecting for feed efficiency with or without including methane, 3) quantify the impact of reducing methane emissions in dairy cattle using the Danish Holstein population as an example. A total of 26,664 CH4 breath records from 647 Danish Holstein cows measured over 7 years in a research farm were analyzed. Records on dry matter intake (DMI), body weight (BW), and energy corrected milk (ECM) were also available. Methane traits were methane concentration (MeC, ppm), methane production (MeP; g/d), methane yield (MeY; g CH4/kg DMI), methane intensity (MeI; g CH4/kg ECM), residual methane concentration (RMeC), residual methane production (RMeP, g/d), and two definitions of residual feed intake with or without including body weight change (RFI1, RFI2). The estimated heritability of MeC was 0.20 ± 0.05 and for MeP, it was 0.21 ± 0.05, whereas heritability estimates for MeY and MeI were 0.22 ± 0.05 and 0.18 ± 0.04, and for the RMeC and RMeP, they were 0.23 ± 0.06 and 0.16 ± 0.02, respectively. Genetic correlations between methane traits ranged from moderate to highly correlated (0.48 ± 0.16–0.98 ± 0.01). Genetic correlations between methane traits and feed efficiency were all positive, ranging from 0.05 ± 0.20 (MeI-RFI2) to 0.76 ± 0.09 (MeP-RFI2). Selection index calculations showed that selecting for feed efficiency has a positive impact on reducing methane emissions’ expected response, independently of the trait used (MeP, RMeP, or MeI). Nevertheless, adding a negative economic value for methane would accelerate the response and help to reach the reduction goal in fewer generations. Therefore, including methane in the breeding goal seems to be a faster way to achieve the desired methane emission reductions in dairy cattle.

Genetic parameters for repeatedly recorded enteric methane concentrations of dairy cows

Animal breeding techniques offer potential to reduce enteric emissions of ruminants to lower the environmental impact of dairy farming. The aim of this study was to estimate the heritability and repeatability of methane (CH₄) concentrations, using the largest data set from long-term repeatedly recorded CH₄ on cows to date, and to evaluate (1) the accuracy of breeding values for different CH₄ traits, including using visits or weekly means, and (2) recording strategies (with varying numbers of records and recorded daughters per sire). The data comprised of long-term recording of CH₄ and carbon dioxide (CO₂), from 1,746 Holstein Friesian cows, on 14 commercial dairy farms throughout the Netherlands. Emissions were recorded in 10- to 35-s intervals, between 64 and 436 d, depending on farms. From each robot visit, CH₄ and CO₂ concentrations were summarized into various traits, averaged per visit and per week: mean, median, mean log, and mean CH₄/CO₂ ratio. Genetic parameters were estimated with animal repeatability models, using a restricted maximum likelihood procedure, and a relationship matrix based on genotypes and pedigree. The heritability was equal for mean and median CH₄ per visit (0.13) but lower for logCH₄ and CH₄/CO₂ (0.07 and 0.01, respectively). Phenotypic and genetic correlations were high (≥0.78) between the CH₄ traits, apart from the genetic correlations with the CH₄/CO₂ trait, which were negative. To achieve a minimum reliability of 50% for the estimated breeding value of a bull, 25 records on mean CH₄, measured on 10 different daughters, were sufficient. Although the heritability and repeatability were higher for weekly (0.32 and 0.68, respectively) than for visit mean CH₄ (0.13 and 0.30, respectively), the reliabilities of estimated breeding values from visit or weekly means were equal; thus, we found no advantage in averaging records to weekly means for genetic evaluations.

Translational multi-omics microbiome research for strategies to improve cattle production and health

Cattle microbiome plays a vital role in cattle growth and performance and affects many economically important traits such as feed efficiency, milk/meat yield and quality, methane emission, immunity and health. To date, most cattle microbiome research has focused on metataxonomic and metagenomic characterization to reveal who are there and what they may do, preventing the determination of the active functional dynamics in vivo and their causal relationships with the traits. Therefore, there is an urgent need to combine other advanced omics approaches to improve microbiome analysis to determine their mode of actions and host-microbiome interactions in vivo. This review will critically discuss the current multi-omics microbiome research in beef and dairy cattle, aiming to provide insights on how the information generated can be applied to future strategies to improve production efficiency, health and welfare, and environment-friendliness in cattle production through microbiome manipulations.

Genome-wide association study for methane emission traits in Danish Holstein cattle

Selecting for lower methane emitting cows requires insight into the most biologically relevant phenotypes for methane emission, which are close to the breeding goal. Several methane phenotypes have been suggested over the last decade. However, the (dis)similarity of their underlying genetic architecture and correlation structures are poorly understood. Therefore, the objective of this study was to test association of SNP and genomic regions through GWAS on 8 CH₄ emission traits in Danish Holstein cattle. The traits studied were methane concentration (MeC; ppm), methane production (MeP ; g/d), 2 definitions of residual methane (RMETc and RMETp: MeC and MeP regressed on metabolic body weight and energy-corrected milk, respectively), 2 definitions of methane intensity (MeI; MeIc = MeC/ECM and MeIp = MeP/ECM); 2 definitions of methane yield per kilogram of dry matter intake (MeY; MeYc = MeC/dry matter intake and MeYp = MeP/dry matter intake). A total of 1,962 cows with genotypes (Illumina BovineSNP50 Chip or Eurogenomic custom SNP chip) and repeated records of the above-mentioned 8 methane traits were analyzed. Strong associations were found with 3 traits (MeC, MeP, and MeYc) on chromosome 13 and with 5 traits (MeC, MeP, MeIp, MeYp, and MeYc) on chromosome 26. For MeIc, MeIp, RMETc, MeYc, and MeYp, some suggestive association signals were identified on chromosome 1. Genomic segments of 1 Mbp (n = 2,525) were tested for their association with these traits, which identified between 33 to 54 significantly associated regions. In a pairwise comparison, MeC and MeP were the traits that shared the highest number of significant segments (17). The same trend was observed when comparing SNP significantly associated with the traits MeC and MeP shared from 23 to 25 SNP (most of which were located in chromosomes 11, 13, and 26). Based on our results on GWAS and genetic correlations, we conclude that MeC is (genetically) more closely linked to MeP than any of the other methane traits analyzed.

Estimating methane coefficients to predict the environmental impact of traits in the Australian dairy breeding program

The dairy industry has been scrutinized for the environmental impact associated with rearing and maintaining cattle for dairy production. There are 3 possible opportunities to reduce emissions through genetic selection: (1) a direct methane trait, (2) a reduction in replacements, and (3) an increase in productivity. Our aim was to estimate the independent effects of traits in the Australian National Breeding Objective on the gross methane production and methane intensity (EI) of the Australian dairy herd of average genetic potential. Based on similar published research, the traits determined to have an effect on emissions include production, fertility, survival, health, and feed efficiency. The independent effect of each trait on the gross emissions produced per animal due to genetic improvement and change in EI due to genetic improvement (intensity value, IV) were estimated and compared. Based on an average Australian dairy herd, the gross emissions emitted per cow per year were 4,297.86 kg of carbon dioxide equivalents (CO₂-eq). The annual product output, expressed in protein equivalents (protein-eq), and EI per cow were 339.39 kg of protein-eq and 12.67 kg of CO₂-eq/kg of protein-eq, respectively. Of the traits included in the National Breeding Objective, genetic progress in survival and feed saved were consistently shown to result in a favorable environmental impact. Conversely, production traits had an unfavorable environmental impact when considering gross emissions, and favorable when considering EI. Fertility had minimal impact as its effects were primarily accounted for through survival. Mastitis resistance only affected IV coefficients and to a very limited extent. These coefficients may be used in selection indexes to apply emphasis on traits based on their environmental impact, as well as applied by governments and stakeholders to track trends in industry emissions. Although initiatives are underway to develop breeding values to reduce methane by combining small methane data sets internationally, alternative options to reduce emissions by utilizing selection indexes should be further explored.

Selective breeding as a mitigation tool for methane emissions from dairy cattle

The global livestock sector, particularly ruminants, contributes substantially to the total anthropogenic greenhouse gases. Management and dietary solutions to reduce enteric methane (CH₄) emissions are extensively researched. Animal breeding that exploits natural variation in CH₄ emissions is an additional mitigation solution that is cost-effective, permanent, and cumulative. We quantified the effect of including CH₄ production in the Dutch breeding goal using selection index theory. The current Dutch national index contains 15 traits, related to milk yield, longevity, health, fertility, conformation and feed efficiency. From the literature, we obtained a heritability of 0.21 for enteric CH₄ production, and genetic correlations of 0.4 with milk lactose, protein, fat and DM intake. Correlations between enteric CH₄ production and other traits in the breeding goal were set to zero. When including CH₄ production in the current breeding goal with a zero economic value, CH₄ production increases each year by 1.5 g/d as a correlated response. When extrapolating this, the average daily CH₄ production of 392 g/d in 2018 will increase to 442 g/d in 2050 (+13%). However, expressing the CH₄ production as CH₄ intensity in the same period shows a reduction of 13%. By putting economic weight on CH₄ production in the breeding goal, selective breeding can reduce the CH₄ intensity even by 24% in 2050. This shows that breeding is a valuable contribution to the whole set of mitigation strategies that could be applied in order to achieve the goals for 2050 set by the EU. If the decision is made to implement animal breeding strategies to reduce enteric CH₄ production, and to achieve the expected breeding impact, there needs to be a sufficient reliability of prediction. The only way to achieve that is to have enough animals phenotyped and genotyped. The power calculations offer insights into the difficulties that will be faced in trying to record enough data. Recording CH₄ data on 100 farms (with on average 150 cows each) for at least 2 years is required to achieve the desired reliability of 0.40 for the genomic prediction.

Breeding for reduced methane emission and feed-efficient Holstein cows: An international response

Selecting for lower methane (CH₄) emitting animals is one of the best approaches to reduce CH₄ given that genetic progress is permanent and cumulative over generations. As genetic selection requires a large number of animals with records and few countries actively record CH₄, combining data from different countries could help to expedite accurate genetic parameters for CH₄ traits and build a future genomic reference population. Additionally, if we want to include CH₄ in the breeding goal, it is important to know the genetic correlations of CH₄ traits with other economically important traits. Therefore, the aim of this study was first to estimate genetic parameters of 7 suggested methane traits, as well as genetic correlations between methane traits and production, maintenance, and efficiency traits using a multicountry database. The second aim was to estimate genetic correlations within parities and stages of lactation for CH₄. The third aim was to evaluate the expected response of economically important traits by including CH₄ traits in the breeding goal. A total of 15,320 methane production (MeP, g/d) records from 2,990 cows belonging to 4 countries (Canada, Australia, Switzerland, and Denmark) were analyzed. Records on dry matter intake (DMI), body weight (BW), body condition score, and milk yield (MY) were also available. Additional traits such as methane yield (MeY; g/kg DMI), methane intensity (MeI; g/kg energy-corrected milk), a genetic standardized methane production, and 3 definitions of residual methane production (g/d), residual feed intake, metabolic BW (MBW), BW change, and energy-corrected milk were calculated. The estimated heritability of MeP was 0.21, whereas heritability estimates for MeY and MeI were 0.30 and 0.38, and for the residual methane traits heritability ranged from 0.13 to 0.16. Genetic correlations between different methane traits were moderate to high (0.41 to 0.97). Genetic correlations between MeP and economically important traits ranged from 0.29 (MY) to 0.65 (BW and MBW), being 0.41 for DMI. Selection index calculations showed that residual methane had the most potential for inclusion in the breeding goal when compared with MeP, MeY, and MeI, as residual methane allows for selection of low methane emitting animals without compromising other economically important traits. Inclusion of residual feed intake in the breeding goal could further reduce methane, as the correlation with residual methane is moderate and elicits a favorable correlated response. Adding a negative economic value for methane could facilitate a substantial reduction in methane emissions while maintaining an increase in milk production.

Residual carbon dioxide as an index of feed efficiency in lactating dairy cows

High feed costs make feed conversion efficiency a desirable target for genetic improvement. Residual feed intake (RFI), calculated as the difference between observed and predicted intake, is a commonly used estimate of feed efficiency. However, determination of feed efficiency in dairy herds is challenging due to difficulties in measuring feed intake of individual animals reliably. Using residual CO₂ (RCO₂) production as an estimate of feed efficiency would allow ranking the cows according to feed efficiency, provided that CO₂ production is closely related to heat production and feed intake. The objective of this study was to evaluate the potential of RCO₂ as an index of feed efficiency using data from respiration calorimetry studies (289 cow per period observations). Heat production was precisely predicted from CO₂ production [root mean square error (RMSE)] adjusted for random effects was 1.5% of observed mean]. Dry matter intake (DMI) was better predicted from energy-corrected milk (ECM) yield and CO₂ production than from ECM yield and body weight in the model (adjusted RSME = 0.92 vs. 1.39 kg/d). Residual CO₂ production estimated as the difference between actual CO₂ production and that predicted from ECM yield, metabolic body weight was closely related to RFI (adjusted RMSE = 0.42) that was calculated as the difference between actual DMI and that predicted from ECM, metabolic body weight, and energy balance (EB). When the cows were categorized in 3 groups of equal sizes on the basis of RCO₂ (low, medium, and high), low RCO₂ cows had lower DMI, RFI, methane production and intensity (g/kg ECM), and heat production, but higher efficiency of metabolizable energy utilization for lactation than high RCO₂ cows. When RFI was predicted from RCO₂, the residuals (observed - predicted) were negatively related to EB and digestibility. Predicting RFI with a 2-variable model based on RCO₂ and digestibility, adjusted RMSE decreased to 0.23 kg/d, and residuals were not significantly related to EB. The cows in low RCO₂ group had a higher energy digestibility than the cows in the high RCO₂ group, and differences in EB were observed between the groups. Error of the model predicting residual ECM production from RCO₂ was 1.41 kg/d. The residuals were positively related to ECM yield and energy digestibility. Predicting residual ECM from RCO₂ and ECM yield decreased adjusted RMSE to 1.07 kg/d, and further to 0.78 kg/d when digestibility was included in the 2-variable model. It is concluded that RCO₂ has a potential for ranking individual cows based on feed efficiency.

Genetic mechanisms underlying feed utilization and implementation of genomic selection for improved feed efficiency in dairy cattle

The economic importance of genetically improving feed efficiency has been recognized by cattle producers worldwide. It has the potential to considerably reduce costs, minimize environmental impact, optimize land and resource use efficiency, and improve the overall cattle industry’s profitability. Feed efficiency is a genetically complex trait that can be described as units of product output (e.g., milk yield) per unit of feed input. The main objective of this review paper is to present an overview of the main genetic and physiological mechanisms underlying feed utilization in ruminants and the process towards implementation of genomic selection for feed efficiency in dairy cattle. In summary, feed efficiency can be improved via numerous metabolic pathways and biological mechanisms through genetic selection. Various studies have indicated that feed efficiency is heritable, and genomic selection can be successfully implemented in dairy cattle with a large enough training population. In this context, some organizations have worked collaboratively to do research and develop training populations for successful implementation of joint international genomic evaluations. The integration of “-omics” technologies, further investments in high-throughput phenotyping, and identification of novel indicator traits will also be paramount in maximizing the rates of genetic progress for feed efficiency in dairy cattle worldwide.

Multitrait genomic prediction of methane emissions in Danish Holstein cattle

In dairy cattle, selecting for lower methane-emitting animals is one of the new challenges of this decade. However, genetic selection requires a large number of animals with records to get accurate estimated breeding values (EBV). Given that CH₄ records are scarce, the use of information on routinely recorded and highly correlated traits with CH₄ has been suggested to increase the accuracy of genomic EBV (GEBV) through multitrait (genomic) prediction. Therefore, the objective of this study was to evaluate accuracies of prediction of GEBV for CH₄ by including or omitting CH₄, energy-corrected milk (ECM), and body weight (BW) as well as genotypic information in multitrait analyses across 2 methods: BLUP and single-step genomic BLUP (SSGBLUP). A total of 2,725 cows with CH₄ concentration in breath (14,125 records), BW (61,667 records), and ECM (61,610 records) were included in the analyses. Approximately 2,000 of these cows were genotyped or imputed to 50K. Ten cross-validation groups were formed by randomly grouping paternal half-sibs. Five scenarios were performed: (1) base scenario with only CH₄ information; (2) without CH₄, but with information from BW, ECM, or BW+ECM only in reference population; (3) without CH₄, but with information from BW, ECM, or BW+ECM in both validation and reference population; (4) with CH₄ information and BW, ECM, or BW+ECM information only in the reference population; and (5) with CH₄ information and BW, ECM, or BW+ECM information in both validation and reference population. As a result, for each method (BLUP, SSGBLUP), 13 sub-scenarios were performed, 1 from scenario 1, and 3 for each of the subsequent 4 scenarios. The average accuracy of GEBV for CH₄ in the base scenario was 0.32 for BLUP and 0.42 for SSGBLUP, and it ranged from 0.10 in scenario 2 to 0.78 in scenario 5 across methods. In terms of bias, the base scenario 1 was unbiased for SSGBLUP; similar results were achieved with scenario 5. Including information on ECM increased the accuracy of GEBV for CH₄ by up to 61%, whereas adding information on both traits (BW and ECM) increased the accuracy by up to 90%. Scenarios that did not include CH₄ in the reference population had the lowest correlations (0.17-0.33) with single-trait CH₄ GEBV, and scenarios with CH₄ in the reference population had the highest correlations (0.41-0.81). Thus, failure to include CH₄ in future reference populations results in predicted CH₄ GEBV, which cannot be used in practical selection. Therefore, recording CH₄ in more animals remains a priority. Finally, multiple-trait genomic prediction using routinely recorded BW and ECM leads to higher prediction accuracies than traditional single-trait genomic prediction for CH₄ and is a viable solution for increasing the accuracies of GEBV for scarcely recorded CH₄ in practice.