Method to establish average relationships among Holstein bull populations over time

Average relationship coefficients between groups of animals were calculated by repeating three calculation steps until convergence was reached. The first step was a stratified sampling of pedigrees. Second, relationships between animals in the sample were calculated, using the group relationships computed in a previous round for unknown parent relationships. Third, results were accumulated by group. Average relationship coefficients calculated in this way can be used for such purposes as to better understand the structure of a population, to more accurately calculate inbreeding coefficients, and to assign unknown parent groups using clustering methods. This method was applied to the pedigree data for the worldwide Holstein sire population. The data was derived from Interbull and North American pedigree databases. Groups were defined by year of birth, sex, and country. Average relationships between US and Canadian bulls increased across time from 0.02 in 1960 to 0.12 in 1999. The increase was continuous except for a plateau in 1992, followed by a slight decrease until 1995 when the trend to increase resumed. Average relationships across time comparing Canadian with European bulls and US with European bulls are quite similar. They increased from about 0 in 1963 to 0.10 in 1999. After a peak of 0.05 in 1966, the relationships dropped to 0.03 in 1976 before starting to increase again. For bulls born in 1996 in the US, Canada, and Germany (representing a typical European country), the average relationships were: 0.116 (CanadaCanada), 0.112 (Canada-USA), 0.1 (Canada-Germany), 0.101 (USA-USA), 0.101 (USAGermany) and 0.084 (Germany-Germany). These results demonstrate the recent dramatic increase in relationships and inbreeding in the worldwide Holstein population, and show the power of this method of calculating relationship coefficients

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, like 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 world-wide to estimate cow body reserves and 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. Body condition score change records were merged with milk MIR spectra recorded on the same week. The data set 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. ΔBCS data in both the first 120 DIM 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 data set and of 3.63*10-3 BCS units in the 305-d data set. 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 data set). 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 data set and of the prediction model used, the combining DIM information with MIR spectral data as prediction variables reduced the RMSE compared with inclusion of DIM alone, albeit the benefit was small (the RMSE from cross-validation was reduced 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 was 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

Review: Opportunities and challenges for the genetic selection of dairy calf disease traits.

Interest in dairy cow health continues to grow as we better understand health's relationship with production potential and animal welfare. Over the past decade, efforts have been made to incorporate health traits into national genetic evaluations. However, they have focused on the mature cow, with calf health largely being neglected. Diarrhoea and respiratory disease comprise the main illnesses with regard to calf health. Conventional methods to control calf disease involve early separation of calves from the dam and housing calves individually. However, public concern regarding these methods, and growing evidence that these methods may negatively impact calf development, mean the dairy industry may move away from these practices. Genetic selection may be a promising tool to address these major disease issues. In this review, we examined current literature for enhancing calf health through genetics and discussed alternative approaches to improve calf health via the use of epidemiological modelling approaches, and the potential of indirectly selecting for improved calf health through improving colostrum quality. Heritability estimates on the observed scale for diarrhoea ranged from 0.03 to 0.20, while for respiratory disease, estimates ranged from 0.02 to 0.24. The breadth in these ranges is due, at least in part, to differences in disease prevalence, population structure, data editing and models, as well as data collection practices, which should be all considered when comparing literature values. Incorporation of epidemiological theory into quantitative genetics provides an opportunity to better determine the level of genetic variation in disease traits, as it accounts for disease transmission among contemporaries. Colostrum intake is a major determinant of whether a calf develops either respiratory disease or diarrhoea. Colostrum traits have the advantage of being measured and reported on a continuous scale, which removes the issues classically associated with binary disease traits. Overall, genetic selection for improved calf health is feasible. However, to ensure the maximum response, first steps by any industry members should focus efforts on standardising recording practices and encouragement of uploading information to genetic evaluation centres through herd management software, as high-quality phenotypes are the backbone of any successful breeding programme

Investigating the potential for genetic selection of dairy calf disease traits using management data.

Genetic selection could be a tool to help improve the health and welfare of calves, however, to date, there is limited research on the genetics of calfhood diseases. This study aimed to understand the current impact of calf diseases, by investigating incidence rates, estimating genetic parameters, and providing industry recommendations to improve calf disease recording practices on farms. Available calf disease data comprised of 69,695 Holstein calf disease records for respiratory problems (RESP) and diarrhea (DIAR), from 62,361 calves collected on 1,617 Canadian dairy herds from 2006 to 2021. Single and multiple trait analysis using both a threshold and linear animal model for each trait were evaluated. Furthermore, each trait was analyzed using 2 scenarios with respect to minimum disease incidence threshold criterion (herd-year incidence of at least 1% and 5%) to highlight the impact of different filtering thresholds on selection potential. Observed scale heritability estimates for RESP and DIAR ranged from 0.02 to 0.07 across analyses, while estimated genetic correlations between the traits ranged from 0.50 to 0.62. Sires were compared based on their estimated breeding value and their diseased daughter incidence rates. On average, calves born to the bottom 10% of sires were 1.8 times more likely to develop RESP and 1.9 times to develop DIAR compared with daughters born to the top 10% of sires. Results from the current study are promising for the inclusion of both DIAR and RESP in Canadian genetic evaluations. However, for effective genetic evaluation we require standardized approaches on data collection and industry outreach to highlight the importance of collecting and uploading this information to herd management software. In particular, it is important that the herd management software is accessible to the national milk recording system to allow for use in national genetic evaluation

Phenotypic analysis of heat stress in Holsteins using test-day production records and NASA POWER meteorological data.

Weather station data and test-day production records can be combined to quantify the effects of heat stress on production traits in dairy cattle. However, meteorological data sets that are retrieved from ground-based weather stations can be limited by spatial and temporal data gaps. The National Aeronautics and Space Administration Prediction of Worldwide Energy Resources (NASA POWER) database provides meteorological data over regions where surface measurements are sparse or nonexistent. The first aim of this study was to determine whether NASA POWER data are a viable alternative resource of weather data for studying heat stress in Canadian Holsteins. The results showed that average, minima, and maxima ambient temperature and dewpoint temperature as well as 4 different types of temperature-humidity index (THI) values from NASA POWER were highly correlated to the corresponding values from weather stations (regression R2 > 0.80). However, the NASA POWER values for the daily average, minima, and maxima wind speed and relative humidity were poorly correlated to the corresponding weather station values (regression R2 = 0.10 to 0.49). The second aim of this study was to quantify the influence of heat stress on Canadian dairy cattle. This was achieved by determining the THI values at which milk, protein, and fat yield started to decline due to heat stress as well as the rates of decline in these traits after the respective thresholds, using segmented polynomial regression models. This was completed for both primiparous and multiparous cows from 5 regions in Canada (Ontario, Quebec, British Columbia, the Prairies, and the Atlantic Maritime). The results showed that all production traits were negatively affected by heat stress and that the patterns of responses for milk, fat, and protein yields to increasing THI differed from each other. We found 3 THI thresholds for milk yield, 1 for fat yield, and 2 for protein yield. All thresholds marked a change in rate of decrease in production yield per unit THI, except for the first milk yield threshold, which marked a greater rate of increase. The first thresholds for milk yield ranged between 47 and 50, the second thresholds ranged between 61 and 69, and the third thresholds ranged between 72 and 76 THI units. The single THI threshold for fat yield ranged between 48 and 55 THI units. Finally, the first and second thresholds ranged between 58 and 62 THI units and 72 and 73 THI units for protein yield, respectively

Associations between feed efficiency and aspects of lactation curves in primiparous Holstein dairy cattle.

Genetic selection for improved feed efficiency in dairy cattle has received renewed attention over the last decade to address the needs of a growing global population. As milk yield is a critical component of feed efficiency metrics in dairy animals, our objective was to evaluate the associations between feed efficiency in primiparous Holstein cattle and parameters of a mathematical model describing individual lactation curves. The Dijkstra lactation curve model was fit to individual lactation records from 34 Holstein heifers with previously estimated measures of feed efficiency. We found that the optimal fit of the Dijkstra model was achieved using daily milk yield records up to 21 d in milk to capture the rise to peak milk yield and using monthly dairy herd improvement records for the remainder of lactation to accurately characterize lactation persistency. In the period of lactation before peak milk yield, improved feed efficiency was associated with a faster increase in daily milk yield over a shorter period of time at the expense of increased mobilization of body reserves; this serves to reinforce the concept that dairy cattle are primarily capital breeders versus income breeders. Feed efficiency in the period following peak lactation, as measured by gross feed efficiency, return over feed costs, and net energy efficiency of lactation, was positively associated with higher peak milk yield. The findings in early lactation suggest that estimates of feed efficiency could be improved by evaluating feed efficiency relative to conception, rather than parturition and lactation, to better account for the energy stored and released from body reserves in capital breeding

Estimating the environmental impact of dairy cattle breeding programs through emission intensity.

A recently developed methodological approach for determining the greenhouse gas emissions impact of national breeding programs was applied to measure the effects of current and future breeding goals on the emission intensity (EI) of the Canadian dairy industry. Emission intensity is the ratio of greenhouse gas outputted in comparison to the product generated. Traits under investigation affected EI by either decreasing the direct emissions yield (i.e. increasing feed performance), changing herd structure (i.e. prolonging herd life) or through the dilution effect of increased production (i.e. increasing fat yield). The intensity value (IV) of each trait, defined as the change in emissions' intensity per unit change in each trait, was calculated for each of the investigated traits. The IV trend of these traits was compared for the current and prospective selection index, as well as for a system with and without quota (the supply management policy designed to prevent overproduction). The overall EI of the average genetic merit Canadian dairy herd per breeding female was 5.07 kg CO2eq/kg protein equivalent output. The annual reduction in EI due to the improvement of production traits was -0.027, -0.018 and -0.006 for fat, protein and milk other solids, respectively. The functional traits, herd life and mastitis resistance, had more modest effects (-0.008 and -0.001, respectively). These results are consistent with international studies that identified traits related to production, survival, health and fertility as having the largest impact on the environmental footprint of dairy cattle. Overall, the dairy industry is becoming more efficient by reducing its EI through selection of environmentally favorable traits, with a 1% annual reduction of EI in Canada

Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins

Background: The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. Methods: In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. Results: In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a - 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. Conclusions: ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations