한국동물자원과학회지 : 51(2) 97~104, 2009 J. Anim. Sci. & Technol. (Kor.) 51(2) 97~104, 2009 국내홀스타인종젖소의선형형질의점수제분석 최태정 * 조광현 * 이기환 * 상병찬 ** 농촌진흥청국립축산과학원 *, 충남대학교낙농학과 ** Analysis for Linear Type Classification Scheme on Holstein Cows in Korea Te Jeong Choi*, Kwang Hyun Cho*, Ki Hwan Lee* and Byeong Chan Sang** National Institute of Animal Science, R.D.A*, Dept. of Dairy Science Chungnam National University** ABSTRACT Complement of test standard, evaluation methods and models are needed to improve national competitiveness and to exchange superior genetic resources through the comparison of genetic evaluation score among nations in dairy cattle. Therefore, this study was conducted for the application of international standard to Korea considering domestic circumstance by changing linearclassification test score system of 50 classes which is currently used in Korea to system of 9 classes which is used in advanced nations of dairy. 15,230 of holstein cow linear type records with first parity records for the fifteen linear type and one total score from 2001 to 2006 and pedigree data which were collected by the Korean Animal Improvement Association were used in this study. Population classified by 9 levels was more normal distributed than 50 levels. Correlation coefficients between 50 and 9 score system showed over 0.98 by each classification scheme. Therefore, the 50 point system can be substituted with 9 point system due to their highly positive correlation. However, scores in all traits were still very contingent on classifier under the 9 point system (p<0.001), and F values between foot angle and front teat attachment showed high fluctuation depending on classifier. It means that subjective opinions of classifier would influence on linear type score as ever even if class scheme transformed to system of 9 class. Therefore, the relevance of transformation to the 9 point system should be assessed after analyses about various environmental factors. (Key words : Linear type traits, Classification scheme, Conformation, Holstein cow). 서론 최근젖소산유능력은과거 10 년간유량연간개량량이 182.9 kg 으로괄목할만한증가를하였으나 ( 젖소유전능력평가보고서, 2007), 경제적수명에도달하지못한채도태되고있는실정으로고능력우가장수하여생애산유량을증대시키는것에대한관심이점차높아지고있다. 따라서, 젖소의경제수명을늘려생애산유량을증대시키기위해서는젖소의신체적기능과밀접한연관을갖고있는기능적체형형질의개량이무엇보다도중요하다고할수있다 (Brotherstone, 1991; 이등, 2006). 이러한이유로현재젖소에대한선형형질의심사는이미여러낙농국가들에서실시되고있으며그평가기록은장수성 (Brotherstone and Hill, 1991; Vollema and Groen, 1997; Weigel 등, 1998; Setati 등, 2004; Tsuruta 등, 2005; Neuenschwander 등, 2005; Vacek 등, 2006), 체중 (Veerkamp and Brotherstone, 1997; Koenen and Groen, 1998), 유방의건강 (Thomas 등, 1984; De Jong and Lansbergen, 1996), 발굽관련문제 (Boelling and Pollott, 1998), 그리고, 분만의용이성 (Dadati 등, 1985; Cue 등, 1990) 등을예측하는데중요하게이용되어져왔다. 국내에서도선형형질에대한자료의수집및심사와유전능력평가가꾸준히이루어지고있으며국내젖소의육종계획의수립시에이를활용하고있다. 우리나라의젖소검정두수는 2006 년말기준으로 20 만두에다가서고있으며농가와유관기관의꾸준한혈통관리로혈통비율도 2006 년 50% 를상회하여유전능력평가를하는데정확도측면에서꾸준히증가해왔다. 그러나국제간의유전평가인 MACE (Multiply Across Country Evaluation) 를통한우리의우수유전자원이교류되기까지는검정기준, 평가방법과모델등에있어많은보완이필요하다. 또한, 이를통하여경쟁력있는젖소를개량하기위해서는낙농선진국에서채택하고있는국제기준방법들을현실에맞게적용 Corresponding author : Kwang Hyun Cho, National Institute of Animal Science, San 9, Eryong-ri, Seonghwan-eup, Cheonan-si, Chungcheongnam-do, 330-801, Ref. of Korea. Tel: +82-41-580-3358, Fax: +82-41-580-3362, E-mail: ckh1219@rda.go.kr - 97 -
해야하는노력이필요하다. 일본의경우 2000 년에체형의국제간유전평가의정확도향상을위해선형식체형심사의평가수치를 1~9 단계로변경하여적용하고있으며 ( 일본가축개량사업국, 2002), 캐나다, 덴마크, 뉴질랜드등의낙농선진국이속한 WHFF (World Holstein Friesian Federation) 에서도 9 단계의선형심사에의한체형심사가실시되고있다. 이러한 9 단계선형심사의도입은선형형질에관한유전평가치의국제간비교시에그정확도를상당히향상시킬수있을것이라고기대된다. 따라서, 우리나라도국제가축기록위원회 (International Committee for Animal Recording, ICAR) 에서추천하고있는국제표준평가방법 (e.g., Hewitt, 2000) 에맞추어 9 단계로변경하는것을적극검토해봐야한다. 이를위해본연구에서는점수변경시에심사자에따른선형형질에대한심사점수의변이의차이를파악하고분산분석및상관분석을통해타당성을검토하여점수제변환을국내현실에맞게적용하기위해실시하였다. 또한, 50 단계심사점수와변환된 9 단계의선형심사점수에영향을미치는심사자, 심사시의연령그리고, 비유단계와같은환경효과에대한분석과검토를추가적으로수행하였다. 1.. 재료및방법 본연구에는한국종축개량협회 (KAIA) 가보유하고있는젖소에대한혈통기록과 2001 년부터 2006 년도까지실시한선형심사기록중최종점수 (final scores) 및 15 개선형형질에대해초산차기록을갖는 150,230 두의암소선형심사기록자료가분석에이용되었다. 선형심사형질들은심사자에의해 1~50 점수단계의범위에서점수화되었으며점수부여기준은미국홀스타인협회의표준평가방법 (Thompson 등, 1981; Holstein Association USA, 2005) 이적용되었다. 선형형질은현행의 50 점제의구분에서국제가축기록위원회 (ICAR) 에서채택하고있는국제표준평가방법을근거로 Table 1 과같이 9 점제로변환하였다. 그리고, 심사자는 4 명, 심사연령 (age at classification) 은 4 단계로, 비유단계 (lactation stage) 는 12 단계로구분하여나누어선형통계분석모형에적용하였다. 2. (1) 자료교정유효자료의선정을위해젖소의생년월일, 심사일, 등록번호, 심사일당시개체의연령, 심사자및선형형질별점수에대한이상치를제거한뒤한개체가 2 개이상의자료를가졌을경우에는마지막심사일에해당하는자료를사용하였으며, 심사자의심사빈도가매우낮고심사자를알수없는자료를제거하였다. 농가, 심사년도, 심사자별 Table 1. A calculating table between 50 and 9 point system 50 point system 9 point system 1 ~ 7 8 ~ 12 13 ~ 17 18 ~ 22 23 ~ 27 28 ~ 32 33 ~ 37 38 ~ 42 43 ~ 50 관측수가 5 이상인개체의심사기록만을이용하였고, 심사시의연령 (age at classification) 은월령으로바꾼뒤정규성검정에따라 15 개월령미만과 80 개월령을초과하는성적은분석에서제외하였고, 2 세이하에서 5 세이상까지총 4 단계로나누어분석에이용하였다. 분만후비유단계 (lactation of stage) 는 1 개월단위로총 12 수준으로분류하여분석모형에적용하였다. (2) 환경요인의효과분석심사점수에대한심사자, 심사시의연령및비유단계의요인의효과를알아보기위해다변량분산분석방법 (multivariate analysis of variance) 이용하여분석을실시하였으며분산분석을위한선형모형은 과같다. 분석은 SAS package 9.1 버전을이용하였고, PROC GLM 절차에서제공되는 4 가지제곱합중 TYPE 제곱합을이용하였다. 최소자승평균치 (LSM) 간유의성검정을위해다음과같은귀무가설을유의수준 5% 로각형질에대한검정을실시하였다. H 0 : LSM(i) = LSM(j) 여기서, LSM(i(j)) : i(j) 번째효과의최소자승평균치 (i j) 이다. Y ijk = μ + classifier i + age j + stage k + e ijk 여기서, Y ijk = i번째심사자의 j번째연령에대한 k번째비유단계의선형심사점수 μ = 전체평균 hyc i = i번째심사자의고정효과 : 4개 age j = j번째연령의고정효과 : 4개 stage k = k번째비유단계의고정효과 : 12개 e ijk = 각개체에미치는잔차효과 1 2 3 4 5 6 7 8 9-98 -
50 점제에서 9 점제로점수변환한후두가지점수간의 상관분석은다음의공식에의거하여 SAS package 의 PROC CORR 절차를이용하여분석하였다. 1.. 결과 Table 2 의자료는 50 점제심사점수를 9 점제로변환한 자료에대한기초통계량과유전력을나타낸것이다. 변환 된점수는 1 에서 9 점의범위를나타냈으며유방의깊이 (UD) 의경우표준편차가 1.63 으로다른형질들에비해높 게나타나자료의변이가크게나타났다. 체격을나타내 는일반외모형질중키 (ST) 의점수가 5.48 ± 1.56 로가장 높게평가되었고, 반면강건성은 3.53 ± 1.29 로가장낮은 평가점수를보였다. 유방관련형질중에서는유방의깊이 (UD) 가 4.13 ± 1.63 로가장높은점수를보였으며, 앞유두 의배열위치 (FTP) 가 3.60 ± 1.35 로가장낮은평가를받은 것으로나타났다. 유전력 (heritability) 에있어서는키, 엉덩 이기울기 (RA), 예각성 (DF), 옆에서본뒷다리 (RLSV) 가 각각 0.42, 0.37, 0.33, 0.32 로서비교적높은유전력을보 이고있고유방관련형질들또한, 0.23 0.29 의높은유전력을보이고있어중요한형질임을알수있었다. 반면발굽기울기 (FA) 는모든형질들중에서 0.15로서가장낮은유전력을나타냈다. 각형질별 50점제집단의분포특성과변환시킨 9점제집단의분포특성은 Fig. 1과같이키 (stature) 와옆에서본뒷다리 (RLSV) 등이 50단계분포보다변환된 9단계자료의분포가더욱정규성을띄고있었으며, 50단계에서는점수대별로분포의폭이작게나타났다. 이는 50단계의심사제도에서심사자간에선호하는점수대의점수편차가크게나타나기때문이었으며, 결국특정선호점수에집중되고그전후의점수의빈도가극단적으로감소하는경향을보이는것으로나타났다. 이결과는 50단계를사용한선형심사에서선형점수의모든부분을유용하게이용하고있지않다는것을보여주는것이다. 비록점수대구분의추정값은세분화할수록더정확하게나오는이치로점수폭을줄이는것이문제점이될수있으나선호심사대에대한분포가아닌정규성을갖는분포로만들기위해서는 9단계로전환이바람직하다고판단된다. 이로써 9단계로의변경은 50단계로심사해온방법에서경계에위치한기준의주관적인개입이감소할것으로예측되며, 9단계점수제를이용함으로써정규분포에근거한폭넓고공정한평가가이루어져유전력의저하를해결하는기능을가질것으로사료된다. 결과적으로젖소에대한선형심사시심사자의넓은점수영역대로인한부담을감소시켜효율적인선형심사가이루어질것으로판단된다. Table 2. Basic statistic and heritability for linear traits transformed by 9 class Traits N Mean ± SD Heritability* Stature (ST) 149,811 5.48 ± 1.56 0.42 Strength (STR) 147,849 3.53 ± 1.29 0.31 Rump angle (RA) 148,945 3.91 ± 1.35 0.37 Thurl width (TW) 149,645 3.88 ± 1.27 0.29 Dairy form (DF) 149,165 4.47 ± 1.52 0.33 Body depth (BD) 149,414 3.97 ± 1.35 0.26 Rear legs side view (RLSV) 150,230 4.86 ± 1.29 0.32 Foot angle (FA) 146,486 3.71 ± 1.33 0.15 Fore udder attachment (FUA) 145,954 3.97 ± 1.48 0.29 Rear udder height (RUH) 147,882 3.82 ± 1.54 0.28 Rear udder width (RUW) 147,831 3.68 ± 1.46 0.23 Udder cleft (UC) 148,294 4.07 ± 1.32 0.24 Udder depth (UD) 146,989 4.13 ± 1.63 0.28 Front teat placement (FTP) 146,708 3.60 ± 1.35 0.26 Front teat length (FTL) 144,530 4.06 ± 1.53 0.26 * Holstein genetic evaluation report in Korea at the second half of 2006. - 99 -
Fig. 1. Distribution of 9 class and 50 class data for linear type traits. 2. Table 3. Analysis of variance for linear traits and final scores by 50 class Source df ST STR BD MS F MS F MS F Cla. 1 3 103269.22 1928.74* 16350.83 472.51* 2622.67 69.59* Age 2 3 167998.41 3137.67* 230423.07 6658.78* 286128.65 7591.69* Stage 3 11 495.88 9.26* 1562.65 45.16* 1764.12 46.81* DF RA TW Cla. 3 60566.54 1778.59* 644.21 9.91* 25822.74 664.72* Age 3 71078.19 2087.28* 846.42 13.02* 246397.57 6342.64* Stage 11 12442.39 365.38* 3718.57 57.19* 784.31 20.19* RLSV FA FUA Cla. 3 7488.36 191.14* 123214.32 2832.71* 100075.41 1652.46* Age 3 40443.90 1032.35* 51890.39 1192.96* 4389.23 72.48* Stage 11 582.66 14.87* 367.55 8.45* 1292.10 21.34* RUH RUW UC Cla. 3 98160.23 1639.87* 80283.59 1798.60* 59437.41 1337.78* Age 3 6065.73 101.33* 121840.43 2729.60* 19678.60 442.92* Stage 11 17298.18 288.98* 18229.15 408.39* 3617.34 81.42* UD FTP FTL Cla. 3 10554.15 196.96* 8602.26 166.80* 12996.47 205.40* Age 3 504054.16 9406.59* 3919.50 76.00* 67174.63 1061.64* Stage 11 6162.68 115.01* 5968.15 115.72* 2538.06 40.11* MS = mean square, Cla. 1 = classifier, Age 2 = age at classification, Stage 3 = lactation stage, * = P < 0.001. - 100-
Table 4. Analysis of variance for linear traits and final scores by 9 class Source df ST STR BD MS F MS F MS F Cla. 3 3482.78 1724.18* 490.30 379.36* 102.04 70.34* age 3 6106.29 3022.98* 8290.68 6414.67* 10444.09 7199.24* stage 11 15.57 7.71* 37.65 29.13* 57.03 39.31* DF RA TW Cla. 3 1902.15 1363.61* 65.90 25.51* 968.65 638.49* age 3 2817.26 2019.63* 42.27 16.36* 9047.54 5963.77* stage 11 477.15 342.06* 152.37 58.98* 30.02 19.79* RLSV FA FUA Cla. 3 229.04 137.97* 3943.13 2388.85* 3466.51 1615.21* age 3 1747.80 1052.87* 1726.50 1045.96* 40.28 18.77* stage 11 23.91 14.4* 11.37 6.89* 59.05 27.51* RUH RUW UC Cla. 3 2926.36 1270.27* 2496.00 1460.84* 1972.15 1148.63* age 3 187.81 81.52* 4308.20 2521.47* 647.43 377.08* stage 11 649.76 282.05* 662.28 387.61* 115.16 67.08* UD FTP FTL Cla. 3 527.83 259.25* 430.54 229.92* 681.59 295.71* age 3 17033.44 8366.32* 64.87 34.64* 2218.69 962.59* stage 11 205.59 100.98* 189.62 101.26* 73.96 32.09* MS = mean square, Cla. 1 = classifier, Age 2 = age at classification, Stage 3 = lactation stage, * = P < 0.001. 50 단계의형질별점수집단과변환한 9 단계의점수집단들의분산분석결과를 Table 3 과 Table 4 에제시하였다. 두점수제의모든선형형질들이심사자, 심사시연령, 비유단계의환경요인과고도의유의성 (p<0.001) 을보이는것으로나타났다. 심사자에따라변이가크게나타나는형질은발굽기울기 (FA), 키 (ST), 앞유방의붙음성 (FUA), 예각성 (DF) 의순이었으며, 변이차이가적게나타나는형질은 엉덩이기울기 (RA), 체심 (BD), 유방의깊이 (UD) 의순이었다. 그리고, 심사시의연령에의해영향을가장많이받는형질은유방의깊이 (UD) 였는데이는송등 (2002) 의연구보고와같은결과였다. 반면가장적게영향을받는형질은엉덩이기울기 (RA) 였다. 또한, 비유단계에따라영향을많이받는형질은뒷유방의너비 (RUW) 와예각성 (DF) 이었으며이는이등 (2006) 의연구보고와같은결과였다. Table 5. Least square means and standard errors for linear traits score according to classifier by 9 class Trait ST STR BD DF RA TW RLSV Cla.* LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE 1 5.23±0.03 3.14±0.02 3.39±0.02 3.21±0.02 4.21±0.03 3.54±0.02 4.61±0.03 2 5.50±0.03 3.44±0.02 3.34±0.03 3.71±0.02 4.11±0.03 3.91±0.03 4.82±0.03 3 5.27±0.03 3.19±0.02 3.22±0.03 2.89±0.02 4.19±0.03 3.42±0.03 4.58±0.03 4 4.31±0.03 3.05±0.02 3.26±0.03 3.23±0.02 4.09±0.03 3.42±0.03 4.61±0.03 Trait FA FUA RUH RUW UC UD FTP FTL Cla.* LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE LSM±SE 1 3.77±0.03 3.43±0.03 3.23±0.03 2.60±0.03 4.16±0.03 4.81±0.03 3.70±0.03 4.01±0.03 2 4.46±0.03 4.09±0.03 4.06±0.03 3.22±0.03 3.81±0.03 4.50±0.03 3.46±0.03 3.79±0.03 3 3.69±0.03 3.86±0.03 3.56±0.03 3.00±0.03 3.83±0.03 4.77±0.03 3.79±0.03 3.75±0.03 4 3.26±0.03 4.44±0.03 3.40±0.03 3.40±0.03 4.65±0.03 4.90±0.03 3.81±0.03 3.53±0.03 * Cla. = Classifier, LSM = Least square mean. - 101-
집단은거의유사한분산의비를보였다. 하지만발굽기울기 (FA) 는 50 단계의점수를 9 단계로전환을해도다른형질들에비해상대적으로심사자에따른차이가많이존재한다는것으로나타났다. Table 6 은 9 단계로전환시킨자료집단과 50 단계집단의자료들에대한심사자, 심사시연령, 비유단계에따른상관분석을실시한결과이다. 분석결과모든형질에서심사자, 심사시의연령, 비유단계에따른상관계수가 0.97 0.99 의범위를보여세분화된집단단계에서그단계를줄여도점수제간에거의변화가없는것으로나타났다.. 고찰 Fig. 2. F value for linear type traits of 9 class and 50 class according to classifier. ST : stature, STR : strength, BD : body depth, DF : dairy form, RA : rump angle, TW : thurl width, RLSV : rear legs side view, FA : foot angel, RUH : rear udder height, RUW : rear udder width, UC : udder cleft, UD : udder depth, FTP : front teat placement, FTL : front teat length. 반면가장영향을적게받는형질은발굽기울기 (FA) 형질이었다. 이는체형형질들이연령과비유단계에유의차를보였다는것을보고한이등 (1994), 이등 (2007), Hayes 와 Mao (1987) 의연구보고와같은결과였다. 3. Table 5 는점수제변환후에도심사자의효과가크게작용하는것으로분석되어심사자에따른추가분석을실시한최소자승평균추정값 (Least squares mean, LSM) 을나타낸표이다. 체심 (BD), 엉덩이기울기 (RA), 앞유두의길이 (FTL) 의경우 1 번째심사자에서높게나타났고, 키 (ST), 강건성 (STR) 과엉덩이나발굽쪽의체형형질에대해서는 2 번째심사자에서높게분석되었다. 예각성 (Angularity/ DF) 은 3 번째심사자에서만높게나타났으며, 특히유방관련형질들은 4 번째심사자의경우에가장높은점수를주는것으로나타났다. Fig. 2 에는형질별 9 단계집단과 50 단계집단의심사자에따른 F 값의차이를나타내었다. 그림에서와같이두 2009 년기준농가의평균도태산차는 2.5 산으로이는우군내에서경제적수명을채우지못하고도태되는소가많다고할수있다. 경제적산차에도달하기위해서는 4 산이상까지갈수있도록번식, 사양, 환경등이조화를이루어야한다. 후대검정우들의선형심사를정확하게심사해야할필요도있지만집단의변화가거의없으며종모우의순위변화가없고집단간의상관관계가높게나타나므로연구의목적에부합하게정규분포에가까운분포와같이폭넓게평가하고유전평가의정확도를높이는방법으로전환해야할것이다. 결론적으로정확한선형심사를바탕으로젖소에대한유전능력평가의정확도를향상시키기위해서는현행의 50 단계의선형심사점수제에서국제표준평가방법인 9 단계로의변환을고려해야할것으로보인다. 본연구결과에서와같이심사자에따라평균점수및편차가다른것으로나타나 9 단계로의변경은가능하지만사전에다른효과들과함께육종가의추정시에모형내보정이나사전보정을해야할것으로판단되었다 (Veerkamp 등, 2002). 즉, 단기간의변경보다는심사자, 연령그리고, 비유단계와같은환경적요인등에대한다양한검토와충분한검정기간을거쳐적합성여부를결정해야할것이다. 그리고, 심사자에따른평가점수의편차를줄이기위해정기적인교육이반드시필요하며, 이를통하여심사척도의표준화가이루어진다면점수제변경후에도심사의정확도가크게향상될것으로사료되었다. 또한, 선형형질에대한점수제의변경이나심사자에따른유전모수추정과같은연구가전무하므로이에대한연구가추가적으로이루어져야할것으로사료된다. Table 6. Correlation analysis for linear type traits according to environmental factors between 9 class and 50 class Traits ST STR BD DF RA TW RLSV FA FUA RUH RUW UC UD FTP TL Classifier 0.98 0.98 0.98 0.98 0.99 0.98 0.98 0.98 0.99 0.99 0.98 0.98 0.99 0.99 0.99 Age* 0.98 0.97 0.98 0.98 0.97 0.98 0.98 0.98 0.99 0.99 0.98 0.98 0.98 0.99 0.99 Lactation stage 0.98 0.98 0.98 0.98 0.97 0.98 0.98 0.98 0.99 0.99 0.98 0.98 0.98 0.99 0.99 * Age at classification. - 102-
. 요약국제간의젖소유전능력평가점수의비교를통한우리의우수유전자원의교류와경쟁력향상을위해서는검정기준, 평가방법과모델등에많은보완이필요하다. 따라서, 본연구는현재국내의 50단계의심사점수를낙농선진국에서채택하고있는 9점제로의변환을통해국제기준을국내현실에맞게적용하기위해실시했다. 분석에는한국종축개량협회 (KAIA) 의젖소에대한혈통기록및 2001에서 2006년도사이에실시된체형측정기록중최종점수와 15개체형형질에대해초산차의기록을갖는 32,487두의암소집단의선형심사기록자료가이용되었다. 9점제로변경한집단이더욱정규성을띄는분포특성을나타내었다. 상관분석결과에서도모든형질에서 50점제와 9점제심사점수간에 0.98 이상의상관계수를나타내었다. 따라서, 50 점제에서 9점제로의변경이가능할것으로판단된다. 하지만, 9점제로의변환후모든형질에서심사자에따라고도의유의성을보였고 (P<0.001), 점수제간 F값은발굽기울기와앞유방의붙음성은심사자간변이가비교적높게나타났다. 이것은 9점제로의변환후에도심사자의주관적인요소가점수에영향을미치는것을의미한다. 따라서, 선형형질점수제변환은환경효과에대한다양한분석을통하여적합성여부를결정해야할것이다.. 인용문헌 1. Boelling, D. and Pollott, G. E. 1998. Locomotion, lameness, hoof and leg traits in cattle: II. Genetic relationships and breeding value. Livest. Prod. Sci. 54:205-215. 2. Brotherstone, S. and Hill, W. G. 1991. Dairy herd life in relation to linear type traits and production. e. Genetic analyses for pedigree and non-pedigree cows. Anim. Prod. 53:298-297. 3. Dadati, E., Kennedy, B. W. and Burnside, E. B. 1985. Relationships between conformation and reproduction in Holstein cows: type and calving performance. J. Dairy Sci. 68:2639-2645. 4. De jong, G. and Lansbergen, L. 1996. Udder health index: selection for mastitis resistance. Proceedings of the International workshop on genetic improvement of functional traits in cattle, Gembloux, Belgium, January 1996. INTERBULL Bull. 12:42-47. 5. Hayes, A. E. and Mao, I. L. 1987. Effect of parity, age and stage of lactation at classification on linear type scores of Holstein cattle. J. Dairy Sci. 70:1898. 6. Hewitt, D. 2000. Type harmonisation working group-report of World Holstein Friesian Federation 5th Classifiers Workshop. Pages 172-177 in 10th World Holstein-Friesian Conference, Sydney. 7. Holstein Association USA. 2005. http//www.holsteinusa.com/ Dec. 12, 2005 accessed. 8. Koenen, E. P. C. and Groen, A. F. 1998. Genetic evaluation of body weight of lactating Holstein heifers using body measurements and conformation traits. J. Dairy Sci. 81:1709-1713. 9. Setati, M. M., Norris, D., Banga, C. B. and Benyi, K. 2004. Relationships between Longevity and Linear Type Traits in Holstein Cattle Population of Southern Africa. Tropical Animal Health and Production, 36(8):807-814. 10. Thompson, J. R., Freeman, A. E., Wilson, D. J., Chapin, C. A., Breger, P. J. and Kuck, A. 1981. Evaluation of a Linear type program in Holsteins. 11. Tomas, C. L., Vinson, W. E., Pearson, R. E., Dicknson, F. N. and Johnson, L. P. 1984. Relationships between linear type scores, objective type measures, and indicators of mastitis. J. Dairy Sci. 67:1281-1292. 12. Neuenschwander, T., Kadarmideem, H. N., Wegmann, S. and de Haas, Y. 2005. Genetics of Parity-Dependant Production Increase and its Relationship with Health, Fertility, Longevity, and Conformation in Swiss Holsteins. J. Dairy. Sci. 88:1540-1551. 13. Tsuruta, S., Misztal, I. and Lawlor, T. J. 2005. Changing Definition of Productive Life in US Holsteins : Effection Genetic Correlations. J. Dairy Sci. 88:1156-1165. 14. Vacek, M., Stipkova, M., Nemcova, E. and Bouska, J. 2006. Relationships between conformation traits and longevity of Holstein cows in the CzechRepublic Czech. J. Anim. Sci. 51(8):327-333. 15. Vollema, A. R. and Groen, A. F. 1997. Genetic correlations between longevity and conformation traits in an upgrading dairy cattle population. J. Dairy Sci. 80:3006-3014. 16. Veerkamp, R. F. and Brotherstone, S. 1997. Genetic correlations between linear type traits, food intake, live weight and condition score in Holstein Friesian dairy cattle. Anim. Sci. 64:358-392. 17. Veerkamp, R. F., Gerritsen, C. L. M., Koenen, E. P. C., Hamoen, A. and Ee Jong, G. 2002. J. Dairy Sci. 85:976-983. 18. Weigel, K. A., Lawlor, T. J., JR., Vanraden, P. M. and Wiggans, G. R. 1998. Use of Linear Type and Production Data to Supplement Early Predicted Transmitting Abilities for Productive Life. J. Dairy Sci. 81:2040-2044. 19. 송치은, 상병찬, 도창희. 2002. 국내홀스타인젖소의선형심 - 103-
사에대한보정계수개발. 한국동물자원과학회지. 44(1):1-12. 20. 이기환. 2008. Holstein종젖소의선형형질에대한유전모수추정에관한연구. 충남대학교석사학위논문. 21. 이광전, 강민식, 조광현, 강민구, 이상구. 1994. 체형형질에대한보정계수개발. 한축지. 36(6):574-583. 22. 이득환, 김은길. 2006. Holstein 젖소의선형심사형질과등급형질에대한환경효과분석. 한국동물자원과학회지. 48(6): 767-776. 23. ( 사 ) 가축개량사업국. 2002. 유용우생애생산성향상기술. 일본 24. 농촌진흥청축산과학원. 2007. 2006년하반기젖소유전능력평가보고서 II. ( 접수일자 2008. 8. 7. / 수정일자 : 2009. 3. 26. / 채택일자 2009. 4. 9.) - 104-