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

Feed efficiency and enteric methane emissions indices are inconsistent with the outcomes of the rumen microbiome composition

The correlation between enteric methane emissions (eME) and feed efficiency (FE) in cattle is linked to the anaerobic fermentation of feedstuffs that occurs in the rumen. Several mathematical indices have been developed to predict feed efficiency and identify low methane emitters in herds. To investigate this, the current study aimed to evaluate the rumen microbial composition in the same group of animals ranked according to six different indices (three indices for FE and three for eME). Thirty-three heifers were ranked into three groups, each consisting of 11 animals, based on FE (feed conversion efficiency - FCE, residual weight gain - RG, and residual feed intake - RFI) and eME indices (production, yield, and intensity). Rumen fluids were collected using a stomach tube and analyzed using 16S rRNA and 18S rRNA, targeting rumen bacteria, archaea, and protozoa. The sequencing analysis revealed that the presence of unique microbial species in the rumen varies across animals ranked by the FE and eME indices. The High RG group harbored 17 unique prokaryotic taxa, while the High FCE group contained only seven. Significant differences existed in the microbial profiles of the animals based on the FE and eME indices. For instance, Raoultibacter was more abundant in the Intermediate RFI group but less so in the Intermediate RG and Intermediate FCE groups. The abundance of Entodinium was higher while Diplodinium was lower in the High FCE group, in contrast to the High RG and High RFI groups. Methanobrevibacter exhibited similar abundances across eME indices. However, the heifers did not demonstrate the same production, yield, and intensity of eME. The present findings underscore the importance of standardizing the FE and eME indices. This standardization is crucial for ensuring consistent and reliable assessments of the composition and function of the rumen microbiome across different herds.

Review: Improving residual feed intake modelling in the context of nutritional- and genetic studies for dairy cattle

The residual feed intake (RFI) model has recently gained popularity for ranking dairy cows for feed efficiency. The RFI model ranks the cows based on their expected feed intake compared to the observed feed intake, where a negative phenotype (eating less than expected) is favourable. Yet interpreting the biological implications of the regression coefficients derived from RFI models has proven challenging. In addition, multitrait modelling of RFI has been proposed as an alternative to the least square RFI in nutrition and genetic studies. To solve the challenge with the biological interpretation of RFI regression coefficients and suggest ways to improve the modelling of RFI, an interdisciplinary effort was required between nutritionists and geneticists. Therefore, this paper aimed to explore the challenges with the traditional least square RFI model and propose solutions to improve the modelling of RFI. In the traditional least square RFI model, one set of fixed effects is used to solve systematic effects (e.g., seasonal effects and age at calving) for traits with different means and variances. Thereby, measurement and model fitting errors can accumulate in the phenotype, resulting in undesirable effects. A multivariate RFI model will likely reduce this problem, as trait-specific fixed effects are used. In addition, regression coefficients for DM intake on milk energy tend to have more biologically meaningful estimates in multitrait RFI models, which indicates a confounding effect between the fixed effects and regression coefficients in the least square RFI model. However, defining precise expectations for regression coefficients from RFI models or sourcing for accurate feed norm coefficients seems difficult, especially if the coefficients are applied to a wide cattle population with varying diets or management systems, for example. To improve multitrait modelling of RFI, we suggest improving the modelling of changes in energy status. Furthermore, a novel method to derive the energy density of the diet and individual digestive efficiency is proposed. Digestive efficiency is defined as the part of the efficiency associated with digestive processes, which primarily reflects the conversion from gross energy to metabolisable energy. We show the model was insensitive to prior values of energy density in feed and that there was individual variation in digestive efficiency. The proposed method needs further development and validation. In summary, using multitrait RFI can improve the accuracy of the ranking of dairy cows' feed efficiency, consequently improving economic and environmental sustainability on dairy farms.

Milk transcriptome biomarker identification to enhance feed efficiency and reduce nutritional costs in dairy ewes

In recent years, rising prices for high-quality protein-based feeds have significantly increased nutrition costs. Consequently, investigating strategies to reduce these expenses and improve feed efficiency (FE) have become increasingly important for the dairy sheep industry. This research investigates the impact of nutritional protein restriction (NPR) during prepuberty and FE on the milk transcriptome of dairy Assaf ewes (sampled during the first lactation). To this end, we first compared transcriptomic differences between NPR and control ewes. Subsequently, we evaluated gene expression differences between ewes with divergent FE, using feed conversion ratio (FCR), residual feed intake (RFI), and consensus classifications of high- and low-FE animals for both indices. Lastly, we assess milk gene expression as a predictor of FE phenotype using random forest. No effect was found for the prepubertal NPR on milk performance or FE. Moreover, at the milk transcriptome level, only one gene, HBB, was differentially expressed between the NPR (n = 14) and the control group (n = 14). Further, the transcriptomic analysis between divergent FE sheep revealed 114 differentially expressed genes (DEGs) for RFI index (high-FERFI = 10 vs low-FERFI = 10), 244 for FCR (high-FEFCR = 10 vs low-FEFCR = 10), and 1 016 DEGs between divergent consensus ewes for both indices (high-FEconsensus = 8 vs low-FEconsensus = 8). These results underscore the critical role of selected FE indices for RNA-Seq analyses, revealing that consensus divergent animals for both indices maximise differences in transcriptomic responses. Genes overexpressed in high-FEconsensus ewes were associated with milk production and mammary gland development, while low-FEconsensus genes were linked to higher metabolic expenditure for tissue organisation and repair. The best prediction accuracy for FE phenotype using random forest was obtained for a set of 44 genes consistently differentially expressed across lactations, with Spearman correlations of 0.37 and 0.22 for FCR and RFI, respectively. These findings provide insights into potential sustainability strategies for dairy sheep, highlighting the utility of transcriptomic markers as FE proxies.

Clustering of countries based on dairy productivity characteristics of Holstein cattle for breeding material selection

Background and Aim

The aim of any breeding process is to create a herd based on certain parameters that reflect an ideal animal vision. Targeted herding involves selecting the source of breeding material to be imported from another country. Therefore, there is a problem in selecting a breeding material importer to rapidly form a uterine canopy with the required properties. The purpose of this study was to evaluate a set of predictive milk productivity traits in Holstein cattle across countries.

Materials and Methods

This research was based on records of 819,358 recorded animals from 28 countries born after January 1, 2018, from open databases. We used the Euclidean metric to construct dendrograms characterizing the similarity of countries according to the complex milk productivity traits of the daughters of bulls. The Ward method was used to minimize intracluster variance when forming clusters and constructing the corresponding diagrams. Principal component analysis was used to reduce dimensionality and eliminate the effect of multicollinearity. The principal components were selected using the Kaiser-Harris criteria.

Results

A ranking of multidimensional complex milk productivity traits in different countries over the past 5 years was performed. A group of leading countries led by the USA was established according to the studied indicators, and the possible reasons for such a division into groups were described.

Conclusion

The pressure of purposeful artificial selection prevails in comparison with the pressure of natural selection concerning milk productivity traits in a certain group of countries, which allows specialists to choose suppliers when buying breeding animals and materials. The findings are based solely on data from recorded animals, which may not represent the entire breed population within each country, especially in regions where record-keeping may be inconsistent. It is expected that further studies will include regional data from large enterprises not part of Interbull, with mandatory verification and validation. An important element of such work is seen as the ability to compare the milk productivity of populations from different countries using a different scale, as well as studying the differentiation of countries by other selection traits of dairy.

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.

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 for feed intake and body weight in dairy cattle using high-throughput 3-dimensional cameras in Danish commercial farms

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

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

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

Genomic evaluation of feed efficiency in US Holstein heifers

There is growing interest in improving feed efficiency traits in dairy cattle. The objectives of this study were to estimate the genetic parameters of residual feed intake (RFI) and its component traits [dry matter intake (DMI), metabolic body weight (MBW), and average daily gain (ADG)] in Holstein heifers, and to develop a system for genomic evaluation for RFI in Holstein dairy calves. The RFI data were collected from 6,563 growing Holstein heifers (initial body weight = 261 ± 52 kg; initial age = 266 ± 42 d) for 70 d, across 182 trials conducted between 2014 and 2022 at the STgenetics Ohio Heifer Center (South Charleston, OH) as part of the EcoFeed program, which aims to improve feed efficiency by genetic selection. The RFI was estimated as the difference between a heifer's actual feed intake and expected feed intake, which was determined by regression of DMI against midpoint MBW, age, and ADG across each trial. A total of 61,283 SNPs were used in genomic analyses. Animals with phenotypes and genotypes were used as training population, and 4 groups of prediction population, each with 2,000 animals, were selected from a pool of Holstein animals with genotypes, based on their relationship with the training population. All traits were analyzed using univariate animal model in DMU version 6 software. Pedigree information and genomic information were used to specify genetic relationships to estimate the variance components and genomic estimated breeding values (GEBV), respectively. Breeding values of the prediction population were estimated by using the 2-step approach: deriving the prediction equation of GEBV from the training population for estimation of GEBV of prediction population with only genotypes. Reliability of breeding values was obtained by approximation based on partitioning a function of the accuracy of training population GEBV and magnitudes of genomic relationships between individuals in the training and prediction population. Heifers had DMI (mean ± SD) of 8.11 ± 1.59 kg over the trial period, with growth rate of 1.08 ± 0.25 kg/d. The heritability estimates (mean ± SE) of RFI, MBW, DMI, and growth rate were 0.24 ± 0.02, 0.23 ± 0.02, 0.27 ± 0.02, and 0.19 ± 0.02, respectively. The range of genomic predicted transmitted abilities (gPTA) of the training population (-0.94 to 0.75) was higher compared with the range of gPTA (-0.82 to 0.73) of different groups of prediction population. Average reliability of breeding values from the training population was 58%, and that of prediction population was 39%. The genomic prediction of RFI provides new tools to select for feed efficiency of heifers. Future research should be directed to find the relationship between RFI of heifers and cows, to select individuals based on their lifetime production efficiencies.

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

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

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

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

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

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

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

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

Simultaneous Bayesian estimation of genetic parameters for curves of weight, feed intake, and residual feed intake in beef cattle

Rates of gain and feed efficiency are important traits in most breeding programs for growing farm animals. The rate of gain (GAIN) is usually expressed over a certain age period and feed efficiency is often expressed as residual feed intake (RFI), defined as observed feed intake (FI) minus expected feed intake based on live weight (WGT) and GAIN. However, the basic traits recorded are always WGT and FI and other traits are derived from these basic records. The aim of this study was to develop a procedure for simultaneous analysis of the basic records and then derive linear traits related to feed efficiency without retorting to any approximation. A bivariate longitudinal random regression model was employed on 13,791 individual longitudinal records of WGT and FI from 2,827 bulls of six different beef breeds tested for their own performance in the period from 7 to 13 mo of age. Genetic and permanent environmental covariance functions for curves of WGT and FI were estimated using Gibbs sampling. Genetic and permanent covariance functions for curves of GAIN were estimated from the first derivative of the function for WGT and finally the covariance functions were extended to curves for RFI, based on the conditional distribution of FI given WGT and GAIN. Furthermore, the covariance functions were extended to include GAIN and RFI defined over different periods of the performance test. These periods included the whole test period as normally used when predicting breeding values for GAIN and RFI for beef bulls. Based on the presented method, breeding values and genetic parameters for derived traits such as GAIN and RFI defined longitudinally or integrated over (parts of) of the test period can be obtained from a joint analysis of the basic records. The resulting covariance functions for WGT, FI, GAIN, and RFI are usually singular but the method presented here does not suffer from the estimation problems associated with defining these traits individually before the genetic analysis. All the results are thus estimated simultaneously, and the set of parameters is consistent.

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.