SJAR Spanish Journal of Agricultural Research SJAR 2171-9292 Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA) 13608 10.5424/sjar/2019171-13608 SHORT COMMUNICATION Estimation of genetic trends for body weight traits in Markhoz goat at different ages Latifi Meysam University of Kurdistan, Faculty of Agriculture, Dept. of Animal Science, Sanandaj, Iran. Razmkabir Mohammad University of Kurdistan, Faculty of Agriculture, Dept. of Animal Science, Sanandaj, Iran. should be addressed to Meysam Latifi: m.latifi@agri.uok.ac.ir 01 03 2019 2019 17 1 10.5424/sjar/2019171-13608 19 06 2018 11 02 2019 © 2019 INIA 2018 This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC-by 4.0) License. Abstract

The objective of the present study was to estimate genetic trends for body weight at different ages in Markhoz goat, including birth weight (BW, n = 4758), weaning weight (WW, n= 3685), 6-month weight (6MW, n = 3420), 9-month weight (9MW, n = 3032) and 12-month weight (12MW, n = 2697). Data and pedigree information were collected from 1992 until 2014 at the Breeding Center of Markhoz goat, Sanandaj, Iran. The GLM procedure of SAS was used for selecting the variables and identifying significant fixed effects in the equation of model. Various animal models were applied for genetic analysis and the best model was determined based on Akaike information criteria (AIC). Breeding values of animals were predicted using Wombat program. Genetic trends were obtained by regressing the average predicted breeding values on birth year for each trait. Based on the best model, direct estimated genetic trends were positive and significance for WW, 6MW, 9MW and 12 MW were 15.51, 26.28, 58.36 and 76.70 g/year, respectively ( p <0.001). Maternal genetic trend for BW and WW were 0.61 and 5.47 g/year, respectively ( p <0.01). The low and moderate generic trends obtained in the present study, indicated the possibility of growth traits improvements through genetic selection at all ages in Markhoz goat.

Additional key words: breeding value; growth traits; animal model. Additional key words: AIC (Akaike information criteria); BW (birth weight); WW (weaning weight); 6MW (6-month weight); 9MW (9-month weight); 12MW (12-month weight). The authors received no specific funding for this work.

Author's contributions: Analysis of data: ML. Both authors wrote, read and approved the final manuscript.

Citation Latifi, M.; Razmkabir, M. (2019). Short communication: Estimation of genetic trends for body weight traits in Markhoz goat at different ages. Spanish Journal of Agricultural Research, Volume 17, Issue 1, e04SC01. https://doi.org/10.5424/sjar/2019171-13608

Competing interests: The authors have declared that no competing interests exist.

Introduction

The Markhoz goat is an indigenous breed originating from Kurdistan Province and west Azerbaijan province in Iran. These goats are small-sized with a wide range in brown, white and black coat color. The range of coat color characters is a unique between goats in the world. Most income for farmers is obtained from sale of kids and meat production ( Rashidi et al ., 2011 , 2015) .

The profitability of goat production for meat and increasing of meat production in goats depend on kid’s weight as the growth performance and production of more kids per doe. Accurate prediction of breeding value of animals is one of the best tools available to maxi­mize response to selection programs and the success of a breeding program can be assessed by examining the actual change in breeding value expressed as a proportion of expected theoretical change of the breeding value mean for the trait under selection ( Jurado et al ., 1994 ). The accuracy of genetic evaluations depends on how well the assumptions of the model match the data ( Wiggans & VanRaden, 1991 ). Selection methods and management can be evaluated by estimating genetic and environmental trend. However, few researches have reported evaluation of genetic trends for body weight traits in goat ( e.g ., Snyman, 2012 ; Hassan et al ., 2013 ; Hasan et al ., 2014 ).

Several studies have been conducted on Markhoz goat including the estimation of genetic parameters or inbreeding effects ( Rashidi et al ., 2008 ; Kheirabadi & Rashidi, 2016 ; Mahmoudi et al ., 2017 ), but there is not any information about genetic trends for body weight traits at different ages in this breed. Therefore, the objective of this study was to determine the best models of genetic analysis for body weight of Markhoz goat kids at different ages and estimation of genetic trends in body weight traits over a 23-year period.

Material and methods Data collection and management

The data used in this study were collected between 1992 and 2014 from Breeding Center of Markhoz goat in Sanandaj, Kurdistan, Iran (1373 m asl and 35 ° 20' N latitude and 47 ° E longitude). The studied traits were birth weight (BW), weaning weight (WW), 6-month weight (6MW), 9-month weight (9MW) and 12-month weight (12MW). Details on the climatic conditions and herd management are presented by Rashidi et al . (2015) . Briefly, in the station, does were mated with selected bucks at about 18 months of age in breeding pens. Each buck was randomly mated to ∼20 does from early October to late November. Kidding was in February and March and kids were weaned until ∼4 months of age. The kids were weighed and ear tagged after birth. The number of records, mean, standard deviation, coefficient of variation, number of sires and dams for the body weight traits are presented in Table 1.

Descriptive statistics for body weight traits in Markhoz goat. Statistical and genetics analysis

The GLM procedure ( SAS, 2001 ) was used for identifying fixed effects that had significant influence on the investigated traits. This was performed on a model including fixed effects of year of birth (23 levels from 1992 to 2014), age of doe (6 classes from 2 to 7 years old), kid’s gender (2 classes of male and female) and birth type (singles, twins and triplets). All of these fixed effects were significant ( p ≤ 0.001) for all body weight traits (BW, WW, 6MW, 9MW and 12MW) and therefore these effects were considered in the model. Age of kid at weighing (in days) was considered as covariate for all analyzed traits except BW.

(Co)variance components and corresponding ge­­ne­tic parameters were estimated using Wombat program ( Meyer, 2013 ). Six different animal models were used for genetics analysis for each trait, by ignoring or including maternal additive genetic effect, permanent environmental effect and covariance between direct-maternal additive genetic effects.

y = X b + Z 1 a + e (1)

y = X b + Z 1 a + Z 3 pe + e (2)

y = X b + Z 1 a + Z 2 m + e cov (a, m) = 0 (3)

y = X b + Z 1 a + Z 2 m + e cov (a, m) = A σ am (4)

y = X b + Z 1 a + Z 2 m + Z 3 pe + e cov (a, m) = 0 (5)

y = X b + Z 1 a + Z 2 m + Z 3 pe + e cov (a, m) = A σ am (6)

where y is a vector of observations of the studied traits; b , a , m , pe and e are vectors of fixed effects, direct additive genetic effects, maternal, permanent environmental effects and the residual effects, respectively. X , Z 1 , Z 2 and Z 3 are incidence matrices relating these effects to the records. A is the additive numerator relationship matrix, and σ am is the covariance between direct additive genetic and maternal additive genetic effects. Assumptions for variance (V) and covariance (Cov) matrices involving random effects were:

V (a) = A σ a , V (m) = A σ m , V (pe) = I d σ pe , V (e) = I n σ e and Cov (a,m) = A σ am

where, I d and I n are identity matrices that have order equal to the number of dams and records, respectively, and σ a , σ m , σ pe and σ e are variance of additive genetic values, maternal additive genetic varian­ce, maternal permanent environmental variance and residual va­riance, respectively.

Model comparison criteria

Akaike information criteria (AIC) were applied to choose the best-fit model amongst all six models ( Akaike, 1974 ). It is defined as AIC = -2 log L + 2 p, where log L is the maximized likelihood and p is the number of parameters in the model. The model with the smallest AIC is considered as the best model for each trait.

Predicted breeding values of individual animals were obtained with Wombat ( Meyer, 2013 ). Genetic trends of the studied traits were obtained by regression of average predicted breeding values on birth year of animals for each trait. Genetic trend analyses were carried out with the regression procedure (Proc Reg) of SAS program ( SAS, 2001 ).

Results and discussion

The most appropriate model for each trait and parameter estimates are given in Table 2. The results obtained in this study indicate that the most appropriate models for BW, WW, 6MW, 9MW and 12 MW were (6), (5), (1), (2) and (1), respectively. Therefore, these models were selected for parameter estimates and prediction of breeding value for each trait.

The most appropriate models, parameter estimates and estimates of genetic trends (g/year) for body weight traits in Markhoz goat.

The average of direct and maternal predicted breeding values by year of birth during the study period (1992-2014) are illustrated in Fig. 1. The estimates of genetic and maternal trends (g/year) for investigated traits are reported in Table 2. As illustrated in Fig. 1, the direct genetic trend of BW between 1992 and 2014 appeared to be constant, but generally showed an increase over time for the other studied traits (WW, 6MW, 9MW and 12MW). The average of annual genetic trends were positive from 1998 onward, and had irregular variation over the years. The annual direct genetic trends were positive and significant ( p <0.001) for all growth traits except BW over the years (Table 2). Fig. 1 shows that the maternal genetic trend appeared to be constant and fluctuated for BW and WW traits, respectively. The annual maternal genetic trends were positive and significant for BW ( p <0.01) and WW ( p <0.001) over the years (Table 2).

Direct and maternal genetic trends of mean breeding values by year of birth observed for body weight traits of Markhoz goat.

Estimates of direct genetic trends for BW, WW, 6MW, 9MW and 12 MW were 1.08, 15.51, 26.28, 58.36 and 76.70 (g/year), respectively. The direct genetic trend for birth weight in the current study (1.08 g/year) was low ( p >0.05) and generally agrees with that reported by Ghavi Hossein-Zadeh (2012) in Moghani sheep (1.63 g/year), Hassani et al . (2009) in Baluchi sheep (0.7 g/year) and Zishiri et al . (2010) in Ile de France sheep (1 g/year). However, higher estimates of genetic trends for BW in various other goat and sheep breeds were reported by Snyman (2012) in Angora goat (4 g/year), Hasan et al . (2014) in Ettawa Grade goat (20 g/year), Kariuki et al . (2010) in Dorper sheep (6 g/year) and Supakorn et al . (2013) in sheep population in Thailand (20 g/year). This difference in genetic trends can be explained by dissimilarity in breeding strategies for different breeds and models used for analysis.

Estimate of direct genetic trend for WW in the pre­sent study was 15.51 g/year ( p <0.001). Also, estimate of direct genetic trend for WW was lower than those reported by Snyman (2012) in Angora goat (56.7 g/year), Hasan et al . (2014) in Ettawa Grade goat (60 g/year) and Shrestha et al . (1996) in Suffolk (23 g/year) and Finnsheep (25 g/year) sheep. Lower estimate (7 g/year) was reported by Gholizadeh & Ghafouri-Kesbi (2015) in Baluchi sheep.

The direct genetic trend estimate for 6MW in the current study (26.28 g/year) was greater than those reported by Hasan et al . (2014) in Ettawa Grade goat (-0.10 g/year), and Gholizadeh & Ghafouri-Kesbi (2015) in Baluchi sheep (4 g/year). Higher estimates of 91 g/year were reported by Mokhtari & Rashidi (2010) in Kermani sheep, 79.4 g/year in Moghani sheep ( Ghavi Hossein-Zadeh, 2012 ) and 72 g/year was reported by Hassani et al . (2009) in Baluchi sheep.

The direct genetic trend value estimated for 9MW (58.36 g/year) in the present study was higher than those of 35 g/year reported by Snyman (2012) in Angora goat at eight month body weight, 7 g/year was reported by Gholizadeh & Ghafouri-Kesbi (2015) in Baluchi sheep and it was lower than those of 81 g/year in Kermani sheep ( Mokhtari & Rashidi, 2010 ) and 66.83 g/year was reported by Ghavi Hossein-Zadeh (2012) in Moghani sheep.

The genetic trend value estimated for 12MW ob­tai­­­ned in the present research (70.76 g/year) was higher than those of -0.30 g/year and 14 g/year reported by Hasan et al . (2014) in Ettawa Grade goat and Gholizadeh & Ghafouri-Kesbi (2015) in Baluchi sheep, respectively.

Also, estimate of genetic trend for 12MW in the present study was lower than the estimates of 156 g/year and 110.2 g/year reported by Mokhtari & Rashidi (2010) and Ghavi Hossein-Zadeh (2012) , respectively. In the present study an increasing trend was observed for direct genetic trend value from birth weight to 12MW. The increasing direct genetic trend value of body weight at later ages might be explained by attributed to the relatively higher additive genetic variation. Also, Shaat et al . (2004) reported that higher genetic trend may be explained by the higher additive genetic variation.

Also, maternal genetic trends for BW (0.61 g/year)

and WW (5.47 g/year) were positive ( p <0.01 and p <0.001, respectively). Results of the current study showed that maternal effects influences on pre-weaning traits in Markhoz goat and this effect need to be considered in selection programs. Maternal genetic trend estimates for BW and WW are rarely reported in literature. The maternal genetic trends for BW in current study were low and generally in agreement with estimates of 1 and 2.36 g/year reported by Snyman (2012) and Ghavi Hossein-Zadeh (2012) , respectively. Estimate of maternal genetic trend for WW was lower than those observed by Snyman (2012) in Angora goat (9 g/year) and Ghavi Hossein-Zadeh (2012) in Moghani sheep (49.2 g/year). The estimate of maternal trends for BW and WW were low and these values were lower than direct trends. These may be due to the larger direct effects on BW and WW than maternal genetic effects in Markhoz goat. Rashidi et al . (2008) reported lower maternal heritabilities estimates than direct heritabilities for BW and WW traits in Markhoz goat.

Results indicate the positive and significant genetic trends for weight traits of Markhoz goat. The low and moderate genetic trends obtained in Markhoz goat indicated that it can be possible the improvement of growth traits through genetic selection at all ages considered.

 Acknowledgement

The authors are grateful to the Breeding Center of Markhoz goat in Sanandaj, Kurdistan, Iran for providing the data and information.

 References Akaike H 1974 A new look at the statistical model identification. IEEE Trans Autom Control 19 716 723 https://doi.org/10.1109/TAC.1974.1100705 Ghavi Hossein-Zadeh N 2012 Bayesian estimates of genetic changes for body weight traits of Moghani sheep using Gibbs sampling. Trop Anim Health Prod 44 531 536 https://doi.org/10.1007/s11250-011-9930-1 Gholizadeh M Ghafouri-Kesbi F 2015 Estimation of genetic parameters for growth-related traits and evaluating the results of a 27-year selection program in Baluchi sheep. Small Rumin Res 130 8 14 https://doi.org/10.1016/j.smallrumres.2015.07.032 Hasan F Jakaria Gunawan A 2014 Genetic and phenotypic parameters of body weight in Ettawa Grade goats. Media Peternakan 37 8 16 https://doi.org/10.5398/medpet.2014.37.1.8 Hassan MR Sultana S Iqbal A Talukder MA 2013 Estimation of heritability, breeding values and genetic trends for growth traits of exotic goat. Int J Nat Sci 3 7 11 Hassani S Deltang Sefidsanghi H Rashidi A Ahani Azari M 2009 Estimation of genetic, phenotypic and environment trends for some growth traits in Baluchi sheep. J Agr Sci Nat Resour 16 126 132 [In Persian]. Jurado JJ Alonso A Alenda R 1994 Selection response for growth in Spanish Merino flock. J Anim Sci 72 1433 1440 https://doi.org/10.2527/1994.7261433x Kariuki CM Ilatsia ED Kosgey IS Kahi AK 2010 Direct and maternal (co)variance components, genetic parameters and annual trends for growth traits of Dorper sheep in semi-arid Kenya. Trop Anim Health Prod 42 473 481 https://doi.org/10.1007/s11250-009-9446-0 Kheirabadi K Rashidi A 2016 Genetic description of growth traits in Markhoz goat using random regression models. Small Rumin Res 144 305 312 https://doi.org/10.1016/j.smallrumres.2016.10.003 Mahmoudi P Rashidi A Razmkabir M 2017 Inbreeding effects on some reproductive traits in Markhoz goats. Anim Prod Sci 58 2178 2183 https://doi.org/10.1071/AN17043 Meyer K 2013 WOMBAT-A program for mixed model analyses by restricted maximum likelihood. Genetics and Breeding Unit, Armidale, 105 Mokhtari MS Rashidi A 2010 Genetic trends estimation for body weights of Kermani sheep at different ages using multivariate animal models. Small Rumin Res 88 23 26 https://doi.org/10.1016/j.smallrumres.2009.11.003 Rashidi A Sheikhahmadi M Rostamzadeh J Shrestha J 2008 Genetic and phenotypic parameter estimates of body weight at different ages and yearling fleece weight in Markhoz goats. As-Aust J Anim Sci 21 1395 1403 Rashidi A Bishop SC Matika O 2011 Genetic parameter estimates for pre-weaning performance and reproduction traits in Markhoz goats. Small Rumin Res 100 100 106 https://doi.org/10.1016/j.smallrumres.2011.05.013 Rashidi A Mokhtari MS Gutiérrez JP 2015 Pedigree analysis and inbreeding effects on early growth traits and greasy fleece weight in Markhoz goat. Small Rumin Res 124 1 8 https://doi.org/10.1016/j.smallrumres.2014.12.011 SAS 2001 SAS /STAT user's Guide: statistics. Release 8.2. SAS Institute Inc., Cary, NC, USA. Shaat I Galal S Mansour H 2004 Genetic trends for lamb weights in flocks of Egyptian Rahmani and Ossimi sheep. Small Rumin Res 51 23 28 https://doi.org/10.1016/S0921-4488(03)00185-8 Shrestha JNB Peters HF Heaney DP Van Vleck LD 1996 Genetic trends over 20 years of selection in the three synthetic Arcotts, Suffolk and Finnish Landrace sheep breeds. 1. Early growth traits. Can J Anim Sci 76 23 34 https://doi.org/10.4141/cjas96-004 Snyman MA 2012 Genetic analysis of body weight in South African Angora kids and young goats. S Afr J Anim Sci 42 146 155 https://doi.org/10.4314/sajas.v42i2.7 Supakorn C Pralomkarn W Anothaisinthawee S 2013 Estimation of genetic parameters and genetic trends for weight and body measurements at birth in sheep populations in Thailand. Songklanakarin J Sci Technol 35 1 10 Wiggans G VanRaden P 1991 Method and effect of adjustment for heterogeneous variance. J Dairy Sci 74 4350 4357 https://doi.org/10.3168/jds.S0022-0302(91)78631-1 Zishiri OT Cloete SWP Olivier JJ Dzama K 2010 Genetic trends in South African terminal sire sheep breeds. S Afr J Anim Sci 40 455 458