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

Large datasets allow estimation of 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 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 DIM, age-parity subclass, trial date, management group, and BW 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, respectively) 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 ECM 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, whereas the sire genomic regressions indicated approximately 35% more DMI for fat than protein. Estimates of annual maintenance in kilograms DMI/kilograms BW per 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 the National Academies for Science, Engineering, and Medicine 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 BW 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.

Genetic analysis of rumination time based on an analysis of 77,697 Israeli dairy cows

Reduction of methane emission may become necessary for sustainable milk production. Several studies indicate a relationship between rumination time and the level of methane emission. The objectives of the current study were to estimate environmental factors affecting daily rumination time in high-yielding dairy cattle, genetic parameters for rumination time across parities, and environmental and genetic correlations between rumination time and economic traits, and to predict the consequence of inclusion of this trait in the Israeli breeding index. The data included more than 30 million daily records from 77,697 Israeli Holstein cows for rumination time and milk production. A lactation measure of daily rumination time per cow was computed as the mean of the residuals from a linear model analysis with rumination time as the dependent variable. The independent variables were parity and the square root, linear, quadradic and inverse of DIM by parity. Because of the shape of the lactation curve for rumination time, separate linear model analyses were performed for records up to 40 DIM and records with >40 DIM. The phenotypic correlation between first- and second-parity lactations for rumination time was almost 0.8, and close to 0.7 for milk. The heritability of lactation rumination time was close to 0.44 for parities 1 to 3. Heritability for milk production decreased from 0.5 in first parity to 0.3 in third parity. For both traits, genetic correlations among parities were all >0.9. Thus, for routine genetic analysis of rumination time, records in the different parities can be considered the same trait. The genetic correlation between rumination time and milk on first parity was 0.25 and increased slightly with increase in parity. Genetic correlations between rumination time, based on the first 40 DIM, were economically unfavorable with retained placenta but economically favorable with metritis, ketosis, and displaced abomasum. Genetic correlations between rumination time and the 9 traits included in the Israeli breeding index (milk, fat, and protein production; SCS; female fertility; herd-life; milk production persistency; calving ease; and calf mortality) were all economically favorable, except for the correlation of 0.17 with SCS. With the current index, daily rumination time with a current mean of 536 min and SD of 90 min is expected to increase by 11 min/d after 10 yr of selection. Inclusion of this trait with a positive index weight equivalent to 10% of the index should increase rumination time by 19 min. All changes in expected gains due to inclusion of rumination time in the index were economically positive, except for fat and SCS. Inclusion of rumination time in the index should result in 1 kg less gain in fat, a miniscule gain of 0.03 for SCS; and gains of 1.5 kg protein, 0.3% female fertility, and 5 d herd-life. Even though the case for a genetic correlation between rumination time and methane emission is still weak, inclusion of this trait in the commercial index may be justified, considering that equipment is now commercially available for routine recording at reasonable cost.

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.

Genetic relationships between behavioral traits and feed efficiency traits in lactating Holstein cows

The evaluation of dairy cow feed efficiency using residual feed intake accounts for known energy sinks. However, behavioral traits may also contribute to the variation in feed efficiency. Our objective was to estimate the heritability and repeatability of behavioral traits and their genetic correlations with feed efficiency and its components in lactating Holstein cows. The first data set consisted of 36,075 daily rumination and lying time records collected using a SMARTBOW ear tag accelerometer (Zoetis, Parsippany, NJ) and 6,371 weekly feed efficiency records of 728 cows from the University of Wisconsin-Madison. The second data set consisted of 59,155 daily activity records, measured as number of steps, recorded by pedometers (AfiAct; S.A.E. Afikim, Kibbutz Afikim, Israel), and 8,626 weekly feed efficiency records of 635 cows from the University of Florida. Feed efficiency and its components included dry matter intake, change in body weight, metabolic body weight, secreted milk energy, and residual feed intake. The statistical models included the fixed effect of cohort, lactation number, and days in milk, and the random effects of animal and permanent environment. Heritability estimates for behavioral traits using daily records were 0.19 ± 0.06 for rumination and activity, and 0.37 ± 0.07 for lying time. Repeatability estimates for behavioral traits using daily data ranged from 0.56 ± 0.02 for activity to 0.62 ± 0.01 for lying time. Both heritability and repeatability estimates were larger when weekly records instead of daily records were used. Rumination and activity had positive genetic correlations with residual feed intake (0.40 ± 0.19 and 0.31 ± 0.22, respectively) while lying time had a negative genetic correlation with this residual feed intake (-0.27 ± 0.11). These results indicate that more efficient cows tend to spend more time lying and less time active. Additionally, less efficient cows tend to eat more and therefore also tend to ruminate longer. Overall, sensor-based behavioral traits are heritable and genetically correlated with feed efficiency and its components and, therefore, they could be used as indicators to identify feed efficient cows within the herd.

Estimation of genetic parameters and single-step genome-wide association studies for eating time and rumination time in Holstein dairy cows

The aims of this study were to estimate genetic parameters and to identify genomic regions associated with eating time (ET) and rumination time (RUT) in Holstein dairy cows. Genetic correlations among ET, RUT, and milk yield traits were also estimated. The data were collected from 2019 to 2022 in 6 dairy herds located in the Walloon Region of Belgium. The dataset consisted of daily ET and RUT records on 284 Holstein cows, from which 41 cows had records only for the first parity (P1), 101 cows had records from both the first and second parities, and 142 cows had records only for the second parity (P2). The number of daily ET and RUT records in the P1 and P2 cows were 18,569 (on 142 cows) and 34,464 (on 243 cows), respectively. Data on 28,994 SNPs located on 29 Bos taurus autosomes (BTA) of 747 animals (435 males) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 20-SNP sliding window (with an average size of 1.52 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Mean (standard deviation; SD) averaged daily ET and RUT were 327.0 (85.66) and 559.4 (77.69) min/d for cows in P1 and 316.0 (82.24) and 574.2 (75.42) min/d for cows in P2, respectively. Mean (standard deviation; SD) heritability estimates for daily ET and RUT were 0.42 (0.09) and 0.45 (0.06) for cows in P1 and 0.45 (0.04) and 0.43 (0.02) for cows in P2, respectively. Mean (SD) daily genetic correlations between daily ET and RUT were 0.27 (0.07) for P1 and 0.34 (0.08) for P2. Genome-wide association analyses identified 6 genomic regions distributed over 5 chromosomes (BTA1, BTA4, BTA11, 2 regions of BTA14, and BTA17) associated with ET or RUT. The findings of this study increase our preliminary understanding of the genetic background of feeding behavior in dairy cows; however, larger datasets are needed to determine whether ET and RUT might have the potential to be used in selection programs.

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.

Heritability and genetic correlations of rumination time with milk-yield and milking traits in Holstein-Friesian cows using an automated milking system

Knowledge of the values of genetic parameters is a prerequisite for conducting a breeding program. This is especially important for rumination, which is considered an indicator of cow's health. Exploring the genetic relations between rumination time, milk yield, and milking traits could make it a valuable tool in dairy cattle breeding strategies. The objective of the research was to estimate heritability, repeatability, and genetic and phenotypic correlations of rumination time (RT), as well as traits associated with milk yield and milking of dairy cows of the Polish Holstein-Friesian breed kept in herds equipped with an automatic milking system. The research takes into consideration daily results for milking in the first lactation and second lactation, from 1 486 cows of the breed milked between 2013 and 2015 year. Cows were housed in 24 free-stall barns and fed a Partial Mixed Ration feed. The barns had an automated milking system (Astronaut A4 - Lely Industry). The cows received a varied dose of the concentrate, either in the milking robot or the feeding station, depending on the level of their milk yield. Our research has shown that RT was a low heritable trait (0.140 ± 0.039) and had a medium repeatability (0.572 ± 0.007). We detected a positive genetic correlation between RT and milk yield (0.341); however, a statistically significant negative relationship was identified between RT and urea content (-0.418) in milk. Estimations of genetic correlations suggest that selecting for higher RT may correspond to reduced urea content in milk. Investigating the genetics aspect of RT and the relationship with milk yield and milking traits may turn this into one of the useful criterion selections for dairy cattle breeding strategies, but should be used carefully. Further analyses on larger data sets and different populations are necessary.

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.

Genetic parameters of behavior traits of beef cattle classified using wearable devices

With the development of wearable devices, it is now possible to monitor livestock behavior 24 h a day. In this study, we estimated the genetic parameters of the daily duration of six behaviors (feeding, moving, lying, standing, ruminating while lying, and ruminating while standing) in beef cattle, automatically classified using wearable devices. The devices were attached to 332 Japanese beef cattle at two stations for approximately 5 months. We compared repeatability, Poisson regression, and random regression models using the deviance information criterion. Poisson regression models were selected for all traits at each station, probably because of the non-normal distribution of the phenotypes. The heritability estimates by the Poisson regression models were moderate at each station: 0.67 and 0.68 for feeding, 0.68 and 0.53 for moving, 0.47 and 0.55 for lying, 0.45 and 0.40 for standing, 0.51 and 0.59 for ruminating while lying, and 0.37 and 0.45 for ruminating while standing. The genetic correlations between these traits were all negative at both stations, whereas the residual correlations showed different directions depending on the station. Although validation studies with larger populations are needed to confirm these findings, this study provides fundamental knowledge of the genetic basis of daily behavior in beef cattle.

Genetic parameter estimates for daily predicted gross feed efficiency and its association with energy-corrected milk in South African Holstein cattle

Genetic parameters for daily predicted gross feed efficiency (pGFE) and energy corrected milk (ECM) in the first three parities of South African Holstein cattle were estimated by repeatability animal models. Data comprised of 11,068 test-day milk production records of 1,575 Holstein cows that calved between 2009 and 2019. Heritability estimates for pGFE were 0.12 ± 0.06, 0.09 ± 0.04 and 0.18 ± 0.05 in early, mid and late lactation, respectively. Estimates were moderate for primiparous (0.21 ± 0.05) and low for multiparous (0.10 ± 0.04) cows. Heritability and repeatability across all lactations were 0.14 ± 0.03 and 0.37 ± 0.03, respectively. Genetic correlations between pGFE in different stages of lactation ranged from 0.87 ± 0.24 (early and mid) to 0.97 ± 0.28 (early and late), while a strong genetic correlation (0.90 ± 0.03) was found between pGFE and ECM, across all lactations. The low to moderate heritability estimates for pGFE suggest potential for genetic improvement of the trait through selection, albeit with a modest accuracy of selection. The high genetic correlation of pGFE with ECM may, however, assist to improve accuracy of selection for feed efficiency by including both traits in multi-trait analyses. These genetic parameters may be used to estimate breeding values for pGFE, which will enable the trait to be incorporated in the breeding objective for South African Holstein cattle.

Individual animal variability in rumination, activity, and lying behavior during the periparturient period of dairy cattle

The aim of the current study was to determine individual animal variability in rumination, activity, and lying behavior during the periparturient period within the context of dairy cattle nutrition, social, and physical environment. Holstein animals (nulliparous = 77, parous = 219) from one sand-bedded, freestall dairy in northwest Wisconsin were enrolled -17 d in milk (DIM, d 0 = calving), when they were fitted with an automated monitoring device (Hi-Tag, SCR Engineers Ltd.). At -11 DIM, animals were fitted with HOBO Pendant G Data Loggers. The HOBO Pendant G Data Loggers were fitted 6 d later because they were set up to collect data for 22 d (d -11 to 11), to avoid constant handling of the animals that could alter their behavior. Prepartum, nulliparous and parous animals were housed separately. Postpartum (1 to 17 ± 3 DIM), primiparous and multiparous cows were commingled. Samples of the total mixed ration were submitted for wet chemistry analysis and determination of physically effective NDF (peNDF). Temperature and humidity data were collected using RH Temp probes (HOBO Pro Series) installed in each of the pens, and the percentages of 30-min intervals within a day with temperature-humidity index ≥68 (PctTHI68) were calculated. Stocking density (cows per stall) during the pre- and postpartum periods were calculated daily. Prepartum data from nulliparous and parous animals were analyzed separately, and postpartum data from primiparous and multiparous animals were analyzed together. Prepartum, nulliparous and parous animals explained 83.9 and 64.5% of the variability in rumination, 70.7 and 60.9% of the variability in activity, and 38.1 and 63.6% of the variability in lying time, respectively. Postpartum, animal explained 49.7, 56.8, and 35.6% of the variability in rumination, activity, and lying time, respectively. Although stocking density, PctTHI68, peNDF, crude protein, and ether extract were associated with the variability in rumination, activity, and lying time, they explained ≤6.6% of the daily variability in these behaviors. We conclude that, within the conditions of the collaborating commercial herd, individual animal is the most important factor explaining daily variability in rumination, activity, and lying time.

Genetic Analysis of Methane Emission Traits in Holstein Dairy Cattle

Genetic selection can be a feasible method to help mitigate enteric methane emissions from dairy cattle, as methane emission-related traits are heritable and genetic gains are persistent and cumulative over time. The objective of this study was to estimate heritability of methane emission phenotypes and the genetic and phenotypic correlations between them in Holstein cattle. We used 1765 individual records of methane emission obtained from 330 Holstein cattle from two Canadian herds. Methane emissions were measured using the GreenFeed system, and three methane traits were analyzed: the amount of daily methane produced (g/d), methane yield (g methane/kg dry matter intake), and methane intensity (g methane/kg milk). Genetic parameters were estimated using univariate and bivariate repeatability animal models. Heritability estimates (±SE) of 0.16 (±0.10), 0.27 (±0.12), and 0.21 (±0.14) were obtained for daily methane production, methane yield, and methane intensity, respectively. A high genetic correlation (rg = 0.94 ± 0.23) between daily methane production and methane intensity indicates that selecting for daily methane production would result in lower methane per unit of milk produced. This study provides preliminary estimates of genetic parameters for methane emission traits, suggesting that there is potential to mitigate methane emission in Holstein cattle through genetic selection.

Individual differences in digesta retention and their relation to chewing in cattle-A pilot investigation

While information on individual differences in digesta mean retention time (MRT) might be interesting when selecting phenotypes for digestive efficiency, MRT measurements are prohibitively labour-intensive for large-scale application. Therefore, more easily measured proxies of MRT might be helpful. We used the opportunity of an experiment applying saliva stimulant in cattle to investigate the effect of different individual chewing behaviour on fluid and particle MRT with a consistent diet. Four non-lactating cattle (670-850 kg body mass [BM]) were used in a 4 × 4 Latin square design, treated with the saliva stimulant pilocarpine in dosages of 0, 1, 2.5 and 5 mg/kg BM per day. The cattle were fed hay with dry matter intake (DMI) assigned according to their metabolic body weight. MRT in the whole gastrointestinal tract (GIT), the reticulorumen (RR) and the distal tract were measured using Co-EDTA, Cr-mordanted fibre and La-mordanted fibre as markers representing fluid, small particles (2 mm) and large particles (1 cm), respectively. The chewing behaviour was measured via noseband pressure sensor and expressed as chewing frequency (chews per time) and chewing intensity (chews per DMI), both for total chewing (ingestion plus rumination) and rumination chewing alone. The animals differed considerably in chewing behaviour and MRT measures. BM did not show a significant effect on chewing behaviour and MRT measures, though it tended to negatively correlated to total chewing intensity. Chewing intensity exerted a significant negative influence on MRT of fluid and particles in the RR, which was not the case for chewing frequency. Chewing frequency showed a significant relationship with MRT of large particles in the GIT. We suggest that chewing behaviour could influence MRT in two ways: (i) by affecting saliva production via the masticatory-salivary reflex and subsequently, the fluid inflow to the RR; (ii) by contributing to particle size reduction. Should the link between chewing behaviour and MRT be corroborated in larger studies, chewing measures, with their large interindividual variation, could emerge as an easy-to-measure proxy for MRT characteristics.

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.

Association of Feed Efficiency, Feeding Rate, and Behaviour with the Milk Performance of Dairy Cows

Identification of the associations of cow feed efficiency with feeding behaviour and milk production is important for supporting recommendations of strategies that optimise milk yield. The objective of this study was to identify associations between measures of feed efficiency, feed intake, feeding rate, rumination time, feeding time, and milk production using data collected from 26 dairy cows during a 3 month period in 2018. Cows averaged (mean ± standard deviation) 2.2 ± 1.7 lactations, 128 ± 40 days in milk, 27.5 ± 5.5 kg/day milk, 1.95 ± 0.69 kg feed/1 kg milk—the measure used to express feed conversion ratio (FCR), 575 ± 72 min/day rumination time, and 264 ± 67 min/day feeding time during the observation period. The coefficient of variation for rumination time (min/d) was 12.5%. A mixed linear model was selected for analyses. The most feed inefficient cows with the highest FCR (≥2.6 kg feed/1 kg milk) showed the lowest milk yield (24.8 kg/day), highest feed intake (78.8 kg), highest feeding rate (0.26 kg/min) and BCS (3.35 point). However, the relative milk yield (milk yield per 100 kg of body weight) was the highest (4.01 kg/day) in the most efficient group with the lowest FCR (≤1.4 kg feed/1 kg milk). Our study showed that the most efficient cows with the lowest FCR (≤1.4 kg feed/1 kg milk) had the highest rumination time (597 min/day; p < 0.05), feeding time (298 min/day; p < 0.05), rumination/activity ratio (4.39; p < 0.05) and rumination/feeding ratio (2.04; p < 0.05). Less active cows (activity time 164 min/day; p < 0.05) were the most efficient cows with the lowest FCR (≤1.4 kg feed/1 kg milk). The behavioural differences observed in this study provide new insight into the association of feed behaviour and feed efficiency with milk performance. Incorporating feeding behaviour into the dry matter intake model can improve its accuracy in the future and benefit breeding programmes.

Herd-level versus animal-level variation in methane emission prediction in grazing dairy cattle

In response to the increased concern over agriculture's contribution to greenhouse gas (GHG) emissions, more detailed assessments of current methane emissions and their variation, within and across individual dairy farms and cattle, are of interest for research and policy development. This assessment will provide insights into possible changes needed to reduce GHG emissions, the nature and direction of these changes, ways to influence farmer behavior and areas to maximize the adoption of emerging mitigation technologies. The objectives of this study were to (1) quantify the variation in enteric fermentation methane emissions within and among seasonal calving dairy farms with the majority of nutritional requirements met through grazed pasture; (2) use this variation to assess the potential of new individual animal emission monitoring technologies and their impact on mitigation policy. We used a large database of cow performance records for milk production and survival from 2 398 herds in New Zealand, and simulation to account for unobserved variation in feed efficiency and methane emissions per unit of feed. Results showed an average of 120 ± 31.4 kg predicted methane (CH₄) per cow per year after accounting for replacement costs, ranging 8.9-323 kg CH₄/cow per year. Whereas milk production, survival and predicted live weight were reasonably effective at predicting both individual and herd average levels of per cow feed intake, substantial within animal variation in emissions per unit of feed reduced the ability of these variables to predict variation in per animal methane output. Animal-level measurement technologies predicting only feed intake but not emissions per unit of feed are unlikely to be effective for advancing national policy goals of reducing dairy farming enteric methane output. This is because farmers seek to profitably utilize all farm feed resources available, so improvements in feed efficiency will not result in the reduction in feed utilization required to reduce methane emissions. At a herd level, average per cow milk production and live weight could form the basis of assigning a farm-level point of obligation for methane emissions. In conclusion, a comprehensive national database infrastructure that was tightly linked to animal identification and movement systems, and captured live weight data from existing farm-level recording systems, would be required to make this effective. Additional policy and incentivization mechanisms would still be required to encourage farmer uptake of mitigation interventions, such as novel feed supplements or vaccines that reduce methane emissions per unit of feed.

Unraveling the Relationship between Milk Yield and Quality at the Test Day with Rumination Time Recorded by a PLF Technology

Simple Summary

Precision livestock farming, by real time monitoring of dairy cows, has the potential to generate a huge amount of data to be used for farm management purposes, as well as in breeding programs. Daily rumination time (RT) recorded by commercial systems is promising in this context because it may be related to individual milk yield and composition. However, it is necessary to assess the ability of sensor data to be used in a predictive model, but also to evaluate and standardize the correct phenotypes, and how they are related to individual variability rather than from other sources. RT data and milk test day (TD) records collected from 691 cows, monitored for thirteen months, were analyzed for the already mentioned goals and to better characterize the effect of high-, medium- and low-level daily RT on milk yield and composition. Our results showed that “animal” in a farm major contributed to the RT total variability, confirming a possible use in breeding program. The higher RT class reported the best productive performance for milk and each solid yield, in spite of a small reduction in their contents, and appears to be related to a higher degree of saturation in the fatty acid profile.

Abstract

The study aimed to estimate the components of rumination time (RT) variability recorded by a neck collar sensor and the relationship between RT and milk composition. Milk test day (TD) and RT data were collected from 691 cows in three farms. Daily RT data of each animal were averaged for 3, 7, and 10 days preceding the TD date (RT_D). Variance component analysis of RT_D, considering the effects of farm, cow, parity, TD date, and lactation phase, showed that a farm, followed by a cow, had major contributions to the total variability. The RT₁₀ variable best performed on TD milk yield and quality records across models by a multi-model inference approach and was adopted to study its relationship with milk traits, by linear mixed models, through a 3-level stratification: low (LRT₁₀ ≤ 8 h/day), medium (8 h/day < MRT₁₀ ≤ 9 h/day), and high (HRT₁₀ > 9 h/day) RT. Cows with HRT₁₀ had greater milk, fat, protein, casein, and lactose daily yield, and lower fat, protein, casein contents, and fat to protein ratio compared to MRT₁₀ and LRT₁₀. Higher percentages of saturated fatty acid and lower unsaturated and monounsaturated fatty acid were found in HRT₁₀, with respect to LRT₁₀ and MRT₁₀ observations.

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.

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.

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.

Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows

Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.

Is rumination time an indicator of methane production in dairy cows?

As long as large-scale recording of expensive-to-measure and labor-consuming traits, such as dry matter intake (DMI) and CH₄ production (CH₄P), continues to be challenging in practical conditions, alternative traits that are already routinely recorded in dairy herds should be investigated. An ideal indicator trait must, in addition to expressing genetic variation, have a strong correlation with the trait of interest. Our aim was to estimate individual level and phenotypic correlations between rumination time (RT), CH₄P, and DMI to determine if RT could be used as an indicator trait for CH₄P and DMI. Data from 343 Danish Holstein cows were collected at the Danish Cattle Research Centre for a period of approximately 3 yr. The data set consisted of 14,890 records for DMI, 15,835 for RT, and 6,693 for CH₄P. Data were divided in primiparous cows only (PC) and all cows (MC), and then divided in lactation stage (early, mid, late, and whole lactation) to analyze the changes over lactation. Linear mixed models, including an animal effect but no pedigree, were used to estimate the correlations among traits. Phenotypic and individual level correlations between RT and both CH₄P and DMI were close to zero, regardless of lactation stage and data set (PC or MC). However, CH₄P and DMI were highly correlated, both across lactation stages and data sets. In conclusion, RT is unsuitable to be used as an indicator trait for either CH₄P or DMI. Our study failed to validate RT as a useful indicator trait for both CH₄P and DMI, but more studies with novel phenotypes can offer different approaches to select and incorporate important yet difficult to record traits into breeding goals and selection indexes.

Rumination time as a potential predictor of common diseases in high-productive Holstein dairy cows

We examined the hypothesis that rumination time (RT) could serve as a useful predictor of various common diseases of high producing dairy cows and hence improve herd management and animal wellbeing. We measured the changes in rumination time (RT) in the days before the recording of diseases (specifically: mastitis, reproductive system diseases, locomotor system issues, and gastroenteric diseases). We built predictive models to assess the association between RT and these diseases, using the former as the outcome variable, and to study the effects of the latter on the former. The average Pseudo-R2 of the fitted models was moderate to low, and this could be due to the fact that RT is influenced by other additional factors which have a greater effect than the predictors used here. Although remaining in a moderate-to-low range, the average Pseudo-R2 of the models regarding locomotion issues and gastroenteric diseases was higher than the others, suggesting the greater effect of these diseases on RT. The results are encouraging, but further work is needed if these models are to become useful predictors.