Genetic components of milk Fourier-transform infrared spectra used to predict breeding values for milk composition and quality traits in dairy goats

The usual practice today is that milk component phenotypes are predicted using Fourier-transform infrared (FTIR) spectra and they are then, together with pedigree information, used in BLUP for calculation of individual estimated breeding values. Here, this is referred to as the indirect prediction (IP) approach. An alternative approach—a direct prediction (DP) method—is proposed, where genetic analyses are directly conducted on the milk FTIR spectral variables. Breeding values of all derived milk traits (protein, fat, fatty acid composition, and coagulation properties, among others) can then be predicted as traits correlated only to the genetic information of the spectra. For the DP, no need exists to predict the phenotypes before calculating breeding values for each of the traits—the genetic analysis is done once for the spectra, and is applicable to all traits derived from the spectra. The aim of the study was to compare the effects of DP and IP of milk composition and quality traits on prediction error variance (PEV) and genetic gain. A data set containing 27,927 milk FTIR spectral observations and milk composition phenotypes (fat, lactose, and protein) belonging to 14,869 goats of 271 herds was used for training and evaluating models. Partial least squares regression was used for calibrating prediction models for fat, protein, and lactose percentages. Restricted maximum likelihood was used to estimate variance components of the spectral variables after principal components analysis was applied to reduce the spectral dimension. Estimated breeding values were predicted for fat, lactose, and protein percentages using DP and IP methods. The DP approach reduced the mean PEV by 3.73, 4.07, and 7.04% for fat, lactose, and protein percentages, respectively, compared with the IP method. Given the reduction in PEV, relative genetic gains were 2.99, 2.78, and 4.85% for fat, lactose, and protein percentages, respectively. We concluded that more accurate estimated breeding values could be found using genetic components of milk FTIR spectra compared with single-trait animal model analyses on phenotypes predicted from the spectra separately. The potential and application is not only limited to milk FTIR spectra, but could also be extended to any spectroscopy techniques implemented in other species and for other traits.     

Prediction accuracy of direct and indirect approaches,and their relationships with prediction ability of calibration models

Milk infrared spectra are routinely used for phenotyping traits of interest through links developed between the traits and spectra. Predicted individual traits are then used in genetic analyses for estimated breeding value (EBV) or for phenotypic predictions using a single-trait mixed model; this approach is referred to as indirect prediction (IP). An alternative approach [direct prediction (DP)] is a direct genetic analysis of (a reduced dimension of) the spectra using a multitrait model to predict multivariate EBV of the spectral components and, ultimately, also to predict the univariate EBV or phenotype for the traits of interest. We simulated 3 traits under different genetic (low: 0.10 to high: 0.90) and residual (zero to high: ±0.90) correlation scenarios between the 3 traits and assumed the first trait is a linear combination of the other 2 traits. The aim was to compare the IP and DP approaches for predictions of EBV and phenotypes under the different correlation scenarios. We also evaluated relationships between performances of the 2 approaches and the accuracy of calibration equations. Moreover, the effect of using different regression coefficients estimated from simulated phenotypes (β_p), true breeding values (β_g), and residuals (β_r) on performance of the 2 approaches were evaluated. The simulated data contained 2,100 parents (100 sires and 2,000 cows) and 8,000 offspring (4 offspring per cow). Of the 8,000 observations, 2,000 were randomly selected and used to develop links between the first and the other 2 traits using partial least square (PLS) regression analysis. The different PLS regression coefficients, such as β_p, β_g, and β_r, were used in subsequent predictions following the IP and DP approaches. We used BLUP analyses for the remaining 6,000 observations using the true (co)variance components that had been used for the simulation. Accuracy of prediction (of EBV and phenotype) was calculated as a correlation between predicted and true values from the simulations. The results showed that accuracies of EBV prediction were higher in the DP than in the IP approach. The reverse was true for accuracy of phenotypic prediction when using β_p but not when using β_g and β_r, where accuracy of phenotypic prediction in the DP was slightly higher than in the IP approach. Within the DP approach, accuracies of EBV when using β_g were higher than when using β_p only at the low genetic correlation scenario. However, we found no differences in EBV prediction accuracy between the β_p and β_g in the IP approach. Accuracy of the calibration models increased with an increase in genetic and residual correlations between the traits. Performance of both approaches increased with an increase in accuracy of the calibration models. In conclusion, the DP approach is a good strategy for EBV prediction but not for phenotypic prediction, where the classical PLS regression-based equations or the IP approach provided better results.     

The effects of heat stress in Italian Holstein dairy cattle

The data set for this study comprised 1,488,474 test-day records for milk, fat, and protein yields and fat and protein percentages from 191,012 first-, second-, and third-parity Holstein cows from 484 farms. Data were collected from 2001 through 2007 and merged with meteorological data from 35 weather stations. A linear model (M1) was used to estimate the effects of the temperature-humidity index (THI) on production traits. Least squares means from M1 were used to detect the THI thresholds for milk production in all parities by using a 2-phase linear regression procedure (M2). A multiple-trait repeatability test-model (M3) was used to estimate variance components for all traits and a dummy regression variable (t) was defined to estimate the production decline caused by heat stress. Additionally, the estimated variance components and M3 were used to estimate traditional and heat-tolerance breeding values (estimated breeding values, EBV) for milk yield and protein percentages at parity 1. An analysis of data (M2) indicated that the daily THI at which milk production started to decline for the 3 parities and traits ranged from 65 to 76. These THI values can be achieved with different temperature/humidity combinations with a range of temperatures from 21 to 36°C and relative humidity values from 5 to 95%. The highest negative effect of THI was observed 4 d before test day over the 3 parities for all traits. The negative effect of THI on production traits indicates that first-parity cows are less sensitive to heat stress than multiparous cows. Over the parities, the general additive genetic variance decreased for protein content and increased for milk yield and fat and protein yield. Additive genetic variance for heat tolerance showed an increase from the first to third parity for milk, protein, and fat yield, and for protein percentage. Genetic correlations between general and heat stress effects were all unfavorable (from −0.24 to −0.56). Three EBV per trait were calculated for each cow and bull (traditional EBV, traditional EBV estimated with the inclusion of THI covariate effect, and heat tolerance EBV) and the rankings of EBV for 283 bulls born after 1985 with at least 50 daughters were compared. When THI was included in the model, the ranking for 17 and 32 bulls changed for milk yield and protein percentage, respectively. The heat tolerance genetic component is not negligible, suggesting that heat tolerance selection should be included in the selection objectives.     

Estimation of genetic parameters for cheese-making traits in Spanish Churra sheep

The global production of sheep milk is growing, and the main industrial use of sheep milk is cheese making. The Spanish Churra sheep breed is one of the most important native dairy breeds in Spain. The present study aimed to estimate genetic parameters for a wide range of traits influencing the cheese-making ability of Churra sheep milk. Using a total of 1,049 Churra ewes, we studied the following cheese-making traits: 4 traits related to milk coagulation properties (rennet coagulation time, curd-firming time, and curd firmness at 30 and 60 min after addition of rennet), 2 traits related to cheese yield (individual laboratory cheese yield and individual laboratory dried curd yield), and 3 traits measuring curd firmness over time (maximum curd firmness, time to attain maximum curd firmness, and syneresis). In addition, a list of milk traits, including the native pH of the milk and several milk production and composition traits (milk yield; the fat, protein, and dried extract percentages; and the somatic cell count), were also analyzed for the studied animals. After discarding the noncoagulating samples (only 3.7%), data of 1,010 ewes were analyzed with multiple-trait animal models by using the restricted maximum likelihood method to estimate (co)variance components, heritabilities, and genetic correlations. In general, the heritability estimates were low to moderate, ranging from 0.08 (for the individual laboratory dried curd yield trait) to 0.42 (for the fat percentage trait). High genetic correlations were found within pairs of related traits (i.e., 0.93 between fat and dried extract percentages, ?0.93 between the log of the curd-firming time and curd firmness at 30 min, 0.70 between individual laboratory cheese yield and individual laboratory dried curd yield, and ?0.94 between time to attain maximum curd firmness and syneresis). Considering all the information provided here, we suggest that in addition to the current consideration of the protein percentage trait for improving cheese yield traits, the inclusion of the pH of milk as a measured trait in the Churra dairy breeding program would represent an efficient strategy for improving the cheese-making ability of milk from this breed.     

The use of Fourier-transform infrared spectroscopy to predict cheese yield and nutrient recovery or whey loss traits from unprocessed bovine milk samples

Cheese yield is an important technological trait in the dairy industry in many countries. The aim of this study was to evaluate the effectiveness of Fourier-transform infrared (FTIR) spectral analysis of fresh unprocessed milk samples for predicting cheese yield and nutrient recovery traits. A total of 1,264 model cheeses were obtained from 1,500-mL milk samples collected from individual Brown Swiss cows. Individual measurements of 7 new cheese yield-related traits were obtained from the laboratory cheese-making procedure, including the fresh cheese yield, total solid cheese yield, and the water retained in curd, all as a percentage of the processed milk, and nutrient recovery (fat, protein, total solids, and energy) in the curd as a percentage of the same nutrient contained in the milk. All individual milk samples were analyzed using a MilkoScan FT6000 over the spectral range from 5,000 to 900 wavenumber × cm⁻¹. Two spectral acquisitions were carried out for each sample and the results were averaged before data analysis. Different chemometric models were fitted and compared with the aim of improving the accuracy of the calibration equations for predicting these traits. The most accurate predictions were obtained for total solid cheese yield and fresh cheese yield, which exhibited coefficients of determination between the predicted and measured values in cross-validation (1-VR) of 0.95 and 0.83, respectively. A less favorable result was obtained for water retained in curd (1-VR = 0.65). Promising results were obtained for recovered protein (1-VR = 0.81), total solids (1-VR = 0.86), and energy (1-VR = 0.76), whereas recovered fat exhibited a low accuracy (1-VR = 0.41). As FTIR spectroscopy is a rapid, cheap, high-throughput technique that is already used to collect standard milk recording data, these FTIR calibrations for cheese yield and nutrient recovery highlight additional potential applications of the technique in the dairy industry, especially for monitoring cheese-making processes and milk payment systems. In addition, the prediction models can be used to provide breeding organizations with information on new phenotypes for cheese yield and milk nutrient recovery, potentially allowing these traits to be enhanced through selection.     

Cheese yield,cheesemaking efficiency,and daily production of 6 breeds of goats

Little is known about the complex process of cheesemaking at the individual level of dairy goats because of the difficulties of producing a high number of model cheeses. The objectives of this work were (1) to study the cheesemaking ability of goat milk; (2) to investigate the variability of cheesemaking-related traits among different farms; (3) to assess the effects of stage of lactation and parity; and (4) to compare 6 breeds of goat (Saanen and Camosciata delle Alpi for the Alpine type; Murciano-Granadina, Maltese, Sarda and Sarda Primitiva for the Mediterranean type) for their cheesemaking ability. For each goat (n = 560) we studied (1) 8 milk quality traits (fat, protein, total solids, casein, lactose, pH, somatic cell score, and bacterial count); (2) 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in curd; (3) 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water); (4) 2 theoretical cheese yield values (fresh cheese and cheese solids) and the related cheesemaking efficiencies; and (5) daily milk yield and 3 daily cheese yield traits (fresh cheese, cheese solids, and water retained in the curd). With respect to individual animal factors, farm was not particularly important for recovery traits or actual and theoretical cheese yield and estimates of efficiency, whereas it highly influenced daily productions. Parity of goats influenced daily cheese production, whereas DIM slightly affected recovery as well as percent and daily cheese yield traits. Breed was the most important source of variation for almost all cheesemaking traits. Compared with those of Alpine type, the 4 Mediterranean breeds had, on average, lower daily milk and cheese productions, greater actual and theoretical cheese yield, and higher recovery of nutrients in the curd. Among Alpine type, Camosciata delle Alpi was characterized by greater nutrients recovery than Saanen. Within the 4 Mediterranean types, the 3 Italians produced much less milk per day, with much more fat and protein and greater recovery traits than the Murciano-Granadina, resulting in greater actual cheese yield. Within the Italian breeds, milk from Sarda and Sarda Primitiva was characterized by lower daily yields, higher protein and fat content, and greater recoveries of nutrients than Maltese goats. These results confirmed the potential of goat milk for cheese production and could be useful to give new possibilities and direction in breeding programs.     

Genetic variability of milk components based on mid-infrared spectral data

The aim of this study was to estimate the genetic parameters of the mid-infrared (MIR) milk spectrum represented by 1,060 data points per sample. The dimensionality of traits was reduced by principal components analysis. Therefore, 46 principal components describing 99.03% of the phenotypic variability were used to create 46 new traits. Variance components were estimated using canonical transformation. Heritability ranged from 0 to 0.35. Twenty-five out of 46 studied traits showed a permanent environment variance greater than genetic variance. Eight traits showed heritability greater than 0.10. Variances of original spectral traits were obtained by back transformation. Heritabilities for each spectral data points ranged from 0.003 to 0.42. In particular, 3 MIR regions showing moderate to high heritability estimates were of potential genetic interest. Heritabilities for specific wave numbers, linked with common milk traits (e.g., lipids, lactose), were similar to those estimated for these traits. This research confirms the genetic variability of the MIR milk spectrum and, therefore, the genetic variation of milk components. The objective of this study was to better understand the genetics of milk composition and, maybe in the future, to select animals to improve milk quality.     

Estimation of genetic parameters for quantitative trait loci for dairy traits in the French Holstein population

A marker-assisted selection program (MAS) has been implemented in dairy cattle in France. The efficiency of such a selection program depends on the use of correct genetic parameters for the marked quantitative trait loci (QTL). Therefore, the objective of this study was to estimate the proportion of genetic variance explained by 4 QTL described in previous studies (these QTL are segregating on chromosomes 6, 14, 20, and 26). Genotypes for 11 markers were available for 3,974 bulls grouped within 54 sire families of the French Holstein population undergoing MAS. The parameters were estimated for 4 QTL and 5 dairy traits: milk, fat and protein yields, and fat and protein percentages. The proportion of genetic variance explained by the QTL ranged from as low as 0.03 to 0.36%. Both lack of marker informativity and poor monitoring of QTL transmission might limit the accuracy of estimation. The QTL explained a larger proportion of genetic variance for milk composition traits. The QTL on chromosome 14 and chromosomes 6 and 20 have their largest influence on fat and protein percentages, respectively. The overall proportions of genetic variance explained by the QTL were 27.0, 30.7, 24.1, 48.2, and 33.6% for milk, fat and protein yields, and fat and protein percentages, respectively. These results clearly indicated that a large part of the genetic variance is explained by a small number of QTL and that their use in MAS might be beneficial for dairy cattle breeding programs.     

Correlations of feed intake predicted with milk infrared spectra and breeding values in the Dutch Holstein population

Feed is a major cost in dairy production, and substantial genetic variation in feed efficiency exists between cows. Therefore, breeders aim to improve feed efficiency of dairy cattle. However, phenotypic data on individual feed intake on commercial farms is scarce, and accurate measurements are very costly. Several studies have shown that information from Fourier-transformed infrared spectra of milk samples (milk infrared, milk IR) can be used to predict phenotypes such as energy balance and energy intake, but this is usually based on small data sets obtained under experimental circumstances. The added value of information from milk IR spectra for estimation of breeding values is unknown. The objectives of this study were (1) to develop prediction equations for dry matter intake (DMI) and residual DMI (rDMI) from milk IR spectra; (2) to apply these for a data set of milk IR spectra from commercial Dutch dairy farms; (3) to estimate genetic parameters for these traits; and (4) to estimate correlations between these predictions and other traits in the breeding goal. We used data from feeding trials where individual feed intake was recorded daily and for which milk IR spectra were determined weekly to develop prediction equations for DMI and rDMI with partial least squares regression. This data set contained over 7,600 weekly averaged DMI records linked with milk IR spectra from 271 cows. The equations were applied for a data set with test day information from 676 Dutch dairy herds with 621,567 records of 78,488 cows. Both milk IR-predicted DMI and rDMI were analyzed with an animal model to obtain genetic parameters and sire effect estimates that could be correlated with breeding values. A partial least squares regression model with 10 components from the milk IR spectra explained around 25% of DMI variation and less than 10% of rDMI variation in the validation set. Nearly all variation in the milk IR spectra was captured by 7 components; additional components contributed marginally to the spectral variation but decreased prediction errors for both traits. Accuracies of predictions of DMI and rDMI from milk IR spectra for a large feeding experiment were 0.47 and 0.26 on average, respectively, with small differences between ration treatments (ranging from 0.43 to 0.55 and from 0.21 to 0.34, respectively) and among lactation stages (ranging from 0.24 to 0.59 and from 0.13 to 0.36, respectively), with the highest prediction accuracies in early lactation. The estimated heritabilities for predicted DMI and rDMI were 0.3 and 0.4, respectively, which suggests genetic potential for both predicted traits. The correlations of sire estimates for milk IR-predicted DMI with official Dutch breeding values were strongest with milk production (0.33), longevity (0.26), and fertility (?0.27), indicating that cows that eat more produce more, live longer, and have poorer fertility. The correlations of sire estimates for predicted DMI and rDMI with the official breeding values for DMI were low (0.14 and 0.03, respectively). This implies that the added value of including milk IR-predicted DMI information in the estimation procedure of breeding values for DMI would be considered insufficient for practical application.     

Short communication: Implementation of a breeding value for heat tolerance in Australian dairy cattle

Excessive ambient temperature and humidity can impair milk production and fertility of dairy cows. Selection for heat-tolerant animals is one possible option to mitigate the effects of heat stress. To enable selection for this trait, we describe the development of a heat tolerance breeding value for Australian dairy cattle. We estimated the direct genomic values of decline in milk, fat, and protein yield per unit increase of temperature-humidity index (THI) using 46,726 single nucleotide polymorphisms and a reference population of 2,236 sires and 11,853 cows for Holsteins and 506 sires and 4,268 cows for Jerseys. This new direct genomic value is the Australian genomic breeding value for heat tolerance (HT ABVg). The components of the HT ABVg are the decline in milk, fat, and protein per unit increase in THI when THI increases above the threshold of 60. These components are weighted by their respective economic values, assumed to be equivalent to the weights applied to milk, fat, and protein yield in the Australian selection indices. Within each breed, the HT ABVg is then standardized to have a mean of 100 and standard deviation (SD) of 5, which is consistent with the presentation of breeding values for many other traits in Australia. The HT ABVg ranged from ?4 to +3 SD in Holsteins and ?3 to +4 SD in Jerseys. The mean reliabilities of HT ABVg among validation sires, calculated from the prediction error variance and additive genetic variance, were 38% in both breeds. The range in ABVg and their reliability suggests that HT can be improved using genomic selection. There has been a deterioration in the genetic trend of HT, and to moderate the decline it is suggested that the HT ABVg should be included in a multitrait economic index with other traits that contribute to farm profit.     

Genetic parameter estimates for body conformation traits using composite index,principal component,and factor analysis

Information about genetic parameters is population specific and it is crucial for designing animal breeding programs and predicting response to selection. This study was carried out to estimate the genetic parameters for 23 body conformation traits of 45,517 Chinese Holstein reared in Eastern China from 1995 to 2017 with the Bayesian inference method using a linear animal mixed model. The methods to integrate these traits included (1) using the composite index from the Dairy Association of China and (2) applying principal component analysis and factor analysis to explore the relationship between the conformation traits. Estimates of heritability using the composite index were low (0.04; feet and legs) to moderate (0.23; body capacity). Strong genetic correlations were observed between the individual body conformation traits. Both principal components (1 to 7; eigenvalues ≥ 1) and latent factors (1 to 7; eigenvalues ≥ 1) explained 60.37% of total variability. Principal component 1 and factor 1 accounted for the traits that are usually associated with milk production. Moderate to low heritability were estimated through multi-trait analysis for principal components (from 0.07 to 0.21) and latent factors (from 0.07 to 0.23). Genetic correlations among the 2 multivariate techniques are typically lower compared with the one existing among the measured traits. Results from these analyses suggest the possibility of using both principal component analysis and factor analysis in morphological evaluation, simplifying the information given by the body conformation traits into new variables that could be useful for the genetic improvement of the Chinese Holstein population. This information could also be used to avoid analyzing large number of correlated traits, thereby improving precision and reducing computation burdens to analyze large and complex data.     

Genetic evaluation and best prediction of lactation persistency

Cows with high persistency tend to produce less milk than expected at the beginning of lactation and more than expected at the end. Best prediction of persistency was calculated as a function of a trait-specific standard lactation curve and a linear regression of test-day deviations on days in milk. Regression coefficients were deviations from a balance point to make yield and persistency phenotypically uncorrelated. The objectives of this study were to calculate (co)variance components and breeding values for best predictions of persistency of milk (PM), fat (PF), protein (PP), and SCS (PSCS) in Holstein cows. Data included 8,682,138 lactations from 4,375,938 cows calving since 1997, and 39,354 sires were evaluated. Sire estimated breeding values (EBV) for PM, PF, and PP were similar and ranged from −0.70 to 0.75 for PM; EBV for PSCS ranged from −0.37 to 0.28. Regressions of sire EBV on birth year were near zero (<0.003) but positive for PM, PF, and PP, and negative for PSCS. Genetic correlations of PM, PF, and PP with PSCS were moderate and favorable, indicating that increasing SCS decreases yield traits, as expected. Genetic correlations among yield and persistency were low to moderate and ranged from −0.09 (PSCS) to 0.18 (PF). This definition of persistency may be more useful than those used in test-day models, which are often correlated with yield. Routine genetic evaluations for persistency are feasible and may allow for improved predictions of yield traits. As calving intervals increase, persistency may have greater value.     

Comparison between genetic parameters of cheese yield and nutrient recovery or whey loss traits measured from individual model cheese-making methods or predicted from unprocessed bovine milk samples using Fourier-transform infrared spectroscopy

Cheese yield is an important technological trait in the dairy industry. The aim of this study was to infer the genetic parameters of some cheese yield-related traits predicted using Fourier-transform infrared (FTIR) spectral analysis and compare the results with those obtained using an individual model cheese-producing procedure. A total of 1,264 model cheeses were produced using 1,500-mL milk samples collected from individual Brown Swiss cows, and individual measurements were taken for 10 traits: 3 cheese yield traits (fresh curd, curd total solids, and curd water as a percent of the weight of the processed milk), 4 milk nutrient recovery traits (fat, protein, total solids, and energy of the curd as a percent of the same nutrient in the processed milk), and 3 daily cheese production traits per cow (fresh curd, total solids, and water weight of the curd). Each unprocessed milk sample was analyzed using a MilkoScan FT6000 (Foss, Hillerød, Denmark) over the spectral range, from 5,000 to 900 wavenumber × cm⁻¹. The FTIR spectrum-based prediction models for the previously mentioned traits were developed using modified partial least-square regression. Cross-validation of the whole data set yielded coefficients of determination between the predicted and measured values in cross-validation of 0.65 to 0.95 for all traits, except for the recovery of fat (0.41). A 3-fold external validation was also used, in which the available data were partitioned into 2 subsets: a training set (one-third of the herds) and a testing set (two-thirds). The training set was used to develop calibration equations, whereas the testing subsets were used for external validation of the calibration equations and to estimate the heritabilities and genetic correlations of the measured and FTIR-predicted phenotypes. The coefficients of determination between the predicted and measured values in cross-validation results obtained from the training sets were very similar to those obtained from the whole data set, but the coefficient of determination of validation values for the external validation sets were much lower for all traits (0.30 to 0.73), and particularly for fat recovery (0.05 to 0.18), for the training sets compared with the full data set. For each testing subset, the (co)variance components for the measured and FTIR-predicted phenotypes were estimated using bivariate Bayesian analyses and linear models. The intraherd heritabilities for the predicted traits obtained from our internal cross-validation using the whole data set ranged from 0.085 for daily yield of curd solids to 0.576 for protein recovery, and were similar to those obtained from the measured traits (0.079 to 0.586, respectively). The heritabilities estimated from the testing data set used for external validation were more variable but similar (on average) to the corresponding values obtained from the whole data set. Moreover, the genetic correlations between the predicted and measured traits were high in general (0.791 to 0.996), and they were always higher than the corresponding phenotypic correlations (0.383 to 0.995), especially for the external validation subset. In conclusion, we herein report that application of the cross-validation technique to the whole data set tended to overestimate the predictive ability of FTIR spectra, give more precise phenotypic predictions than the calibrations obtained using smaller data sets, and yield genetic correlations similar to those obtained from the measured traits. Collectively, our findings indicate that FTIR predictions have the potential to be used as indicator traits for the rapid and inexpensive selection of dairy populations for improvement of cheese yield, milk nutrient recovery in curd, and daily cheese production per cow.     

Use of test-day records beyond three hundred five days for estimation of three hundred five-day breeding values for production traits and somatic cell score of Canadian Holsteins

The Canadian Test-Day Model includes test-day (TD) records from 5 to 305 d in milk (DIM). Because 60% of Canadian Holstein cows have at least one lactation longer than 305 d, a significant number of TD records beyond 305 DIM could be included in the genetic evaluation. The aim of this study was to investigate whether TD records beyond 305 DIM could be useful for estimation of 305-d estimated breeding value (EBV) for milk, fat, and protein yields and somatic cell score. Data were 48,638,184 TD milk, fat, and protein yields and somatic cell scores from the first 3 lactations of 2,826,456 Canadian Holstein cows. All production traits were preadjusted for the effect of pregnancy. Subsets of data were created for variance-component estimation by random sampling of 50 herds. Variance components were estimated using Gibbs sampling. Full data sets were used for estimation of breeding values. Three multiple-trait, multiple-lactation random regression models with TD records up to 305 DIM (M305), 335 DIM (M335), and 365 DIM (M365) were fitted. Two additional models (M305a and M305b) used TD records up to 305 DIM and variance components previously estimated by M335 and M365, respectively. The effects common to all models were fixed effects of herd × test-date and DIM class, fixed regression on DIM nested within region × age × season class, and random regressions for additive genetic and permanent environmental effects. Legendre polynomials of order 6 and 4 were fitted for fixed and random regressions, respectively. Rapid increase of additive genetic and permanent environmental variances at extremes of lactations was observed with all 3 models. The increase of additive genetic and permanent environmental variances was at earlier DIM with M305, resulting in greater variances at 305 DIM with M305 than with M335 and M365. Model M305 had the best ability to predict TD yields from 5 through 305 DIM and less error of prediction of 305-d EBV than M335 and M365. Model M335 had smaller change of 305-d EBV of bulls over the period of 7 yr than did M305 and M365. Model M305a had the least error of prediction and change of 305-d EBV from all models. Therefore, the use of TD records of Holstein cows from 5 through 305 DIM and variance components estimated using records up to 335 DIM is recommended for the Canadian Test-Day Model.     

Correlation between breeding values for milk fatty acids and Nordic Total Merit index traits for Danish Holstein and Danish Jersey

Milk fatty acid composition is gaining interest in the Danish dairy industry both to develop new dairy products and as a management tool. To be able to implement milk fatty acid (FA) composition in the breeding program, it is important to know the correlations with the traits in the breeding goal. To estimate these correlations, we measured milk fat composition in Danish Holstein (DH) and Danish Jersey (DJ) cattle breeds using mid-infrared spectroscopy. Breeding values were estimated for specific FA and for groups of FA. Correlations with the estimated breeding values (EBV) underlying the Nordic Total Merit index (NTM) were calculated within breed. For both DH and DJ, we showed that FA EBV had moderate correlations with the NTM and production traits. For both DH and DJ, the correlation of FA EBV and NTM were in the same direction, except for C16:0 (0 in DH, 0.23 in DJ). A few correlations differed between DH and DJ. The correlation between claw health index and C18:0 was negative in DH (−0.09) but positive in DJ (0.12). In addition, some correlations were not significant in DH but were significant in DJ. The correlations between udder health index and long-chain FA, trans FA, C16:0, and C18:0 were not significant in DH (−0.05 to 0.02), but were significant in DJ (−0.17, −0.15, 0.14, and −0.16, respectively). For both DH and DJ, the correlations between FA EBV and nonproduction traits were low. This implies that it is possible to breed for a different fat composition in the milk without affecting the nonproduction traits in the breeding goal.     

Effects of breed,farm intensiveness,and cow productivity level on cheese-making ability predicted using infrared spectral data at the population level

Fourier-transform infrared (FTIR) spectra collected during milk recording schemes at population level can be used for predicting novel traits of interest for farm management, cows' genetic improvement, and milk payment systems. The aims of this study were as follows. (1) To predict cheese yield traits using FTIR spectra from routine milk recordings exploiting previously developed calibration equations. (2) To compare the predicted cheese-making abilities of different dairy and dual-purpose breeds. (3) To analyze the effects of herds' level of intensiveness (HL) and of the cow's level of productivity (CL). (4) To compare the patterns of predicted cheese yields with the patterns of milk composition in different breeds to discern the drivers of cheese-making efficiency. The major sources of variation of FTIR predictions of cheese yield ability (fresh cheese or cheese solids produced per unit milk) of individual milk samples were studied on 115,819 cows of 4 breeds (2 specialized dairy breeds, Holstein and Brown Swiss, and 2 dual-purpose breeds, Simmental and Alpine Grey) from 6,430 herds and exploiting 1,759,706 FTIR test-day spectra collected over 7 yr of milk sampling. Calibration equations used were previously developed on 1,264 individual laboratory model cheese procedures (cross-validation R² 0.85 and 0.95 for fresh and solids cheese yields, respectively). The linear model used for statistical analysis included the effects of parity, lactation stage, year of calving, month of sampling, HL, CL, breed of cow, and the interactions breed × HL and breed × CL. The HL and CL stratifications (5 classes each) were based on average daily secretion of milk net energy per cow. All effects were highly significant. The major conclusions were as follows. (1) The FTIR-based prediction of cheese yield of milk goes beyond the knowledge of fat and protein content, partially explaining differences in cheese-making ability in different cows, breeds and herds. (2) Differences in cheese yields of different breeds are only partially explained by milk fat and protein composition, and less productive breeds are characterized by a higher milk nutrient content as well as a higher recovery of nutrients in the cheese. (3) High-intensive herds not only produce much more milk, but the milk has a higher nutrient content and a higher cheese yield, whereas within herds, compared with less productive cows, the more productive cows have a much greater milk yield, milk with a greater content of fat but not of protein, and a moderate improvement in cheese yield, differing little from expectations based on milk composition. Finally, (4) the effects of HL and CL on milk quality and cheese-making ability are similar but not identical in different breeds, the less productive ones having some advantage in terms of cheese-making ability. We can obtain FTIR-based prediction of cheese yield from individual milk samples retrospectively at population level, which seems to go beyond the simple knowledge of milk composition, incorporating information on nutrient retention ability in cheese, with possible advantages for management of farms, genetic improvement of dairy cows, and milk payment systems.     

An empirical method for prediction of cheese yield

Theoretical cheese yield can be estimated from the milk fat and casein or protein content of milk using classical formulae, such as the VanSlyke formula. These equations are reliable predictors of theoretical or actual yield based on accurately measured milk fat and casein content. Many cheese makers desire to base payment for milk to dairy farmers on the yield of cheese. In small factories, however, accurate measurement of fat and casein content of milk by either chemical methods or infrared milk analysis is too time consuming and expensive. Therefore, an empirical test to predict cheese yield was developed which uses simple equipment (i.e., clinical centrifuge, analytical balance, and forced air oven) to carry out a miniature cheese making, followed by a gravimetric measurement of dry weight yield. A linear regression of calculated theoretical versus dry weight yields for milks of known fat and casein content was calculated. A regression equation of y = 1.275x + 1.528, where y is theoretical yield and x is measured dry solids yield (r2 = 0.981), for Cheddar cheese was developed using milks with a range of theoretical yield from 7 to 11.8%. The standard deviation of the difference (SDD) between theoretical cheese yield and dry solids yield was 0.194 and the coefficient of variation (SDD/mean x 100) was 1.95% upon cross validation. For cheeses without a well-established theoretical cheese yield equation, the measured dry weight yields could be directly correlated to the observed yields in the factory; this would more accurately reflect the expected yield performance. Payments for milk based on these measurements would more accurately reflect quality and composition of the milk and the actual average recovery of fat and casein achieved under practical cheese making conditions.     

Assessing the effect of pregnancy stage on milk composition of dairy cows using mid-infrared spectra

Changes in milk production traits (i.e., milk yield, fat, and protein contents) with the pregnancy stage are well documented. To our knowledge, the effect of pregnancy on the detailed milk composition has not been studied so far. The mid-infrared (MIR) spectrum reflects the detailed composition of a milk sample and is obtained by a nonexhaustive and widely used method for milk analysis. Therefore, this study aimed to investigate the effect of pregnancy on milk MIR spectrum in addition to milk production traits (milk yield, fat, and protein contents). A model including regression on the number of days pregnant was applied on milk production traits (milk yield, fat, and protein contents) and on 212 spectral points from the MIR spectra of 9,757 primiparous Holstein cows from Walloon herds. Effects of pregnancy stage were expressed on a relative scale (effect divided by the squared root of the phenotypic variance); this allowed comparisons between effects on milk traits and on 212 spectral points. Effect of pregnancy stage on production traits were in line with previous studies indicating that the model accounted well for the pregnancy effect. Trends of the relative effect of the pregnancy stage on the 212 spectral points were consistent with known and observed effect on milk traits. The highest effect of the pregnancy was observed in the MIR spectral region from 968 to 1,577 cm^?1. For some specific wavenumbers, the effect was higher than for fat and protein contents in the beginning of the pregnancy (from 30 to 90 or 120 d pregnant). In conclusion, the effect of early pregnancy can be observed in the detailed milk composition through the analysis of the MIR spectrum of bovine milk. Further analyses are warranted to explore deeply the use of MIR spectra of bovine milk for breeding and management of dairy cow pregnancy.     

Standardization of milk infrared spectra for the retroactive application of calibration models

The objective of this study was to standardize the infrared spectra obtained over time and across 2 milk laboratories of Canada to create a uniform historical database and allow (1) the retroactive application of calibration models for prediction of fine milk composition; and (2) the direct use of spectral information for the development of indicators of animal health and efficiency. Spectral variation across laboratories and over time was inspected by principal components analysis (PCA). Shifts in the PCA scores were detected over time, leading to the definition of different subsets of spectra having homogeneous infrared signal. To evaluate the possibility of using common equations on spectra collected by the 2 instruments and over time, we developed a standardization (STD) method. For each subset of data having homogeneous infrared signal, a total of 99 spectra corresponding to the percentiles of the distribution of the absorbance at each wavenumber were created and used to build the STD matrices. Equations predicting contents of saturated fatty acids, short-chain fatty acids, and C18:0 were created and applied on different subsets of spectra, before and after STD. After STD, bias and root mean squared error of prediction decreased by 66% and 32%, respectively. When calibration equations were applied to the historical nonstandardized database of spectra, shifts in the predictions could be observed over time for all investigated traits. Shifts in the distribution of the predictions over time corresponded to the shifts identified by the inspection of the PCA scores. After STD, shifts in the predicted fatty acid contents were greatly reduced. Standardization reduced spectral variability between instruments and over time, allowing the merging of milk spectra data from different instruments into a common database, the retroactive use of calibrations equations, or the direct use of the spectral data without restrictions.     

Milk metagenomics and cheese-making properties as affected by indoor farming and summer highland grazing

The study of the complex relationships between milk metagenomics and milk composition and cheese-making efficiency as affected by indoor farming and summer highland grazing was the aim of the present work. The experimental design considered monthly sampling (over 5 mo) of the milk produced by 12 Brown Swiss cows divided into 2 groups: the first remained on a lowland indoor farm from June to October, and the second was moved to highland pastures in July and then returned to the lowland farm in September. The resulting 60 milk samples (2 kg each) were used to analyze milk composition, milk coagulation, curd firming, and syneresis processes, and to make individual model cheeses to measure cheese yields and nutrient recoveries in the cheese. After DNA extraction and Illumina Miseq sequencing, milk microbiota amplicons were also processed by means of an open-source pipeline called Quantitative Insights Into Microbial Ecology (Qiime2, version 2018.2; https://qiime2.org). Out of a total of 44 taxa analyzed, 13 bacterial taxa were considered important for the dairy industry (lactic acid bacteria, LAB, 5 taxa; and spoilage bacteria, 4) and for human (other probiotics, 2) and animal health (pathogenic bacteria, 2). The results revealed the transhumant group of cows transferred to summer highland pastures showed an increase in almost all the LAB taxa, bifidobacteria, and propionibacteria, and a reduction in spoilage taxa. All the metagenomic changes disappeared when the transhumant cows were moved back to the permanent indoor farm. The relationships between 17 microbial traits and 30 compositional and technological milk traits were investigated through analysis of correlation and latent explanatory factor analysis. Eight latent factors were identified, explaining 75.3% of the total variance, 2 of which were mainly based on microbial traits: pro-dairy bacteria (14% of total variance, improving during summer pasturing) and pathogenic bacteria (6.0% of total variance). Some bacterial traits contributed to other compositional-technological latent factors (gelation, udder health, and caseins).