Estimation of Genetic Parameters for Live Weight in South African Holstein Cattle ()
1. Introduction
The Holstein is the most popular dairy cattle breed in South Africa. Presently, about 53,151 Holstein cows are participating in the South African National Milk Recording Scheme [1] . The average milk production of these cows per 305 day lactation is 7441 kg compared to Ayrshire, Guernsey and Jersey which produce 6072 kg, 5570 kg and 5187 kg respectively [1] .
Continued decline in dairy herd profitability poses a serious challenge to the South African dairy industry. Due to this problem, the number of milk producers in South Africa has declined from 3899 in January 2007 to 2083 in September 2013 [2] . There is, therefore, a pressing need to improve the profitability of milk production in South African dairy herds. In order to achieve such improvement, from a genetic improvement perspective, all economically important traits should be included in the selection objective. Live weight is among the traits of economic importance in dairy production, as it affects dairy herd profitability through its effect on cow maintenance costs and income from beef (culled cows). Cows with increased live weight tend to have higher maintenance requirements, which are not compensated for by higher revenue from the sale of beef (culled cows), resulting in negative economic value [3] - [5] . In South African dairy cattle, a kilogram increase in live weight was recently reported to result in a decrease in profit of up to ZAR 7.49 [6] . In the past, little or no emphasis was placed on non-production traits, such as live weight, in most breeding objectives worldwide [7] . In South Africa, genetic trends show that dairy cattle have been selected mainly for increased yield and improved type [8] [9] . Such narrow breeding objectives hamper improvement in total economic merit. In recent years, there has been a shift towards more balanced, profit-focused breeding objectives, globally, by including all economically relevant traits [7] .
The aim of this study was to estimate genetic parameters for live weight in South African Holstein cattle. Such parameters are a prerequisite to incorporating live weight into the selection objective of South African Holstein cattle.
2. Materials and Methods
Two herds that routinely record live weight data were identified by South African National Milk Recording Scheme technicians. The first herd was Groenpunt Correctional Services Dairy Farm, situated in Vereeniging, Gauteng Province, which has an average of 55 lactating cows. The second herd was Elsenburg Experimental Station, Western Cape Department of Agriculture, which is situated in Elsenburg, Western Cape Province with an average of 162 lactating cows. Only live weights of lactating Holstein cows were used. Cows were weighed daily after the morning milking using weighing scales approved by the National Livestock Improvement Scheme. Cows at Groenpunt Correctional Services Dairy Farm were fed 32 kg of a total mixed ration per cow per day.The diet consisted of lucerne hay, yellow maize, silage, brewers grain, weeping love grass and salt. After milking, cows are provided an average of 9 kg of dairy meal per cow per day depending on production. Moreover due to different production stages cows are housed separately in sheds.
The feeding program and management for Elsenburg Experimental Station was as described by [10] . Live weight data were merged with milk production and pedigree data from the Integrated Registration and Genetics Information System of South Africa (Intergis). Live weights collected between 1 and 305 days of lactation were selected. Cows were aged between 24 and 96 months and ranged from lactation 1 to 8. Data editing was carried out using the Statistical Analysis System software [11] . The original data set consisted of 29,079 live weight records of 195 cows. Herd-year-season of weighing (HYS) was defined as the contemporary group. Data were edited to remove records without live weight, birth date, calving date or identification number. After editing, the remaining data set had 9843 live weight records on 94 cows.
2.1. Environmental Factors Affecting Live Weight
An analysis of variance (ANOVA) was conducted to determine environmental factors affecting live weight; such effects need to be fitted in the model for variance component estimation. The General Linear Models (GLM) procedure of the Statistical Analysis System software [11] was used and the following environmental factors were tested for:
・ Herd-year-season of weighing;
・ Parity;
・ Age at weighing;
・ Days in milk at weighing;
・ Quadratic effects for days in milk;
・ Two way interactions between all the above effects.
2.2. Estimation of Variance Components for Live Weight
Variance components for live weight were estimated by the Restricted Maximum Likelihood (REML) procedure, using the ASREML programme [12] . The following single trait repeatability animal model was used:
(1)
where:
y = vector of observations of live weight;
X = incidence matrix relating fixed effects to observations;
b = vector of fixed effects influencing live weight;
Z = incidence matrix relating random animal additive genetic effects to observations;
u = vector of animal additive genetic effects;
W = incidence matrix relating random permanent environmental effects to observations;
pe = vector of permanent environmental effects, to account for effects influencing the repeated live weight records;
e = is the vector of random residual effects.
Random animal additive genetic effects (a) were assumed to have the distribution, where A is the additive genetic relationship matrix and is the animal additive genetic variance. Residual effects were assumed to be normally, independently and identically distributed with mean 0 and variance (i.e. , and cov(a, e) = 0. Permanent environmental effects were assumed to be distributed with, where I is an identity matrix, is the variance due to permanent environmental effects and cov(pe, e) = 0 and cov(a, pe) = 0.
Thus, the (co)variance structure of the random effects was assumed to be as follows:
(2)
2.3. Calculation of Phenotypic Variance, Heritability and Repeatability
The phenotypic variance of live weight was calculated as the sum of additive genetic variance, permanent environmental variance and residual variance (Falconer & Mackay, 1996) :
(3)
The heritability of live weight (h2) was calculated as the ratio of animal additive genetic variance to phenotypic variance (Falconer & Mackay, 1996):
(4)
The repeatability of live weight (r) was calculated as the sum of additive genetic variance and permanent environmental variance to the phenotypic variance (Falconer & Mackay, 1996):
(5)
3. Results and Discussion
An overall mean of 570 ± 0.80 kg (Table 1) was obtained for live weight across parities. The same mean was obtained for Holstein cattle in Ireland [13] . Another study in Ireland [14] reported a higher mean for pre-calving Holstein cows. Higher means of 587 kg at parturition and 641 kg in third parity, respectively, were also found in United States Holstein cattle [15] [16] . Danish Friesian cattle were also heavier at 593 kg [17] . Using data from
Table 1. Summary of the descriptive statistics.
Elsenburg Experimental Farm, one of the herds used in the current study, but from earlier years, [10] obtained a lower mean of 533 kg. The current study analysed data on eight parities, whereas [10] only focused on the first
three parities. The higher mean live weight for the current study may partly be due to a correlated increase in live weight with increase in yield traits in this herd, over time. There has been a positive genetic trend in yield traits in recent years in the herd [18] and a high positive correlation is known to exist between live weight and yield traits in dairy cattle [10] . A high positive correlation is known to exist between live weight and yield traits in dairy cattle [10] .
Studies conducted in New Zealand and the Netherlands reported lower mean live weights of 404 kg and 520 kg, respectively [19] [20] in Friesian cattle. The low live weight for New Zealand Friesians is probably due to sustained negative selection on live weight. New Zealand is one of the few countries that include live weight, with negative emphasis, in the selection objective [21] .
3.1. Environmental Factors Influencing Live Weight
Environmental effects significantly affecting live weight (P < 0.001) were herd-year-season of weighing, parity, age at weighing, days in milk and quadratic effects of days in milk as illustrated in Table 2. Several other studies [10] [13] [17] [19] [22] [23] have also reported parity to be a significant source of variation in live weight.
Herd-year-season of weighing corrects for the effects associated with yearly and seasonal climatic variation and management differences among herds [24] . Age at weighing also had a significant effect (P < 0.001) on live weight in a study on Holstein cows in Finland [22] . Alawneh [19] also observed a significant effect (P < 0.05) of days in milk on individual live weights record between 0 - 100 days of lactation.
3.2. Variance Component Estimates
Estimating variance components for live weight in the South African Holstein cattle population forms part of an objective for this study. Such estimates are a prerequisite to incorporating live weight into the selection objective of South African Holstein cattle. Estimates of variance components, heritability and repeatability for live weight are presented in Table 3. Live weight had high heritability (0.74) and repeatability (0.86) estimates. The high repeatability indicates a high influence of permanent environmental effects on repeated measures of live weight in South African Holstein cattle. Toshniwal et al. [25] obtained a slightly lower repeatability of 0.76 in United States Holstein cows.
The heritability was similar to estimates reported in Norway and the Netherlands, which ranged between 0.71 and 0.88 [26] [27] . Such a high estimate of heritability implies that there is scope for significant response to selection on live weight in the South African Holstein cattle population. Veerkamp et al. [28] , Berry et al. [29] and Muller et al. [10] found moderate heritabilities (0.60 - 0.61) in South African, Irish and Netherlands Holstein cows, respectively. Fairly lower heritability estimates, between 0.40 and 0.50, were reported by [23] [25] [30] [31] . Distinctly lower heritabilities of 0.24 - 0.36 were found in other studies [15] [32] - [34] . Differences in heritability estimates among studies may, in part, be due to the variation in stages at which live weight was measured. For example, [28] used live weights in parity 1 and [10] used live weights measured in lactations 1 to 3. Veerkamp and Brotherstone [30] measured live weight at 26 weeks of lactation. Toshniwal et al. [25] measured daily live weight at 1 - 26 weeks and [31] predicted live weight from type traits at 29 - 33 weeks of lactation. Heritability is also a parameter that varies from one population to another.
4. Conclusion
Results of this study may assist in including live weight in the selection objective for South African Holstein cattle. Herd-year-season of weighing, parity, age at weighing, days in milk and quadratic effects for days in milk significantly influence live weight in the two South African herds used in the current study. These factors should, therefore, be accounted for when analysing live weight data. The high estimate of heritability of live weight
Table 2. Environmental factors influencing live weight in South African Holstein cattle.
*Significant (p < 0.001); NS, Not significant.
Table 3. Variance components, heritability & repeatability of live weight.
Heritability; Additive genetic variance; permanent environmental variance; phenotypic variance; repeatability; residual (error) variance; SE = standard error.
obtained in the present study indicates that there is potential for significant rate of genetic change if selection is applied on live weight in this population. Variance components, on which this heritability estimate is based, can
be used to calculate estimated breeding values for live weight; thus making it possible to include live weight in the in breeding objective for South African Holstein cattle. It should be noted, however, that these estimates are based on only two herds and therefore are not representative of the whole population. There is therefore a need to build on the results of this study by collecting and analysing live weight from a much larger proportion of South African Holstein herds.
Acknowledgements
The authors would like to acknowledge the Agricultural Research Council and National Research Foundation for providing financial support to this study. We wish to thankDr CJC Muller (Elsenburg Experimental Station) and staff at Groenpunt Correctional Services Dairy Farm for kindly providing us with the live weight records and other relevant information.
NOTES
*Corresponding author.