تجزیه و تحلیل ماندگاری در بره های نژاد ماکویی با استفاده از توابع خطر نسبی

نوع مقاله : مقاله پژوهشی

نویسندگان

1 گروه علوم دامی، دانشکده کشاورزی، دانشگاه تبریز، تبریز، آذربایجان شرقی، ایران

2 گروه علوم دامی دانشگاه تبریز

چکیده

زمینه مطالعاتی و هدف: به توانایی حیوان برای حفظ بقاء در دوره پرورش، ماندگاری گفته می‌شود. از ویژگی‌های مهم تجزیه و تحلیل صفات ماندگاری، مشاهدات ناتمام (سانسور شده) است. یعنی رخداد مورد مطالعه لزوماً در زمان ثبت مشاهدات اتفاق نمی‌افتد و رخدادها برای برخی افراد اتفاق می‌افتد ولی برای برخی دیگر مشاهده‌ای برای آن رخداد وجود ندارد. آنالیز ماندگاری به روش‌های خطی، آستانه‌ای و توابع خطر نسبی انجام می‌گیرد. مدل‌های خطی کاربرد آسان‌تری دارند ولی به علت غیرخطی بودن صفت بقاء دارای اریب هستند. مدل‌های آستانه‌ای دقت بیشتری داشته و تخمین‌های ژنتیکی بالاتری برآورد می‌کنند. اما مدل‌های خطر نسبی به دلیل استفاده از داده‌های سانسورشده، دقت بالاتری از مدل‌های قبلی دارند. از اینرو، هدف از مطالعه حاضر، برآورد پارامترهای ژنتیکی و تاثیر عوامل غیر ژنتیکی بر روی ماندگاری بره‌های گوسفند نژاد ماکویی تا سن یک و دوسالگی با استفاده از تابع خطر نسبی تحت توزیع ویبول بود. روش کار: 2332 رکورد ثبت شده طی سال‌های 1389-1365 در ایستگاه اصلاح نژاد گوسفند ماکوئی مورد استفاده قرار گرفت. بررسی عوامل غیر ژنتیکی مؤثر بر ماندگاری(جنسیت، تیپ تولد، سال تولد، ماه تولد، وزن تولد، وزن ازشیرگیری، سن ازشیرگیری و سن مادر در هنگام زایش)، با رویه‌ی LIFEREG نرم افزار آماری SAS نسخه 9.4 انجام گردید. برای آنالیز ماندگاری و بررسی پارامترهای ژنتیکی صفت ماندگاری از نرم افزار SURVIVAL KIT استفاده شد. نتایج: تمامی اثرات ثابت به جز سن مادر و سن شیرگیری دارای اثر معنی‌داری بر ماندگاری بره‌ها بودند(05/0>P). نسبت خطر بره‌های نر برای ماندگاری تا سن یک‌سالگی 17/2و دوسالگی 61/2 برابر ماده‌ها بود. نسبت خطر حذف دوقلوها و سه‌قلوها بالاتر از یک‌قلوها بود (22/1-24/1 برای دوقلوها و 49/1-46/1 برای سه‌قلوها). با افزایش وزن تولد و وزن شیرگیری خطر حذف ابتدا کاهش و سپس افزایش می‌یابد. وراثت پذیری برای صفات بقاء تا سن یک و دوسالگی به ترتیب 062/0 و 079/0 برآورد شد. نتیجه گیری نهایی: نظر به اینکه وراثت پذیری در سطح پایینی برآورد شد، کنترل عوامل محیطی و مدیریت، سهم مهمی بر بهبود ماندگاری بره‌های نژاد ماکویی خواهد داشت.

کلیدواژه‌ها


عنوان مقاله [English]

Survival Analysis of Makuie Sheep Breed’s Lambs Using Proportional Hazard Functions

نویسندگان [English]

  • mahdi nezhadali 1
  • Sadegh Alijani 1
  • arash javanmard 2
  • ali hosseinkhani 2
1 department of animal science, faculty of agriculture , university of tabriz, tabriz, east Azerbaijan, iran
2 Assistant Professor, Department of Animal Sciences, University of Tabriz
چکیده [English]

Introduction: The ability of an animal to survive in the breeding period is called the survival trait. One of the most important features of survival analysis is unfinished (censored) observations. In other words, the event under study does not necessarily occur at the time of recording the observations. Thus, survival observations and data are generally incomplete, and events occur for some people but for others there is no observation for that event. In general, nonparametric regression methods (Kaplan-Meier estimator, Cox hazard ratio model) and parametric methods are used for this type of analysis more than other methods. Survival and longevity analysis is performed by linear, threshold and hazard ratio functions methods. Linear models are easier to use, but due to the nonlinearity of the survival trait, they are skewed, and the estimates of heritability with this method are low and close to zero. Threshold models are more accurate and obtain higher genetic estimates. In the analysis of hazard ratio method, the accuracy is higher than previous models due using the more information (use of censored data), but it is time consuming and its computational load is higher. Therefore, the aim of the present study was to estimate the genetic parameters and the effect of fixed factors on the survival of Makuie sheep lambs up to one and two years of age using the proportional hazard function under Weibull distribution.
Material and methods: In this study, 2332 records of Makuie lambs that collected between 1986 to 2010 in Makuie sheep breeding station, were used. In the first step, fixed factors affecting longevity (gender, type of birth, year of birth, month of birth, birth weight, weaning weight, weaning age and age of dam) Was investigated by LIFEREG procedure of SAS statistical software version 9.4 under Weibull distribution. Then significant effects were used in the survival analysis with the proportional hazard model. Survival up to one and two years of age were analyzed using the Weibull hazard model by SURVIVAL KIT software, version 6.12. SURVIVAL KIT software is able to provide components of genetic variance and can estimates breeding values for animals using mixed models. Survival from post-weaning to one year (for analysis up to one year of age) and two years (for analysis up to two years of age) were assessed by day.
Results and discussions: In the analysis of survival from weaning to the one year of age, 38.12% of observations and in the study of survival up to the two years of age, 13.68% of the observations were censored. According to the study, the year of birth had a significant effect on the survival of lambs. The year of birth affects the amount of rainfall and ultimately the availability of food resources by affecting the pasture, and thus changes the quality and amount of nutrition. survival up to one year of age Lambs born in 2004 and for survival up to two years of age lambs born in 1994 had the lowest risk of elimination. Lambs born in 2010 had the highest risk of elimination. Scrutiny the effect of birth month on the survival of lambs, it was found that for survival up to one year of age, lambs born in April had a higher risk of elimination. However, for survival up to two years of age, lambs born in May had the highest risk of elimination and the amount of risk for May compared to March was 1.96. Lambs born in May reach autumn at a young age and during a period of rapid growth. Because the studied flock is kept semi-densely, access to feed for these lambs is reduced and they are removed from the flock due to their low weight and insufficient growth rate. Birth weight is one of the most important components affecting the survival rate and directly affects the elimination rate. Threshold weights have the highest elimination risk ratio and the weight of 3.5-4.7 kg has the lowest elimination risk ratio for survival up to one and two years of age. There is a curved relationship between elimination risk and birth weight, so with increasing birth weight to medium weight (3.5-4.5 kg), the risk ratio decreases and then increases. Mortality at low weights is due to the inability to control body temperature and high mortality at high birth weights is due to infertility. Weaning weight is an important factor influencing survival after weaning. Similar to birth weight, weaning weight has a curved relationship with survival. So that with increasing the weaning weight from 10 kg to 19 kg, the risk of elimination for survival to one and two years of age decreases from 1.8 to 0.97 and 0.93, respectively, and then with increasing the weaning weight, risk ratio reaches to 2.53 and 2.38, respectively. Lack of growth-induced production encourages the elimination of low-weight weaned lambs. So the elimination risk ratio is high at low weights. But as the weaning weight increases, the desire to keep it for later periods will be greater, and thus the removal pressure will be lower. The average weight of weaning in the herd is equal to 19.01 kg, the results also show that the lowest elimination risk ratio and consequently the highest survival for weaning weight is 18-20 kg. Therefore, in order to increase the survival and reduce the elimination risk ratio, management and corrective programs can be designed to select the optimal weaning weight. For survival up to one and two years of age, the risk of elimination for male lambs was higher than females. In the period between weaning until the age of one year, the elimination risk ratio for male lambs was 2.17 relation to females and for survival to the age of two years, the elimination risk ratio of males was 2.61 Compared to females. This difference in elimination rate can be attributed to sex-related factors, sex characteristics that may not have yet been discovered, systemic diseases, different management practices for males and females, and high elimination rate for males due to surplus fattening. There are three types of births in Makuie sheep (singles, twins and triplets). Lambs born with twins and triplets have lower birth weight and weaning weight than single lambs. Due to the above topics, low birth weight and low weaning weight, reduces the survival and, increase the risk of elimination. Another reason for the reduced survival of multiple births is the limitation in dam milk production, which can be due to poor maternal genetics or inadequate nutritional intake. Heritability of survival to one and two years of age for Makuie sheep lambs according to the censorship ratio of 38.12% and 13.68% (respectively) was 0.062 and 0.079, respectively.
Conclusions: Since heritability was estimated at a low level, control of environmental factors and management will play an important role in improving the survival of Makuie lambs.

کلیدواژه‌ها [English]

  • Heritability Estimation
  • survival triat
  • Causes of Elimination
  • Makui Sheep
  • Weibull Model
  • Hazard risk ratio
  • Elimination
Atkins K, 1980. The comparative productivity of five ewe breeds: Lamb growth and survival. Animal Production Science 20: 272-279.
Bagheri M, 2010. Investigation of factors which affect lori-Bakhtiari lambs mortality from birth to weaning. Veterinary Journal (Pajouhesh & Sazandegi) No 85: 8-15 (In Persian).
Bahri Binabaj F, Tahmoorespour M, Aslami Nejad A A, and Vatankhah M, 2013. Determination of the best birth weight for survival to weaning in Karakul breed lambs. Animal Production Research, Vol. 2, No. 1: 53- 63 (In Persian).
Barazandeh A, Moghbeli S M, Vatankhah M and Hossein-Zadeh N G, 2012. Lamb survival analysis from birth to weaning in Iranian Kermani sheep. Tropical Animal Health and Production 44: 929-934.
Chniter M, Hammadi M, Khorchani T, Krit R, Lahsoumi B, Sassi M B and Hamouda M B ,2011. Phenotypic and seasonal factors influence birth weight, growth rate and lamb mortality in D’man sheep maintained under intensive management in Tunisian oases. Small Ruminant Research 99(2): 166-170.
Do Prado Paim T, Da Silva A F, Martins R F S, Borges B O, Lima P D M T, Cardoso C C and McManus C ,2013. Performance, survivability and carcass traits of crossbred lambs from five paternal breeds with local hair breed Santa Inês ewes. Small Ruminant Research 112(1): 28-34.
Ducrocq V and Casella G ,1996. A Bayesian analysis of mixed survival models mixed model, variance component analysis, proportional hazards model. Genetics Selection Evolution. BioMed Central, 1996, 28 (6): 505-529.
Ducrocq V, Sölkner J and Mészáro G, 2010. Survival Kit V6—A software package for survival analysis. 9th World Congr. Genet. Appl. Livest. Prod., Leipzig, Germany. Communication, (0232).
Getachew T, Gizaw S, Wurzinger M, Haile A, Rischkowsky B, Okeyo A, Sölkner J and Mészáros G ,2015. Survival analysis of genetic and non-genetic factors influencing ewe longevity and lamb survival of Ethiopian sheep breeds. Livestock Science 176: 22-32.
Hatcher S, Atkins K and Safari E ,2010. Lamb survival in Australian Merino sheep: a genetic analysis. Journal of Animal Science 88: 3198-3205.
Jafari S, Hashemi A, Manafiazar G, Darvishzadeh R, Razzagzadeh S and Farhadian M ,2012. Genetic analysis of growth traits in Iranian Makuiee sheep breed. Italian Journal of Animal Science 11(1): 98-102
Kalbfleisch J D and Prentice R L ,2011. The statistical analysis of failure time data (Vol. 360). John Wiley & Sons.INC., Publication.
Kern G, Kemper N, Traulsen I, Henze C, Stamer E and Krieter J, 2010. Analysis of different effects on longevity in four sheep breeds of northern Germany. Small Ruminant Research 90(1): 71-74.
Mandal A, Prasad H, Kumar A, Roy R and Sharma N ,2007. Factors associated with lamb mortalities in Muzaffarnagari sheep. Small Ruminant Research 71: 273-279.
Maxa J, Sharifi A, Pedersen J, Gauly M, Simianer H and Norberg E ,2009. Genetic parameters and factors influencing survival to twenty-four hours after birth in Danish meat sheep breeds. Journal of Animal Science 87: 1888-1895.
Mohammadi H, Shahrebabak MM, Sharebabak HM and Vatankhah M ,2012. Estimation of genetic parameters of reproductive traits in Zandi sheep using linear and threshold models. Czech Journal of Animal Science 57: 382-388.
Nezhadali M, Alijani S and Javanmard A, 2020. Comparison of Non-Linear Models to Describe of Growth Pattern in Makuiee Sheep. Research On Animal Production 11(27): 88-94.
Riggio V, Finocchiaro R and Bishop S ,2008. Genetic parameters for early lamb survival and growth in Scottish Blackface sheep. Journal of Animal Science 86: 1758-1764.
Sawalha R M, Conington J, Brotherstone S and Villanueva B ,2007. Analyses of lamb survival of Scottish Blackface sheep. Animal-Cambridge University Press-, 1(1): 151.
Seasakhti D, Vatankhah M, Merzaei H R and Yousef Ellahi M ,2010. Estimates of some environmental factors and genetic parameters on Lori-Bakhtiari lambs’ survival. Animal Science Journal 84: 65-70. (In Persian)
Southey B R, Rodriguez-Zas S.L and Leymaster K A ,2001. Survival analysis of lamb mortality in a terminal sire composite population. Journal of Animal Science 79(9): 2298-2306.
 Valizadeh R ,2010. Iranian sheep and goat industry at a glance. Ferdowsi University of Mashhad .1–9.
Vatankhah M ,2013. Estimation of the genetic parameters for survival rate in Lori-Bakhtiari lambs using linear and Weibull proportional hazard models. Journal of Agricultural Science and Technology 15: 1133-1143.
Vostrý L and Milerski M ,2013. Genetic and non-genetic effects influencing lamb survivability in the Czech Republic. Small Ruminant Research 113(1): 47-54.
Yazdi M H, Visscher P M, Ducrocq V and Thompson R ,2002. Heritability, reliability of genetic evaluations and response to selection in proportional hazard models. Journal of Dairy Science 85(6): 1563-1577.