Immune cells phenotype and bioenergetic measures in CD4+ T cells differ between high and low feed efficient dairy cows

Identifying sources of variance that contribute to residual feed intake (RFI) can aid in improving feed efficiency. The objectives of this study were to investigate immune cells phenotype and bioenergetic measures in CD4+ T cells in low feed efficient (LE) and high feed efficient (HE) dairy cows. Sixty-four Holstein cows were enrolled at 93 ± 22 days in milk (DIM) and monitored for 7 weeks to measure RFI. Cows with the highest RFI (LE; n = 14) or lowest RFI (HE; n = 14) were selected to determine immune cells phenotype using flow cytometry. Blood was sampled in the same LE and HE cows at 234 ± 22 DIM to isolate peripheral blood mononuclear cells, followed by magnetic separation of CD4+ T lymphocytes using bovine specific monoclonal antibodies. The metabolic function of isolated CD4+ T lymphocytes was evaluated under resting and activated states. An increased expression of CD62L+ cells within CD8+ T lymphocytes and CD21+ B lymphocytes was observed in HE cows compared to LE cows. CD4+ T lymphocytes of HE cows exhibited an increased mitochondrial and glycolytic activity in resting and activated states compared to LE cows. These data suggest that immune cells in HE cows exhibit an increased metabolic function, which might influence nutrient partitioning and utilization and serve as a source of variation in feed efficiency that warrants future investigation.

Estimates of genetic parameters for rumination time, feed efficiency, and methane production traits in first-lactation Holstein cows

The large-scale recording of traits such as feed efficiency (FE) and methane emissions (ME) for use in genetic improvement programs is complex, costly, and time-consuming. Therefore, heritable traits that can be continuously recorded in dairy herds and are correlated with FE and ME traits could provide useful information for genetic evaluation. Rumination time has been suggested to be associated with FE, methane production (MeP; ME in g/d), and production traits at the phenotypic level. Therefore, the objective of this study was to investigate the genetic relationships among rumination time (RT), FE, methane and production traits using 7,358 records from 656 first-lactation Holstein cows. The estimated heritabilities were moderate for RT (0.45 ± 0.14), MeP (0.36 ± 0.12), milk yield (0.40 ± 0.08), fat yield (0.29 ± 0.06), protein yield (0.32 ± 0.07), and energy-corrected milk (0.28 ± 0.07), but were low and nonsignificant for FE (0.15 ± 0.07), which was defined as the residual of the multiple linear regression of DMI on energy-corrected milk and metabolic body weight. A favorable negative genetic correlation was estimated between RT and MeP (-0.53 ± 0.24), whereas a positive favorable correlation was estimated between RT and energy-corrected milk (0.49 ± 0.11). The estimated genetic correlation of RT with FE (-0.01 ± 0.17) was not significantly different from zero but showed a trend of a low correlation with dry matter intake (0.21 ± 0.13). These results indicate that RT is genetically associated with MeP and milk production traits, but high standard errors indicate that further analyses should be conducted to verify these findings when more data for RT, MeP, and FE become available.

Assessing different cross-validation schemes for predicting novel traits using sensor data: an application to dry matter intake and residual feed intake using milk spectral data

Feed efficiency is important for economic profitability of dairy farms; however, recording daily dry matter intakes (DMI) is expensive. Our objective was to investigate the potential use of milk mid-infrared (MIR) spectral data to predict proxy phenotypes for DMI based on different cross-validation schemes. We were specifically interested in comparisons between a model that included only MIR data (Model M1), a model that incorporated different energy sink predictors, such as body weight, body weight change, and milk energy (Model M2), and an extended model that incorporated both energy sinks and MIR data (Model M3). Models M2 and M3 also included various cow level variables (stage of lactation, age at calving, parity) such that any improvement in model performance from M2 to M3, whether through a smaller root mean squared error (RMSE) or a greater squared predictive correlation (R2), could indicate a potential benefit of MIR to predict residual feed intake. The data used in our study originated from a multi-institutional project on the genetics of feed efficiency in US Holsteins. Analyses were conducted on 2 different trait definitions based on different period lengths: averaged across weeks vs. averaged across 28-d. Specifically, there were 19,942 weekly records on 1,812 cows across 46 experiments or cohorts and 3,724 28-d records on 1,700 cows across 43 different cohorts. The cross-validation analyses involved 3 different k-fold schemes. First, a 10-fold cow-independent cross-validation was conducted whereby all records from any one cow were kept together in either training or test sets. Similarly, a 10-fold experiment-independent cross-validation kept entire experiments together whereas a 4-fold herd-independent cross-validation kept entire herds together in either training or test sets. Based on cow-independent cross-validation for both weekly and 28-d DMI, adding MIR predictors to energy sinks (Models M3 vs M2) significantly (P < 10-10) reduced average RMSE to 1.59 kg and increased average R2 to 0.89. However, adding MIR to energy sinks (M3) to predict DMI either within an experiment-independent or herd-independent cross-validation scheme seemed to demonstrate no merit (P > 0.05) compared with an energy sink model (M2) for either R2 or RMSE (respectively, 0.68 and 2.55 kg for M2 in herd-independent scheme). We further noted that with broader cross-validation schemes, i.e., from cow-independent to experiment-independent to herd-independent schemes, the mean and slope bias increased. Given that proxy DMI phenotypes for cows would need to be almost entirely generated in herds having no DMI or training data of their own, herd-independent cross-validation assessments of predictive performance should be emphasized. Hence, more research on predictive algorithms suitable for broader cross-validation schemes and a more earnest effort on calibration of spectrophotometers against each other should be considered.

Dry matter intake in US Holstein cows: exploring the genomic and phenotypic impact of milk components and body weight composite

Large data sets allow estimating feed required for individual milk components or body maintenance. Phenotypic regressions are useful for nutrition management, but genetic regressions are more useful in breeding programs. Dry matter intake (DMI) records from 8,513 lactations of 6,621 Holstein cows were predicted from phenotypes or genomic evaluations for milk components and body size traits. The mixed models also included days in milk, age-parity subclass, trial date, management group, and body weight change during 28- and 42-d feeding trials in mid-lactation. Phenotypic regressions of DMI on milk (0.014 ± 0.006), fat (3.06 ± 0.01), and protein (4.79 ± 0.25) were much less than corresponding genomic regressions (0.08 ± 0.03, 11.30 ± 0.47, and 9.35 ± 0.87) or sire genomic regressions multiplied by 2 (0.048 ± 0.04, 6.73 ± 0.94, and 4.98 ± 1.75). Thus, marginal feed costs as fractions of marginal milk revenue were higher from genetic than phenotypic regressions. According to the energy-corrected milk formula, fat production requires 69% more DMI than protein production. In the phenotypic regression, it was estimated that protein production requires 56% more DMI than fat. However, the genomic regression for the animal showed a difference of only 21% more DMI for protein compared with fat, while the sire genomic regressions indicated approximately 35% more DMI for fat than protein. Estimates of annual maintenance in kg DMI / kg body weight/lactation were similar from phenotypic regression (5.9 ± 0.14), genomic regression (5.8 ± 0.31), and sire genomic regression multiplied by 2 (5.3 ± 0.55) and are larger than those estimated by NASEM (2021) based on NEL equations. Multiple regressions on genomic evaluations for the 5 type traits in body weight composite (BWC) showed that strength was the type trait most associated with body weight and DMI, agreeing with the current BWC formula, whereas other traits were less useful predictors, especially for DMI. The Net Merit formula used to weight different genetic traits to achieve an economically optimal overall selection response was revised in 2021 to better account for these estimated regressions. To improve profitability, breeding programs should select smaller cows with negative residual feed intake that produce more milk, fat, and protein.

Longitudinal modeling of residual carbon dioxide and residual feed intake in the Nordic Red dairy cattle

Feed utilization efficiency is an important trait in dairy production playing a significant role in reducing feed costs and lowering methane emission. One of the metrics used to measure feed efficiency in dairy cows is residual feed intake (RFI). This metric requires routine measurement of feed intake. Since there is a positive high correlation between heat production and carbon dioxide (CO2) production on the one hand and heat production and efficiency on the other hand, residual carbon dioxide (RCO2) might be a useful metric to improve feed efficiency. The objectives of this study were to model the trajectories of RCO2 and RFI as well as to estimate their repeatabilities and correlations at different stages of lactation. Daily CO2 output and feed intake were recorded from 46 primiparous Nordic Red dairy cows using two Greenfeed Emissions Monitoring™ systems from 2 to 305 days in milk (DIM). Edited data comprised 5 995 daily averages. To calculate predicted values of CO2 and DM intake (DMI), prediction models were developed by fitting multiple regression models to observations. Subsequently, RCO2 and RFI were calculated by subtracting predicted values of CO2 and DMI from their corresponding actual observations. A random regression bivariate model was fitted to estimate repeatabilities and animal correlations within lactation at different DIMs between RCO2 and RFI traits. The model fitted included fixed effects of year-month of recording, lactation month, fixed regressions as well as random regressions for the animal effect. The residual variance was considered to be heterogeneous. Repeatabilities and animal correlations of RCO2 and RFI between selected DIM (for every 30 DIM i.e., 6, 36,…, 246 and 276) were calculated. Repeatability of RCO2 was high at the beginning of lactation (0.72 at DIM 6) and decreased around the peak of milk production (0.27 at DIM 96) and again increased gradually toward the end of lactation. Similarly, RFI also had high repeatability at the beginning (0.86 at DIM 6); however, it decreased in mid-lactation (0.37 at DIM 156) and then increased toward the end of lactation. Animal correlations between RCO2 and RFI were moderate to high on the same DIM and ranged from 0.37 to 0.88. Overall, we found that animals with higher CO2 production than expected also consume more DMI than expected, but the moderate correlation between RCO2 and RFI found in this study calls for more research to assess the potential of RCO2 to become a new feed efficiency metric.

Estimation of genetic parameters for feed efficiency traits using random regression models in dairy cattle

Feed efficiency has become an increasingly important research topic in recent years. As feed costs rise and the environmental impacts of agriculture become more apparent, improving the efficiency with which dairy cows convert feed to milk is increasingly important. However, feed intake is expensive to measure accurately on large populations, making the inclusion of this trait in breeding programs difficult. Understanding how the genetic parameters of feed efficiency and traits related to feed efficiency vary throughout the lactation period is valuable to gain understanding into the genetic nature of feed efficiency. This study used 121,226 dry matter intake (DMI) records, 120,500 energy-corrected milk (ECM) records, and 98,975 metabolic body weight (MBW) records, collected on 7,440 first-lactation Holstein cows from 6 countries (Canada, Denmark, Germany, Spain, Switzerland, and the United States), from January 2003 to February 2022. Genetic parameters were estimated using a multiple-trait random regression model with a fourth-order Legendre polynomial for all traits. Weekly phenotypes for DMI were re-parameterized using linear regressions of DMI on ECM and MBW, creating a measure of feed efficiency that was genetically corrected for ECM and MBW, referred to as genomic residual feed intake (gRFI). Heritability (SE) estimates varied from 0.15 (0.03) to 0.29 (0.02) for DMI, 0.24 (0.01) to 0.29 (0.03) for ECM, 0.55 (0.03) to 0.83 (0.05) for MBW, and 0.12 (0.03) to 0.22 (0.06) for gRFI. In general, heritability estimates were lower in the first stage of lactation compared with the later stages of lactation. Additive genetic correlations between weeks of lactation varied, with stronger correlations between weeks of lactation that were close together. The results of this study contribute to a better understanding of the change in genetic parameters across the first lactation, providing insight into potential selection strategies to include feed efficiency in breeding programs.

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.

An artificial intelligence approach of feature engineering and ensemble methods depicts the rumen microbiome contribution to feed efficiency in dairy cows

Genetic selection has remarkably helped U.S. dairy farms to decrease their carbon footprint by more than doubling milk production per cow over time. Despite the environmental and economic benefits of improved feed and milk production efficiency, there is a critical need to explore phenotypical variance for feed utilization to advance the long-term sustainability of dairy farms. Feed is a major expense in dairy operations, and their enteric fermentation is a major source of greenhouse gases in agriculture. The challenges to expanding the phenotypic database, especially for feed efficiency predictions, and the lack of understanding of its drivers limit its utilization. Herein, we leveraged an artificial intelligence approach with feature engineering and ensemble methods to explore the predictive power of the rumen microbiome for feed and milk production efficiency traits, as rumen microbes play a central role in physiological responses in dairy cows. The novel ensemble method allowed to further identify key microbes linked to the efficiency measures. We used a population of 454 genotyped Holstein cows in the U.S. and Canada with individually measured feed and milk production efficiency phenotypes. The study underscored that the rumen microbiome is a major driver of residual feed intake (RFI), the most robust feed efficiency measure evaluated in the study, accounting for 36% of its variation. Further analyses showed that several alpha-diversity metrics were lower in more feed-efficient cows. For RFI, [Ruminococcus] gauvreauii group was the only genus positively associated with an improved feed efficiency status while seven other taxa were associated with inefficiency. The study also highlights that the rumen microbiome is pivotal for the unexplained variance in milk fat and protein production efficiency. Estimation of the carbon footprint of these cows shows that selection for better RFI could reduce up to 5 kg of diet consumed per cow daily, potentially reducing up to 37.5% of CH4. These findings shed light that the integration of artificial intelligence approaches, microbiology, and ruminant nutrition can be a path to further advance our understanding of the rumen microbiome on nutrient requirements and lactation performance of dairy cows to support the long-term sustainability of the dairy community.

Consistency of dry matter intake in Holstein cows: Heritability estimates and associations with feed efficiency

Resilience can be defined as the capacity to maintain performance or bounce back to normal functioning after a perturbation, and studying fluctuations in daily feed intake may be an effective way to identify resilient dairy cows. Our goal was to develop new phenotypes based on daily dry matter intake (DMI) consistency in Holstein cows, estimate genetic parameters and genetic correlations with feed efficiency and milk yield consistency, and evaluate their relationships with production, longevity, health, and reproduction traits. Data consisted of 397,334 daily DMI records of 6,238 lactating Holstein cows collected from 2007 to 2022 at 6 research stations across the United States. Consistency phenotypes were calculated based on the deviations from expected daily DMI for individual cows during their respective feeding trials, which ranged from 27 to 151 d in duration. Expected values were derived from different models, including simple average, quadratic and cubic quantile regression with a 0.5 quantile, and locally estimated scatterplot smoothing (LOESS) regression with span parameters 0.5 and 0.7. We then calculated the log of variance (log-Var-DMI) of daily deviations for each model as the consistency phenotype. Consistency of milk yield was also calculated, as a reference, using the same methods (log-Var-Milk). Genetic parameters were estimated using an animal model, including lactation, days in milk and cohort as fixed effects, and animal as random effect. Relationships between log-Var-DMI and traits currently considered in the US national genetic evaluation were evaluated using Spearman's rank correlations between sires' breeding values. Heritability estimates for log-Var-DMI ranged from 0.11 ± 0.02 to 0.14 ± 0.02 across models. Different methods (simple average, quantile regressions, and LOESS regressions) used to calculate log-Var-DMI yielded very similar results, with genetic correlations ranging from 0.94 to 0.99. Estimated genetic correlations between log-Var-DMI and log-Var-Milk ranged from 0.51 to 0.62. Estimated genetic correlations between log-Var-DMI and feed efficiency ranged from 0.55 to 0.60 with secreted milk energy, from 0.59 to 0.63 with metabolic body weight, and from 0.26 to 0.31 with residual feed intake (RFI). Relationships between log-Var-DMI and the traits in the national genetic evaluation were moderate and positive correlations with milk yield (0.20 to 0.21), moderate and negative correlations with female fertility (-0.07 to -0.20), no significant correlations with health and longevity, and favorable correlations with feed efficiency (-0.23 to -0.25 with feed saved and 0.21 to 0.26 with RFI). We concluded that DMI consistency is heritable and may be an indicator of resilience. Cows with lower variation in the difference between actual and expected daily DMI (more consistency) may be more effective in maintaining performance in the face of challenges or perturbations, whereas cows with greater variation in observed versus expected daily DMI (less consistency) are less feed efficient and may be less resilient.

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Repeatabilities of individual measures of feed intake and body weight on in-house commercial dairy cattle using a 3-dimensional camera system

In this study a 3-dimensional (3D) camera system was set up to measure individual feed intake of dairy cows in a commercial in-house setting. The system was developed to identify the cows while eating, predict body weight based on the curvature of the back of the cow, and quantify the amount of feed eaten by the cow at each visit of eating. The identification of the cow was based on recognizing the patterns, colors, and curvatures of the back from a reference database obtained in a corridor after milking, where images were taken of all cows with a simultaneous reading of the electronic ear tag. Body weight is predicted using the curvatures of the back of the cow. Feed intake is quantified as the difference in surface of the feed a cow can reach before and after a visit is initiated. This estimate is in liters but converted to kilograms, using the density of the feed in the specific herd. A total of 9,142 cows were measured in 19 herds across 3 breeds: Jersey (2,513 cows), Red Dairy Cattle (2,813 cows), and Holstein (3,816 cows). Mean daily feed intake was higher for Red Dairy Cattle (61.72 kg) and Holstein (64.59 kg) than for Jersey (55.74 kg). Repeatability estimates for daily feed intake as a weekly average was 0.62, 0.65, and 0.63 for Jersey, Red Dairy, and Holstein cattle, respectively. Mean body weight was higher for Red Dairy (647.9 kg) and Holstein (683.8 kg) than for Jersey (469.6 kg). Repeatability estimates for body weight as a weekly average was 0.83, 0.85, and 0.88 for Jersey, Red Dairy, and Holstein, respectively. The perspectives in having such records available is huge both for the farmer and for the dairy industry. The records can both be used for improving management in farms on an individual cow level and herd level, but also for genetic evaluation and selection as well as testing feeding regimens. Feed intake can be measured on an individual level using a 3D camera system.

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.

Liver and Muscle Transcriptomes Differ in Mid-Lactation Cows Divergent in Feed Efficiency in the Presence or Absence of Supplemental Rumen-Protected Choline

Improving dairy cow feed efficiency is critical to the sustainability and profitability of dairy production, yet the underlying mechanisms that contribute to individual cow variation in feed efficiency are not fully understood. The objectives of this study were to (1) identify genes and associated pathways that are altered in cows with high- or low-residual feed intake (RFI) using RNA sequencing, and (2) determine if rumen-protected choline supplementation during mid-lactation would influence performance or feed efficiency. Mid-lactation (134 ± 20 days in milk) multiparous Holstein cows were randomly assigned to either supplementation of 0 g/d supplementation (CTL; n = 32) or 30 g/d of a rumen-protected choline product (RPC; 13.2 g choline ion; n = 32; Balchem Corp., New Hampton, NY, USA). Residual feed intake was determined as dry matter intake regressed on milk energy output, days in milk, body weight change, metabolic body weight, and dietary treatment. The 12 cows with the highest RFI (low feed efficient; LE) and 12 cows with the lowest RFI (high feed efficient; HE), balanced by dietary treatment, were selected for blood, liver, and muscle analysis. No differences in production or feed efficiency were detected with RPC supplementation, although albumin was greater and arachidonic acid tended to be greater in RPC cows. Concentrations of β-hydroxybutyrate were greater in HE cows. Between HE and LE, 268 and 315 differentially expressed genes in liver and muscle tissue, respectively, were identified through RNA sequencing. Pathway analysis indicated differences in cell cycling, oxidative stress, and immunity in liver and differences in glucose and fatty acid pathways in muscle. The current work indicates that unique differences in liver and muscle post-absorptive nutrient metabolism contribute to sources of variation in feed efficiency and that differences in amino acid and fatty acid oxidation, cell cycling, and immune function should be further examined.

Opportunities and limits of commercial farm data to study the genetic determinism of feed efficiency throughout lactation in dairy sheep

The collective economic and environmental interest of the whole dairy sheep sector is to reduce feed costs and the negative impact of milk production on the environment. Thus, this study focused on the characterisation and genetic selection potential of feed efficiency in the Lacaune breed. Estimates for feed efficiency in dairy ewes are limited, mainly due to a lack of individual feed intake measurements in the sheepfold or in the pasture. We estimated the genetic parameters for two approximated (not entirely based on individual data) feed efficiency traits (lactation feed conversion ratio (LFCR) and residual energy intake (REI)) and daily milk yield (DMY) at different stages of lactation and throughout lactation. The accuracy of the efficiency traits was first evaluated on samples from Lacaune dairy ewes that were monitored individually, especially for their feed intake. Then, feed efficiency estimation methods were applied on eight commercial farms corresponding to 4 680 Lacaune dairy ewes over two milk lactations (30 854 records). Animals were collectively (for a large part of feed intake) or individually (for milk performance and dynamics of body fat reserves) monitored at different lactation stages. The heritabilities of LFCR and REI were estimated over lactations at 0.10 ± 0.01 and 0.11 ± 0.01, respectively. High genetic correlations were observed between the two efficiency traits and milk production traits, with a genetic correlation between LFCR and DMY of 0.74 ± 0.04 and between REI and DMY of -0.79 ± 0.04. A strong influence of environmental factors such as farm, year of milk production and lactation stage affected the genetic link between REI and milk production traits. Efficiency values observed in early lactation when animals were bred in the sheepfold were less genetically correlated with values obtained later in lactation when animals were grass-fed. However, individual characterisation of feed efficiency remains difficult due to the collective feeding context in dairy ewe farms.

Estimating breeding values for feed efficiency in dairy cattle by regression on expected feed intake

The efficiency with which a dairy cow utilises feed for the various physiological and metabolic processes can be evaluated by metrics that contrast realised feed intake with expected feed intake. In this study, we presented a new metric - regression on expected feed intake (ReFI). This metric is based on the idea of regressing DM intake (DMI) on expected DMI using a random regression model, where energy requirement formulations are applied for the calculation of expected DMI covariables. We compared this new metric with the metrics residual feed intake (RFI) and genetic residual feed intake (gRFI), by applying them on 18 581 feed efficiency records from 654 primiparous Nordic Red dairy cows. We estimated variance components for the three metrics and their respective genetic correlations with intake and production traits. In addition, we examined the phenotypes of superior cows. With ReFI, we estimated for feed efficiency a higher genetic variation (4.7%) and heritability (0.23) compared to applying RFI or gRFI. The ReFI metric was genetically uncorrelated with DMI and negatively correlated within energy-corrected milk (ECM), whereas the RFI metric was genetically positively correlated with DMI and metabolic BW. The gRFI metric was genetically positively correlated with DMI and uncorrelated with energy sink traits. Overall, the estimated SE were large. The ReFI metric resulted in a different ranking of cows compared to those based on RFI or gRFI and was superior in selecting the most efficient animals. When the selection was based on ReFI breeding values, then the 10% most efficient cows produced 12.3% more ECM per unit metabolisable energy intake, whereas the corresponding values were only 4.3 or 5.9% when using RFI or gRFI breeding values, respectively. Based on ReFI, superior cows had also higher milk production, whereas based on RFI or gRFI milk production either decreased or was unaffected, respectively. The superiority of the ReFI metric in selecting efficient cows was due to a better modelling of the expected feed intake. The ReFI metric simplified modelling of feed utilisation efficiency in dairy cattle and resulted in breeding values that are equal to percentages of feed saved.

Impact of parity differences on residual feed intake estimation in Holstein cows

Residual feed intake (RFI) has been used as a measure of feed efficiency in farm animals. In lactating dairy cattle, RFI is typically obtained as the difference between dry matter intake observations and predictions from regression on known energy sinks, and effects of parity, days in milk, and cohort. The impact of parity (lactation number) on the estimation of RFI is not well understood, so the objectives of this study were to (1) evaluate alternative RFI models in which the energy sinks (metabolic body weight, body weight change, and secreted milk energy) were nested or not nested within parity, and (2) estimate variance components and genetic correlations for RFI across parities. Data consisted of 72,474 weekly RFI records of 5,813 lactating Holstein cows collected from 2007 to 2022 in 5 research stations across the United States. Estimates of heritability, repeatability, and genetic correlations between weekly RFI for parities 1, 2, and 3 were obtained using bivariate repeatability animal models. The nested RFI model showed better goodness of fit than the nonnested model, and some partial regression coefficients of dry matter intake on energy sinks were heterogeneous between parities. However, the Spearman's rank correlation between RFI values calculated from nested and nonnested models was equal to 0.99. Similarly, Spearman's rank correlation between the RFI breeding values from these 2 models was equal to 0.98. Heritability estimates for RFI were equal to 0.16 for parity 1, 0.19 for parity 2, and 0.22 for parity 3. Repeatability estimates for RFI across weeks within parities were high, ranging from 0.51 to 0.57. Spearman's rank correlations of sires' breeding values were 0.99 between parities 1 and 2, 0.91 between parities 1 and 3, and 0.92 between parities 2 and 3. We conclude that nesting energy sinks within parity when computing RFI improves model goodness of fit, but the impact on the estimated breading values appears to be minimal.

Animal board invited review: Genomic-based improvement of cattle in response to climate change

Climate change brings challenges to cattle production, such as the need to adapt to new climates and pressure to reduce greenhouse emissions (GHG). In general, the improvement of traits in current breeding goals is favourably correlated with the reduction of GHG. Current breeding goals and tools for increasing cattle production efficiency have reduced GHG. The same amount of production can be achieved by a much smaller number of animals. Genomic selection (GS) may offer a cost-effective way of using an efficient breeding approach, even in low- and middle-income countries. As climate change increases the intensity of heatwaves, adaptation to heat stress leads to lower efficiency of production and, thus, is unfavourable to the goal of reducing GHG. Furthermore, there is evidence that heat stress during cow pregnancy can have many generation-long lowering effects on milk production. Both adaptation and reduction of GHG are among the difficult-to-measure traits for which GS is more efficient and suitable than the traditional non-genomic breeding evaluation approach. Nevertheless, the commonly used within-breed selection may be insufficient to meet the new challenges; thus, cross-breeding based on selecting highly efficient and highly adaptive breeds may be needed. Genomic introgression offers an efficient approach for cross-breeding that is expected to provide high genetic progress with a low rate of inbreeding. However, well-adapted breeds may have a small number of animals, which is a source of concern from a genetic biodiversity point of view. Furthermore, low animal numbers also limit the efficiency of genomic introgression. Sustainable cattle production in countries that have already intensified production is likely to emphasise better health, reproduction, feed efficiency, heat stress and other adaptation traits instead of higher production. This may require the application of innovative technologies for phenotyping and further use of new big data techniques to extract information for breeding.

Impact of epistasis effects on the accuracy of predicting phenotypic values of residual feed intake in U. S Holstein cows

The impact of genomic epistasis effects on the accuracy of predicting the phenotypic values of residual feed intake (RFI) in U.S. Holstein cows was evaluated using 6215 Holstein cows and 78,964 SNPs. Two SNP models and seven epistasis models were initially evaluated. Heritability estimates and the accuracy of predicting the RFI phenotypic values from 10-fold cross-validation studies identified the model with SNP additive effects and additive × additive (A×A) epistasis effects (A + A×A model) to be the best prediction model. Under the A + A×A model, additive heritability was 0.141, and A×A heritability was 0.263 that consisted of 0.260 inter-chromosome A×A heritability and 0.003 intra-chromosome A×A heritability, showing that inter-chromosome A×A effects were responsible for the accuracy increases due to A×A. Under the SNP additive model (A-only model), the additive heritability was 0.171. In the 10 validation populations, the average accuracy for predicting the RFI phenotypic values was 0.246 (with range 0.197-0.333) under A + A×A model and was 0.231 (with range of 0.188-0.319) under the A-only model. The average increase in the accuracy of predicting the RFI phenotypic values by the A + A×A model over the A-only model was 6.49% (with range of 3.02-14.29%). Results in this study showed A×A epistasis effects had a positive impact on the accuracy of predicting the RFI phenotypic values when combined with additive effects in the prediction model.

Prediction of dry matter intake and gross feed efficiency using milk production and live weight in first-parity Holstein cows

Direct measurement of dry matter intake (DMI) presents a major challenge in estimating gross feed efficiency (GFE) in dairy cattle. This challenge can, however, be resolved through the prediction of DMI and GFE from easy-to-measure traits such as milk production (i.e. milk yield, energy-corrected milk (ECM), butterfat, protein, lactose) and live weight (LW). The main objective of this study was, therefore, to investigate the feasibility of predicting dry matter intake and gross feed efficiency for first-parity Holstein cows using milk production traits and LW. Data comprised of 30 daily measurements of DMI and milk production traits, and 25 daily LW records of a group of 100 first-parity Holstein cows, fed a total mixed ration. Gross feed efficiency was calculated as kg ECM divided by kg DMI. The initial step was to estimate correlations of milk production traits and LW with DMI and GFE, to identify the best potential predictors of DMI and GFE. Subsequently, a forward stepwise regression analysis was used to develop models to predict DMI and GFE from LW and milk production traits, followed by within-herd validations. Means for DMI, butterfat yield (BFY) and LW were 21.91 ± 2.77 kg/day, 0.95 ± 0.14 kg/day and 572 ± 15.58 kg/day, respectively. Mean GFE was 1.32 ± 0.22. Dry matter intake had positive correlations with milk yield (MY) (r = 0.32, p < 0.001) and LW (r = 0.76, p < 0.0001) and an antagonistic association with butterfat percent (BFP) (r =  - 0.55, p < 0.001). On the other hand, GFE was positively associated with MY (r = 0.36, p < 0.001), BFP (r = 0.53, p < 0.001) and BFY (r = 0.83, p < 0.0001), and negatively correlated with LW (r =  - 0.23, p > 0.05). Dry matter intake was predicted reliably by a model comprising of only LW and MY (R² = 0.79; root mean squared error (RMSE) = 1.05 kg/day). A model that included BFY, MY and LW had the highest ability to predict GFE (R² = 0.98; RMSE = 0.05). Live weight and BFY were the main predictor traits for DMI and GFE, respectively. The best models for predicting DMI and GFE were as follows: DMI (kg/day) =  - 54.21 - 0.192 × MY (kg/day) + 0.146 × LW (kg/day) and GFE (kg/day) = 4.120 + 0.024 × MY (kg/day) + 1.000 × BFY (kg/day) - 0.008 × LW (kg/day). Thus, daily DMI (kg/day) and GFE can be reliably predicted from LW and milk production traits using these developed models in first-parity Holstein cows. This presents a big promise to generate large quantities of data of individual cow DMI and GFE, which can be used to implement genetic improvement of feed efficiency.

Sharing of either phenotypes or genetic variants can increase the accuracy of genomic prediction of feed efficiency

Background

Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows.

Results

GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (r_g) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase.

Conclusions

The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00749-z.

Defining a nitrogen efficiency index in Holstein cows and assessing its potential effect on the breeding program of bulls

The purposes of this study were (1) to explore the relationship between 3 milk mid-infrared predicted features including nitrogen intake (NINT), milk true protein N (MTPN), and milk urea-N yield (MUNY); (2) to integrate these 3 features into an N efficiency index (NEI) and analyses approximate genetic correlations between the NEI and 37 traits (indices) of interest; and (3) to assess the potential effect of including the NEI into breeding programs of bulls. The edited data were 1,043,171 test-day records on 342,847 cows in 1,931 herds and 143,595 test-day records on 53,660 cows in 766 herds used for estimating breeding values (EBV) and variance components, respectively. The used records were within 5 to 50 d in milk. The records were grouped into primiparous and multiparous. The genetic parameters for the included mid-infrared features and EBV of the animals included in the pedigree were estimated using a multiple-trait repeatability animal model. Then, the EBV of the NINT, MTPN, MUNY were integrated into the NEI using a selection index assuming weights based on the N partitioning. The approximate genetic correlations between the NEI and 37 traits of interest were estimated using the EBV of the selected bulls. The bulls born from 2011 to 2014 with NEI were selected and the NEI distribution of these bulls having EBV for the 8 selected traits (indices) was checked. The heritability and repeatability estimates for NINT, MTPN, and MUNY ranged from 0.09 to 0.13, and 0.37 to 0.65, respectively. The genetic and phenotypic correlations between NINT, MTPN, and MUNY ranged from -0.31 to 0.87, and -0.02 to 0.42, respectively. The NEI ranged from -13.13 to 12.55 kg/d. In total, 736 bulls with reliability ≥0.50 for all included traits (NEI and 37 traits) and at least 10 daughters distributed in at least 10 herds were selected to investigate genetic aspects of the NEI. The NEI had positive genetic correlations with production yield traits (0.08-0.46), and negative genetic correlations with the investigated functional traits and indices (-0.71 to -0.07), except for the production economic index and functional type economic index. The daughters of bulls with higher NEI had lower NINT and MUNY, and higher MTPN. Furthermore, 26% of the bulls (n = 50) with NEI born between 2011 to 2014 had higher NEI and global economic index than the average in the selected bulls. Finally, the developed NEI has the advantage of large-scale prediction and therefore has the potential for routine application in dairy cattle breeding in the future.

Evaluation of feed efficiency traits in different Hereford populations and their effect on variance component estimation

Context Residual feed intake is a relevant trait for beef cattle, given the positive impact on reducing feeding costs and greenhouse gas emissions. The lack of large databases is a restriction when estimating accurate genetic parameters for dry matter intake (DMI) and residual feed intake (RFI), and combining different data sets could be an alternative to increase the amount of data and achieve better estimations. Aim The main objective was to compare Uruguayan data (URY; 780 bulls) and Canadian data (CAN; 1597 bulls), and to assess the adequacy of pooling both data sets (ALL) for the estimation of genetic parameters for DMI and RFI. Methods Feed intake and growth traits phenotypes in both data sets were measured following the same protocols established by the Beef Improvement Federation. Pedigree connections among data sets existed, but were weak. Performance data were analysed for each data set, and individual partial regression coefficients for each energy sink on DMI were obtained and compared. Univariate and multivariate variance components were estimated by the restricted maximum likelihood (REML) for DMI, RFI and their energy sinks traits (average daily gain, metabolic mid weight and back fat thickness). Key results There were some differences in phenotypic performance among data (P < 0.01); however, no differences (P > 0.1) were observed for phenotypic values of RFI between sets. Heritability estimates for DMI were 0.42 (URY), 0.41 (CAN) and 0.45 for ALL data, whereas heritability estimates for RFI were 0.34 (URY), 0.20 (CAN) and 0.25 for ALL data. The results obtained indicate selection on reducing RFI could lead to a decrease in DMI, without compromising other performance traits, as genetic correlations between RFI, growth and liveweight were low or close to 0 (−0.12–0.07). Conclusions As genetic parameters were similar between national data sets (URY, CAN), pooling data (ALL) provided more accurate parameter estimations, as they presented smaller standard deviations, especially in multivariate analysis. Implications Parameters estimated here may be used in international or national genetic evaluation programs.

Estimates of genetic parameters for feeding behavior traits and their associations with feed efficiency in Holstein cows

Residual feed intake (RFI) is commonly used to measure feed efficiency but individual intake recording systems are needed. Feeding behavior may be used as an indicator trait for feed efficiency using less expensive precision livestock farming technologies. Our goal was to estimate genetic parameters for feeding behavior and the genetic correlations with feed efficiency in Holstein cows. Data consisted of 75,877 daily feeding behavior records of 1,328 mid-lactation Holstein cows in 31 experiments conducted from 2009 to 2020 with an automated intake recording system. Feeding behavior traits included number of feeder visits per day, number of meals per day, duration of each feeder visit, duration of each meal, total duration of feeder visits, intake per visit, intake per meal [kg of dry matter (DM)], feeding rate per visit, and feeding rate per meal (kg of DM per min). The meal criterion was estimated as 26.4 min, which means that any pair of feeder visits separated by less than 26.4 min were considered part of the same meal. The statistical model included lactation and days in milk as fixed effects, and experiment-treatment, animal, and permanent environment as random effects. Genetic parameters for feeding behavior traits were estimated using daily records and weekly averages. Estimates of heritability for daily feeding behavior traits ranged from 0.09 ± 0.02 (number of meals; mean ± standard error) to 0.23 ± 0.03 (feeding rate per meal), with repeatability estimates ranging from 0.23 ± 0.01 (number of meals) to 0.52 ± 0.02 (number of feeder visits). Estimates of heritability for weekly averages of feeding behavior traits ranged from 0.19 ± 0.04 (number of meals) to 0.32 ± 0.04 (feeding rate per visit), with repeatability estimates ranging from 0.46 ± 0.02 (duration of each meal) to 0.62 ± 0.02 (feeding rate per visit and per meal). Most of the feeding behavior measures were strongly genetically correlated, showing that with more visits or meals per day, cows spend less time in each feeder visit or meal with lower intake per visit or meal. Weekly averages for feeding behavior traits were analyzed jointly with RFI and its components. Number of meals was genetically correlated with milk energy (0.48), metabolic body weight (-0.27), and RFI (0.19). Duration of each feeder visit and meal were genetically correlated with milk energy (0.43 and 0.44, respectively). Total duration of feeder visits per day was genetically correlated with DM intake (0.29), milk energy (0.62), metabolic body weight (-0.37), and RFI (0.20). Intake per visit and meal were genetically correlated with DM intake (0.63 and 0.87), milk energy (0.47 and 0.69), metabolic body weight (0.47 and 0.68), and RFI (0.31 and 0.65). Feeding rate was genetically correlated with DM intake (0.69), metabolic body weight (0.67), RFI (0.47), and milk energy (0.21). We conclude that measures of feeding behavior could be useful indicators of dairy cow feed efficiency, and individual cows that eat at a slower rate may be more feed efficient.

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.

Associations of Wearable Sensor Measures With Feed Intake, Production Traits, Lactation, and Environmental Parameters Impacting Feed Efficiency in Dairy Cattle

Feed efficiency is an important trait to dairy production because of its impact on sustainability and profitability. Measuring individual cow feed intake on commercial farms would be unfeasibly costly at present. Thus, developing cheap and portable indicators of feed intake would be highly beneficial for genetic selection and precision feeding management tools. Given the growing use of automated sensors on dairy farms, the objective of this study was to determine the relationship between measurements recorded from multiple wearable sensors and feed intake. A total of three different wearable sensors were evaluated for their association with dry mater intake (DMI). The sensors measured activity (sensors = 3), rumination (sensors = 1), ear temperature (sensors = 1), rumen pH (sensors = 1) and rumen temperature (sensors = 1). A range of 56–340 cows with assorted sensors from 24 to 313 days in milk (DIM) were modeled to evaluate associations with DIM, parity, and contemporary group (CG; comprised of pen and study cohort). Models extending upon these variables included known energy sinks (i.e., milk production, milk fat/protein and metabolic body weight), to characterize the association of sensors measures and DMI. Statistically significant (i.e., P &lt; 0.05) regression coefficients for individual sensor measures with DMI ranged from 9.01E-07 to −3.45 kg DMI/day. When integrating all measures from a single sensor in a model, estimated regression coefficients ranged 8.83E-07 to −3.48 kg DMI/day. Significant associations were also identified for milk production traits, parity, DIM and CG. Associations tended to be highest for timepoints around the time of feeding and when multiple measurements within a sensor were integrated in a single model. The findings of this study indicate sensor measures are associated with feed intake and other energy sink traits and variables impacting feed efficiency. This information would be helpful to improve feed and feeding efficiency on commercial farms as proxy measurements for feed intake.

Effect of diet energy level and genomic residual feed intake on bred Holstein dairy heifer growth and feed efficiency

The objective of this study was to determine growth, feed intake, and feed efficiency of postbred dairy heifers with different genomic residual feed intake (RFI) predicted as a lactating cow when offered diets differing in energy density. Postbred Holstein heifers (n = 128, ages 14-20 mo) were blocked by initial weight (high, medium-high, medium-low, and low) with 32 heifers per block. Each weight block was sorted by RFI (high or low) to obtain 2 pens of heifers with high and low genomically predicted RFI within each block (8 heifers per pen). Low RFI heifers were expected to have greater feed efficiency than high RFI heifers. Dietary treatments consisted of a higher energy control diet based on corn silage and alfalfa haylage [HE; 62.7% total digestible nutrients, 11.8% crude protein, and 45.6% neutral detergent fiber; dry matter (DM) basis], and a lower energy diet diluted with straw (LE; 57.0% total digestible nutrients, 11.7% crude protein, and 50.1% neutral detergent fiber; DM basis). Each pen within a block was randomly allocated a diet treatment to obtain a 2 × 2 factorial arrangement (2 RFI levels and 2 dietary energy levels). Diets were offered in a 120-d trial. Dry matter intake by heifers was affected by diet (11.0 vs. 10.0 kg/d for HE and LE, respectively) but not by RFI or the interaction of RFI and diet. Daily gain was affected by the interaction of RFI and diet, with low RFI heifers gaining more than high RFI heifers when fed LE (0.94 vs. 0.85 kg/d for low and high RFI, respectively), but no difference for RFI groups when fed HE (1.16 vs. 1.19 kg/d for low and high RFI, respectively). Respective feed efficiencies were improved for low RFI compared with high RFI heifers when fed LE (10.6 vs. 11.8 kg of feed DM/kg of gain), but no effect of RFI was found when fed HE (9.4 vs. 9.5 kg of DM/kg of gain for high and low RFI, respectively). No effect of RFI or diet on first-lactation performance through 150 DIM was observed. Based on these results, the feed efficiency of heifers having different genomic RFI may be dependent on diet energy level, whereby low RFI heifers utilized the LE diet more efficiently. The higher fiber straw (LE) diet controlled intake and maintained more desirable heifer weight gains. This suggests that selection for improved RFI in lactating cows may improve feed efficiency in growing heifers when fed to meet growth goals of 0.9 to 1.0 kg of gain/d.

Opportunities to Harness High-Throughput and Novel Sensing Phenotypes to Improve Feed Efficiency in Dairy Cattle

Simple Summary

Sensors, routinely collected on-farm tests, and other repeatable, high-throughput measurements can provide novel phenotype information on a frequent basis. Information from these sensors and high-throughput measurements could be harnessed to monitor or predict individual dairy cow feed intake. Predictive algorithms would allow for genetic selection of animals that consume less feed while producing the same amount of milk. Improved monitoring of feed intake could reduce the cost of milk production, improve animal health, and reduce the environmental impact of the dairy industry. Moreover, data from these information sources could aid in animal management (e.g., precision feeding and health detection). In order to implement tools, the relationship of measurements with feed intake needs to be established and prediction equations developed. Lastly, consideration should be given to the frequency of data collection, the need for standardization of data and other potential limitations of tools in the prediction of feed intake. This review summarizes measurements of feed efficiency, factors that may impact the efficiency and feed consumption of an animal, tools that have been researched and new traits that could be utilized for the prediction of feed intake and efficiency, and prediction equations for feed intake and efficiency presented in the literature to date.

Abstract

Feed for dairy cattle has a major impact on profitability and the environmental impact of farms. Sustainable dairy production relies on continued improvement in feed efficiency as a way to reduce costs and nutrient loss from feed. Advances in breeding, feeding and management have led to the dilution of maintenance energy and thus more efficient dairy cattle. Still, many additional opportunities are available to improve individual animal feed efficiency. Sensing technologies such as wearable sensors, image-based and high-throughput phenotyping technologies (e.g., milk testing) are becoming more available on commercial farm. The application of these technologies as indicator traits for feed intake and efficiency related traits would be advantageous to provide additional information to predict and manage feed efficiency. This review focuses on precision livestock technologies and high-throughput phenotyping in use today as well as those that could be developed in the future as possible indicators of feed intake. Several technologies such as milk spectral data, activity, rumen measures, and image-based phenotypes have been associated with feed intake. Future applications will depend on the ability to repeatably measure and calibrate these data across locations, so that they can be integrated for use in predicting and managing feed intake and efficiency on farm.

Circulating Metabolites Indicate Differences in High and Low Residual Feed Intake Holstein Dairy Cows

Selection for more feed efficient dairy cows is key to improving sustainability and profitability of dairy production; however, underlying mechanisms contributing to individual animal feed efficiency are not fully understood. The objective of this study was to identify circulating metabolites, and pathways associated with those metabolites, that differ between efficient and inefficient Holstein dairy cows using targeted metabolite quantification and untargeted metabolomics. The top and bottom fifteen percent of cows (n = 28/group) with the lowest and highest residual feed intake in mid-lactation feed efficiency trials were grouped retrospectively as high-efficient (HE) and low-efficient (LE). Blood samples were collected for quantification of energy metabolites, markers of hepatic function, and acylcarnitines, in addition to a broader investigation using untargeted metabolomics. Short-chain acylcarnitines, C3-acylcarnitine, and C4-acylcarntine were lower in HE cows (n = 18/group). Untargeted metabolomics and multivariate analysis identified thirty-nine differential metabolites between HE and LE (n = 8/group), of which twenty-five were lower and fourteen were higher in HE. Pathway enrichment analysis indicated differences in tryptophan metabolism. Combined results from targeted metabolite quantification and untargeted metabolomics indicate differences in fatty acid and amino acid metabolism between HE and LE cows. These differences may indicate post-absorptive nutrient use efficiency as a contributor to individual animal variation in feed efficiency.

Combining datasets in a dynamic residual feed intake model and comparison with linear model results in lactating Holstein cattle

A new method to estimate residual feed intake (RFI) was recently developed based on a multi-trait random regression model. This approach deals with the dynamic nature of the lactation, which is in contrast with classical linear approaches. However, an issue remains: pooling data across sites and years, which implies dealing with different (and sometimes unknown) diet energy contents. This will be needed for genomic evaluation. In this study, we tested whether merging two individual datasets into a larger one can lead to valuable results in comparison to analysing them on their own with the multi-trait random regression model. Three datasets were defined: the first one with 1 063 lactations, the second one with 205 lactations from a second farm and the third one combining the data of the two first datasets (1 268 lactations). The model was applied to the three datasets to estimate individual RFI as well as variance components and correlations between the four traits included in the model (fat and protein corrected milk production, BW, feed intake and body condition score), and a fixed month-year-farm effect was used to define the contemporary group. The variance components and correlations between animal effects of the four traits were very similar irrespective of the dataset used with correlations higher than 0.94 between the different datasets. The RFI estimates for animals from their single farm only were also very similar (r > 0.95) to the ones computed from the merged dataset (Dataset 3). This highlights that the contemporary group correction in the model adequately accounts for differences between the two feeding environments. The dynamic model can thus be used to produce RFI estimates from merged datasets, at least when animals are raised in similar systems. In addition, the 205 lactations from the second farm were also used to estimate the RFI with a linear approach. The RFI estimated by the two approaches were similar when the considered period was rather short (r = 0.85 for RFI for the first 84 days of lactation) but this correlation weakened as the period length grew (r = 0.77 for RFI for the first 168 days of lactation). This weakening in correlations between the two approaches when increasing the used time-period reflects that only the dynamic model permits the regression coefficients to evolve in line with the physiological changes through the lactation. The results of this study enlarge the possibilities of use for the dynamic RFI model.

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.

An alternative interpretation of residual feed intake by phenotypic recursive relationships in dairy cattle

There has been increasing interest in residual feed intake (RFI) as a measure of net feed efficiency in dairy cattle. Residual feed intake phenotypes are obtained as residuals from linear regression encompassing relevant factors (i.e., energy sinks) to account for body tissue mobilization. By rearranging the single-trait linear regression, we showed a causal RFI interpretation underlying the linear regression for RFI. It postulates recursive effects in energy allocation from energy sinks on dry matter intake, but the feedback or simultaneous effects are nonexistent. A Bayesian recursive structural equation model was proposed for directly predicting RFI and energy sinks and estimating relevant genetic parameters simultaneously. A simplified Markov chain Monte Carlo algorithm was described. The recursive model is asymptotically equivalent to one-step linear regression for RFI, yet extends the analytical capacity to multiple-trait analysis.

Residual energy intake, energy balance, and liability to diseases: Genetic parameters and relationships in German Holstein dairy cows

Residual energy intake (REI) is an often-suggested trait for direct selection of dairy cows for feed efficiency. Cows with lower REI seem to be more efficient but are also in a more severe negative energy balance (EB), especially in early lactation. A negative EB leads to a higher liability to diseases. Due to this fact, this study aims to investigate the genetic relationship between REI and liability to diseases. Health and production data were recorded from 1,370 German Holstein dairy cows from 8 research farms over a period of 2 yr. We calculated 2 phenotypes for REI that considered the following energy sinks: milk energy content, metabolic body weight, body weight change, body condition score, and body condition score change. Genetic parameters were estimated with threshold or linear random regression models from days in milk (DIM) 1 to 305. Heritabilities for REI, EB, and all diseases ranged from 0.12 to 0.39, 0.15 to 0.31, and 0.09 to 0.20, respectively. Genetic correlations between selected DIM for REI and EB were higher for adjacent DIM than for more distant DIM. Pearson correlation coefficients between estimated breeding values (EBV) for REI and EB varied between 0.47 and 0.81; they were highest in mid lactation. Correlations between EBV for all diseases and REI as well as EB were negative, with lowest values in early lactation. Within the first 50 DIM, proportions of diseased days for cows with lowest EBV for REI were almost twice as high as for cows with highest EBV for REI. In conclusion, selecting dairy cows for lower REI should be treated with caution because of an unfavorable relationship with liability to diseases, especially in early lactation.

Importance of considering body weight change in response to dietary protein deficiency in lactating dairy cows

Our objective was to quantify the contribution of body weight (BW) change to the overall response of lactating dairy cows to a shortage of dietary protein. Lactating Holstein cows (n = 166; 92 primiparous, 74 multiparous) with initial milk yield of 41 ± 10 kg/d were fed high and low-protein diets in 7 blocks. Blocks were repeated in the same crossover design with periods of 28 to 35 d. Production of 69 of the 166 cows (42 primiparous, 27 multiparous) was also measured in late lactation. Low-protein diets were 14% crude protein (CP) in peak lactation and 13% CP in late lactation and were formulated to contain adequate rumen-degradable protein to maintain rumen function but inadequate rumen undegradable protein for the average cow in this study. High-protein diets were 18% CP in peak lactation and 16% CP in late lactation and contained extra expeller soybean meal to meet metabolizable protein requirements. Body weight changes were used to predict body energy and protein changes, which were added to milk components to calculate total captured energy and protein. Fixed effects of diet, parity, treatment sequence nested in each block, treatment period nested in block, interaction of diet and parity, and the random effects of block and cow nested within block were included in the model to compare cow responses to diets within each lactation stage. In peak lactation, reducing protein from 18 to 14% resulted in estimated daily losses of 2.9 Mcal of milk energy, 2.2 Mcal of body energy, 127 g of milk protein, and 16 g of body protein. Therefore, BW loss accounted for 43% of the decrease in captured energy and 11% of the decrease in captured protein when cows were fed deficient protein. In late lactation, BW loss accounted for 51% of the decrease in captured energy and 14% of the decrease in captured protein when cows were fed deficient protein. We suggest that BW change should be considered when assessing cow responses to changes in dietary protein.

Breeding goal traits accounting for feed intake capacity and roughage or concentrate intake separately

Current breeding tools aiming to improve feed efficiency use definitions based on total dry matter intake (DMI); for example, residual feed intake or feed saved. This research aimed to define alternative traits using existing data that differentiate between feed intake capacity and roughage or concentrate intake, and to investigate the phenotypic and genetic relationships among these traits. The data set contained 39,017 weekly milk yield, live weight, and DMI records of 3,164 cows. The 4 defined traits were as follows: (1) Feed intake capacity (FIC), defined as the difference between how much a cow ate and how much she was expected to eat based on diet satiety value and status of the cow (parity and lactation stage); (2) feed saved (FS), defined as the difference between the measured and the predicted DMI, based on the regression of DMI on milk components within experiment; (3) residual roughage intake (RRI), defined as the difference between the measured and the predicted roughage intake, based on the regression of roughage intake on milk components and concentrate intake within experiment; and (4) residual concentrate intake (RCI), defined as the difference between the measured and the predicted concentrate intake, based on the regression of concentrate intake on milk components and roughage intake within experiment. The phenotypic correlations were -0.72 between FIC and FS, -0.84 between FS and RRI, and -0.53 between FS and RCI. Heritability of FIC, FS, RRI, and RCI were estimated to be 0.21, 0.12, 0.15, and 0.03, respectively. The genetic correlations were -0.81 between FS and FIC, -0.96 between FS and RRI, and -0.25 between FS and RCI. Concentrate intake and RCI had low heritability. Genetic correlation between DMI and FIC was 0.98. Although the defined traits had moderate phenotypic correlations, the genetic correlations between DMI, FS, FIC, and RRI were above 0.79 (in absolute terms), suggesting that these traits are genetically similar. Therefore, selecting for FIC is expected to simply increase DMI and RRI, and there seems to be little advantage in separating concentrate and roughage intake in the genetic evaluation, because measured concentrate intake was determined by the feeding system in our data and not by the genetics of the cow.

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.

Assessment of the Relationship between Postpartum Health and Mid-Lactation Performance, Behavior, and Feed Efficiency in Holstein Dairy Cows

The objective of this study was to investigate the relationships between postpartum health disorders and mid-lactation performance, feed efficiency, and sensor-derived behavioral traits. Multiparous cows (n = 179) were monitored for health disorders for 21 days postpartum and enrolled in a 45-day trial between 50 to 200 days in milk, wherein feed intake, milk yield and components, body weight, body condition score, and activity, lying, and feeding behaviors were recorded. Feed efficiency was measured as residual feed intake and the ratio of fat- or energy-corrected milk to dry matter intake. Cows were classified as either having hyperketonemia (HYK; n = 72) or not (n = 107) and grouped by frequency of postpartum health disorders: none (HLT; n = 94), one (DIS; n = 63), or ≥2 (DIS+; n = 22). Cows that were diagnosed with HYK had higher mid-lactation yields of fat- and energy-corrected milk. No differences in feed efficiency were detected between HYK or health status groups. Highly active mid-lactation time was higher in healthy animals, and rumination time was lower in ≥4th lactation cows compared with HYK or DIS and DIS+ cows. Differences in mid-lactation behaviors between HYK and health status groups may reflect the long-term impacts of health disorders. The lack of a relationship between postpartum health and mid-lactation feed efficiency indicates that health disorders do not have long-lasting impacts on feed efficiency.

Comparison of methods to predict feed intake and residual feed intake using behavioral and metabolite data in addition to classical performance variables

Predicting dry matter intake (DMI) and feed efficiency by leveraging the use of data streams available on farm could aid efforts to improve the feed efficiency of dairy cattle. Residual feed intake (RFI) is the difference between predicted and observed feed intake after accounting for body size, body weight change, and milk production, making it a valuable metric for feed efficiency research. Our objective was to develop and evaluate DMI and RFI prediction models using multiple linear regression (MLR), partial least squares regression, artificial neural networks, and stacked ensembles using different combinations of cow descriptive, performance, sensor-derived behavioral (SMARTBOW; Zoetis), and blood metabolite data. Data were collected from mid-lactation Holstein cows (n = 124; 102 multiparous, 22 primiparous) split equally between 2 replicates of 45-d duration with ad libitum access to feed. Within each predictive approach, 4 data streams were added in sequence: dataset M (week of lactation, parity, milk yield, and milk components), dataset MB (dataset M plus body condition score and metabolic body weight), dataset MBS (dataset MB plus sensor-derived behavioral variables), and dataset MBSP (dataset MBS plus physiological blood metabolites). The combination of 4 datasets and 4 analytical approaches resulted in 16 analyses of DMI and RFI, using variables averaged within cow across the study period. Additional models using weekly averaged data within cow and study were built using all predictive approaches for datasets M, MB, and MBS. Model performance was assessed using the coefficient of determination, concordance correlation coefficient, and root mean square error of prediction. Predictive models of DMI performed similarly across all approaches, and models using dataset MBS had the greatest model performance. The best approach-dataset combination was MLR-dataset MBS, although several models performed similarly. Weekly DMI models had the greatest performance with MLR and partial least squares regression approaches. Dataset MBS models had incrementally better performance than datasets MB and M. Within each approach-dataset combination, models with DMI averaged over the study period had slightly greater model performance than DMI averaged weekly. Predictive performance of all RFI models was poor, but slight improvements when using MLR applied to dataset MBS suggest that rumination and activity behaviors may explain some of the variation in RFI. Overall, similar performance of MLR, compared with machine learning techniques, indicates MLR may be sufficient to predict DMI. The improvement in model performance with each additional data stream supports the idea of integrating data streams to improve model predictions and farm management decisions.

Effects of Incorporating Dry Matter Intake and Residual Feed Intake into a Selection Index for Dairy Cattle Using Deterministic Modeling

The inclusion of feed efficiency in the breeding goal for dairy cattle has been discussed for many years. The effects of incorporating feed efficiency into a selection index were assessed by indirect selection (dry matter intake) and direct selection (residual feed intake) using deterministic modeling. Both traits were investigated in three ways: (1) restricting the trait genetic gain to zero, (2) applying negative selection pressure, and (3) applying positive selection pressure. Changes in response to selection from economic and genetic gain perspectives were used to evaluate the impact of including feed efficiency with direct or indirect selection in an index. Improving feed efficiency through direct selection on residual feed intake was the best scenario analyzed, with the highest overall economic response including favorable responses to selection for production and feed efficiency. Over time, the response to selection is cumulative, with the potential for animals to reduce consumption by 0.16 kg to 2.7 kg of dry matter per day while maintaining production. As the selection pressure increased on residual feed intake, the response to selection for production, health, and fertility traits and body condition score became increasingly less favorable. This work provides insight into the potential long-term effects of selecting for feed efficiency as residual feed intake.

Invited review: Disentangling residual feed intake-Insights and approaches to make it more fit for purpose in the modern context

Residual feed intake (RFI) is an increasingly used trait to analyze feed efficiency in livestock, and in some sectors such as dairy cattle, it is one of the most frequently used traits. Although the principle for calculating RFI is always the same (i.e., using the residual of a regression of intake on performance predictors), a wide range of models are found in the literature, with different predictors, different ways of considering intake, and more recently, different statistical approaches. Consequently, the results are not easily comparable from one study to another as they reflect different biological variabilities, and the relationship between the residual (i.e., RFI) and the underlying true efficiency also differs. In this review, the components of the RFI equation are explored with respect to the underlying biological processes. The aim of this decomposition is to provide a better understanding of which of the processes in this complex trait contribute significantly to the individual variability in efficiency. The intricacies associated with the residual term, as well as the energy sinks and the intake term, are broken down and discussed. Based on this exploration as well as on some recent literature, new forms of the RFI equation are proposed to better separate the efficiency terms from errors and inaccuracies. The review also considers the time period of measurement of RFI. This is a key consideration for the accuracy of the RFI estimation itself, and also for understanding the relationships between short-term efficiency, animal resilience, and long-term efficiency. As livestock production moves toward sustainable efficiency, these considerations are increasingly important to bring to bear in RFI estimations.

A new method to estimate residual feed intake in dairy cattle using time series data

In dairy, the usual way to measure feed efficiency is through the residual feed intake (RFI) method. However, this method is, in its classical form, a linear regression, which, by construction, does not take into account the evolution of the RFI components across time, inducing approximations in the results. We present here a new approach that incorporates the dynamic dimension of the data. Using a multitrait random regression model, the correlations between milk, live weight, DM intake (DMI) and body condition score (BCS) were investigated across the lactation. In addition, at each time point, by a matrix regression on the variance-covariance matrix and on the animal effects from the three predictor traits, a predicted animal effect for intake was estimated, which, by difference with the actual animal effect for intake, gave a RFI estimation. This model was tested on historical data from the Aarhus University experimental farm (1 469 lactations out of 740 cows). Correlations between animal effects were positive and high for milk and DMI and for weight and DMI, with a maximum mid-lactation, stable across time at around 0.4 for weight and BCS, and slowly decreasing along the lactation for milk and weight, DMI and BCS, and milk and BCS. At the Legendre polynomial coefficient scale, the correlations were estimated with a high accuracy (averaged SE of 0.04, min = 0.02, max = 0.05). The predicted animal effect for intake was always extremely highly correlated with the milk production and highly correlated with BW for the most part of the lactation, but only slightly correlated with BCS, with the correlation becoming negative in the second half of the lactation. The estimated RFI possessed all the characteristics of a classical RFI, with a mean at zero at each time point and a phenotypic independence from its predictors. The correlation between the averaged RFI over the lactation and RFI at each time point was always positive and above 0.5, and maximum mid-lactation (>0.9). The model performed reasonably well in the presence of missing data. This approach allows a dynamic estimation of the traits, free from all time-related issues inherent to the traditional RFI methodology, and can easily be adapted and used in a genetic or genomic selection context.

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.

Residual Feed Intake in Dairy Ewes: An Evidence of Intraflock Variability

Simple Summary

Few, if any, reference is available in residual feed intake in dairy sheep. In this study, carried out during more than two months with French Lacaune dairy ewes in mid-lactation, we demonstrated an intraflock variability in feed efficiency determined by, beyond litter size and daily milking frequency, evident differences between the individuals in their efficiency of using the available total mixed rations.

Abstract

This study examined the intraflock variability of feed efficiency in dairy ewes, through monitoring residual feed intakes (RFI). Primiparous lactating ewes (n = 43; 57.7 ± 0.91 kg body weight [BW] at lambing), representative of a French Lacaune dairy flock, were allocated in an equilibrated 2 × 2 factorial design experiment, lasting for 63 days during mid-lactation and combining 2 litter sizes (singletons, SING or twins, TWIN) and 2 daily milking frequencies (once, ONE or twice, TWO). Weaning occurred, and milking started, at 35 days after lambing (DIM). Ewes were individually fed a diet based on ryegrass silage, local hay, and supplements. Individual dry matter intake (DMI) was recorded daily and further used to evaluate (and compare) differences in RFI between ewes at 42, 49, 56, 63, 70, 77, 84, 91, and 98. Average individual RFI were calculated weekly since the first week (i.e., 35–42 DIM). Total (BW) and metabolic (BW^0.75) body weight, body condition score BCS, milk yield, and plasma non-esterified fatty acids NEFA were monitored weekly. Differences in DMI were mainly due to the lactation stage and litter size and were 11% higher in ewes with TWIN compared to SING. This was positively correlated to milk yield and consistent with differences in RFI which varied due to litter size and to the milking frequency × lactation stage interaction. Ewes that lambed SING showed higher feed efficiency (−0.08 ± 0.018 vs. 0.13 ± 0.014 kg DM/ewe/d of RFI in SING vs. TWIN, respectively), whereas there were no differences in BW or BCS. Milking frequency did not affect DMI but milk yields were higher in TWO, which was related to a higher feed efficiency in this group (0.115 ± 0.016 vs. −0.07 ± 0.016 kg DM/ewe/d of RFI in ONE vs. TWO, respectively). Average RFI was affected (p < 0.0001) by the ewe, thus allowing a ranking among individuals to be established. High (n = 22) or low (n = 21) feed efficiency ewes averaged −0.17 ± 0.09 or 0.18 ± 0.09 kg DM/d RFI, respectively. Estimates of RFI were not correlated to the individual milk production potential. Even if no differences in BW, BW^0.75, or BCS were detected, high-efficiency ewes mobilized 1.5 times their body reserves (0.30 vs. 0.20 mmol NEFA/L of plasma) when compared to the low-efficiency group. The observed intraflock variability in feed efficiency of this dairy ewes’ flock was affected by litter size and milking frequency but also by evident differences between individuals’ physiologies.