Genetic parameters and principal components analysis of breeding value for birth and weaning weight in Egyptian buffalo

Objective The objectives of the current study were to study the main environmental factors affecting birth weight (BW) and weaning weight (WW), estimate variance components, genetic parameters and genetic trend and to evaluate the variability and relationships among breeding value of BW and WW using principal components analysis (PCA). Methods A total of 16,370 records were collected from 8,271 buffalo calves. Genetic parameters and breeding values were estimated using a bivariate animal model which includes direct, maternal and permanent maternal effects. These estimates were standardized and used in PCA. Results The direct heritability estimates were 0.06 and 0.41 for BW and WW, respectively whereas direct maternal heritability values were 0.03 and 0.14, respectively. Proportions of variance due to permanent environmental effects of dam were 0.455 and 0.280 for BW and WW respectively. The genetic correlation between BW and WWs was weak approaching zero, but the maternal correlation was 0.26. The first two principal components (PC1 and PC2) were estimated utilizing the standardized breeding values according to Kaiser method. The total variance explained by the first two PCs was 71.17% in which 45.91% and 25.25% were explained by PC1 and PC2, respectively. The direct breeding values of BW were related to PC2 but those of WW and maternal breeding values of BW and WWs were associated with PC1. Conclusion The results of genetic parameters and PCA indicate that BW and WWs were not genetically correlated and improving growth traits of Egyptian buffaloes could be achieved using WW without any adverse effect by BW.


INTRODUCTION
Water buffalo are important source of meat, high quality milk and skin as well as other by products, not only in Egypt but also in many countries worldwide. Thus, buffalos contribute significantly in agriculture economy and human wellbeing in these regions. The total num ber of buffaloes worldwide is about 195.1 million heads, in which 97% are raised in Asia [1] and 2% (~4 million head) are reared in Egypt [1] mainly for milk and meat production. Egyptian buffaloes, known as river buffalo (Bubalus bubalis) contributed 39.2% and 45.2% of national total meat and milk production, respectively [1]. As 96% of buffaloes popula tion in Egypt are owned by small holders, buffaloes play an essential role in rural economy and have challenges in establishing genetic improvement program.
Calf birth weight (BW) is one of the main breeding goals of buffalo breeders, because of its positive relationship with production and reproduction performance. Although, MMI Salem and AMS Amin [2] showed that the stillbirth incidence was similar in either light or heavy weights at birth in the Egyptian buffalo, strong genetic correlations between BW and other growth perfor mance traits including weaning and yearling weights were reported in Pakistani NiliRavi buffalo [3] and Surti buffalo [4]. in addition, strong genetic correlations between BW and several production and reproduction traits such as milk yield, peak daily milk yield, lactation length, age at first calving and calving interval were reported [5]. Accurate estimations of variance components and genetic parameters such as (i.e. heritability and genetic correlation) for Egyptian buffalo growth performance traits are important to design suitable genetic improvement programs. Principal component analy sis (PCA) is a multivariate data technique used to reduce the size of a set of variables by removing repeated informa tion while maintaining the maximum variancecovariance structure of these variables [6]. The PCA has been used in the animal breeding and genetics field to reduce the size of the direct additive genetic covariance matrix in multiple trait models [6,7], and to study the relationship among predicted breeding value [8] as well as population and family structure using genomic data. Thus, the objectives of the current study were to: i) detect the main environmental factors affecting BW [7] and weaning weight (WW); ii) estimate variance components, genetic parameters and genetic trend; iii) eval uate the variability and relationships among breeding value of BW and WWs using PCA .

Dataset
Data were collated from five buffalo experimental herds (El Nattafe ElGadid, ElNattafe ElKadim, Nubariya, Sids and Gimeza) belonging to the Animal Production Research In stitute (APRI), Agriculture research center (ARC), Ministry of Agriculture and Land Reclamation, Egypt. A total of 16,370 records were available on 8,327 buffalo calves in which 8,271 and 8,099 records for BW and WW, respectively, were in cluded in the analysis. These calves were born between 1980 and 2018 with pedigree having 188 sires and 2,211 dams (Table 1).

Statistical analyses
The systematic environmental effects on growth performance traits were examined as fixed effects using least squares meth ods implemented in general linear model procedure of SAS [9]. These fixed effects included the effects of season of calving (4 seasons), year of calving (39 years), sex (male and female), herd (5 herds), gestation length (GL; 4 levels), weight of dam (WD; 5 levels) and parity (14 parities). The effect of GL had 4 classes; between 295 and 325 days with increment of 15 days. Whereas the WD had 5 classes between 300 and 500 kg with increment of 50 kg. The linear model was fitted as follows: Where, Y ijklmnop , the phenotypic record of BW or WW; μ, the effect of the intercept; A i , the fixed effect of ith season of calving; B j , the fixed effect of jth year of calving; C k , the fixed effect of kth sex; D l , the fixed effect of lth herd; G m , the fixed effect of mth gestation length; W n , the fixed effect of nth WD; P o , the fixed effect of oth parity and e ijklmnop , random residual assumed to be independent normally distributed with mean zero and variance σ 2 e . The significant fixed effects were used to form contemporary groups, which were included in ge netic parameters analysis.
Variance components, heritability and breeding values were estimated using the bivariate mixed animal model us ing Wombat software [10]. The model can be described in the matrix notation follows: And the assumptions of the variances were as: fixed effect of n th WD; Po, the fixed effect of o th parity and eijklmnop, random re be independent normally distributed with mean zero and variance σ 2 e. Th effects were used to form contemporary groups, which were included in g analysis.
Variance components, heritability and breeding values were estimated mixed animal model using Wombat software [10]. The model can be descr notation follows: And the assumptions of the variances were as: Where y, a vector of BW or WW observations; b, a vec tor of fixed effects with an incidence matrix X; a, a vector of random animal effects with incidence matrix Z; m, a vector of random dam effects with incidence matrix S; c, a vector of random permanent environmental effects of dam with incidence matrix W; and e, a vector of random residual effects. A is the numerator relationship matrix between animals, I is an identity matrix, σ 2 a is direct additive genetic variance, σ 2 m is maternal additive genetic variance, σ 2 am is the direct × maternal genetic covariance, σ 2 c is maternal permanent environmental variance, and σ 2 e is residual variance. Direct and maternal genetic trends were estimated by regressing yearly mean estimates of breeding value of animal and dam on year of calving.

Principal component analysis
The principal component (PC) could be used as index to eval uate animals for multiple traits [6]. First, in order to avoid the effect of traits' different scaling and magnitudes, the estimated direct and maternal breeding values obtained for BW and WW were standardized for mean zero and unit variance. The standardization was performed according to Boligon et al [6] as follows: z i = (x i -  is the mean of ith trait, and s i is the corre sponding standard deviation. Then, the standardized estimated breeding values (SEBVs) were used in the PCA performed using FactoMineR package [11] in R software. Generally, the PCA reduces the information contained in estimated di rect and maternal breeding values for BWs and WWs in fewer orthogonal latent variables, PCs, with minimal loss of information or variability [6]. In the current study, the PCs explained the highest percentage of variance (i.e with values greater than one eigenvalues) were used. Using the SEBVs in this analysis, each PC can generate a new value, called a principal component score, which is the result of the sum of the SEBVs of each trait weighed by their respec tive standardized score coefficient (SSC) [12]. Thus, the PC could be used to evaluate animals for multiple traits. The SSC of each SEBV in each principal component were ob tained using the following formula:

RESULTS AND DISCUSSION 161
Where PCS jl is the principal component score for the lth animal in the jth principal component, SSC ij is the SSC for EBVs of ith trait in the jth PC and SEBV il is the standardised estimated breeding value of the ith trait for the lth animal.

Environmental factors
The current study assessed the nongenetic factors affected the growth performance (i.e. BW and WW) in Egyptian buffalo ( Table 2). The sex of calf significantly affected (p< 0.0001) BW and WW (p<0.05) where male calves were heavier at birth and weaning than female calves. These dif ferences due to sex effects were in agreement with reported results by Pandya et al [4] in Surti buffalo. These differenc es are probably due to physiological mechanisms associated with sexual endocrine systems which cause the differences between masculinity and feminist characteristics and the natural buildup of male and female bodies [13]. Season of calving also had a significant (p<0.05) effect on both BW and WW where the heavier BW and WW buffaloes calves born in winter compared to those born in other seasons confirming the significant effect of season on buffalo calves' growth performance reported in India [4] and Egypt [14]. Such positive effect of winter calving on BW and WW may be due to the appropriate climatic conditions especially tem perature, and feeding regimes. Furthermore,, the effect of year of calving on BW had significant effect as in the results stated in other studies [4,14]. The effect of year of calving on BW reflects year to year variability with respect to feeding and management practices and climate changes. Moreover, in the current study, the effect of herd was significant (p< 0.0001) on BW but it did show significant effect on WW.
The strong effect of herd may be attributed to differences in management practices among herds during the dam gesta tion period. Also, Gestation length had a significant effect on BW (p<0.05) but it did not have effect on WW. The calves having gestation length between 260 to 310 days had the heaviest BW, followed by those having gestation length longer than 325 days whereas the calves having gestation length below 296 days were the lightest BW. These results were in agreement with what reported in other studies on beef cattle by Jamrozik and Miller [15], where they found a positive genetic correlation between BW and GL. The current study found a significant (p<0.0001) effect of parity on BW whereas there was no effect for parity on WW. Calf BW was lower for primiparous buffaloes than those of multiparous. Fur thermore, BW of calves increased with parity number. This positive relationship has reported in other studies [3,13]. The parity number and age at parturition are reproductive traits which are associated with the physiological status of the fe male buffaloes and positively effect on calves' BWs. The effect of parity on BW may be due to the maturity status of the dams buffaloes as in the late parities, they have high body capacity compatible with better development of fetus [16].
In the current study we also found that, the WD had the same trend for BW, confirming the results reported in pre vious studies [17].

Genetic parameters
The direct heritability estimates were 0.06 and 0.41 for BW and WW, respectively whereas direct maternal heritability values were 0.03 and 0.14, respectively (Table 3). These heri tability values are in the range reported for direct heritability of BW (0.05 to 0.188) and maternal heritability (0.03 to 0.349) on other buffalo populations [4,18,19]. The low direct and maternal heritability values for BW due to maternal environ mental variance being more than 7 folds of direct heritability, can be an indication that genetic selection for BW could be effective, when improving dam feeding and management practices during last phase of pregnancy might play an essen tial role to increase BW of calves [20]. In Brazilian buffaloes calves, Malhado et al [19] reported low heritability estimates 0.09 and 0.03 for direct and maternal heritabilities, respec tively. Also, in Colombian buffalo calves, Bolívar et al [21] found low direct and maternal heritability estimates of 0.05 and 0.01 respectively. On the other hand, Pandya et al [4] and Gupta et al [18] obtained heritability values of 0.188 and 0.349, respectively, for BW of Surti Murrah buffalo. Although, the heritability estimates are population specific, the high esti mates reported in Surti buffalo were probably attributed to the use of animal model without accounting for maternal ef fects which may overestimate the direct heritability estimates. In this study, we found a negative correlation between direct and maternal (r am ) effects for BW (-0.16) which was reported by other studies [22]. This negative genetic correlation could be due to another negative correlation between dam and calf resulted from an adverse effect of high nutrition during early growth of calf [22]. In the current study, high (0.41) and moderate (0.14) di rect and maternal heritability estimates were obtained for WW, respectively indicating that the genetic selection to im prove WW is possible with good response of maternal ability [23]. Such lower estimate of maternal heritability than direct , maternal heritability; C 2 , proportion of phenotypic variance due maternal permanent environmental variance; h 2 t , total heritability; r am , genetic correlation between direct × maternal effect; r g , direct genetic correlation between birth and weaning weights; r m , maternal genetic correlation between birth and weaning weights; r c , maternal permanent environmental correlation between birth and weaning weights.
heritability indicated that the direct genetic effect has large influence on WWs in buffaloes [21]. The current heritability estimates for WW are in agreement with estimated reported in Brazilian buffalo. Nonetheless, the current estimates were higher than direct heritability estimate reported by Ashmawy and ElBramony [24] (0.19) in Egyptian buffalo, and Pandya et al [4] (0.17) in Surti buffalo. Also, in Colombian buffalo, found direct and maternal heritability values of 0.16 and 0.04, respectively [8]. The variation in WW due to maternal envi ronment was moderate but smaller than its value for BW, indicating that the permanent environmental influence of dam on WW has a carryover effect from birth to weaning [25]. The current results showed that the proportion of phe notypic variance due maternal permanent environmental variance (C 2 ) was higher than h 2 m for BW and WW, indicat ing that maternal environmental effect was dominant until weaning [23]. Estimates of C 2 and h 2 m for different beef breeds were close to those obtained in the current study [26]. The directmaternal genetic correlation (r am ) was negative for WW and higher than its value for BW. This negative r am was also found in Australian beef cattle [22] and Spanish beef cattle [27]. The negative genetic correlation may be due to a nega tive correlation between dam and calf resulted from an adverse effect of the high nutrition during early growth of calf [22].
A low genetic correlation between BW and WW (0.05) found in the current study indicated that selection to improve WW in buffalo may not affect BW in a breeding objectives that include both traits. Such relationship is very important to avoid dystocia which may occur for calves with high BW. High estimates of genetic correlation between BW and WW of Egyptian buffalo calves were reported in other studies [18] using animal model which did not account for maternal ef fects. The estimates genetic correlation between growth traits in Beef cattle are usually moderate to high [3,25,27], but others were low and negative such as those by ElSaied et al [28]. The weak correlation between BW and WW may be a de crease in correlation response because of a decrease in selection intensity and population size [29]. Thus, the low genetic cor relation in current study may be attributed to low selection intensity in Egyptian buffalo populations compared to beef cattle. In this study, the maternal genetic r m and maternal permanent environmental r c correlations were 0.262 and 0.003 respectively which were lower than estimated values 0.33 and 0.64 reported in Japanese beef cattle [25].
Remarkable fluctuations in breeding values of calves and dams, and negative direct and maternal genetic trend in some years were found (Figure 1). Averages of direct and maternal genetic trend were -0.003, 0.002, -0.020, and 0.001 kg/yr for BW and WW, respectively indicating that there was no spe cific breeding programs have been applied to improve growth traits in the studied buffalo herds. Therefore, there was no genetic improvement in either trait. Other studies showed improvement in BW and WW where the direct and mater nal genetic trends for BW and WW of Brazilian buffalo were 0.006, -0.01, 0.23, and 0.05 kg/yr, respectively [19]. None theless, a negative direct genetic trend was -0.03 kg/yr was reported for water buffalo in Brazil [30].

Principal components
The total variance explained by PCs was 71.17% in which 45.91% and 25.25% were explained by PC1 and PC2, respec tively (Table 4). These results indicated that the first two PCs out of 4 showed eigenvalues greater than 1. Therefore, the first two components were sufficient to explain the most of variation in breeding value estimates for early growth traits    in the current buffalo calves herds. The variation in breeding values for growth traits have been studied by Boligon et al [31] where they found the first three PCs explained 93.62% of total additive genetic variation of growth traits in Nellore cattle. In Colombian buffalo, three PCs greater explained 65.87% of the variance of breeding value for growth traits [8]. The twodimensional graph of PC1 vs PC2 ( Figure 2) and the correlation between BV, MV of both traits and PC (Table 5) showed that the direct breeding values of BW were correlated with PC2 and those for direct and maternal breed ing values of BW and WW were correlated with PC1. PC1 explained most of direct and maternal breeding value variation and was correlated with the direct and maternal breeding values of WW and maternal breeding values of BW. Thus, it could be considered as genetic index to WW and also as maternal index for BW and WW. Only 25.25% of di rect and maternal breeding value variation was explained by PC2 which was correlated with direct breeding values of BW and considered as genetic index related BW. Therefore, se lection to improve WW could be applied using PC1 separate from selection to improve BW in Egyptian buffalo. These re sults supported the reported genetic correlation between BW and WW in Table 3. The twodimensional graph of PC1 vs PC2 illustrated in Figure 2, showed a negative correlation between maternal and genetic breeding value of BW and WW supported the directmaternal genetic correlations reported in Table 3. A study, on Canchim cattle, reported similar trends where the first two principal components explained 73.37% of total additive genetic variance, and the PC1 was consid ered as genetic index for reproduction traits as well as PC2 was a genetic index correlated with body weight [12].
The PC scores for each animal in the first two principal components were calculated as: The weights of these indices were SSC for each SEBV. The larger the absolute value of the SSC, the greater the relative importance. The 2 PC scores allowing capture the main in formation in predicted breeding value for BW and WW. Animals with higher PCS1 could be used to improve WW, and animals with higher PCS2 could be used to improve BWs. Thus, selection with previous index could be used in genetic selection programs of Egyptian buffalo. As previ ously mentioned regarding the use of PCA is very practical  technique to reduce the data dimensionality especially with limited number of observations in the dependent variables under the lacking of national recording system in Egypt. Also, under Egyptian breeding programs, it has been diffi cult to develop the economic indices and if developed they were rarely used because of the difficulty in economic weights for each trait as the Egyptian market is not stable. Thus, we can consider the PCSs as new traits which can be included in new breeding objectives and used as selection criteria un der Egypt condition. Meanwhile, in other studies the PCSs can be used to reduces the genetic variance components accounting for selection bias and correlated response in cor related traits [32].

CONCLUSION
All environmental factors affect BW, and sex, season of calv ing and herd affected on WW. The results indicated low heritability value for BW, high heritability value for WW and low genetic correlation between BW and WW. More over, 2 PCs components were sufficient to capture the breeding value variation of the studied traits. Finally PC1 and PC2 could be considered as new traits representing WW and BW, respectively. These results suggested improving growth traits of Egyptian buffaloes could be achieved using WW without any adverse effect by BW.