Livestock Research for Rural Development 32 (3) 2020 LRRD Search LRRD Misssion Guide for preparation of papers LRRD Newsletter

Citation of this paper

Phenotypic and genetic parameters of milk traits of Fogera cattle in Ethiopia

Demelash Kassahun1,2, Damitie Kebede2, Aynalem Haile3, Adebabay Kebede2, Assemu Tesfa2, Mekonen Tilahun2, Addisu Bitew2, Tewodros Bimrow2, Eyasu Lakew2, Mulugeta Meseret2, Zeleke Tesema4 and Mengistie Taye2,5

1 Andassa Livestock Research Center, P.O. Box 27, Bahir Dar, Ethiopia
2 College of Agriculture and Environmental Sciences, Bahir Dar University, P O Box 5501, Bahir Dar, Ethiopia
mengistietaye@gmail.com
3 International Center for Agricultural Research in the Dry Areas (ICARDA), P O Box 5689, Addis Ababa, Ethiopia
4 Sirinka Agricultural Research Center, P O Box 74, Sirinka, Ethiopia
5 Biotechnolgy Research Institute, Bahir Dar University, P O Box 79, Bahir Dar, Ethiopia

Abstract

This study was conducted to estimate the production performance and the genetic parameters for milk traits of Fogera cattle at Andassa Livestock Research center. Data collected on milk production traits from 2000 to 2017 from Fogera cattle herds maintained at Andassa Livestock Research Center were used. The General Linear Model (GLM) procedure of SAS was used to analyze the least-squares mean for phenotypic traits, and Wombat software was used to estimate genetic parameters and breeding values. The performance milk traits analyzed include daily milk yield (DMY), lactation milk yield (LMY), and lactation length (LL). Dam parity, calving year and season were the main fixed effects considered for performance estimation of the milk traits.

The overall means for DMY, LMY, and LL obtained were 1.98±0.40 liters, 482±0.02 liters, and 238±0.3 days. The direct heritability estimates for DMY, LMY, and LL of Fogera cows were 0.33±0.27, 0.27±0.001, and 0.20±0.23. The repeatability estimates for DMY, LMY and LL were 0.33, 0.55, and 0.48. The phenotypic correlation for DMY-LMY, DMY-LL, and LMY-LL were 0.49, 0.48, and 0.55, while the genetic correlations were 0.53, 0.46, and 0.53 The mean breeding values for DMY, LMY and LL were 0.016 liters, 18.77. liters and 4.30 days. LMY increased by 0.86 liters/year while the DMY and LL decreased by 0.0008 liters/year and 0.17 days/year. Both, heritability and repeatability estimates for milk production traits indicate the presence of moderate variation within the studied cattle population and thus selection could be used as the best tool of genetic improvement of Fogera cattle for milk production.

Keywords: estimated breeding value, genetic trend, heritability, lactation length


Introduction

With its diverse and suitable agro-ecological zones, Ethiopia has been recognized for its huge and diverse (35 cattle breeds) genetic resources ( http://dagris.ilri.cgiar.org/). Despite this, the contribution of the sector to the livelihood of farm households and the country is below what is expected due to several factors including high prevalence of disease, feed shortage, and low productive potential of local breeds. With this regard, the government of Ethiopia has initiated breed improvement programs through selective breeding.

Fogera cattle is Zenga cattle adapted to the Fogera plains around Lake Tana with characteristics of: docile use for draft, meat and milk; short stumpy pointed horns; hump ranging from thoracic to cevico-thoracic; folded to moderately large dewlap; and black-and-white or black-and-gray coat color (Rege and Tawah 1999). Fogera cattle, because it is selected for environmental adaptation, its production and productivity is low. The Andassa Livestock Research center has been working on the improvement of Fogera cattle through selective breeding.

The potential for genetic improvement of a trait largely depends on genetic variation existing in the population of the study. Understanding the potential of indigenous cattle and identification of responsible factors for milk production is essential to improve the milk production performance. Previously, the milk production performances of Fogera cattle under station management conditions was reported by Bitew et al (2010). However, information on the genetic parameter estimates of milk production traits of Fogera cattle is lacking. Therefore, the aim of this study was to estimate phenotypic, genetic parameters and breeding values for milk traits of Fogera cattle at Andassa Livestock Research Center.


Material and methods

Description of the study area and herd management

Milk production data of Fogera cattle collected from 2000 to 2017 from Andassa Livestock Research Center in Ethiopia were used. The center is located in Amhara National Regional State at 11029’ north latitude and 37029’ east longitude, at an elevation of 1730 meters above sea level. The average annual rainfall is 1150 mm, and the temperature ranges from 6.5 OC to 30 OC.

Cattle in the center were allowed to graze natural pasture for 8 hours, and during the dry season, they were provided with hay harvested from the natural pasture. Animals were watered from Andassa river. The mating system was natural, a bull being assigned to 40-50 Fogera cows based on pedigree information. Herd groups were managed based on age, sex, breed and physiological status of the animal. Calves after birth had free access to suckle their dams for 4 days to consume enough colostrum. Then calves were separated from their dams and allowed to partially suckle (two teats) at milking times until weaning age of 8 months. Cows were partially hand-milked (two teats) twice a day at about 6:00 AM in the morning and 4:00 PM in the afternoon. The total milk yield produced per day was calculated by doubling the daily milk yield.

Data source and management

Milk data (daily milk yield, lactation milk yield, and lactation length) collected and recorded in a herd book from 2000 to 2017 were used. The data were filtered and organized using Microsoft Excel. Records with no pedigree information were removed during data cleaning.

Statistical analysis

The data were analyzed using the General Linear Model (GLM) procedures of the Statistical Analysis System (SAS 2004). Dam parity, calving year, and season were the main fixed effects considered. WOMBAT software (Meyer 2012) was used to estimate the genetic parameters of heritability, repeatability, genetic and phenotypic correlation.

The statistical model used for the analysis of least squares means of milk traits was:

Yijk = µ+Yi+Sj+Pk+e ijk

Where:

Yijk = The observation on daily milk yield (DMY), lactation milk yield (LMY), and lactation length (LL)

µ = Overall mean,

Yi = fixed effect of ith calving year (2000-2017)

Sj = fixed effect of jth calving season (dry, wet)

Qk = fixed effect of kth parity (1-7)

eijk = effect of random errors associated with each observation.

The model used to estimate variance components and heritability of milk traits was as follow:

Y = Xb + Za + Wpe + e

Where, Y is the vector of observations of the traits (DMY, LMY, and LL), b is the vector of fixed effects, a is the vector of solutions for the coefficients of direct animal (additive) genetic effects, pe is the vector of solution for permanent environmental effects and e is was the vector of residual effects. X, Z, W are the correspondent incidence matrices of the fixed effects, additive genetic and permanent environmental random effects, respectively.

The heritability and repeatability were computed as follows:

h2 = σa2/ σp2 and r = σa 2 + σpe2/ σp2

Where: h2 is heritability, r is repeatability, σa2 is direct additive genetic variance, σpe2 is permanent environmental variance related to repeated records and σp2 is phenotypic variance (Falconer and Mackay 1996). The genetic trends for milk traits were estimated by averaging the estimated breeding value estimated by year of birth. The pedigree information and the descriptive statistics are presented in Table 1 and Table 2, respectively.

Table 1. Information related to pedigree, traits analyzed and sample sizes

Traits

DMY

LMY

LL

Number of animals used

282

277

256

Animals after pruning

179

209

198

Animals without recording

37

42

41

Unknown Dams

128

53

136

Unknown sire

49

139

52

Dams with records and progeny

22

31

25

Animals with an unknown sire with records

12

11

11

Animals with an unknown dam with records

91

97

95

Animals with both parents unknown

48

8

49

Progeny per sire

26

29

38

Progeny per dam

33

44

28

Animals with maternal grandsire

26

39

31

Animals with maternal grand??

4

16

9

DMY, daily milk yield; LMY, lactation milk yield; LL, lactation length



Table 2. Descriptive statistics of the data set

Traits

Units

N

Mean

SD

Min.

Max.

Range

Daily milk yield

Liters

652

1.98

0.60

0.34

5.20

4.86

Lactation milk yield

Liters

652

489

184

43.9

1298

1253

Lactation length

Days

652

243

72.79

106

539

433

SD, standard deviation


Results and discussion

The effect of non-genetic factors on milk production traits
Daily milk yield

The overall mean of daily milk yield of Fogera cattle was 1.98±0.2 liters (Table 3). This result is higher than the result (1.7±0.1 liters) reported for Boran cattle under on-station management (Haile et al 2009), 1.8±0.84 liters for Kereyu cattle (Garoma et al 2014) and 1.6±0.5 liters for highland Zebu around west Gojjam under on-farm conditions (Minuye et al 2016). On the other hand, higher result for daily milk yield was reported for Barca cattle breed (2.98±0.69 liters) by Tadesse and Dessie (2003), for Arsi (2.7 liters) and for Zebu cattle (2.8 liters) by Kiwuwa et al (1983). The mean here obtained is also lower than the reports (2.8±0.15 liters/day) for Boran cattle (Demeke et al 2004) and 2.2±0.5 liters/day for local cattle breed at Beddelle areas (Gelmessa et al 2013).

Parity had a significant effect on daily milk production of Fogera cattle ( p<0.001). This result is in line with the previous results of Bitew et al (2010) for Fogera cattle and Demeke et al (2004) for Ethiopian Boran cattle. Cows with lower parities had lower daily milk yield than those of higher parities. This might be because the size of udder and teats increased with the maturity of the cows and subsequently the milk production capacity increased with parity. However, this result contradicts the study of Bayou et al (2015) who noted that parity of cow had a non-significant effect on daily milk yield of Boran cattle.

The year of calving had an effect on daily milk yield (p<0.001). A similar result was reported by Bitew et al (2010) and Demeke et al (2004).

The overall mean lactation milk yield was 489±7.50 liters (Table 3). This figure is lower when compared with results for the other indigenous cattle breeds of Ethiopia. For instance, it was much lower than the average yield reported by Kiwuwa et al (1983) for Arsi cattle (809 liters) and for Zebu cattle (929 liters) and also lower than the report of Demeke et al (2004) for Boran cattle (529±65 liters). The present result is lower than 507±39 liters for Boran cattle (Haile et al 2009) and 672±196 litters for Zebu*HF crosses (Tadesse and Dessie 2003). The lower lactation milk yield for Fogera cattle may be due to the reason that milking persists only for eight months of lactation as the weaning age of Fogera calves is at eight months.

Parity had a significant effect on the lactation milk yield (Table 3). The finding of this study showed a consistent increase in milk yield with increasing yield up to the fifth parity and declined thereafter. This result was supported by the previous reports of Bitew et al ( 2010), Demeke et al (2004) and Bayou et al (2015). Udder development increases as the age of cows advances till maturity so that the mammary gland of cows with higher parities can produce more milk than in the case of cows of lower parities.

Lactation length

The overall mean lactation length was 243±0.3 days (Table 3). The lactation length of indigenous cattle ranged from 193±6 days for Boran (Demeke et al 2004) to 303 days for Zebu (Kiwuwa et al 1983). Lactation length of 305 days is usually accepted as a benchmark in modern dairy farms. However, the lactation length of Fogera cattle was lower than most of the other indigenous breeds due to the weaning age of Fogera calves at 8 months.

Parity had no effect on the lactation length. This result contradicts with the results of Demeke et al (2004) and Bitew et al (2010). However, cows that calved in the dry season had shorter lactation length than those calved during the wet season, which is in with the results of Bitew et al (2010) and Bayou et al (2015). This could be associated with the availability of feed resources i.e. the availability of feed is better during the wet season than dry season.

The variation of lactation length with the different years might be explained by the variation of annual rainfall, which directly or indirectly is associated with the availability of feeds.

Table 3. Least squares means and standard errors of milk yield (LSM±SE)

Sources of
variation

LMY (Liters)

LL (Days)

DMY (Liters)

N

LSM±SE

N

LSM±SE

N

LSM±SE

Overall

654

489 ±7.50

654

243±2.97

654

1.98±0.02

CV

32.29

25.97

28.22

Parity

617

***

617

NS

609

***

1

137

431±16.5c

137

242±8.03

145

1.83±0.04c

2

146

489±14.9abc

146

248±6.38

130

1.94±0.04bc

3

130

521±15.0a

130

229±5.10

130

2.09±0.06ab

4

90

521±20.3ab

90

244±7.40

90

2.11±0.06ab

5

70

552±22.3a

70

243±7.26

70

2.21±0.06a

6

32

511±30.8bc

32

242±9.70

32

2.02±0.08ab

7

12

442±35.4c

12

212±2.20

12

2.13±0.12a

Season

625

NS

**

NS

Dry

213

483±13.1

213

228±5.31

213

2.06±0.04

Wet

412

508±9.20

399

246±3.56

399

2.04±0.02

Year

654

***

654

***

654

***

2003

13

834±65.9a

13

348±19.4a

13

2.42±0.15a

2004

91

610±20.4b

91

298±12.0b

91

1.79±0.05cd

2005

33

550±38.0bcd

33

293±15.6b

33

2.13±0.08ab

2006

19

571±57.9bc

19

294±22.5b

19

1.96±0.10bcd

2007

26

369±21.4h

26

211±9.41cd

26

1.72±0.08d

2008

72

376±14.5gh

72

210±4.11cd

72

1.78±0.70cd

2009

16

439±25.0fgh

16

170±3.38e

16

1.85±0.09cd

2010

55

458±24.3def

55

212±6.21cd

55

2.15±0.09ab

2011

33

430±24.5efgh

33

208±8.19cd

33

2.15±0.15ab

2012

39

445±25.8efg

39

224±6.04cd

39

1.99±0.09bc

2013

29

451±22.4def

29

228±5.48cd

29

2.01±0.09abc

2014

60

415±28.0fgh

60

202±8.66de

60

2.07±0.11abc

2015

60

453±17.9def9

60

219±5.68cd

60

2.11±0.06ab

2016

47

528±19.7bcde

47

223±4.78cd

47

2.36±0.07a

2017

61

499±24.5efg

61

223±6.98cd

53

2.20±0.00abc

N = the number of observations, *** P < 0.001, **P < 0.01, means with in the same column with different letters are significantly different, LMY= lactation milk yield, LL= lactation length, DMY= daily milk yield; SE = Standard error; LSM = Least squares means, NS = Non-significant

Variance component and heritability estimates for milk yield traits

The direct heritability estimates for DMY, LL and LMY were 0.33±0.27, 0.20±0.23 and 0.27±0.001, respectively (Table 4). The current estimate for milk yield is lower than the value (0.34±0.13) reported for Ethiopian Boran cattle by Gebreyohannis et al (2014) and that reported by Haile et al (2009). The present direct heritability estimate for lactation length is lower than that (0.26) of Ethiopian Boran cattle reported by Haile et al (2009). However, the heritability of lactation length for Ethiopian Boran and Horo breeds reported by (Gebreyohannes et al 2013) was lower than the present study. The observed moderate level of heritability estimates for milk traits indicates that there is enough variation within the study population and confirms the possibility of genetic improvement through selection.

Table 4. Variance components and heritability estimates for milk traits

Traits

DMY

LL

LMY

σ2a

0.17 ±0.15

1118±1312

0.20±0.0005

σ2pe

0.001±0.31

2692±1110

0.20±0.0003

σ2e

0.34±0.14

1776±1086

0.31±0.0003

σ2p

0.51±0.006

5587±719

0.72±0.0002

h2a

0.33±0.27

0.20±0.23

0.27±0.001

r

0.33

0.48

0.55

Log L

-13838524.8

-678

-36.8

LMY= lactation milk yield, LL= lactation length, DMY= daily milk yield

Repeatability (r) of milk traits

The repeatability estimates for milk yield (Table 4) are lower than reported by Haile et al (2009) for Ethiopian Boran cattle and higher than reported by Gebreyohannes et al (2013) for Horro and their crossbreds. The repeatability estimates of lactation length (Table 4) are lower than the result (0.13 ± 0.17) reported for Mpwapwa cattle (Chawala et al 2017) and that (0.46) noted by Haile et al (2010) for Boran cattle. The estimated repeatability of lactation milk yield and lactation length was highly repeatable while daily milk yield was moderately repeatable. Therefore, selection and culling of animals based on their lactation milk yield performance could be carried out with a reasonable degree of accuracy.

Phenotypic and genotypic correlation among milk traits

Moderate phenotypic and genetic correlations were observed between daily milk yield and lactation milk yield (Table 5). Haile et al (2009) reported higher phenotypic correlation between daily milk yield and lactation length (0.78±0.12); and daily milk yield and lactation milk yield (0.57±0.23) of Ethiopian Boran cattle. The authors also reported a higher genetic correlation (0.57) between daily milk yield and lactation milk yield than the current estimates. The current genetic correlation estimates indicate that selection for one of the traits could have a moderate effect on the other traits through correlated responses.

Table 5. Phenotypic (above diagonal) and genetic correlations (below diagonal) for DMY, LMY, and LL

Parameters

DMY

LMY

LL

DMY

*

0.49± 0.000

0.48 ± 0.003

LMY

0.53 ± 0.001

*

0.5 5± 0.000

LL

0.46 ±0.007

0.53 ± 0.001

*

LMY= lactation milk yield, LL= lactation length, DMY= daily milk yield

Breeding value estimates for milk production trait

The estimated mean breeding values for daily and lactation milk yield and lactation length were 0.02 liters, 10.3 liters and 6.64 days, respectively (Figures 1 and 2). Means of estimates of breeding value by year of birth for the studied traits showed the opposite trend; the highest estimated mean breeding value (0.068 liters) for daily milk yield was observed in 2013 whereas the lowest estimated mean breeding value (-0.02) was observed in 2005. The breeding value for lactation milk yield ranged from -25.2 to 64.1 liters while the breeding value for lactation length was between 2.0 and 11.7 days. The high ranges of breeding values for milk production traits maybe because of the existence of high and low producing cows with the absence of culling and due to non-efficient breeding programs.

Figure 1. Estimated genetic trend for daily and lactation milk yield of Fogera cows at Andassa Livestock Research Center


Figure 2. The estimated genetic trend for lactation length of Fogera cows at Andassa Livestock Research Center

Lactation milk yield increased by 0.86 liters/year while the DMY and LL genetically decreased by 0.0008 liters/year and 0.17 days/year, respectively. The decreasing trend for DMY and LL implies inefficiency in selection based on phenotypic performance and absence of proper culling of unproductive cows. In addition, environmental stresses (heat stress, low quality and quantity of nutrition) and lack of accurate pedigree structure of the herds could be the possible reasons. Researchers and experts should think how to reverse the trend through an appropriate and accurate method of selection.


Conclusions


Acknowledgment

The authors would like to thank Andassa Livestock Research Center and Amhara Region Agricultural Research Institute (ARARI) for financing this research work.


References

Bayou E, Haile A, Gizaw S and Mekasha Y 2015 Evaluation of non-genetic factors affecting calf growth, reproductive performance and milk yield of traditionally managed Sheko cattle in southwest Ethiopia. Springer Plus, 4(1), 568.

Bitew A, Taye M, Kebede A, Mekuriaw G, Tassew A, Mulugeta T and Goshu G 2010 Milk yield and calf growth performance of cattle under partial suckling system at Andassa Livestock Research Centre, North West Ethiopia. Livestock Research for Rural Development. Volume 22, Article #136. http://www.lrrd.org/lrrd22/8/bite22136.htm

Chawala A R, Banos G, Komwihangilo D M, Peters A and Chagunda M G G 2017 Phenotypic and genetic parameters for selected production and reproduction traits of Mpwapwa cattle in low-input production systems. South African Journal of Animal Science, Volume 47, Number 3: 307-319.

Demeke S, Neser F W C and Schoeman S J 2004 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. Journal of Animal Breeding and Genetics, 121(3), 163-175.

Falconer D S and Mackay T F C 1996 Introduction to quantitative genetics. 1996. Harlow, Essex, UK: Longmans Green, 3.

Galmessa U, Dessalegn J, Tola A, Prasad S and Kebede L M 2013 Dairy production potential and challenges in western Oromia milk value chain, Oromia, Ethiopia. Journal of Agriculture and Sustainability, 2(1): 1-21.

Garoma S 2014 Reproductive and productive performance of Kereyu Sanga cattle in Fentalle District of Oromia Region, Ethiopia. Journal of Cell and Animal Biology, 8(2), 28-33.

Gebreyohannes G, Koonawootrittriron S, Elso M and Suwanasopee T 2014 Genotype by environment interaction effect on lactation pattern and milk production traits in an Ethiopian dairy cattle population. Nature and Science, 48, 38-51.

Gebreyohannes G, Koonawootrittriron S, Elzo M A and Suwanasopee T 2013 Variance components and genetic parameters for milk production and lactation pattern in an Ethiopian multibreed dairy cattle population. Asian-Australasian Journal of Animal Sciences, 26(9), 1237-1246

Haile A, Joshi B K, Ayalew W, Tegegne A and Singh A 2009 Genetic evaluation of Ethiopian Boran cattle and their crosses with Holstein Friesian in central Ethiopia: milk production traits. Animal, 3(4), 486-493.

Kiwuwa G H, Trail J C, Kurtu M Y, Worku G, Anderson F M and Darkins J 1983 Crossbred dairy cattle productivity in Arsi region, Ethiopia. ILCA Research Report No. 11, Addis Ababa, Ethiopia. 14-18.

Meyer K 2012 WOMBAT, A program for mixed model analyses by restricted maximum likelihood. User notes. Animal Genetics and Breeding Unit, University of New England Armidale, Australia.

Minuye M 2016 Milk production and reproductive performance of local and crossbreed dairy cows in selected districts of west Gojjam zone, Amhara region, Ethiopia. MSc. Thesis. Bahir Dar University, Ethiopia. p. 42-61.

Rege J E O and Tawah C L 1999 The state of African cattle genetic resources II. Geographical distribution, characteristics and uses of present-day breeds and strains. Journal of Animal Genetic Resources, 26: 1-25.

SAS 2004 SAS user's guide version 9.1: Statistics. Cary, NC: SAS Institute Inc.

Tadesse M and Dessie T 2003 Milk production performance of Zebu, Holstein Friesian and their crosses in Ethiopia. Livestock Research for Rural Development, Volume 15, Article #26. http://www.lrrd.org/lrrd15/3/Tade153.htm

Tesfa A 2015 Genetic and nongenetic parameter estimates for performance traits of Fogera cattle. M.Sc. Thesis, Bahir Dar University, Bahir Dar.


Received 9 December 2019; Accepted 27 January 2020; Published 2 March 2020

Go to top