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.

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.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing.

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Breeding for reduced methane emission and feed-efficient Holstein cows: An international response.

Selecting for lower methane (CH4) emitting animals is one of the best approaches to reduce CH4 given that genetic progress is permanent and cumulative over generations. As genetic selection requires a large number of animals with records and few countries actively record CH4, combining data from different countries could help to expedite accurate genetic parameters for CH4 traits and build a future genomic reference population. Additionally, if we want to include CH4 in the breeding goal, it is important to know the genetic correlations of CH4 traits with other economically important traits. Therefore, the aim of this study was first to estimate genetic parameters of 7 suggested methane traits, as well as genetic correlations between methane traits and production, maintenance, and efficiency traits using a multicountry database. The second aim was to estimate genetic correlations within parities and stages of lactation for CH4. The third aim was to evaluate the expected response of economically important traits by including CH4 traits in the breeding goal. A total of 15,320 methane production (MeP, g/d) records from 2,990 cows belonging to 4 countries (Canada, Australia, Switzerland, and Denmark) were analyzed. Records on dry matter intake (DMI), body weight (BW), body condition score, and milk yield (MY) were also available. Additional traits such as methane yield (MeY; g/kg DMI), methane intensity (MeI; g/kg energy-corrected milk), a genetic standardized methane production, and 3 definitions of residual methane production (g/d), residual feed intake, metabolic BW (MBW), BW change, and energy-corrected milk were calculated. The estimated heritability of MeP was 0.21, whereas heritability estimates for MeY and MeI were 0.30 and 0.38, and for the residual methane traits heritability ranged from 0.13 to 0.16. Genetic correlations between different methane traits were moderate to high (0.41 to 0.97). Genetic correlations between MeP and economically important traits ranged from 0.29 (MY) to 0.65 (BW and MBW), being 0.41 for DMI. Selection index calculations showed that residual methane had the most potential for inclusion in the breeding goal when compared with MeP, MeY, and MeI, as residual methane allows for selection of low methane emitting animals without compromising other economically important traits. Inclusion of residual feed intake in the breeding goal could further reduce methane, as the correlation with residual methane is moderate and elicits a favorable correlated response. Adding a negative economic value for methane could facilitate a substantial reduction in methane emissions while maintaining an increase in milk production.

Development and comparison of RNA-sequencing pipelines for more accurate SNP identification: practical example of functional SNP detection associated with feed efficiency in Nellore beef cattle.

BACKGROUND: Optimization of an RNA-Sequencing (RNA-Seq) pipeline is critical to maximize power and accuracy to identify genetic variants, including SNPs, which may serve as genetic markers to select for feed efficiency, leading to economic benefits for beef production. This study used RNA-Seq data (GEO Accession ID: PRJEB7696 and PRJEB15314) from muscle and liver tissue, respectively, from 12 Nellore beef steers selected from 585 steers with residual feed intake measures (RFI; n = 6 low-RFI, n = 6 high-RFI). Three RNA-Seq pipelines were compared including multi-sample calling from i) non-merged samples; ii) merged samples by RFI group, iii) merged samples by RFI and tissue group. The RNA-Seq reads were aligned against the UMD3.1 bovine reference genome (release 94) assembly using STAR aligner. Variants were called using BCFtools and variant effect prediction (VeP) and functional annotation (ToppGene) analyses were performed. RESULTS: On average, total reads detected for Approach i) non-merged samples for liver and muscle, were 18,362,086.3 and 35,645,898.7, respectively. For Approach ii), merging samples by RFI group, total reads detected for each merged group was 162,030,705, and for Approach iii), merging samples by RFI group and tissues, was 324,061,410, revealing the highest read depth for Approach iii). Additionally, Approach iii) merging samples by RFI group and tissues, revealed the highest read depth per variant coverage (572.59 ± 3993.11) and encompassed the majority of localized positional genes detected by each approach. This suggests Approach iii) had optimized detection power, read depth, and accuracy of SNP calling, therefore increasing confidence of variant detection and reducing false positive detection. Approach iii) was then used to detect unique SNPs fixed within low- (12,145) and high-RFI (14,663) groups. Functional annotation of SNPs revealed positional candidate genes, for each RFI group (2886 for low-RFI, 3075 for high-RFI), which were significantly (P < 0.05) associated with immune and metabolic pathways. CONCLUSION: The most optimized RNA-Seq pipeline allowed for more accurate identification of SNPs, associated positional candidate genes, and significantly associated metabolic pathways in muscle and liver tissues, providing insight on the underlying genetic architecture of feed efficiency in beef cattle.

Targeted genotyping to identify potential functional variants associated with cholesterol content in bovine milk.

High blood cholesterol concentration, mainly caused by high dietary cholesterol, is a potential risk factor for human health. Dairy products are important sources of human dietary cholesterol intake. Therefore, monitoring bovine milk cholesterol concentration is important for human health benefit. Genetic selection for improvement of cow milk cholesterol content requires understanding of the genetics of milk cholesterol. For this purpose, we performed analyses of additive and dominance effects of 126 potentially functional SNPs within 43 candidate genes with milk cholesterol content [expressed as mg of cholesterol in 100 g of fat (CHL_fat) or in 100 mg of milk (CHL_milk)]. The additive and dominance effects of SNPs rs380643365 in AGPAT1 (P = 0.04) and rs134357240 in SOAT1 (P = 0.035) genes associated significantly with CHL_fat. Moreover, five (rs109326954 and rs523413537 in DGAT1, rs109376747 in LDLR, rs42781651 in FAM198B and rs109967779 in ACAT2) and four (rs137347384 in RBM19, rs109376747 in LDLR, rs42016945 in PPARG and rs110862179 in SCAP) SNPs were significantly associated with CHL_milk (P < 0.05) based on additive and dominance effect analyses respectively. Rs109326954 and rs523413537 in DGAT1 explained a considerable portion of the phenotypic variance of CHL_milk (7.54 and 6.84% respectively), and might be useful in selection programs for reduced milk cholesterol content. Several significantly associated SNPs were in genes (such as ACAT2 and LDLR) involved in cholesterol metabolism in the liver or cholesterol transport, suggesting multiple mechanisms regulating milk cholesterol content. Nine and seven SNPs identified by additive or dominance effect analyses associated significantly with milk yield and fat yield respectively. Further analyses are required to better understand the consequences of these variants and their potential use in genomic selection of the studied traits.

Short communication: Time-dependent genetic parameters and single-step genome-wide association analyses for predicted milk fatty acid composition in Ayrshire and Jersey dairy cattle.

Milk fat content and fatty acid (FA) composition have great economic value to the dairy industry as they are directly associated with taste and chemical-physical characteristics of milk and dairy products. In addition, consumers' choices are not only based on the nutritional aspects of food, but also on products known to promote better health. Milk FA composition is also related to the metabolic status and physiological stages of cows and thus can also be used as indicator for other novel traits of interest (e.g., metabolic diseases and methane yield). Genetic selection is a promising alternative to manipulate milk FA composition. In this study, we aimed to (1) estimate time-dependent genetic parameters for 5 milk FA groups (i.e., short-chain, medium-chain, long-chain, saturated, and unsaturated) predicted based on milk mid-infrared spectroscopy, for Canadian Ayrshire and Jersey breeds, and (2) conduct a time-dependent, single-step genome-wide association study to identify genomic regions, candidate genes, and metabolic pathways associated with milk FA. We analyzed 31,709 test-day records of 9,648 Ayrshire cows from 268 herds, and 34,341 records of 11,479 Jersey cows from 883 herds. The genomic database contained a total of 2,330 Ayrshire and 1,019 Jersey animals. The average daily heritability ranged from 0.18 (long-chain FA) to 0.34 (medium-chain FA) in Ayrshire, and from 0.25 (long-chain and unsaturated FA) to 0.52 (medium-chain and saturated FA) in Jersey. Important genomic regions were identified in Bos taurus autosomes BTA3, BTA5, BTA12, BTA13, BTA14, BTA16, BTA18, BTA20, and BTA21. The proportion of the variance explained by 20 adjacent SNP ranged from 0.71% (saturated FA) to 1.11% (long-chain FA) in Ayrshire, and from 0.70% (unsaturated FA) to 3.09% (medium-chain FA) in Jersey cattle. Important candidate genes and pathways were also identified, such as the PTK2 and TRAPPC9 genes, associated with milk fat percentage, and HMGCS, FGF10, and C6 genes, associated with fertility traits and immune response. Our findings on the genetic parameters and candidate genes contribute to a better understanding of the genetic architecture of milk FA composition in Ayrshire and Jersey dairy cattle.

Genomic analyses for predicted milk fatty acid composition throughout lactation in North American Holstein cattle.

Milk fat composition has important implications in the nutritional and processing properties of milk. Additionally, milk fat composition is associated with cow physiological and health status. The main objectives of this study were (1) to estimate genetic parameters for 5 milk fatty acid (FA) groups (i.e., short-chain, medium-chain, long-chain, saturated, and unsaturated) predicted from milk infrared spectra using a large data set; (2) to predict genomic breeding values using a longitudinal single-step genomic BLUP approach; and (3) to conduct a single-step GWAS aiming to identify genomic regions, candidate genes, and metabolic pathways associated with milk FA, and consequently, to understand the underlying biology of these traits. We used 629,769 test-day records of 201,465 first-parity Holstein cows from 6,105 herds. A total of 8,865 genotyped (Illumina BovineSNP50K BeadChip, Illumina, San Diego, CA) animals were considered for the genomic analyses. The average daily heritability ranged from 0.24 (unsaturated FA) to 0.47 (medium-chain and saturated FA). The reliability of the genomic breeding values ranged from 0.56 (long-chain fatty acid) to 0.74 (medium-chain fatty acid) when using the default τ and ω scaling parameters, whereas it ranged from 0.58 (long-chain fatty acid) to 0.73 (short-chain fatty acid) when using the optimal τ and ω values (i.e., τ = 1.5 and ω = 0.6), as defined in a previous study in the same population. Relevant chromosomal regions were identified in Bos taurus autosomes 5 and 14. The proportion of the variance explained by 20 adjacent single nucleotide polymorphisms ranged from 0.71% (saturated FA) to 15.12% (long-chain FA). Important candidate genes and pathways were also identified. In summary, our results contribute to a better understanding of the genetic architecture of predicted milk FA in dairy cattle and reinforce the relevance of using genomic information for genetic analyses of these traits.

Symposium review: Multiple-trait single-step genomic evaluation for hoof health.

Hoof lesions represent an important issue in modern dairy herds, with reported prevalence in different countries ranging from 40 to 70%. This high prevalence of hoof lesions has both economic and social consequences, resulting in increased labor expenses and decreasing animal production, longevity, reproduction, health, and welfare. Therefore, a key goal of dairy herds is to reduce the incidence of hoof lesions, which can be achieved both by improving management practices and through genetic selection. The Canadian dairy industry has recently released a hoof health sub-index. This national genetic evaluation program for hoof health was achieved by creating a centralized data collection system that routinely transfers data recorded by hoof trimmers into a coherent and sustainable national database. The 8 most prevalent lesions (digital dermatitis, interdigital dermatitis, interdigital hyperplasia, heel horn erosion, sole hemorrhage, sole ulcer, toe ulcer, and white line lesion) in Canada are analyzed with a multiple-trait model using a single-step genomic BLUP method. Estimated genomic breeding values for each lesion are combined into a sub-index according to their economic value and prevalence. In addition, data recorded within this system were used to create an interactive management report for dairy producers by Canadian DHI, including the prevalence of lesions on farm, their trends over time, and benchmarks with provincial and national averages.

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.

Genome-wide association analysis for ß-hydroxybutyrate concentration in Milk in Holstein dairy cattle.

BACKGROUND: Ketosis in dairy cattle has been shown to cause a high morbidity in the farm and substantial financial losses to dairy farmers. Ketosis symptoms, however, are difficult to identify, therefore, the amount of ketone bodies (mainly ß-hydroxybutyric acid, BHB) is used as an indicator of subclinical ketosis in cows. It has also been shown that milk BHB concentrations have a strong correlation with ketosis in dairy cattle. Mid-infrared spectroscopy (MIR) has recently became a fast, cheap and high-throughput method for analyzing milk components. The aim of this study was to perform a genome-wide association study (GWAS) on the MIR-predicted milk BHB to identify genomic regions, genes and pathways potentially affecting subclinical ketosis in North American Holstein dairy cattle. RESULTS: Several significant regions were identified associated with MIR-predicted milk BHB concentrations (indicator of subclinical ketosis) in the first lactation (SCK1) and second and later lactations (SCK2) in Holstein dairy cows. The strongest association was located on BTA6 for SCK1 and BTA14 on SCK2. Several SNPs on BTA6 were identified in regions and variants reported previously to be associated with susceptibility to ketosis and clinical mastitis in Jersey and Holstein dairy cattle, respectively. One highly significant SNP on BTA14 was found within the DGAT1 gene with known functions on fat metabolism and inflammatory response in dairy cattle. A region on BTA6 and three SNPs on BTA20 were found to overlap between SCK1 and SCK2. However, a novel region on BTA20 (55-63 Mb) for SCK2 was also identified, which was not reported in previous association studies. Enrichment analysis of the list of candidate genes within the identified regions for MIR-predicted milk BHB concentrations yielded molecular functions and biological processes that may be involved in the inflammatory response and lipid metabolism in dairy cattle. CONCLUSIONS: The results of this study confirmed several SNPs and genes identified in previous studies as associated with ketosis susceptibility and immune response, and also found a novel region that can be used for further analysis to identify causal variations and key regulatory genes that affect clinical/ subclinical ketosis.

Genetic and genomic analyses of embryo production in dairy cattle.

The Canadian dairy industry has been using invivo and invitro assisted reproductive technologies to produce embryos. Technological improvements have helped increase the number and quality of embryos produced, but genetic and genomic tools for improving these traits have yet to be assessed for the Canadian Holstein population. Genetic parameters and a genome-wide association study were performed in Canadian Holstein for the total number of embryos (NE) and the number of viable embryos (VE). Results showed potential for genetic selection for both NE and VE, with heritability estimates (± s.e.) of approximately 0.15±0.01. Genetic correlations between the number of embryos produced using different procedures (invivo and invitro) suggested that a similar number of embryos should be expected from a donor regardless of the procedure used. A region on chromosome 11 of the bovine genome was found to be significantly associated with the number of embryos, indicating a potential regulatory role of this region on embryo production. Overall, these findings are of interest for the Canadian dairy industry because they provide useful information for breeders that are interested in producing embryos from the elite donors in their herds or in the population using assisted reproductive technologies.

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.

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.

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0.

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.