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 composed 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 effect 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, standardized approaches on data collection and industry outreach to highlight the importance of collecting and uploading this information to herd management software is required. 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.

Estimates of genetic parameters for rumination time, feed efficiency, and methane production traits in first-lactation Holstein cows.

The large-scale recording of traits such as feed efficiency (FE) and methane emissions (ME) for use in genetic improvement programs is complex, costly, and time-consuming. Therefore, heritable traits that can be continuously recorded in dairy herds and are correlated with FE and ME traits could provide useful information for genetic evaluation. Rumination time has been suggested to be associated with FE, methane production (MeP; ME in g/d), and production traits at the phenotypic level. Therefore, the objective of this study was to investigate the genetic relationships among rumination time (RT), FE, methane and production traits using 7,358 records from 656 first-lactation Holstein cows. The estimated heritabilities were moderate for RT (0.45 ± 0.14), MeP (0.36 ± 0.12), milk yield (0.40 ± 0.08), fat yield (0.29 ± 0.06), protein yield (0.32 ± 0.07), and energy-corrected milk (0.28 ± 0.07), but were low and nonsignificant for FE (0.15 ± 0.07), which was defined as the residual of the multiple linear regression of DMI on energy-corrected milk and metabolic body weight. A favorable negative genetic correlation was estimated between RT and MeP (-0.53 ± 0.24), whereas a positive favorable correlation was estimated between RT and energy-corrected milk (0.49 ± 0.11). The estimated genetic correlation of RT with FE (-0.01 ± 0.17) was not significantly different from zero but showed a trend of a low correlation with dry matter intake (0.21 ± 0.13). These results indicate that RT is genetically associated with MeP and milk production traits, but high standard errors indicate that further analyses should be conducted to verify these findings when more data for RT, MeP, and FE become available.

Estimation of genetic parameters for feed efficiency traits using random regression models in dairy cattle.

Feed efficiency has become an increasingly important research topic in recent years. As feed costs rise and the environmental impacts of agriculture become more apparent, improving the efficiency with which dairy cows convert feed to milk is increasingly important. However, feed intake is expensive to measure accurately on large populations, making the inclusion of this trait in breeding programs difficult. Understanding how the genetic parameters of feed efficiency and traits related to feed efficiency vary throughout the lactation period is valuable to gain understanding into the genetic nature of feed efficiency. This study used 121,226 dry matter intake (DMI) records, 120,500 energy-corrected milk (ECM) records, and 98,975 metabolic body weight (MBW) records, collected on 7,440 first-lactation Holstein cows from 6 countries (Canada, Denmark, Germany, Spain, Switzerland, and the United States), from January 2003 to February 2022. Genetic parameters were estimated using a multiple-trait random regression model with a fourth-order Legendre polynomial for all traits. Weekly phenotypes for DMI were re-parameterized using linear regressions of DMI on ECM and MBW, creating a measure of feed efficiency that was genetically corrected for ECM and MBW, referred to as genomic residual feed intake (gRFI). Heritability (SE) estimates varied from 0.15 (0.03) to 0.29 (0.02) for DMI, 0.24 (0.01) to 0.29 (0.03) for ECM, 0.55 (0.03) to 0.83 (0.05) for MBW, and 0.12 (0.03) to 0.22 (0.06) for gRFI. In general, heritability estimates were lower in the first stage of lactation compared with the later stages of lactation. Additive genetic correlations between weeks of lactation varied, with stronger correlations between weeks of lactation that were close together. The results of this study contribute to a better understanding of the change in genetic parameters across the first lactation, providing insight into potential selection strategies to include feed efficiency in breeding programs.

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.

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.

Genetic evaluation of heat tolerance in Holsteins using test-day production records and NASA POWER weather data.

Heat stress is a prominent issue in livestock production, even for intensively housed dairy herds in Canada. Production records and meteorological data can be combined to assess heat tolerance in dairy cattle. The overall aim of this study was to evaluate the possibility of genetic evaluation for heat tolerance in Canadian dairy cattle. The 2 specific objectives were (1) to estimate the genetic parameters for milk, fat, and protein yield for Holsteins while accounting for high environmental heat loads, and (2) to determine if a genotype-by-environment interaction causes reranking of top-ranked sires between environments with low and high heat loads. A repeatability test-day model with a heat stress function was used to evaluate the genetic merit for milk, fat, and protein yield under heat stress and at thermal comfort for first parity in 5 regions in Canada. The heat stress function for each trait was defined using a specific temperature-humidity index (THI) threshold. The purpose of this function was to quantify the level of heat stress that was experienced by the dairy cattle. The estimated genetic correlation between the general additive genetic effect and the additive effect on the slope of the change in the trait phenotype for milk, fat, and protein yield ranged from -0.16 to -0.30, -0.20 to -0.44, and -0.28 to -0.42, respectively. These negative correlations imply that there is an antagonistic relationship between sensitivity to heat stress and level of production. The heritabilities for milk, fat, and protein yield at 15 units above the THI threshold ranged from 0.15 to 0.27, 0.11 to 0.15, and 0.11 to 0.15, respectively. Finally, the rank correlations between the breeding values from a repeatability model with no heat stress effect and the breeding values accounting for heat stress for the 100 top-ranked bulls indicated possible interaction between milk production traits and THI, resulting in substantial reranking of the top-ranked sires in Canada, especially for milk yield. This is the first study to implement weather data from the NASA POWER database in a genetic evaluation of heat tolerance in dairy cattle. The NASA POWER database is a novel alternative meteorological resource that is potentially more reliable and consistent and with broader coverage than weather station data increasing the number of animals that could be included in a heat stress evaluation.

The effect of progesterone concentrations during superovulation of Holstein heifers in a randomized trial.

The aim of this study was to evaluate the effect of different progesterone (P4) concentrations during the follicular growth on the intensity of estrous expression, ovarian response to the superovulatory treatment, and embryo production and quality in superovulated heifers. A total of 63 Holstein heifers were randomly assigned into 2 experimental groups: Low P4 (n = 31) and High P4 (n = 32). Animals received a pre-synchronization protocol followed by a protocol of superovulation that included the allocated P4 treatment. Activity was monitored continuously by an automated activity monitor, and estrus characteristics (maximum intensity and duration) were recorded. Embryo collection was performed 7 d post artificial insemination (AI). Embryos were counted and graded from good or excellent (1) to degenerated (4). The outcomes of interest were: number and diameter of follicles at the time of AI, ovulation success (confirmed 7 d post-AI), time to estrus event, maximum intensity and duration of estrus, number and quality of embryos. Data were analyzed according to the type of outcome variable using logistic, linear, or Poisson regression models. A total of 105 embryos (High P4: n = 42; Low P4: n = 63) were graded for quality. Different P4 levels did not affect the maximum intensity (High P4 = 497.8 ± 23.9%; Low P4 = 542.2 ± 23.5%) or the duration (High P4 = 13.5 ± 1.5 h; Low P4 = 14.3 ± 1.4 h) of estrus. Heifers in the High P4 treatment had greater number of follicles at time of AI (High P4 = 16.6 ± 1.6 follicles; Low P4 = 13.9 ± 1.2 follicles), but with smaller diameter (High P4 = 11.3 ± 0.1 mm; Low P4 = 12.0 ± 0.1 mm) compared with Low P4. High P4 heifers tended to have better embryo quality compared with Low P4 heifers (odds ratio = 1.98; 95% CI = 0.90-4.35). High P4 heifers had less embryos than Low P4 heifers, but this was modified by the CIDR (intravaginal implant of P4) removal to estrus interval (interval 0-21 h: mean ratio = 1.15, 95% CI = 0.42-1.87; interval 22-46 h: mean ratio = 0.58, 95% CI = 0.27-0.96). Although estrous expression was not associated with embryo quality, as the duration and the maximum intensity of estrous expression increased, the number of embryos recovered 7 d post-AI increased (duration: mean ratio = 1.04; 95% CI = 1.03-1.05; maximum intensity: mean ratio = 1.50; 95% CI = 1.42-1.58). In conclusion, P4 during the follicular growth, and intensity of estrus, are playing a role in regulating the quality and the number of embryos produced by superovulated heifers. This study was supported by contributions from Resilient Dairy Genome Project and the Natural Sciences and Engineering Research Council.

Machine learning classification of breeding protocol descriptions from Canadian Holsteins.

Dairy farmers are motivated to ensure cows become pregnant in an optimal and timely manner. Although timed artificial insemination (TAI) is a successful management tool in dairy cattle, it masks an animal's innate fertility performance, likely reducing the accuracy of genetic evaluations for fertility traits. Therefore, separating fertility traits based on the recorded management technique involved in the breeding process or adding the breeding protocol as an effect to the model can be viable approaches to address the potential bias caused by such management decisions. Nevertheless, there is a lack of specificity and uniformity in the recording of breeding protocol descriptions by dairy farmers. Therefore, this study investigated the use of 8 supervised machine learning algorithms to classify 1,835 unique breeding protocol descriptions from 981 herds into the following 2 classes: TAI or other than TAI. Our results showed that models that used a stacking classifier algorithm had the highest Matthews correlation coefficient (0.94 ± 0.04, mean ± SD) and maximized precision and recall (F1-score = 0.96 ± 0.03) on test data. Nonetheless, their F1-scores on test data were not different from 5 out of the other 7 algorithms considered. Altogether, results presented herein suggest machine learning algorithms can be used to produce robust models that correctly identify TAI protocols from dairy cattle breeding records, thus opening the opportunity for unbiased genetic evaluation of animals based on their natural fertility.

Occurrence and greater intensity of estrus in recipient lactating dairy cows improve pregnancy per embryo transfer.

The aim of this study was to determine the association between occurrence and intensity of estrous expression with pregnancy success in recipient lactating dairy cows subjected to embryo transfer (ET). Two observational studies were conducted. Holstein cows were synchronized using the same timed ET protocol, based on estradiol and progesterone in both experiments. At 9 d after the end of the timed ET protocol only animals that had ovulated were implanted with a 7-d embryo [experiment 1 (Exp. 1); n = 1,401 ET events from 1,045 cows, and experiment 2 (Exp. 2); n = 1,147 ET events from 657 cows]. Embryos were produced in vivo (Exp. 1 and Exp. 2) and in vitro (only Exp. 2), then transferred to recipient cows as fresh or frozen-thawed. Pregnancy was confirmed at 29 and 58 d after the end of timed ET protocol. In Exp. 1, animals had their estrous expression monitored through a tail chalk applied on the tail head of the cows and evaluated daily for chalk removal (no estrus: 100% of chalk remaining; estrus: <50% of chalk remaining). In Exp. 2, cows were continuously monitored by a leg-mounted automated activity monitor. Estrous expression was quantified using the relative increase in physical activity at estrus in relation to the days before estrus. Estrous expression was classified as no estrus [<100% relative increase in activity (RI)], weak intensity (100-299% RI), and strong intensity (≥300% RI). Data were analyzed by analysis of variance using mixed linear regression models (GLIMMIX) in SAS (SAS Institute Inc.). A total of 65.2% (914/1,401) and 89.2% (1,019/1,142) of cows from Exp. 1 and Exp. 2, respectively, displayed estrus at the end of the ovulation synchronization protocol. In Exp. 1, cows expressing estrus before to ET had greater pregnancy per ET than those that did not [41.0 ± 2.3% (381/914) vs. 31.5 ± 2.9% (151/487), respectively]. Similarly, in Exp. 2, cows classified in the strong intensity group had greater pregnancy per ET compared with cows in the weak intensity and no estrus groups [41.3 ± 2.2% (213/571) vs. 32.7 ± 2.7% (115/353) vs. 11.3 ± 3.5% (26/218), respectively]. There was no effect of ET type on pregnancy per ET in Exp. 1. However, in Exp. 2, cows that received an in vivo-produced embryo, either fresh or frozen, had greater pregnancy per ET compared with cows that received in vitro-produced embryo. Cows receiving embryos in the early blastocyst and blastocyst stage had greater fertility compared with cows receiving embryos in the morula stage. There was an interaction between the occurrence of estrus and the stage of embryo development on pregnancy per ET, cows which displayed estrus and received a morula or early blastocyst had greater pregnancy per ET than cows that did not display estrus. In conclusion, the occurrence and the intensity of estrous expression improved pregnancy per ET in recipient lactating dairy cows and thus could be used as a tool to assist in the decision making of reproduction strategies in dairy farms.

Effect of synchronized breeding on genetic evaluations of fertility traits in dairy cattle.

Estrus detection has become more difficult over the years due to decreases in the estrus expression of high-producing dairy cows, and increased herd sizes and animal density. Through the use of hormonal synchronization protocols, also known as timed artificial insemination (TAI) protocols, it is possible to alleviate some of the challenges associated with estrus detection. However, TAI masks cows' fertility performance, resulting in an unfair comparison of treated animals and innately fertile animals. Consequently, genetically inferior and superior cows show similar phenotypes, making it difficult to distinguish between them. As genetic programs rely on the collection of accurate phenotypic data, phenotypes collected on treated animals likely add bias to genetic evaluations. In this study, to assess the effect of TAI, the rank correlation of bulls for a given trait using only TAI records were compared with the same trait using only heat detection records. A total of 270,434 records from 192,539 animals split across heifers, first and second parity cows were analyzed for the traits: calving to first service, first service to conception, and days open. Results showed large reranking across all traits and parities between bulls compared based on either having only TAI records or only heat detection records, suggesting that a bias does indeed exist. Large reranking was also observed for both the heat detection and TAI groups among the top 100 bulls in the control group, which included all records. Furthermore, breeding method was added to the model to assess its effect on bull ranking. However, there were only minor changes in the rank correlations between scenario groups. Therefore, more complex methods to account for the apparent bias created by TAI should be investigated; for this, the method by which these data are collected needs to be improved through creating a standardized way of recording breeding codes. Though the results of this study suggest the presence of bias within current fertility evaluations, additional research is required to confirm the findings of this study, including looking at high-reliability bulls specifically, to determine if the levels of reranking remain. Future studies should also aim to understand the potential genetic differences between the fertility traits split via management technology, possibly in a multiple-trait analysis.

Potential effects of hormonal synchronized breeding on genetic evaluations of fertility traits in dairy cattle: A simulation study.

About 30% of producers use hormone protocols to synchronize ovulation and perform timed artificial insemination (AI) in Canada. Days from calving to first service (CTFS) and first service to conception (FSTC) become masked phenotypes leading to biased genetic evaluations of cows for these fertility traits. The objectives of this study were to (1) demonstrate and quantify the potential amount of bias in genetic evaluations, and (2) find a procedure that could remove the bias. Simulation was used for both objectives. The proposed solution was to identify cows that have been treated by hormone protocols, make their CTFS and FSTC missing, and perform a multiple trait analysis including traits that have high genetic correlations with CTFS and FSTC, and which are not affected by the hormone protocols themselves. A total of 12 scenarios (S1-S12) were tested, changing the percentage of herds and cows that were randomly selected to be under timed AI. Cows that were given hormone protocols had CTFS of 86 d and FSTC of 0, which were used in genetic evaluation. Four criteria were used to indirectly measure the presence of bias: (1) the correlation between true (TBV) and estimated (EBV) breeding values (accuracy); (2) the differences in the mean EBV of top 25, 50, and 75 sires; (3) changes in correlation between TBV and EBV rankings; and (4) the changes in mean EBV over the simulated generations. All criteria changed unfavorably and proportionally to the increased use of timed AI. The accuracy within each class of animals (cows, dams, or sires) decreased proportionally with increased use of timed AI, varying from 0.32 (S12) to 0.52 (S1) for bull EBV for CTFS. The average EBV of the top sires (best 25, 50, 75, or 100 sires) approached population average EBV values when increasing the number of treated animals. The sire rank correlation between EBV and TBV within simulated scenarios was smaller for scenarios with more synchronized animals, going from 0.38 (S12) to 0.67 (S1). The long-term use of hormonal synchronized cows clearly decreased the mean EBV over generations in the population for CTFS and FSTC. The inclusion of genetically correlated traits in a multiple trait model was effective in removing the bias due to the presence of hormonal synchronized cows. However, given the constraints within the simulation, it is important that further investigation with real data is conducted to determine the true effect of including timed AI records within genetic evaluations of fertility traits in dairy cattle.

Estimated genetic parameters for all genetically evaluated traits in Canadian Holsteins.

Genetic improvement is a crucial tool to deal with the increasing demand for high quality, sustainably produced dairy. Breeding programs are based on genetic parameters, such as heritability and genetic correlations, for economically important traits in a population. In this study, we estimated population genetic parameters and genetic trends for 67 traits evaluated on heifers and first-lactation Canadian Holstein cows. The data consisted of approximately 500,000 records with pedigree information collected from 1980 to 2019. Genetic parameters were estimated using bivariate linear animal models under a Bayesian approach. Analyses for the 67 traits resulted in 2,211 bivariate combinations, from which the estimated genetic parameters are reported here. The most highly heritable traits were fat percent (0.66) and protein percent (0.69), followed by stature (0.47). Lowest heritabilities (0.01) were observed for disease-related traits, such as lameness and toe ulcer, and calf survival. The genetic correlations between gestation length, calf size, and calving ease measured on both heifer and cows were close to unity. On the other hand, traits such as body condition score and pin width, cystic ovaries and sole ulcer, rear teat placement, and toe ulcer were genetically unrelated. This study reports genetic parameters that have not been previously published for Canadian Holstein cows, and provides updates of those previously estimated. These estimates are useful for building new indexes, updating existing selection indexes, and for predicting correlated responses due to inclusion of novel traits in the breeding programs.

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.

Discovering lethal alleles across the turkey genome using a transmission ratio distortion approach.

Deviation from Mendelian inheritance expectations (transmission ratio distortion, TRD) has been observed in several species, including the mouse and humans. In this study, TRD was characterized in the turkey genome using both allelic (specific- and unspecific-parent TRD) and genotypic (additive- and dominance-TRD) parameterizations within a Bayesian framework. In this study, we evaluated TRD for 23 243 genotyped Turkeys across 56 393 autosomal SNPs. The analyses included 500 sires, 2013 dams and 11 047 offspring (trios). Three different haplotype sliding windows of 4, 10 and 20 SNPs were used across the autosomal chromosomes. Based on the genotypic parameterizations, 14 haplotypes showed additive and dominance TRD effects highlighting regions with a recessive TRD pattern. In contrast, the allelic model uncovered 12 haplotype alleles with the allelic TRD pattern which showed an underrepresentation of heterozygous offspring in addition to the absence of homozygous animals. For regions with the allelic pattern, only one particular region showed a parent-specific TRD where the penetrance was high via the dam, but low via the sire. The gene set analysis uncovered several gene ontology functional terms, Reactome pathways and several Medical Subject Headings that showed significant enrichment of genes associated with TRD. Many of these gene ontology functional terms (e.g. mitotic spindle assembly checkpoint, DRM complex and Aneuploidy), Reactome pathways (e.g. Mismatch repair) and Medical Subject Headings (e.g. Adenosine monophosphate) are known to be related to fertility, embryo development and lethality. The results of this study revealed potential novel candidate lethal haplotypes, functional terms and pathways that may enhance breeding programs in Turkeys through reducing mortality and improving reproduction rate.

The dynamic behavior of feed efficiency in primiparous dairy cattle.

A variety of measures are currently used to evaluate feed efficiency in livestock, with some of the most popular being the linear regression of nutrient intake on the various energy sinks of the cohort being evaluated. The adoption of novel on-farm technologies has greatly increased the breadth and depth of data being collected on a daily basis, which creates the opportunity for the implementation of more advanced feed efficiency models. The objectives of the current study were to leverage data collected using high-throughput biometric technologies to evaluate changes of early-lactation feed efficiency on a daily basis and to develop a novel measure of feed efficiency based on the NRC net energy model. Daily records of nutrient intake, fat- and protein-corrected milk yield, body weight, and body condition score were used to model 7 different indices of feed efficiency (gross feed efficiency, return over feed costs, residual feed intake, residual net energy intake, and net energy efficiencies of lactation, maintenance, and body reserve flux) in primiparous Holstein heifers (n = 40) from 0 to 150 d in milk. We observed that feed efficiency was dynamic and fluctuated over the course of the experimental period regardless of the efficiency measure used, but most measures followed similar trends in changes of efficiency on a day-to-day basis. Of the measures evaluated, return over feed costs would be most suited for management purposes due to practicality and simplicity, whereas other measures, such as the efficiency of net energy flux to and from body reserves, present potential targets for inclusion in breeding goals. Net energy efficiency of lactation was estimated to have higher repeatability compared with residual feed or energy intake, which may reduce the time required to estimate whole-lactation efficiency. Because decreased net energy balance was strongly correlated with increased efficiency for all indices evaluated, future work investigating the associations of feed efficiency with the magnitude and duration of negative energy balance is warranted.

Symposium review: The choice and collection of new relevant phenotypes for fertility selection.

In dairy production, high fertility contributes to herd profitability by achieving greater production and maintaining short calving intervals. Improved management practices and genetic selection have contributed to reversing negative trends in dairy cow fertility, but further progress is still required. Phenotypes included in current genetic evaluations are largely interval and binary traits calculated from insemination and calving date records. Several indicator traits such as calving, health, variation in body condition score, and longevity traits also apply to genetic improvement of fertility. Several fertility traits are included in the selection indices of many countries, but for improved selection, the development of novel phenotypes that more closely describe the physiology of reproduction and limit management bias could be more effective. Progesterone-based phenotypes can be determined from milk samples to describe the heritable interval from calving to corpus luteum activity, as well as additional measures of cow cyclicity. A fundamental component of artificial insemination practices is the observation of estrus. Novel phenotypes collected on estrous activity could be used to select for cows clearly displaying heat, as those cows are more likely to be inseminated at the right time and therefore have greater fertility performance. On-farm technologies, including in-line milk testing and activity monitors, may allow for phenotyping novel traits on large numbers of animals. Additionally, selection for improved fertility using traditional traits could benefit from refined and accurate recording and implementation of parameters such as pregnancy confirmation and reproductive management strategy, to differentiate embryonic or fetal loss, and to ensure selection for reproductive capability without producer intervention. Opportunities exist to achieve genetic improvement of reproductive efficiency in cattle using novel phenotypes, which is required for long-term sustainability of the dairy cattle population and industry.

Genetics and genomics of reproductive disorders in Canadian Holstein cattle.

In Canada, reproductive disorders known to affect the profitability of dairy cattle herds have been recorded by producers on a voluntary basis since 2007. Previous studies have shown the feasibility of using producer-recorded health data for genetic evaluations. Despite low heritability estimates and limited availability of phenotypic information, sufficient genetic variation has been observed for those traits to indicate that genetic progress, although slow, can be achieved. Pedigree- and genomic-based analyses were performed on producer-recorded health data of reproductive disorders, including retained placenta (RETP), metritis (METR), and cystic ovaries (CYST) using traditional BLUP and single-step genomic BLUP. Genome-wide association studies and functional analyses were carried out to unravel significant genomic regions and biological pathways, and to better understand the genetic mechanisms underlying RETP, METR, and CYST. Heritability estimates (posterior standard deviation in parentheses) were 0.02 (0.003), 0.01 (0.004), and 0.02 (0.003) for CYST, METR, and RETP, respectively. A moderate to strong genetic correlation of 0.69 (0.102) was found between METR and RETP. Averaged over all traits, sire proof reliabilities increased by approximately 11 percentage points with the incorporation of genomic data using a multiple-trait linear model. Biological pathways and associated genes underlying the studied traits were identified and will contribute to a better understanding of the biology of these 3 health disorders in dairy cattle.

Estimating the effect of the deleterious recessive haplotypes AH1 and AH2 on reproduction performance of Ayrshire cattle.

The effects of 2 deleterious recessive haplotypes on reproduction performance of Ayrshire cattle, Ayrshire Haplotype 1 (AH1) and Ayrshire Haplotype 2 (AH2), were investigated in Canadian Ayrshire cattle. We calculated their phenotypic effects on stillbirth (SB) rate and 56-d nonreturn rate (NRR) by estimating the interaction of service sire carrier status with maternal grandsire carrier status using the official Canadian evaluation models for those 2 traits. The interaction term included 9 subclasses for the 3 possible statuses of each bull: haplotype carrier, noncarrier, or not genotyped. For AH1, 394 carriers and 1,433 noncarriers were available, whereas 313 carriers and 1,543 noncarriers were available for the AH2 haplotype. The number of matings considered for SB was 34,312 for heifers (first parity) and 115,935 for cows (later parities). For NRR, 49,479 matings for heifers and 160,528 for cows were used to estimate the haplotype effects. We observed a negative effect of AH1 on SB rates, which was 2.0% higher for matings of AH1-carrier sires to dams that had an AH1-carrier sire; this effect was found for both heifers and cows. However, AH1 had small, generally nonsignificant effects on NRR. The AH2 haplotype had a substantial negative effect on NRR, with 5.1% more heifers and 4.0% more cows returning to service, but the effects on SB rates were inconsistent and mostly small effects. Our results validate the harmful effects of AH1 and AH2 on reproduction traits in the Canadian Ayrshire population. This information will be of great interest for the dairy industry, allowing producers to make mating decisions that would reduce reproductive losses.

Use of a single-step approach for integrating foreign information into national genomic evaluation in Holstein cattle.

The use of multi-trait across-country evaluation (MACE) and the exchange of genomic information among countries allows national breeding programs to combine foreign and national data to increase the size of the training populations and potentially increase accuracy of genomic prediction of breeding values. By including genotyped and nongenotyped animals simultaneously in the evaluation, the single-step genomic BLUP (GBLUP) approach has the potential to deliver more accurate and less biased genomic evaluations. A single-step genomic BLUP approach, which enables integration of data from MACE evaluations, can be used to obtain genomic predictions while avoiding double-counting of information. The objectives of this study were to apply a single-step approach that simultaneously includes domestic and MACE information for genomic evaluation of workability traits in Canadian Holstein cattle, and compare the results obtained with this methodology with those obtained using a multi-step approach (msGBLUP). By including MACE bulls in the training population, msGBLUP led to an increase in reliability of genomic predictions of 4.8 and 15.4% for milking temperament and milking speed, respectively, compared with a traditional evaluation using only pedigree and phenotypic information. Integration of MACE data through a single-step approach (ssGBLUPIM) yielded the highest reliabilities compared with other considered methods. Integration of MACE data also helped reduce bias of genomic predictions. When using ssGBLUPIM, the bias of genomic predictions decreased by half compared with msGBLUP using domestic and MACE information. Therefore, the reliability and bias of genomic predictions for both traits improved substantially when a single-step approach was used for evaluation compared with a multi-step approach. The use of a single-step approach with integration of MACE information provides an alternative to the current method used in Canadian genomic evaluations.

Invited review: Determination of large-scale individual dry matter intake phenotypes in dairy cattle.

Feed efficiency has been widely studied in many areas of dairy science and is currently seeing renewed interest in the field of breeding and genetics. A critical part of determining how efficiently an animal utilizes feed is accurately measuring individual dry matter (DM) intake. Currently, multiple methods are used to measure feed intake or determine the DM content of that feed, resulting in different levels of accuracy of measurement. Furthermore, the scale at which data need to be collected for use in genetic analyses makes some methodologies impractical. This systematic review aims to provide an overview of the current methodologies used to measure both feed intake in ruminants and DM content of feedstuffs, current methods to predict individual DM intake, and applications of large-scale intake measurements. Overall, advances in milk spectral data analysis present a promising method of estimating individual DM intake on a herd scale with further validation of prediction models. Although measurements of individual feed intake rely on the same underlying principle, the methods selected are largely dictated by the costs of capital, labor, and necessary analyses. Finally, DM methodologies were synthesized into a comprehensive protocol for use in a variety of feedstuffs.