초록
본 연구는 축산연구소 한우시험장에서 출생한 한우 암소로부터 시간적인 간격을 두고 조사된 체중측정 기록에 대해 비선형의 성장곡선 모형을 적용하여 추정된 성장곡선 모수의 유전적인 경향을 평가하기 위해 실시하였다. 본 연구에서 성장곡선 모수들의 유전력 추정은 단형질 모형과 다형질 모형으로 분석하였으며 단형질 모형의 경우 선형모형은 출생년도-계절과 어미소의 나이의 효과가 포함된 동기우 집단을 고정효과로 상가적 개체유전효과를 임의효과로 하는 Model I과 Model I에 최종 체중측정시의 일령을 일차식 공변이로 추가시킨 Model II 등 두 가지 분석모형을, 그리고 다형질 모형의 경우 출생년도-계절과 어미소 나이의 효과를 고정효과로 하는 Model I과 Model I에 최종 측정시 일령을 공변이로 추가시킨 Model II 등 두 가지 분석모형을 이용하였는데, 단형질 모형의 Model I을 이용하여 추정된 성장곡선 모수 중 성숙체중의 유전력은 모형별로 0.09~0.22의 범위였으며, 성장비는 0.07~0.13의 범위였고, 성숙률은 0.05~0.07의 범위였다. 그리고 Model II를 이용하였을 때는 모형별로 성숙체중이 0.12~0.28, 성장비가 0.07~0.13의 범위였으며 성숙률은 0.12로 Gompertz 모형이나, Von Bertalanffy 모형 그리고 Logistic 모형이 모두 같았다. 한편 다형질 모형의 Model I을 이용하여 추정된 성장곡선모수 중 성숙체중의 유전력은 모형별로 0.09~0.17의 범위였으며, 성장비는 0.07~ 0.13의 범위였고, 성숙률은 0.06으로 세모형이 같았다. 그리고 Model II를 이용하였을 때는 성숙체중은 0.10~0.23, 성장비는 0.00~0.01, 성숙률은 0.06~0.11의 범위였다. 본 연구에서 추정된 성장곡선 모수들의 유전력은 외국의 육우에서 보고되는 유전력보다 낮았으며 한우수소에서 보고된 것과 유사한 결과였다. 그리고 Model II는 성숙체중과 성숙률의 유전력이 Model I보다 크게 추정되어 최종 측정시 일령을 공변이로 첨가할 경우 성숙체중과 성숙률의 상가적유전분산의 크기를 증가시키는 결과를 얻었다. 각 월령별 실측체중과 각 성장곡선 모형에 적합시켜 추정한 월령별 체중들에 대해서는 단형질모형을 이용하여 유전력을 추정하였는데 분석에 이용된 선형모형은 출생년도-계절과 어미소의 나이의 효과가 포함된 동기우 집단을 고정효과로 상가적 개체유전효과를 임의효과로 하는 Model I이었다. 실측체중의 경우 24개월령 체중만 0.52로 한우에 대한 타 연구자들의 결과에 비해 높았고 그 외의 월령별 체중은 타 연구자들의 결과 범위에 포함되는 성적이었다. 각 성장곡선모형으로 적합시켜 구한 생시체중의 유전력은 Gom- pertz 모형이 0.08, Von Bertalanffy 모형이 0.08 그리고 Logistic 모형이 0.06으로서 실측된 생시체중의 유전력 0.27에 비해 높게 나타났다. 그리고 실측체중의 경우 24개월령 체중의 유전력이 0.52, 36개월령 체중의 유전력이 0.32로서 36개월령의 유전력이 24개월령의 유전력에 비해 낮아지는데 적합체중의 경우에는 36개월령 체중의 유전력과 24개월령 체중의 유전력의 차이가 없거나(Gompertz 모형), 오히려 36개월령 체중이 24개월령 체중에 비해 유전력 추정치가 높아지고 있다(Von Bertalanffy 모형, Logistic 모형). 이렇게 적합체중에서 생시의 유전력이 낮아지거나 실측체중의 경우처럼 24개월령 체중보다 36개월령 체중의 유전력이 낮아지지 않는 것은 본 연구에 이용된 각 성장모형들이 생시체중을 실측체중보다 높게 추정하고 36개월령 체중을 낮게 추정하기 때문인 것으로 판단된다. 본 연구 결과로 볼 때 성장곡선 모형으로 추정된 월령별 체중들간에 유전력의 차이가 나타나 한우 암소의 성장예측을 위한 성장곡선의 사용은 중요하게 다루어져야 할 것으로 사료되며, 성장곡선 모수들에 대한 유전능력을 예측하여 한우 암소집단에 대한 선발과 도태의 기준으로 활용한다면 암소의 육용형 개량에 도움이 될 것으로 사료된다.
The objective of this study was to estimate genetic variances of growth curve parameters in Hanwoo cows. The data used in this study were records from 1,083 Hanwoo cows raised at Hanwoo Experiment Station, National Livestock Research Institute(NLRI). First evaluation model(Model I) fit year-season of birth and age of dam as fixed effects and second model(Model II) added age at the final weight as a linear covariate to Model I. Heritability estimates of A, b and k from Gompertz model were 0.22, 0.11 and 0.07 using modelⅠ and 0.28, 0.11 and 0.12 using modelⅡ. Those from Von Bertalanffy model were 0.22, 0.11 and 0.07 using modelⅠ, 0.28, 0.11 and 0.12 using modelⅡ. Heritability estimates of A, b and k from Logistic model were 0.14, 0.07 and 0.05 using modelⅠ, 0.18, 0.07 and 0.12 using modelⅡ. Heritability estimates of A from Gompertz model were higher than those from Von Bertalanffy model or Logistic model in both model Ⅰand model Ⅱ. Heritability estimates of b from Logistic model were higher than those from Gompertz model or Von Bertalanffy model in both modelⅠand model Ⅱ. Heritability estimates of birth weight, weaning weight, 3 month weight, 6 month weight, 9 month weight, 12 month weight, 18 month weight, 24 month weight, 36 month weight were after linear age adjustment 0.27, 0.11, 0.19, 0.14, 0.16, 0.23, 0.52 and 0.32, respectively. Heritability estimates of birth weight, weaning weight, 3 month weight, 6 month weight, 9 month weight and 24 month weight fit by Gompertz model were larger than those estimated from linearly adjusted data. Heritability estimates of 12 month weight, 18 month weight and 36 month weight fit by Von Bertalanffy model were larger than those estimated from linearly adjusted data. In the multitrait analyses for parameters from Gompertz model, genetic and phenotypic correlations between A and k parameters were -0.47 and -0.67 using modelⅠand -0.56 and -0.63 using model Ⅱ. Those between the A and b parameters were 0.69 and 0.34 using modelⅠand 0.72 and 0.37 using model Ⅱ. Those between the b and k parameters were -0.26 and 0.01 using modelⅠand -0.30 and 0.01 using model Ⅱ. In the multitrait analyses for parameters from Von Bertalanffy model, genetic and phenotypic correlations between A and k parameters were -0.49 and -0.67 suing model Ⅰ and -0.57 and -0.70 using modelⅡ. Those between the A and b parameters were 0.61 and 0.33 using modelⅠ and 0.60 and 0.30 using model Ⅱ. Those between the b and k parameters were -0.20 and 0.02 using modelⅠ and 0.16 and 0.00 using modelⅡ. In the multitrait analyses for parameters from Logistic model, genetic and phenotypic correlations between A and k parameters were -0.43 and -0.67 using model Ⅰ and -0.50 and -0.63 using modelⅡ. Those between the A and b parameters were 0.47 and 0.22 using modelⅠ and 0.38 and 0.24 using modelⅡ. Those between the b and k parameters were -0.09 and 0.02 using model Ⅰ and -0.02 and 0.13 using model Ⅱ.