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.

Comparison of greenhouse gas emissions from sheep measured using both respiration and portable accumulation chambers

Methane (CH4) is a potent gas produced by ruminants, and new measurement techniques are required to generate large datasets suitable for genetic analysis. One such technique are portable accumulation chambers (PAC), a short-term sampling method. The objectives of the current study were to explore the relationship between CH4 and carbon dioxide (CO2) output measured using both PAC and respiration chambers (RC) in growing lambs, and separately investigate the relationship among CH4, CO2 and measured ad libitum DM intake (DMI). Methane, CO2 and DMI were measured on 30 Suffolk and 30 Texel ewe lambs (age 253 ± 12 days) using the RC and PAC sequentially. The experiment was conducted over a 14-day period, with DMI measured from days 1 to 14; measurements in RC were conducted from days 10 to 12, while measurements in PAC were taken twice, the day immediately prior to the lambs entering the RC (day 9; PAC Pre-RC) and on the day lambs exited the RC (day 13; PAC Post-RC). Greater CH4 and CO2 output was measured in the RC than in the PAC (P < 0.01); similarly mean CH4 yield was greater when measured in the RC (15.39 ± 0.452 g CH4/kg DMI) compared to PAC (8.01 ± 0.767 g CH4/kg DMI). A moderate correlation of 0.37 was found between CH4 output measured in PAC Pre-RC and the RC, the corresponding regression coefficient of CH4 output measured in the RC regressed on CH4 output measured in PAC Pre-RC was close to unity (0.74; SE 0.224). The variance of CH4 and CO2 output within the measurement technique did not differ from each other (P > 0.05). Moderate to strong correlations were found between CH4 and CO2 per kg of live weight and CH4 and CO2 yield. Results from this study highlight the suitability of PAC as a ranking tool to rank animals based on their gaseous output when compared to the RC. However, repeated measurements separated by several days may be beneficial if precise rankings are required. Given the close to unity regression coefficient of CH4 output measured in the RC regressed on CH4 output measured in PAC Pre-RC suggests that PAC could also be potentially used to estimate absolute CH4 output; however, further research is required to substantiate this claim. When DMI is unknown, CH4 and CO2 per kg of live weight are a suitable alternative to the measurement of CH4 and CO2 yield.

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.

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.

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.

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.

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.

Animal board invited review: Practical applications of genomic information in livestock

Access to high-dimensional genomic information in many livestock species is accelerating. This has been greatly aided not only by continual reductions in genotyping costs but also an expansion in the services available that leverage genomic information to create a greater return-on-investment. Genomic information on individual animals has many uses including (1) parentage verification and discovery, (2) traceability, (3) karyotyping, (4) sex determination, (5) reporting and monitoring of mutations conferring major effects or congenital defects, (6) better estimating inbreeding of individuals and coancestry among individuals, (7) mating advice, (8) determining breed composition, (9) enabling precision management, and (10) genomic evaluations; genomic evaluations exploit genome-wide genotype information to improve the accuracy of predicting an animal's (and by extension its progeny's) genetic merit. Genomic data also provide a huge resource for research, albeit the outcome from this research, if successful, should eventually be realised through one of the ten applications already mentioned. The process for generating a genotype all the way from sample procurement to identifying erroneous genotypes is described, as are the steps that should be considered when developing a bespoke genotyping panel for practical application.

Genetic (co)variance components for slaughter traits in a multi-breed sheep population

Carcass value is one of the main contributors to revenue in meat sheep enterprises, while age at slaughter is also a major component to the cost of production. Despite the contribution of such traits to overall profit, little is actually known on the extent of exploitable genetic variability in the traits that govern carcass value (i.e. carcass weight, carcass conformation, carcass fat) and age at slaughter, especially independent of each other. The objective of the present study was to estimate genetic (co)variances for and among carcass weight, carcass conformation, carcass fat, kill-out percentage and age at slaughter as well as their genetic (co)variances with traits measured earlier in life. Data consisted of slaughter records from 15 714 lambs, with 12 630 of these lambs having at least one live weight measure. The heritability (SE) of carcass weight, carcass conformation, carcass fat, kill-out percentage, and age at slaughter was 0.14 (0.02), 0.19 (0.02), 0.08 (0.01), 0.22 (0.03), and 0.16 (0.02), respectively. The maternal heritability for age at slaughter was 0.07 (0.02); no maternal genetic influence was found on any of the other slaughter traits. The coefficient of genetic variation for carcass weight and age at slaughter was 3 and 8%, respectively. The correlations between the direct genetic effects for live weight throughout life, and carcass weight were weak up to weaning but were strong (0.83) thereafter. The correlation between the direct genetic effects of birth weight and age at slaughter was zero, but varied from -0.91 to -0.56 between live weight measured later in life and age at slaughter. Results demonstrate significant exploitable genetic variability in a range of slaughter traits with the prediction of genetic merit for carcass traits and age at slaughter being possible using live weight measures taken on live animals. For example, the accuracy of selection for slaughter traits (comprising of age at slaughter, carcass conformation and carcass fat) from weaning weight records available on 100 progeny was 0.37; when slaughter data were also available for 10 progeny, the accuracy of selection increased to 0.56.

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.

Evaluating the use of statistical and machine learning methods for estimating breed composition of purebred and crossbred animals in thirteen cattle breeds using genomic information

Introduction: The ability to accurately predict breed composition using genomic information has many potential uses including increasing the accuracy of genetic evaluations, optimising mating plans and as a parameter for genotype quality control. The objective of the present study was to use a database of genotyped purebred and crossbred cattle to compare breed composition predictions using a freely available software, Admixture, with those from a single nucleotide polymorphism Best Linear Unbiased Prediction (SNP-BLUP) approach; a supplementary objective was to determine the accuracy and general robustness of low-density genotype panels for predicting breed composition. Methods: All animals had genotype information on 49,213 autosomal single nucleotide polymorphism (SNPs). Thirteen breeds were included in the analysis and 500 purebred animals per breed were used to establish the breed training populations. Accuracy of breed composition prediction was determined using a separate validation population of 3,146 verified purebred and 4,330 two and three-way crossbred cattle. Results: When all 49,213 autosomal SNPs were used for breed prediction, a minimal absolute mean difference of 0.04 between Admixture vs. SNP-BLUP breed predictions was evident. For crossbreds, the average absolute difference in breed prediction estimates generated using SNP-BLUP and Admixture was 0.068 with a root mean square error of 0.08. Breed predictions from low-density SNP panels were generated using both SNP-BLUP and Admixture and compared to breed prediction estimates using all 49,213 SNPs (representing the gold standard). Breed composition estimates of crossbreds required more SNPs than predicting the breed composition of purebreds. SNP-BLUP required ≥3,000 SNPs to predict crossbred breed composition, but only 2,000 SNPs were required to predict purebred breed status. The absolute mean (standard deviation) difference across all panels <2,000 SNPs was 0.091 (0.054) and 0.315 (0.316) when predicting the breed composition of all animals using Admixture and SNP-BLUP, respectively compared to the gold standard prediction. Discussion: Nevertheless, a negligible absolute mean (standard deviation) difference of 0.009 (0.123) in breed prediction existed between SNP-BLUP and Admixture once ≥3,000 SNPs were considered, indicating that the prediction of breed composition could be readily integrated into SNP-BLUP pipelines used for genomic evaluations thereby avoiding the necessity for a stand-alone software.

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.

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.

Association between the prion protein genotype and animal performance traits in a large multibreed sheep population

Genetic susceptibility to scrapie, a fatal disease of sheep and goats, is modulated by polymorphisms in the prion protein (PrP). Neither the frequency of the PrP genotypes nor their association with animal performance has been investigated in a large multibreed Irish sheep population. Scrapie genotypes were available on 16 416 animals; the breeds represented included purebred Belclare (733), Charollais (333), Suffolk (739), Texel (1 857), Vendeen (191), and crossbreds (12 563). Performance data on lambing, lamb and ewe performance as well as health traits were available. The association between alternative approaches of describing the PrP genotype (i.e. 15 individually called PrP genotypes, five genotype classes representing susceptibility to scrapie, or number of ARR haplotypes) and animal performance were quantified using animal linear mixed models. All 15 of the possible scrapie genotypes were detected, although the frequency differed by breed. The frequency of the five PrP haplotypes in the entire population were 0.70 (ARR), 0.15 (ARQ), 0.11 (ARH), 0.02 (AHQ) and 0.01 (VRQ); the most susceptible haplotype (VRQ) was only detected in purebred Texels and crossbreds. No association was detected between the PrP genotype of either the animal or dam and any of the lambing traits (i.e. lambing difficulty score, perinatal mortality and birth weight). With the exception of ultrasound muscle depth, no association between the PrP genotype and any of the lamb performance traits (i.e. lamb BW and carcass) was observed. Lambs carrying the category four PrP genotype (i.e. ARR/VRQ) had 1.20 (SE = 0.45) mm, 1.38 (SE = 0.12) mm, 1.47 (S = 0.25) mm shallower ultrasound muscle depth relative to lambs of the less susceptible scrapie categories of 1, 2, 3, respectively (P < 0.05). Nonetheless, no association between PrP genotype and lamb carcass conformation, the ultimate end goal of producers, was detected. Ewe litter size, body condition score or lameness did not differ by PrP genotype of the ewe (P > 0.05). For ewe mature BW, ARH/VRQ ewes differed from most other ewe PrP genotypes and were, on average, 3.79 (SE = 1.66) kg heavier than ARR/ARR genotype ewes. Lamb dag score differed by dam PrP genotype (P < 0.05), although the differences were small. Results from this study show that scrapie is segregating within the Irish sheep population, but the PrP genotype was not associated with most traits investigated and, where associations were detected, the biological significance was minimal. This suggests minimal impact of selection on PrP genotype on performance, at least for the traits investigated in the present study.

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.

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.

Commercial beef farms excelling in terminal and maternal genetic merit generate more gross profit

Validation of beef total merit breeding indexes for improving performance and profitability has previously been undertaken at the individual animal level; however, no herd-level validation of beef genetic merit and profit has been previously investigated. The objective of the present study was to quantify the relationship between herd profitability and both herd-average terminal and maternal genetic merit across 1,311 commercial Irish beef herds. Herd-level physical and financial performance data were available from a financial benchmarking tool used by Irish farmers and their extension advisors. Animal genetic merit data originated from the Irish Cattle Breeding Federation who undertake the national beef and dairy genetic evaluations. Herd-average genetic merit variables included the terminal index of young animals, the maternal index of dams, and the terminal index of service sires. The herds represented three production systems: 1) cow-calf to beef, 2) cow-calf to weanling/yearling, and 3) weanling/yearling to beef. Associations between herd financial performance metrics and herd average genetic merit variables were quantified using a series of linear mixed models with year, production system, herd size, stocking rate, concentrate input, and the two-way interactions between production system and herd size, stocking rate, and concentrate input included as nuisance factors. Herd nested within the county of Ireland (n = 26) was included as a repeated effect. Herds with young cattle excelling in terminal index enjoyed greater gross and net profit per hectare (ha), per livestock unit (LU), and per kg net live-weight output. The change in gross profit per LU per unit change in the terminal index of young animals was €1.41 (SE = 0.23), while the respective regression coefficient for net profit per LU was €1.37 (SE = 0.30); the standard deviation of the terminal index is €37. Herd-average dam maternal index and sire terminal index were both independently positively associated with gross profit per ha and gross profit per LU. Each one unit increase in dam maternal index (standard deviation of €38) was associated with a €1.40 (SE = 0.48) and €0.76 (SE = 0.29) greater gross profit per ha and per LU, respectively. Results from the present study at the herd-level concur with previous validation studies at the individual animal level thus instilling further confidence among stakeholders as to the expected improvement in herd profitability with improving genetic merit.

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.

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.

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.

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.

Short communication: Differences in genetic merit for visually-assessed body condition score materialises as phenotypic differences in tactile-based body condition score in commercial dairy cows

Body condition score (BCS) is a known risk factor for cow health and well-being. Many different BCS scales and systems for assessment exist;while the scales used for assessing BCS vary, differences in how BCS is assessed (i.e., visual versus visual plus tactile) and the extent of training and experience of the assessor (i.e., professional linear classifiers versus producers) also contributes to the underlying variability. Registered dairy cows globally are routinely assessed for linear type traits which describe biological extremes in the morphological attributes; BCS and a correlated trait angularity are within this suite of traits assessed. These linear-type data are used to generate estimates of genetic merit (predicted transmitting ability), but how these estimates manifest themselves as phenotypic differences when assessed by producers on commercial multiparous cows has never been quantified. To evaluate this, 58440 phenotypic BCS records from 48823 lactations in 38608 cows were used. Associations were undertaken using linear mixed models relating phenotypic BCS to genetic merit after accounting for nuisance factors. Differences in genetic merit for either BCS or angularity (assessed visually by professionals on a 1 to 9 scale just once during lactation in primiparous registered cows) translated to phenotypic difference in BCS (assessed by producers using both tactile and visual assessment on a 1 to 5 scale across lactation in commercial dairy cows). The partial correlation between test phenotypic BCS and genetic merit for either BCS or angularity was 0.13 and 0.10, respectively. Based on the model coefficients estimated in the present study, the mean expected difference in phenotypic BCS on a 1 to 5 scale between the top and bottom 10% on genetic merit for BCS or angularity was 0.28 and 0.31 units, respectively. Results from the present study clearly provide confidence that genetic merit for BCS or angularity based on a single visual assessment in primiparous cows is useful to breed for cows of better body condition, irrespective of stage of lactation or parity.

Are Type and Screen Samples Routinely Necessary Before Laparoscopic Cholecystectomy?

AIMS

Type and screen (T&S) samples are routinely requested before all laparoscopic cholecystectomies (LCs) at our centre despite the low reported risk of major vascular injury and peri-operative transfusion. Our retrospective case series aimed to identify local transfusion need to inform policy.

METHODS

Emergency and elective LC performed at a single tertiary centre between March 2014 and October 2016 (30 months) were analysed. This included all patients aged ≥ 16, and procedures converted to open where LC was the primary procedure. Peri-operative complications and transfusion data were obtained from electronic records.

RESULTS

In total, 1002 consecutive patients met inclusion criteria; 12 patients were transfused during index admission (1.20%). No patients required emergency transfusion or had major vascular injuries. Despite local policy, 106 patients (10.6%) did not have a valid T&S sample prior to their procedure. Transfused patients were more likely to be emergency admissions (n = 10/12). The most common indications for transfusion were pre-operative anaemia (n = 7/12) and septic coagulopathy (n = 2/12).

CONCLUSIONS

Peri-operative transfusions at our centre were low. No patients required intra-operative blood transfusions dependent on a pre-operative T&S sample. Patients requiring transfusion were predictable from their pre-operative clinical status. We propose that a highly selective T&S policy is safe and can reduce costs.

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.

Formulation of a decision support tool incorporating both genetic and non-genetic effects to rank young growing cattle on expected market value

While breeding indexes exist globally to identify candidate parents of the next generation, fewer tools exist that provide guidance on the expected monetary value of young animals. The objective of the present study was therefore to develop the framework for a cattle decision-support tool which incorporates both the genetic and non-genetic information of an animal and, in doing so, better predict the potential market value of an animal, whatever the age. Two novel monetary indexes were constructed and their predictive ability of carcass value was compared to that of the Irish national Terminal breeding index, typical of other terminal indexes used globally. A constructed Harvest index was composed of three carcass-related traits [i.e., 1) carcass weight, 2) carcass conformation and 3) carcass fat, each weighted by their respective economic value] and aimed at purchasers of animals close to harvest; the second index, termed the Calf index, also included docility and feed intake (weighted by their respective economic value), thus targeting purchasers of younger calves for growing (and eventually harvesting). Genetic and non-genetic fixed and random effect model solutions from the Irish national genetic evaluations underpinned all indexes. The two novel indexes were formulated using three alternative estimates of an animal's total merit for comparative purposes: 1) an index based solely on the animal's breed solutions, 2) an index which also included within-breed animal differences, and 3) an index which, as well as considering additive and non-additive genetic effects, also included non-genetic effects (referred to as production values [PVs]). As more information (i.e., within breed effects and subsequently non-genetic effects) was included in the total merit estimate, the correlations strengthened between the two proposed indexes and the animal's calculated carcass market value; the correlation coefficients almost doubled in strength when total merit was based on PV-based estimates as compared to the breed solutions alone. Including phenotypic live-weight data, collected during the animal's life, strengthened the predictive ability of the indexes further. Based on the results presented, the proposed indexes may fill the void in decision support when purchasing or selling cattle. In addition, given the dynamic nature of indexes, they have the potential to be updated in real-time as information becomes available.

Inter-animal genetic variability exist in organoleptic properties of prime beef meat

The objective of the present study was to estimate genetic parameters for four organoleptic traits in beef meat, namely tenderness, juiciness, flavour and chewiness using data from 5380 young crossbred progeny of 748 different sires. As well as using the mean animal sensory score across all panellists for a given trait, other aggregate functions such as the median and modal values were also investigated. The heritability (SE) of mean tenderness, juiciness, flavour and chewiness was 0.16 (0.04), 0.14 (0.04), 0.11 (0.03) and 0.21 (0.06), respectively; heritability estimates for the other aggregate values of these traits were generally lower. All genetic correlations between tenderness, juiciness and flavour were positive (0.52 to 0.68) while the genetic correlations between these three traits with chewiness were all negative varying from -0.95 to -0.48. Weak genetic correlations (≤|0.16|) were evident between the sensory traits and all of carcass weight, conformation and subcutaneous fat cover.

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.

Short communication: Animal-level factors associated with whether a dairy female is mated to a dairy or beef bull

When serving a female, the producer must decide whether to mate her to a dairy or beef bull. Tools assisting in this decision could be a useful component of the decision process. A database of 2,283,100 artificial inseminations from 806,725 dairy females was used to investigate what factors were associated with servicing a given female to a beef bull. The probability of being inseminated with a beef bull increased with each service and as the breeding season progressed. An older cow had greater odds of being served with a beef bull, as did cows that calved later in the year, had recently experienced dystocia, were a longer time calved, or were of a poor overall genetic merit compared with herdmates. Cows with low somatic cell count in the previous lactation compared with herdmates were less likely to be mated to a beef bull, as were cows that yielded relatively higher milk solids in the previous lactation. Relative to a first-parity cow, the odds of a fifth-parity cow being mated to a beef bull were 1.35, whereas those of a tenth-parity cow were 2.11. The odds of a female in the worst 10% for total genetic merit being mated to a beef bull were 2.90 times those of a female in the top 10%. Although dystocia was associated with the likelihood of being mated to a beef bull, the actual likelihood did not vary much by level of dystocia experienced. Relative to the first service, the odds of the third and fifth services being to a beef bull were 2.23 and 3.71, respectively. These probability estimates can form the back-end system supporting decisions on mating type for a female within a sire mating advice system but also in risk analysis of farm management.

Short communication: The beef merit of the sire mated to a dairy female affects her subsequent performance

Much of the research to date on dairy × beef matings has focused only on the greater revenue attainable from these beef-cross calves. The objective of the present study was to quantify the mean effect on cow performance following the birth of calves differing in beef merit; all calves were born without calving assistance. Beef merit in the present study was based on the breed of the sire but also its genetic merit for carcass weight and conformation. The cross-sectional study used up to 346,765 calving events from 230,255 Holstein-Friesian cows in 3,604 herds. Performance traits of interest were those associated with milk production, including somatic cell count, as well as female reproductive performance. Sire breed was associated with all yield traits, somatic cell count, and both pregnancy rate and the interval from calving to first service; no association existed with either submission rate or number of services. Relative to a Holstein-Friesian sire, the mean 305-d milk yield (in kg) was 45.22 (standard error, SE = 4.0), 62.0 (SE = 36.8), 65.4 (SE = 9.6), 101.1 (SE = 31.6), 36.7 (SE = 4.9), 51.5 (SE = 10.7), 53.3 (SE = 31.5), and 43.3 (SE = 23.4) less for cows that gave birth to Angus-, Aubrac-, Beligan Blue-, Charolais-, Hereford-, Limousin-, Saler-, or Simmental-sired calves, respectively. Service sire accounted for only 1% of the phenotypic variation in all 3 milk production traits when fitted as a random effect in the model. The regression coefficients of phenotypic milk, fat, and protein yields on sire (of calf) predicted transmitting ability for carcass weight were -1.84 (SE = 0.17), -0.10 (SE = 0.01), and -0.08 kg (SE = 0.01), respectively. The respective regression coefficients on sire (of calf) predicted transmitting ability for carcass conformation (scale of 1 to 15; 1 = poor and 15 = excellent) were -23.46 (SE = 1.81), -1.20 (SE = 0.08), and -1.05 units (SE = 0.06). The biological significance of the sire breed effects or the measure of sire genetic merit on the reproductive traits was either not different from zero or biologically small. Although statistically significant associations existed between sire beef merit and both milk and reproductive performance of the mate, the actual size of the associations was biologically small.

Heteropaternal superfecundation frequently occurs in multiple-bearing mob-mated sheep

Heteropaternal superfecundation may be defined as the fertilisation of two or more ova during the same oestrus cycle as a result of more than one coital act from different males; this results in foetuses being born in the same litter of the same age but different paternity. Heteropaternal superfecundation is more likely to occur in poly-ovulatory species like sheep; moreover, female sheep are often mob-mated with several rams concurrently, thus providing an opportunity for a given female to be served by multiple males during the same oestrus cycle. The objective of the present study was to determine the frequency of heteropaternal superfecundation in six sheep flocks where most of the ewes, lambs and rams were genotyped. A total of 685 multiple-birth litters were available where the sire, dam and all lambs were genotyped. Of the 539 pairs of twins included in the analysis, 160 (i.e. 30%) were sired by two different rams. Of the 137 sets of triplets included in the analysis, 73 (i.e. 53%) were sired by more than one ram. Of the nine sets of quadruplets, eight were sired by two rams with the remaining litter being mono-paternal. The overall incidence of heteropaternal superfecundation among litters was therefore 35%. Given that the incidence of multiple births in these flocks was 65%, heteropaternal superfecundation is expected to be relatively common in sheep; this is especially true as all but two of the litter-mates were polyzygotic. Genotyping of progeny is one practical solution to identity such individuals.

On-farm net benefit of genotyping candidate female replacement cattle and sheep

The net benefit from investing in any technology is a function of the cost of implementation and the expected return in revenue. The objective of the present study was to quantify, using deterministic equations, the net monetary benefit from investing in genotyping of commercial females. Three case studies were presented reflecting dairy cows, beef cows and ewes based on Irish population parameters; sensitivity analyses were also performed. Parameters considered in the sensitivity analyses included the accuracy of genomic evaluations, replacement rate, proportion of female selection candidates retained as replacements, the cost of genotyping, the sire parentage error rate and the age of the female when it first gave birth. Results were presented as an annualised monetary net benefit over the lifetime of an individual, after discounting for the timing of expressions. In the base scenarios, the net benefit was greatest for dairy, followed by beef and then sheep. The net benefit improved as the reliability of the genomic evaluations improved and, in fact, a negative net benefit of genotyping was less frequent when the reliability of the genomic evaluations was high. The impact of a 10% point increase in genomic reliability was, however, greatest in sheep, followed by beef and then dairy. The net benefit of genotyping female selection candidates reduced as replacement rate increased. As genotyping costs increased, the net benefit reduced irrespective of the percentage of selection candidates kept, the replacement rate or even the population considered. Nonetheless, the association between the genotyping cost and the net benefit of genotyping differed by the percentage of selection candidates kept. Across all replacement rates evaluated, retaining 25% of the selection candidates resulted in the greatest net benefit when genotyping cost was low but the lowest net benefit when genotyping cost was high. Genotyping breakeven cost was non-linearly associated with the percentage of selection candidates retained, reaching a maximum when 50% of selection candidates were retained, irrespective of replacement rate, genomic reliability or the population. The genotyping breakeven cost was also non-linearly associated with replacement rate. The approaches outlined within provide the back-end framework for a decision support tool to quantify the net benefit of genotyping, once parameterised by the relevant population metrics.

Observed progeny performance validates the benefit of mating genetically elite beef sires to dairy females

While several studies in cattle have confirmed the improved performance achievable from selection on total merit indexes, these studies have solely been confined to specific-purpose beef or dairy total merit indexes. Validation studies of total merit indexes used to select beef sires for use on dairy females are lacking. The objective here was to fill this void by quantifying the performance of beef × dairy progeny where the sire excels in either a total merit index encompassing calving performance and beef performance traits (dairy-beef index; DBI) or excels in a subindex based solely on calving performance (CLV); for comparative purposes, these beef × dairy progeny were also compared with dairy × dairy progeny. A total of 123,785 calving records from 101,773 dairy cows calving in 3,065 dairy herds were used; of these, 48,875 progeny also had carcass information. The beef sires were stratified into 5 equally sized groups based separately on their DBI or CLV. Linear and threshold mixed models were used to compare calving and carcass performance of all 3 sire genotypes. Of the 415 sires that ranked in the highest of the 5 strata on the CLV subindex, only 52% of them ranked in the highest stratum for the DBI. The percentage of primiparae requiring any assistance at calving was 2 to 3 percentage units greater for the higher DBI sires relative to both the higher CLV beef sires and the dairy sires (not ranked on anything); no difference existed in multiparae. The extent of calving difficulty in primiparae was, however, less in higher DBI beef sires relative to both the higher CLV beef sires and the dairy sires, although the differences were biologically small. Perinatal mortality was greatest in the beef sires relative to the dairy sires, but no difference existed between the high CLV or high DBI beef sires. No difference in progeny gestation length was evident between the high DBI or high CLV beef sires, although both were >2 d longer than progeny from dairy sires. The higher DBI sires produced progeny with heavier, more conformed carcasses relative to the progeny from both high CLV beef sires and dairy sires. No differences existed between the progeny of the beef sires ranked highly on the CLV versus those ranked highly on the DBI for the probability of achieving the specification for carcass weight (between 270 and 380 kg) or fat score; the higher DBI animals, however, had a 4 to 10% greater probability of achieving the minimum carcass conformation required. In all instances, the beef sires had a greater probability of achieving all specifications relative to the progeny from the dairy sires with the difference for conformation being particularly large. Results indicate that more balanced progeny can be generated using a DBI, helping meet the requirements of both dairy and beef producers. Ignoring market failure across sectors, using higher DBI sires could increase dairy herd profit by 3 to 5% over and above the status quo approach to selection in dairy (i.e., CLV subindex).

Choice of artificial insemination beef bulls used to mate with female dairy cattle

Understanding the preferences of dairy cattle producers when selecting beef bulls for mating can help inform beef breeding programs as well as provide default parameters in mating advice systems. The objective of the present study was to characterize the genetic merit of beef artificial insemination (AI) bulls used in dairy herds, with particular reference to traits associated with both calving performance and carcass merit. The characteristics of the beef AI bulls used were compared with those of the dairy AI bulls used on the same farms. A total of 2,733,524 AI records from 928,437 females in 5,967 Irish dairy herds were used. Sire predicted transmitting ability (PTA) values and associated reliability values for calving performance and carcass traits based on national genetic evaluations from prior to the insemination were used. Fixed effects models were used to relate both genetic merit and the associated reliability of the dairy and beef bulls used on the farm with herd size, the extent of Holstein-Friesian × Jersey crossbreeding adopted by the herd, whether the herd used a technician insemination service or do-it-yourself, and the parity of the female mated. The mean direct calving difficulty PTA of the beef bulls used was 1.85 units higher than that of the dairy bulls but with over 3 times greater variability in the beef bulls. This 1.85 units equates biologically to an expectation of 1.85 more dystocia events per 100 dairy cows mated in the beef × dairy matings. The mean calving difficulty PTA of the dairy AI bulls used reduced with increasing herd size, whereas the mean calving difficulty PTA of the beef AI bulls used increased as herd size increased from 75 cows or fewer to 155 cows; the largest herds (>155 cows) used notably easier-calving beef bulls, albeit the calving difficulty PTA of the beef bulls was 3.33 units versus 1.67 units for the dairy bulls used in these herds. Although we found a general tendency for larger herds to use dairy AI bulls with lower reliability, this trend was not obvious in the beef AI bulls used. Irrespective of whether dairy or beef AI bulls were considered, herds that operated more extensive Holstein-Friesian × Jersey crossbreeding (i.e., more than 50% crossbred cows) used, on average, easier calving, shorter gestation-length bulls with lighter expected progeny carcasses of poorer conformation. Mean calving difficulty PTA of dairy bulls used increased from 1.39 in heifers to 1.79 in first-parity cows and to 1.82 in second-parity cows, remaining relatively constant thereafter. In contrast, the mean calving difficulty PTA of the beef bulls used increased consistently with cow parity. Results from the present study demonstrate a clear difference in the mean acceptable genetic merit of beef AI bulls relative to dairy AI bulls but also indicates that these acceptable limits vary by herd characteristics.

Reaffirmation of known major genes and the identification of novel candidate genes associated with carcass-related metrics based on whole genome sequence within a large multi-breed cattle population

Background

The high narrow sense heritability of carcass traits suggests that the underlying additive genetic potential of an individual should be strongly correlated with both animal carcass quality and quantity, and therefore, by extension, carcass value. Therefore, the objective of the present study was to detect genomic regions associated with three carcass traits, namely carcass weight, conformation and fat cover, using imputed whole genome sequence in 28,470 dairy and beef sires from six breeds with a total of 2,199,926 phenotyped progeny.

Results

Major genes previously associated with carcass performance were identified, as well as several putative novel candidate genes that likely operate both within and across breeds. The role of MSTN in carcass performance was re-affirmed with the segregating Q204X mutation explaining 1.21, 1.11 and 5.95% of the genetic variance in carcass weight, fat and conformation, respectively in the Charolais population. In addition, a genomic region on BTA6 encompassing the NCAPG/LCORL locus, which is a known candidate locus associated with body size, was associated with carcass weight in Angus, Charolais and Limousin. Novel candidate genes identified included ZFAT in Angus, and SLC40A1 and the olfactory gene cluster on BTA15 in Charolais. Although the majority of associations were breed specific, associations that operated across breeds included SORCS1 on BTA26, MCTP2 on BTA21 and ARL15 on BTA20; these are of particular interest due to their potential informativeness in across-breed genomic evaluations. Genomic regions affecting all three carcass traits were identified in each of the breeds, although these were mainly concentrated on BTA2 and BTA6, surrounding MSTN and NCAPG/LCORL, respectively. This suggests that although major genes may be associated with all three carcass traits, the majority of genes containing significant variants (unadjusted p-value < 10^− 4) may be trait specific associations of small effect.

Conclusions

Although plausible novel candidate genes were identified, the proportion of variance explained by these candidates was minimal thus reaffirming that while carcass performance may be affected by major genes in the form of MSTN and NCAPG/LCORL, the majority of variance is attributed to the additive (and possibly multiplicative) effect of many polymorphisms of small effect.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6071-9) contains supplementary material, which is available to authorized users.

A breeding index to rank beef bulls for use on dairy females to maximize profit

The desire to increase profit on dairy farms necessitates consideration of the revenue attainable from the sale of surplus calves for meat production. However, the generation of calves that are expected to excel in efficiency of growth and carcass merit must not be achieved to the detriment of the dairy female and her ability to calve and re-establish pregnancy early postcalving without any compromise in milk production. Given the relatively high heritability of many traits associated with calving performance and carcass merit, and the tendency for many of these traits to be moderately to strongly antagonistic, a breeding index that encompasses both calving performance and meat production could be a useful tool to fill the void in supporting decisions on bull selection. The objective of the present study was to derive a dairy-beef index (DBI) framework to rank beef bulls for use on dairy females with the aim of striking a balance between the efficiency of valuable meat growth in the calf and the subsequent performance of the dam. Traits considered for inclusion in this DBI were (1) direct calving difficulty; (2) direct gestation length; (3) calf mortality; (4) feed intake; (5) carcass merit reflected by carcass weight, conformation, and fat and the ability to achieve minimum standards for each; (6) docility; and (7) whether the calf was polled. Each trait was weighted by its respective economic weight, most of which were derived from the analyses of available phenotypic data, supplemented with some assumptions on costs and prices. The genetic merit for a range of performance metrics of 3,835 artificial insemination beef bulls from 14 breeds ranked on this proposed DBI was compared with an index comprising only direct calving difficulty and gestation length (the 2 generally most important characteristics of dairy farmers when selecting beef bulls). Within the Angus breed (i.e., the beef breed most commonly used on dairy females), the correlation between the DBI and the index of genetic merit for direct calving difficulty plus gestation length was 0.74; the mean of the within-breed correlations across all other breeds was 0.87. The ranking of breeds changed considerably when ranked based on the top 20 artificial insemination bulls excelling in the DBI versus excelling in the index of calving difficulty and gestation length. Dairy breeds ranked highest on the index of calving difficulty and gestation length, whereas the Holstein and Friesian breeds were intermediate on the DBI; the Jersey breed was one of the poorest breeds on DBI, superior only to the Charolais breed. The results clearly demonstrate that superior carcass and growth performance can be achieved with the appropriate selection of beef bulls for use on dairy females with only a very modest increase in collateral effect on cow performance (i.e., 2-3% greater dystocia expected and a 6-d-longer gestation length).

A mating advice system in dairy cattle incorporating genomic information

This study investigated the effects of alternative mating programs that incorporate genomic information on expected progeny herd performance and inbreeding, as well as methods to include un-genotyped animals in such mating programs. A total of 54,535 Holstein-Friesian cattle with imputed high-density genotypes (547,650 SNP after edits) were available. First, to quantify the accuracy of imputing un-genotyped animals (often an issue in populations), a sub-population of 729 genotyped animals had their genotypes masked, and their allele dosages were imputed, using linear regression exploiting information on genotyped relatives. The reference population for imputation included all genotyped animals, excluding the 729 selected animals and their sires, dams, and grandsires, and had either (1) their sires' genotypes, (2) their dams' genotypes (3) both their sires' and their dams' genotypes, or (4) both their sires' and maternal grandsires' genotypes introduced into the reference population. The correlations between true genotypes and the imputed allele dosages ranged from 0.58 (sire only) to 0.68 (both sire and dam). A herd of 100 cows was then simulated (1,000 replicates) from the sub-population of 729 imputed animals. The top 10 bulls from the genotyped population, based on their total genetic merit index (TMI) were selected to be used as sires. Three mating allotment methods were investigated: (1) random mating, (2) sequential mating based on maximizing only the expected TMI of the progeny, and (3) linear programming to maximize a generated index constructed to maximize genetic merit and minimize expected progeny inbreeding as well as intra- and inter-progeny variability in genetic merit. Relationships among candidate parents were calculated using either the pedigree relationship matrix or the genomic relationship matrix; the latter was constructed using either the true genotypes of both parents or the true genotypes of the sire plus the imputed allele dosages of the dam. Using the genomic co-ancestry estimates resulted in lower average herd expected genomic inbreeding levels compared with using the pedigree-based co-ancestry estimates. Additionally, if the dams were not genotyped, using their imputed allele dosages also resulted in lower average herd expected inbreeding levels compared with using the pedigree co-ancestry estimates. The inter-progeny coefficient of variation for selected traits, milk and fertility, estimated breeding values were reduced by 12 to 65% using the linear programing method compared with sequential mating.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

Risk factors associated with animal mortality in pasture-based, seasonal-calving dairy and beef herds

Animal mortality is indicative of animal health and welfare standards, which are of growing concern to the agricultural industry. The objective of the present study was to ascertain risk factors associated with mortality at multiple life stages in pasture-based, seasonal-calving dairy and beef herds. Males and females were stratified into seven life stages based on age (0 to 2 d, 3 to 7 d, 8 to 30 d, 31 to 182 d, 183 to 365 d, 366 to 730 d, and 731 to 1,095 d) whereas females with ≥1 calving event were further stratified into five life stages based on cow parity number (1, 2, 3, 4, and 5). Mortality was defined as whether an animal died during each life stage; only animals that either survived the entire duration or died during a life stage were considered. The data, following edits, consisted of 4,404,122 records from 1,358,712 animals. Multivariable logistic regression was used to estimate the logit of the probability of mortality in each life stage separately. The odds of a young animal (i.e., aged ≤ 1,095 d) dying was generally greater if veterinary assistance was required at their birth relative to no assistance (odds ratio [OR]: 3.10 to 31.85), if the animal was a twin relative to a singleton (OR: 1.46 to 2.31) or if the animal was male relative to female (OR: 1.14 to 6.15). Moreover, the odds of a cow (i.e., females with ≥1 calving event) dying were greater when she required veterinary assistance at calving (OR: 2.69 to 7.55) compared with a cow that did not require any assistance, if she produced twin relative to singleton progeny (OR: 1.59 to 2.03) or male relative to female progeny (OR: 1.09 to 1.20). Additionally, the odds of a first or second parity cow dying when she herself had received veterinary assistance at birth were only 0.63 to 0.66 times that of a cow that was provided no assistance at birth. For both young animals and cows, the odds of dying generally increased with herd size, whereas animals residing in expanding herds had lower odds of dying. Results from the present study indicate that the risk factors associated with mortality in pasture-based, seasonal-calving herds are similar to those reported in literature in confinement, nonseasonal-calving herds. Moreover, the present study identifies that these risk factors are similar in both dairy and beef herds, yet the magnitude of the association often differs and also changes with life stage.

Genetic variability in the humoral immune response to bovine herpesvirus-1 infection in dairy cattle and genetic correlations with performance traits

Bovine herpesvirus-1 (BoHV-1) is a viral pathogen of global significance that is known to instigate several diseases in cattle, the most notable of which include infectious bovine rhinotracheitis and bovine respiratory disease. The genetic variability in the humoral immune response to BoHV-1 has, to our knowledge, not ever been quantified. Therefore, the objectives of the present study were to estimate the genetic parameters for the humoral immune response to BoHV-1 in Irish female dairy cattle, as well as to investigate the genetic relationship between the humoral immune response to BoHV-1 with milk production performance, fertility performance, and animal mortality. Information on antibody response to BoHV-1 was available to the present study from 2 BoHV-1 sero-prevalence research studies conducted between the years 2010 to 2015, inclusive; after edits, BoHV-1 antibody test results were available on a total of 7,501 female cattle from 58 dairy herds. National records of milk production (i.e., 305-d milk yield, fat yield, protein yield, and somatic cell score; n = 1,211,905 milk-recorded cows), fertility performance (i.e., calving performance, pregnancy diagnosis, and insemination data; n = 2,365,657 cows) together with animal mortality data (i.e., birth, farm movement, death, slaughter, and export events; n = 12,853,257 animals) were also available. Animal linear mixed models were used to quantify variance components for BoHV-1 as well as to estimate genetic correlations among traits. The estimated genetic parameters for the humoral immune response to BoHV-1 in the present study (i.e., heritability range: 0.09 to 0.16) were similar to estimates previously reported for clinical signs of bovine respiratory disease in dairy and beef cattle (i.e., heritability range: 0.05 to 0.11). Results from the present study suggest that breeding for resistance to BoHV-1 infection could reduce the incidence of respiratory disease in cattle while having little or no effect on genetic selection for milk yield or milk constituents (i.e., genetic correlations ranged from -0.13 to 0.17). Moreover, even though standard errors were large, results also suggest that breeding for resistance to BoHV-1 infection may indirectly improve fertility performance while also reducing the incidence of mortality in older animals (i.e., animals >182 d of age). Results can be used to inform breeding programs of potential genetic gains achievable for resistance to BoHV-1 infection in cattle.

Heritability estimates of meat sensory characteristics are a function of the number of panellists and their inter-correlations

The objective of the present study was to quantify, using simulated data, the impact on estimated heritability of varying the number of panellists and their inter-correlations using meat sensory tenderness in cattle as an example. Estimated parameters from actual sensory-based tenderness scores from 9 individual panellists on 1252 beef cattle were used to parameterise the simulation. A single "tenderness score" for each of 10 panellists was simulated for 15,000 cattle. Heritability estimates were calculated for each of the 10 panellists individually as well as the mean score per animal for all n combinations of panellists. Heritability estimates improved with increasing number of panellists in line with expectations from a deterministic equation. The increase in heritability was due to a reduction in the residual variance, albeit the rate of reduction in residual variance declined with each additional panellist included in the calculated mean tenderness score. Results highlight the importance of reporting the number of panellist scores per animal as well as their inter-correlations in sensory-based studies.

The effect of non-steroidal anti-inflammatory drugs on severity of acute pancreatitis and pancreatic necrosis

Introduction Acute pancreatitis (AP) is a common emergency presentation and can be disabling. There is significant morbidity and mortality associated with AP, and it places a considerable burden on the healthcare system. Non-steroidal anti-inflammatory drugs (NSAIDs) have been shown to have a protective effect in some elective contexts. This retrospective study aimed to evaluate the effect of NSAIDs on the course of AP and the severity of the disease. Methods A retrospective analysis was carried out of 324 patients admitted as an emergency with a diagnosis of AP to two UK hospitals. Patients were divided into two groups: those already taking NSAIDs for other co-morbidities and those not taking NSAIDs. Variables compared included: admission to a high dependency or intensive care unit; pancreatic necrosis; pseudocyst development; need for surgery; serum inflammatory markers; modified early warning scores on days 1, 3 and 5; length of stay; and mortality. Results Patients not taking NSAIDs were more likely to have a C-reactive protein level of ≥150mg/l (p=0.007). Patients in the NSAID group experienced less pancreatic necrosis (p=0.019) and lower rates of pseudocyst formation (p=0.010). Other variables showed no difference between the two groups, specifically length of stay and mortality. Conclusions Routine NSAID use may exert a protective effect on the development of AP, its severity, and complications. Therapeutic use of NSAIDs in acute presentations with pancreatitis should be further evaluated.

Symposium review: Breeding a better cow-Will she be adaptable?

Adaption is a process that makes an individual or population more suited to their environment. Long-term adaptation is predicated on ample usable genetic variation. Evolutionary forces influencing the extent and dynamics of genetic variation in a population include random drift, mutation, recombination, selection, and migration; the relative importance of each differs by population (i.e., drift is likely to be more influential in smaller populations) and number of generations exposed to selection (i.e., mutation is expected to contribute substantially to genetic variability following many generations of selection). The infinitesimal model, which underpins most genetic and genomic evaluations, assumes that each quantitative trait is controlled by an infinitely large number of unlinked and non-epistatic loci, each with an infinitely small effect. Under the infinitesimal model, selection is not expected to noticeably alter the allele frequencies, despite a potential substantial change in the population mean; the exception is in the first few generations of selection when genetic variance is expected to decline, after which it stabilizes. Despite the common use of the heritability statistic in quantitative genetics as a descriptor of adaption or response to selection, it is arguably the coefficient of genetic variation that is more informative to gauge adaptation potential and should, therefore, always be cited in such studies; for example, the heritability of fertility traits in dairy cows is generally low, yet the coefficient of genetic variation for most traits is comparable to many other performance traits, thus supporting the observed rapid genetic gain in fertility performance in dairy populations. Empirical evidence from long-term selection studies, across a range of animal and plant species, fails to support the premise that selection will deplete genetic variability. Even after 100 yr (synonymous with 100 generations) of selection in corn for high protein or oil content, there appears to be no obvious plateauing in the response to selection. Although populations in several selection experiments did reach a selection limit after multiple generations of directional selection, this does not equate to an exhaustion of genetic variance; such a declaration is supported by the observed rapid responses to reverse selection once implemented in long-term selection studies. New technologies such as genome-wide enabled selection and genome editing, as well as having the potential to accelerate genetic gain, could also increase the genetic variation, or at least reduce the erosion of genetic variance over time. In conclusion, there is no evidence, either theoretical or empirical, to indicate that dairy cow breeding programs will be unable to adapt to evolving challenges and opportunities, at least not because of an absence of ample genetic variability.

Genetic co-variance functions for live weight, feed intake, and efficiency measures in growing pigs

The objective of the present study was to estimate genetic co-variance parameters pertaining to live weight, feed intake, and 2 efficiency traits (i.e., residual feed intake and residual daily gain) in a population of pigs over a defined growing phase using Legendre polynomial equations. The data set used consisted of 51,893 live weight records and 903,436 feed intake, residual feed intake (defined as the difference between an animal's actual feed intake and its expected feed intake), and residual daily gain (defined as the difference between an animal's actual growth rate and its expected growth rate) records from 10,201 growing pigs. Genetic co-variance parameters for all traits were estimated using random regression Legendre polynomials. Daily heritability estimates for live weight ranged from 0.25 ± 0.04 (d 73) to 0.50 ± 0.03 (d 122). Low to moderate heritability estimates were evident for feed intake, ranging from 0.07 ± 0.03 (d 66) to 0.25 ± 0.02 (d 170). The estimated heritability for residual feed intake was generally lower than those of both live weight and feed intake and ranged from 0.04 ± 0.01 (d 96) to 0.17 ± 0.02 (d 159). The heritability for feed intake and residual feed intake increased in the early stages of the test period and subsequently sharply declined, coinciding with older ages. Heritability estimates for residual daily gain ranged from 0.26 ± 0.03 (d 188) to 0.42 ± 0.03 (d 101). Genetic correlations within trait were strongest between adjacent ages but weakened as the interval between ages increased; however, the genetic correlations within all traits tended to strengthen between the extremes of the trajectory. Moderate to strong genetic correlations were evident among live weight, feed intake, and the efficiency traits, particularly in the early stage of the trial period (d 66 to 86), but weakened with age. Results from this study could be implemented into the national genetic evaluation for pigs, providing comprehensive information on the profile of growth and efficiency throughout the growing period of the animal's life, thus helping producers identify genetically superior animals.

Impact of alternative definitions of contemporary groups on genetic evaluations of traits recorded at lambing

The objective of this study was to quantify the impact of alternative contemporary group definitions for lambing traits on genetic evaluations in the Irish multibreed sheep population. Three lambing traits were considered for analysis: lambing difficulty, birth weight, and survival. Eight alternative contemporary group definitions were investigated for each lambing trait; all contemporary groups were formed within flock of lambing and included (flock by) week of lambing, week of lambing by litter size (i.e., singles vs. multiples), 2-wk interval (i.e., fortnight) of lambing, fortnight of lambing by litter size, month of lambing, and month of lambing by litter size or were based on an optimized algorithm that creates contemporary groups based on animals from the same flock that are born in close proximity of date. Three alternative scenarios were modeled for each of the lambing traits using the contemporary group definitions: the first scenario (termed Current Scenario) represented the editing criteria currently employed in the Irish national genetic evaluations; the second scenario (No Restriction Scenario) removed any restriction on number of records per contemporary group, and the final scenario (Variation Scenario) included only data from contemporary groups with some variability in the dependent variable. Variance components and EBV for each of the 3 lambing traits were estimated using linear mixed models. The direct heritability estimates ranged from 0.09 ± 0.02 to 0.29 ± 0.02 for lambing difficulty, 0.11 ± 0.01 to 0.24 ± 0.01 for birth weight, and 0.05 ± 0.02 to 0.10 ± 0.02 for lamb survival. Irrespective of lambing trait, greater estimated accuracy of the sire EBV was achieved with the No Restriction Scenario. Results for the ability to predict future lambing characteristics, based on only the direct and maternal EBV, revealed that the area under the receiver operator characteristic curve for the dichotomized lambing assistance phenotype varied from 0.56 to 0.66; a lambing event predicted to be in the worst 10% risk category of a difficult lambing on the basis of genetic merit alone was 5.48 times (95% CI: 3.94 to 7.61; < 0.001) more likely to require assistance at lambing compared to a lambing event in the best 10%. Results show that the use of contemporary groups formed over short time periods, coupled with moderate editing of the data, yielded superior predictions for all lambing traits.

Genetic and nongenetic factors associated with milk color in dairy cows

Milk color is one of the sensory properties that can influence consumer choice of one product over another and it influences the quality of processed dairy products. This study aims to quantify the cow-level genetic and nongenetic factors associated with bovine milk color traits. A total of 136,807 spectra from Irish commercial and research herds (with multiple breeds and crosses) were used. Milk lightness (Lˆ*), red-green index (aˆ*) and yellow-blue index (bˆ*) were predicted for individual milk samples using only the mid-infrared spectrum of the milk sample. Factors associated with milk color were breed, stage of lactation, parity, milking-time, udder health status, pasture grazing, and seasonal calving. (Co)variance components for Lˆ*,aˆ*, and bˆ* were estimated using random regressions on the additive genetic and within-lactation permanent environmental effects. Greater bˆ* value (i.e., more yellow color) was evident in milk from Jersey cows. Milk Lˆ* increased consistently with stage of lactation, whereas aˆ* increased until mid lactation to subsequently plateau. Milk bˆ* deteriorated until 31 to 60 DIM, but then improved thereafter until the end of lactation. Relative to multiparous cows, milk yielded by primiparae was, on average, lighter (i.e., greater Lˆ*), more red (i.e., greater aˆ*), and less yellow (i.e., lower bˆ*). Milk from the morning milk session had lower Lˆ*,aˆ*, and bˆ* Heritability estimates (±SE) for milk color varied between 0.15 ± 0.02 (30 DIM) and 0.46 ± 0.02 (210 DIM) for Lˆ*, between 0.09 ± 0.01 (30 DIM) and 0.15 ± 0.02 (305 DIM) for aˆ*, and between 0.18 ± 0.02 (21 DIM) and 0.56 ± 0.03 (305 DIM) for bˆ* For all the 3 milk color features, the within-trait genetic correlations approached unity as the time intervals compared shortened and were generally <0.40 between the peripheries of the lactation. Strong positive genetic correlations existed between bˆ* value and milk fat concentration, ranging from 0.82 ± 0.19 at 5 DIM to 0.96 ± 0.01 at 305 DIM and confirming the observed phenotypic correlation (0.64, SE = 0.01). Results of the present study suggest that breeding strategies for the enhancement of milk color traits could be implemented for dairy cattle populations. Such strategies, coupled with the knowledge of milk color traits variation due to nongenetic factors, may represent a tool for the dairy processors to reduce, if not eliminate, the use of artificial pigments during milk manufacturing.