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.

Genetic analysis of milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation.

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.

Exploring a Bayesian sparse factor model-based strategy for the genetic analysis of thousands of MIR-spectra traits for animal breeding.

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).

Genetic parameters for mid-infrared-spectroscopy-predicted mastitis phenotypes and related traits.

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.

Defining valid breeding goals for animal breeds.

BACKGROUND: The objective of any valid breeding program is to increase the suitability of a breed for its future purposes. The approach most often followed in animal breeding for optimizing breeding goals assumes that the sole desire of the owners is profit maximization. As this assumption is often violated, a generalized approach is needed that does not rely on this assumption. RESULTS: The generalized approach is based on the niche concept. The niche of a breed is a set of environments in which a small population of the breed would have a positive population growth rate. Its growth rate depends on demand from prospective consumers and supply from producers. The approach involves defining the niche that is envisaged for the breed and identifying the trait optima that maximize the breed's adaptation to its envisaged niche within the set of permissible breeding goals. The set of permissible breeding goals is the set of all potential breeding goals that are compatible with animal welfare and could be reached within the planning horizon of the breeding program. In general, the breed's adaptation depends on the satisfaction of the producers with the animals and on the satisfaction of the consumers with the products produced by the animals. When consumers buy live animals, then the breed needs to adapt to both the environments provided by the producers, and the environments provided by the consumers. The profit function is replaced by a more general adaptedness function that measures the breed's adaptation to its envisaged niche. CONCLUSIONS: The proposed approach coincides with the traditional approach if the producers have the sole desire to maximize their income, and if consumer preferences are well reflected by the product prices. If these assumptions are not met, then the traditional approach to breeding goal optimization is unlikely to result in a valid breeding goal. Using the example of companion breeds, this paper shows that the proposed approach has the potential to fill the gap.

Integrating heterogeneous across-country data for proxy-based random forest prediction of enteric methane in dairy cattle.

Direct measurements of methane (CH4) from individual animals are difficult and expensive. Predictions based on proxies for CH4 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 CH4 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 CH4 emissions from proxies in dairy cows, and assess effects of imputing missing data points on prediction accuracy. Data on CH4 emissions and proxies for CH4 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 CH4, 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 CH4 if coupled with state-of-the-art techniques for imputation and advanced ML algorithms for predictive modeling.

Genomic prediction and genetic correlations estimated for milk production and fatty acid traits in Walloon Holstein cattle using random regression models.

The aims of this study were to: (1) estimate genetic correlation for milk production traits (milk, fat and protein yields and fat and protein contents) and fatty acids (FA: C16:0, C18:1 cis-9, LCFA, SFA, and UFA) over days in milk, (2) investigate the performance of genomic predictions using single-step GBLUP (ssGBLUP) based on random regression models (RRM), and (3) identify the optimal scaling and weighting factors to be used in the construction of the H matrix. A total of 302 684 test-day records of 63.875 first lactation Walloon Holstein cows were used. Positive genetic correlations were found between milk yield and fat and protein yield (rg from 0.46 to 0.85) and between fat yield and milk FA (rg from 0.17 to 0.47). On the other hand, negative correlations were estimated between fat and protein contents (rg from -0.22 to -0.59), between milk yield and milk FA (rg from -0.22 to -0.62), and between protein yield and milk FA (rg from -0.11 to -0.19). The selection for high fat content increases milk FA throughout lactation (rg from 0.61 to 0.98). The test-day ssGBLUP approach showed considerably higher prediction reliability than the parent average for all milk production and FA traits, even when no scaling and weighting factors were used in the H matrix. The highest validation reliabilities (r2 from 0.09 to 0.38) and less biased predictions (b1 from 0.76 to 0.92) were obtained using the optimal parameters (i.e., ω = 0.7 and α = 0.6) for the genomic evaluation of milk production traits. For milk FA, the optimal parameters were ω = 0.6 and α = 0.6. However, biased predictions were still observed (b1 from 0.32 to 0.81). The findings suggest that using ssGBLUP based on RRM is feasible for the genomic prediction of daily milk production and FA traits in Walloon Holstein dairy cattle.

The pattern of linkage disequilibrium in Dual-Purpose Belgian Blue cattle.

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.

Genome-wide association for milk urea concentration in Dual-Purpose Belgian Blue cows.

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.

Development of a genomic tool for breed assignment by comparison of different classification models: Application to three local cattle breeds.

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.

Random regression test-day models to describe milk production and fatty acid traits in first lactation Walloon Holstein cows.

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.

Investigation on heroin and cocaine quality in Luxembourg.

BACKGROUND: Heroin and cocaine are among the most dangerous illicit drugs available and their presence on the market is increasing. These facts have led to the investigation of the quality of heroin and cocaine samples seized in Luxembourg by police and customs but also collected at the national supervised drug consumption facilities. METHODS: Samples obtained from 2019 to 2020 were analyzed to determine their composition and content using GC-MS, HPLC-UV and LC-Q-ToF. The statistical evaluation of concentration changes depending on the source of collection is based on an ANOVA single factor test and a two-tailed t test. RESULTS: Results showed important differences between seizure and collection sources. For both drugs, customs samples had significantly higher concentrations than police samples and the latter had significantly higher concentrations than samples from drug consumption facilities, whereas for heroin two cutting steps were identified, for cocaine samples only one appears to occur on the local market. Indeed, cocaine samples seized by police consisted of a mixture of low and high concentration samples. CONCLUSION: The results show that extensive adulteration with pharmacological active and inactive compounds takes place at local levels, which, however, are different for heroin and cocaine. This knowledge on variability of quality of drugs should be considered in the elaboration of drug and harm prevention strategies.

Improving robustness and accuracy of predicted daily methane emissions of dairy cows using milk mid-infrared spectra.

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.

Short communication: Investigation of the temporal relationships between milk mid-infrared predicted biomarkers and lameness events in later lactation.

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.

Additional considerations to the use of single-step genomic predictions in a dominance setting.

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 ).

Challenges and priorities for modelling livestock health and pathogens in the context of climate change.

Climate change has the potential to impair livestock health, with consequences for animal welfare, productivity, greenhouse gas emissions, and human livelihoods and health. Modelling has an important role in assessing the impacts of climate change on livestock systems and the efficacy of potential adaptation strategies, to support decision making for more efficient, resilient and sustainable production. However, a coherent set of challenges and research priorities for modelling livestock health and pathogens under climate change has not previously been available. To identify such challenges and priorities, researchers from across Europe were engaged in a horizon-scanning study, involving workshop and questionnaire based exercises and focussed literature reviews. Eighteen key challenges were identified and grouped into six categories based on subject-specific and capacity building requirements. Across a number of challenges, the need for inventories relating model types to different applications (e.g. the pathogen species, region, scale of focus and purpose to which they can be applied) was identified, in order to identify gaps in capability in relation to the impacts of climate change on animal health. The need for collaboration and learning across disciplines was highlighted in several challenges, e.g. to better understand and model complex ecological interactions between pathogens, vectors, wildlife hosts and livestock in the context of climate change. Collaboration between socio-economic and biophysical disciplines was seen as important for better engagement with stakeholders and for improved modelling of the costs and benefits of poor livestock health. The need for more comprehensive validation of empirical relationships, for harmonising terminology and measurements, and for building capacity for under-researched nations, systems and health problems indicated the importance of joined up approaches across nations. The challenges and priorities identified can help focus the development of modelling capacity and future research structures in this vital field. Well-funded networks capable of managing the long-term development of shared resources are required in order to create a cohesive modelling community equipped to tackle the complex challenges of climate change.

Unified method to integrate and blend several, potentially related, sources of information for genetic evaluation.

BACKGROUND: A condition to predict unbiased estimated breeding values by best linear unbiased prediction is to use simultaneously all available data. However, this condition is not often fully met. For example, in dairy cattle, internal (i.e. local) populations lead to evaluations based only on internal records while widely used foreign sires have been selected using internally unavailable external records. In such cases, internal genetic evaluations may be less accurate and biased. Because external records are unavailable, methods were developed to combine external information that summarizes these records, i.e. external estimated breeding values and associated reliabilities, with internal records to improve accuracy of internal genetic evaluations. Two issues of these methods concern double-counting of contributions due to relationships and due to records. These issues could be worse if external information came from several evaluations, at least partially based on the same records, and combined into a single internal evaluation. Based on a Bayesian approach, the aim of this research was to develop a unified method to integrate and blend simultaneously several sources of information into an internal genetic evaluation by avoiding double-counting of contributions due to relationships and due to records. RESULTS: This research resulted in equations that integrate and blend simultaneously several sources of information and avoid double-counting of contributions due to relationships and due to records. The performance of the developed equations was evaluated using simulated and real datasets. The results showed that the developed equations integrated and blended several sources of information well into a genetic evaluation. The developed equations also avoided double-counting of contributions due to relationships and due to records. Furthermore, because all available external sources of information were correctly propagated, relatives of external animals benefited from the integrated information and, therefore, more reliable estimated breeding values were obtained. CONCLUSIONS: The proposed unified method integrated and blended several sources of information well into a genetic evaluation by avoiding double-counting of contributions due to relationships and due to records. The unified method can also be extended to other types of situations such as single-step genomic or multi-trait evaluations, combining information across different traits.

Development and validation of a classification model for boar taint detection in pork fat samples.

This study aims to characterize a complete volatile organic compound profile of pork neck fat for boar taint prediction. The objectives are to identify specific compounds related to boar taint and to develop a classification model. In addition to the well-known androstenone, skatole and indole, 10 other features were found to be discriminant according to untargeted volatolomic analyses were conducted on 129 samples using HS-SPME-GC×GC-TOFMS. To select the odor-positive samples among the 129 analyzed, the selection was made by combining human nose evaluations with the skatole and androstenone concentrations determined using UHPLC-MS/MS. A comparison of the data of the two populations was performed and a statistical model analysis was built on 70 samples out of the total of 129 samples fully positive or fully negative through these two orthogonal methods for tainted prediction. Then, the model was applied to the 59 remaining samples. Finally, 7 samples were classified as tainted.

Inversion of a part of the numerator relationship matrix using pedigree information.

BACKGROUND: In recent theoretical developments, the information available (e.g. genotypes) divides the original population into two groups: animals with this information (selected animals) and animals without this information (excluded animals). These developments require inversion of the part of the pedigree-based numerator relationship matrix that describes the genetic covariance between selected animals (A22). Our main objective was to propose and evaluate methodology that takes advantage of any potential sparsity in the inverse of A22 in order to reduce the computing time required for its inversion. This potential sparsity is brought out by searching the pedigree for dependencies between the selected animals. Jointly, we expected distant ancestors to provide relationship ties that increase the density of matrix A22 but that their effect on A22-1 might be minor. This hypothesis was also tested. METHODS: The inverse of A22 can be computed from the inverse of the triangular factor (T-1) obtained by Cholesky root-free decomposition of A22. We propose an algorithm that sets up the sparsity pattern of T-1 using pedigree information. This algorithm provides positions of the elements of T-1 worth to be computed (i.e. different from zero). A recursive computation of A22-1 is then achieved with or without information on the sparsity pattern and time required for each computation was recorded. For three numbers of selected animals (4000; 8000 and 12 000), A22 was computed using different pedigree extractions and the closeness of the resulting A22-1 to the inverse computed using the fully extracted pedigree was measured by an appropriate norm. RESULTS: The use of prior information on the sparsity of T-1 decreased the computing time for inversion by a factor of 1.73 on average. Computational issues and practical uses of the different algorithms were discussed. Cases involving more than 12 000 selected animals were considered. Inclusion of 10 generations was determined to be sufficient when computing A22. CONCLUSIONS: Depending on the size and structure of the selected sub-population, gains in time to compute A22-1 are possible and these gains may increase as the number of selected animals increases. Given the sequential nature of most computational steps, the proposed algorithm can benefit from optimization and may be convenient for genomic evaluations.

The use of a genomic relationship matrix for breed assignment of cattle breeds: comparison and combination with a machine learning method.

To develop a breed assignment model, three main steps are generally followed: 1) The selection of breed informative single nucleotide polymorphism (SNP); 2) The training of a model, based on a reference population, that allows to classify animals to their breed of origin; and 3) The validation of the developed model on external animals i.e., that were not used in previous steps. However, there is no consensus in the literature about which methodology to follow for the first step, nor about the number of SNP to be selected. This can raise many questions when developing the model and lead to the use of sophisticated methodologies for selecting SNP (e.g., with iterative algorithms, partitions of SNP, or combination of several methods). Therefore, it may be of interest to avoid the first step by the use of all the available SNP. For this purpose, we propose the use of a genomic relationship matrix (GRM), combined or not with a machine learning method, for breed assignment. We compared it with a previously developed model based on selected informative SNP. Four methodologies were investigated: 1) The PLS_NSC methodology: selection of SNP based on a partial least square-discriminant analysis (PLS-DA) and breed assignment by classification based on the nearest shrunken centroids (NSC) method; 2) Breed assignment based on the highest mean relatedness of an animal to the reference populations of each breed (referred to mean_GRM); 3) Breed assignment based on the highest SD of the relatedness of an animal to the reference populations of each breed (referred to SD_GRM) and 4) The GRM_SVM methodology: the use of means and SD of the relatedness defined in mean_GRM and SD_GRM methodologies combined with the linear support vector machine (SVM), a machine learning method used for classification. Regarding mean global accuracies, results showed that the use of mean_GRM or GRM_SVM was not significantly different (Bonferroni corrected P > 0.0083) than the model based on a reduced SNP panel (PLS_NSC). Moreover, the mean_GRM and GRM_SVM methodology were more efficient than PLS_NSC as it was faster to compute. Therefore, it is possible to bypass the selection of SNP and, by the use of a GRM, to develop an efficient breed assignment model. In routine, we recommend the use of GRM_SVM over mean_GRM as it gave a slightly increased global accuracy, which can help endangered breeds to be maintained. The script to execute the different methodologies can be accessed on: https://github.com/hwilmot675/Breed_assignment.

Breed assignment models generally rely on three main steps: 1) Selection of markers that allow to distinguish the breeds under study; 2) Development of a classification model that assigns each animal to its breed of origin; and 3) Validation of the developed model with new animals, to verify that the developed model is not overfitted. The first step often raises several questions about the methodology to select the best markers or about the number of markers to select. That is why it can be interesting to avoid this first step and to use an appropriate methodology that performs similarly without the need for single nucleotide polymorphism (SNP) selection. In this study, we developed different methodologies based on the genomic relationship matrix (GRM), combined or not with a machine learning method, to assign animals to their breed of origin. The results showed that the model based on a GRM combined with a machine learning method showed equivalent percentage of correct assignment to a previously developed model relying on SNP selection while being substantially faster to compute. It is therefore possible to assign animals to their breed by the use of a GRM and to bypass the first step of selection of SNP.