Genomic prediction of lactation curves for milk,fat, protein,and somatic cell score in Holstein cattle

Application of random regression models (RRM) in a 2-step genomic prediction might be a feasible way to select young animals based on the complete pattern of the lactation curve. In this context, the prediction reliability and bias of genomic estimated breeding value (GEBV) for milk, fat, and protein yields and somatic cell score over days in milk (DIM) using a 2-step genomic approach were investigated. In addition, the effect of including cows in the training and validation populations was investigated. Estimated breeding values for each DIM (from 5 to 305 d) from the first 3 lactations of Holstein animals were deregressed and used as pseudophenotypes in the second step. Individual additive genomic random regression coefficients for each trait were predicted using RRM and genomic best linear unbiased prediction and further used to derive GEBV for each DIM. Theoretical reliabilities of GEBV obtained by the RRM were slightly higher than theoretical reliabilities obtained by the accumulated yield up to 305 d (P305). However, validation reliabilities estimated for GEBV using P305 were higher than for GEBV using RRM. For all traits, higher theoretical and validation reliabilities were estimated when incorporating genomic information. Less biased GEBV estimates were found when using RRM compared with P305, and different validation reliability and bias patterns for GEBV over time were observed across traits and lactations. Including cows in the training population increased the theoretical reliabilities and bias of GEBV; nonetheless, the inclusion of cows in the validation population does not seem to affect the regression coefficients and the theoretical reliabilities. In summary, the use of RRM in 2-step genomic prediction produced fairly accurate GEBV over the entire lactation curve for all analyzed traits. Thus, selecting young animals based on the pattern of lactation curves seems to be a feasible alternative in genomic selection of Holstein cattle for milk production traits.     

Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.     

Genome-wide association for milk production traits and somatic cell score in different lactation stages of Ayrshire,Holstein, and Jersey dairy cattle

We performed genome-wide association analyses for milk, fat, and protein yields and somatic cell score based on lactation stages in the first 3 parities of Canadian Ayrshire, Holstein, and Jersey cattle. The genome-wide association analyses were performed considering 3 different lactation stages for each trait and parity: from 5 to 95, from 96 to 215, and from 216 to 305 d in milk. Effects of single nucleotide polymorphisms (SNP) for each lactation stage, trait, parity, and breed were estimated by back-solving the direct breeding values estimated using the genomic best linear unbiased predictor and single-trait random regression test-day models containing only the fixed population average curve and the random genomic curves. To identify important genomic regions related to the analyzed lactation stages, traits, parities and breeds, moving windows (SNP-by-SNP) of 20 adjacent SNP explaining more than 0.30% of total genetic variance were selected for further analyses of candidate genes. A lower number of genomic windows with a relatively higher proportion of the explained genetic variance was found in the Holstein breed compared with the Ayrshire and Jersey breeds. Genomic regions associated with the analyzed traits were located on 12, 8, and 15 chromosomes for the Ayrshire, Holstein, and Jersey breeds, respectively. Especially for the Holstein breed, many of the identified candidate genes supported previous reports in the literature. However, well-known genes with major effects on milk production traits (e.g., diacylglycerol O-acyltransferase 1) showed contrasting results among lactation stages, traits, and parities of different breeds. Therefore, our results suggest evidence of differential sets of candidate genes underlying the phenotypic expression of the analyzed traits across breeds, parities, and lactation stages. Further functional studies are needed to validate our findings in independent populations.     

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

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

Single-step genome-wide association for longitudinal traits of Canadian Ayrshire,Holstein, and Jersey dairy cattle

Estimating single nucleotide polymorphism (SNP) effects over time is essential to identify and validate candidate genes (or quantitative trait loci) associated with time-dependent variation of economically important traits and to better understand the underlying mechanisms of lactation biology. Therefore, in this study, we aimed to estimate time-dependent effects of SNP and identifying candidate genes associated with milk (MY), fat (FY), and protein (PY) yields, and somatic cell score (SCS) in the first 3 lactations of Canadian Ayrshire, Holstein, and Jersey breeds, as well as suggest their potential pattern of phenotypic effect over time. Random regression coefficients for the additive direct genetic effect were estimated for each animal using single-step genomic BLUP, based on 2 random regression models: one considering MY, FY, and PY in the first 3 lactations and the other considering SCS in the first 3 lactations. Thereafter, SNP solutions were obtained for random regression coefficients, which were used to estimate the SNP effects over time (from 5 to 305 d in lactation). The top 1% of SNP that showed a high magnitude of SNP effect in at least 1 d in lactation were selected as relevant SNP for further analyses of candidate genes, and clustered according to the trajectory of their SNP effects over time. The majority of SNP selected for MY, FY, and PY increased the magnitude of their effects over time, for all breeds. In contrast, for SCS, most selected SNP decreased the magnitude of their effects over time, especially for the Holstein and Jersey breeds. In general, we identified a different set of candidate genes for each breed, and similar genes were found across different lactations for the same trait in the same breed. For some of the candidate genes, the suggested pattern of phenotypic effect changed among lactations. Among the lactations, candidate genes (and their suggested phenotypic effect over time) identified for the second and third lactations were more similar to each other than for the first lactation. Well-known candidate genes with major effects on milk production traits presented different suggested patterns of phenotypic effect across breeds, traits, and lactations in which they were identified. The candidate genes identified in this study can be used as target genes in studies of gene expression.     

Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo

Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.     

Genetic and environmental relationships between body condition score and milk production traits in Canadian Holsteins

The objective of this research was to estimate genetic parameters of first-lactation body condition score (BCS), milk yield, fat percentage (Fat%), protein percentage (Prot%), somatic cell score (SCS), milk urea nitrogen (MUN), lactose percentage (Lact%), and fat to protein ratio (F:P) using multiple-trait random regression animal models. Changes in covariances between BCS and milk production traits on a daily basis have not been investigated before and could be useful for determining which BCS estimated breeding values (EBV) might be practical for selection in the future. Field staff from Valacta milk recording agency (Sainte-Anne-de-Bellevue, QC, Canada) collected BCS from Québec herds several times per cow throughout the lactation. Average daily heritabilities and genetic correlations among the various traits were similar to literature values. On an average daily basis, BCS was genetically unfavorably correlated with milk yield (i.e., increased milk yield was associated with lower body condition). The unfavorable genetic correlation between BCS and milk yield became stronger as lactation progressed, but was equivalent to zero for the first month of lactation. Favorable genetic correlations were found between BCS with Prot%, SCS, and Lact% (i.e., greater BCS was associated with greater Prot%, lower SCS, and greater Lact%). These correlations were strongest in early lactation. On an average daily basis, BCS was not genetically correlated with Fat% or MUN, but was negatively correlated with F:P. Furthermore, BCS at 5 and 50 d in milk (DIM) had the most favorable genetic correlations with milk production traits over the lactation (at 5, 50, 150, and 250 DIM). Thus, early lactation BCS EBV shows potential for selection. Regardless, this study showed that the level of association BCS has with milk production traits is not constant over the lactation. Simultaneous selection for both BCS and milk production traits should be considered, mainly due to the unfavorable genetic correlation between BCS with milk yield.     

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.     

Relationships between milk yield and somatic cell score in Canadian Holsteins from simultaneous and recursive random regression models

Multiple-trait random regression animal models with simultaneous and recursive links between phenotypes for milk yield and somatic cell score (SCS) on the same test day were fitted to Canadian Holstein data. All models included fixed herd test-day effects and fixed regressions within region-age at calving-season of calving classes, and animal additive genetic and permanent environmental regressions with random coefficients. Regressions were Legendre polynomials of order 4 on a scale from 5 to 305 d in milk (DIM). Bayesian methods via Gibbs sampling were used for the estimation of model parameters. Heterogeneity of structural coefficients was modeled across (the first 3 lactations) and within (4 DIM intervals) lactation. Model comparisons in terms of Bayes factors indicated the superiority of simultaneous models over the standard multiple-trait model and recursive parameterizations. A moderate heterogeneous (both across- and within-lactation) negative effect of SCS on milk yield (from −0.36 for 116 to 265 DIM in lactation 1 to −0.81 for 5 to 45 DIM in lactation 3) and a smaller positive reciprocal effect of SCS on milk yield (from 0.007 for 5 to 45 DIM in lactation 2 to 0.023 for 46 to 115 DIM in lactation 3) were estimated in the most plausible specification. No noticeable differences among models were detected for genetic and environmental variances and genetic parameters for the first 2 regression coefficients. The curves of genetic and permanent environmental variances, heritabilities, and genetic and phenotypic correlations between milk yield and SCS on a daily basis were different for different models. Rankings of bulls and cows for 305-d milk yield, average daily SCS, and milk lactation persistency remained the same among models. No apparent benefits are expected from fitting causal phenotypic relationships between milk yield and SCS on the same test day in the random regression test-day model for genetic evaluation purposes.     

Genetic correlations between the cumulative pseudo-survival rate,milk yield,and somatic cell score during lactation in Holstein cattle in Japan using a random regression model

Trends in genetic correlations between longevity, milk yield, and somatic cell score (SCS) during lactation in cows are difficult to trace. In this study, changes in the genetic correlations between milk yield, SCS, and cumulative pseudo-survival rate (PSR) during lactation were examined, and the effect of milk yield and SCS information on the reliability of estimated breeding value (EBV) of PSR were determined. Test day milk yield, SCS, and PSR records were obtained for Holstein cows in Japan from 2004 to 2013. A random subset of the data was used for the analysis (825 herds, 205,383 cows). This data set was randomly divided into 5 subsets (162–168 herds, 83,389–95,854 cows), and genetic parameters were estimated in each subset independently. Data were analyzed using multiple-trait random regression animal models including either the residual effect for the whole lactation period (H0), the residual effects for 5 lactation stages (H5), or both of these residual effects (HD). Milk yield heritability increased until 310 to 351 d in milk (DIM) and SCS heritability increased until 330 to 344 DIM. Heritability estimates for PSR increased with DIM from 0.00 to 0.05. The genetic correlation between milk yield and SCS increased negatively to under ?0.60 at 455 DIM. The genetic correlation between milk yield and PSR increased until 342 to 355 DIM (0.53–0.57). The genetic correlation between the SCS and PSR was ?0.82 to ?0.83 at around 180 DIM, and decreased to ?0.65 to ?0.71 at 455 DIM. The reliability of EBV of PSR for sires with 30 or more recorded daughters was 0.17 to 0.45 when the effects of correlated traits were ignored. The maximum reliability of EBV was observed at 257 (H0) or 322 (HD) DIM. When the correlations of PSR with milk yield and SCS were considered, the reliabilities of PSR estimates increased to 0.31–0.76. The genetic parameter estimates of H5 were the same as those for HD. The rank correlation coefficients of the EBV of PSR between H0 and H5 or HD were greater than 0.9. Additionally, the reliabilities of EBV of PSR of H0 were similar to those for H5 and HD. Therefore, the genetic parameter estimates in H0 were not substantially different from those in H5 and HD. When milk yield and SCS, which were genetically correlated with PSR, were used, the reliability of PSR increased. Estimates of the genetic correlations between PSR and milk yield and between PSR and SCS are useful for management and breeding decisions to extend the herd life of cows.     

Genome-wide association study for selected cheese-making properties in Dual-Purpose Belgian Blue cows

This study aimed to estimate genetic parameters and identify genomic region(s) associated with selected cheese-making properties (CMP) in Dual-Purpose Belgian Blue (DPBB) cows. Edited data were 46,301 test-day records of milk yield, fat percentage, protein percentage, casein percentage, milk calcium content (CC), coagulation time (CT), curd firmness after 30 min from rennet addition (a30), and milk titratable acidity (MTA) collected from 2014 to 2020 on 4,077 first-parity (26,027 test-day records), and 3,258 second-parity DPBB cows (20,274 test-day records) distributed in 124 herds in the Walloon Region of Belgium. Data of 28,266 SNP, located on 29 Bos taurus autosomes (BTA) of 1,699 animals were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 25 consecutive SNPs (with an average size of ～2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Heritability estimates for the included CMP ranged from 0.19 (CC) to 0.50 (MTA), and 0.24 (CC) to 0.41 (MTA) in the first and second parity, respectively. The genetic correlation estimated between CT and a30 varied from ?0.61 to ?0.41 and from ?0.55 to ?0.38 in the first and second lactations, respectively. Negative genetic correlations were found between CT and milk yield and composition, while those estimated between curd firmness and milk composition were positive. Genome-wide association analyses results identified 4 genomic regions (BTA1, BTA3, BTA7, and BTA11) associated with the considered CMP. The identified genomic regions showed contrasting results between parities and among the different stages of each parity. It suggests that different sets of candidate genes underlie the phenotypic expression of the considered CMP between parities and lactation stages of each parity. The findings of this study can be used for future implementation and use of genomic evaluation to improve the cheese-making traits in DPBB cows.     

Genomic-based genetic parameters for milkability traits derived from automatic milking systems in North American Holstein cattle

The number of dairy farms adopting automatic milking systems (AMS) has considerably increased around the world aiming to reduce labor costs, improve cow welfare, increase overall performance, and generate a large amount of daily data, including production, behavior, health, and milk quality records. In this context, this study aimed to (1) estimate genomic-based variance components for milkability traits derived from AMS in North American Holstein cattle based on random regression models; and (2) derive and estimate genetic parameters for novel behavioral indicators based on AMS-derived data. A total of 1,752,713 daily records collected using 36 milking robot stations and 70,958 test-day records from 4,118 genotyped Holstein cows were used in this study. A total of 57,600 SNP remained after quality control. The daily-measured traits evaluated were milk yield (MY, kg), somatic cell score (SCS, score unit), milk electrical conductivity (EC, mS), milking efficiency (ME, kg/min), average milk flow rate (FR, kg/min), maximum milk flow rate (FRM, kg/min), milking time (MT, min), milking failures (MFAIL), and milking refusals (MREF). Variance components and genetic parameters for MY, SCS, ME, FR, FRM, MT, and EC were estimated using the AIREMLF90 software under a random regression model fitting a third-order Legendre orthogonal polynomial. A threshold Bayesian model using the THRGIBBS1F90 software was used for genetically evaluating MFAIL and MREF. The daily heritability estimates across days in milk (DIM) ranged from 0.07 to 0.28 for MY, 0.02 to 0.08 for SCS, 0.38 to 0.49 for EC, 0.45 to 0.56 for ME, 0.43 to 0.52 for FR, 0.47 to 0.58 for FRM, and 0.22 to 0.28 for MT. The estimates of heritability (± SD) for MFAIL and MREF were 0.02 ± 0.01 and 0.09 ± 0.01, respectively. Slight differences in the genetic correlations were observed across DIM for each trait. Strong and positive genetic correlations were observed among ME, FR, and FRM, with estimates ranging from 0.94 to 0.99. Also, moderate to high and negative genetic correlations (ranging from −0.48 to −0.86) were observed between MT and other traits such as SCS, ME, FR, and FRM. The genetic correlation (± SD) between MFAIL and MREF was 0.25 ± 0.02, indicating that both traits are influenced by different sets of genes. High and negative genetic correlations were observed between MFAIL and FR (−0.58 ± 0.02) and MFAIL and FRM (−0.56 ± 0.02), indicating that cows with more MFAIL are those with lower FR. The use of random regression models is a useful alternative for genetically evaluating AMS-derived traits measured throughout the lactation. All the milkability traits evaluated in this study are heritable and have demonstrated selective potential, suggesting that their use in dairy cattle breeding programs can improve dairy production efficiency in AMS.     

Short communication: Genetic parameters of production traits in Chinese Holsteins using a random regression test-day model

The objective of this study was to estimate genetic parameters of production traits in the first 3 parities in Chinese Holsteins. Data were a random sample of complete herds (109,005 test-day records of 9,706 cows from 54 herds) extracted from the original data set, which included 362,304 test-day records of 30,942 Holstein cows from 105 herds. A test-day animal model with multiple-trait random regression and the Gibbs sampling method were used for parameter estimation. Regression curves were modeled using Legendre polynomials of order 4. The multiple-trait analysis included milk, fat, and protein yield, and somatic cell score (SCS). Average daily heritabilities ranged between 0.222 and 0.346 for the yield traits and between 0.092 and 0.187 for SCS. Heritabilities were higher in the third lactation for all traits. Within-parity genetic correlations were very high among the yield traits (>0.806) and were close to zero between SCS and yield traits, especially for first-parity cows. Results were similar to previous literature estimates from studies that used the same model as applied to this study. The estimates found in this study will be used to perform breeding value estimation for national genetic evaluations in Chinese Holsteins.     

Genetic correlation patterns between somatic cell score and protein yield in the Italian Holstein-Friesian population

Genetic parameters for somatic cell score (SCS) in the Italian Holstein-Friesian population were estimated addressing the pattern of genetic correlation with protein yield in different parities (first, second, and third) and on different days in milk within each parity. Three approaches for parameter estimation were applied using random samples of herds from the national database of the Italian Holstein Association. Genetic correlations for lactation measures (305-d protein yield and lactation SCS) were positive in the first parity (0.31) and close to zero in the second (0.01) and third (0.09) parities. These results indicated that larger values of SCS were genetically associated with increased production. The second and third sets of estimates were based on random regression test-day models, modeling the shape of lactation curve with the Wilmink function and fourth-order Legendre polynomials, respectively. Genetic correlations from both random regression models showed a specific pattern associated with days in milk within and across parities. Estimates varied from positive to negative in the first and second parity, and from null to negative in the third parity. Patterns were similar for both random regression models. The average overall correlation between SCS and protein yield was zero or slightly positive in the first lactation and ranged from zero to negative in later lactations. Correlation estimates differed by parity and stage of lactation. They also demonstrated the dubiousness of applying a single genetic correlation measure between SCS and protein in setting selection strategies. Differences in magnitude and the sign of genetic correlations between SCS and yields across and within parities should be accounted for in selection schemes.     

Best prediction of yields for long lactations

Lactation records of any reasonable length now can be processed with the selection index method known as best prediction (BP). Previous prediction programs were limited to the 305-d standard used since 1935. Best prediction was implemented in 1998 to calculate lactation records in USDA genetic evaluations, replacing the test interval method used since 1969 to calculate lactation records. Best prediction is more complex but also more accurate, particularly when testing is less frequent. Programs were reorganized to output better graphics, give users simpler access to options, and provide additional output, such as BP of daily yields. Test-day data for 6 breeds were extracted from the national dairy database, and lactation lengths were required to be ≥500 d (Ayrshire, Milking Shorthorn) or ≥800 d (all others). Average yield and SD at any day in milk (DIM) were estimated by fitting 3-parameter Wood's curves (milk, fat, protein) and 4-parameter exponential functions (somatic cell score) to means and SD of 15- (≤300 DIM) and 30-d (>300 DIM) intervals. Correlations among TD yields were estimated using an autoregressive matrix to account for biological changes and an identity matrix to model daily measurement error. Autoregressive parameters (r) were estimated separately for first (r = 0.998) and later parities (r = 0.995). These r values were slightly larger than previous estimates due to the inclusion of the identity matrix. Correlations between traits were modified so that correlations between somatic cell score and other traits may be nonzero. The new lactation curves and correlation functions were validated by extracting TD data from the national database, estimating 305-d yields using the original and new programs, and correlating those results. Daily BP of yield were validated using daily milk weights from on-farm meters in university research herds. Correlations ranged from 0.900 to 0.988 for 305-d milk yield. High correlations ranged from 0.844 to 0.988 for daily yields, although correlations were as low as 0.015 on d 1 of lactation, which may be due to calving-related disorders that are not accounted for by BP. Correlations between 305-d yield calculated using 50-d intervals from 50 to 250 DIM and 305-yield calculated using all TD to 500 DIM increased as TD data accumulated. Many cows can profitably produce for >305 DIM, and the revised program provides a flexible tool to model these records.     

Climate sensitivity of milk production traits and milk fatty acids in genotyped Holstein dairy cows

The aim of this study was the evaluation of climate sensitivity via genomic reaction norm models [i.e., to infer cow milk production and milk fatty acid (FA) responses on temperature-humidity index (THI) alterations]. Test-day milk traits were recorded between 2010 and 2016 from 5,257 first-lactation genotyped Holstein dairy cows. The cows were kept in 16 large-scale cooperator herds, being daughters of 344 genotyped sires. The longitudinal data consisted of 47,789 test-day records for the production traits milk yield (MY), fat yield (FY), and protein yield (PY), and of 20,742 test-day records for 6 FA including C16:0, C18:0, saturated fatty acids (SFA), unsaturated fatty acids (UFA), monounsaturated fatty acids (MUFA), and polyunsaturated fatty acids (PUFA). After quality control of the genotypic data, 41,057 SNP markers remained for genomic analyses. Meteorological data from the weather station in closest herd distance were used for the calculation of maximum hourly daily THI. Genomic reaction norm models were applied to estimate genetic parameters in a single-step approach for production traits and FA in dependency of THI at different lactation stages, and to evaluate the model stability. In a first evaluation strategy (New_sire), all phenotypic records from daughters of genotyped sires born after 2010 were masked, to mimic a validation population. In the second strategy (New_env), only daughter records of the new sires recorded in the most extreme THI classes were masked, aiming at predicting sire genomic estimated breeding values (GEBV) under heat stress conditions. Model stability was the correlation between GEBV of the new sires in the reduced data set with respective GEBV estimated from all phenotypic data. Among all test-day production traits, PY responded as the most sensitive to heat stress. As observed for the remaining production traits, genetic variances were quite stable across THI, but genetic correlations between PY from temperate climates with PY from extreme THI classes dropped to 0.68. Genetic variances in dependency of THI were very similar for C16:0 and SFA, indicating marginal climatic sensitivity. In the early lactation stage, genetic variances for C18:0, MUFA, PUFA, and UFA were significantly larger in the extreme THI classes compared with the estimates under thermoneutral conditions. For C18:0 and MUFA, PUFA, and UFA in the middle THI classes, genetic correlations in same traits from the early and the later lactation stages were lower than 0.50, indicating strong days in milk influence. Interestingly, within lactation stages, genetic correlations for C18:0 and UFA recorded at low and high THI were quite large, indicating similar genetic mechanisms under stress conditions. The model stability was improved when applying the New_env instead of New_sire strategy, especially for FA in the first stage of lactation. Results indicate moderately accurate genomic predictions for milk traits in extreme THI classes when considering phenotypic data from a broad range of remaining THI. Phenotypically, thermal stress conditions contributed to an increase of UFA, suggesting value as a heat stress biomarker. Furthermore, the quite large genetic variances for UFA at high THI suggest the consideration of UFA in selection strategies for improved heat stress resistance.     

Processing characteristics of dairy cow milk are moderately heritable

Milk processing attributes represent a group of milk quality traits that are important to the dairy industry to inform product portfolio. However, because of the resources required to routinely measure such quality traits, precise genetic parameter estimates from a large population of animals are lacking for these traits. Milk processing characteristics considered in the present study—rennet coagulation time, curd-firming time, curd firmness at 30 and 60 min after rennet addition, heat coagulation time, casein micelle size, and milk pH—were all estimated using mid-infrared spectroscopy prediction equations. Variance components for these traits were estimated using 136,807 test-day records from 5 to 305 d in milk (DIM) from 9,824 cows using random regressions to model the additive genetic and within-lactation permanent environmental variances. Heritability estimates ranged from 0.18 ± 0.01 (26 DIM) to 0.38 ± 0.02 (180 DIM) for rennet coagulation time; from 0.26 ± 0.02 (5 DIM) to 0.57 ± 0.02 (174 DIM) for curd-firming time; from 0.16 ± 0.01 (30 DIM) to 0.56 ± 0.02 (271 DIM) for curd firmness at 30 min; from 0.13 ± 0.01 (30 DIM) to 0.48 ± 0.02 (271 DIM) for curd firmness at 60 min; from 0.08 ± 0.01 (17 DIM) to 0.24 ± 0.01 (180 DIM) for heat coagulation time; from 0.23 ± 0.02 (30 DIM) to 0.43 ± 0.02 (261 DIM) for casein micelle size; and from 0.20 ± 0.01 (30 DIM) to 0.36 ± 0.02 (151 DIM) for milk pH. Within-trait genetic correlations across DIM weakened as the number of days between compared intervals increased but were mostly >0.4 except between the peripheries of the lactation. Eigenvalues and associated eigenfunctions of the additive genetic covariance matrix for all traits revealed that at least the 80% of the genetic variation among animals in lactation profiles was associated with the height of the lactation profile. Curd-firming time and curd firmness at 30 min were weakly to moderately genetically correlated with milk yield (from 0.33 ± 0.05 to 0.59 ± 0.05 for curd-firming time, and from ?0.62 ± 0.03 to ?0.21 ± 0.06 for curd firmness at 30 min). Milk protein concentration was strongly genetically correlated with curd firmness at 30 min (0.84 ± 0.02 to 0.94 ± 0.01) but only weakly genetically correlated with milk heat coagulation time (?0.27 ± 0.07 to 0.19 ± 0.06). Results from the present study indicate the existence of exploitable genetic variation for milk processing characteristics. Because of possible indirect deterioration in milk processing characteristics due to selection for greater milk yield, emphasis on milk processing characteristics is advised.     

A model for the genetic evaluation of number of clinical mastitis cases per lactation in Czech Holstein cows

Cases of mastitis from 9,550 lactations of 6,242 cows were recorded on 5 farms in the Czech Republic from 1996 to 2008. The number of clinical mastitis (CM) cases per cow adjusted to a lactation length of 305 d was analyzed with 4 linear single-trait animal models and one 3-trait model, which also included lactation mean somatic cell score (SCS) and 305-d milk yield. Factors included in the model of choice were parity, combined effect of herd and a 2-yr calving period, calving season, permanent environmental effect of the cow, and additive genetic effect of the cow. From both the single-trait and multiple-trait models, estimated heritability of number of CM cases was 0.11 (±0.015 for the multiple-trait model). Permanent environmental effects accounted for approximately one-third of the phenotypic variance. Heritability estimates for lactation mean SCS and 305-d milk yield were 0.17 ± 0.019 and 0.25 ± 0.011, respectively, and genetic correlations of these traits with number of CM cases were 0.80 ± 0.059 and 0.34 ± 0.079, respectively. Genetic evaluation of the number of CM cases in Czech Holsteins could be carried out including data from all parities using a 3-trait animal model with SCS and milk yield as additional traits.     

Genomic loci and genetic parameters for uterine diseases in first-parity Holstein cows and associations with milk production and fertility

Based on the clinical stage (e.g., vaginal discharge) and bacterial species, several forms of uterine diseases (UD) exist and can be classified as different traits [i.e., different stages of endometritis (EM) and metritis (MET)], which may differ in their genetic background and causal physiological mechanisms. Consequently, the present study aimed to study (1) the effect of UD on 305-d lactation and fertility, (2) the estimation of heritabilities for UD traits using pedigree- and SNP-based relationships, and (3) genome-wide associations to detect significant SNP markers and to infer candidate genes for UD traits. The data set contained herd manager and veterinarian recorded UD traits of 14,810 first-lactating genotyped Holstein cows from 63 large-scale contract herds. Binary defined UD traits (healthy or diseased) according to the clinical stage were endometritis catarrhalis (EM I), endometritis mucopurulenta (EM II), endometritis purulenta (EM III), pyometra (EM IV), endometritis (EM_SOD; superordinate diagnosis = no specific clinical stage defined), and MET. The binary defined trait UDall included all EM and MET diagnoses. The prevalence of UDall was 26.7%. The effect of UD on 305-d lactation and fertility was estimated via linear and generalized linear mixed models. We applied linear single-trait animal models and threshold models to estimate pedigree- and SNP-based heritabilities for UD traits, and bivariate linear models for genetic correlation estimations between UDall with 305-d lactation and fertility traits. A diagnosis for UDall had significant unfavorable effects on the female fertility traits calving interval, interval from calving to first service, days open, and nonreturn rate after 90 d, but was unrelated to 305-d lactation records for production traits milk yield, protein yield, and fat yield. Heritabilities for UDall and EM stages were close to zero, displaying maximal values of 0.05 for pedigree and 0.07 for SNP-based relationship matrices. For MET, pedigree- and SNP-based heritabilities were <0.001 and 0.07, respectively. Genetic correlations ranged from 0.20 to 0.31 between UDall with 305-d milk, protein, and fat yield, and from 0.17 to 0.40 with fertility traits. The GWAS revealed 5 SNP on bovine chromosomes (BTA) 1, 8, 10, 23 for UDall, 5 SNP on BTA 26 for EM I, 1 SNP on BTA 19 for EM II, 4 SNP on BTA 2, 18, 20, 25 for EM III, and 4 SNP on BTA 4, 16, 20 for EM IV above the significance threshold. For EM_SOD, we identified 15 significantly associated SNP on 4 chromosomes, and 4 significant SNP on BTA 3, 20, 22, 28 for MET. Marker associations for UD traits were annotated to 24 potential candidate genes using the ENSEMBL database. Six of these genes were previously reported to be involved in uterine defense mechanisms or in endometritis. Further detected genes contribute to immune response mechanisms during bacterial infections. Different SNP significantly influenced different UD stages, explaining the inter-individual variations in clinical severity of uterine infections.     

Test-day genetic analysis of condition score and heart girth in Holstein Friesian cows.

This study aimed to estimate heritability for condition score and heart girth using a test-day model, to investigate the genetic relationships between condition score, heart girth, and milk yield traits and to analyze the genetic relationships of condition score and heart girth measured at different stages of lactation. Cows from 25 dairy herds were scored for body condition and measured for heart girth at 3-mo intervals for 2 yr. Approximately 5000 test-day observations on condition score, heart girth, and milk fat and protein yield from 1344 Italian Friesian cows were analyzed using two approaches: 1) repeated observations for a trait were considered repeated measurements of the same trait; 2) observations for a trait collected in different stages of lactation (dry period, 1 to 75, 76 to 130, 131 to 210, and 211 to 300 DIM) were treated as different traits. (Co)variance components and related parameters were estimated using REML multiple-trait procedures and animal models with unequal design for different traits. Heritability estimates for fat and protein test-day yield and for test-day condition score and heart girth were 0.22, 0.18, 0.29, and 0.33, respectively. Condition score was negatively correlated with yield traits and positively correlated with heart girth, whereas genetic relationships between heart girth and milk yield traits were negligible. Heritability estimates were 0.27 for condition score recorded in the first half of lactation (1 to 75 and 76 to 130 DIM), 0.36 for condition score in the second half of lactation (131 to 210 and 211 to 300 DIM) and 0.32 for condition score recorded on dry cows. Genetic correlations between condition scores measured in different lactation stages were generally high (0.85 or more), with the exception of the relationships between the first and the last stage of lactation (0.74) and between the first half of lactation and the dry period (0.7). Heritability estimates for heart girth in different lactation stages ranged from 0.31 to 0.40, and genetic correlations between high girth measured in different lactation stages were higher than 0.80.