Using milk mid-infrared spectroscopy to estimate cow-level nitrogen efficiency metrics.

Minimizing pollution from the dairy sector is paramount; one potential cause of such pollution is excess nitrogen. Nitrogen pollution contributes to a deterioration in water quality as well as an increase in both eutrophication and greenhouse gases. It is therefore essential to minimize the loss of nitrogen from the sector, including excretion from the cow. Breeding programs are one potential strategy to improve the efficiency with which nitrogen is used by dairy cows but relies on routine access to individual cow information on how efficiently each cows uses the nitrogen it ingests. A total of 3,497 test-day records for individual cow nitrogen efficiency metrics along with milk yield and the associated milk spectra were used to investigate the ability of milk infrared spectral data to predict these nitrogen traits; both traditional partial least squares regression and neural networks were used in the prediction process. The data originated from 4 farms across 11 years. The nitrogen traits investigated were nitrogen intake, nitrogen use efficiency, and nitrogen balance. Both nitrogen use efficiency and nitrogen balance were calculated considering nitrogen intake, nitrogen in milk, nitrogen in the conceptus, nitrogen used for the growth, nitrogen stored in the reserves, and nitrogen mobilized from the reserves. Irrespective of the nitrogen-related trait being investigated, the best prediction from 4-fold cross-validation were achieved using neural networks that considered both the morning and evening milk spectra along with milk yield, parity, and days in milk in the prediction process. The coefficient of determination in the cross-validation was 0.61, 0.74, and 0.58 for nitrogen intake, nitrogen use efficiency, and nitrogen balance, respectively. In a separate series of validation approaches, the calibration and validation was stratified by herd (n = 4) and separately by year. For these scenarios, partial least squares regression generated more accurate predictions compared with neural networks; the coefficient of determination was always lower than 0.29 and 0.60 when validation was stratified by herd and year, respectively. Therefore, if the variability of the data being predicted in the validation data sets is similar to that in the data used to develop the predictions, then nitrogen-related traits can be predicted with reasonable accuracy. In contrast, where the variability of the data that exists in the validation data set is poorly represented in the calibration data set, then poor predictions will ensue.

Genetic covariance components for measures of nitrogen utilization in grazing dairy cows.

Improved nitrogen utilization of dairy production systems should improve not only the economic output of the systems but also the environmental metrics. One strategy to improve efficiency is through breeding programs. Improving a trait through breeding is conditional on the presence of exploitable genetic variability. Using a database of 1,291 deeply phenotyped grazing dairy cows, the genetic variability for 2 definitions of nitrogen utilization was studied: nitrogen use efficiency (i.e., nitrogen output in milk and meat divided by nitrogen available) and nitrogen balance (i.e., nitrogen available less nitrogen output in milk and meat). Variance components for both variables were estimated using animal repeatability linear mixed models. Genetic variability was detected for both nitrogen utilization metrics, even though their heritability estimates were low (<0.10). Validation of genetic evaluations revealed that animals divergent for nitrogen use efficiency or nitrogen balance indeed differed phenotypically, further demonstrating that breeding for improved nitrogen efficiency should result in a shift in the population mean toward better efficiency. Nitrogen use efficiency and nitrogen balance were not genetically correlated with each other (<|0.28|), and neither metric was correlated with milk urea nitrogen (<|0.12|). Nitrogen balance was unfavorably correlated with milk yield, showing the importance of including the nitrogen utilization metrics in a breeding index to improve nitrogen utilization without negatively impacting milk yield. In conclusion, improvement of nitrogen utilization through breeding is possible, even if more nitrogen utilization phenotypic data need to be collected to improve the selection accuracy considering the low heritability estimates.

The association between calf birth weight and the postcalving performance of its dairy dam in the absence of dystocia.

The objective of the study was to quantify the association between the birth weight of a calf and the subsequent performance of its dairy dam in the absence of any recorded calving assistance. A total of 11,592 lactation records from 4,549 spring-calving dairy cows were used. The association between a series of quantitative cow performance metrics (dependent variable) and calf birth weight (independent variable) was determined using linear mixed models; logistic regression was used where the dependent variable was binary. Nuisance factors in the models were calf sex, heterosis coefficient of both the cow and calf, dry period length immediately before the birth of the calf, cow age at calving relative to the median cow age per parity, breed proportion of the cow, cow live weight between 100 and 200 d of lactation relative to the mean cow weight per parity, and contemporary group. Calf birth weight was included in the model as either a continuous or a categorical variable. Primiparous and multiparous cows were analyzed separately. Mean (SD) calf birth weight was 36.2 (6.8) kg. In primiparous cows, calf birth weight was associated with milk yield in the first 60 d of lactation, calving to first service interval, calving body weight (BW), and both nadir BW and body condition score (BCS). In multiparous cows, calf birth weight was associated with total milk, fat, and protein yield in the first 60 and 305 d of lactation, peak milk yield, total milk solids, both calving and nadir BW, and BCS loss from calving to nadir. Relative to primiparous cows that gave birth to calves weighing 34 to 37 kg (i.e., population mean), their contemporaries who gave birth to calves that weighed 15 to 29 kg produced 9.82 kg more milk in the first 60 d of lactation, had a 2-d shorter interval to first service, and were 8.08 kg and 5.51 kg lighter at calving and nadir BW, respectively; the former was also 0.05 units lower in BCS (5-point scale, 1 = emaciated and 5 = obese) at nadir. Relative to multiparous cows that gave birth to calves that were 34 to 37 kg birth weight, multiparous cows that gave birth to calves that were 15 to 29 kg yielded 59.63 kg, 2.44 kg, and 1.76 kg less milk, fat, and protein, respectively, in the first 60 d of lactation; produced 17.69 kg less milk solids throughout the 305-d lactation; and were also 10.49 kg lighter at nadir and lost 0.01 units more BCS to nadir. In a separate series of analyses, sire breed was added to the model as a fixed effect with and without calf birth weight. When calf birth weight was not adjusted for, 60-d milk yield for multiparous cows who gave birth to calves sired by a traditional beef breed (i.e., Angus, Hereford) produced 59.63 kg more milk than multiparous cows who gave birth to calves sired by a Holstein-Friesian. Hence, calf birth weight is associated with some subsequent dam performance measures; however, where associations do exist, the effect is biologically small.

Predicting methane emissions of individual grazing dairy cows from spectral analyses of their milk samples.

Data on the enteric methane emissions of individual cows are useful not just in assisting management decisions and calculating herd inventories but also as inputs for animal genetic evaluations. Data generation for many animal characteristics, including enteric methane emissions, can be expensive and time consuming, so being able to extract as much information as possible from available samples or data sources is worthy of investigation. The objective of the present study was to attempt to predict individual cow methane emissions from the information contained within milk samples, specifically the spectrum of light transmittance across different wavelengths of the mid-infrared (MIR) region of the electromagnetic spectrum. A total of 93,888 individual spot measures of methane (i.e., individual samples of an animal's breath when using the GreenFeed technology) from 384 lactations on 277 grazing dairy cows were collapsed into weekly averages expressed as grams per day; each weekly average coincided with a MIR spectral analysis of a morning or evening individual cow milk sample. Associations between the spectra and enteric methane measures were performed separately using partial least squares regression or neural networks with different tuning parameters evaluated. Several alternative definitions of the enteric methane phenotype (i.e., average enteric methane in the 6 d preceding or 6 d following taking the milk sample or the average of the 6 d before and after the milk sample, all of which also included the enteric methane emitted on the day of milk sampling), the candidate model features (e.g., milk yield, milk composition, and milk MIR) as well as validation strategy (i.e., cross-validation or leave-one-experimental treatment-out) were evaluated. Irrespective of the validation method, the prediction accuracy was best when the average of the milk MIR from the morning and evening milk sample was used and the prediction model was developed using neural networks; concurrently including milk yield and days in milk in the prediction model generated superior predictions relative to just the spectral information alone. Furthermore, prediction accuracy was best when the enteric methane phenotype was the average of at least 20 methane spot measures across a 6-d period flanking each side of the milk sample with associated spectral data. Based on the strategy that achieved the best accuracy of prediction, the correlation between the actual and predicted daily methane emissions when based on 4-fold cross-validation varied per validation stratum from 0.68 to 0.75; the corresponding range when validated on each of the 8 different experimental treatments focusing on alternative pasture grazing systems represented in the dataset varied from 0.55 to 0.71. The root mean square error of prediction across the 4-folds of cross-validation was 37.46 g/d, whereas the root mean square error averaged across all folds of leave-one-treatment-out was 37.50 g/d. Results suggest that even with the likely measurement errors contained within the MIR spectrum and gold standard enteric methane phenotype, enteric methane can be reasonably well predicted from the infrared spectrum of milk samples. What is yet to be established, however, is whether (a) genetic variation exists in this predicted enteric methane phenotype and (b) selection on estimates of genetic merit for this phenotype translate to actual phenotypic differences in enteric methane emissions.

Impact on prime animal beef merit from breeding solely for lighter dairy cows.

As the proportion of prime carcasses originating from dairy herds increases, the focus is shifting to the beef merit of the progeny from dairy herds. Several dairy cow total merit indexes include a negative weight on measures of cow size. However, there is a lack of knowledge on the impact of genetic selection, solely for lighter or smaller-sized dairy cows, on the beef performance of their progeny. Therefore, the objective of this study was to quantify the genetic correlations among cow size traits (i.e., cow body weight (BW), cow carcass weight (CW)), cow body condition score (BCS), cow carcass conformation (CC), and cow carcass fat cover (CF), as well as the correlations between these cow traits and a series of beef performance slaughter-related traits (i.e., CW, CC, CF, and age at slaughter (AS)) in their progeny. After data editing, there were 52,950 cow BW and BCS records, along with 57,509 cow carcass traits (i.e., CW, CC, and CF); carcass records from 346,350 prime animals along with AS records from 316,073 prime animals were also used. Heritability estimates ranged from moderate to high (0.18 to 0.62) for all cow and prime animal traits. The same carcass trait in cows and prime animals were strongly genetically correlated with each other (0.76 to 0.85), implying that they are influenced by very similar genomic variants. Selecting exclusively for cows with higher BCS (i.e., fatter) will, on average, produce more conformed prime animals carcasses, owing to a moderate genetic correlation (0.30) between both traits. Genetic correlations revealed that selecting exclusively for lighter BW or CW cows will, on average, result in lighter prime animal carcasses of poor CC, while also delaying slaughter age. Nonetheless, selective breeding through total merit indexes should be successful in breeding for smaller dairy cows, and desirable prime animal carcass traits concurrently, because of the non-unity genetic correlations between the cow and prime animal traits; this will help to achieve a more ethical, environmentally sustainable, and economically viable dairy-beef industry.

Mean breed performance of the progeny from beef-on-dairy matings.

Gains through breeding can be achieved through a combination of both between-breed and within-breed selection. Two suites of traits of particular interest to dairy producers when selecting beef bulls for mating to dairy females are calving-related attributes and the expected value of the subsequent calf, the latter usually being a function of expected carcass value. Estimated breed effects can be informative, particularly in the absence of across-breed genetic evaluations. The objective of the present study was to use a large national database of the progeny from beef-on-dairy matings to estimate the mean breed effects of the used beef sires. Calving performance (i.e., gestation length, calving difficulty score, and perinatal morality) as well as calf value were investigated; a series of slaughter-related traits (i.e., carcass metrics and age at slaughter) of the prime progeny were also investigated. Phenotypic data on up to 977,037 progeny for calving performance, 79,903 for calf price and 103,175 for carcass traits (including dairy × dairy progeny for comparative purposes) were used; sire breeds represented were Holstein-Friesian, Angus, Aubrac, Belgian Blue, Charolais, Hereford, Limousin, Salers, and Simmental. Large interbreed differences existed. The mean gestation length of male calves from beef sires varied from 282.3 d (Angus) to 287.4 d (Limousin) which were all longer than the mean of 280.9 d for Holstein-Friesian sired male calves. Relative to a Holstein-Friesian sire, the odds of dystocia varied from 1.43 (Angus) to 4.77 (Belgian Blue) but, once adjusted for both the estimated maternal genetic merit of the dam and direct genetic merit of the calf for calving difficulty, the range in odds ratios shrunk. A difference of €125.4 existed in calf sale price between the progeny of the different beef breeds investigated which represented over twice the residual standard deviation in calf price within the day of sale-Angus was the cheapest while Charolais calves were, on average, the most expensive calves. Mean carcass weight of steers, not adjusted for age at slaughter or carcass fat, varied from 327.1 kg (Angus) to 363.2 kg (Belgian Blue) for the beef breeds with the mean carcass weight of Holstein-Friesian steer progeny being 322.4 kg. Belgian Blues had, on average, the best carcass conformation with the Herefords and Angus having the worst of all beef breeds. Angus and Hereford steers were slaughtered the youngest of all beef breeds but just 9 d younger than the average of all other beef breeds yet 24 d younger than Holstein-Friesian sired progeny. Clear breed differences in calving and carcass performance exist among beef breeds mated to dairy females. Those breeds excelling in calving performance were not necessarily the best for carcass merit.

Usefulness of mid-infrared spectroscopy as a tool to estimate body condition score change from milk samples in intensively fed dairy cows.

Directly measuring individual cow energy balance is not trivial. Other traits such as body condition score (BCS) and BCS change (ΔBCS) can, however, be used as an indicator of cow energy status. Body condition score is a metric used worldwide to estimate cow body reserves, but the estimation of ΔBCS was, until now, conditional on the availability of multiple BCS assessments. The aim of the present study was to estimate ΔBCS from milk mid-infrared (MIR) spectra and days in milk (DIM) in intensively fed dairy cows using statistical prediction methods. Daily BCS was interpolated from cubic splines fitted through the BCS records and daily ΔBCS was calculated from these splines. The ΔBCS records were merged with milk MIR spectra recorded on the same week. The dataset comprised 37,077 ΔBCS phenotypes across 9,403 lactations from 6,988 cows in 151 herds based in Quebec, Canada. Partial least squares regression (PLSR) and a neural network (NN) were then used to estimate ΔBCS from (1) MIR spectra only, (2) DIM only, or (3) MIR spectra and DIM together. The ΔBCS data in both the first 120 and 305 DIM of lactation were used to develop the estimates. Daily ΔBCS had a standard deviation of 4.40 × 10-3 BCS units in the 120-d dataset and of 3.63 × 10-3 BCS units in the 305-d dataset. A 4-fold cross-validation was used to calibrate and test the prediction equations. External validation was also conducted using more recent years of data. Irrespective of whether based on the first 120 or 305 DIM, or when MIR spectra only, DIM only or MIR spectra and DIM were jointly used as prediction variables, NN produced the lowest root mean square error (RMSE) of cross-validation (1.81 × 10-3 BCS units and 1.51 × 10-3 BCS units, respectively, using the 120-d and 305-d dataset). Relative to predictions for the entire 305 DIM, the RMSE of cross-validation was 15.4% and 1.5% lower in the first 120 DIM when using PLSR and NN, respectively. Predictions from DIM only were more accurate than those using just MIR spectra data but, irrespective of the dataset and of the prediction model used, combining DIM information with MIR spectral data as prediction variables reduced the RMSE compared with the inclusion of DIM alone, albeit the benefit was small (the RMSE from cross-validation reduced by up to 5.5% when DIM and spectral data were jointly used as model features instead of DIM only). However, when predicting extreme ΔBCS records, the MIR spectral data were more informative than DIM. Model performance when predicting ΔBCS records in future years was similar to that from cross-validation demonstrating the ability of MIR spectra of milk and DIM combined to estimate ΔBCS, particularly in early lactation. This can be used to routinely generate estimates of ΔBCS to aid in day-to-day individual cow management.

Estimation of body condition score change in dairy cows in a seasonal calving pasture-based system using routinely available milk mid-infrared spectra and machine learning techniques.

Body condition score (BCS) is a subjective estimate of body reserves in cows. Body condition score and its change in early lactation have been associated with cow fertility and health. The aim of the present study was to estimate change in BCS (ΔBCS) using mid-infrared spectra of the milk, with a particular focus on estimating ΔBCS in cows losing BCS at the fastest rate (i.e., the cows most of interest to the producer). A total of 73,193 BCS records (scale 1 to 5) from 6,572 cows were recorded. Daily BCS was interpolated from cubic splines fitted through the BCS records, and subsequently used to calculate daily ΔBCS. Body condition score change records were merged with milk mid-infrared spectra recorded on the same week. Both morning (a.m.) and evening (p.m.) spectra were available. Two different statistical methods were used to estimate ΔBCS: partial least squares regression and a neural network (NN). Several combinations of variables were included as model features, such as days in milk (DIM) only, a.m. spectra only and DIM, p.m. spectra only and DIM, and a.m. and p.m. spectra as well as DIM. The data used to estimate ΔBCS were either based on the first 120 DIM or all 305 DIM. Daily ΔBCS had a standard deviation of 1.65 × 10-3 BCS units in the 305 DIM data set and of 1.98 × 10-3 BCS units in the 120 DIM data set. Each data set was divided into 4 sub-data sets, 3 of which were used for training the prediction model and the fourth to test it. This process was repeated until all the sub-data sets were considered as the test data set once. Using all 305 DIM, the lowest root mean square error of validation (RMSEV; 0.96 × 10-3 BCS units) and the strongest correlation between actual and estimated ΔBCS (0.82) was achieved with NN using a.m. and p.m. spectra and DIM. Using the 120 DIM data, the lowest RMSEV (0.98 × 10-3 BCS units) and the strongest correlation between actual and estimated ΔBCS (0.87) was achieved with NN using DIM and either a.m. spectra only or a.m. and p.m. spectra together. The RMSEV for records in the lowest 2.5% ΔBCS percentile per DIM in early lactation was reduced up to a maximum of 13% when spectra and DIM were both considered in the model compared with a model that considered just DIM. The performance of the NN using DIM and a.m. spectra only with the 120 DIM data was robust across different strata of farm, parity, year of sampling, and breed. Results from the present study demonstrate the ability of mid-infrared spectra of milk coupled with machine learning techniques to estimate ΔBCS; specifically, the inclusion of spectral data reduced the RMSEV over and above using DIM alone, particularly for cows losing BCS at the fastest rate. This approach can be used to routinely generate estimates of ΔBCS that can subsequently be used for farm decisions.

Early detection of subclinical mastitis in lactating dairy cows using cow-level features.

Subclinical mastitis in cows affects their health, well-being, longevity, and performance, leading to reduced productivity and profit. Early prediction of subclinical mastitis can enable dairy farmers to perform interventions to mitigate its effect. The present study investigated how well predictive models built using machine learning techniques can detect subclinical mastitis up to 7 d before its occurrence. The data set used consisted of 1,346,207 milk-day (i.e., a day when milk was collected on both morning and evening) records spanning 9 yr from 2,389 cows producing on 7 Irish research farms. Individual cow composite milk yield and maximum milk flow were available twice daily, whereas milk composition (i.e., fat, lactose, protein) and somatic cell count (SCC) were collected once per week. Other features describing parity, calving dates, predicted transmitting ability for SCC, body weight, and history of subclinical mastitis were also available. The results of the study showed that a gradient boosting machine model trained to predict the onset of subclinical mastitis 7 d before a subclinical case occurs achieved a sensitivity and specificity of 69.45 and 95.64%, respectively. Reduced data collection frequency, where milk composition and SCC were recorded only every 15, 30, 45, and 60 d was simulated by masking data, to reflect the frequency of recording of this data on commercial dairy farms in Ireland. The sensitivity and specificity scores reduced as recording frequency reduced with respective scores of 66.93 and 80.43% when milk composition and SCC were recorded just every 60 d. Results demonstrate that models built on data that could be recorded routinely available on commercial dairy farms, can achieve useful predictive ability of subclinical mastitis even with reduced frequency of milk composition and SCC recording.

Cow-level factors associated with nitrogen utilization in grazing dairy cows using a cross-sectional analysis of a large database.

Reducing nitrogen pollution while maintaining milk production is a major challenge of dairy production. One of the keys to delivering on this challenge is to improve the efficiency of how dairy cows use nitrogen. Thus, estimating the nitrogen utilization of lactating grazing dairy cows and exploring the association between animal factors and productivity with nitrogen utilization are the first steps to understanding the nitrogen utilization complex in dairy cows. Nitrogen utilization metrics were derived from milk and body weight records from 1,291 grazing dairy cows of multiple breeds and crossbreeds; all cows had sporadic information on nitrogen intake concurrent with information on nitrogen sinks (and other nitrogen sources, such as body tissue mobilization). Several nitrogen utilization metrics were investigated, including nitrogen use efficiency (nitrogen output as products such as milk and meat divided by nitrogen intake) and nitrogen excreted (nitrogen intake less the nitrogen output as products such as milk and meat). In the present study, a primiparous Holstein-Friesian used, on average, 20.6% of the nitrogen it ate, excreting the surplus as feces and urine, representing 402 g of nitrogen per day. Intercow variability existed, with a between-cow standard deviation of 0.0094 for nitrogen use efficiency and 24 g of nitrogen per day for nitrogen excretion. As lactation progressed, nitrogen use efficiency declined and nitrogen excretion increased. Nevertheless, nitrogen use efficiency improved (i.e., decreased) from first to second parity, even though it did not improve from second to third parity or greater. Furthermore, nitrogen excretion continued to increase from first to third parity or greater. Nitrogen use efficiency and nitrogen excretion were negatively correlated (-0.56 to -0.40), signifying that dairy cows who partition more of the ingested nitrogen into products such as milk and meat, on average, also excrete less nitrogen. Milk urea nitrogen was, at best, weakly correlated with nitrogen use efficiency and nitrogen excretion; the correlations were between -0.01 and 0.06. In conclusion, several cow-level factors such as parity, stage of lactation, and breed were associated with the range of different nitrogen efficiency metrics investigated; moreover, even after accounting for such effects, 4.8% to 6.3% of the remaining variation in the nitrogen use efficiency and nitrogen balance metrics were attributable to intercow differences.

Re-assessing the importance of linear type traits in predicting genetic merit for survival in an aging Holstein-Friesian dairy cow population.

The cumulative improvement achieved in the genetic merit for reproductive performance in dairy populations will likely improve dairy cow longevity; therefore, it is time to reassess whether linear type traits are still suitable predictors of survival in an aging dairy cow population. The objective of the present study was therefore to estimate the genetic correlations between linear type traits and survival from one parity to the next and, in doing so, evaluate if those genetic correlations change with advancing parity. After edits, 152,894 lactation survival records (first to ninth parity) were available from 52,447 Holstein-Friesian cows, along with linear type trait records from 52,121 Holstein-Friesian cows. A series of bivariate random regression models were used to estimate the genetic covariances between survival in different parities and each linear type trait. Heritability estimates for survival per parity ranged from 0.02 (SE = 0.004; first parity) to 0.05 (SE = 0.01; ninth parity). Pairwise genetic correlations between survival among different parities varied from 0.42 (first and ninth parity) to 1.00 (eighth to ninth parity), with the strength of these genetic correlations being inversely related to the interval between the compared parities. The genetic correlations between survival and the individual linear type traits varied across parities for 9 of the 20 linear type traits examined, but the correlations with only 3 of these linear type traits strengthened as the cows aged; these 3 traits were rear udder height, teat length, and udder depth. Given that linear type traits are frequently scored in first parity and are genetically correlated with survival in older parities, they may be suitable early predictors of survival, especially for later parity cows. Additionally, the direction of the genetic correlations between survival and rear udder height, teat length, and udder depth did not change between parities; hence, selection for survival in older parities using these linear type traits should not hinder genetic improvement for survival in younger parities.

Exploiting genetic variability in the trajectory of lactation yield and somatic cell score with each progressing parity.

The inclusion of reproductive performance in dairy cow breeding schemes has resulted in a cumulative improvement in genetic merit for reproductive performance; this improvement should manifest in longer productive lives through a reduced requirement for involuntary culling. Nonetheless, the average length of dairy cow productive life has not changed in most populations, suggesting that risk factors for culling, especially in older cows, are possibly more associated with lower yield or high somatic cell score (SCS) than compromised reproductive performance. The objective of the present study was to understand the dynamics of lactation yields and SCS in dairy cows across parities and, in doing so, quantify the potential to alter this trajectory through breeding. After edits, 3,470,520 305-d milk, fat, and protein yields, as well as milk fat and protein percentage and somatic cell count records from 1,162,473 dairy cows were available for analysis. Random regression animal models were used to identify the parity in which individual cows reached their maximum lactation yields, and highest average milk composition and SCS; also estimated from these models were the (co)variance components for yield, composition, and SCS per parity across parities. Estimated breeding values for all traits per parity were calculated for cows reaching ≥fifth parity. Of the cows included in the analyses, 91.0%, 92.2%, and 83.4% reached maximum milk, fat, and protein yield in fifth parity, respectively. Conversely, 95.9% of cows reached their highest average fat percentage in first parity and 62.9% of cows reached their highest average protein percentage in third parity. In contrast to both milk yield and composition traits, 98.4% of cows reached their highest average SCS in eighth parity. Individual parity estimates of heritability for milk yield traits, milk composition, and SCS ranged from 0.28 to 0.44, 0.47 to 0.69, and 0.13 to 0.23, respectively. The strength of the genetic correlations per trait among parities was inversely related to the interval between the parities compared; the weakest genetic correlation was 0.67 (standard error = 0.02) between milk yield in parities 1 and 8. Eigenvalues and eigenfunctions of the additive genetic covariance matrices for all investigated traits revealed potential to alter the trajectory of parity profiles for milk yield, milk composition, and SCS. This was further demonstrated when evaluating the trajectories of animal estimated breeding values per parity.

Are subjectively scored linear type traits suitable predictors of the genetic merit for feed intake in grazing Holstein-Friesian dairy cows?

Measuring dry matter intake (DMI) in grazing dairy cows using currently available techniques is invasive, time consuming, and expensive. An alternative to directly measuring DMI for use in genetic evaluations is to identify a set of readily available animal features that can be used in a multitrait genetic evaluation for DMI. The objectives of the present study were thus to estimate the genetic correlations between readily available body-related linear type traits and DMI in grazing lactating Holstein-Friesian cows, but importantly also estimate the partial genetic correlations between these linear traits and DMI, after adjusting for differences in genetic merit for body weight. Also of interest was whether the predictive ability derived from the estimated genetic correlations materialized upon validation. After edits, a total of 8,055 test-day records of DMI, body weight, and milk yield from 1,331 Holstein-Friesian cows were available, as were chest width, body depth, and stature from 47,141 first lactation Holstein-Friesian cows. In addition to considering the routinely recorded linear type traits individually, novel composite traits were defined as the product of the linear type traits as an approximation of rumen volume. All linear type traits were moderately heritable, with heritability estimates ranging from 0.27 (standard error = 0.14) to 0.49 (standard error = 0.15); furthermore, all linear type traits were genetically correlated (0.29 to 0.63, standard error 0.14 to 0.12) with DMI. The genetic correlations between the individual linear type traits and DMI, when adjusted for genetic differences in body weight, varied from -0.51 (stature) to 0.48 (chest width). These genetic correlations between DMI and linear type traits suggest linear type traits may be useful predictors of DMI, even when body weight information is available. Nonetheless, estimated genetic merit of DMI derived from a multitrait genetic evaluation of linear type traits did not correlate strongly with actual DMI in a set of validation animals; the benefit was even less if body weight data were also available.

Concordance rate in cattle and sheep between genotypes differing in Illumina GenCall quality score.

Proper quality control of data prior to downstream analyses is fundamental to ensure integrity of results; quality control of genomic data is no exception. While many metrics of quality control of genomic data exist, the objective of the present study was to quantify the genotype and allele concordance rate between called single nucleotide polymorphism (SNP) genotypes differing in GenCall (GC) score; the GC score is a confidence measure assigned to each Illumina genotype call. This objective was achieved using Illumina beadchip genotype data from 771 cattle (12 428 767 genotypes in total post-editing) and 80 sheep (1 557 360 SNPs genotypes in total post-editing) each genotyped in duplicate. The called genotype with the lowest associated GC score was compared to the genotype called for the same SNP in the same duplicated animal sample but with a GC score of >0.90 (assumed to represent the true genotype). The mean genotype concordance rate for a GC score of <0.300, 0.300-0.549, and ≥0.550 in the cattle (sheep in parenthesis) was 0.9467 (0.9864), 0.9707 (0.9953), and 0.9994 (0.99997) respectively; the respective allele concordance rate was 0.9730 (0.9930), 0.9849 (0.9976), and 0.9997 (0.99998). Hence, concordance eroded as the GC score of the called genotype reduced, albeit the impact was not dramatic and was not very noticeable until a GC score of <0.55. Moreover, the impact was greater and more consistent in the cattle population than in the sheep population. Furthermore, an impact of GC score on genotype concordance rate existed even for the same SNP GenTrain value; the GenTrain value is a statistical score that depicts the shape of the genotype clusters and the relative distance between the called genotype clusters.

Invited review: Beef-on-dairy-The generation of crossbred beef × dairy cattle.

Because a growing proportion of the beef output in many countries originates from dairy herds, the most critical decisions about the genetic merit of most carcasses harvested are being made by dairy producers. Interest in the generation of more valuable calves from dairy females is intensifying, and the most likely vehicle is the use of appropriately selected beef bulls for mating to the dairy females. This is especially true given the growing potential to undertake more beef × dairy matings as herd metrics improve (e.g., reproductive performance) and technological advances are more widely adopted (e.g., sexed semen). Clear breed differences (among beef breeds but also compared with dairy breeds) exist for a whole plethora of performance traits, but considerable within-breed variability has also been demonstrated. Although such variability has implications for the choice of bull to mate to dairy females, the fact that dairy females themselves exhibit such genetic variability implies that "one size fits all" may not be appropriate for bull selection. Although differences in a whole series of key performance indicators have been documented between beef and beef-on-dairy animals, of particular note is the reported lower environmental hoofprint associated with beef-on-dairy production systems if the environmental overhead of the mature cow is attributed to the milk she eventually produces. Despite the known contribution of beef (i.e., both surplus calves and cull cows) to the overall gross output of most dairy herds globally, and the fact that each dairy female contributes half her genetic merit to her progeny, proxies for meat yield (i.e., veal or beef) are not directly considered in the vast majority of dairy cow breeding objectives. Breeding objectives to identify beef bulls suitable for dairy production systems are now being developed and validated, demonstrating the financial benefit of using such breeding objectives over and above a focus on dairy bulls or easy-calving, short-gestation beef bulls. When this approach is complemented by management-based decision-support tools, considerable potential exists to improve the profitability and sustainability of modern dairy production systems by exploiting beef-on-dairy breeding strategies using the most appropriate beef bulls.

Contribution of genetic variability to phenotypic differences in on-farm efficiency metrics of dairy cows based on body weight and milk solids yield.

Milk solids per kilogram of body weight (BW) is growing in popularity as a measure of dairy cow lactation efficiency. Little is known on the extent of genetic variability that exist in this trait but also the direction and strength of genetic correlations with other performance traits. Such genetic correlations are important to know if producers are to consider actively selecting cows excelling in milk solids per kilogram of BW. The objective of the present study was to use a large data set of commercial Irish dairy cows to quantify the extent of genetic variability in milk solids per kilogram of BW and related traits but also their genetic and phenotypic inter-relationships. Mid-lactation BW and body condition score (BCS), along with 305-d milk solids yield (i.e., fat plus protein yield) were available on 12,413 lactations from 11,062 cows in 85 different commercial dairy herds. (Co)variance components were estimated using repeatability animal linear mixed models. The genetic correlation between milk solids and body weight was only 0.05, which when coupled with the observed large genetic variability in both traits, indicate massive potential to select for both traits in opposite directions. The genetic correlations between both milk solids and BW with BCS; however, need to be considered in any breeding strategy. The genetic standard deviation, heritability, and repeatability of milk solids per kilogram of BW was 0.08, 0.37, and 0.57, respectively. The genetic correlation between milk solids per kilogram of BW with milk solids, BW, and BCS was 0.62, -0.75, and -0.41, respectively. Therefore, based on genetic regression, each increase of 0.10 units in genetic merit for milk solids per kilogram of BW is expected to result in, on average, an increase in 16.1 kg 305-d milk solids yield, a reduction of 25.6 kg of BW and a reduction of 0.05 BCS units (scale of 1-5 where 1 is emaciated). The genetic standard deviation (heritability) for 305-d milk solids yield adjusted phenotypically to a common BW was 27.3 kg (0.22). The genetic correlation between this adjusted milk solids trait with milk solids, BW, and BCS was 0.91, -0.12, and -0.26, respectively. Once also adjusted phenotypically to a common BCS, the genetic standard deviation (heritability) for milk solids adjusted phenotypically to a common BW was 26.8 kg (0.22) where the genetic correlation with milk solids, BW and BCS was 0.91, -0.21, and -0.07, respectively. The genetic standard deviation (heritability) of BW adjusted phenotypically for differences in milk solids was 35.3 kg (0.61), which reduced to 33.2 kg when also phenotypically adjusted for differences in BCS. Results suggest considerable opportunity exists to change milk solids yield independent of BW, and vice versa. The opportunity is reduced slightly once also corrected for differences in BCS. Inter-animal BCS differences should be considered if selection on such metrics is contemplated.

Cross-sectional analyses of a national database to determine if superior genetic merit translates to superior dairy cow performance.

Various studies have validated that genetic divergence in dairy cattle translates to phenotypic differences; nonetheless, many studies that consider the breeding goal, or associated traits, have generally been small scale, often undertaken in controlled environments, and they lack consideration for the entire suite of traits included in the breeding goal. Therefore, the objective of the present study was to fill this void, and in doing so, provide producers with confidence that the estimated breeding values (EBV) included in the breeding goal do (or otherwise) translate to desired changes in performance among commercial cattle; an additional outcome of such an approach is the identification of potential areas for improvements. Performance data on 536,923 Irish dairy cows (and their progeny) from 13,399 commercial spring-calving herds were used. Association analyses between the cow's EBV of each trait included in the Irish total merit index for dairy cows (which was derived before her own performance data accumulated) and her subsequent performance were undertaken using linear mixed models; milk production, fertility, calving, maintenance (i.e., liveweight), beef, health, and management traits were all considered in the analyses. Results confirm that excelling in EBV for individual traits, as well as on the total merit index, generally delivers superior phenotypic performance; examples of the improved performance for genetically elite animals include a greater yield and concentration of both milk fat and milk protein, despite a lower milk volume, superior reproductive performance, better survival, improved udder and hoof health, lighter cows, and fewer calving complications; all these gains were achieved with minimal to no effect on the beef merit of the dairy cow's progeny. The associated phenotypic change in each performance trait per unit change in its respective EBV was largely in line with the direction and magnitude of expectation, the exception being for calving interval. Per unit change in calving interval EBV, the direction of phenotypic response was as anticipated but the magnitude of the response was only half of what was expected. Despite the deviation from expectation between the calving interval EBV and its associated phenotype, a superior total merit index or a superior fertility EBV was indeed associated with an improvement in all detailed fertility performance phenotypes investigated. Results substantiate that breeding is a sustainable strategy of improving phenotypic performance in commercial dairy cattle and, by extension, profit.

Prediction of genetic merit for live weight and body condition score in dairy cows using routinely available linear type and carcass data.

Accurate estimates of genetic merit for both live weight and body condition score (BCS) could be useful additions to both national- and herd-breeding programs. Although recording live weight and BCS is not technologically arduous, data available for use in routine genetic evaluations are generally lacking. The objective of the present study was to explore the usefulness of routinely recorded data, namely linear type traits (which also included BCS but only assessed visually) and carcass traits in the pursuit of genetic evaluations for both live weight and BCS in dairy cows. The data consisted of on-farm records of live weight and BCS (assessed using both visual and tactile cues) from 33,242 dairy cows in 201 commercial Irish herds. These data were complemented with information on 6 body-related linear type traits (i.e., stature, angularity, chest width, body depth, BCS, and rump width) and 3 cull cow carcass measures (i.e., carcass weight, conformation, and fat cover) on a selection of these animals plus close relatives. (Co)variance components were estimated using animal linear mixed models. The genetic correlation between the type traits stature, angularity, body depth, chest width, rump width, and visually-assessed BCS with live weight was 0.68, -0.28, 0.43, 0.64, 0.61, and 0.44, respectively. The genetic correlation between angularity and BCS measured on farm (based on both visual and tactile appraisal) was -0.79; the genetic and phenotypic correlation between BCS assessed visually as part of the linear assessment with BCS assessed by producers using both tactile and visual cues was 0.90 and 0.27, respectively. The genetic (phenotypic) correlation between cull cow carcass weight and live weight was 0.81 (0.21), and the genetic (phenotypic) correlation between cull cow carcass fat cover and BCS assessed on live cows was 0.44 (0.12). Estimated breeding values (EBV) for live weight and BCS in a validation population of cows were generated using a multitrait evaluation with observations for just the type traits, just the carcass traits, and both the type traits and carcass traits; the EBV were compared with the respective live weight and BCS phenotypic observations. The regression of phenotypic live weight on its EBV from the multitrait evaluations was 1.00 (i.e., the expectation) when the EBV was generated using just linear type trait data, but less than 1 (0.83) when using just carcass data. However, the regression changed across parities and stages of lactation. The partial correlation (after adjusting for contemporary group, parity by stage of lactation, heterosis, and recombination loss) between phenotypic live weight and EBV for live weight estimated using the 3 different scenarios (i.e., type only, carcass only, type plus carcass) ranged from 0.38 to 0.43. Although the prediction of phenotypic BCS from its respective EBV was relatively good when using just the linear type trait data (regression coefficient of 0.83 with a partial correlation of 0.22), the predictive ability of BCS EBV based on just carcass data was poor and should not be used. Overall, linear type trait data are a useful source of information to predict live weight and BCS with minimal additional predictive value from also including carcass data. Nonetheless, in the absence of linear type trait data, information on carcass traits can be useful in predicting genetic merit for mature cow live weight. Prediction of cow BCS from cow carcass data is not recommended.

Genetic and nongenetic factors associated with lactation length in seasonal-calving, pasture-based dairy cows.

Lactation yield estimates standardized to common lactation lengths of 270-d or 305-d equivalents are commonly used in management decision support tools and dairy cow genetic evaluations. The use of such measurements to quantify the (genetic) merit of individual cows fails to penalize cows that do not reach the standardized lactation length, or indeed reward cows that lactate for more than the standardized lactation length. The objective of the present study was to quantify the genetic and nongenetic factors associated with lactation length in seasonal-calving, pasture-based dairy cows. A total of 616,350 lactation length records from 285,598 Irish cows were used. Linear mixed models were used to quantify the associations between lactation length and calving month, parity, age at calving, previous dry period length, calving difficulty score, heterosis, recombination loss, breed, and herd size, as well as to estimate the genetic and residual variance components of lactation length. The median lactation length in the edited data set was 288 d, with 27% of cows achieving lactations of at least 305 d. Relative to cows calving in January, the lactations of cow calving in February, March, or April was, on average, 4.2, 12.7, and 21.9 d shorter, respectively. The lactation length of a first parity cow was, on average, 7.8, 8.6, and 8.4 d shorter than that of second, third, and fourth parity cows, respectively. Norwegian Red and Montbéliarde cows had, on average, a 4.7- and 1.6-d shorter lactation than Holstein-Friesian cows, respectively. The heritability estimate, coefficient of genetic variation, and repeatability estimate of lactation length were 0.02, 1.2%, and 0.04, respectively. Based on the genetic standard deviation for lactation length estimated in the present study (3.3 d), cows ranked in the top 20% for genetic merit for lactation length would be expected to have lactations 9.2 d longer than cows in the bottom 20%, demonstrating exploitable genetic variability. Given the vast array of genetic and nongenetic factors associated with lactation length, an approach which combines improved management practices and selective breeding may be an efficient and effective strategy to lengthen lactations.

Predicting cow milk quality traits from routinely available milk spectra using statistical machine learning methods.

Numerous statistical machine learning methods suitable for application to highly correlated features, as those that exist for spectral data, could potentially improve prediction performance over the commonly used partial least squares approach. Milk samples from 622 individual cows with known detailed protein composition and technological trait data accompanied by mid-infrared spectra were available to assess the predictive ability of different regression and classification algorithms. The regression-based approaches were partial least squares regression (PLSR), ridge regression (RR), least absolute shrinkage and selection operator (LASSO), elastic net, principal component regression, projection pursuit regression, spike and slab regression, random forests, boosting decision trees, neural networks (NN), and a post-hoc approach of model averaging (MA). Several classification methods (i.e., partial least squares discriminant analysis (PLSDA), random forests, boosting decision trees, and support vector machines (SVM)) were also used after stratifying the traits of interest into categories. In the regression analyses, MA was the best prediction method for 6 of the 14 traits investigated [curd firmness at 60 min, αS1-casein (CN), αS2-CN, κ-CN, α-lactalbumin, and ß-lactoglobulin B], whereas NN and RR were the best algorithms for 3 traits each (rennet coagulation time, curd-firming time, and heat stability, and curd firmness at 30 min, ß-CN, and ß-lactoglobulin A, respectively), PLSR was best for pH, and LASSO was best for CN micelle size. When traits were divided into 2 classes, SVM had the greatest accuracy for the majority of the traits investigated. Although the well-established PLSR-based method performed competitively, the application of statistical machine learning methods for regression analyses reduced the root mean square error compared with PLSR from between 0.18% (κ-CN) to 3.67% (heat stability). The use of modern statistical machine learning methods for trait prediction from mid-infrared spectroscopy may improve the prediction accuracy for some traits.