Estimation of genetic parameters of the productive and reproductive traits in Ethiopian Holstein using multi-trait models.
During the last decades breeding program of dairy herds in Ethiopia had focused on meeting the ever-increasing demand for milk and milk products of the ever increasing population growth in Ethiopia and also contribute to poverty alleviation. Genetic improvement of the indigenous cattle had prioritised on crossbreeding backed by within country selection program on Holstein Friesian herds at Holetta bull dam farm and other follower herds to ensure young bull replacement for artificial insemination services in the country. But, a breeding plan and implementation of selection program lacks detailed genetic parameters analysis on important milk production and reproduction traits like lactation milk yield (LMY), lactation length (LL), age at first calving (AFC), days open (DO), and calving intervals (CI) etc. Therefore, an accurate estimation of genetic parameters in tropical herds is important for planning and implementing efficient breeding programs .
Genetic parameter estimates of locally available exotic and crossbred dairy cattle in Ethiopia have so far been carried out using univariate models. Evaluations based on single trait models does not account for the covariances among traits. Therefore, using estimates from univariate models that do not account for the relationship among traits may result in inaccurate estimates of breeding values that will leads to selection bias. To date multivariate analyses are generally used to select the best replacment stock so that the effect of selection for milk production would be known to affect reproduction traits in the anticipated direction. In additon, using multivariate analyses has a great importance in providing reliable and unbiased esimates of geneic parameters . Variance covariance estimates from a multivariate data analysis in the dairy herd helps to assess the magnitude of genetic correlations among breeding goal traits that enables the setting up of total merit index for the accurate evaluation of the genetic merit of animals in the herd. In general, information on the estimates of genetic parameters of productive and reproductive traits of Holstein Friesian dairy cattle particularly those with multivariate animal models are scarce in tropical environments. Therefore, the objective of this study was to estimate genetic parameters for productive and reproductive traits fitting multi- trait model that would contribute to the design of the national dairy genetic improvement programs in tropical dairy production systems.
MATERIALS AND METHODS
Animals and phenotypes
Data from 1,567 animals of Holstein Friesian born to eight sire origins and 136 sires, 1,431 cows that calved between 1979 and 2013 were used. The data was edited in such a way that incomplete lactation record on lactation milk yield less than 1,000 kg and lactation lengths less than 220 days and greater than 450 days were excluded to avoid bias Ayied et al . In addition, all lactation length records between 21 and 320 days were standardized to 305-days milk yield using the projection procedures of Rege . Lactations greater than 321 days were adjusted by a regression model that considered coefficient of determination and variance inflation factor for each independent variable in the model as described by Khan et al . All animals with AFC between 18 to 60 months were considered as normal expected age range and those animals with records greater than 60 months age for AFC were considered in the terminal category, similarly calving interval between 300 and 900 days and days open between 21 and 500 days were included in the final analysis . Finally a total of 3,552, 3,733, 2,938, 1,125,2,764, and 2,773 records on LMY, 305 days milk yield (305-d MY), LL, AFC, CI, and DO, respectively were used in the final analysis of the mixed models fitted to generate the estimates (Table 1).
Seasons of calving was classified into three groups based on weather and climatic conditions of the area. These included June to September as long rainy season, March to May as short rainy season and October to February as dry season . Further, lactations were classified into three parities as 1, 2, and 3+. Party three and above were all pooled together as parity three due to very few number of observations in later lactations as well as due to the higher correlation between 3rd and later lactations. Sires used in the farm were assigned into eight groups based on their source or country of origin. i.e. those imported from Cuba, Finland, Kenya, Israel, Italy, United States of America, and the remaining recruited in Holetta and unknown group from Ethiopia.
Estimation of variance and (co)variance components were made using the multivariate mixed model program package (DMU) using average information restricted maximum likelihood algorithm . To estimate variance and (co)variance components and correlations, single and multi-trait animal models were used for the different traits depending on the nature of the data. A full multi-trait animal model for six traits was initially tried, but it had convergence problems. Other multiple-trait analyses were tried by reducing the number of traits step wisely one at a time from 6 but only the 2 trait and 3 trait multi-trait models resulted in logical estimates. Preliminary analysis was performed to identify the fixed effects which should be considered in the genetic analysis. The presence of any significant differences was checked by using Tukey Kramer multiple comparison tests. Fixed effects which were significant (p<0.05) were fitted in to the model to estimate the genetic parameters. Two different models, model 1 and 2 as shown below were used for estimation of variance components. In general, the fixed effects included in these models were season of calving, year of calving/birth, origin of sire and parity whilst permanent environment, animal genetic and residuals as random effects. The description of the models in matrix notation was:
Model 1 is used for univariate analysis for AFC
Y = Xb + Za + e
Where, Y is the vector of observation; b includes vector of fixed effects; a is the vector of solutions for the coefficients of direct animal (additive) genetic random effects; e is the vector of residual effects; and X and Z are incidence matrices of fine fired effects and additive genetic random effects, respectively.
Model 2 was used for univariate and multivariate (two trait and three trait) analysis of productive and reproductive traits (LMY, LL, 305-d MY, Cl, and D O) fitting permanent environmental effort due to repeated records per cow.
Y = Xb + Za + Wpe + e
Where, Y is the vector of observations; b includes vector of fixed effects; a is the vector of solutions for the coefficients of direct animal (additive) genetic random effects; pe is the vector of solution for permanent environmental effects; e is the vector of residual effects; and X, Z, W, are the correspondent incidence matrices of the fixed effects and additive genetic and permanent environmental random effects, respectively.
In the two trait analysis we considered 305-d MY as a major trait to evaluate the milk yield of dairy cows as in most modern dairy farms, a lactation length of 305 days commonly accepted as a standard. This standard allows for ideal calving every 12 months with a 60-day dry period. Therefore, we fitted lactation length and reproductive traits with 305-d MY and presented the heritability estimates of all traits in the two trait model. We expected a different h2 estimate for 305-d MY. However, we found that all the estimates were similar and we just picked the estimates from one of the runs. Likewise, in the three trait model analysis we considered 305-d MY as a major trait and AFC as second trait based on its [h.sup.2] estimate (genetic variability) from the univariate analysis and economic value of the traits as indicated by Emana . So we have made a combined analysis of 305-d MY and AFC with LL, CI, and DO and we got h2 estimates for all traits just like two trait analysis. In both two and three trait we used LMY instead of 305-d MY merely to test its association with production and fertility traits.
Estimation of genetic parameters
Heritability was estimated as the ratio of the additive genetic variance to total phenotypic variance; and repeatability, as the ratio of the sum of the additive genetic variance and permanent environmental variance to phenotypic variance .
[h.sup.2] = [[sigma].sup.2.sub.a]/[[sigma].sup.2.sub.a] + [[sigma].sup.2.sub.pe] + [[sigma].sup.2.sub.e] r = [[sigma].sup.2.sub.a] + [[sigma].sup.2.sub.pe]/[[sigma].sup.2.sub.a] + [[sigma].sup.2.sub.pe] + [[sigma].sup.2.sub.e]
Where, [h.sup.2] is heritability; r is repeatability; [[sigma].sup.2.sub.a] is direct additive genetic variance; [[sigma].sup.2.sub.pe] is permanent environmental variance related to repeated racerds and [[sigma].sup.2.sub.e] is residual variance.
The genetic correlations were estimated as the ratio of additive covariance of the two traits (trait i and j) to the products of additive genetic standard deviations of the two traits whilst phenotypic correlations were estimated similarly as ratio of the sum of the genetic and environmental covariances (phenotypic covariances) to the products of phenotypic standard deviations of the two traits.
[r.sub.g] = [[sigma].sub.aij]/[square root of [[sigma].sup.2.sub.ai][[sigma].sup.2.sub.aj]] [r.sub.p] = [[sigma].sub.pij]/[square root of [[sigma].sup.2.sub.pi][[sigma].sup.2.sub.pj]]
Where, [r.sub.g] is genetic correlations; [r.sub.p] is phenotypic correlation; [[sigma].sup.2.sub.ai] is additive genetic variance for trait i; [[sigma].sup.2.sub.aj] is additive genetic variance for trait j; [[sigma].sup.2.sub.pi] is phenotypic variance for trait i, [[sigma].sup.2.sub.pj] is phenotypic variance for trait j; [[sigma].sub.aij] is additive genetic covariance between trait i and j; and [[sigma].sub.pij] is phenotypic covariance between trait i and j.
RESULTS AND DISCUSSION
The heritability estimates for LMY in the present study from single and three trait models were similar (0.17 [+ or -] 0.04) whilst in the two traitmodelwas0.25 [+ or -] 0.02 (Table2). These values of heritability for LMY indicated presence of modest genetic variability among the small population in Holetta. The increment in the estimates in two trait analyses for LMY and 305-d MY record of an animal under the study indicated the presence of high correlation among the traits that amplified the heritability estimates. The current estimate from the univariate model is comparable to those estimates reported for Holstein Friesian dairy cattle [11-13]. On the other hand, it is higher than those reported for first lactation Holstein Friesians in Egypt .
The heritability estimates for adjusted 305-d MY from univariate, 2 trait and 3 trait analysis were 0.15 [+ or -] 0.04, 0.16 [+ or -] 0.04, and 0.15 [+ or -] 0.04, respectively (Table 2). These values of heritability for 305-d MY indicate presence of modest genetic variability among the small population in the study. However, unlike the estimates for LMY the heritability estimates in univraite, two and three traits models did not show much change, this could be due to adjustment made on lactation days that triggers weak correlation to other traits. The current estimates from all models were comparable with 0.11 for Egyptian Holstein Friesian dairy cattle . Similarly, closely heritability estimate of 0.19 reported for the same breed in a temperate environment . However, higher estimates were reported by several other authors from different production environments [4,16,17]. Such differences in the heritability estimates are expected and are mainly results of size of dataset, genetic variation with in population, management and environmental conditions and the methods used for parameter estimation. In addition, we believed that a long term biased phenotypic selection and involuntary culling of animals from the existing herd might have contributed to the decreased hereditary differences and genetic variability within the population. Although estimate of heritability for LMY and adjusted 305-d MY indicates that still there is an opportunity for improving these traits through selection, future improvement of the herd could be meaningful if the current phenotypic selection is changed to selection based on genetic merit combined with improved herd management as well as introducing of new animals from other source to increase genetic variability within the herd.
The heritability estimates for lactation length ranged from 0.03 [+ or -] 0.03 to 0.08 [+ or -] 0.03 (Table 2) and these estimates correspond to those reported by Ojango and Pollott . It is slightly higher than heritability estimates of 0.003 [+ or -] 0.08 in Sudan using paternal half-sib correlation . These low heritability estimates obtained in this study could be attributable to the high residual variance. In the study the heritability estimates of LL shows slight improvement from univariate model to 3 traits model. This enhanced accuracy is due to the ability of the model to account for relationship between the traits and better connections in the data due to residual covariance between traits . Moreover, multivariate models gave better reliabilities than univariate models in the prediction of genetic merit which is mainly a result of the model's ability to use extra information from correlated traits .
Heritability estimates for AFC were 0.47, 0.38, and 0.38 from univariate, 2 trait and 3 trait models, respectively (Table 2). These high heritability estimates for AFC from this study may be explained by weak selection made on this trait in contrast to the intense selection practiced for milk production and presence of genetic variability in the small herd as the result of the different sources of sires being used in the study period. The current estimates were comparable with estimates of Rege  and Ojango and Pollott  who reported 0.47 and 0.38, respectively for the same breed in Kenya. Differences between the estimates of heritability obtained in the study and estimates from other countries are most likely caused by management and climate differences affecting genetic and environmental variances and methods of estimation . Our study suggests that improvement of this trait through selection in line with several other authors who have come to the same conclusion [4,22].
Heritability estimates for calving interval were 0.11 [+ or -] 0.04, 0.11 [+ or -] 0.04, and 0.13 [+ or -] 0.04 from univariate, 2 trait and 3 trait models, respectively (Table 2). These estimates for CI is still within the literature estimates for the traits, indicating the existence of certain level of genetic diversity. The present result is lower than estimates of 0.17 for Iraqi Friesian cows . However, it is higher than the estimates reported for the same breed in other tropical countries [17,23,24]. This relatively low heritability estimates for CI could be explained by large environmental variance. Therefore, improvements in nutrition and reproductive management should lead to a considerable decrease in length of CI than making purely genetic selection alone. However, in the study, multi-trait analysis with milk production trait and AFC indicates slight improvement in the estimates of heritability of calving interval. Incorporating information on correlated traits may enhance the accuracy of selection and reduce selection bias made by univariate analysis . Therefore, a multiple trait model selection including milk production with fertility traits could improve the selection accuracy for selecting the right animals with good fertility traits [15,25].
Heritability estimates for days open in the same study were 0.09 [+ or -] 0.03, 0.09 [+ or -] 0.03, and 0.10 [+ or -] 0.04 from univariate, 2 trait and 3 trait models, respectively (Table 2). The heritability estimates from all models were similar and coincide with most literatures. The heritability estimates for DO is also indicative of the presence of genetic basis that deserve selection or culling of an animal. The current estimates were close to 0.082 reported by Zambrano and Echeverri , but higher than 0.042 that was reported in temperate zone . On the contrary, higher estimates of days open have been reported in Egypt and Iraq for the same breed [26,13].
In the current study the repeatability estimate for LMY and adjusted 305-d MY is moderate (Table 2). For LMY the repeatability estimates were 0.39 and 0.45 form univariate and 2 trait models, respectively, which are comparable with that reported by Ojango and Pollott , but higher than the reported value of 0.26 [+ or -] 0.08 for the same breed in Pakistan . As for adjusted 305-d MY the repeatability estimate were 0.42 and 0.41, these values were in close agreement with 0.43 reported by Dematawewa and Berger . Nevertheless, Rege  and Kadarmideen et al  have reported higher repeatability estimates of 0.55 and 0.58, respectively. These moderate repeatability estimates in the present study indicated that selection of dairy cattle for milk yield based on early performance records is reliable.
Regarding lactation length the repeatability estimate was very low which correspond with Ojango and Pollott . However, slightly higher estimate of 0.19 was reported for the same breed . The low repeatability estimate implied that variation in lactation length was more a result of variation in management and feeding in the given environment of a particular lactation rather than of factors associated genetic makeup of the cow [17,12].
The repeatability estimates for calving interval and days open from all models were very low (Table 2). Estimates for CI and DO were comparable with Rege  who reported 0.25 for the same breed in Kenya, but higher than estimates reported by Zambrano and Echeverri . The low repeatability estimates obtained in this study suggest that CI and DO are strongly influenced by temporary environmental factors, this is mainly because of the complex nature of reproductive traits, difficulties in detection of oestrus, and various other managerial and nutritional factors.
Genetic correlations between milk production traits were high. The highest genetic correlation was 0.99 [+ or -] 0.01 between LMY and 305-d MY and the lowest was 0.73 [+ or -] 0.19 between 305-d MY and LL (Table 3). This is comparable with Ahmad et al . The highest genetic correlation between milk productions traits were result from pleiotropy . This means that genetic improvement of one milk production trait could result in a correlated response in the correlated trait. Generally the very high correlation between these two traits is indicative of the fact that using 305-d MY can be sufficient to meet selection for lactation milk yield.
The genetic correlation between 305-d MY and AFC was -0.24 [+ or -] 0.11 (Table 3). This is in close agreement with -0.22 [+ or -] 0.08 reported in Ethiopia . On the contrary, a relatively higher and positive (0.54) genetic correlation was reported for Holstein Friesian dairy cattle in Kenya . The negative genetic correlation between 305-d MY and AFC indicates that genes which affect the 305-d MY trait positively will also be responsible for an early onset of puberty among the current population found at Holetta bull dam herd.
Similarly genetic correlation between 305-d MY and CI found to be -0.10 [+ or -] 0.15. This result is accords with-0.03 correlation estimate of Friesian cattle in Kenya . However, higher and positive genetic correlations of 0.54 [+ or -] 0.01 was reported for the same breed in Ethiopia . Regarding genetic correlation between 305-d MY and DO, it was close to zero (-0.02 [+ or -] 0.16). Positive and lower estimates were reported in Turkey  and Egypt . On the other hand, higher and positive estimates of genetic correlation was observed for Ethiopian Holstein . The negative genetic correlation of 305-d MY and fertility traits indicates that improvement of fertility is possible to certain degree in this herd as milk yield increase. However, since the standard error of the estimates were high, the results needs to be interpreted with caution. The higher standard errors obtained from pair wise combinations of CI and DO with 305-d MY could be due to the smaller proportion of CI and DO records in the dataset.
In the current study, positive genetic correlations between various reproductive traits were observed. This positive genetic correlation between reproductive traits indicates that improvement of one trait has positive impact on the other trait. The genetic correlation between CI and DO was 0.99 [+ or -] 0.01 (Table 3). This is in close agreement to a nearly perfect genetic correlation of 0.99  and 0.98 . The positive and strong genetic correlation in the present study indicate the evidence of common genetic and physiological mechanisms controlling those traits. In the tropics in general and in the Holetta bull dam herd in particular the recording system is still at its early stage of development. As a result, recording of calving interval is relatively easier than days open. Therefore, evaluation based on calving interval would be easy to handle by the available human resource and recording infrastructure than evaluation based on days open.
In this study phenotypic correlations between milk productions traits were high (Table 3). The highest phenotypic correlation was 0.84 [+ or -] 0.01 between 305-d MY and LMY and the lowest was 0.05 [+ or -] 0.02 between 305-d MY and LL. This explains that these traits are influenced by more or less same genetic and environmental factors. The highest phenotypic correlation among LMY and 305-d MY was in agreement with the estimate reported for crossbreeds dairy cattle in Egypt [26,27].
Phenotypic correlations between 305-d MY with fertility traits (AFC, CI, and DO) were in general negative (Table 3). The phenotypic correlation between 305-d MY and AFC in present study was higher than Ahmad et al  for crossbreed cattle (-0.02) and -0.02 for first parity Holstein Friesian . Conversely, positive phenotypic correlation was reportedfor Jersey cattle in Ethiopia . In the current study the phenotypic correlation between CI and DO is 0.98 [+ or -] 0.01 (Table 3). This strong and positive phenotypic correlations is in agreement with results of Kadarmideen et al , Ghiasi et al , and Zambrano and Echeverri  who reported 0.95, 0.95, and 1.0 for Holstein Friesian cattle in different production environment, respectively.
Estimates of genetic parameters are essential strategic step in launching a breeding programs which could be used as a tool for implementation of selective breed improvement in a population. Therefore, the correct estimation of genetic parameters that guarantee for accurate prediction of genetic merit of an individual is essential. The heritability estimates for all production and reproduction traits indicated that the major part of variation for those traits were due to environmental factors. Although the heritability estimate for LMY and 305-d MY were relatively low, there is still some room for genetic improvement in the herd. The repeatability estimate for lactation length and all reproductive traits were low indicating that these traits are mainly influenced by managerial and temporary environmental effects. The genetic correlations observed for milk production and reproductive traits ranged from moderate to high. Multi-trait analysis also showed slight improvement in heritability estimates of CI and LL than univariate model analysis. This indicates that selection with multitrait model could be better than univariate model due to the correction of selection bias by multivariate model as a result of using the information from the correlated traits. Therefore, in the future verification of some lower estimates of certain traits including large set of data and fitting multivariate model for both productive and reproductive traits is recommended.
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
The authors are grateful for the technical and financial support by Jimma University, Ethiopia and the Natural Resource Institute of Finland (Luke). The National Artificial Insemination Centre (NAIC), Ethiopia and particularly the management Holetta dairy farm are duly acknowledged for allowing us to use their farm record.
[1.] Choudhary V, Kothekar MD, Raheja KL, et al. Genetic evaluation of first lactation traits in Sahiwal cattle using restricted maximum likelihood technique. Asian-Australas J Anim Sci 2003;16:639-43.
[2.] Unalan A, Cebeci Z. Estimation of genetic parameters and correlations for the first three lactation milk yields in Holstein Friesian cattle by the REML method. Turk J Vet Anim Sci 2004;28:1043-9.
[3.] Ayied YA, Jadoa JA, Abdulrada JA. Heritabilties and breeding values of production and reproduction traits of Holstein cattle in Iraq. J Basrah Res Sci 2011;37:66-70.
[4.] Rege JEO. Genetic analysis of reproductive and productive performance of Friesian cattle in Kenya. J Anim Breed Genet 1991;108: 412-23.
[5.] Khan MS, Rehman ZU, Ahmad S, Hassan FU, Rehman MSU. Prediction of standard lactation milk yield from completed lactations of longer duration in cattle and buffaloes. Pakistan Vet J 2011;32:122-4.
[6.] Ajili N, Rekik A, Gara B, Bouraoui R. Relationships among milk production, reproductive traits, and herd life for Tunisian Holstein Friesian cows. African J Agric Res 2007;2:047-051.
[7.] Sendros D, Neser FWC, Schoeman SJ. Estimates of genetic parameters for Boran, Friesian, and crosses of Friesian and Jersey with the Boran cattle in the tropical highlands of Ethiopia: milk production traits and cow weight. J Anim Breed Genet 2004;121:163-75.
[8.] Madsen P, Jensen J. DMU A package for analysing multivariate mixed models. Version 6, release 5.2. Tjele Denmark: Center for Quantitative Genetics and Genomics Dept. of Molecular Biology and Genetics, University of Aarhus Research Centre; 2013.
[9.] Emana G. Economic index weight age for dairy traits preference in central highlands of Ethiopia [master's thesis]. Jimma, Ethiopia: Jimma University; 2015.
[10.] Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. London, England; Longman group; 1996. p.463.
[11.] Gader AZ, Ahmed MK A, Musa LMA, Peters KJ. Milk yield and reproductive performance of Friesian cows under Sudan tropical conditions. Arch Tierz 2007;50:155-64.
[12.] Usman T, Guo G, Suhail SM, et al. Performance traits study of Holstein Friesian cattle under subtropical conditions. J Anim Plant Sci 2012; 22:92-5.
[13.] Hermiz HN, Juma KH, Saadi SK, Aldoori TS. Genetic parameters of production, reproduction and growth traits of Holstein cows. Dirasat Agric Sci 2005;32:157-62.
[14.] Hammoud MH, Salem MMI. The genetic evaluation of some first lactation traits of Holstein cows in Egypt. Alex J Agric Res 2013; 58:1-8.
[15.] Dematawewa CMB, Berger PJ. Genetic and phenotypic parameters for 305-day yield, fertility, and survival in Holsteins. J Dairy Sci 1998; 81:2700-9.
[16.] Kadarmideen HN, Thompson R, Simm G. Linear and threshold model genetic parameters for disease, fertility and milk production in dairy cattle. In UK dairy cattle. Anim Sci 2000;71:411-9.
[17.] Ojango JM, Pollott GE. Genetics of milk yield and fertility traits in Holstein Friesian cattle on large scale Kenyan farms. J Anim Sci 2001;79:1742-50.
[18.] Eid II, Elsheikh MO, Yousif IA. Estimation of genetic and non-genetic parameters of Friesian Cattle under hot climate. J Agric Sci 2012;4: 95-102.
[19.] Thompson R, Meyer K. A review of theoretical aspects in the estimation of breeding values for multi-trait selection. Livest Prod Sci 1986; 15:299-313.
[20.] Strabel T, Misztal I, Bertrand JK. Approximation of reliabilities of multiple trait models with maternal effects. J Anim Sci 2001;79: 833-9.
[21.] Montaldo HH, Castillo-Juarez H, Valencia-Posadas M, Cienfuegos-Rivas EG, Ruiz-Lopez FJ. Genetic and environmental parameters for milk production, udder health, and fertility traits in Mexican
Holstein cows. J Dairy Sci 2010;93:2168-75.
[22.] Makgahlela ML, Banga CB, Norris D, Dzama K, Ng JW Genetic correlations between female fertility and production traits in South African Holstein cattle. S Afr J Anim Sci 2007;37:180-8.
[23.] Mhamdi N, Aloulou R, Brar SK, Bouallegue M, Ben-Hamouda M. Phenotypic and genetic parameters of reproductive traits in Tunisian Holstein cows. Biotechnol Anim Husb 2010;26:297-307.
[24.] Zambrano JC, Echeverri J. Genetic and environmental variance and covariance Parameters for some Reproductive traits of Holstein and Jersey cattle in Antioquia (Colombia). R Bras Zootec 2014;43:132-9.
[25.] Kadarmideen HN, Thompson R, Coffey MP, Kossaibati MA. Genetic parameters and evaluations from single and multiple trait analysis of dairy cow fertility and milk production. Livest Prod Sci 2003;81: 183-95.
[26.] Shalaby NA, Oudah EZM, Abdel-Momin M. Genetic analysis of some productive and reproductive traits and sire evaluation imported and locally born Friesian cattle raise in Egypt. Pakistan J Boil Sci 2001;4:893-901.
[27.] Ahmad M, Van der Werf JHJ, Javed K. Genetic and phenotypic correlations for some economic traits in dairy cattle. Pakistan Vet J 2001; 21:81-6.
[28.] Yosef T. Genetic and non-genetic analysis of fertility and production traits in Holetta and Adaa Berga dairy herds [master's thesis]. Alemaya, Ethiopia: Alemaya University; 2006.
[29.] Atil H, khattab AS, Yakupoglu C. Genetic analysis for milk traits in different herds of Holstein Friesian cattle in Turkey. J Biol Sci 2001; 1:737-41.
[30.] Ghiasi H, Pakdel A, Nejati-Javaremi A, et al. Genetic variance components for female fertility in Iranian Holstein cows. Livest Sci 2011; 139:277-80.
Wondossen Ayalew (1) *, Mohammed Aliy (2), Enyew Negussie (3)
* Corresponding Author: Wondossen Ayalew Tel: +251-113220169, Fax: +251-113220041, E-mail: firstname.lastname@example.org
(1) Wolkite University College of Agriculture and Natural Resource, PO. Box 07, Wolkite, Ethiopia
(2) Jimma University College of Agriculture and Veterinary Medicine, P.O. Box, 307 Jimma, Ethiopia
(3) Natural Resources Institute Finland (Luke), Biometrical Genetics, Jokioinen 31600, Finland
Submitted Mar 13, 2017; Revised May 3, 2017; Accepted Jun 13, 2017
Table 1. Statistical description on the data set Variables N Mean SD Minimum Maximum LL 2,938 315.04 48.54 220 450 LMY 3,552 3,574.69 1,288.69 1,001 9,800 305-d MY 3,733 3,504.02 1,222.56 1,004 9,301 AFC 1,125 39.26 9.8 18 60 CI 2,764 473.57 124.32 301 897 DO 2,773 184.54 107.55 25 500 N, number of observations; SD, standard deviation; 305-d MY, 305 days milk yield; LMY, lactation milk yield; LL, lactation length; AFC, age at first calving; CI, calving interval; DO, days open. Table 2. Estimates of permanent environmental ([[sigma].sup.2.sub.pe]), additive genetic ([[sigma].sup.2.sub.a]), and residual variances (a2e), heritability (h2) and repeatability (r) of production and reproduction traits using univariate, two and tree trait models Models Traits [[sigma]. [[sigma]. [[sigma]. sup.2. sup.2. sup.2. sub.pe] sub.a] sub.e] Univariate 305-d MY 287,808.6 160,216.1 646,915.1 model LMY 294,985 229,963 830,492 LL 181.8 60.3 1,834.3 AFC - 22.7 25.9 CI 1,526.5 1,483.02 10,965 DO 1,020.0 982.6 8,571.5 Two traits 305-d MY 291,396.8 195,089.6 738,703.6 model LMY 174,216.2 273,808.6 644,915.1 LL 116.5 136.6 1,834.4 AFC - 18.2 29.7 CI 1,526.0 1,678.2 10,927.8 DO 1,006.5 1,107.0 8,558.8 Three 305-d MY 287,808.6 160,216.1 646,915.1 traits model LMY 229,781 294,639 830,914 LL 116.5 166.6 1,804.4 AFC 1,767.7 1,436.8 1,0954.7 CI - 32.5 52.3 DO 1,080.2 1,017.2 8,579.2 Models Traits [h.sup.2] r [+ or -] SE [+ or -] SE Univariate 305-d MY 0.15 [+ or -] 0.04 0.42 [+ or -] 0.02 model LMY 0. 17 [+ or -] 0.04 0.39 [+ or -] 0.02 LL 0.03 [+ or -] 0.03 0.12 [+ or -] 0.02 AFC 0.47 [+ or -] 0.06 - CI 0.11 [+ or -] 0.04 0.22 [+ or -] 0.02 DO 0.09 [+ or -] 0.03 0. 19 [+ or -] 0.02 Two traits 305-d MY 0.16 [+ or -] 0.04 0.41 [+ or -] 0.02 model LMY 0.25 [+ or -] 0.02 0.45 [+ or -] 0.02 LL 0.06 [+ or -] 0.03 0.12 [+ or -] 0.02 AFC 0.38 [+ or -] 0.27 -- CI 0.11 [+ or -] 0.04 0.23 [+ or -] 0.02 DO 0.09 [+ or -] 0.03 0.20 [+ or -] 0.02 Three 305-d MY 0.15 [+ or -] 0.04 -- traits model LMY 0.17 [+ or -] 0.04 -- LL 0.08 [+ or -] 0.03 -- AFC 0.38 [+ or -] 0.25 -- CI 0.13 [+ or -] 0.04 -- DO 0.10 [+ or -] 0.04 -- SE, standard error; 305-d MY, 305 days milk yield; LMY, lactation milk yield; LL, lactation length; AFC, age at first calving; CI, calving interval; DO, days open. Table 3. Genetic (above the diagonal) and phenotypic (below the diagonal) correlations between milk production and fertility traits using two trait analysis Traits LMY 305-d MY LMY - 0.99 [+ or -] 0.01 305-d MY 0.84 [+ or -] 0.01 -- LL 0.51 [+ or -] 0.02 0.05 [+ or -] 0.02 AFC -0.41 [+ or -] 0.11 -0.24 [+ or -] 0.11 CI 0.15 [+ or -] 0.02 -0.12 [+ or -] 0.02 DO 0.17 [+ or -] 0.02 -0.11 [+ or -] 0.02 Traits LL AFC LMY 0.91 [+ or -] 0.08 -0.41 [+ or -] 0.11 305-d MY 0.74 [+ or -] 0.19 -0.24 [+ or -] 0.11 LL -- -0.72 [+ or -] 0.38 AFC 0.08 [+ or -] 0.04 -- CI 0.39 [+ or -] 0.01 0.11 [+ or -] 0.03 DO 0.42 [+ or -] 0.02 0.10 + 0.02 Traits CI DO LMY 0.16 [+ or -] 0.12 0.25 [+ or -] 0.13 305-d MY -0.10 [+ or -] 0.13 -0.02 [+ or -] 0.14 LL 0.99 [+ or -] 0.10 0.98 [+ or -] 0.19 AFC 0.16 [+ or -] 0.12 0.06 [+ or -] 0.13 CI -- 0.99 [+ or -] 0.01 DO 0.98 [+ or -] 0.01 -- SE, standard error; 305-d MY, 305 days milk yield; LMY, lactation milk yield; LL, lactation length; AFC, age at first calving; CI, calving interval; DO, days open.
|Printer friendly Cite/link Email Feedback|
|Author:||Ayalew, Wondossen; Aliy, Mohammed; Negussie, Enyew|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Nov 1, 2017|
|Previous Article:||Tracing the breeding farm of domesticated pig using feature selection (Sms scrofa).|
|Next Article:||Effects of various weaning times on growth performance, rumen fermentation and microbial population of yellow cattle calves.|