Analysis of (co) variance components and estimation of breeding value of growth and production traits in Dahlem Red chicken using pedigree relationship in an animal model

Variance and covariance components of growth and production traits were analyzed employing REML animal model to assess the Dahlem Red (PD-3) chicken population for direct additive genetic, maternal effects and to estimate the estimated breeding value (EBV), genetic parameters, genetic trends and rate of inbreeding (ΔF) utilizing seven generation’s data. The generation and hatch had significant (P≤0.01) effect on the body weight at 0 day (BW0), 2 (BW2), 4 (BW4) and 6 weeks (BW6) and shank length at six weeks of age (SL6). The average least squares means (LSM) for BW6 and SL6 were 273.93±0.62 g and 53.97±0.05 mm, respectively. All the production traits were significantly (P≤0.01) influenced by generation and hatch. The average LSM for age at sexual maturity (ASM), egg production up to 40 weeks (EP40) and egg mass up to 40 weeks (EM40) were 168.82±0.25 d, 72.60±0.41 eggs and 4.21±0.07 kg, respectively. Model 5 with additive direct, maternal genetic, maternal permanent environmental and residual variance components was the best for BW0, BW2 and BW4 based on the AIC values obtained in WOMBAT. Model 4 was the best model for BW6, SL6, ASM, EP40 and EM40 with additive direct, maternal permanent environmental and residual variance components. Maternal effects were higher during early age, decreased with age, and remained present until 20 weeks of age. The heritability (h2) estimates were low to moderate in magnitude for all the growth traits and ranged from 0.02±0.03 to 0.19±0.03. The maternal heritability was high at hatch (0.35±0.06), decreased gradually until 4th week (0.02±0.01) and ceased afterwards. The heritabilities of EP40 (0.11±0.03) and EM40 (0.12±0.04) were low. The direct additive genetic correlations (ra) between BW2, BW4, BW6 and SL6 were high and positive (P≤ 0.05). The additive genetic and maternal permanent environmental correlation between EP40 and EM40 were high and positive (P≤ 0.05). The EBV of EM40 was significant (P≤ 0.05) with 0.48 kg/generation in PD-3 chicken at the end of the seventh generation. The EBV of EP40 showed an increasing trend with a genetic gain of 1.87 eggs per generation. The average inbreeding coefficient of the population was 0.019 and average ΔF was 0.007 over the last seven generations of selection. The EBV trends for primary and associated traits showed linear trends in the desired direction and negligible inbreeding.

Introduction Poultry production has grown exponentially in structure and operation from a traditional backyard activity into a major commercial agri based industry over the years. Development of high yielding layer (310-340 eggs) and broiler (2.4-2.6 kg at 6 weeks) varieties together with standardized practices on nutrition, housing, management and disease control [1] have contributed to spectacular growth rates in the egg (8.5% per annum) and broiler production (7.8% per annum) in India [2]. Backyard poultry contributes about 17% (17.09 billion) to the total egg production (103.32 billion) of the country and chicken meat is 50% (4.06 MMT) of the total meat production (8.10 MMT) in India [2]. India ranks 3 rd in egg production and 5 th in meat production in the world, but the availability of chicken products is far below the recommended levels. The per capita availability is about 79 eggs and 3.2 kg chicken meat per annum against the recommended level of 180 eggs and 10.25 kgs meat as per the Indian Council of Medical Research. Backyard poultry farming with improved chicken varieties such as Vanaraja, Gramapriya, Srinidhi, and Rajasree is promoted by the Government as it is one of the proven technologies for poverty alleviation and nutritional security in rural and tribal areas of the country. A wide gap exists between the availability and requirement of chicken products, thus providing an opportunity to expand the poultry farming in the commercial as well as the backyard sectors. Since commercial poultry has inherent difficulties for expansion in rural and tribal areas, the backyard poultry has bright scope in these areas. Promoting the backyard poultry farming among the rural masses is going to increase its share in total poultry production apart from providing healthy animal protein in terms of eggs and meat. Dahlem Red (PD-3) is the female parent of Gramapriya, an improved backyard chicken layer variety propagated across the country. The improvement in egg production, egg mass, growth and other economic traits in the terminal cross is being achieved through the application of selection in the parent lines.
The estimated breeding value (EBV) determines the rate of improvement in the primary trait of selection. Variability in the population is the basis for any genetic improvement program in livestock and poultry. Genetic improvement is determined by the response to selection in the primary as well as other correlated traits of economic importance. Thorough knowledge on the genetic basis, gene action and the environmental effect is essential for planning breeding programs for genetic improvement of the desired traits in poultry. Selection response of trait depends on the population size, gene frequency, mutation, allelic fixation, random drift and physiological limits [3]. The magnitude and direction of correlated responses play a major role in developing effective breeding strategies for improved productivity, more so among positively correlated traits due to linkage and pleiotropy [4]. Partition of variance into every possible source and their interaction effects reduces the error variance and minimizes over-weightage to some components, leading to precision in the estimates. Animal model accommodates additional components of variance, which makes the estimates less biased by confounding environmental variation and precise modeling [5]. The traditional models ignore the maternal and permanent environmental effects in chicken leading to overestimation of additive genetic variance resulting in high h 2 estimates. The maternal effect is defined as the situation where the phenotype of the offspring is determined not only by the environment and its genotype, but also by the genotype and environment of their mother. A maternal effect is a situation wherein traits are influenced by the genotype and environment of the mother. Maternal effects are important in the development and expression of the economic traits due to genetic or environmental differences between dams or by the combination of genetic or environmental differences [6,7]. The egg weight is a major factor contributing to maternal effects for body weight, however, others like hatch weight, incubation conditions, nutrition, etc. may also play a role. Correlations of egg weight with hatching weight and with the subsequent body weight of chicks may also reflect maternal effects. The inclusion of maternal effects in the model reduces the bias and increases the precision of genetic parameters [8]. Many authors estimated the genetic parameters without considering the maternal effects in chicken [9][10][11][12][13][14][15][16]. Some publications with additive, maternal and permanent effects on various economic traits using animal model [4,17] and diallel analysis [10,18] were reported in chicken. However, studies with large data sets over the generations utilizing an animal model in PD-3 line are limited from India. The analysis of data using traditional models overestimated the genetic parameters and lead to inadequate inferences affecting the long term breeding experiment especially in India.
Therefore, the present study was taken up with an aim to estimate the breeding value, genetic trends of economic traits, genetic (additive and maternal) and non-genetic environmental effects and performance of important economic traits in (PD-3 line using Average Information Restricted Maximum Likelihood (AIREML) animal model.

Materials and methods
The study was carried out at the experimental poultry farm of ICAR-Directorate of Poultry Research, Hyderabad, Telangana, India. Hyderabad is in Deccan plateau placed between 172 3' N and 78˚28' E at 500 m from mean sea level. The climate is hot tropical with maximum temperature ranging from 20˚C in winter to 45˚C in summer.

Ethic and biosafety statement
The experiment was approved by the Institutional Animal Ethics Committee (IAEC/DPR/17/ 4). The manuscript does not contain clinical studies and patient data.

Population and management
PD-3 population was selected over the generations [19,20] for higher egg production and egg mass. It is the female parent line of Gramapriya, a popular brown egg layer suitable for rural poultry farming developed by ICAR-Directorate of Poultry Research. PD-3 line was developed from Dahlem Red breed. The work on the development of PD-3 line as a female line for production of egg type rural poultry was started in 1998-99. A pedigreed population of this line was produced for the first time (S-0) in 2000-01. After recording the performance for 6 generations, a random mating program using pooled semen for two generations was followed. The variability in the existing population was increased subsequently with the addition of Dahlem Red population from Anand Agriculture University, Anand, Gujarat, India. The base population was stabilized for two generations with random inter se mating and the third generation was considered as the base population with the initiation of pedigreed mating and selection. The population was selected for higher 40 week egg mass (EM40) during the last seven generations.
About 3500 to 4000 healthy chicks were produced in 7-8 hatches in each generation in a pedigreed mating with 50 sires and 250 dams. These chicks were wing banded and reared on deep litter in an open sided poultry house under standard management. About 600 females and 200 males were housed in individual cages in each generation at 12-14 weeks of age. The chicks were fed ad-libitum layer starter diet (2800 Kcal: ME and 18%: CP) based on maize-soybean meal up to 6 weeks of age in each generation, with layer grower ration (2700 Kcal: ME and 18% CP) until 16 weeks of age and with layer breeder ration (2650 Kcal: ME, 16.50% CP and Calcium: 3.5%) until 40 weeks of age. The chicks were vaccinated against Marek's disease (MD) on the 1 st day; Newcastle disease (ND) on the 5 th and 14 th days; infectious bursal disease (IBD) on the 24 th day and fowl pox (FP) on the 42 nd day, ND with R 2 B at the 9 th week, infectious coryza at the 12 th week and ND, IBD and infectious bronchitis (IB) combined vaccine at the 18 th and the 40 th weeks of age.

Data and traits
Data on growth and production traits of 16,292 PD-3 birds descended from 336 sires and 1076 dams collected over seven generations (S1 to S7) from 2012-13 to 2018-19 were utilized in the study. In each generation, juvenile body weight at 0 day (BW0), 2 nd (BW2), 4 th (BW4) and 6 th (BW6) week of age; shank length at 6 th week (SL6); adult body weight at 20 th (BW20) and 40 th (BW40) weeks were measured. The body weight was measured to 0.1g accuracy using digital balance, while shank length was measured to the nearest of 0.01mm accuracy using digital Vernier Calipers. Age at sexual maturity (ASM); egg weight at 28 (EW28) and 40 (EW40) weeks; egg production (EP40) and egg mass (EM40) up to 40 weeks of age were recorded. The weight of eggs was recorded using a digital balance to an accuracy of 0.01g. The data were not recorded on BW0 in S-3 and S-7; BW2 in S-6 and hatch 8; BW40 in S-7 generation. The descriptive statistics of growth and production data are presented in Table 1.

Statistical analysis
Variance and covariance components were estimated by average information restricted maximum likelihood (AIREML) fitting an animal model [21]. Data were first analyzed by least squares analysis of variance (SPSS 12) to identify the fixed effects to be included in the model [22]. Two statistical models were used for identifying the significant effects on the traits. For BW0, BW2, BW4, BW6, SL6, ASM, BW20, BW40, EP40 and EM40, the statistical model included the fixed effect of generations (7 levels) and hatch number (8 levels). Only significant effects (P�0.05) were included in the models, which were subsequently used for the genetic analysis. Convergence was assumed when change of value of the natural logarithm of the likelihood function in two consecutive iterations was lower than 5× 10 −4 . Univariate animal models were fitted to estimate (co)variance components for all the traits. The six single-trait models Table 1. Characteristics of data on growth and production traits in the PD-3 line. used to estimate variance components were defined as follows: Where y is the vector of records; β, a, m, pe and ε are vectors of fixed, direct additive genetic, maternal additive genetic, permanent environmental effects of the dam, and residual effects, respectively; with association matrices X, Z a, Z m and Zpe ; A is the numerator relationship matrix between animals; and σ am is the covariance between additive direct and maternal genetic effects. Assumptions for variance (V) and covariance (Cov) matrices involving random effects were V(a) = Aσ 2 a , V(m) = Aσ 2 m, V(c) = Iσ 2 c , V(e) = Iσ 2 e , and Cov(a,m) = Aσ am, where I is an identity matrix and σ 2 a, σ 2 m, σ 2 c and σ 2 e are additive direct (a), additive maternal (m), maternal permanent environmental (c) and residual variances (e), respectively. The heritability (h 2 ) was calculated as h 2 = σ 2 a /σ 2 p, where σ 2 p = σ 2 a + σ 2 m + σ am + σ 2 c + σ 2 e . The m 2 , r am and c 2 were calculated as m 2 = σ 2 m/ σ 2 p , r am = σ am / σ 2 p c 2 = σ 2 c/ σ 2 p, respectively. The total heritability (h 2 t ), was calculated using the formula h 2 t = h 2 + 0.5m 2 + 1.5mr am h [23]. Akaike's Information Criterion (AIC) was used for selecting the best model among the tested models [24]. The model with lowest AIC value was chosen as most appropriate model and used to study the genetic parameters.
The best models from the single trait analyses were combined with appropriate (co)variance between random effects in the model for the bivariate analysis. The best model identified for a specific trait using likelihood ratio test was only used for the bivariate analysis. The estimates of genetic parameters i.e., genetic, phenotypic and environmental correlations between different economic traits were obtained by AIREML fitting an animal model. To test the significance of the genetic covariance, the full model was compared with the model in which genetic covariance was zero (COV A = 0). Significance of maternal permanent environmental covariance was also tested accordingly (COV C = 0). [25]. Significance of phenotypic correlations was tested by the hypothesis test to decide whether the value of the correlation coefficient was significantly different from zero [26].
The EBV obtained from the best single trait model suited for each trait was used to plot the genetic trend. The genetic trend was estimated by regression of the EBV of the females that contributed in each generation for the trait under selection (EM40) and the other production traits under study [26]. The EBV of contributed individuals was utilized for estimating the genetic trends of growth traits.

Growth and production traits
The least squares means (LSM) for body weight up to six weeks of age and SL6 are presented in Table 2. Generation and hatch had significant (P�0.01) effect on the body weight (BW0, BW2, BW4 and BW6) and shank length (SL6) in PD-3 line. The BW6 and SL6, the important juvenile traits showed an increasing trend over the generations. Hatch 1 recorded significantly (P�0.01) higher body weights and shank length. The interaction between the generation and hatch was significant (P�0.01). The LSM for interaction effects of growth traits is presented in Table 3. The overall average LSM for BW6 and SL6 were 273.93±0.62 g and 53.97±0.05 mm, respectively.
The LSM for ASM, BW20, BW40, EW28, EW40, EP40 and EM40 are presented in Table 4. All the production traits were significantly (P�0.01) influenced by the generation, hatch and their interaction. The average LSM for ASM, EP40 and EM40 were 168.82±0.25 d, 72.60±0.41 eggs and 4.21±0.07 kg, respectively, in PD-3 line. The ASM showed decreasing trend, whereas egg production and egg mass showed an increasing trend over the generations. The performance of various traits was higher in early hatches compared to later hatches ( Table 4). The LSM for interaction effects for production traits are presented in Table 5. The performance of different traits was higher in the later generations in early hatches in general with few exceptions.

(co)Variance components
The appropriate model for different economic traits was selected based on the AIC values obtained in WOMBAT [21]. The AIC values of models for traits analyzed are presented in Table 6. The estimates of (co) variance components with respect to additive, maternal permanent environmental and residual effects arrived by employing the best model for growth traits

PLOS ONE
(co) variance components and breeding value of economic traits in Dahlem Red chicken

PLOS ONE
(co) variance components and breeding value of economic traits in Dahlem Red chicken are presented in Table 7. Model 5 was the best for BW0, BW2 and BW4 and model 4 for BW6 and SL6. The estimates of various models are presented in S1 Table in S1 Material. The variance was partitioned into additive direct, maternal genetic, maternal permanent environmental and residual variance for growth traits up to 4 weeks of age. BW6 and SL6 had all the effects, except maternal genetic variance. The appropriate models for production traits were model 4 for ASM, EP40 and EM40; model 1 for EW28 and EW40; model 5 for BW20 and model 2 for BW40, respectively ( Table 8). The estimates of various models are presented in S2 Table in S1 Material. The ASM, EP40 and EM40 had additive direct, maternal permanent environmental and residual variance components, respectively. BW20 had the source of variation from additive direct, maternal genetic, maternal permanent environmental and residual components, while BW40 had similar sources, except for maternal permanent environmental variance. EW28 and EW40 had only additive genetic and residual components in Model 1.

Heritability
The heritability (h 2 ) estimates were low to moderate in magnitude for all the growth traits and ranged from 0.02±0.03 to 0.19 ±0.03 (Table 7). The trend of additive and maternal heritability estimates for growth and production traits over the age was depicted in Figs 1 and 2, respectively. The maternal heritability (m 2 ) at birth was high in magnitude (0.35±0.06), which reduced gradually by the end of the 4 th week (Fig 1). The c 2 was also higher at birth and gradually reduced with age. The ASM was moderately heritable (0.19±0.04). The heritability of EP40 and EM40 was low with 0.11±0.03 and 0.12±0.04 h 2 estimates, respectively. BW20 and BW40  Fig 2).

Correlations
The direct additive genetic correlation (r a ) between BW0 and other juvenile traits was less in magnitude and positive, while maternal genetic (r g ) and maternal permanent environmental (r c ) correlation were significantly higher; phenotypic correlation (r p ) was less in magnitude ( Table 9). The r g was observed till 4 weeks of age and r c continued up to 20 weeks of age. The r a between BW2, BW4, BW6 and SL6 was highly correlated with significant (P� 0.05) positive association. The r c was observed up to 20 weeks of age with significant high correlation coefficients between the body weights and shank length. The r p was higher among BW2, BW4, BW6 and SL6 with a positive association. The r a between BW6 and adult body weights (BW20 and BW40) was high in magnitude. The correlation coefficients from different components between production traits are presented in Table 10. The r a between ASM and BW20 was negative as heavier birds matured early. ASM and BW40 had a low degree of positive association. ASM had significantly (P� 0.05) higher negative correlation (r a , r c and r p ) with EP40 and EM40 (Table 10). The r a and r p between ASM and egg weights (EW28 and EW40) were of low magnitude and positive in Table 8. Estimates of (co)variance components and genetic parameters for production traits of PD-3 line.  direction. The r a between EP40 and BW40 was negative, while between EP40 and BW20 was positive with less magnitude. EP40 and EM40 had highly significant (P� 0.05) positive association with r a , r c and r p with. The EP40 and egg weight (EW28 and EW40) had negative association, while EM40 and egg weight had a positive association for both r a and r p . The EM40 had a positive association with all the production traits, except ASM. The adult body weight (BW20 and BW40) and egg weight (EW28 and EW40) had highly significant (P� 0.05) positive r a and low magnitude r p. The BW20 and BW40 had highly significant (p� 0.05) positive association for r a and r g and r p components.

Breeding value and genetic gain
The average EBV for BW6 and SL6, the important juvenile traits, was 17.65 g and 1.33 mm, respectively during the S-7 generation. The genetic trend of BW6 and SL6 showed significant (P < 0.05) and linear increment over the generations. The EBV of EM40 (0.48 kg) at the end of seventh generation significantly (P� 0.05) increased linearly as a response to selection in PD-3 line. The average genetic and phenotypic response was 0.08 and 0.11 kg per generation for EM40, the primary trait of selection (Figs 3 and 4). ASM decreased significantly with an increment of one day in each generation. The breeding value of EP40 showed an increasing trend with a genetic gain of 1.87 eggs per generation. The EW28 and EW40 increased linearly with an average genetic gain of 0.09 and 0.04 g, respectively per generation.

Inbreeding
The inbreeding was almost negligible in the population. The average inbreeding coefficient of the population was 0.019 and the average ΔF was 0.007 over the last seven generations of selection. The inbreeding coefficient of the population increased gradually from 0.00 to 0.04 in the 7 th generation. The trend of inbreeding and rate of inbreeding is presented in Fig 5.

Discussion
The effects of generation, hatch and their interaction were significant (P� 0.05) on body weight and shank length ( Table 2) with an increasing linear trend over the generations and

PLOS ONE
(co) variance components and breeding value of economic traits in Dahlem Red chicken decreasing trend among the hatches. Significant effect of generation and hatch on body weight and shank length similar to the present findings was reported in PD-1 line [4,12], in PB-2 line [17] and in White Leghorn layer strains [14]. The higher body weight in the early phase of life is important as it helps the bird to attain maturity at an early age and aid indirectly in increasing the egg production. Higher body weight in a female parent may contribute to the higher body weight in the terminal cross, which was desirable in a rural chicken variety. The LSM was significantly higher either in the first or second hatch, which may be due to the variations in environmental factors like maternal effects, hatching conditions and rolling reactions etc. The LSM of ASM significantly (P� 0.01) reduced over the generations as a correlated response to the selection for EM40 in PD-3 line. The increased body weight over the generations may have resulted in reduced ASM in a positive direction. Selection for body weight resulted in the reduction of ASM though direct selection was not practiced for this trait [27]. The significant (P� 0.01) increase in adult body weight (BW20 and BW40) might be attributed to the management, environment and feed restriction schedule followed during the grower phase to maintain the target body weight at laying. The body weight and egg production had a negative relationship in PD-3 line in the present study. The negative association between body weight and egg production in chicken was documented by many authors [4,7,11,20,28]. Higher early body weight is of great importance for the onset of egg production as heavier birds mature early and lay more number of eggs. The egg production (EP40) has shown a significant increasing trend with fluctuations across the generations (Table 4) as a correlated response to selection for higher egg mass as both are positively correlated traits. The EW28 showed a gradually increasing trend, while EW40 showed a fluctuating trend with the egg weight stabilized between 53 and 55 g. Selection for higher egg mass might be the possible reason for maintaining egg weight without any reduction, in spite of the increased EP40 as selection based on EM40 takes care of both egg number and egg weight [20]. The EM40 showed an increasing trend with exceptions during 2 nd and 5 th generation where lower egg mass was observed, which was due to the reduced egg production in the generation. However, a significant positive trend was observed for EM40 over the generations (Table 4).
Model 5 was appropriate for juvenile body weight (BW0, BW2 and BW4) with a, m, c, r and p effects and Model 4 for BW6 and SL6 with a, c, r and p effects. Maternal effects, both m and c were observed up to 20 weeks of age in PD-3 (Table 8). Jasouri et al. (2017) observed maternal genetic effects on body weight up to 12 weeks of age in dual purpose chicken [29]. The reduction in maternal effects as age increased was reported in broiler crosses [9] and commercial broilers [30]. In the present study, the maternal genetic effect continued till 40 weeks of age, contrary to the earlier findings (Table 8). Maternal environment effects are pre and post-ovipositional, the later may be pre and post hatch effects which influence egg quality traits and chick weight at hatch [31]. Maternal effects are due to the contribution of maternal additive, dominant gene effects and the differences between favorable and unfavorable allele frequencies [9]. The maternal effects were observed on growth traits up to six weeks of age in PD-1 line [4] and PB-2 broiler female line [17], contrary to the present findings. The continued existence of maternal effects in PD-3 line beyond juvenile phase might be due to the breed characteristic as Dahlem Red breed is one of the excellent female parent lines. The inclusion of one of the maternal effects either m, or am or c in the model could be enough to adjust the variation that occurred due to all the maternal effects [32]. The importance of "m" effects for early body weight was demonstrated, though the magnitude was less than that of "c" effects [7,8] in different chicken populations.
Model 4 was the appropriate model for ASM, EP40 and EM40 with effects of a and c, while model 1 was the best for EW28 and EW40 with direct additive genetic effect. Similar to the present findings, maternal effects in addition to the direct additive effect for the egg production and egg weight [29,33,34] were reported in chicken. However, some authors reported that model 1 without or with negligible maternal effects was suitable for egg production traits [8,35,36], contrary to the present study. A comprehensive review of the literature revealed that additive genetic effect increased and maternal effects decreased with age for production traits [37][38][39].
The estimates of h 2 determined using the animal model was more precise due to partitioning of the variance and covariance into all possible sources. The h 2 estimates of low magnitude with high precision and accuracy using the animal model were documented in the literature [5,7,8,29]. The variation in h 2 estimates over the time was due to the activation of various genes responsible for egg production [39] and persistency [40]. Higher h 2 estimates for body weight (0.24-0.47), ASM (0.36), egg number (0.17), egg mass (0.16) and egg weight (0.32-0.43) were observed in Mazandaran chicken [39] due to ignoring the maternal effects compared to the present estimates. The lower estimates in the present study might be due to the prolonged selection for higher egg mass, which reduced the variability in the population and also due to the inclusion of additional sources of variance in the model. Higher magnitude of h 2 estimates was also reported for egg production in Horro chicken from Ethiopia [41] in commercial layer line [42] and in Korean chicken [43].
The heritability estimated using the animal model was lesser as compared to the Henderson variance component analysis model. Maternal effects account for only small portion (2-8%) of the variance in chicken, but ignoring these effects results in biasness and over estimation of genetic parameters [35,36,44] leading to inadequate inferences hampering the breeding goals [6,34,45]. Similar to the findings of the present study, h 2 estimates for growth and egg production traits were reported in PD-1 [4]; in PB-2 [17] from India using the REML including the maternal effects in the model. The precise and accurate h 2 estimates in the present study were due to the inclusion of maternal effects in the model, which will aid the breeder in planning the breeding program with reliable prediction.
Genetic correlation plays a significant role in the success of the breeding experiment, as selection for one trait improves the performance in other traits as a positive correlated response, whereas negative association depresses the performance [7,11]. In addition, r a , r c continue to exist up to 20 weeks of age (Table 9) indicating the existence of prolonged maternal effects as PD-3 line is a female line, contrary to the findings in other lines wherein the maternal effects ceased at an early age [4,17]. The factors like hatch conditions, egg size and uterus have considerable influence on early body weight leading to the high effect of c, which reduces or becomes negligible at later stages of life [4]. However, in this study, effect of c prolonged up to 40 weeks of age. The r g between growth traits (Body weight and Shank length) was high in magnitude in PD-3 line in the present study, similar to the reports in PD-1 line [4] in PB-2 broiler line [17]; in Mazandaran chicken [39]; in Thai native chicken [46].
The correlation (r a , r c and r p ) between ASM and EP40 was significant with a high degree of negative association in PD-3 line. A similar significant inverse relationship between ASM and egg production in chicken was reported by many authors [4,11,17]. The correlation between ASM and egg weight was positive as the birds matured late laid heavier eggs. The association between ASM and BW20 was negative and was in the desirable direction as heavier birds matured earlier. A similar association between ASM and egg weight; and ASM and BW20 were reported by many authors [11,14]. The association (r a ) between BW40 and EP40 was negative ( Table 10). The inverse relationship between body weight and egg production was an established fact in chicken [7,11,28]. The EM40 had a direct positive association with EP40, EW28 and EW40, while had an antagonistic relationship with ASM. The direct positive association was justified as egg mass is the product of egg number and egg weight, and the higher number of eggs with better egg weight resulted in higher egg mass. The selection for egg mass maintains both egg number and egg weight at an optimum level, which is commonly practiced in breeder lines [20]. The early matured bird recorded higher egg mass with an inverse relationship as these birds produced more number of eggs with optimum egg weight leading to higher egg mass. The direction, magnitude and precision of correlation enable the breeder to fix the favourable traits in selection for improvement of the traits simultaneously [4]. Multi trait selection with adequate importance to each trait is the better option for optimal growth and egg production in chicken as these traits have antagonistic relationships making it more complex.
The average EBV of EM40, the primary trait of selection recorded a significant linear trend over the generations indicating the positive selection response in the population with significant improvement (Fig 3). The positive genetic trend was observed in BW6 and SL6 as correlated responses to selection, though the selection was not practiced for growth traits in this line. Similar observations of indirect response in correlated traits were reported in other chicken lines [4,17,20]. The direct positive genetic trend for BW6 and SL6 was observed in PD-6, a rural male parent line [47]. The EBV of ASM reduced linearly over the generations as an indirect response to selection for egg mass. The EBV of EP 40 showed a linear increasing trend with an average genetic gain of 1.87 eggs per generation (Fig 3), while the phenotypic trend showed fluctuations (Fig 4). The positive genetic and phenotypic gains in the population might be due to the effect of selection for egg mass, which improved both egg production and egg weight. The EW40 showed a linear increasing trend with fluctuations over the generations as selection for egg mass maintained the egg weight in the population. The EBV and the genetic trend observed in all the production traits was due to the correlated response to selection for egg mass as the direct selection was not practiced for those traits.
The ΔF in the population was very low (0.007), which may be due to the ideal breeding strategy followed in the PD-3 line, wherein close relatives were not included in the mating plan. The number of sires and dams utilized and contributed was higher than those required for maintaining effective population size. The inbreeding coefficient of the population at the end of the 7 th generation was 0.019. The inbreeding was not set in PD-3 chicken or negligible after 7 generations of selection. The delayed onset of inbreeding in pedigreed breeding populations similar to the present study was reported in pureline chicken [20,47]. The PD-3 population was in ideal condition with respect to genetic architecture after seven generations; however, the inbreeding coefficient may increase in further generations of selection.

Conclusion
The EBV of primary and associated traits was positive with a linear trend in the desired direction indicating the efficacy of selection practiced in PD-3 line with negligible inbreeding over the seven generations. Partitioning of variance into additive, maternal genetic, maternal genetic permanent environmental and residual sources improved the precision of the BV and genetic parameters. The accuracy in EBV and the genetic parameters aids in selecting the suitable breeding strategy for genetic improvement of PD-3 line. The study concluded that the positive genetic gains with respect to economic traits in PD-3 line lead to improvement in terminal crosses intended for backyard poultry farming.