Genetic parameters for daily milk somatic cell score and relationships with yield traits of primiparous Holstein cattle in Iran

Background Despite the importance of relationships between somatic cell score (SCS) and currently selected traits (milk, fat and protein yield) of Holstein cows, there was a lack of comprehensive literature for it in Iran. Therefore we tried to examine heritabilities and relationships between these traits using a fixed-regression animal model and Bayesian inference. The data set consisted of 1,078,966 test-day observations from 146,765 primiparous daughters of 1930 sires, with calvings from 2002 to 2013. Results Marginal posterior means of heritability estimates for SCS (0.03 ± 0.002) were distinctly lower than those for milk (0.204 ± 0.006), fat (0.096 ± 0.004) and protein (0.147 ± 0.005) yields. In the case of phenotypic correlations, the relationships between production and SCS were near zero at the beginning of lactation but become increasingly negative as days in milk increased. Although all environmental correlations between production and SCS were negative (−0.177 ± 0.007, −0.165 ± 0.008 and −0.152 ± 0.007 between SCS and milk, fat, and protein yield, respectively), slightly antagonistic genetic correlations were found; with posterior mean of relationships ranging from 0.01 ± 0.039 to 0.11 ± 0.036. This genetic opposition was distinctly higher for protein than for fat. Conclusion Although small, the positive genetic correlations suggest some genetic antagonism between desired increased milk production and reduced SCS (i.e., single-trait selection for increased milk production will also increase SCS).


Background
Selection has traditionally focused on production traits. Today it is generally accepted that undesirable genetic associations (i.e., positive) exist between production and health disorders, including mastitis [1][2][3]. Consequently, it is believed that selection for increased milk production would increase susceptibility to both clinical-and subclinical-mastitis [1,2,4]. Mastitis is the most common and the costly disease of dairy cattle (for review, see Heringstad et al. [2]). In addition to the cost of veterinary treatments, decrease in milk yield quality, early culling, extra labor, increased disease risk in the future and replacement costs, mastitis contributes to consumer concerns [2]. Reduced animal welfare, economic losses, and poor milk quality, consequently, form incentives to reduce the incidence of mastitis.
The SCC (concentration of somatic cells) of cow's milk is an economically important trait in dairy cattle because it can be used both as an indicator of mastitis and as a measure of response to infection [2]. Studies [1,5,6] suggest that genetic correlation of clinical mastitis and SCC is moderate to high, implying that selection for low SCC will reduce the incidence of mastitis. In addition, high SCC milk is undesirable for processors because it reduces the shelf-life of consumer milk [7] and diminishes the quality and quantity of milk protein, thereby reducing cheese yields [8]. SCC also has an economic value in their own right in industries where there are milk price penalties for milk supplied for processing with very high counts. High SCC is not only associated with udder health and milk losses, but also negatively affects the longevity [9] and fertility [10] of dairy cows. Thus, incorporation of SCC in selection decisions seems desirable. Commonly, SCC is log-transformed to SCS (standing for somatic cell score).
Possibilities for simultaneous improvement in SCS and yield traits are dependent on the relationship between the traits. Early estimates of the genetic correlation between these traits, especially from the Iran, were few and based on small data sets. Therefore, the objective of the current study is to estimate heritabilities and genetic, phenotypic and environmental relationships between SCS and currently selected traits (milk, fat and protein yield) on test-day basis using fixed regression animal models (FR-AM) and Bayesian inference in the first parity for Holstein cows.

Methods
A total of 1,078,966 test-day SCC and production traits (milk, fat, and protein yield) records from 146,765 first lactations of Holstein cows recorded from 2002 to 2013 were analyzed. The cows were descendants of 1930 sires (average 76 cows/sire), were distributed over 268 herds (average 547 cows/herd), with age at first calving ranging from 18 to 36 months (in 19 classes). The data were obtained from the Animal Breeding Center of Iran. The test-day yield was limited to ranges of 3 to 75 kg of milk yield, 1.5 to 8 % fact content, 1 to 7 % protein content, and 1 (× 10 3 cells/mL) to 500 (× 10 3 cells/mL) SCC; allowable length of lactation was from 5 to 305 d. Testday SCC were converted to SCS using a base 2 logarithmic function: SCS = log 2 (SCC/100) + 3 [11]. Only data of animals with at least five animal-observations were kept. A total of 15,517 contemporary groups (herd within year-month of test day) were formed. All known pedigree information of the cows was traced back as far as possible, resulting in relationship matrix of 273,078 animals for the final data set. The structure of the data set after editing is summarized in Table 1.
The analyses were conducted using a multi-trait FR-AM. The model involved the direct additive, permanent environmental and residual random effects. The fixed part of the lactation curve was modeled with a herd by year-month of test-day effect (HYM) and Legendre polynomials with six coefficients for the age-season of calving effect (AS; 76 classes) and with five coefficients for the herd-year of calving effect (HY; 1802 classes). The equation presentation of the model is given by: where y ijklmn was the nth test-day record of the mth cow for a trait i (milk, fat, protein yield, or SCS); HYM ij was the jth herd-year-month of test-day effect for a trait i; β ikq was the qth fixed regression coefficient for a trait i specific to the kth age-season of calving class; δ ilq was the qth fixed regression coefficient for a trait i specific to the lth herd-year of calving class; q was the number of covariates; z q (t) was a vector of covariates of size q describing the shape of lactation curve of fixed regressions evaluated at t days in milk; α im and ρ im , respectively, were the additive genetic and permanent environmental effects of cow m for a trait i; and e was the residual. The (co)variance components were estimated by Bayesian inference using the Gibbs sampler of the GIBBS3F90 program [12]. A chain length of 200,000 cycles was established, with a burn-in period of 10,000 cycles, and a sampling interval of ten cycles, corresponding to 19,000 samples for subsequent analysis. Convergence of Gibbs chains was monitored by visual inspections of trace plots. The software R [13] was employed for drawing samples from posterior distributions of parameters.

Result
Phenotypic correlations between yield traits and SCS on daily basis are in Fig. 1. The relationships were near zero at the beginning of lactation but become increasingly negative as days in milk increased. Average phenotypic correlations between SCS and milk, fat, and protein yields were −0.16, −0.06, and −0.09, respectively.

Genetic parameters
Means and posterior standard errors of distributions for variances of random effects are presented in Table 2. Lower genetic variability was observed for protein and fat yields, respectively. Additive genetic variances were always lower than the other variances for all traits. Realizations of heritabilities (mean and 95 % pointwise credible interval) from a Gibbs chain with 190,000 cycles are  Fig. 2. For brevity, only trace plots and marginal posterior densities of heritability for SCS and milk yield are shown. Heritabilities for other two production traits (i.e., fat and protein yield) were intermediate to those for SCS and milk yield and were therefore omitted from Fig. 2. The plots indicate that the algorithm mixed well, in spite of differences among traits. In particular, the mixing of the Gibbs sampler was slightly worse for milk, compared to SCS. Means of the posterior densities of heritability for milk, fat and protein yields were 0.204, 0.096, and 0.147, respectively, while Monte Carlo Standard Error ranged from 0.004 to 0.006. Mean of the marginal posterior density of heritability for SCS was low and equal to 0.03, with a 95 % Bayesian credibility region ranging from 0.026 to 0.034 ( Fig. 2; upper panel).

Relationships between traits
Posterior estimates of the genetic correlation between the milk production traits were moderate to high (Fig. 3), with the posterior mean (standard error) varied between 0.62 (0.014) and 0.90 (0.004). The highest genetic correlation was between milk and protein, whereas the correlation between milk and fat was the lowest. Estimated posterior environmental correlations were all high, about 0.85 to 0.97, and again the highest estimates were between milk and protein (Fig. 3).
Posterior estimates of the genetic correlation between SCS and production traits were low and often symmetric (Fig. 4). Means of the marginal posterior density (highest posterior density of 95 %) of genetic correlations were all positive, but small and averaged 0.07 (0.006 to 0.139), 0.01 (−0.066 to 0.086), and 0.11 (0.033 to 0.171) for milk, fat, and protein, respectively. Posterior distributions of the environmental correlations between SCS and yield traits had a larger mean (−0.15 to −0.18) and a lower standard error (0.007 to 0.008) than corresponding genetic correlations; milk and protein had the largest and smallest correlation, respectively. All these relationships (with the exception of the genetic correlation between SCS and fat yield) can be considered statistically significant, because their 95 % credible intervals not included zero.

Discussion
Phenotypic correlations between SCS and yield traits were all negative, ranging from 0.0 to −0.24, in agreement with the literature finding that yield traits decreased under mastitis conditions. However phenotypic relationships seemed stranger with advancing days in milk; therefore, the effects of SCS on yields of milk, fat and protein were larger in the late than in the beginning of lactation. This result in contrary to results of most  studies [14,15], which have shown larger correlations during the first 30 d than later. Generally, the negative correlations between yield and SCS are consistent with the deleterious effect of poor mammary health on production [16]. Across production traits, milk yield showed the highest heritability estimates (>0.2). In the earlier study of the Iranian population conducted by Kheirabadi et al. [17], obtained by a random regression test-day model, the daily heritabilities for milk yield were from 0.05 to 0.26 and they were also higher than those for fat (0.03 to 0.12) and protein (0.04 to 0.24). Similar conclusions about heritabilities of milk traits of Holstein cattle in Iranian were reported by Razmkabir [18]. Overall, these results suggest that production traits have enough genetic variation to develop breeding programs. In contrast heritability of SCS was considerable smaller than those for production traits, which in good agreement with other studies [14,15]. Posterior mean of heritability for milk yield assessed by Wu et al. [19] in the primiparous Norwegian Red cows was 0.18, very close to our results (0.204 ± 0.006), while posterior mean of heritability for SCS was significantly higher (0.12) than our results (0.03 ± 0.002). Haile-Mariam et al. [20] obtained a heritability of 0.07 when test-day records are analysed as repeated measurements. Kheirabadi and Alijani [21] also estimated small values of SCS daily heritabilities for primiparous Holstein cows in Iran (0.03 to 0.07).
As a consequence, standard errors of the estimates of genetic correlations between yield traits were negligible (<0.02). These results imply that yield traits are influenced by similar genes. This corroborates the idea that increasing milk yield by selection would have a positive effect on the other economic traits (referring to fat and protein yields). It seems clear from results in this and other studies [6,17,22,23] that for any pair of traits, genetic correlations were consistently lower than corresponding environmental correlations; this is due to the fact that environmental effects among traits were more correlated than additive genetic effects.
Previous results from studies [6,19,23] reporting undesirable genetic correlations of SCS with production traits were confirmed. Mrode and Swanson [1] found, for first lactation, a weighted average genetic correlation between SCC and milk, fat, and protein yields of 0.14 (SE 0.04 to 0.05). Although small, the positive genetic correlations suggest some genetic antagonism between desired increased milk production and reduced SCS (i.e., single-trait selection for increased milk production will also increase SCS). Presumably there is, then, some genetic antagonism between increased milk yield and deteriorate incidence of mastitis as genetic correlations between SCC and mastitis are high [1,5,6]. The magnitude of the correlations, however, varied slightly between traits. In agreement with our results, Carlén et al. [6] found the lowest correlation between SCS and fat (0.17 ± 0.04) and the highest between SCS and protein (0.23 ± 0.04) in first-lactation Holstein cows. In general, the relationships between SCS, milk and protein obtained in the present study were in agreement with the results of Yazgan et al. (2010), who found average genetic correlations of SCS with milk and protein equal to 0.07 and 0.14, respectively. Therefore, these results suggest that an increase in somatic cells occurs with an increase of protein content. As shown in Fig. 4, the posterior distributions of relationships between SCS and production traits for the environmental effect were sharper (i.e., less variation), and posterior means had larger absolute values than the corresponding genetic correlation; this reflected the decline

Conclusions
As is usually found in the literature, the genetic correlations between SCS and yield traits were positive, reflecting an antagonism between production and mastitis resistance, whereas the phenotypic correlations were negative, reflecting an unfavorable effect of mastitis on production. This emphasizes the need for including udder health traits in the breeding goal. Results of the present study indicate the presence of low additive genetic variation for SCS in this population of Holstein cows. Therefore, a relatively low genetic progress will be expected following selection. Despite the low estimates of heritability for SCS, the utilization of SCS in order to improve the resistance to mastitis, and in consequence, an increase the profitability of production system, is recommended.