A. Fleming

A 100-Year Review: Identification and genetic selection of economically important traits in dairy cattle

Over the past 100 yr, the range of traits considered for genetic selection in dairy cattle populations has progressed to meet the demands of both industry and society. At the turn of the 20th century, dairy farmers were interested in increasing milk production; however, a systematic strategy for selection was not available. Organized milk performance recording took shape, followed quickly by conformation scoring. Methodological advances in both genetic theory and statistics around the middle of the century, together with technological innovations in computing, paved the way for powerful multitrait analyses. As more sophisticated analytical techniques for traits were developed and incorporated into selection programs, production began to increase rapidly, and the wheels of genetic progress began to turn. By the end of the century, the focus of selection had moved away from being purely production oriented toward a more balanced breeding goal. This shift occurred partly due to increasing health and fertility issues and partly due to societal pressure and welfare concerns. Traits encompassing longevity, fertility, calving, health, and workability have now been integrated into selection indices. Current research focuses on fitness, health, welfare, milk quality, and environmental sustainability, underlying the concentrated emphasis on a more comprehensive breeding goal. In the future, on-farm sensors, data loggers, precision measurement techniques, and other technological aids will provide even more data for use in selection, and the difficulty will lie not in measuring phenotypes but rather in choosing which traits to select for.

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.

Genetic analysis of groups of mid-infrared predicted fatty acids in milk

The objective of this study was to investigate genetic variability of mid-infrared predicted fatty acid groups in Canadian Holstein cattle. Genetic parameters were estimated for 5 groups of fatty acids: short-chain (4 to 10 carbons), medium-chain (11 to 16 carbons), long-chain (17 to 22 carbons), saturated, and unsaturated fatty acids. The data set included 49,127 test-day records from 10,029 first-lactation Holstein cows in 810 herds. The random regression animal test-day model included days in milk, herd-test date, and age-season of calving (polynomial regression) as fixed effects, herd-year of calving, animal additive genetic effect, and permanent environment effects as random polynomial regressions, and random residual effect. Legendre polynomials of the third degree were selected for the fixed regression for age-season of calving effect and Legendre polynomials of the fourth degree were selected for the random regression for animal additive genetic, permanent environment, and herd-year effect. The average daily heritability over the lactation for the medium-chain fatty acid group (0.32) was higher than for the short-chain (0.24) and long-chain (0.23) fatty acid groups. The average daily heritability for the saturated fatty acid group (0.33) was greater than for the unsaturated fatty acid group (0.21). Estimated average daily genetic correlations were positive among all fatty acid groups and ranged from moderate to high (0.63-0.96). The genetic correlations illustrated similarities and differences in their origin and the makeup of the groupings based on chain length and saturation. These results provide evidence for the existence of genetic variation in mid-infrared predicted fatty acid groups, and the possibility of improving milk fatty acid profile through genetic selection in Canadian dairy cattle.

Prediction of milk fatty acid content with mid-infrared spectroscopy in Canadian dairy cattle using differently distributed model development sets

The fatty acid profile of milk is a prevailing issue due to the potential negative or positive effects of different fatty acids to human health and nutrition. Mid-infrared spectroscopy can be used to obtain predictions of otherwise costly fatty acid phenotypes in a widespread and rapid manner. The objective of this study was to evaluate the prediction of fatty acid content for the Canadian dairy cattle population from mid-infrared spectral data and to compare the results produced by altering the partial least squares (PLS) model development set used. The PLS model development sets used to develop the predictions were reference fatty acids expressed as (1) grams per 100 g of fatty acid, (2) grams per 100 g of milk, (3) the natural logarithmic transform of grams per 100 g of milk, and (4) subsets of samples randomly selected by removing excess records around the mean to present a more uniform distribution, repeated 10 times. Gas chromatography measured fatty acid concentration and spectral data for 2,023 milk samples of 373 cows from 4 breeds and 44 herds were used in the model development. The coefficient of determination of cross-validation (Rcv2) increased when fatty acids were expressed on a per 100 g of milk basis compared with on a per 100 g of fat basis for all examined fatty acids. The logarithmic transformation used to create a more Gaussian distribution in the development set had little effect on the prediction accuracy. The individual fatty acids C12:0, C14:0, C16:0, C18:0, C18:1n-9 cis, and saturated, monounsaturated, unsaturated, short-chain, medium-chain, and long-chain fatty acid groups had (Rcv2) greater than 0.70. When model development was performed with subsets of the original samples, slight increases in (Rcv2) values were observed for the majority of fatty acids. The difference in (Rcv2) between the top- and bottom-performing prediction equation across the different subsets for a single predicted fatty acid was on average 0.055 depending on which samples were randomly selected to be used in the PLS model development set. Predictions for fatty acids with high accuracies can be used to monitor fatty acid contents for cows in milk recording programs and possibly for genetic evaluation.

Variation in fat globule size in bovine milk and its prediction using mid-infrared spectroscopy

The objectives of this study were to investigate the sources of variation in milk fat globule (MFG) size in bovine milk and its prediction using mid-infrared (MIR) spectroscopy. Mean MFG size was measured in 2,076 milk samples from 399 Ayrshire, Brown Swiss, Holstein, and Jersey cows, and expressed as volume moment mean (D[4,3]) and surface moment mean (D[3,2]). The mid-infrared spectra of the samples and milk performance data were also recorded during routine milk recording and testing. The effects of breed, herd nested within breed, days in milk, season, milking period, age at calving, parity, and individual animal on the variation observed in MFG size were investigated. Breed, herd nested within breed, days in milk, season, and milking period significantly affected mean MFG size. Milk fat globule size was the largest at the beginning of lactation and subsequently decreased. Milk samples with the smallest MFG on average came from Holstein cows, and those with the largest were from Jersey and Brown Swiss cows. Partial least squares regression was used to predict MFG size from MIR spectra of samples with a calibration data set containing 2,034 and 2,032 samples for D[4,3] and D[3,2], respectively. Coefficients of determination of cross validation for D[4,3] and D[3,2] prediction models were 0.51 and 0.54, respectively. The associated ratio of performance deviation values were 1.43 and 1.48 for D[4,3] and D[3,2], respectively. With these models, individual mean MFG size could not be accurately predicted, but results may be sufficient to screen samples for having either small or large MFG on average. Significant but low correlations of D[4,3] and D[3,2] with milk fat yield were estimated (0.16 and 0.21, respectively). Significant and moderate Pearson correlation coefficients for fat percent with D[4,3] and D[3,2] were assessed (0.34 and 0.36, respectively). This correlation was greater between milk fat percentage and predicted MFG size than with measured MFG size with coefficients of 0.47 and 0.49 for D[4,3] and D[3,2], respectively. The MIR prediction equations are potentially overusing the correlation between fat and MFG size and exploiting the strong relationship between the MIR spectra and total milk fat. However, the predictions of MFG size are able to determine variation in mean globule size beyond what would be achieved just by looking at the correlation with fat production.

Genetic correlations of mid-infrared-predicted milk fatty acid groups with milk production traits

The objective of this research was to estimate the genetic correlations between milk mid-infrared-predicted fatty acid groups and production traits in first-parity Canadian Holsteins. Contents of short-chain, medium-chain, long-chain, saturated, and unsaturated fatty acid groupings in milk samples can be predicted using mid-infrared spectral data for cows enrolled in milk recording programs. Predicted fatty acid group contents were obtained for 49,127 test-day milk samples from 10,029 first-parity Holstein cows in 810 herds. Milk yield, fat and protein yield, fat and protein percentage, fat-to-protein ratio, and somatic cell score were also available for these test days. Genetic parameters were estimated for the fatty acid groups and production traits using multiple-trait random regression test day models by Bayesian methods via Gibbs sampling. Three separate 8- or 9-trait analyses were performed, including the 5 fatty acid groups with different combinations of the production traits. Posterior standard deviations ranged from <0.001 to 0.01. Average daily genetic correlations were negative and similar to each other for the fatty acid groups with milk yield (-0.62 to -0.59) and with protein yield (-0.32 to -0.25). Weak and positive average daily genetic correlations were found between somatic cell score and the fatty acid groups (from 0.25 to 0.36). Stronger genetic correlations with fat yield, fat and protein percentage, and fat-to-protein ratio were found with medium-chain and saturated fatty acid groups compared with those with long-chain and unsaturated fatty acid groups. Genetic correlations were very strong between the fatty acid groups and fat percentage, ranging between 0.88 for unsaturated and 0.99 for saturated fatty acids. Daily genetic correlations from 5 to 305 d in milk with milk, protein yield and percentage, and somatic cell score traits showed similar patterns for all fatty acid groups. The daily genetic correlations with fat yield at the beginning of lactation were decreasing for long-chain and unsaturated fatty acid groups and increasing for short-chain fatty acids. Genetic correlations between fat percentage and fatty acids were increasing at the beginning of lactation for short- and medium-chain and saturated fatty acids, but slightly decreasing for long-chain and unsaturated fatty acid groups. These results can be used in defining fatty acid traits and breeding objectives.

Phenotypic investigation of fine milk components in bovine milk and their prediction using mid-infrared spectroscopy

This study aimed to examine the phenotypic variation observed in fine milk components, the use of mid-infrared (MIR) spectroscopy to predict these components, and the correlations with other milk p...

Heritabilities of measured and mid-infrared predicted milk fat globule size, milk fat and protein percentages, and their genetic correlations

The objective of this study was to estimate the heritability of milk fat globule (MFG) size and mid-infrared (MIR) predicted MFG size in Holstein cattle. The genetic correlations between measured and predicted MFG size with milk fat and protein percentage were also investigated. Average MFG size was measured in 1,583 milk samples taken from 254 Holstein cows from 29 herds across Canada. Size was expressed as volume moment mean (D[4,3]) and surface moment mean (D[3,2]). Analyzed milk samples also had average MFG size predicted from their MIR spectral records. Fat and protein percentages were obtained for all test-day milk samples in the cows lactation. Univariate and bivariate repeatability animal models were used to estimate heritability and genetic correlations. Moderate heritabilities of 0.364 and 0.466 were found for D[4,3] and D[3,2], respectively, and a strong genetic correlation was found between the 2 traits (0.98). The heritabilities for the MIR-predicted MFG size were lower than those estimated for the measured MFG size at 0.300 for predicted D[4,3] and 0.239 for predicted D[3,2]. The genetic correlation between measured and predicted D[4,3] was 0.685; the correlation was slightly higher between measured and predicted D[3,2] at 0.764, likely due to the better prediction accuracy of D[3,2]. Milk fat percentage had moderate genetic correlations with both D[4,3] and D[3,2] (0.538 and 0.681, respectively). The genetic correlation between predicted MFG size and fat percentage was much stronger (greater than 0.97 for both predicted D[4,3] and D[3,2]). The stronger correlation suggests a limitation for the use of the predicted values of MFG size as indicator traits for true average MFG size in milk in selection programs. Larger samples sizes are required to provide better evidence of the estimated genetic parameters. A genetic component appears to exist for the average MFG size in bovine milk, and the variation could be exploited in selection programs.