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.

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.

Relationship between feed efficiency and resilience in dairy ewes subjected to acute underfeeding

Selection of dairy sheep based on production levels has caused a loss of rusticity, which might compromise their future resilience to nutritional challenges. Although refocusing breeding programs toward improved feed efficiency (FE) is expected, more-efficient ewes also seem to be more productive. As a first step to examine the relationship between FE and resilience in dairy sheep, in this study we explored the variation in the response to and the recovery from an acute nutritional challenge in high-yielding Assaf ewes phenotypically divergent for FE. First, feed intake, milk yield and composition, and body weight changes were recorded individually over a 3-wk period in a total of 40 sheep fed a total mixed ration (TMR) ad libitum. Data were used to calculate their FE index (FEI, defined as the difference between the actual and predicted intake estimated through net energy requirements for maintenance, production, and weight change). The highest and lowest FE ewes (H-FE and L-FE groups, respectively; 10 animals/group) were selected and then subjected to the nutritional challenge (i.e., withdrawing the TMR and limiting their diet only to the consumption of straw for 3 d). Afterward, sheep were fed again the TMR ad libitum. Temporal patterns of variation in performance traits, and ruminal fermentation and blood parameters were examined. A good consistency between FEI, residual feed intake, and feed conversion ratio was observed. Results supported that H-FE were more productive than L-FE sheep at similar intake level. Average time trends of milk yield generated by a piecewise model suggest that temporal patterns of variation in this trait would be related to prechallenge production level (i.e., H-FE presented quicker response and recovery than L-FE). Considering all studied traits, the overall response to and recovery from underfeeding was apparently similar or even better in H-FE than in L-FE. This would refute the initial hypothesis of a poorer resilience of more-efficient sheep to an acute underfeeding. However, the question remains whether a longer term feed restriction might impair the ability of H-FE ewes to maintain or revert to a high-production status, which would require further research.

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

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

Implementation of Feed Efficiency in Iranian Holstein Breeding Program

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

An evaluation of detailed animal characteristics influencing the lactation production efficiency of spring-calving, pasture-based dairy cattle

Selection for feed efficiency, the ratio of output (e.g., milk yield) to feed intake, has traditionally been limited on commercial dairy farms by the necessity for detailed individual animal intake and performance data within large animal populations. The objective of the experiment was to evaluate the effects of individual animal characteristics (animal breed, genetic potential, milk production, body weight (BW), daily total dry matter intake (TDMI), and energy balance) on a cost-effective production efficiency parameter calculated as the annual fat and protein (milk solids) production per unit of mid-lactation BW (MSperBW_lact). A total of 1,788 individual animal intake records measured at various stages of lactation (early, mid, and late lactation) from 207 Holstein-Friesian and 200 Jersey × Holstein-Friesian cows were used. The derived efficiency traits included daily kilograms of milk solids produced per 100 kg of BW (dMSperBW_int) and daily kilograms of milk solids produced per kilogram of TDMI (dMSperTDMI). The TDMI per 100 kg of BW was also calculated (TDMI/BW_int) at each stage of lactation. Animals were subsequently either ranked as the top 25% (Heff) or bottom 25% (Leff) based on their lactation production efficiency (MSperBW_lact). Dairy cow breed significantly affected animal characteristics over the entire lactation and during specific periods of intake measurements. Jersey crossbred animals produced more milk, based on a lower TDMI, and achieved an increased intake per kilogram of BW. Similarly, Heff produced more milk over longer lactations, weighed less, were older, and achieved a higher TDMI compared with the Leff animals. Both Jersey × Holstein-Friesian and Heff cows achieved superior production efficiency due to lower maintenance energy requirements, and consequentially increased milk solids production per kilogram of BW and per kilogram of TDMI at all stages of lactation. Indeed, within breed, Heff animals weighed 20 kg less and produced 15% more milk solids over the total lactation than Leff. In addition, Heff achieved increased daily milk solids yield (+0.16 kg) and milk solids yield per kilogram of TDMI (+ 0.23 kg/kg DM) during intake measurement periods. Moreover, the strong and consistently positive correlations between MSperBW_lact and detailed production efficiency traits (dMSperBW_int, dMSperTDMI) reported here demonstrate that MSperBW_lact is a robust measure that can be applied within commercial grazing dairy systems to increase the selection intensity for highly efficient animals.

Genetic analysis of production traits and body size measurements and their relationships with metabolic diseases in German Holstein cattle

This study sheds light on the genetic complexity and interplay of production, body size, and metabolic health in dairy cattle. Phenotypes for body size-related traits from conformation classification (130,166 animals) and production (101,562 animals) of primiparous German Holstein cows were available. Additionally, 21,992, 16,641, and 7,096 animals were from herds with recordings of the metabolic diseases ketosis, displaced abomasum, and milk fever in first, second, and third lactation. Moreover, all animals were genotyped. Heritabilities of traits and genetic correlations between all traits were estimated and GWAS were performed. Heritability was between 0.240 and 0.333 for production and between 0.149 and 0.368 for body size traits. Metabolic diseases were lowly heritable, with estimates ranging from 0.011 to 0.029 in primiparous cows, from 0.008 to 0.031 in second lactation, and from 0.037 to 0.052 in third lactation. Production was found to have negative genetic correlations with body condition score (BCS; -0.279 to -0.343) and udder depth (-0.348 to -0.419). Positive correlations were observed for production and body depth (0.138-0.228), dairy character (DCH) (0.334-0.422), and stature (STAT) (0.084-0.158). In first parity cows, metabolic disease traits were unfavorably correlated with production, with genetic correlations varying from 0.111 to 0.224, implying that higher yielding cows have more metabolic problems. Genetic correlations of disease traits in second and third lactation with production in primiparous cows were low to moderate and in most cases unfavorable. While BCS was negatively correlated with metabolic diseases (-0.255 to -0.470), positive correlations were found between disease traits and DCH (0.269-0.469) as well as STAT (0.172-0.242). Thus, the results indicate that larger and sharper animals with low BCS are more susceptible to metabolic disorders. Genome-wide association studies revealed several significantly associated SNPs for production and conformation traits, confirming previous findings from literature. Moreover, for production and conformation traits, shared significant signals on Bos taurus autosome (BTA) 5 (88.36 Mb) and BTA 6 (86.40 to 87.27 Mb) were found, implying pleiotropy. Additionally, significant SNPs were observed for metabolic diseases on BTA 3, 10, 14, 17, and 26 in first lactation and on BTA 2, 6, 8, 17, and 23 in third lactation. Overall, this study provides important insights into the genetic basis and interrelations of relevant traits in today's Holstein cattle breeding programs, and findings may help to improve selection decisions.

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.

Are subjectively scored linear type traits suitable predictors of the genetic merit for feed intake in grazing Holstein-Friesian dairy cows?

Measuring dry matter intake (DMI) in grazing dairy cows using currently available techniques is invasive, time consuming, and expensive. An alternative to directly measuring DMI for use in genetic evaluations is to identify a set of readily available animal features that can be used in a multitrait genetic evaluation for DMI. The objectives of the present study were thus to estimate the genetic correlations between readily available body-related linear type traits and DMI in grazing lactating Holstein-Friesian cows, but importantly also estimate the partial genetic correlations between these linear traits and DMI, after adjusting for differences in genetic merit for body weight. Also of interest was whether the predictive ability derived from the estimated genetic correlations materialized upon validation. After edits, a total of 8,055 test-day records of DMI, body weight, and milk yield from 1,331 Holstein-Friesian cows were available, as were chest width, body depth, and stature from 47,141 first lactation Holstein-Friesian cows. In addition to considering the routinely recorded linear type traits individually, novel composite traits were defined as the product of the linear type traits as an approximation of rumen volume. All linear type traits were moderately heritable, with heritability estimates ranging from 0.27 (standard error = 0.14) to 0.49 (standard error = 0.15); furthermore, all linear type traits were genetically correlated (0.29 to 0.63, standard error 0.14 to 0.12) with DMI. The genetic correlations between the individual linear type traits and DMI, when adjusted for genetic differences in body weight, varied from -0.51 (stature) to 0.48 (chest width). These genetic correlations between DMI and linear type traits suggest linear type traits may be useful predictors of DMI, even when body weight information is available. Nonetheless, estimated genetic merit of DMI derived from a multitrait genetic evaluation of linear type traits did not correlate strongly with actual DMI in a set of validation animals; the benefit was even less if body weight data were also available.

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.

Contribution of genetic variability to phenotypic differences in on-farm efficiency metrics of dairy cows based on body weight and milk solids yield

Milk solids per kilogram of body weight (BW) is growing in popularity as a measure of dairy cow lactation efficiency. Little is known on the extent of genetic variability that exist in this trait but also the direction and strength of genetic correlations with other performance traits. Such genetic correlations are important to know if producers are to consider actively selecting cows excelling in milk solids per kilogram of BW. The objective of the present study was to use a large data set of commercial Irish dairy cows to quantify the extent of genetic variability in milk solids per kilogram of BW and related traits but also their genetic and phenotypic inter-relationships. Mid-lactation BW and body condition score (BCS), along with 305-d milk solids yield (i.e., fat plus protein yield) were available on 12,413 lactations from 11,062 cows in 85 different commercial dairy herds. (Co)variance components were estimated using repeatability animal linear mixed models. The genetic correlation between milk solids and body weight was only 0.05, which when coupled with the observed large genetic variability in both traits, indicate massive potential to select for both traits in opposite directions. The genetic correlations between both milk solids and BW with BCS; however, need to be considered in any breeding strategy. The genetic standard deviation, heritability, and repeatability of milk solids per kilogram of BW was 0.08, 0.37, and 0.57, respectively. The genetic correlation between milk solids per kilogram of BW with milk solids, BW, and BCS was 0.62, -0.75, and -0.41, respectively. Therefore, based on genetic regression, each increase of 0.10 units in genetic merit for milk solids per kilogram of BW is expected to result in, on average, an increase in 16.1 kg 305-d milk solids yield, a reduction of 25.6 kg of BW and a reduction of 0.05 BCS units (scale of 1-5 where 1 is emaciated). The genetic standard deviation (heritability) for 305-d milk solids yield adjusted phenotypically to a common BW was 27.3 kg (0.22). The genetic correlation between this adjusted milk solids trait with milk solids, BW, and BCS was 0.91, -0.12, and -0.26, respectively. Once also adjusted phenotypically to a common BCS, the genetic standard deviation (heritability) for milk solids adjusted phenotypically to a common BW was 26.8 kg (0.22) where the genetic correlation with milk solids, BW and BCS was 0.91, -0.21, and -0.07, respectively. The genetic standard deviation (heritability) of BW adjusted phenotypically for differences in milk solids was 35.3 kg (0.61), which reduced to 33.2 kg when also phenotypically adjusted for differences in BCS. Results suggest considerable opportunity exists to change milk solids yield independent of BW, and vice versa. The opportunity is reduced slightly once also corrected for differences in BCS. Inter-animal BCS differences should be considered if selection on such metrics is contemplated.

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.

Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle

Feed efficiency and energy balance are important traits underpinning profitability and environmental sustainability in animal production. They are complex traits, and our understanding of their underlying biology is currently limited. One measure of feed efficiency is residual feed intake (RFI), which is the difference between actual and predicted intake. Variation in RFI among individuals is attributable to the metabolic efficiency of energy utilization. High RFI (H_RFI) animals require more energy per unit of weight gain or milk produced compared with low RFI (L_RFI) animals. Energy balance (EB) is a closely related trait calculated very similarly to RFI. Cellular energy metabolism in mitochondria involves mitochondrial protein (MiP) encoded by both nuclear (NuMiP) and mitochondrial (MtMiP) genomes. We hypothesized that MiP genes are differentially expressed (DE) between H_RFI and L_RFI animal groups and similarly between negative and positive EB groups. Our study aimed to characterize MiP gene expression in white blood cells of H_RFI and L_RFI cows using RNA sequencing to identify genes and biological pathways associated with feed efficiency in dairy cattle. We used the top and bottom 14 cows ranked for RFI and EB out of 109 animals as H_RFI and L_RFI, and positive and negative EB groups, respectively. The gene expression counts across all nuclear and mitochondrial genes for animals in each group were used for differential gene expression analyses, weighted gene correlation network analysis, functional enrichment, and identification of hub genes. Out of 244 DE genes between RFI groups, 38 were MiP genes. The DE genes were enriched for the oxidative phosphorylation (OXPHOS) and ribosome pathways. The DE MiP genes were underexpressed in L_RFI (and negative EB) compared with the H_RFI (and positive EB) groups, suggestive of reduced mitochondrial activity in the L_RFI group. None of the MtMiP genes were among the DE MiP genes between the groups, which suggests a non-rate limiting role of MtMiP genes in feed efficiency and warrants further investigation. The role of MiP, particularly the NuMiP and OXPHOS pathways in RFI, was also supported by our gene correlation network analysis and the hub gene identification. We validated the findings in an independent data set. Overall, our study suggested that differences in feed efficiency in dairy cows may be linked to differences in cellular energy demand. This study broadens our knowledge of the biology of feed efficiency in dairy cattle.

Genetic parameters for methane emission traits in Australian dairy cows

Methane is a greenhouse gas of high interest to the dairy industry, with 57% of Australia's dairy emissions attributed to enteric methane. Enteric methane emissions also constitute a loss of approximately 6.5% of ingested energy. Genetic selection offers a unique mitigation strategy to decrease the methane emissions of dairy cattle, while simultaneously improving their energy efficiency. Breeding objectives should focus on improving the overall sustainability of dairy cattle by reducing methane emissions without negatively affecting important economic traits. Common definitions for methane production, methane yield, and methane intensity are widely accepted, but there is not yet consensus for the most appropriate method to calculate residual methane production, as the different methods have not been compared. In this study, we examined 9 definitions of residual methane production. Records of individual cow methane, dry matter intake (DMI), and energy corrected milk (ECM) were obtained from 379 animals and measured over a 5-d period from 12 batches across 5 yr using the SF₆ tracer method and an electronic feed recording system, respectively. The 9 methods of calculating residual methane involved genetic and phenotypic regression of methane production on a combination of DMI and ECM corrected for days in milk, parity, and experimental batch using phenotypes or direct genomic values. As direct genomic values (DGV) for DMI are not routinely evaluated in Australia at this time, DGV for FeedSaved, which is derived from DGV for residual feed intake and estimated breeding value for bodyweight, were used. Heritability estimates were calculated using univariate models, and correlations were estimated using bivariate models corrected for the fixed effects of year-batch, days in milk, and lactation number, and fitted using a genomic relationship matrix. Residual methane production candidate traits had low to moderate heritability (0.10 ± 0.09 to 0.21 ± 0.10), with residual methane production corrected for ECM being the highest. All definitions of residual methane were highly correlated phenotypically (>0.87) and genetically (>0.79) with one another and moderately to highly with other methane candidate traits (>0.59), with high standard errors. The results suggest that direct selection for a residual methane production trait would result in indirect, favorable improvement in all other methane traits. The high standard errors highlight the importance of expanding data sets by measuring more animals for their methane emissions and DMI, or through exploration of proxy traits and combining data via international collaboration.

Bayesian estimation of genetic variance and response to selection on linear or ratio traits of feed efficiency in dairy cattle

This study aimed to estimate genetic parameters of the linear trait genetic residual feed intake (RFI) and the ratio traits feed conversion ratio (FCR) and feed conversion efficiency (FCE) along with dry matter intake (DMI) and energy sink traits such as energy-corrected milk (ECM), body weight (BW), body condition score (BCS), and BW change (BWC) across different weeks in the first lactation of Danish Holstein cows. A second objective was to conduct a Bayesian analysis of direct and correlated superiority of the selected group when selecting on genetic RFI, FCR, or FCE. Feed intake and energy sink traits were recorded during wk 1 to 44 of lactation on 847 primiparous Danish Holstein cows. A Bayesian multivariate random regression animal model was used to analyze DMI, ECM, BW, and BCS in different weeks of lactation. Genetic RFI was obtained by conditioning DMI on ECM, BW, BCS, and BWC using genetic partial regression coefficients. The posterior distribution of the breeding values for FCR and FCE was derived from the posterior distribution of functions of "fixed" environmental effects and random additive genetic effects on DMI and ECM. Genetic superiority of the selected group was defined as the difference in additive genetic mean of the selected top individuals expected to be potential parents, and the total population after integrating genetic trends out of the posterior distribution of selection responses. Posterior means of heritability of genetic RFI ranged from 0.10 to 0.15, genetic variance of FCR and FCE ranged from 2.13 × 10^-3 to 3.2 × 10^-3 (kg² DMI/kg² ECM) and 6.11 × 10^-3 to 2.4 × 10^-2 (kg² ECM/kg² DMI), respectively. Selection against RFI showed a direct response of -1.01 to -2.23 kg/d RFI and correlated responses of -0.031 to -0.056 kg/kg for FCR, 0.104 to 0.160 kg/kg for FCE, and -0.316 to -1.057 kg/d for DMI in different weeks of lactation. Selection against RFI had no significant effect on production traits but selection for ratio traits reduced BW and BCS. Posterior means of genetic correlation between DMI and ratio traits were low. In conclusion, the Bayesian procedure allowed us to estimate genetic RFI without the need for separate multiple regression analysis and considered the non-normal posterior distribution of ratio traits. Selection against genetic RFI might be an effective means to improve feed efficiency compared with ratio traits for feed efficiency in dairy cattle.

Symposium review: Decomposing efficiency of milk production and maximizing profit

The dairy industry has focused on maximizing milk yield, as it is believed that this maximizes profit mainly through dilution of maintenance costs. Efficiency of milk production has received, until recently, considerably less attention. The most common method to determine biological efficiency of milk production is feed efficiency (FE), which is defined as the amount of milk produced relative to the amount of nutrients consumed. Economic efficiency is best measured as income over feed cost or gross margin obtained from feed investments. Feed efficiency is affected by a myriad of factors, but overall they could be clustered as follows: (1) physiological status of the cow (e.g., age, state of lactation, health, level of production, environmental conditions), (2) digestive function (e.g., feeding behavior, passage rate, rumen fermentation, rumen and hindgut microbiome), (3) metabolic partitioning (e.g., homeorhesis, insulin sensitivity, hormonal profile), (4) genetics (ultimately dictating the 2 previous aspects), and (5) nutrition (e.g., ration formulation, nutrient balance). Over the years, energy requirements for maintenance seem to have progressively increased, but efficiency of overall nutrient use for milk production has also increased due to dilution of nutrient requirements for maintenance. However, empirical evidence from the literature suggests that marginal increases in milk require progressively greater marginal increases in nutrient supply. Thus, the dilution of maintenance requirements associated with increases in production is partially overcome by a progressive diminishing marginal biological response to incremental energy and protein supplies. Because FE follows the law of diminishing returns, and because marginal feed costs increase progressively with milk production, profits associated with improving milk yield might, in some cases, be considerably lower than expected.

Factors associated with the weight of individual primal cuts and their inter-relationship in cattle

Input parameters for decision support tools are comprised of, amongst others, knowledge of the associated factors and the extent of those associations with the animal-level feature of interest. The objective of the present study was to quantify the association between animal-level factors with primal cut yields in cattle and to understand the extent of the variability in primal cut yields independent carcass weight. The data used consisted of the weight of 14 primal carcass cuts (as well as carcass weight, conformation, and fat score) on up to 54,250 young cattle slaughtered between the years 2013 and 2017. Linear mixed models, with contemporary group of herd-sex-season of slaughter as a random effect, were used to quantify the associations between a range of model fixed effects with each primal cut separately. Fixed effects in the model were dam parity, heterosis coefficient, recombination loss, a covariate per breed representing the proportion of Angus, Belgian Blue, Charolais, Jersey, Hereford, Limousin, Simmental, and Holstein-Friesian and a three-way interaction between whether the animal was born in a dairy or beef herd, sex, and age at slaughter, with or without carcass weight as a covariate in the mixed model. The raw correlations among all cuts were all positive varying from 0.33 (between the bavette and the striploin) to 0.93 (between the topside and knuckle). The partial correlation among cuts, following adjustment for differences in carcass weight, varied from -0.36 to 0.74. Age at slaughter, sex, dam parity, and breed were all associated (P < 0.05) with the primal cut weight. Knowledge of the relationship between the individual primal cuts, and the solutions from the models developed in the study, could prove useful inputs for decision support systems to increase performance.

Energy profiling of dairy cows from routine milk mid-infrared analysis

The balance of body energy within and across lactations can have health and fertility consequences for the dairy cow. This study aimed to create a large calibration data set of dairy cow body energy traits across the cow's productive life, with concurrent milk mid-infrared (MIR) spectral data, to generate a prediction tool for use in commercial dairy herds. Detailed phenotypic data from 1,101 Holstein Friesian cows from the Langhill research herd (SRUC, Scotland) were used to generate energy balance (EB) and effective energy intake (EI), both in megajoules per day. Pretreatment of spectral data involved standardization to account for drift over time and machine. Body energy estimates were aligned with their spectral data to generate a prediction of these traits based on milk MIR spectroscopy. After data edits, partial least squares analysis generated prediction equations with a coefficient of determination from split sample 10-fold cross validation of 0.77 and 0.75 for EB and EI, respectively. These prediction equations were applied to national milk MIR spectra on over 11 million animal test dates (January 2013 to December 2016) from 4,453 farms. The predictions generated from these were subject to phenotypic analyses with a fixed regression model highlighting differences between the main dairy breeds in terms of energy traits. Genetic analyses generated heritability estimates for EB and EI ranging from 0.12 to 0.17 and 0.13 to 0.15, respectively. This study shows that MIR-based predictions from routinely collected national data can be used to generate predictions of dairy cow energy turnover profiles for both animal management and genetic improvement of such difficult and expensive-to-record traits.

High-density genome-wide association study for residual feed intake in Holstein dairy cattle

Improving feed efficiency (FE) of dairy cattle may boost farm profitability and reduce the environmental footprint of the dairy industry. Residual feed intake (RFI), a candidate FE trait in dairy cattle, can be defined to be genetically uncorrelated with major energy sink traits (e.g., milk production, body weight) by including genomic predicted transmitting ability of such traits in genetic analyses for RFI. We examined the genetic basis of RFI through genome-wide association (GWA) analyses and post-GWA enrichment analyses and identified candidate genes and biological pathways associated with RFI in dairy cattle. Data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States. Of these cows, 3,555 were genotyped and were imputed to a high-density list of 312,614 SNP. We used a single-step GWA method to combine information from genotyped and nongenotyped animals with phenotypes as well as their ancestors' information. The estimated genomic breeding values from a single-step genomic BLUP were back-solved to obtain the individual SNP effects for RFI. The proportion of genetic variance explained by each 5-SNP sliding window was also calculated for RFI. Our GWA analyses suggested that RFI is a highly polygenic trait regulated by many genes with small effects. The closest genes to the top SNP and sliding windows were associated with dry matter intake (DMI), RFI, energy homeostasis and energy balance regulation, digestion and metabolism of carbohydrates and proteins, immune regulation, leptin signaling, mitochondrial ATP activities, rumen development, skeletal muscle development, and spermatogenesis. The region of 40.7 to 41.5 Mb on BTA25 (UMD3.1 reference genome) was the top associated region for RFI. The closest genes to this region, CARD11 and EIF3B, were previously shown to be related to RFI of dairy cattle and FE of broilers, respectively. Another candidate region, 57.7 to 58.2 Mb on BTA18, which is associated with DMI and leptin signaling, was also associated with RFI in this study. Post-GWA enrichment analyses used a sum-based marker-set test based on 4 public annotation databases: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, Reactome pathways, and medical subject heading (MeSH) terms. Results of these analyses were consistent with those from the top GWA signals. Across the 4 databases, GWA signals for RFI were highly enriched in the biosynthesis and metabolism of amino acids and proteins, digestion and metabolism of carbohydrates, skeletal development, mitochondrial electron transport, immunity, rumen bacteria activities, and sperm motility. Our findings offer novel insight into the genetic basis of RFI and identify candidate regions and biological pathways associated with RFI in dairy cattle.

Determining the economic value of daily dry matter intake and associated methane emissions in dairy cattle

Feed represents a substantial proportion of production costs in the dairy industry and is a useful target for improving overall system efficiency and sustainability. The objective of this study was to develop methodology to estimate the economic value for a feed efficiency trait and the associated methane production relevant to Canada. The approach quantifies the level of economic savings achieved by selecting animals that convert consumed feed into product while minimizing the feed energy used for inefficient metabolism, maintenance and digestion. We define a selection criterion trait called Feed Performance (FP) as a 1 kg increase in more efficiently used feed in a first parity lactating cow. The impact of a change in this trait on the total lifetime value of more efficiently used feed via correlated selection responses in other life stages is then quantified. The resulting improved conversion of feed was also applied to determine the resulting reduction in output of emissions (and their relative value based on a national emissions value) under an assumption of constant methane yield, where methane yield is defined as kg methane/kg dry matter intake (DMI). Overall, increasing the FP estimated breeding value by one unit (i.e. 1 kg of more efficiently converted DMI during the cow's first lactation) translates to a total lifetime saving of 3.23 kg in DMI and 0.055 kg in methane with the economic values of CAD $0.82 and CAD $0.07, respectively. Therefore, the estimated total economic value for FP is CAD $0.89/unit. The proposed model is robust and could also be applied to determine the economic value for feed efficiency traits within a selection index in other production systems and countries.

Reliability of breeding values for feed intake and feed efficiency traits in dairy cattle: When dry matter intake recordings are sparse under different scenarios

Currently, routine recordings of dry matter intake (DMI) in commercial herds are practically nonexistent. Recording DMI from commercial herds is a prerequisite for the inclusion of feed efficiency (FE) traits in dairy cattle breeding goals. To develop future on-farm phenotyping strategies, recording strategies that are low cost and less demanding logistically and that give relatively accurate estimates of the animal's genetic merit are therefore needed. The objectives of this study were (1) to estimate genetic parameters for daily DMI and FE traits and use the estimated parameters to simulate daily DMI phenotypes under different DMI recording scenarios (SCN) and (2) to use the simulated data to estimate for different scenarios the associated reliability of estimated breeding value and accuracies of genomic prediction for varying sizes of reference populations. Five on-farm daily DMI recording scenarios were simulated: once weekly (SCN1), once monthly (SCN2), every 2 mo (SCN3), every 3 mo (SCN4), and every 4 mo (SCN5). To estimate reliability of estimated breeding values, DMI and FE observations and true breeding values were simulated based on variance components estimated for daily observations of Nordic Red cows. To emulate realistic on-farm recording, 5 data set replicates, each with 36,037 DMI and FE records, were simulated for real pedigree and data structure of 789 Holstein cows. Observations for the 5 DMI recording scenarios were generated by discarding data in a step-wise manner from the full simulated data per the scenario's definitions. For each of these scenarios, reliabilities were calculated as correlation between the true and estimated breeding values. Variance components and genetic parameters were estimated for daily DMI, residual feed intake (RFI), and energy conversion efficiency (ECE) fitting the random regression model. Data for variance components were from 227 primiparous Nordic Red dairy cows covering 8 to 280 d in milk. Lactation-wise heritability for DMI, RFI, and ECE was 0.33, 0.12, and 0.32, respectively, and daily heritability estimates during lactation ranged from 0.18 to 0.45, 0.08 to 0.32, and 0.08 to 0.45 for DMI, RFI, and ECE, respectively. Genetic correlations for DMI between different stages of lactation ranged from -0.50 to 0.99. The comparison of different on-farm DMI recording scenarios indicated that adopting a less-frequent recording scenario (SCN3) gave a similar level of accuracy as SCN1 when 17 more daughters are recorded per sire over the 46 needed for SCN1. Such a strategy is less demanding logistically and is low cost because fewer observations need to be collected per animal. The accuracy of genomic predictions associated with the 5 recording scenarios indicated that setting up a relatively larger reference population and adopting a less-frequent DMI sampling scenario (e.g., SCN3) is promising. When the same reference population size was considered, the genomic prediction accuracy of SCN3 was only 5.0 to 7.0 percentage points lower than that for the most expensive DMI recording strategy (SCN1). We concluded that DMI recording strategies that are sparse in terms of records per cow but with slightly more cows recorded per sire are advantageous both in genomic selection and in traditional progeny testing schemes when accuracy, logistics, and cost implications are considered.

Little genetic variability in resilience among cattle exists for a range of performance traits across herds in Ireland differing in Fasciola hepatica prevalence

It is anticipated that in the future, livestock will be exposed to a greater risk of infection from parasitic diseases. Therefore, future breeding strategies for livestock, which are generally long-term strategies for change, should target animals adaptable to environments with a high parasitic load. Covariance components were estimated in the present study for a selection of dairy and beef performance traits over herd-years differing in Fasciola hepatica load using random regression sire models. Herd-year prevalence of F. hepatica was determined by using F. hepatica-damaged liver phenotypes which were recorded in abattoirs nationally. The data analyzed consisted up to 83,821 lactation records from dairy cows for a range of milk production and fertility traits, as well as 105,054 young animals with carcass-related information obtained at slaughter. Reaction norms for individual sires were derived from the random regression coefficients. The heritability and additive genetic standard deviations for all traits analyzed remained relatively constant as herd-year F. hepatica prevalence gradient increased up to a prevalence level of 0.7; although there was a large increase in heritability and additive genetic standard deviation for milk and fertility traits in the observed F. hepatica prevalence levels >0.7, only 5% of the data existed in herd-year prevalence levels >0.7. Very little rescaling, therefore, exists across differing herd-year F. hepatica prevalence levels. Within-trait genetic correlations among the performance traits across different herd-year F. hepatica prevalence levels were less than unity for all traits. Nevertheless, within-trait genetic correlations for milk production and carcass traits were all >0.8 for F. hepatica prevalence levels between 0.2 and 0.8. The lowest estimate of within-trait genetic correlations for the different fertility traits ranged from -0.03 (SE = 1.09) in age of first calving to 0.54 (SE = 0.22) for calving to first service interval. Therefore, there was reranking of sires for fertility traits across different F. hepatica prevalence levels. In conclusion, there was little or no genetic variability in sensitivity to F. hepatica prevalence levels among cattle for milk production and carcass traits. But, some genetic variability in sensitivity among dairy cows did exist for fertility traits measured across herds differing in F. hepatica prevalence.