Search

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]

Quantifying the level of linkage disequilibrium (LD), non‐random association of alleles at two or more loci, is important to determine the number of markers needed for genomic selection. The aims of this study were to evaluate the extent of LD in Dual‐Purpose Belgian Blue (DPBB) and to compare the level of LD in DPBB with that of Walloon Holstein. Data of 28,427 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA), of 639 DPBB and 398 Holstein bulls were used. The level of LD between pairwise SNPs separated by up to 10 Mb was evaluated, separately for each breed, using the squared correlation of the alleles at two loci. The analysis of molecular variance showed that the percentage of variation within populations (85.48%) was higher than between populations (14.52%). However, permutation tests showed a significant genetic differentiation between the two studied populations (p <.01). The average LD found between adjacent SNP pairs in DPBB (0.16 (SD = 0.22)) was generally lower than in Holstein (0.23 (SD = 0.27)). The proportion of SNPs in useful LD (r2 > 0.30) within a genomic distance of ≤0.10 Mb between SNPs was 18.58% and 28.23% in DPBB and Holstein bulls, respectively. In both breeds, the effective population size decreased over generations; however, the decline was greater in DPBB than that in Holstein. Based on results, it can be concluded that at least 68,000 SNPs are needed for implementing genomic selection in DPBB cattle with enough accuracy. [ABSTRACT FROM AUTHOR]

Assignment of individual cattle to a specific breed can often not rely on pedigree information. This is especially the case for local breeds for which the development of genomic assignment tools is required to allow individuals of unknown origin to be included to their herd books. A breed assignment model can be based on two specific stages: (a) the selection of breed‐informative markers and (b) the assignment of individuals to a breed with a classification method. However, the performance of combination of methods used in these two stages has been rarely studied until now. In this study, the combination of 16 different SNP panels with four classification methods was developed on 562 reference genotypes from 12 cattle breeds. Based on their performances, best models were validated on three local breeds of interest. In cross‐validation, 14 models had a global cross‐validation accuracy higher than 90%, with a maximum of 98.22%. In validation, best models used 7,153 or 2,005 SNPs, based on a partial least squares‐discriminant analysis (PLS‐DA) and assigned individuals to breeds based on nearest shrunken centroids. The average validation sensitivity of the first two best models for the three local breeds of interest were 98.33% and 97.5%. Moreover, results reported in this study suggest that further studies should consider the PLS‐DA method when selecting breed‐informative SNPs. [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]

Recent publications indicate that single‐step models are suitable to estimate breeding values, dominance deviations and total genetic values with acceptable quality. Additive single‐step methods implicitly extend known number of allele information from genotyped to non‐genotyped animals. This theory is well derived in an additive setting. It was recently shown, at least empirically, that this basic strategy can be extended to dominance with reasonable prediction quality. Our study addressed two additional issues. It illustrated the theoretical basis for extension and validated genomic predictions to dominance based on single‐step genomic best linear unbiased prediction theory. This development was then extended to include inbreeding into dominance relationships, which is a currently not yet solved issue. Different parametrizations of dominance relationship matrices were proposed. Five dominance single‐step inverse matrices were tested and described as C1, C2, C3, C4 and C5. Genotypes were simulated for a real pig population (n = 11,943 animals). In order to avoid any confounding issues with additive effects, pseudo‐records including only dominance deviations and residuals were simulated. SNP effects of heterozygous genotypes were summed up to generate true dominance deviations. We added random noise to those values and used them as phenotypes. Accuracy was defined as correlation between true and predicted dominance deviations. We conducted five replicates and estimated accuracies in three sets: between all (S1), non‐genotyped (S2) and inbred non‐genotyped (S3) animals. Potential bias was assessed by regressing true dominance deviations on predicted values. Matrices accounting for inbreeding (C3, C4 and C5) best fit. Accuracies were on average 0.77, 0.40 and 0.46 in S1, S2 and S3, respectively. In addition, C3, C4 and C5 scenarios have shown better accuracies than C1 and C2, and dominance deviations were less biased. Better matrix compatibility (accuracy and bias) was observed by re‐scaling diagonal elements to 1 minus the inbreeding coefficient (C5). [ABSTRACT FROM AUTHOR]

Genetic improvement of udder health in dairy cows is of high relevance as mastitis is one of the most prevalent diseases. Since it is known that the heritability of mastitis is low and direct data on mastitis cases are often not available in large numbers, auxiliary traits, such as somatic cell count (SCC) are used for the genetic evaluation of udder health. In previous studies, models to predict clinical mastitis based on mid‐infrared (MIR) spectral data and a somatic cell count‐derived score (SCS) were developed. Those models can provide a probability of mastitis for each cow at every test‐day, which is potentially useful as an additional auxiliary trait for the genetic evaluation of udder health. Furthermore, MIR spectral data were used to estimate contents of lactoferrin, a glycoprotein positively associated with immune response. The present study aimed to estimate heritabilities (h2) and genetic correlations (ra) for clinical mastitis diagnosis (CM), SCS, MIR‐predicted mastitis probability (MIRprob), MIR + SCS‐predicted mastitis probability (MIRSCSprob) and lactoferrin estimates (LF). Data for this study were collected within the routine milk recording and health monitoring system of Austria from 2014 to 2021 and included records of approximately 54,000 Fleckvieh cows. Analyses were performed in two datasets, including test‐day records from 5 to 150 or 5 to 305 days in milk. Prediction models were applied to obtain MIR‐ and SCS‐based phenotypes (MIRprob, MIRSCSprob, LF). To estimate heritabilities and genetic correlations bivariate linear animal models were applied for all traits. A lactation model was used for CM, defined as a binary trait, and a test‐day model for all other continuous traits. In addition to the random animal genetic effect, the fixed effects year‐season of calving and parity‐age at calving and the random permanent environmental effect were considered in all models. For CM the random herd‐year effect, for continuous traits the random herd‐test day effect and the covariate days in milk (linear and quadratic) were additionally fitted. The obtained genetic parameters were similar in both datasets. The heritability found for CM was expectedly low (h2 = 0.02). For SCS and MIRSCSprob, heritability estimates ranged from 0.23 to 0.25, and for MIRprob and LF from 0.15 to 0.17. CM was highly correlated with SCS and MIRSCSprob (ra = 0.85 to 0.88). Genetic correlations of CM were moderate with MIRprob (ra = 0.26 and 0.37) during 150 and 305 days in milk, respectively and low with LF (h2 = 0.10 and 0.11). However, basic selection index calculations indicate that the added value of the new MIR‐predicted phenotypes is limited for genetic evaluation of udder health. [ABSTRACT FROM AUTHOR]