ناهمگني اجزاي واريانس پروتئين شير در سطوح مختلف توليد گله- سال و تاثير آن بر پارامترهاي ژنتيكي و ارزش اصلاحي برآورد شده گاوهاي هلشتاين ايران

Heterogeneity of Variance Components for Milk Protein Yield in Different Levels of Herd-year Production and Its Effects on Genetic Parameters and Estimated Breeding Value of Iranian Holsteins

در اين مطالعه از ركوردهاي 305 روز مقدار پروتئين مربوط به سه شكم زايش اول گاوهاي هلشتاين ايران به منظور بررسي تاثير ناهمگني اجزاي واريانس بر پارامترهاي ژنتيكي، ارزش اصلاحي برآورد شده و تغيير رتبه گاوهاي نر و ماده ممتاز استفاده شد. داده ها طي سال هاي 1362-1392 توسط مركز اصلاح نژاد دام و بهبود توليدات دامي كشور جمع آوري شده بودند و براي دوره هاي شيردهي اول، دوم و سوم به ترتيب 141670، 115395 و 82529 ركورد وجود داشت. داده ها بر اساس سطوح مختلف توليد گله-سال به سه دسته بالا، متوسط و پايين گروه بندي شدند. ناهمگني واريانس ها با آزمون بارتلت بررسي شد و نتايج آن در هر سه دوره شيردهي معني دار (0/01˂p) بود كه نشان دهنده ناهمگن بودن واريانس ها بين سطوح توليدي است. به منظور برطرف نمودن احتمالي و يا كاهش ناهمگني واريانس ها از روش هاي تبديل جذري و باكس كاكس و همچنين روش پيش تصحيح ويشر استفاده شد. پارامترهاي ژنتيكي، وراثت پذيري ها و اجزاي واريانس با استفاده از مدل حيواني و نرم افزار VCE برآورد شدند. برآوردهاي وراثت پذيري در روش پيش تصحيحي ويشر نسبت به ساير روش ها بالاتر بود كه ممكن است به دليل تخمين صحيح تر اثرات ژنتيكي افزايشي در نتيجه در نظر گرفتن ناهمگني واريانس ها باشد. همبستگي هاي اسپيرمن و همچنين نسبت حيوانات انتخاب شده قبل و بعد از تبديل داده ها بر اساس ارزش هاي اصلاحي برآورد شدند. برآوردهاي همبستگي وقتي كه تمامي حيوانات نر و ماده ممتاز در نظر گرفته شدند بالا بودند. ناهمگني واريانس تاثير قابل ملاحظه اي بر تغيير رتبه و انتخاب 5%گاوهاي نر و 1% ماده هاي ممتاز داشت. روش هاي تبديل باكس-كاكس، جذري و پيش تصحيح ويشر باعث شدند به ترتيب 4%، 19% و 10% گاوهاي نر و 10، 21 و 7% ماده هاي برتر در مقايسه با سناريوي واريانس همگن از ليست انتخاب خارج شوند. نتايج مطالعه حاضر نشان داد كه واريانس ها بين سطوح مختلف توليد گله-سال همگن نبوده و ممكن است روي رتبه بندي و انتخاب گاوهاي ممتاز تاثير گذارد.

This study was carried out to investigate different data transformation methods on homogeneity and heterogeneity of variance components. Data included 305-day lactation records for protein yield from the first three lactations of Iranian Holstein cows collected from 1983 to 2014 by the Animal Breeding Center and Promotion of Animal Products of Iran. Data included 141670 records for 1st lactation, 115395 records for 2nd lactation and 82529 records for 3rd lactation. Records were categorized to 3 classes according to the average of herd-year production. For testing the heterogeneity of variance components Bartlet test was used and it was significant among all three lactations. A pre-correction method and two different data transformation methods including Box-Cox and Square root were used to correct for heterogeneity of variance. Genetic parameter and heritability estimates were estimated by VCE program, under an animal model. Spearman correlations and proportion of animals selected before and after data transformation were also estimated. Application of the Visscher adjustment method resulted in slightly higher heritabilities, which may be due to the more accurate estimation of additive genetic effects when heterogeneity is considered. Heterogeneity of variance had a significant effect on re-ranking and selection of 5% top sires and 1% of top dams. Pre-correction, Box-Cox and Square root method caused a proportion of 4%, 19% and 10% of top sires and 10%, 21% and 7% of top dams, respectively, to be excluded from selection when compared to the homogenous variance scenario. The results of this research indicate that the variance between different levels of herd-year production is not homogeneous and may influence the ranking and genetic evaluation of top cows