Search

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

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

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

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., (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/).)

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 (h 2 ) and genetic correlations (r a ) 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 (h 2  = 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 (r a  = 0.85 to 0.88). Genetic correlations of CM were moderate with MIRprob (r a  = 0.26 and 0.37) during 150 and 305 days in milk, respectively and low with LF (h 2  = 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., (© 2024 The Authors. Journal of Animal Breeding and Genetics published by John Wiley & Sons Ltd.)

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., (© 2023. The Author(s).)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of the American Society of Animal Science.)

Mass spectrometry (MS)-based techniques, including liquid chromatography coupling, shotgun lipidomics, MS imaging, and ion mobility, are widely used to analyze lipids. However, with enhanced separation capacity and an optimized chemical derivatization approach, comprehensive two-dimensional gas chromatography (GC×GC) can be a powerful tool to investigate some groups of small lipids in the framework of lipidomics. This study describes the optimization of a dedicated two-stage derivatization and extraction process to analyze different saturated and unsaturated fatty acids in plasma by two-dimensional gas chromatography-time-of-flight mass spectrometry (GC×GC-TOFMS) using a full factorial design. The optimized condition has a composite desirability of 0.9159. This optimized sample preparation and chromatographic condition were implemented to differentiate between positive (BT) and negative (UT) boar-tainted pigs based on fatty acid profiling in pig serum using GC×GC-TOFMS. A chemometric screening, including unsupervised (PCA, HCA) and supervised analysis (PLS-DA), as well as univariate analysis (volcano plot), was performed. The results suggested that the concentration of PUFA ω-6 and cholesterol derivatives were significantly increased in BT pigs, whereas SFA and PUFA ω-3 concentrations were increased in UT pigs. The metabolic pathway and quantitative enrichment analysis suggest the significant involvement of linolenic acid metabolism.

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., (© 2022 John Wiley & Sons Ltd.)

This research aims to develop a predictive model to discriminate milk produced from a cattle diet either based on grass or not using milk mid-infrared spectrometry and the month of testing (an indirect indicator of the feeding ration). The dataset contained 3,377,715 spectra collected between 2011 and 2021 from 2449 farms and 3 grazing traits defined following the month of testing. Records from 30% of the randomly selected farms were kept in the calibration set, and the remaining records were used to validate the models. Around 90% of the records were correctly discriminated. This accuracy is very good, as some records could be erroneously assigned. The probability of belonging to the GRASS modality allowed confirmation of the model's ability to detect the transition period even if the model was not trained on this data. Indeed, the probability increased from the spring to the summer and then decreased. The discrimination was mainly explained by the changes in the milk fat, mineral, and protein compositions. A hierarchical clustering from the averaged probability per farm and year highlighted 12 groups illustrating different management practices. The probability of belonging to the GRASS class could be used in a tool counting the number of grazing days.

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 ( r g from 0.46 to 0.85) and between fat yield and milk FA ( r g from 0.17 to 0.47). On the other hand, negative correlations were estimated between fat and protein contents ( r g from -0.22 to -0.59), between milk yield and milk FA ( r g from -0.22 to -0.62), and between protein yield and milk FA ( r g from -0.11 to -0.19). The selection for high fat content increases milk FA throughout lactation ( r g 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 ( r 2 from 0.09 to 0.38) and less biased predictions ( b 1 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 ( b 1 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.

Monitoring for mastitis on dairy farms is of particular importance, as it is one of the most prevalent bovine diseases. A commonly used indicator for mastitis monitoring is somatic cell count. A supplementary tool to predict mastitis risk may be mid-infrared (MIR) spectroscopy of milk. Because bovine health status can affect milk composition, this technique is already routinely used to determine standard milk components. The aim of the present study was to compare the performance of models to predict clinical mastitis based on MIR spectral data and/or somatic cell count score (SCS), and to explore differences of prediction accuracies for acute and chronic clinical mastitis diagnoses. Test-day data of the routine Austrian milk recording system and diagnosis data of its health monitoring, from 59,002 cows of the breeds Fleckvieh (dual purpose Simmental), Holstein Friesian and Brown Swiss, were used. Test-day records within 21 days before and 21 days after a mastitis diagnosis were defined as mastitis cases. Three different models (MIR, SCS, MIR + SCS) were compared, applying Partial Least Squares Discriminant Analysis. Results of external validation in the overall time window (-/+21 days) showed area under receiver operating characteristic curves (AUC) of 0.70 when based only on MIR, 0.72 when based only on SCS, and 0.76 when based on both. Considering as mastitis cases only the test-day records within 7 days after mastitis diagnosis, the corresponding areas under the curve were 0.77, 0.83 and 0.85. Hence, the model combining MIR spectral data and SCS was performing best. Mastitis probabilities derived from the prediction models are potentially valuable for routine mastitis monitoring for farmers, as well as for the genetic evaluation of the trait udder health.

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., (© 2022 Wiley-VCH GmbH.)

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

In recent years, cannabis and hemp-based products have become increasingly popular for recreational use, edibles, beverages, health care products, and medicines. The rapid detection and differentiation of phytocannabinoids is, therefore, essential to assess the potency and the therapeutic and nutritional values of cannabis cultivars. Here, we implemented SpiderMass technology for in vivo detection of cannabidiolic acid (CBDA) and ∆ 9 -tetrahydrocannabinolicacid (∆ 9 -THCA), and other endogenous organic plant compounds, to access distribution gradients within the plants and differentiate between cultivars. The SpiderMass system is composed of an IR-laser handheld microsampling probe connected to a mass spectrometer through a transfer tube. The analysis was performed on different plant organs from freshly cultivated cannabis plants in only a few seconds. SpiderMass analysis easily discriminated the two acid phytocannabinoid isomers via MS/MS, and the built statistical models differentiated between four cannabis cultivars. Different abundancies of the two acid phytocannabinoids were found along the plant as well as between different cultivars. Overall, these results introduce direct analysis by SpiderMass as a compelling analytical alternative for rapid hemp analysis.

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 (r 2  > 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., (© 2021 Wiley-VCH GmbH.)

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., (© 2021 Wiley-VCH GmbH.)

Measuring the mineral composition of milk is of major interest in the dairy sector. This study aims to develop and validate robust multi-breed and multi-country models predicting the major minerals through milk mid-infrared spectrometry using partial least square regressions. A total of 1281 samples coming from five countries were analyzed to obtain spectra and in ICP-AES to measure the mineral reference contents. Models were built from records coming from four countries ( n = 1181) and validated using records from the fifth country, Austria ( n = 100). The importance of including local samples was tested by integrating 30 Austrian samples in the model while validating with the remaining 70 samples. The best performances were achieved using this second set of models, confirming the need to cover the spectral variability of a country before making a prediction. Validation root mean square errors were 54.56, 63.60, 7.30, 59.87, and 152.89 mg/kg for Na, Ca, Mg, P, and K, respectively. The built models were applied on the Walloon milk recording large-scale spectral database, including 3,510,077. The large-scale predictions on this dairy herd improvement database provide new insight regarding the minerals' variability in the population, as well as the effect of parity, stage of lactation, breeds, and seasons.

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., (© 2021. The Author(s).)

Background: A robust proxy for estimating methane (CH 4 ) 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 CH 4 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 CH 4 were developed using a combined dataset including daily CH 4 measurements (n = 1089; g d -1 ) collected using the SF 6 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 SF 6 dataset predicted the expected pattern in CH 4 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: R 2 = 0.68 and standard error = 57 g CH 4 d -1 )., Conclusions: The models developed accounted for more of the observed variability in CH 4 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., (© 2020 Society of Chemical Industry.)

We predicted dry matter intake of dairy cows using parity, week of lactation, milk yield, milk mid-infrared (MIR) spectrum, and MIR-based predictions of bodyweight, fat, protein, lactose, and fatty acids content in milk. The dataset comprised 10,711 samples of 534 dairy cows with a geographical diversity (Australia, Canada, Denmark, and Ireland). We set up partial least square (PLS) regressions with different constructs and a one-hidden-layer artificial neural network (ANN) using the highest contribution variables. In the ANN, we replaced the spectra with their projections to the 25 first PLS factors explaining 99% of the spectral variability to reduce the model complexity. Cow-independent 10 × 10-fold cross-validation (CV) achieved the best performance with root mean square errors (RMSE CV ) of 3.27 ± 0.08 kg for the PLS regression and 3.25 ± 0.13 kg for ANN. Although the available data were significantly different, we also performed a country-independent validation (CIV) to measure the models' performance fairly. We found RMSE CIV varying from 3.73 to 6.03 kg for PLS and 3.69 to 5.08 kg for ANN. Ultimately, based on the country-independent validation, we discussed the developed models' performance with those achieved by the National Research Council's equation.

Knowing the body weight (BW) of a cow at a specific moment or measuring its changes through time is of interest for management purposes. The current work aimed to validate the feasibility of predicting BW using the day in milk, parity, milk yield, and milk mid-infrared (MIR) spectrum from a multiple-country dataset and reduce the number of predictors to limit the risk of over-fitting and potentially improve its accuracy. The BW modeling procedure involved feature selections and herd-independent validation in identifying the most interesting subsets of predictors and then external validation of the models. From 1849 records collected in 9 herds from 360 Holstein cows, the best performing models achieved a root mean square error (RMSE) for the herd-independent validation between 52 ± 2.34 kg to 56 ± 3.16 kg, including from 5 to 62 predictors. Among these models, three performed remarkably well in external validation using an independent dataset (N = 4067), resulting in RMSE ranging from 52 to 56 kg. The results suggest that multiple optimal BW predictive models coexist due to the high correlations between adjacent spectral points.

The use of abnormal milk mid-infrared (MIR) spectrum strongly affects prediction quality, even if the prediction equations used are accurate. So, this record must be detected after or before the prediction process to avoid erroneous spectral extrapolation or the use of poor-quality spectral data by dairy herd improvement (DHI) organizations. For financial or practical reasons, adapting the quality protocol currently used to improve the accuracy of fat and protein contents is unfeasible. This study proposed three different statistical methods that would be easy to implement by DHI organizations to solve this issue: the deletion of 1% of the extreme high and low predictive values (M1), the deletion of records based on the Global-H (GH) distance (M2), and the deletion of records based on the absolute fat residual value (M3). Additionally, the combinations of these three methods were investigated. A total of 346,818 milk samples were analyzed by MIR spectrometry to predict the contents of fat, protein, and fatty acids. Then, the same traits were also predicted externally using their corresponded standardized MIR spectra. The interest in cleaning procedures was assessed by estimating the root mean square differences (RMSDs) between those internal and external predicted phenotypes. All methods allowed for a decrease in the RMSD, with a gain ranging from 0.32% to 41.39%. Based on the obtained results, the "M1 and M2" combination should be preferred to be more parsimonious in the data loss, as it had the higher ratio of RMSD gain to data loss. This method deleted the records based on the 2% extreme predictions and a GH threshold set at 5. However, to ensure the lowest RMSD, the "M2 or M3" combination, considering a GH threshold of 5 and an absolute fat residual difference set at 0.30 g/dL of milk, was the most relevant. Both combinations involved M2 confirming the high interest of calculating the GH distance for all samples to predict. However, if it is impossible to estimate the GH distance due to a lack of relevant information to compute this statistical parameter, the obtained results recommended the use of M1 combined with M3. The limitation used in M3 must be adapted by the DHI, as this will depend on the spectral data and the equation used. The methodology proposed in this study can be generalized for other MIR-based phenotypes.

Boar taint is an unpleasant odor in male pig meat, mainly caused by androstenone, skatole, and indole, which are deposited in the fat tissue. Piglet castration is the most common practice to prevent boar taint. However, castration is likely to be banished in a few years due to animal welfare concerns. Alternatives to castration, such as genetic selection, have been assessed. Androstenone and skatole have moderate to high heritability, which makes it feasible to select against these compounds. This review presents the latest results obtained on genetic selection against boar taint, on correlation with other traits, on differences in breeds, and on candidate genes related to boar taint. QTLs for androstenone and skatole have been reported mainly on chromosomes 6, 7, and 14. These chromosomes were reported to contain genes responsible for synthesis and degradation of androstenone and skatole. A myriad of work has been done to find markers or genes that can be used to select animals with lower boar taint. The selection against boar taint could decrease performance of some reproduction traits. However, a favorable response on production traits has been observed by selecting against boar taint. Selection results have shown that it is possible to reduce boar taint in few generations. In addition, modifications in diet and environment conditions could be associated with genetic selection to reduce boar taint. Nevertheless, costs to measure and select against boar taint should be rewarded with incentives from the market; otherwise, it would be difficult to implement genetic selection.

Phenotypes related to feed efficiency were predicted from records easily acquired by breeding organizations. A total of 461,036 and 354,148 records were collected from the first and second parity Holstein cows. Equations were applied to the milk mid-infrared spectra to predict the main milk components and coupled with animal characteristics to predict the body weight (pBW). Dry matter intake (pDMI) was predicted from pBW using the National Research Council (NRC) equation. The consumption index (pIC) was estimated from pDMI and fat, and protein corrected milk. All traits were modeled using single trait test-day models. Descriptive statistics were within the expected range. Milk yield, pDMI, and pBW were phenotypically positively related (r ranged from 0.08 to 0.64). As expected, pIC was phenotypically negatively correlated with milk yield (-0.77 and -0.80 for the first and second lactation) and slightly positively correlated with pBW (0.16 and 0.07 for the first and second lactation). Later, parity cows seemed to have a better feed efficiency as they had a lower pIC. Although the prediction accuracy was moderate, the observed behaviors of studied traits by year, stage of lactation, and parity were in agreement with the literature. Moreover, as a genetic component was highlighted (heritability around 0.18), it would be interesting to realize a genetic evaluation of these traits and compare the obtained breeding values with the ones estimated for sires having daughters with reference feed efficiency records.

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

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 C 1 , C 2 , C 3 , C 4 and C 5 . 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 (S 1 ), non-genotyped (S 2 ) and inbred non-genotyped (S 3 ) animals. Potential bias was assessed by regressing true dominance deviations on predicted values. Matrices accounting for inbreeding (C 3 , C 4 and C 5 ) best fit. Accuracies were on average 0.77, 0.40 and 0.46 in S 1 , S 2 and S 3 , respectively. In addition, C 3 , C 4 and C 5 scenarios have shown better accuracies than C 1 and C 2 , 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 (C 5 )., (© 2019 Blackwell Verlag GmbH.)

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 ( h 2 ) 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., Competing Interests: The authors declare that they have no conflict of interest., (Copyright: © 2019 Marwa Grayaa et al.)

The aim of this study was to evaluate the effect of UGT1A1 polymorphisms on Raltegravir (RAL) and its metabolite RAL-glucuronide trough plasma concentrations ([RAL]plasma and [RAL-glu]plasma) and on the metabolic ratio (MR): [RAL-glu]plasma/[RAL]plasma. UGT1A1 genotyping was performed on 96 patients. 44% (n = 42) were homozygous UGT1A1*1/*1 while 50% (n = 48) and 6% (n = 6) were UGT1A1*28 and UGT1A1*36 carriers, respectively. The median concentration and interquartile range (IQR) of [RAL]plasma were 88.5 ng/ml (41.0-236), 168 ng/ml (85.8-318) and 92.5 ng/ml (36.4-316) for UGT1A1*1/*1, UGT1A1*28 and UGT1A1*36 carriers, respectively. Only the difference between UGT1A1*1/*1 and *28 carriers was statistically significant (p = 0.022). The median MR (IQR) were 5.8 (3-10), 2.9 (1.6-5.3) and 3.2 (1.7-5.9) for UGT1A1*1/*1, UGT1A1*28 and UGT1A1*36 carriers, respectively. Only the difference between UGT1A1*1/*1 and *28 carriers was statistically significant (p = 0.004) with an allele-dependent effect: UGT1A1*28 homozygous having lower MR than heterozygous carriers who show lower MR compared to *1/*1. Except for the sensation of fatigue, this PK effect did not correlate with clinical adverse events or biological abnormalities. In Conclusion, we demonstrate that UGT1A1*28 polymorphism has a significant impact on RAL metabolism: UGT1A1*28 carriers being characterized by higher [RAL]plasma and lower MR.

Electronic cigarette use has raised concern worldwide regarding potential health risks and its position in tobacco cessation strategies. As part of any toxicity assessment, the chemical characterization of e-liquids and their related vapors are among fundamental data to be determined. Considering the lack of available reference methods, we developed and validated several analytical procedures in order to conduct a multicomponent analysis of six e-liquid refills and their resultant vapor emissions (generated by a smoking machine), and compared them with tobacco smoke. We combined several techniques including gas-chromatography, high and ultra-performance liquid chromatography and inductively coupled plasma with mass spectrometry or ultraviolet and flame ionization detection in order to identify the main e-liquid constituents (propylene glycol, glycerol and nicotine), as well as multiple potentially harmful components (trace elements, polycyclic aromatic hydrocarbons (PAHs), pesticides and carbonyl compounds). Regarding propylene glycol, glycerol and nicotine concentrations, the six tested e-liquids comply with the advertised composition and contain only traces of pollutants. Noticeable lower concentrations of trace elements (≤3.4 pg/mL puff), pesticides (

Through centuries of both natural and artificial selection, a variety of local cattle populations arose with highly specific phenotypes. However, the intensification and expansion of scale in animal production systems led to the predominance of a few highly productive cattle breeds. The loss of local populations is often considered irreversible and with them specific qualities and rare variants could be lost as well. Over these last years, the interest in these local breeds has increased again leading to increasing efforts to conserve these breeds or even revive lost populations, e.g. through the use of crosses with similar breeds. However, the remaining populations are expected to contain crossbred individuals resulting from introgressions. They are likely to carry exogenous genes that affect the breed's authenticity on a genomic level. Using the revived Campine breed as a case study, 289 individuals registered as purebreds were genotyped on the Illumina BovineSNP50. In addition, genomic information on the Illumina BovineHD and Illumina BovineSNP50 of ten breeds was available to assess the current population structure, genetic diversity, and introgression with phenotypically similar and/or historically related breeds. Introgression with Holstein and beef cattle genotypes was limited to only a few farms. While the current population shows a substantial amount of within-breed variation, the majority of genotypes can be separated from other breeds in the study, supporting the re-establishment of the Campine breed. The majority of the population is genetically close to the Deep Red (NL), Improved Red (NL) and Eastern Belgium Red and White (BE) cattle, breeds known for their historical ties to the Campine breed. This would support an open herdbook policy, thereby increasing the population size and consequently providing a more secure future for the breed.

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., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Reducing the frequency of milk recording would help reduce the costs of official milk recording. However, this approach could also negatively affect the accuracy of predicting daily yields. This problem has been investigated in numerous studies. In addition, published equations take into account milking intervals (MI), and these are often not available and/or are unreliable in practice. The first objective of this study was to propose models in which the MI was replaced by a combination of data easily recorded by dairy farmers. The second objective was to further investigate the fatty acids (FA) present in milk. Equations to predict daily yield from AM or PM data were based on a calibration database containing 79,971 records related to 51 traits [milk yield (expected AM, expected PM, and expected daily); fat content (expected AM, expected PM, and expected daily); fat yield (expected AM, expected PM, and expected daily; g/day); levels of seven different FAs or FA groups (expected AM, expected PM, and expected daily; g/dL milk), and the corresponding FA yields for these seven FA types/groups (expected AM, expected PM, and expected daily; g/day)]. These equations were validated using two distinct external datasets. The results obtained from the proposed models were compared to previously published results for models which included a MI effect. The corresponding correlation values ranged from 96.4% to 97.6% when the daily yields were estimated from the AM milkings and ranged from 96.9% to 98.3% when the daily yields were estimated from the PM milkings. The simplicity of these proposed models should facilitate their use by breeding and milk recording organizations.

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.

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.

In Greece, seven native horse breeds have been identified so far. Among these, the Skyros pony is outstanding through having a distinct phenotype. In the present study, the aim was to assess genetic diversity in this breed, by using different types of genetic loci and available genealogical information. Its relationships with the other Greek, as well as foreign, domestic breeds were also investigated. Through microsatellite and pedigree analysis it appeared that the Skyros presented a similar level of genetic diversity to the other European breeds. Nevertheless, comparisons between DNA-based and pedigree-based results revealed that a loss of genetic diversity had probably already occurred before the beginning of breed registration. Tests indicated the possible existence of a recent bottleneck in two of the three main herds of Skyros pony. Nonetheless, relatively high levels of heterozygosity and Polymorphism Information Content indicated sufficient residual genetic variability, probably useful in planning future strategies for breed conservation. Three other Greek breeds were also analyzed. A comparison of these with domestic breeds elsewhere, revealed the closest relationships to be with the Middle Eastern types, whereas the Skyros itself remained isolated, without any close relationship, whatsoever.