مطالعه توزيع آماري اثرات QTL بر صحت ارزش¬هاي اصلاحي ژنومي برآورد شده با روش Bayesian

Study of QTL Effects Distribution on Accuracy of Genomic Breeding values Estimated Using Bayesian Method

هدف تحقیق بررسی اثر دو توزیع گاما و بتا به طور جداگانه بر صحت ارزیابی ژنتیكی بود. به این منظور ژنومی متشكل از 10 كروموزوم هر یك به طول 200 سانتی¬مورگان شبیه¬سازی شد. نشانگرها بر روی ژنوم با فواصل 2/0 سانتی¬مورگان طراحی شدند و جایگاه¬های ژنی مؤثر بر صفات بر روی ژنوم با توزیع تصادفی و تعداد متغیر شبیه¬سازی شدند و تنها اثرات افزایشی ژن¬ها در نظر گرفته شد. در ابتدا جمعیت پایه¬ای از حیوانات با اندازه مؤثر 100 شبیه¬سازی شد و این ساختار جمعیتی برای 50 نسل با آمیزش تصادفی ادامه یافت تا عدم تعادل پیوستگی بین نشانگرها و QTL ایجاد شود. پس از این اندازه جمعیت به 1000 فرد در نسل 51 (نسل مرجع) گسترش یافت. در نسل مرجع براساس اطلاعات ژنومیك و فنوتیپی اثرات نشانگری محاسبه گردید. در نسل 52 تا 57 (نسل هدف) ارزش اصلاحی محاسبه شد. نتایج نشان داد كه در همه توزیع¬ها هر چه از جمعیت مرجع دور می¬شویم، با افزایش نسل¬ها از نسل 51 به 57 صحت برآورد ارزش اصلاحی ژنومی كاهش می‏یابد. همچنین وراثت¬پذیری بالا (2/0) نسبت به وراثت¬پذیری پایین (05/0) در تعداد QTL مشابه از صحت بیشتری برخوردار می¬باشد. در مقایسه توزیع¬ها در وراثت¬پذیری پایین، با تعداد 10 QTL، توزیع گاما 2، با تعداد 20 QTL، توزیع گاما 1 و با تعداد 50 و 100 QTL توزیع بتا برآوردها صحت بیشتری در دو روش Lasso و Ridge داشتند؛ در وراثت¬پذیری بالا، با تعداد 50 و 100 QTL توزیع گاما 2 در دو روش Lasso و Ridge برتری داشت. به طوركلی توزیع گاما باعث افزایش صحت برآورد ارزش اصلاحی ژنومی شد.

Introduction Genetic evaluation and estimation of breeding value are one of the most fundamental elements of breeding programmes for genetic improvement. Recently, genomic selection has become an efficient method to approach this aim. The accuracy of estimated Genomic breeding value is the most important factor in genomic selection. Different studies have been performed addressing the factors affecting the accuracy of estimated Genomic breeding value. The aim of this study was to evaluate the effect of beta and gamma distributions on the accuracy of genetic evaluation. Materials and Methods A genome consisted of 10 chromosomes with 200 cm length was simulated. Markers were spaced on 0.2 cm intervals and different numbers of QTL with random distribution were simulated. Only additive gene effects were considered. The base population was simulated with an effective size of 100 animals and this structure continued up to generation 50 to creating linkage disequilibrium between the markers and QTL. The population size was increased to 1000 animals in generation 51 (reference generation). Marker effects were calculated from the genomic and phenotypic information. Genomic breeding value was computed in generations 52 to 57 (training generation). Effects of gamma 1 distribution (shape=0.4, scale=1.66), gamma 2 distribution (shape=0.4, scale=1) and beta distribution (shape1=3.11, shape2=1.16) were studied in the reference and training groups. The heritability values were 0.2 and 0.05. Results and Discussion The results showed that accuracy of genomic breeding value reduced with passing generation (from 51 to 57) for two gamma distributions and beta distribution; this decrease may be due to two factors: recombination has negative impact on accuracy of genomic breeding value and selection reduces genetic variance as the number of generations increases. Accuracy of genomic estimated breeding value increased as the heritability increased so that the high heritability had more accuracy than low heritability in same QTL number. Number and distribution of genes is an important factor for accuracy of estimated breeding value. Duncan test was conducted by SPSS software. Results illustrated that there was no significant difference between the different distributions. Comparing accuracy of estimated breeding value showed that in the low heritability scenario with 10 and 20 QTL, gamma distribution 2 and gamma distribution 1 performed well, respectively, whilst in 50 and 100 QTL scenario, beta distribution was superior in both Lasso and Ridge methods. In the high heritability scenario with 50, 100 QTL gamma distributions 2 were superior in both Lasso and Ridge methods. With four QTL (10, 20, 50 & 100), in high heritability scenario, estimated genomic breeding value was often increased by increasing the number of QTL. This may be due to increasing linkage disequilibrium between markers and QTL. In general, the gamma distribution led to the increased accuracy of the estimations in both Lasso and Ridge methods. Conclusion Marker density, method to estimate marker effects, QTL distribution, number of QTL, number of generations and trait heritability are some effective factors on accuracy of estimated genomic breeding value. The accuracy of estimated genomic breeding value is output of these factors and the distribution of genes is an important factor for accuracy of estimated genomic breeding value. We can conclude that, accuracy is reduced with increasing number of generations from base population to training population while the accuracy of estimated genomic breeding value is increased when breeding value of the reference group is used in lieu of the phenotypic records. In addition, accuracy of estimated genomic breeding value is enhanced by increasing heritability, so that, between three the distributions simulated in high heritability scenario, gamma 2 distribution increased accuracy of the estimates. Although, the size and distribution of QTL effects still greatly influence the effectiveness of the genomic prediction methods, but as a suggestion, models of genetic variation of genomic assessment should be considered, since a method of estimating breeding value may have (or produce) a better estimation with a specific model.