Search

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]

In the context of genetic research, but even more for routine genetic evaluations, a recurrent problem is the difficulty to obtain measurements for certain required phenotypes. This can be caused by limited facilities, difficulties in sampling, expensive analyses, limited staff relative to number of animals, and many other reasons. It may therefore be necessary to select which animals to phenotype while maintaining the phenotyping quality and avoiding bias or omitting some animals of interest. This implies developing strategies to do the phenotyping in a non-biased way using a priori available minimum knowledge. An example could be predicting post weaning growth of cattle based on earlier information as birth weight and weaning weight to optimize the use of feed intake or methane measurement facilities. In pigs, similar situations might occur for growth or emissions. However, another important example is boar taint (BT) which shows that not only sequential weighing, but also partial measurements of BT may require attention. BT is defined as urinary and fecal and associated with 2 major compounds, androstenone (AND) and skatole (SKA). The classical method of odor discrimination, "Human Nose Score" (HNS), is notoriously unreliable. Even if odor is known to be complex, adding AND and SKA to HNS would help reveal better the boar taint condition (i.e., phenotype). However, the measurements of AND and SKA or other compounds in fat are expensive. Therefore, it would be important to optimize these measurements by selecting relevant fat samples taken at slaughterhouses. We tested the following optimization strategy. The objective was to obtain a good distribution of a posteriori measured AND and SKA in the selected samples given initial asymmetric populations with positive HNS (3 to 10% declared positives). In a first step, the strategy pooled 2 available HNS for 997 across 22 contemporary groups (CG) and added limited chemical measurements (218/997) available for older samples. The data were obtained from a boar progeny testing scheme; therefore, CG were linked through connecting sires from artificial insemination having sons in different CG. A multi-trait BLUP model was used adjusting for metabolic weight using a 4-generation pedigree. A total of 15% of the animals in the last 6 contemporary groups were selected using first positive HNS (17/300), then keeping the highest and lowest predictions and maintaining diversity by not repeating boars from the same sire-family (i.e., when a family was already present, the next animal was taken). A posteriori chemical analyses were performed on selected samples to compare predicted and observed classifications in tainted and untainted. Research showed an accuracy of 0.73 (Table 1), leading to 58% tainted and 42% untainted compared with the expected ratio. Improvements of this strategy are currently in development, including using genomics. [ABSTRACT FROM AUTHOR]

Evaluating genetically heat stress (HS) of cattle is still a major challenge, also because of the need to have access to HS sensitive phenotypes relating their reaction to temperature and humidity as described by the temperature-humidity index (THI). A solution could be to add the information of data from sensors. However, activity (ACT) and rumination (RUM) sensors reflect no "production" [i.e., milk yields (MY), milk composition, but also intake or growth] related aspects which is suboptimal. Therefore, the objectives of this study were thus to assess the possible usefulness for adding sensors data for HS genetic evaluation based on a dairy situation. But the context described can be easily enlarged and include beef cattle and pigs. Strategy was illustrated using daily SenseHub collar records that were obtained from October 2019 to July 2022 in six dairy herds in Wallonia in Southern Belgium. A total of 453 Walloon Holstein cows were followed during this period for activity time. Meteorological data and a HS sensitive "production" trait (here MY) were also obtained from 2015 to 2022 for 1,740 cows from the same herds. In total 32,154 MY, 130,867 ACT and 130,848 RUM records were available. The thresholds at which the different traits start to be affected by HS were estimated at a THI of 63, 64 and 66, respectively, for MY, ACT and RUM. A three-trait random regression reaction norm model was fitted with these different thresholds. Estimated heritability values were 0.19 ± 0.03, 0.14 ± 0.06 and 0.21 ± 0.05 respectively for MY, ACT and RUM using AIREML. Ratios of variance associated to THI and constant genetic effect were 0.002, 0.02 and 0.04 indicating a stronger reaction by RUM but also ACT than MY to THI. However, it is uncertain if the reactions of these sensor traits are reflected in the production traits. Estimated genetic correlations of regressions on THI for ACT and RUM showed values of 0.51 ± 0.70 and 0.08 ± 0.54 with MY and 0.77 ± 0.69 between ACT and RUM. The first eigenvector of the correlation matrix explained 65% of the total variability and standardized coefficients of this eigenvector were 0.42, 0.69 and 0.59 for MY, ACT and RUM. A total of 50% of the variance associated with the regression on THI for MY could be explained by the regressions on THI for ACT and RUM. Adding the constant genetic effects for ACT and RUM only marginally increased to 51% the ability to explain the effect of THI on MY. Based on these results, reactions of different indicator traits, here sensor data, may reflect a different reaction than production and therefore economically important traits. Additional research, also beyond dairy cattle, should continue to address this topic that is very important for the validation of sensor data detecting HS. [ABSTRACT FROM AUTHOR]

With the rapid development of animal phenomics and deep phenotyping, we can get thousands of traditional but also molecular phenotypes per individual. However, there is still a lack of exploration regarding how to handle this huge amount of data in the context of animal breeding, presenting a challenge that we are likely to encounter more and more in the future. This study aimed to (1) explore the use of the Mega-scale linear mixed model (MegaLMM), a factor model-based approach, able to simultaneously estimate (co)variance components and genetic parameters in the context of thousands of milk traits, hereafter called thousand-trait (TT) models; (2) compare the phenotype values and genomic breeding values (u) predictions for focal traits (i.e., traits that are targeted for prediction, compared with secondary traits that are helping to evaluate), from single-trait (ST) and TT models, respectively; (3) propose a new approximate method of estimated genomic breeding values (U) prediction with TT models and MegaLMM. 3,421 milk mid-infrared (MIR) spectra wavepoints (called secondary traits) and 3 focal traits [average fat percent (Fat), average methane (CH4), and average somatic cell score (SCS)] collected on 3,302 first-parity Holstein cows were used. The 3,421 milk MIR wavepoints traits were composed of 311 wavepoints in 11 classes (months in lactation). Genotyping information of 564,439 SNP was available for all animals and was used to calculate the genomic relationship matrix. The MegaLMM was implemented in the framework of the Bayesian sparse factor model and solved through Gibbs sampling (Markov chain Monte Carlo). The heritabilities of the studied 3,421 milk MIR wavepoints gradually increased and then decreased in units of 311 wavepoints throughout the lactation. The genetic and phenotypic correlations between the first 311 wavepoints and the other 3,110 wavepoints were low. The accuracies of phenotype predictions from the ST model were lower than those from the TT model for Fat (0.51 vs. 0.93), CH4 (0.30 vs. 0.86), and SCS (0.14 vs. 0.33). The same trend was observed for the accuracies of u predictions: Fat (0.59 vs. 0.86), CH4 (0.47 vs. 0.78), and SCS (0.39 vs. 0.59). The average correlation between U predicted from the TT model and the new approximate method was 0.90. The new approximate method used for estimating U in MegaLMM will enhance the suitability of MegaLMM for applications in animal breeding. This study conducted an initial investigation into the application of thousands of traits in animal breeding and showed that the TT model is beneficial for the prediction of focal traits (phenotype and breeding values), especially for difficult-to-measure traits (e.g., CH4)., (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/).)