Next Article in Journal
Rule Discovery in Milk Content towards Mastitis Diagnosis: Dealing with Farm Heterogeneity over Multiple Years through Classification Based on Associations
Next Article in Special Issue
Detection and Visualization of Heterozygosity-Rich Regions and Runs of Homozygosity in Worldwide Sheep Populations
Previous Article in Journal
Canola Meal versus Soybean Meal as Protein Supplements in the Diets of Lactating Dairy Cows Affects the Greenhouse Gas Intensity of Milk
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Aspects of Somatic Cell Count in Holstein Dairy Cows in Iran

1
Department of Animal Science, Shiraz University, Shiraz 7155713876, Iran
2
Department of Farm Animal Health, University of Utrecht, Yalelaan 7, 3584 CL Utrecht, The Netherlands
*
Author to whom correspondence should be addressed.
Animals 2021, 11(6), 1637; https://doi.org/10.3390/ani11061637
Submission received: 23 April 2021 / Revised: 26 May 2021 / Accepted: 28 May 2021 / Published: 1 June 2021
(This article belongs to the Special Issue Phenotypic and Genotypic Characterization of Farm Animals)

Abstract

:

Simple Summary

Holstein Friesian is the dominant breed of modern dairy cattle in Iran, therefore the development of a genetic evaluation system for economically important traits for this population is needed. Genetic parameters for SCS and its genetic correlations with production traits were estimated in Iranian Holstein cows. The moderate heritability for SCS and its low negative genetic correlations with yield traits indicate that genetic selection for decreasing SCS would have a relatively medium genetic progress with no necessarily antagonistic effects on lactation performance. The estimates found in this study can be considered as the first step to include SCS in the national genetic evaluations in Iranian Holsteins using a multiple-trait, multiple-lactation random regression model.

Abstract

The aim of this study was to estimate the genetic parameters of somatic cell count (SCC) and its relationship with production traits in the first three parities in Iranian Holstein dairy cows. Data were 1,891,559 test-day records of SCC, milk yield, and milk compositions on 276,217 lactations on 147,278 cows distributed in 134 herds. The number of test-day records in the first, second and third parities were 995,788 (on 147,278 cows), 593,848 (on 85,153 cows), and 301,923 (on 43,786 cows), respectively. Test-day SCCs were transformed to somatic cell scores (SCS). A random regression test-day animal model through four-trait three-lactation was used to estimate variance components for test-day records of SCS and lactation traits were included. Gibbs sampling was used to obtain marginal posterior distributions for the various parameters using a single chain of 200,000 iterates in which the first 50,000 iterates of each chain were regarded as a burn-in period. The mean heritability estimates for SCS (0.15 to 0.18) were lower than those for milk yield (0.36 to 0.38), fat yield (0.30 to 0.31), protein yield (0.31 to 0.32), fat percentage (0.21 to 0.25), and protein percentage (0.21 to 0.22). Low negative genetic correlations ranging from −0.05 to −0.30 were found between SCS and yield traits (milk, fat, and protein yields). The genetic correlation found between SCS and fat percentage was close to zero, however, a low positive genetic correlation ranging from 0.12 to 0.17 was found between SCS and protein percentage. Based on the results, it can be concluded that genetic selection for decreasing SCS would also increase lactation yield. The estimates found in this study can be used to perform breeding value estimations for national genetic evaluations in Iranian Holsteins using a multiple-trait, multiple-lactation random regression model.

1. Introduction

There are about seven million cattle in Iran, of which some 700,000 are Holstein distributed in more than 2000 herds. The herd sizes vary from 20 to 3000 dairy cows with a mean 305-day milk production of 7200 kg [1,2]. Although the use of selection indices has been a common practice during the last decades [3,4], in Iran, bull and cow selection is mainly based on the estimated breeding value of milk yield. The only nationally devised genetic index in use includes the production traits and longevity [2]. However, milk price in Iran is based on bulk tank somatic cell count (SCC) thereby affecting revenues. Somatic cell score (SCS), the average of the log of test-day SCC, could be included in the economic breeding index as a proxy for udder health. However, accurate genetic parameters for the prevailing population are required for any trait to be included in a selection index. Different models (e.g., repeatability model, random regression, or fixed regression test-day model) are used for the national genetic evaluation of SCS in Interbull member countries [5]. The first random regression test-day (RR-TDM) model for national genetic evaluation of dairy cattle was adopted at the beginning of the 21st century [6]. Since then, RR-TD models have become the models of choice for genetic evaluation of production traits in dairy cattle [7,8]. The RR-TD model was subsequently adopted for national genetic evaluation of dairy breeds, and most of the countries that participated in the international evaluations for production traits of Holstein bulls used a test-day model for national genetic evaluations, while others used a lactation animal model [5]. Holstein Friesian is the dominant breed of modern dairy cattle in Iran, therefore the development of a genetic evaluation system for economically important traits for this population is needed. The objective of the present study was to estimate genetic parameters for SCS and its genetic correlations with production traits in the first three parities in Iranian Holstein cows using a multiple-trait, multiple-lactation random regression test-day model.

2. Materials and Methods

2.1. Data

Data used in this study were test-day records from Holstein cows collected between January 2009 to December 2018 by the Animal Breeding Center of Iran (Karaj, Iran, 35°49′ N and 50°59′ E). All evaluated herds were purebred Holsteins, managed under conditions similar to those used in most developed countries, and were under official performance and pedigree recording. The diet, fed as a total mixed ration (TMR), mainly consisted of corn silage, alfalfa hay, barley grain, fat powder, beet pulp, and feed additives. Monthly milk recording was performed by trained technicians of the Iranian Animal Breeding Center, according to the guidelines of the International Committee for Animal Recording [9]. Cows with missing birth dates, calving dates, or lactation numbers were excluded. Only records from the first three parities that had data for SCC, milk yield, and milk compositions on a given test-day remained in the dataset. Records from days in milk (DIM) lower than 5 and greater than 365 were eliminated. Within cows, if lactation number three was present, lactations number one and two were also present, and if lactation number two was present, lactation number one was also present. Age at the first calving (AFC) was calculated as the difference between birth date and calving date at the first parity and restricted to the range of 640 to 1200 d. Edited data consisted of 1,891,559 test-day records of milk yield (MY), fat yield (FY), protein yield (PY), fat percentage (FP), protein percentage (PP), and somatic cell count (SCC) on 276,217 lactations on 147,278 cows in 134 herds. The number of test-day records in the first, second, and third parity were 995,788 (on 147,278 cows), 593,848 (on 85,153 cows), and 301,923 (on 43,786 cows), respectively. Full pedigree records were included from 386,955 females and 4219 males. Test-day SCCs were transformed into somatic cell scores (SCS) based on the following equation: SCS = log 2 ( SCC / 100 , 000 ) + 3 .

2.2. Model

The following random regression (RR) test-day animal model through four-trait three-lactation was used to estimate variance components for test-day records of MY, FY, PY, FP, PP, and SCS:
y ijklm =   µ   +   HTDp i + b = 0 4 AS j ø b ( t ) + b = 0 2 a l ø b ( t ) + b = 0 2 pe l ø b ( t ) +   e ijklm
where y ijklm is the test-day record (MY, FY, PY, FP, PP, and SCS) on DIM m of cow l in parity k, belonging to ith class of HTDp, and jth class of AS, HTDp i is the fixed effect of ith class of herd-test-day parity, AS k is the fixed effect of age-season of calving defined as the following: age at calving (five intervals of age at calving were created for each parity) × season of calving (four seasons: winter from Jan–Mar, spring from Apr–Jun, summer from Jul–Sep and autumn from Oct–Dec), b = 0 2 a l ø b ( t ) and b = 0 2 pe l ø b ( t ) are, respectively, the random regression coefficients of additive and permanent environment effects, e ijklm is the residual random effect. Regression curves were modeled using Legendre polynomials of order 2. Residual random effects were assumed to be normally distributed with a mean of 0.0. Homogeneity of residual variance was tested for SCS by computing the standard deviation (SD) of observed residuals (difference between observed and estimated values) for each DIM in the first three parities. The (co)variance components were estimated by Bayesian inference using the Gibbs sampler of the GIBBS3F90 program [10]. Gibbs sampling was used to obtain marginal posterior distributions for the various parameters using a single chain of 200,000 iterates. The first 50,000 iterates of each chain were regarded as a burn-in period to allow sampling from the proper marginal distributions. The length of the burn-in period was determined by visually inspecting plots of sample values across rounds. Daily heritability was defined as the ratio of genetic variance to the sum of genetic, permanent environmental, and residual variances at a given DIM.

3. Results

Descriptive statistics of the dataset used for parameter estimation for all traits included are shown in Table 1. The average level of SCS increased with increasing parity having means (SD) of 2.39 (1.97), 2.99 (2.06), and 3.37 (2.11) in the first, second, and third parities, respectively. The average daily milk yield was 31.6 kg with 3.36% of fat and 3.09% of protein in the first parity. The corresponding values for the second and third parities were 34.7 kg, 3.40%, and 3.13% and 35.8 kg, 3.42%, 3.10%. Figure 1 shows trajectories of SCS along with those for MY, FY, PY, FP, and PP by DIM for each parity. Lower SCS was found at the time of the milk yield peak, then increased for the remainder of the lactation. The lactation curves for MY showed that the first parity was more persistent but the second and third parities had higher peak yields to decline more thereafter compared to the first parity. Lactation curves for FY and PY had the same pattern as those for MY. However, in the latter part of the lactation, average first-parity yield traits (MY, FY, and PY) were higher compared to higher parities.

3.1. Variance Components

The phenotypic, permanent environmental, genetic, and residual variances by DIM for SCS in each parity are presented in Figure 2. Phenotypic, permanent environmental, and residual variance found for SCS increased with increasing parity. The average phenotypic variances by DIM for the first, second, and third parities were, respectively, 3.59, 3.80, and 4.02. The average permanent environmental variances found for SCS by DIM for the first, second, and third parities were, respectively, 0.75, 0.93, and 1.09. The corresponding values for the residual variance were, respectively, 2.10, 2.17, and 2.22. The genetic variance found for SCS in the first parity (0.74) was higher than those found for the second (0.70) and third parities (0.71). A high difference in phenotypic and permanent environmental variances was found between parities at the beginning of the lactation. This difference decreased with increasing DIM and disappeared at the end of the lactation. The observed genetic variance found in the first parity compared to those found in later parities was lower at early lactation, but higher in late lactation. The pattern of genetic variance by DIM found in the second parity was similar to that found for the third parity.

3.2. Heritabilities

The patterns of estimated heritabilities on a daily basis for SCS in the first three parities are represented in Figure 2. The results showed that h2 of SCS in the first parity compared to those for higher parities was lower in early lactation, but higher in mid and late lactation. The heritability pattern found in the second parity was the same as that found in the third parity. The mean (range) heritability estimates across lactation for SCS were 0.18 (0.11 to 0.20), 0.15 (0.11 to 0.17), and 0.15 (0.12 to 0.17), respectively, in the first, second, and third parities.

3.3. Relationships between Traits

The average phenotypic, genetic, and permanent environmental correlations between SCS and production traits are presented in Table 2, Table 3 and Table 4, respectively. The patterns of phenotypic and genetic correlations estimated on a daily basis between SCS and production traits included are presented in Figure 3. Weak negative phenotypic correlations ranging from −0.04 to −0.27 were found between SCS and yield traits (MY, FY, PY); however, weak positive phenotypic correlations ranging from 0.03 to 0.08 were found between SCS and FP and PP. The mean (range) genetic correlations between SCS and MY in the first three parities were, respectively, −0.21 (−0.27 to −0.04), −0.28 (−0.47 to −0.17), and −0.30 (−0.45 to −0.27). The corresponding values found between SCS and FY were −0.10 (−0.13 to −0.03), −0.14 (−0.23 to −0.08), and −0.16 (−0.22 to −0.11), and those found between SCS and PY were −0.05 (−0.09 to 0.11), −0.13 (−0.23 to −0.04), and −0.15 (−0.22 to −0.10). The mean (range) genetic correlations between SCS with PP in the first three parities were, respectively, 0.17 (0.05 to 0.25), 0.12 (0.07 to 0.17), and 0.13 (0.09 to 0.18). The corresponding values found between SCS and FP were −0.01 (−0.06 to 0.02), 0.00 (−0.03 to 0.05), and −0.00 (−0.4 to 0.06). For primiparous cows, the lowest genetic correlation between SCS and MY was found in the mid-lactation. The genetic correlation between SCS and MY in primiparous cows was lower at the mid-lactation compared to correlations found at the beginning and end of the lactation. The genetic correlation between SCS and MY in multiparous cows was high at the beginning of lactation, decreased with increasing DIM, and was the lowest at the end of lactation.

4. Discussion

Lower SCS was found at peak milk yield, then increased for the remainder of the lactation which is in close agreement with previous studies [8,11]. The average level of SCS increased with increasing parity which agrees with previous studies [11,12,13,14]. An increase in the SCS with increasing parity can be attributed to the fact that older cows have more opportunity for exposure to mastitis-causing pathogens [15,16,17]. The average heritabilities estimated for SCS in this study were 0.18, 0.15, and 0.15 for the first, second, and third parity, respectively, which are slightly higher than those reported in previous studies [8,18,19]. Miglior et al. [8] reported that the average daily heritabilities for SCS in Chinese Holstein were 0.09, 0.15, and 0.19, for the first, second, and third parity, respectively. Haile-Mariam et al. [12] reported that average daily heritability for SCS in Australian dairy cattle ranged from 0.09 to 0.11 in the first three parities. The results showed that average phenotypic correlations between SCS and yield traits ranged from −0.04 to −0.27. The genetic correlations found between milk yield and SCS ranged from −0.21 to −0.30, in the first three parities. It is well documented that the correlation between milk yield and SCS varies both within and across parities [11]. The negative relationships found between SCS and yield traits are consistent with the deleterious effect of poor mammary health on production [20]. Franzoi et al. [21] reported that an increase in SCC impairs lactation performance. An elevated SCS may indicate a bacterial infection of the udder which could have an adverse effect on milk production explaining the negative relationship between SCS and MY found in this study. In addition, higher MY may cause a decrease in SCS at the same level of infection through dilution effects, assuming that the SCS is proportional to the degree of infection. Miglior et al. [8] reported that the genetic correlations between milk yield and SCS were −0.12, −0.20, and −0.36, respectively, in the first, second, and third parities in Chinese Holstein cows. Koivula et al. [22] reported a positive genetic correlation between SCC and yield traits in the first parity, and near zero in the second parity in Ayrshire dairy cows in Finland.

5. Conclusions

Genetic parameters were estimated for SCS in Iranian Holstein dairy cows using a random regression test-day model. SCS had a low negative correlation with the production traits, but its magnitude increased in later parities. The moderate heritability for SCS and its low negative genetic correlations with yield traits indicates that genetic selection for decreasing SCS would have a relatively medium genetic progress with no necessarily antagonistic effects on lactation performance. The estimates found in this study can be considered as the first step to include SCS in the national genetic evaluations in Iranian Holsteins using a multiple-trait, multiple-lactation random regression model.

Author Contributions

Conceptualization, H.A. and M.H.; methodology, H.A.; formal analysis, H.A.; data curation, H.A.; writing—original draft preparation, H.A.; writing—review and editing, M.H.; supervision, H.A.; project administration, H.A.; funding acquisition, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the vice chancellor of research of Shiraz University (grant number: 97GR1M155290, 2014).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Relevant information supporting the results not presented here is provided in supplemental files (https://github.com/hadiatashi/Atashi-et-al-2021-SCS/SCS.rar) (26 May 2021).

Acknowledgments

The author greatly appreciates the cooperation of the Animal Breeding Center of Iran (Karaj, Iran) in providing the data.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Atashi, H.; Zamiri, M.J.; Sayyadnejad, M.B.; Akhlaghi, A. Trends in the reproductive performance of Holstein dairy cows in Iran. Trop. Anim. Health Prod. 2012, 44, 2001–2006. [Google Scholar] [CrossRef] [PubMed]
  2. Sadeghi-Sefidmazgi, A.; Moradi-Shahrbabak, M.; Nejati-Javaremi, A.; Shadparvar, A. Estimation of economic values in three breeding perspectives for longevity and milk production traits in Holstein dairy cattle in Iran. Ital. J. Anim. Sci. 2009, 8, 359–375. [Google Scholar] [CrossRef] [Green Version]
  3. Miglior, F.; Muir, B.; Van Doormaal, B. Selection indices in Holstein cattle of various countries. J. Dairy Sci. 2005, 88, 1255–1263. [Google Scholar] [CrossRef]
  4. Van Raden, P. Selection of dairy cattle for lifetime profit. In Proceedings of the 7th World Congress on Genetics Applied to Livestock Production, Montpellier, France, 19–23 August 2002; pp. 127–130. [Google Scholar]
  5. Mark, T.; Fikse, W.; Emanuelson, U.; Philipsson, J. International genetic evaluations of Holstein sires for milk somatic cell and clinical mastitis. J. Dairy Sci. 2002, 85, 2384–2392. [Google Scholar] [CrossRef]
  6. Schaeffer, L.; Jamrozik, J.; Kistemaker, G.; Van Doormaal, J. Experience with a test-day model. J. Dairy Sci. 2000, 83, 1135–1144. [Google Scholar] [CrossRef]
  7. Lidauer, M.H.; Pösö, J.; Pedersen, J.; Lassen, J.; Madsen, P.; Mäntysaari, E.A.; Nielsen, U.S.; Eriksson, J.-Å.; Johansson, K.; Pitkänen, T. Across-country test-day model evaluations for Holstein, Nordic Red Cattle, and Jersey. J. Dairy Sci. 2015, 98, 1296–1309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Miglior, F.; Gong, W.; Wang, Y.; Kistemaker, G.; Sewalem, A.; Jamrozik, J. Genetic parameters of production traits in Chinese Holsteins using a random regression test-day model. J. Dairy Sci. 2009, 92, 4697–4706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. ICAR. ICAR (International Committee for Animal Recording). International Agreement of Recording Practices. 2012. Available online: http://aberdeenangus.ro/wp-content/uploads/2014/03/ICAR.pdf (accessed on 26 May 2021).
  10. Aguilar, I.; Tsuruta, S.; Masuda, Y.; Lourenco, D.; Legarra, A.; Misztal, I. BLUPF90 suite of programs for animal breeding with focus on genomics. In Proceedings of the World Congress on Genetics Applied to Livestock Production, Auckland, New Zealand, 10–15 February 2018; pp. 11–16. [Google Scholar]
  11. Jamrozik, J.; Bohmanova, J.; Schaeffer, L. Relationships between milk yield and somatic cell score in Canadian Holsteins from simultaneous and recursive random regression models. J. Dairy Sci. 2010, 93, 1216–1233. [Google Scholar] [CrossRef] [PubMed]
  12. Haile-Mariam, M.; Goddard, M.; Bowman, P. Estimates of genetic parameters for daily somatic cell count of Australian dairy cattle. J. Dairy Sci. 2001, 84, 1255–1264. [Google Scholar] [CrossRef]
  13. Alam, M.; Cho, C.; Choi, T.; Park, B.; Choi, J.; Choy, Y.; Lee, S.; Cho, K. Estimation of genetic parameters for Somatic Cell Scores of Holsteins using multi-trait lactation models in Korea. Asian-Australas. J. Anim. Sci. 2015, 28, 303. [Google Scholar] [CrossRef] [PubMed]
  14. Martins, A.; Silvestre, A.; Petim-Batista, M.; Colaço, J. Somatic cell score genetic parameter estimates of dairy cattle in Portugal using fractional polynomials. J. Anim. Sci. 2011, 89, 1281–1285. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, Z.; Li, X.; Yang, F.; Luo, J.; Wang, X.; Liu, L.; Li, H. Influences of season, parity, lactation, udder area, milk yield, and clinical symptoms on intramammary infection in dairy cows. J. Dairy Sci. 2016, 99, 6484–6493. [Google Scholar] [CrossRef] [PubMed]
  16. Abebe, R.; Hatiya, H.; Abera, M.; Megersa, B.; Asmare, K. Bovine mastitis: Prevalence, risk factors and isolation of Staphylococcus aureus in dairy herds at Hawassa milk shed, South Ethiopia. BMC Vete. Res. 2016, 12, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Girma, S.; Mammo, A.; Bogele, K.; Sori, T.; Tadesse, F.; Jibat, T. Study on prevalence of bovine mastitis and its major causative agents in West Harerghe zone, Doba district, Ethiopia. J. Vet. Med. Anim. Health 2012, 4, 116–123. [Google Scholar]
  18. Costa, A.; Lopez-Villalobos, N.; Visentin, G.; De Marchi, M.; Cassandro, M.; Penasa, M. Heritability and repeatability of milk lactose and its relationships with traditional milk traits, somatic cell score and freezing point in Holstein cows. Animal 2019, 13, 909–916. [Google Scholar] [CrossRef]
  19. Soumri, N.; Carabaño, M.J.; González-Recio, O.; Bedhiaf-Romdhani, S. Genetic parameters of somatic cell scores using random regression test-day models with Legendre polynomials in Tunisian dairy cattle. Livest. Sci. 2020, 241, 104178. [Google Scholar] [CrossRef]
  20. Rajala-Schultz, P.; Gröhn, Y.; McCulloch, C.; Guard, C. Effects of clinical mastitis on milk yield in dairy cows. J. Dairy Sci. 1999, 82, 1213–1220. [Google Scholar] [CrossRef]
  21. Franzoi, M.; Manuelian, C.; Penasa, M.; De Marchi, M. Effects of somatic cell score on milk yield and mid-infrared predicted composition and technological traits of Brown Swiss, Holstein Friesian, and Simmental cattle breeds. J. Dairy Sci. 2020, 103, 791–804. [Google Scholar] [CrossRef] [PubMed]
  22. Koivula, M.; Mäntysaari, E.; Negussie, E.; Serenius, T. Genetic and phenotypic relationships among milk yield and somatic cell count before and after clinical mastitis. J. Dairy Sci. 2005, 88, 827–833. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Lactation curve for daily milk yield (MY), SCS, fat yield (FY), protein yield (PY), fat percentage (FP), and protein percentage (PP) in the first three parities.
Figure 1. Lactation curve for daily milk yield (MY), SCS, fat yield (FY), protein yield (PY), fat percentage (FP), and protein percentage (PP) in the first three parities.
Animals 11 01637 g001
Figure 2. Phenotypic (A), permanent environmental (B), additive genetic (C), residual variance (D), and heritability (E) estimated for SCS over days in milk in the first three parities.
Figure 2. Phenotypic (A), permanent environmental (B), additive genetic (C), residual variance (D), and heritability (E) estimated for SCS over days in milk in the first three parities.
Animals 11 01637 g002
Figure 3. Phenotypic and genetic correlations estimated between SCS and productive traits over days in milk in the first three parities. Phenotypic correlation in the first parity (A), second parity (B), and third parity (C). Genetic correlation in the first parity (D), second parity (E), and third parity (F). SCSMY, SCSFY, SCSPY, SCSFP, and SCSPP, are abbreviations for somatic cell score and milk yield (SCSMY), fat yield (SCSFY), protein yield (SCSPY), fat percentage (SCSFP), and protein percentage (SCSPP), respectively.
Figure 3. Phenotypic and genetic correlations estimated between SCS and productive traits over days in milk in the first three parities. Phenotypic correlation in the first parity (A), second parity (B), and third parity (C). Genetic correlation in the first parity (D), second parity (E), and third parity (F). SCSMY, SCSFY, SCSPY, SCSFP, and SCSPP, are abbreviations for somatic cell score and milk yield (SCSMY), fat yield (SCSFY), protein yield (SCSPY), fat percentage (SCSFP), and protein percentage (SCSPP), respectively.
Animals 11 01637 g003
Table 1. Descriptive statistics for milk yield and compositions for the first three parities in Holstein dairy cows in Iran.
Table 1. Descriptive statistics for milk yield and compositions for the first three parities in Holstein dairy cows in Iran.
First ParitySecond ParityThird Parity
TraitsMeanSDCV (%)MeanSDCV (%)MeanSDCV (%)
Milk yield31.67.62434.710.43035.811.332
Fat yield 110.493.343211.634.343712.094.7239
Protein yield 29.702.382410.733.102910.963.3531
Fat percentage3.360.83253.400.88263.420.8926
Protein percentage3.090.39133.130.41133.100.4113
SCS 32.391.97822.992.06683.372.1163
1 Daily fat yield (hg), 2 daily protein yield (hg), 3 SCS = log2 (SCC/100,000) + 3.
Table 2. Mean (range) phenotypic correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
Table 2. Mean (range) phenotypic correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
First ParitySecond ParityThird Parity
Milk yield−0.19 (−0.00 to −0.00)−0.25 (−0.00 to −0.00)−0.27 (−0.00 to −0.00)
Fat yield−0.04 (−0.00 to −0.00)−0.08 (−0.00 to −0.00)−0.09 (−0.00 to −0.00)
Protein yield−0.06 (−0.00 to −0.00)−0.10 (−0.00 to −0.00)−0.11 (−0.00 to −0.00)
Fat percentage0.05 (−0.00 to −0.00)0.04 (−0.00 to −0.00)0.03 (−0.00 to −0.00)
Protein percentage0.08 (−0.00 to −0.00)0.08 (−0.00 to −0.00)0.08 (−0.00 to −0.00)
Table 3. Mean (range) genetic correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
Table 3. Mean (range) genetic correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
First ParitySecond ParityThird Parity
Milk yield−0.21 (−0.27 to −0.04)−0.28 (−0.47 to −0.17)−0.30 (−0.45 to −0.27)
Fat yield−0.10 (−0.13 to −0.03)−0.14 (−0.23 to −0.08)−0.16 (−0.22 to −0.11)
Protein yield−0.05 (−0.09 to 0.11)−0.13 (−0.23 to −0.04)−0.15 (−0.22 to −0.10)
Fat percentage−0.01 (−0.06 to 0.02)0.00 (−0.03 to 0.05)0.00 (−0.0 to 0.09)
Protein percentage0.17 (0.05 to 0.25)0.12 (0.07 to 0.17)0.13 (0.09 to 0.18)
Table 4. Mean (range) permanent environmental correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
Table 4. Mean (range) permanent environmental correlations between SCS and production traits in the first three parities of Iranian Holstein cows.
First ParitySecond ParityThird Parity
Milk yield−0.16(−0.29 to −0.10)−0.21(−0.42 to −0.10)−0.24(−0.42 to −0.10)
Fat yield−0.04(−0.11 to −0.00)−0.10(−0.21 to −0.01)−0.11(−0.21 to −0.03)
Protein yield−0.06(−0.11 to −0.02)−0.09(−0.21 to −0.02)−0.10(−0.19 to −0.02)
Fat percentage0.06(0.02 to 0.09)0.02(0.00 to 0.04)0.02(−0.01 to 0.05)
Protein percentage0.05(0.03 to 0.08)0.07(0.05 to 0.08)0.03(0.05 to 0.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Atashi, H.; Hostens, M. Genetic Aspects of Somatic Cell Count in Holstein Dairy Cows in Iran. Animals 2021, 11, 1637. https://doi.org/10.3390/ani11061637

AMA Style

Atashi H, Hostens M. Genetic Aspects of Somatic Cell Count in Holstein Dairy Cows in Iran. Animals. 2021; 11(6):1637. https://doi.org/10.3390/ani11061637

Chicago/Turabian Style

Atashi, Hadi, and Miel Hostens. 2021. "Genetic Aspects of Somatic Cell Count in Holstein Dairy Cows in Iran" Animals 11, no. 6: 1637. https://doi.org/10.3390/ani11061637

APA Style

Atashi, H., & Hostens, M. (2021). Genetic Aspects of Somatic Cell Count in Holstein Dairy Cows in Iran. Animals, 11(6), 1637. https://doi.org/10.3390/ani11061637

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop