ارزيابي ژنتيكي توليد شير در گاوهاي شيري نژاد هلشتاين استان خراسان با استفاده از مدل تابعيت تصادفي تكه‌اي

Genetic Evaluation of Milk Production in Holstein Dairy Cattle of Khorasan Province Using a Spline Random Regression Model

در اين تحقيق، تخمين پارامترها و ارزيابي ژنتيكي گاو ‌هاي هلشتاين استان خراسان براي صفت شير بر اساس يك مدل دام روز آزمون با تابعيت تصادفي تكه‌اي ( Cubic spline random regression ) انجام شد. داده‌هاي مورد استفاده 32854 ركورد روز آزمون شير دو و سه بار دوشش در روز مربوط به 3842 رأس گاو هلشتاين شكم اول (فرزندان 466 رأس گاو نر) در 125 گله بود كه در بين سال‌هاي 1380 تا 1384 زايش داشتند. در مدل دام روز - آزمون، اثر عوامل ثابت محيطي گله- سال و ماه ركوردگيري- دفعات دوشش، متغير‌هاي كمكي خطي و درجه دوم سن زايش اول، درصد ژن هلشتاين (متغير كمكي خطي)، اثرات تصادفي ژنتيكي افزايشي و محيط دائمي گاو قرار داده شد. در مدل مزبور به‌منظور در نظر گرفتن تفاوت‌هاي ژنتيكي و محيطي شكل منحني شيردهي بين گاوها از تابعيت تصادفي تكه‌اي توان سوم استفاده شد. ميانگين پسين ارزش اصلاحي پيش بيني شده گاوها با روش آماري بيزين (از طريق نمونه‌گيري گيبس با 100000 زنجيره با نرم افزار RRGIBBS ) محاسبه شد. ميانگين ارزش اصلاحي پيش بيني شده شير 305 روز براي گاوهاي داراي ركورد برابر با 52/90 كيلوگرم (0/05>P) بود. هم‌بستگي بين ارزش اصلاحي ماه‌هاي شيردهي با افزايش فاصله بين آنها كاهش نشان داد . بيشترين هم‌بستگي رتبه‌اي بين ماه‌هاي هشتم و نهم (0/998) و كمترين هم‌بستگي رتبه‌اي بين ماه‌هاي اول و دهم (0/312) بود. ارزش اصلاحي ماه‌هاي اول و ششم به‌ ترتيب كمترين (0/553) و بيشترين (0/990) هم‌بستگي رتبه‌اي را با ارزش اصلاحي شير 305 روز داشت. روند ژنتيكي برآورد شده شير 305 روز بر حسب سال تولد گاوها 17/75 كيلوگرم در سال بود كه به لحاظ آماري تفاوت معني‌دار از صفر نداشت ( 0/165= p value).

This study was undertaken to estimate genetic trend and parameters of Holstein cattle in Khorasan province for milk yield using a spline random regression test day animal model. A total of 32854 monthly test day milk records (twice and thrice a day milking) obtained from 3842 Holstein heifers (progeny of 466 sires) distributed in 125 herds and calved from 2001 to 2005 was used to predict breeding value of individual animals. In the model, fixed effects of herd including year-month of recording, milking times, age at calving (linear and quadratic covariables), Holstein gene percentage (linear covariable) as well as random effects of additive genetic and permanent environment were studied. To take account of the shape of the lactation curve at genetic and environmental levels, cubic spline polynomials were also included in the test day model. Bayesian method by applying Gibbs sampling technique (100000 chains applying RRGIBBS software) was utilized to obtain posterior means of predicted breeding value of animals for milk yield at individual month of lactation. The results showed that mean of breeding value for 305-day milk yield was 52.90 kg (p<0.05). Spearman rank correlations between predicted breeding values at different months of lactation decreased as the interval between them increased. The highest and lowest rank correlations were found between months 8 and 9 (0.998) and between months 1 and 10 (0.312), respectively. Predicted breeding value of 305-day milk had the lowest and highest rank correlations with predicted breeding value at months 1 (0.553) and 6 (0.990), respectively. Regression analysis of average predicted breeding value of progenies in their birth year showed that the amount of genetic trend for 305 day milk yield was 17.75 kg per year, statistically no different from zero (p value=0.165).