ISRN Veterinary ScienceVolume 2013 (2013), Article ID 470585, 5 pageshttp://dx.doi.org/10.1155/2013/470585

Research Article

Estimation of Genetic and Phenotypic Parameters for Production Traits and Somatic Cell Count for Jersey Dairy Cattle in Zimbabwe

Edward Missanjo,¹ Venancio Imbayarwo-Chikosi,² and Tinyiko Halimani²

¹Malawi College of Forestry and Wildlife, Private Bag 6, Dedza, Malawi²Department of Animal Science, Faculty of Agriculture, University of Zimbabwe, P.O. Box MP 167, Mount Pleasant, Harare, Zimbabwe

Received 27 May 2013; Accepted 24 June 2013

Academic Editors: M. H. Kogut and S. Whisnant

Copyright © 2013 Edward Missanjo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Genetic and phenotypic parameters for production traits and somatic cell count (SCC) for Jersey dairy cattle in Zimbabwe were estimated. A total of 10986 lactation records were obtained from Zimbabwe Livestock Identification Trust, with cows calving in the period from 1996 to 2008. An ASReml program fitting an animal model was used for the analyses. Heritability estimates for milk yield, fat yield, protein yield, fat percentage, protein percentage, and Log₁₀SCC were 0.30, 0.32, 0.33, 0.42, 0.44, and 0.08, respectively. The corresponding repeatability estimates were 0.39, 0.38, 0.39, 0.49, 0.51, and 0.16, respectively. The genetic and phenotypic correlations between different production traits ranged from −0.86 to 0.95 and from −0.88 to 0.98, respectively. The genetic and phenotypic correlations between production traits and Log₁₀SCC were weak almost nonsignificantly differentl from zero. The results imply that milk traits for Jersey dairy cattle in Zimbabwe are more heritable. Therefore, these traits may be preferred by breeders as selection criteria for development of effective genetic improvement programme.

1. Introduction

Knowledge of genetic parameters is the basis of sound livestock improvement programmes. Estimates of heritabilities and genetic correlations are essential population parameters required in animal breeding research and in design and application of practical animal breeding programmes. Genetic parameters are a characteristic of the population in which they were estimated and may change overtime due to selection and management decisions [1].

The Jersey is one of the dairy breeds found in Zimbabwe. It is the second most important breed of dairy after Holstein-Friesian. It is a small breed of dairy cattle, with mature live weight ranging from 360 to 540 kg. Originally bred on the British Channel Island of Jersey, the breed is popular for the high butterfat content of its milk and the low maintenance costs incurred due to its lower bodyweight, as well as its genial disposition. In addition, Jerseys can thrive on locally produced food. They can tolerate high temperatures; heifers mature more quickly than those of other breeds and can be mated at 13 to 15 months; they produce an average herd milk production of 3500 to 5000 kg [2].

Jersey is an ideal breed for crossbreeding with Bos indicus to produce a hardy, disease-tolerant, dairy-type cow which does not need a high plane of nutrition to produce reasonable milk yield and is suited to dairying in the communal areas [1]. However, there has been a rapid decrease in average milk yield per cow for Zimbabwean Jersey from 5000 kg to 3900 kg from the year 2000 to 2010 [3]. Furthermore, there is paucity of information on production traits of Jersey cattle in Zimbabwe as most researches have been concentrated on Holstein-Friesian breed [4]. In addition, there have been massive changes in the dairy farming industry since 2000 due to the land reform programme. This may have resulted in the conversion of use of some dairy farms into crop agriculture; a subsequent decrease in herd sizes; and changes to a number of management and production parameters within the herds. Imports of semen and replacement stock were also affected. This may have altered variance components and hence the genetic parameters for yield in most dairy herds including the Jersey herds.

It thus becomes imperative to genetically evaluate the Jersey production traits as a way of increasing productivity through selection. The objective of this study was to analyse lactation records in order to fit an appropriate model for evaluation of the Zimbabwean Jersey dairy cattle regarding production traits and somatic cell count (SCC). Genetic and phenotypic parameters were also estimated. This led to the development of a selection index [5].

2. Materials and Methods

2.1. Environment

Zimbabwe is located in Southern Africa in the tropical savannah region. The total land area is 390,759 km², and it is divided into five agroecological regions. Rainfall patterns and crop production progressively deteriorate from region I to V. However, livestock production including dairying is practised in all the regions. In the regions with low rainfall, dairying is assisted by production of drought-resistant fodder crops. Most dairy farms are located within 40 km of the major cities and towns [6].

3. Data and Data Edits

The standard 305-day milk production records of pure bred Jersey were obtained from Zimbabwe Livestock Identification Trust (LIT). The file contained production records of 25 herds, with cows calving during the period 1996 to 2008. The following production traits were recorded: 305-day milk, fat and protein yields; milk yield over the complete lactation period; length of lactation; calving interval; days dry; days open and somatic cell count. Identity of herd, cow, dam, sire, date of birth, and date of calving were also recorded.

The records were edited using Statistical Analysis of Systems version 9.1.3 [7] and Pedigree Viewer version 6.4 b [8]. Edits included removal of herds with no information on fat, protein, or milk yields and records where calving date, days dry, days open, date of birth, somatic cell count, or length of lactation had not been recorded. Calving interval classes of 20-day interval were formed from 300 to 460 days. The first class had a calving interval of zero, and this represented the first lactation cows only. Calving intervals below 300 days were deleted and those above 460 days were pooled. Further, those with 305-day milk yield less than 3000 kg and greater than 6000 kg were deleted. Days dry classes of 15-day interval from 0 to 150 days were also formed. Dry periods greater than 150 days were pooled. Age at calving was classified within parity as shown in Table 1. All duplicate records were deleted in order to remove bias from the data. Further edits imposed a restriction of at least five daughters per sire. This was done to improve estimation of the sire’s breeding values and the accuracy of the estimated heritability since the biggest influence on the accuracy of estimation of heritability is the number of daughters per sire [4]. The restriction also ensured connectedness in the data and reduced bias due to preferential treatment of daughters of a particular sire in a group [1]. Then, four subdata sets were created. Data set 1 was for all lactations. This was for those records of cows with at least two consecutive lactations. Data sets 2 and 3 were for those records of cows in lactations 1 and 2, respectively. Data set 4 was for those records of cows in lactation 3 and above. The characteristics of the data sets for the traits analysed are given in Tables 2 and 3.

Table 1: Restriction on age at calving by lactation.

Lactation
Age at calving (month)

Minimum
Maximum

1
22
36

2
34
48

3+
46
126

Table 2: Characteristics of the data sets.

Parity

All
1
2
3+

Records
6725
1057
1117
4551

Sires
135
133
135
132

Dams
1093
1005
1054
1001

Herds
21
21
21
21

Paternal family
175

Maternal family
14

Means and standard deviations in parenthesis

Milk (kg)
4468 (533)
4420 (527)
4528 (528)
4464 (535)

Fat (kg)
199 (36)
194 (35)
223 (41)
194 (33)

Protein (kg)
158 (25)
155 (25)
159 (25)
158 (26)

Fat%
4.42 (0.37)
4.35 (0.28)
4.89 (0.39)
4.33 (0.30)

Protein%
3.52 (0.15)
3.50 (0.16)
3.50 (0.14)
3.53 (0.15)

Log₁₀SCC
5.54 (0.56)
5.32 (0.57)
5.38 (0.59)
5.50 (0.52)

Table 3: Pedigree structure means by Tier.

Tier
1
2
3
4

Individuals
6.00
277.00
64 43.00
11 60.00

Milk (kg)
—
44 64.84
45 39.21
45 38.83

Fat (kg)
—
198.57
205.45
211.13

Protein (kg)
—
157.78
160.88
161.53

Fat%
—
4.42
4.49
4.63

Protein%
—
3.52
3.53
3.54

Log₁₀SCC
—
5.45
5.34
4.91

Tier 1: grand-grand-grandparents.

Tier 2: grand-grandparents.

Tier 3: grandparents.

Tier 4: parents.

4. Statistical Analysis

Multivariate analyses to estimate variance components, heritability, predicted breeding values (EBVs), and genetic and phenotypic corrections were done using animal models and the ASReml program developed by [9].

Two models were fitted to the data. Model 1 (a repeatability animal model) was for data sets 1 and 4. The model is described as follows: where   is the observed value for all traits (305-day milk yield; fat and protein yields);  is the overall mean common to all observations;  is the fixed effect of herd year of calving-season of calving ;  is the fixed effect of classes of calving interval ;  is the fixed effect of classes of days dry ;  and are the linear and quadratic regression coefficients, respectively, on age at calving in months;  is the random animal effect of cow ;  is the random permanent environmental effect of cow ; and  is the random residual effects, .

Model 2 was used for data sets 2 and 3. This model was similar to model 1. The only difference is that for data set 2 fixed effects of the previous calving interval and days dry were excluded because these values were just the same (zero) since cows had calved for the first time. The random effects of the permanent environmental effects were also excluded, since there were no repeat records within parity. For data set 3, the permanent environmental effect was excluded but the fixed effects of the previous calving interval and days dry were included in the mixed animal model.

5. Results and Discussion

Variance components, heritability, and repeatability estimates are given in Table 4.

Table 4: Variance components and genetic parameters (±s.e) for production and Log₁₀SCC traits by parity.

Parity
Trait

All
Milk (kg)
86263
25532
173158

0.39

Fat (kg)
428
78
822

0.38

Protein (kg)
201
36
371

0.39

Fat%
0.0587
0.0095
0.0710

0.49

Protein%
0.0105
0.0016
0.0117

0.51

Log₁₀SCC
0.0252
0.0224
0.2500

0.16

1
Milk (kg)
106632
—
172145

—
—

Fat (kg)
482
—
750

—
—

Protein (kg)
247
—
386

—
—

Fat%
0.0364
—
0.0405

—
—

Protein%
0.0122
—
0.0126

—
—

Log₁₀SCC
0.0357
—
0.2844

—
—

2
Milk (kg)
73587
—
205573

—
—

Fat (kg)
533
—
1122

—
—

Protein (kg)
207
—
399

—
—

Fat%
0.0614
—
0.0910

—
—

Protein%
0.0087
—
0.0122

—
—

Log₁₀SCC
0.0473
—
0.3039

—
—

3+
Milk (kg)
61399
18852
206258

0.28

Fat (kg)
245
77
779

0.29

Protein (kg)
139
53
409

0.32

Fat%
0.0282
0.0081
0.0515

0.41

Protein%
0.0081
0.0023
0.0136

0.43

Log₁₀SCC
0.0420
0.0350
0.1952

0.28

: additive genetic variance, : permanent environmental variance, : residual variance, : heritability, : permanent environment, : repeatability, and SCC: somatic cell count.

Estimates of the additive genetic variances for yield traits in multiparous cows were dramatically smaller than those estimates from first lactation cows. As a result, heritability estimates for yield traits with repeated records were lower. The reason behind this is that in multiparous cows a repeatability animal model was used. One of the assumptions behind the use of repeatability models in the present study is that all lactations are genetically the same trait. This implies that the same set of genes exerts a common influence on the phenotypic expression of the first and later lactations [2]. This assumption is supported by observations that the genetic correlations among records are close to unity [10]. Conversely, changes in the additive genetic variances for the percentage traits in all records were small compared to those from the first lactation cows. Apparently, Jersey cows are able to systematically maintain milk composition throughout their lives [2].

The heritability and repeatability estimates for yield traits were somewhat lower than the average estimates of such estimates given by [11, 12], but heritability and repeatability estimates for protein and protein percentage are reported here for the first time for this breed. Banga [11] and Makuza et al. [12] reported heritability and repeatability estimates of 0.54, 0.56 and 0.36, 0.38 for milk yield in Zimbabwean Jersey cows, respectively. This is a little too high possibly due to method of estimation which was Harvey’s least square mixed models by [11] and due to the small sample sizes by [12]. However, the results obtained in this study were slightly higher than those reported by [2] working with Jersey cattle population in a subtropical environment. Roman et al. [2] reported heritability estimates of 0.18, 0.14, and 0.10 for milk, fat, and protein yields, respectively. Differences in the estimates are expected as a result of differences in populations, in estimation methods, and in the mathematical models employed and because of sampling errors [13]. Small heritabilities (less than 0.10) mean that regardless of genetic evaluations and selection methods used, the genetic gain would be relatively small [12].

The heritability estimates for all records for SCC obtained in this study were not significantly different from zero. However, heritability estimates for SCC were lower in the first lactation and increases the later lactations, but they were also not significantly different from zero. This suggests that little genetic progress can be made by selection of this trait. High values of SCC are highly associated with mastitis [14]; therefore, the low heritability value of SCC implies that substantial improvement has been made of Jersey herds in Zimbabwe from mastitis. As a result, this has led to good quality milk.

The heritability estimates for all production traits are highest during the first lactation and decreases towards the later lactations. This is in agreement with the results reported by various researchers [2, 15, 16]. This result suggests a sizeable genetic maternal effect in the first lactation which drops with the second lactation and disappears by the third lactation [17].

The phenotypic and genetic correlations of 305-day production and SCC variables are given in Table 5. The genetic correlations of milk with yields of protein and fat and of protein yield with protein were very high, especially for milk yield with protein yield. Phenotypic correlations were higher than the corresponding genetic correlations. This means that an increase in protein yield can be achieved through selection of milk or fat yields. The correlation between protein yield and protein percentage in this study indicates that it was negative, but in general, this correlation is positive [18]. The cause for this is the high positive correlation between protein yield and milk yield and the high negative correlation between protein percentage and milk yield. Similar trend was also observed in the first lactation and later lactations. However, the correlations were generally high in the first two lactations than in later lactations and all records. The reason behind this is the low heritability estimates in later lactations and all records [19].

Table 5: Genetic (below diagonal) and phenotypic (above diagonal) correlations and their standard errors between production traits and Log₁₀SCC across and within lactations.

Parity
Trait
Milk (kg)
Fat (kg)
Protein (kg)
Fat%
Protein%
Log₁₀SCC

Milk (kg)

0.93 (0.00)
0.98 (0.00)
0.62 (0.21)
−0.78 (0.02)
−0.04 (0.03)

Fat (kg)
0.91 (0.00)

0.93 (0.00)
0.86 (0.01)
−0.78 (0.02)
−0.01 (0.03)

All
Protein (kg)
0.95 (0.00)
0.91 (0.00)

0.64 (0.21)
−0.88 (0.01)
−0.03 (0.03)

Fat%
−0.60 (0.21)
0.83 (0.01)
−0.61 (0.21)

0.59 (0.30)
−0.04 (0.03)

Protein%
−0.75 (0.02)
−0.74 (0.02)
−0.86 (0.01)
0.55 (0.30)

−0.01 (0.03)

Log₁₀SCC
−0.03 (0.03)
−0.03 (0.03)
−0.01 (0.03)
−0.03 (0.03)
−0.01 (0.03)

Milk (kg)

0.98 (0.00)
0.98 (0.00)
0.95 (0.00)
−0.89 (0.01)
−0.01 (0.04)

Fat (kg)
0.96 (0.00)

0.98 (0.00)
0.98 (0.00)
−0.94 (0.00)
−0.01 (0.04)

1
Protein (kg)
0.96 (0.00)
0.92 (0.00)

0.98 (0.00)
−0.95 (0.00)
−0.01 (0.04)

Fat%
−0.94 (0.00)
0.95 (0.00)
−0.94 (0.00)

0.98 (0.00)
−0.02 (0.04)

Protein%
−0.85 (0.01)
−0.89 (0.01)
−0.92 (0.00)
0.97 (0.00)

−0.01 (0.04)

Log₁₀SCC
−0.38 (0.32)
−0.59 (0.31)
−0.42 (0.32)
−0.63 (0.37)
−0.44 (0.27)

Milk (kg)

0.98 (0.00)
0.98 (0.00)
0.85 (0.01)
−0.89 (0.01)
−0.01 (0.04)

Fat (kg)
0.92 (0.00)

0.97 (0.00)
0.93 (0.00)
−0.87 (0.01)
−0.01 (0.04)

2
Protein (kg)
0.92 (0.00)
0.97 (0.00)

0.83 (0.01)
−0.94 (0.00)
−0.01 (0.04)

Fat%
−0.85 (0.01)
0.89 (0.01)
−0.76 (0.02)

0.75 (0.02)
−0.05 (0.07)

Protein%
−0.87 (0.01)
−0.89 (0.01)
−0.88 (0.01)
0.72 (0.02)

−0.01 (0.04)

Log₁₀SCC
−0.49 (0.36)
−0.63 (0.42)
−0.57 (0.32)
−0.69 (0.37)
−0.56 (0.32)

Milk (kg)

0.95 (0.00)
0.98 (0.00)
0.62 (0.21)
−0.75 (0.02)
−0.05 (0.07)

Fat (kg)
0.92 (0.00)

0.97 (0.00)
0.82 (0.01)
−0.81 (0.01)
−0.04 (0.07)

3+
Protein (kg)
0.96 (0.00)
0.92 (0.00)

0.68 (0.21)
−0.85 (0.01)
−0.04 (0.07)

Fat%
−0.59 (0.32)
0.78 (0.02)
−0.63 (0.21)

0.73 (0.02)
−0.01 (0.04)

Protein%
−0.72 (0.02)
−0.78 (0.02)
−0.81 (0.01)
0.68 (0.21)

−0.01 (0.04)

Log₁₀SCC
−0.19 (0.04)
−0.18 (0.04)
−0.18 (0.04)
−0.14 (0.04)
−0.13 (0.04)

This study also shows that the genetic and phenotypic correlations between production traits and SCC were weak almost nonsignificantly, different from zero. Juozaitiene et al. [20] and Němcová et al. [21] also reported very weak genetic and phenotypic correlations, while [22, 23] reported absence of genetic correlation between production traits and SCC.

6. Conclusion

Variance due to permanent environmental effects was an important source of variation. The milk traits for Jersey dairy cattle in Zimbabwe were found to be more heritable than the SCC, and the correlation among the milk traits was significantly high. The correlation between the production traits and SCC was weak almost nonsignificantly different from zero. Therefore, these production traits may be preferred by breeders as selection criteria for development of effective genetic improvement programme.

Acknowledgments

The authors thank Mr. N. B. Mware and his colleagues of Zimbabwe Livestock Identification Trust (LIT) for providing them with data that was used in this study.

References

V. E. Imbayarwo-Chikosi, Dairy Cattle Genetics and Breeding Module, Faculties of Agriculture and Veterinary, University of Zimbabwe, Harare, Zimbabwe, 2010.

R. M. Roman, C. J. Wilcox, and F. G. Martin, “Estimates of repeatability and heritability of productive and reproductive traits in a herd of Jersey cattle,” Genetics and Molecular Biology, vol. 23, no. 1, pp. 113–119, 2000.

N. B. Mware, “Dairy herd improvement annual herd averages 1999 to 2010 report,” Tech. Rep., Livestock Identification Trust, Harare, Zimbabwe, 2010.

K. Kunaka and S. M. Makuza, “Genetic and environmental trends for milk traits for Zimbabwean Holstein-Friesian population,” Pakistan Journal of Biological Sciences, vol. 8, no. 7, pp. 1011–11015, 2005.

E. Missanjo, V. Imbayarwo-Chikosi, and T. Halimani, “A proposed selection index for Jersey Cattle in Zimbabwe,” ISRN Veterinary Science, vol. 2013, Article ID 148030, 3 pages, 2013.

USDA, Agro-Climatic Zones in Zimbabwe, United States Department of Agriculture, 2009.

SAS 9.1.3 Qualification ToolS User's Guide, SAS Institute, Cary, NC, USA, 2004.

B. Kinghorn and S. Kinghorn, Pedigree Viewer, University of New England, Biddeford, Me, USA, 2010.

A. R. Gilmour, B. R. Cullis, S. J. Welham, and R. Thompson, ASReml Reference Manual, NSW Agriculture Biometrical Bulletin 3, 2002.

T. Druet, F. Jaffrézic, and V. Ducrocq, “Estimation of genetic parameters for test day records of dairy traits in the first three lactations,” Genetics Selection Evolution, vol. 37, no. 3, pp. 257–271, 2005.

C. Banga, “Genetic parameters for milk production traits in Jersey cattle,” Zimbabwe Journal of Agricultural Research, vol. 30, no. 1, pp. 45–48, 1992.

S. M. Makuza, V. Muchenje, and S. Chiyanike, “Genetic evaluation of grade, appendix and pedigree cow classes in Holstein, Jersey and crossbred dairy breeds in Zimbabwe,” in Proceedings of the Regional Conference, Lilongwe, Malawi, 2000.

A. Sattar and R. H. Mirza, “Haematological parameters in exotic cows during gestation and lactation under subtropical conditions,” Pakistan Veterinary Journal, vol. 29, no. 3, pp. 129–132, 2009.

E. M. Missanjo, V. E. Imbayarwo-Chikosi, and T. E. Halimani, “Genetic trends production and somatic cell count for Jersey cattle in Zimbabwe born 1994 to 2005,” Tropical Animal Health and Production, vol. 44, pp. 1921–11925, 2012.

J. M. Abdallah and B. T. McDaniel, “Genetic parameters and trends of milk, fat, days open, and body weight after calving in North Carolina experimental herds,” Journal of Dairy Science, vol. 83, no. 6, pp. 1364–1370, 2000.

I. Kaya, Y. Akbaş, and C. Uzmay, “Estimation of breeding values for dairy cattle using test-day milk yields,” Turkish Journal of Veterinary and Animal Sciences, vol. 27, no. 2, pp. 459–464, 2003.

A. Hashemi, I. Bernousi, S. R. Zadeh, and M. Ranjbari, “Estimates of genetic parameters of productive traits in Holstein-native crossbreds in W. Azarbaijan Province-Iran,” Journal of Animal and Veterinary Advances, vol. 8, no. 5, pp. 917–919, 2009.

V. Pantelić, M. M. Petrović, S. Aleksić, L. Sretenović, D. Ostojić-Andrić, and Ž. Novaković, “Investigation of the genetic correlation between milk and fertility traits of first calving cows of Simmental breed,” Biotechnology in Animal Husbandry, vol. 24, no. 5-6, pp. 1–18, 2008.

M. L. Makgahlela, B. E. Mostert, and C. B. Banga, “Genetic relationships between calving interval and linear type traits in South African Holstein and Jersey cattle,” South African Journal of Animal Sciences, vol. 39, no. 1, pp. 90–92, 2009.

V. Juozaitiene, A. Juozaitis, and R. Micikeviciene, “Relationship between somatic cell count and milk production or morphological traits of udder in Black-and-White cows,” Turkish Journal of Veterinary and Animal Sciences, vol. 30, no. 1, pp. 47–51, 2006.

E. Němcová, M. Štípková, L. Zavadilová, J. Bouška, and M. Vacek, “The relationship between somatic cell count, milk production and six linearly scored type traits in Holstein cows,” Czech Journal of Animal Science, vol. 52, no. 12, pp. 437–446, 2007.

T. Ikonen, S. Morri, A.-M. Tyrisevä, O. Ruottinen, and M. Ojala, “Genetic and phenotypic correlations between milk coagulation properties, milk production traits, somatic cell count, casein content, and pH of milk,” Journal of Dairy Science, vol. 87, no. 2, pp. 458–467, 2004.

M. Koivula, E. A. Mäntysaari, E. Negussie, and T. Serenius, “Genetic and phenotypic relationships among milk yield and somatic cell count before and after clinical mastitis,” Journal of Dairy Science, vol. 88, no. 2, pp. 827–833, 2005.