مقایسه روش های Gblup و بیزی در برآورد ارزش های اصلاحی ژنومی با معماری های مختلف ژنتیکی

انتخاب ژنومی با استفاده از نشانگرهای با تراکم بالا، بخصوص snpهایی که کل ژنوم را پوشش می دهند و اغلب در عدم تعادل پیوستگی با qtlهای مجاور خود قرار دارند، ارزش ژنتیکی کل را پیش بینی می کند. این پژوهش با هدف بررسی تاثیر عواملی چون تراکم نشانگرها، تعداد qtl، وراثت پذیری صفت و نوع توزیع اثر qtl بر صحت برآورد ارزش های اصلاحی ژنومی با استفاده از داده های شبیه سازی در گوسفند انجام گرفت. به همین منظور ژنومی متشکل از سه کروموزوم، هر یک به طول 100 سانتی مورگان با سه مقدار وراثت پذیری برای صفت مورد بررسی (0.1، 0.3 و 0.5) و سه پنل نشانگری (500، 1000 و 1500) در سه سطح تعداد qtl (50، 100 و 150) با دو اثر توزیع یکنواخت و گاما برای qtl شبیه سازی شدند. صحت ارزش های اصلاحی ژنومی برآورد شده با استفاده از پنج روش gblup، بیزa، بیزb، بیزc و بیز lasso مورد مقایسه قرار گرفتند. نتایج این تحقیق نشان داد که هر چه تراکم نشانگر و وراثت پذیری صفت افزایش یافته و تعداد qtl موثر بر صفت کمتر باشد، صحت ارزش اصلاحی برآورد شده بالاتر خواهد بود. در بین روش های آماری زمانی که تعداد qtl موثر بر صفت پایین است و توزیع اثر qtl گاما در نظر گرفته شد، پیش بینی ارزش های اصلاحی ژنومی با روش بیزb عملکرد بهتری داشت.

A Comparison of GBLUP and Bayesian Methods in Prediction of Genomic Breeding Values under Different Genetic Architectures

Introduction[1] Genomic Selection (GS) has been proved to be a powerful tool for estimating genetic values in livestock breeding. Newly developed sequencing technologies have dramatically reduced the cost of genotyping and significantly increased the scale of genotype data that used for GS. The estimation of breeding values in order to select the best animals as parents of the next generation is the main goal of animal breeding programs. Traditional methods of genetic evaluation were performed using a combination of phenotypic and pedigree information to produce estimated breeding values. Most simulation studies of genomic selection (GS) methods have considered genetic architectures in which the number and relative magnitudes of quantitative trait loci (QTL) have varied. Among the Bayesian methods, those using markerspecific shrinkage of effects (e.g., BayesA or BayesB of or the Bayesian LASSO are commonly used in animal breeding applications. The Bayesian methods proposed differ in the way of looking at the variances of parameters. In classical livestock breeding methods, selection for important economic traits using pedigree information with individual phenotypic records was performed and best Linear Prediction of Breeding Values (BLUP) is achieved. In genome selection, genomic breeding values of all individuals can be predicted with high accuracy using a linear model. Various factors can be affecting the accuracy of genomic breeding values. Therefore, the present study aimed to evaluate the accuracy of estimating genomic breeding values in different genetic architectures including different distributions of gene effects, different numbers of QTL, different levels of heritability and different marker densities using GBLUP and Bayesian methods including Bayes A, Bayes B, Bayes C and Bayes LASSO. In addition to comparing the performance of different methods in different genetic architectures, a marker density and QTL numbers were introduced for simulation programs of sheep populations. Materials and Methods To create a basic population (G0), 100 heads of livestock, including 50 males and 50 females, were considered. The frequency of primary alleles for singlenucleotide polymorphisms in the basal generation was considered to be 0.5. To create the first generation (G1), the parents were randomly selected from the males and females of the G0 generation. Parental gametes were simulated based on the assumption of disconnection imbalance using the Halldan location function method, and then randomly generated gametes were randomly selected and mixed to create a new generation of G1 generation. A genome with a length of 300 cM was simulated and 500, 1000 and 1500 SNPs were equally spaced over the chromosome. Three different numbers of QTL (50, 100 and 150) were considered and QTLs were uniformly distributed over the chromosome. One hundred individuals, including 50 males and 50 females, were simulated for the base population. The first generation structure was followed through to the 50th generation of random mating to make linkage disequilibrium populations. Generation 51 was assumed as a training population and the other generations (52 to 60) as validation populations. Five methods, GBLUP, Bayes A, Bayes B, Bayes C and Bayesian LASSO, were used to estimate genomic breeding values. Results and Discussion In all five methods, the accuracy of genomic values decreased as the number of QTLs increased from 50 to 150. The reason for this can be attributed to the limited amount of genetic variance distributed over many QTLs. Also predicting accuracy of all five methods increased with increasing marker density. Results showed that increasing marker density at low (0.1) and high (0.5) heritability levels, increased genomic accuracy but increasing at moderate heritability (0.3) traits did not affect the accuracy of genomic evaluation. Accuracy of genomic breeding values in the gamma distribution provides better gene effects to uniform distributions. Conclusion The results showed that factors such as marker density, QTL numbers, distribution QTL effect and trait heritability were effective in estimating the accuracy of genomic breeding values. In high heritability traits, the higher markers density and lower QTL numbers, leading to increase accuracy of estimating genomic breeding values. In genomic studies, if the trait is affected by a small number of QTLs, estimation of breeding values by Bayes B method can yield a more favorable result. Marker densities did not affect the accuracy of genomic evaluation in traits of moderate heritability, and since most of the economic traits in native species of sheep are moderate heritability, 500 to 1000 markers can be used to estimate breeding values in simulation programs.