Search

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., (The Authors. Published by Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

Direct measurements of methane (CH 4 ) from individual animals are difficult and expensive. Predictions based on proxies for CH 4 are a viable alternative. Most prediction models are based on multiple linear regressions (MLR) and predictor variables that are not routinely available in commercial farms, such as dry matter intake (DMI) and diet composition. The use of machine learning (ML) algorithms to predict CH 4 emissions from across-country heterogeneous data sets has not been reported. The objectives were to compare performances of ML ensemble algorithm random forest (RF) and MLR models in predicting CH 4 emissions from proxies in dairy cows, and assess effects of imputing missing data points on prediction accuracy. Data on CH 4 emissions and proxies for CH 4 from 20 herds were provided by 10 countries. The integrated data set contained 43,519 records from 3,483 cows, with 18.7% missing data points imputed using k-nearest neighbor imputation. Three data sets were created, 3k (no missing records), 21k (missing DMI imputed from milk, fat, protein, body weight), and 41k (missing DMI, milk fat, and protein records imputed). These data sets were used to test scenarios (with or without DMI, imputed vs. nonimputed DMI, milk fat, and protein), and prediction models (RF vs. MLR). Model predictive ability was evaluated within and between herds through 10-fold cross-validation. Prediction accuracy was measured as correlation between observed and predicted CH 4 , root mean squared error (RMSE) and mean normalized discounted cumulative gain (NDCG). Inclusion of DMI in the model improved within and between-herd prediction accuracy to 0.77 (RMSE = 23.3%) and 0.58 (RMSE = 31.9%) in RF and to 0.50 (RMSE = 0.327) and 0.13 (RMSE = 42.71) in MLR, respectively than when DMI was not included in the predictive model. When missing DMI records were imputed, within and between-herd accuracy increased to 0.84 (RMSE = 18.5%) and 0.63 (RMSE = 29.9%), respectively. In all scenarios, RF models out-performed MLR models. Results suggest routinely measured variables from dairy farms can be used in developing globally robust prediction models for CH 4 if coupled with state-of-the-art techniques for imputation and advanced ML algorithms for predictive modeling., (© 2022, The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).)

This study reports on the exploration of temporal relationships between milk mid-infrared predicted biomarkers and lameness events. Lameness in dairy cows is an issue that can vary greatly in severity and is of concern for both producers and consumers. Metabolic disorders are often associated with lameness. However, lameness can arise weeks or even months after the metabolic disorder, making the detection of causality difficult. We already use mid-infrared technology to predict major milk components, such as fat or protein, during routine milk recording and for milk payment. It was recently shown that this technology can also be used to predict novel biomarkers linked to metabolic disorders in cows, such as oleic acid (18:1 cis-9), β-hydroxybutyrate, acetone, and citrate in milk. We used these novel biomarkers as proxies for metabolic issues. Other studies have explored the possibility of using mid-infrared spectra to predict metabolic diseases and found it (potentially) usable for indicating classes of metabolic problems. We wanted to explore the possible relationship between mid-infrared-based metabolites and lameness over the course of lactation. In total, data were recorded from 6,292 cows on 161 farms in Austria. Lameness data were recorded between March 2014 and March 2015 and consisted of 37,555 records. Mid-infrared data were recorded between July and December 2014 and consisted of 9,152 records. Our approach consisted of fitting preadjustments to the data using fixed effects, computing pair-wise correlations, and finally applying polynomial smoothing of the correlations for a given biomarker at a certain month in lactation and the lameness events scored on severity scale from sound or non-lame (lameness score of 1) to severely lame (lameness score of 5) throughout the lactation. The final correlations between biomarkers and lameness scores were significant, but not high. However, for the results of the present study, we should not look at the correlations in terms of absolute values, but rather as indicators of a relationship through time. When doing so, we can see that metabolic problems occurring in mo 1 and 3 seem more linked to long-term effects on hoof and leg health than those in mo 2. However, the quantity (only 1 pair-wise correlation exceeded 1,000 observations) and the quality (due to limited data, no separation according to more metabolic-related diseases could be done) of the data should be improved., (The Authors. Published by FASS Inc. and Elsevier Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).)

The objectives of this study were to estimate genetic parameters and identify genomic regions associated with milk urea concentration (MU) in Dual‐Purpose Belgian Blue (DPBB) cows. The data were 29,693 test‐day records of milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FP), protein percentage (PP) and MU collected between 2014 and 2020 on 2498 first parity cows (16,935 test‐day records) and 1939 second‐parity cows (12,758 test‐day records) from 49 herds in the Walloon Region of Belgium. Data of 28,266 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), on 1699 animals (639 males and 1060 females) were used. Random regression test‐day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method using a single chain of 100,000 iterations after a burn‐in period of 20,000. SNP solutions were estimated using a single‐step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by windows of 25 consecutive SNPs (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. The mean (SD) of MU was 22.89 (10.07) and 22.35 (10.07) mg/dl for first and second parity, respectively. The mean (SD) heritability estimates for daily MU were 0.18 (0.01) and 0.22 (0.02), for first and second parity, respectively. The mean (SD) genetic correlations between daily MU and MY, FY, PY, FP and PP were −0.05 (0.09), −0.07 (0.11), −0.03 (0.13), −0.05 (0.08) and −0.03 (0.11) for first parity, respectively. The corresponding values estimated for second parity were 0.02 (0.10), −0.02 (0.09), 0.02 (0.08), −0.08 (0.06) and −0.05 (0.05). The genome‐wide association analyses identified three genomic regions (BTA2, BTA3 and BTA13) associated with MU. The identified regions showed contrasting results between parities and among different stages within each parity. This suggests that different groups of candidate genes underlie the phenotypic expression of MU between parities and among different lactation stages within a parity. The results of this study can be used for future implementation and use of genomic evaluation to reduce MU in DPBB cows. [ABSTRACT FROM AUTHOR]

We investigated the use of different Legendre polynomial orders to estimate genetic parameters for milk production and fatty acid (FA) traits in the first lactation Walloon Holstein cows. The data set comprised 302,684 test‐day records of milk yield, fat and protein contents, and FAs generated by mid‐infrared (MIR) spectroscopy, C16:0 (palmitic acid), C18:1 cis‐9 (oleic acid), LCFAs (long‐chain FAs), SFAs (saturated FAs) and UFAs (unsaturated FAs) were studied. The models included random regression coefficients for herd‐year of calving (h), additive genetic (a) and permanent environment (p) effects. The selection of the best random regression model (RRM) was based on the deviance information criterion (DIC), and genetic parameters were estimated via a Bayesian approach. For all analysed random effects, DIC values decreased as the order of the Legendre polynomials increased. Best‐fit models had fifth‐order (degree 4) for the p effect and ranged from second‐ to fifth‐order (degree 1–4) for the a and h effects (LEGhap: LEG555 for milk yield and protein content; LEG335 for fat content and SFA; LEG545 for C16:0 and UFA; and LEG535 for C18:1 cis‐9 and LCFA). Based on the best‐fit models, an effect of overcorrection was observed in early lactation (5–35 days in milk [DIM]). On the contrary, third‐order (LEG333; degree 2) models showed flat residual trajectories throughout lactation. In general, the estimates of genetic variance tended to increase over DIM, for all traits. Heritabilities for milk production traits ranged from 0.11 to 0.58. Milk FA heritabilities ranged from low‐to‐high magnitude (0.03–0.56). High Spearman correlations (>0.90 for all bulls and >0.97 for top 100) were found among breeding values for 155 and 305 DIM between the best RRM and LEG333 model. Therefore, third‐order Legendre polynomials seem to be most parsimonious and sufficient to describe milk production and FA traits in Walloon Holstein cows. [ABSTRACT FROM AUTHOR]

BACKGROUND: A robust proxy for estimating methane (CH4) emissions of individual dairy cows would be valuable especially for selective breeding. This study aimed to improve the robustness and accuracy of prediction models that estimate daily CH4 emissions from milk Fourier transform mid‐infrared (FT‐MIR) spectra by (i) increasing the reference dataset and (ii) adjusting for routinely recorded phenotypic information. Prediction equations for CH4 were developed using a combined dataset including daily CH4 measurements (n = 1089; g d−1) collected using the SF6 tracer technique (n = 513) and measurements using respiration chambers (RC, n = 576). Furthermore, in addition to the milk FT‐MIR spectra, the variables of milk yield (MY) on the test day, parity (P) and breed (B) of cows were included in the regression analysis as explanatory variables. RESULTS: Models developed based on a combined RC and SF6 dataset predicted the expected pattern in CH4 values (in g d−1) during a lactation cycle, namely an increase during the first weeks after calving followed by a gradual decrease until the end of lactation. The model including MY, P and B information provided the best prediction results (cross‐validation statistics: R2 = 0.68 and standard error = 57 g CH4 d−1). CONCLUSIONS: The models developed accounted for more of the observed variability in CH4 emissions than previously developed models and thus were considered more robust. This approach is suitable for large‐scale studies (e.g. animal genetic evaluation) where robustness is paramount for accurate predictions across a range of animal conditions. © 2020 Society of Chemical Industry [ABSTRACT FROM AUTHOR]

Genetic parameters were estimated for first lactation survival defined as a binary trait (alive or dead to second calving) and the curve shape traits of milk yield, fat and protein percentages using information from 25 981 primiparous Tunisian Holsteins. For each trait, shape curves (i.e. peak lactation, persistency), level of production adjusted to 305 days in milk (DIMs) for total milk yield (TMY), and average fat (TF %) and protein (TP %) percentages were defined. Variance components were estimated with a linear random regression model under three bivariate animal models. Production traits were modelled by fixed herd × test-day (TD) interaction effects, fixed classes of 25 DIMs × age of calving × season of calving interaction effects, fixed classes of pregnancy, random environment effects and random additive genetic effects. Survival was modelled by fixed herd × year of calving interaction effects and age of calving × season of calving interaction effects, random permanent environment effects, and random additive genetic effects. Heritability (h2) estimates were 0.03 (±0.01) for survival and 0.23 (±0.01), 0.31 (±0.01) and 0.31 (±0.01) for TMY, TF % and TP %, respectively. Genetic correlations between survival and TMY, TF % and TP % were 0.26 (±0.08), -0.24 (±0.06) and -0.13 (±0.06), respectively. Genetic correlations between survival and persistency for fat and protein percentages were -0.35 (±0.09) and -0.19 (±0.09), respectively. Cows that had higher persistencies for fat and protein percentages were more likely not to survive. [ABSTRACT FROM AUTHOR]

Copyright of Biotechnologie, Agronomie, Societe et Environnement is the property of Les Presses Agronomiques de Gembloux and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Milk citrate is regarded as an early biomarker of negative energy balance in dairy cows during early lactation and serves as a suitable candidate phenotype for genomic selection due to its wide availability across a large number of cows through milk mid-infrared spectra prediction. However, its genetic background is not well known. Therefore, the objectives of this study were to (1) analyze the genetic parameters of milk citrate; (2) identify genomic regions associated with milk citrate; and (3) analyze the functional annotation of candidate genes and quantitative trait loci (QTL) related to milk citrate in Walloon Holstein cows. In total, 134,517 test-day milk-citrate phenotypes (mmol/L) collected within the first 50 d in milk on 52,198 Holstein cows were used. These milk-citrate phenotypes, predicted by milk mid-infrared spectra, were divided into 3 traits according to the first (citrate1), second (citrate2), and third to fifth parity (citrate3+). Genomic information for 566,170 SNPs was available for 4,479 animals. A multiple-trait repeatability model was used to estimate genetic parameters. A single-step GWAS was used to identify candidate genes for citrate and post-GWAS analysis was done to investigate the relationship and function of the identified candidate genes. The heritabilities estimated for citrate1, citrate2, and citrate3+ were 0.40, 0.37, and 0.35, respectively. The genetic correlations among the 3 traits ranged from 0.98 to 0.99. The genomic correlations among the 3 traits were also close to 1.00 across the genomic regions (1 Mb) in the whole genome, which means that citrate can be considered as a single trait in the first 5 parities. In total, 603 significant SNPs located on 3 genomic regions (chromosome 7, 68.569–68.575 Mb; chromosome 14, 0.15–1.90 Mb; and chromosome 20, 54.00–64.28 Mb), were identified to be associated with milk citrate. We identified 89 candidate genes including GPT , ANKH , PPP1R16A , and 32 QTL reported in the literature related to the identified significant SNPs. These identified QTL were mainly reported associated with milk fatty acids and metabolic diseases in dairy cows. This study suggests that milk citrate in Holstein cows is highly heritable and has the potential to be used as an early proxy for the negative energy balance of Holstein cows in a breeding objective. [ABSTRACT FROM AUTHOR]