Analysis of genetic diversity in five Iranian sheep population using microsatellites markers

Document Type : Research Paper

Authors

Abstract

Indigenous races in each country are as a national capital and strategic product in the economy and prosperity of the country that maintenance of these races are very valuable. Due to the importance of maintaining diversity and identification of genetic resources in indigenous breeds, five breeds of Iranian sheep population (Arabi, Arman, Dalagh, Karakul and Lory) using four microsatellite markers (McMA2, McMA26, OarHH35 and BM6444 ) was evaluated in terms of genetic diversity. Genomic DNA was extracted from 225 blood samples by optimized salting-out DNA extraction procedure. The PCR reactions were successfully performed with all primers. The results showed that all loci were quite polymorph. Hardy-Weinberg equilibrium test using chi-square test showed that some various combinations of loci-population were in a state of Hardy-Weinberg disequilibrium (P<0.05). The highest number of alleles was observed in loci of McMA2, McMA26 and OarHH35 of Dalagh, Arabic, Lory and Karakul populations respectively (12 alleles) and the lowest number of alleles for OarHH35 in Dalagh population. The maximum expected heterozygosity in combination locus-population of loci McMA2، McMA26، OarHH35، BM6444 were 0.885, 0.9, 0.88 and 0.87 respectively. With respect to the studied loci, it is concluded that all studied sheep populations have wide genetic diversity.

Keywords


بررسی تنوع ژنتیکی پنج جمعیت گوسفند ایرانی با استفاده از نشانگرهای ریزماهواره­ای

 

محمدتقی واجدابراهیمی*1، محمدرضا محمدآبادی2، علی اسماعیلی زاده 3

 

1 دانشجوی کارشناسی ارشد ژنتیک و اصلاح دام گروه علوم دامی، دانشگاه شهید باهنر کرمان

2 دانشیار گروه علوم دامی، دانشگاه شهید باهنر کرمان

3 استاد گروه علوم دامی، دانشگاه شهید باهنر کرمان

تاریخ دریافت: 05/07/1393، تاریخ پذیرش: 01/09/1394

 

چکیده

نژادهای بومی در هر کشور به‌عنوان یک سرمایه ملی و محصولی استراتژیک در اقتصاد و رونق آن کشور محسوب می‌شوند که حفظ و نگهداری این نژادها بسیار ارزشمند است. با توجه به اهمیت حفظ تنوع در نژادهای بومی و شناسایی ذخایر ژنتیکی، پنج جمعیت از نژادهای گوسفند ایرانی (عربی، آرمان، دالاق، قره­گل و لری) با استفاده از 4 نشانگر ریزماهواره­ای (McMA2، McMA26، OarHH35 و BM6444) از نظر تنوع ژنتیکی مورد بررسی قرار گرفت. از تعداد 225 نمونه خون DNA ژنومی به روش استخراج نمکی بهینه‌شده، به دست آمد. واکنش‌های PCR به ‌خوبی انجام شد. نتایج نشان داد که تمامی جایگاه­ها کاملاً چندشکل بودند. آزمون تعادل هاردی–وینبرگ به روش آزمون مربع کای نشان داد که برخی ترکیبات مختلف جایگاه­ها-جمعیت در حالت عدم تعادل قرار داشتند (05/0P<). بیشترین تعداد آلل مشاهده ‌شده در جایگاه­های McMA2، McMA26 و OarHH35 به ترتیب مربوط به جمعیت­های دالاق، عربی، لری و قره­گل (12 آلل) و کمترین تعداد از آن جایگاه OarHH35 در جمعیت دالاق بود. حداکثر هتروزیگوسیتی مورد انتظار در حالت ترکیب جمعیت-جایگاه برای جایگاه­های McMA2، McMA26، OarHH35، BM6444 به ترتیب 885/0، 9/0، 88/0و 87/0 بودند. در مجموع می­توان نتیجه گرفت که جمعیت­های گوسفند مورد بررسی در این پژوهش با توجه به جایگاه­های مورد مطالعه، از تنوع ژنتیکی نسبتاً بالایی برخوردارند.

واژه­های کلیدی: گوسفند، نشانگرهای ریزماهواره­ای، تنوع ژنتیکی، چندشکلی.



مقدمه

 

نژادهای بومی در هر کشور به‌عنوان یک سرمایه ملی و محصولی استراتژیک در اقتصاد و رونق آن کشور محسوب می­شوند و نیز به دلیل شدت انتخاب پایین، تعداد زیاد پرورش‌دهندگان و محدودیت استفاده از تلقیح مصنوعی، معمولا از تنوع ژنتیکی بالایی برخوردارند. از طرفی باگذشت زمان و کسب آگاهی بیشتر نسبت به اهمیت صفات مختلف، نیازهای جدیدی مطرح می­شود که متخصصین اصلاح نژاد را بر آن می­دارد که از ذخیره ژنی دام­های بومی استفاده نمایند. این مسئله، بخصوص با افزایش تولید محصولات دامی و تولید محصولات پیش‌بینی‌نشده در آینده، لزوم حفظ تنوع ژنتیکی در دام­های بومی را الزامی ساخته است (Frankham, 1994) چرا که یک گونه بدون تنوع ژنتیکی کافی قادر به سازگاری با تغییرات محیطی و مبارزه با انگل­ها نیست (Askari et al., 2011). در ایران بیش از 50 میلیون رأس گوسفند، شامل 27 نژاد و اکوتیپ وجود دارد (Zamani et al., 2013) که باید تنوع ژنتیکی آنها محاسبه شده و سپس در حفظ این نژادهای بومی تلاش شود. حفاظت باید بر اساس دانش عمیقی از منابع ژنتیکی نژادهای خاص باشد، لذا تلاش برای شناسایی و تعیین خصوصیات ژنتیکی نژادهای بومی و محلی بسیار اهمیت دارد (Shojaei et al., 2010). استفاده از نشانگرهای مولکولی در سال­های اخیر جهت تعیین تنوع ژنتیکی بین جمعیت­ها و حیوانات حفاظت‌شده، کاربرد گسترده­ای یافته است. میزان چندشکلی به‌دست‌آمده از این نشانگرهای ژنتیکی، یکی از پارامترهای قابل ارزیابی برای مطالعه جمعیت­های مختلف و درک تفاوت­های ژنتیکی بین جمعیت­هاست (Davis et al., 2002). در این میان، ریزماهواره‌ها به سبب برتری‌هایی که دارند، برای مطالعات ژنتیک جمعیت نسبت به سایر نشانگرها ارجح هستند. انجمن بین المللی ژنتیک حیوانی[1] ریزماهواره­ها را به عنوان بهترین نشانگر جهت تعیین تنوع ژنتیکی گونه­های حیوانی معرفی کرده است (Mohammadifar and Mohammadabadi, 2011). بر اساس بررسی­های ثبت شده توسط FAO، 66 درصد کل مطالعات تعیین فاصله ژنتیکی با استفاده از ریزماهواره­ها انجام گرفته است و 70 درصد پژوهشگران ریزماهواره­ها را انتخاب کرده­اند (Mohammadifar and Mohammadabadi, 2011).

در مطالعه­ای Esmail-Khanian et al (2007) تنوع ژنتیکی 156 رأس گوسفند بلوچی ایران را با استفاده از 19 جایگاه ریزماهواره­ای بررسی کردند. آن­ها گزارش کردند که بر اساس آزمون تعادل هاردی-وینبرگ این جمعیت در بیشتر جایگاه­ها انحراف معنی‌داری از حالت تعادل نشان داد (05/0P<). دامنه هتروزیگوسیتی برای این جایگاه­ها بین 1/0 تا 93/0 متفاوت بود. درمجموع آن­ها نشان دادند که بیشتر جایگاه­های موردمطالعه دارای چندشکلی بوده و می­توان از آن‌ها برای مطالعات بعدی استفاده کرد. در پژوهشی دیگر Banabazi et al (2007) تنوع ژنتیکی درون و بین پنج جمعیت گوسفند ایرانی (سنجابی، کردی کردستان، کردی خراسان، مهربان و مغانی) را با استفاده از شش جفت آغازگر ریزماهواره­ای بررسی کردند و نشان دادند که تمامی جایگاه­ها چندشکل بودند، تمامی جمعیت­های موردمطالعه در جایگاه­های ریزماهواره­ای موردبررسی، در تعادل هاردی-وینبرگ بوده و بیشترین و کمترین هتروزیگوسیتی مورد انتظار در همه جایگاه­ها برآورد گردید که به ترتیب مربوط به جمعیت مهربان (847/0) و کردی خراسان (744/0) بود.

میزان چندشکلی برخی از نشانگرهای ریزماهواره­ای­ پیوسته به ژن دوقلوزایی FecB را Mohammadi and Saberivand (2007) در گوسفندان قزل بررسی کردند. آنها میزان چندشکلی 3 جفت پرایمر ریزماهواره­ای متصل به ژن چندقلوزایی FecB در 300 گوسفند نژاد قزل را مطالعه کردند. جایگاه OarHH55  و OarHH35 به ترتیب با 6 و 9 آلل کمترین و بیشترین آلل را تولید کردند. میزان هتروزیگوسیتی و میزان شاخص شانون برای جفت پرایمرهای OarHH55، OarHH35 وBM1329  به ترتیب (9333/0 و 7461/1)، (96/0 و 0786/2) و (91/0 و 5698/1) بود و نتایج به‌دست‌آمده نشان داد که هر سه جایگاه مورد مطالعه از حالت تعادل هاردی-وینبرگ انحراف داشتند.

برای مطالعه 120 گوسفند زندی Nanekarani et al (2010) از 15 نشانگر ریزماهواره­ای استفاده کردند. آن­ها بیان داشتند که همه جایگاه­های مورد مطالعه چندشکل بودند و به خاطر فزونی هتروزیگوت­ها در همه جایگاه­ها انحراف از تعادل مشاهده شد (01/0P<). میانگین محتوای چندشکلی در جمعیت به ازای تمام جایگاه­ها، برابر 808/0 بود که نشان دهنده تغییرپذیری بسیار بالا در جمعیت موردمطالعه بود. نتایج پارامترهای تنوع در این مطالعه نشان داد که تنوع ژنتیکی بالایی در نژاد زندی وجود دارد. در مطالعه­ای به ‌منظور بررسی مقدار کارایی و سودمندی نشانگرهای ریزماهواره­ای ویژه کروموزوم Y گاوی در تعیین تنوع ژنتیکی در گوسفند، وضعیت تکثیر 17 نشانگر ریزماهواره­ای ویژه کروموزوم Y گاوی در گوسفند کرمانی مورد بررسی قرار گرفت (Mohammadifar and Mohammadabadi, 2011). در مجموع 102 آلل در این جایگاه‌ها شناسایی شدند که بیشترین و کمترین باندهای مشاهده‌شده به ترتیب 16 و 5 باند بود. نتایج این مطالعه نشان داد که می­توان از ریزماهواره­های ویژه کروموزوم Y گاوی به‌منظور بررسی تنوع ژنتیکی، فراوانی ژنی، هتروزیگوسیتی و فاصله ژنتیکی در گوسفندان استفاده کرد. در پژوهشیTomasco et al. (2002) چندشکلی 10 جایگاه ریزماهواره­ای در گله­های گوسفند اروگوئه را مورد ارزیابی قرار دادند. در این بررسی جایگاه­ها چندشکلی بالایی نشان دادند و تعداد آلل گزارش‌شده برای این 10 جایگاه 7 تا 15 آلل و شاخص PIC بین 63/0 و 87/0 متغیر بود. نتایج نشان داد که این سری جایگاه­ها برای آزمون ارتباطات ژنتیکی مؤثر می­باشند. در مطالعه دیگری میزان تنوع ژنتیکی نژاد گوسفند نجدی[2] عربستان با 19 نشانگر ریزماهواره­ای اندازه­گیری شد و درمجموع 173 آلل مشاهده گردید. به‌طور متوسط هتروزیگوسیتی قابل‌مشاهده، محتوای اطلاعات چندشکلی (PIC) و شاخص شانون به ترتیب 67/0، 71/0 و 69/1 به دست آمد و نتایج نشان داد که این نژاد از تنوع ژنتیکی مناسبی برخوردار است و ازلحاظ هم­خونی درخطر نیست و می­توان از تنوع ژنتیکی درون این نژاد برای پیشرفت روند ژنتیکی استفاده کرد (Musthafa Muneeb et al., 2012). در مطالعه­ای دیگر Jakaria et al. (2012) با استفاده از 18 جایگاه ریزماهواره­ای تنوع ژنتیکی جمعیت گوسفندهای اندونزی را بررسی کردند. در این پژوهش 180 آلل از 17 جایگاه شناسایی شدند. متوسط هتروزیگوسیتی مشاهده‌شده و مورد انتظار به ترتیب 5749/0 و 6896/0 بود. میانگین تمایز ژنتیکی برای هم‌خونی در بین و داخل جمعیت‌ها و همچنین میانگین تمایز ژنتیکی جمعیت گوسفندان به ترتیب 1006/0، 1647/0 و 0712/0 بود. بر اساس نتایج این پژوهش، جمعیت گوسفندهای اندونزی به یک برنامه حفاظت ژنتیکی و اصلاحی نیاز دارند. با توجه به این ‌که نژادهای عربی (47 رأس)، آرمان (46 رأس)، دالاق (44 رأس)، قره­گل (43 رأس) و لری (45 رأس) ایران تاکنون با نشانگرهای ریزماهواره­ای مورد مطالعه قرار نگرفته­اند، لذا هدف تحقیق حاضر مطالعه تنوع درون و بین جمعیتی موجود در این پنج نژاد گوسفند ایرانی به کمک 4 جایگاه ریزماهواره­ای بود.

 

مواد و روش‌ها

در مطالعه حاضر از 225 گوسفند شامل نژادهای عربی (47 رأس)، آرمان (46 رأس)، دالاق (44 رأس)، قره­گل (43 رأس) و لری (45 رأس) از مراکز اصلاح نژاد و جهاد کشاورزی نقاط مختلف کشور، نمونه‌های خون کامل از سیاهرگ وداج و با استفاده از لوله CBC/5ml  حاوی محلول ضد انعقاد K2/K3 EDTAجمع‌آوری و تهیه گردید. استخراج DNA از نمونه‌های خون کامل به روش بهینه‌شده و تغییریافته استخراج نمکی (Abadi et al., 2009) انجام گردید که برای تجزیه سلول­های قرمز از بافر شستشو و نیز برای تجزیه اسیدهای نوکلئیک از بافر لایسیز، استفاده گردید. مشخصات جایگاه­های موردمطالعه در جدول 1 آورده شده است. واکنش­های PCR با استفاده از 4 جفت آغازگر ریزماهواره­ای McMA2، McMA26، OarHH35 و BM6444 (در حجم نهایی 25 میکرولیتر) انجام شد (جدول 2). در این پژوهش از دستگاه ترموسایکلر peqSTAR 96X ساخت آلمان، با بلوک 96 تایی برای تیوب‌های PCR به حجم 2/0 میلی‌لیتر، استفاده شد. برنامه دمایی برای آغازگرها موردمطالعه به‌صورت واسرشته سازی اولیه در دمای 94 درجه سانتی­گراد به مدت 5 دقیقه و به دنبال آن 30 چرخه شامل واسرشته سازی در دمای 95 درجه سانتی­گراد به مدت 30 ثانیه، اتصال آغازگرها با توجه به دمای اتصال بهینه‌شده برای هر جایگاه و به مدت 45 ثانیه، بسط اولیه در دمای 72 درجه سانتی­گراد به مدت 30 ثانیه و سرانجام بسط نهایی در دمای 72 درجه سانتی­گراد به مدت 5 دقیقه بهینه شد. فرآورده­های PCR به‌دست‌آمده روی ژل پلی‌آکریل‌آمید 8 درصد الکتروفورز گردید و نمایان‌سازی آلل­ها به روش رنگ­آمیزی نیترات نقره انجام شد (Promega, 1993). آزمون نسبت درست نمایی برای تعادل هاردی- وینبرگ (05/0P<) در هریک از جایگاه­ها، همچنین تنوع درون جمعیت به‌صورت متوسط هتروزیگوسیتی مورد انتظار و معیارهای چندشکلی ازجمله تعداد آلل مشاهده‌ شده و مؤثر با استفاده از نرم‌افزار PopGene (Yeh et al., 1999) انجام شد. محتوای اطلاعات چندشکلی هریک از جایگاه­ها نیز با استفاده از نرم‌افزار Excel و فرمول PIC= N(1-∑Pi^2 )/N-1 به دست آمدند(Buchanan and Thue., 1998).

 

نتایج و بحث

نتایج حاصل از الکتروفورز ژل آگارز 8/0 درصد نشان داد که کیفیت DNA های استخراج‌شده با استفاده از روش نمکی بهینه‌شده بسیار مناسب و تقریباً تمامی باندهای نمونه‌ها شفاف و بدون کشیدگی بودند (شکل 1).

واکنش‌های PCR با چهار جفت آغازگر انجام گردید، تمامی جایگاه‌های موردمطالعه طبق برنامه­های دمایی به‌دست‌آمده برای آن جایگاه، تکثیر شدند و به‌منظور تائید تکثیر قطعات موردنظر، محصولات PCR با استفاده از الکتروفورز روی ژل آگارز 1% مورد ارزیابی قرار گرفتند (شکل2). پس از تائید تکثیر قطعات موردنظر، الکتروفورز روی ژل آکریل‌آمید محصولات PCR انجام گرفت و سپس ژل‌های آکریل‌آمید رنگ‌آمیزی شدند. شکل 3 نمونه‌­ای از الگوهای باندی به‌دست‌آمده را برای جایگاه‌ BM6444 در جمعیت عربی نشان می‌دهد.

 

 

جدول 1- مشخصات جایگاه‌های انتخاب‌شده برای تحقیق حاضر.

Table 1- Characteristics of selected microsatellites in present study.

منبع

Reference

شماره دسترسی[3]**

Accession Numbers

آغازگرها

(5´-3´)

Primers

نام جایگاه و

موقعیت کروموزومی*

Marker and Chromosome

Maddox et al. (2000)

AF098773

TCA CCC AAC AAT CAT GAA AC

TTA AAT CGA GTG TGA ATG GG

McMA2

(13)

Maddox et al. (2000)

AF098961

TCT CTG CTT TCC AGC CTT ATT C

AGA GCT TTT AGG ACA GCC ACC

McMA26

(18)

Henry et al.

(1993)

 

L12554

AAT TGC ATT CAG TAT CTT TAA ACA TCT GGC/  ATG AAA ATA TAA AGA GAA TGA ACC ACA CGG

OarHH35

(4)

Buchanan et al.

(1994)

G18444

CTC TGG GTA CAA CAC TGA GTC C

TAG AGA GTT TCC CTG TCC ATC C

BM6444

(2)

* اعداد داخل پرانتز شماره کروموزومی است که ریزماهواره مربوطه بر روی آن واقع‌شده است.

** این شماره‌ها کُد دسترسی به هر یک از این ریزماهواره‌ها در بانک ژن NCBI با آدرس اینترنتی http://www.ncbi.nlm.nih.gov:80/entrez/viewer.cgi است.

* The numbers in parenthesis is chromosomes of marker that is located on it.

** The accession number of the markers in the NCBI and the URL is http://www.ncbi.nlm.nih.gov:80/entrez/viewer.cgi.

 

         

 

 

جدول 2- غلظت اجزاء مختلف PCR.

Table 2- The concentration of PCR components.

اجزاء واکنش

Components

حجم مورداستفاده در واکنش (میکرولیتر)

volume reaction (ml)

غلظت استوک

Stoke concentration

DNA

2

50-100 ng / µl

H2o

متغیر( changeable)

بافر (Buffer)

2.5

10X

dNTPs

0.5

10 mM

MgCl2

1

1.5 Mm

Primer Forward

1

pM 10

Primer Reverse

1

pM 10

TaqDNA Polymerase

0.3

5 Unit/ µl

حجم نهائی واکنش(Total)

25


 

 

شکل 1- الکتروفورز چند نمونه DNA استخراج‌شده گوسفند روی ژل آگارز 8/0 درصد.

Figure 1- Electerophoresis of some sheep DNA samples on 0.8% Agarose gel.

 

شکل 2- نمونه­ای از الکتروفورز محصولات PCR جایگاه OarHH35 در جمعیت آرمان روی ژل آگارز 1%.   L: DNA Ladder 50 bp، 1-7: محصولات PCR

Figure 2- A sample of electrophoresis of OarHH35 locus PCR products in Arman population on the 1% Agarose gel. L: DNA Ladder 50 bp, Lanes 1-7: PCR products

 

شکل 3- نمونه­ای از الکتروفورز محصولات PCR جایگاه BM6444 در جمعیت عربی روی ژل پلی‌آکریل‌آمید 8%.   L: DNA Ladder 50 bp، 1-14: محصولات PCR.

Figure 3- A sample of electrophoresis of BM6444 locus PCR products in Arabi population on the 8% Polyacrylamide gel. L: DNA Ladder 50 bp, Lanes 1-14: PCR products.

جدول 3- تعداد آلل واقعی (Na)، تعداد آلل مؤثر (Ne) و هتروزیگوسیتی مورد انتظار (He) برای ترکیبات مختلف جایگاه-جمعیت و برای هر جمعیت.

Table 3- The number of actual alleles (Na), effective alleles (Ne) and expected heterozygosity (He) for different combinations of locus-population and for each population.

جایگاه (Marker)

 

جمعیت (Population)

McMA2

BM6444

McMA26

OarHH35

Na

Ne

He

Na

Ne

He

Na

Ne

He

Na

Ne

He

عربی (Arabi)

11

8.2

0.88

11

7.7

0.87

12

10.3

0.9

9

7

0.85

آرمان (Arman)

10

7.88

0.875

9

7.1

0.86

11

7.9

0.87

9

7.1

0.86

دالاق (Dalagh)

12

8.56

0.885

10

8.01

0.87

11

7.7

0.87

8

5.9

0.83

قره­گل (Karakul)

11

7.73

0.87

10

7.29

0.86

11

8.04

0.87

12

8.6

0.88

لری (Lory)

10

7.27

0.86

9

7.14

0.86

12

8.42

0.88

9

7.19

0.85

 

 

 

جدول 4- معیار شاخص Shannon () و محتوای اطلاعات چندشکلی (PIC) برای ترکیبات مختلف جایگاه-جمعیت و برای هر جمعیت.

Table 4- Shannon index () and polymorphism information conten (PIC) for different combinations of locus-population and for each population.

   جایگاه(Marker)

 

جمعیت (Population)

McMA2

BM6444

McMA26

OarHH35

I*

PIC

I

PIC

I

PIC

I

PIC

عربی (Arabi)

2.266

0.88

2.23

0.88

2.396

0.9

2.065

0.87

آرمان (Arman)

2.18

0.87

2.074

0.86

2.242

0.87

2.07

0.86

دالاق (Dalagh)

2.311

0.88

2.184

0.87

2.198

0.87

1.926

0.83

قره­گل  (Karakul)

2.22

0.87

2.142

0.86

2.225

0.87

2.316

0.88

لری (Lory)

2.137

0.86

2.075

0.86

2.30

0.88

2.05

0.85

* شاخص اطلاعات شانون (Shannon andWeaver., 1949)



 

 

 

نتایج حاصل از آزمون نسبت درست نمایی برای جایگاه McMA2 نشان داد که این جایگاه فقط در جمعیت­های عربی از حالت تعادل هاردی-وینبرگ انحراف دارد. نتایج حاصل از آزمون مربع کای نیز نشان داد که جمعیت گوسفند آرمان در تمامی جایگاه­ها، به‌جز McMA26 انحراف معنی‌داری از تعادل هاردی-وینبرگ دارد (05/0P<). علت این انحرافات را می­توان احتمالاً ناشی از کوچک بودن اندازه نمونه‌ها، انتخاب و مهاجرت دانست. هنگامی‌که در هر جایگاه تمامی جمعیت‌ها باهم در نظر گرفته شد، با هر دو آزمون نیز انحراف معنی‌داری از تعادل هاردی-واینبرگ در بیشتر جایگاه­ها دیده شد (05/0P<)، علت چنین امری را نیز می‌توان تکامل غیر هم‌جهتی دانست که در جمعیت‌های مختلف برای یک جایگاه خاص در طول زمان در اثر تفاوت‌های جغرافیایی مناطق پراکنش آن‌ها روی‌داده است. این امر موجب می‌گردد جایگاه موردنظر برای کل جمعیت‌های ادغام‌شده از تعادل هاردی-واینبرگ انحراف داشته باشد. نرخ بالای جهش در ریزماهواره­ها و ایجاد آلل­های جدید و وجود آلل­های نول در برخی از نشانگرها نیز از عوامل مهم انحراف از تعادل است. چنین مشاهداتی در مطالعات Nanekarani et al. (2010)، Sun et al.  (2010)، Mohammadifar and Mohammadabadi (2011) ، Bhatia and Arora  (2007) و Crispim et al.  (2014) که روی گوسفندان دیگر انجام گردیده است نیز گزارش‌شده است. تنها Buchanan et al.  (1994) هیچ انحرافی از تعادل هاردی- واینبرگ را گزارش نکردند. به‌طورکلی انحراف از تعادل را می­توان به علت وجود آمیزش­های خویشاوندی در بین گونه­ها، نزدیک هم بودن تعداد محدودی از آلل­ها، ناکافی بودن تعداد نمونه­ها و خطای نمونه‌برداری دانست (Chauhan et al., 2007).

مقایسه تعداد، نوع و دامنه اندازه آلل‌های حاصل از مطالعه حاضر تفاوت‌ها و شباهت‌هایی با مطالعات Sun et al.  (2010)، Mohammadifar and Mohammadabadi (2011) ، Bhatia and Arora  (2007)، Crispim et al.  (2014)، Nanekarani et al. (2010) و Esmail-Khanian et al. (2007) دارد که این شباهت­ها و تفاوت­ها را نمی­توان با دلایل خاصی تفسیر نمود زیرا محیط، مواد و روش­های استفاده ‌شده در این مطالعه با مطالعات دیگر متفاوت بوده و همچنین نشانگرهای ریزماهواره­ای­ چندشکلی زیادی دارند که این امر نیز طبیعی است. دامنه اندازه آلل‌ها برای جایگاه‌های McMA2، McMA26، BM6444 و OarHH35 به ترتیب 202-155، 218-184، 177-78 و 164-87 جفت باز به دست آمد. تعداد واقعی آلل (Na) و تعداد آلل مؤثر (Ne) برای ترکیبات مختلف جایگاه-جمعیت، هر جمعیت به ازای تمامی جایگاه‌ها در جدول 3 خلاصه‌ شده است. همان‌طور که در این جدول دیده می‌شود. بیشترین تعداد آلل مشاهده‌شده در جایگاه­های McMA2، McMA26 و OarHH35 به ترتیب مربوط به جمعیت­های دالاق، عربی، لری و قره­گل (12 آلل) و کمترین تعداد از آن جایگاه OarHH35 در جمعیت دالاق بود. کمترین تعداد آلل مؤثر مربوط به جایگاه OarHH35 در جمعیت دالاق (9/5) و بیشترین مقدار آن مربوط به جایگاه McMA26 در جمعیت عربی (3/10) بود، درضمن بدلیل کوچک بودن اندازة نمونه و چند شکلی فراوان ریزماهواره ها نبایستی به دنبال آلل­های  اختصاصی گشت.

در مورد جایگاه McMA2، Mohammadifar and Mohammadabadi (2011)، Nanekarani et al. (2010) و Esmail-Khanian et al. (2007) تعداد آلل یکسان و دامنه آللی تقریباً مشابهی را گزارش کردند، اما Maddox et al. (2000) تعداد 17 آلل را برای این جایگاه گزارش کردند که خیلی بیشتر از تعداد آلل به‌دست‌آمده برای این جایگاه در مطالعه حاضر بود. تعداد آلل مشاهده‌ شده و دامنه آللی مربوط به جایگاه         McMA26  در مطالعه حاضر با گزارش Esmail-Khanian et al. (2007) مشابه بود، اما Qanbari et al. (2007)  در مطالعه­ای که روی نژاد افشاری داشتند، این جایگاه را تک شکل گزارش کردند. در پژوهشی Nanekarani et al. (2010) در بررسی تنوع ژنتیکی گوسفندان نژاد زندی تعداد 11 آلل و Maddox et al. (2000)  تعداد 17 آلل را برای این جایگاه گزارش کردند. در بین جایگاه­های موردمطالعه، این جایگاه کمترین دامنه آللی را نیز داشت (bp34). در مورد جایگاه BM6444، تعداد آلل مشاهده‌شده در این جایگاه متفاوت با مطالعات قبلی بود و نیز با داشتن دامنه آللی 99 جفت باز، بیشترین دامنه را در بین جایگاه­های موردمطالعه داشت. در بررسی تنوع ژنتیکی جمعیت گوسفند بلوچی Esmail-Khanian et al. (2007) تعداد 15 آلل را برای این جایگاه گزارش کردند که تقریباً مشابه با تعداد آلل مشاهده‌شده در گزارش حاضر (13 آلل) بود. در بررسی چندشکلی برخی از نشانگرهای ریزماهواره­ای مرتبط با ژن اینهیبین در گوسفند سنجابی Solimani et al. (2011) تعداد آلل متفاوتی را (4 آلل) نسبت به مطالعه حاضر و مشاهدات قبلی گزارش کردند. جایگاه OarHH35 نسبت به مطالعات گذشته تعداد آلل بیشتری مشاهده شد (12 آلل). در مطالعه Mohammadi and Saberivand (2007) روی تعیین چندشکلی برخی از نشانگرهای ریزماهواره­ای پیوسته به ژن دوقلوزایی FecB در گوسفند قزل، تعداد 9 آلل را در دامنه آللی 130-170 گزارش کردند و Solimani et al. (2011) نیز در بررسی چندشکلی برخی از نشانگرهای ریزماهواره­ای مرتبط با ژن اینهیبین در گوسفند سنجابی، تعداد 3 آلل با دامنه آللی 113-163 جفت باز را برای این جایگاه گزارش کردند.

تنوع درون جمعیتی یا تنوع ژنی به‌صورت هتروزیگوسیتی مورد انتظار نااُریب  Nei در هر جایگاه و برای هر جمعیت در جدول (3) خلاصه ‌شده است. دامنه He نیز بین 83/0-9/0 بود که کمترین مقدار مربوط به جایگاه OarHH35 در جمعیت دالاق و بیشترین مقدار آن مربوط به جایگاه­های McMA26 در جمعیت عربی بود.

هنگامی‌که برای هر جمعیت هتروزیگوسیتی مورد انتظار در سطح تمامی جایگاه‌ها و به‌صورت میانگین هتروزیگوسیتی مورد انتظار جایگاه‌ها، محاسبه شد، بیشترین و کمترین مقدار He به ترتیب برای جمعیت‌های عربی (875/0) و لری (862/0) و با دامنه محدود 013/0 به دست آمد. البته میانگین هتروزیگوسیتی مورد انتظار برای جمعیت قره­گل (87/0) نیز تنها تفاوت اندکی با مقدار حداکثر در جمعیت عربی داشت. دامنه انحراف معیار (St. Dev) هتروزیگوسیتی مورد انتظار نااُریب برای تمامی چهار جایگاه در جمعیت‌های مختلف بین 0132/0و 0493/0 است که حاکی از تفاوت اندک میان آن‌ها است. نتایج، حاکی از تنوع بالا در جمعیت­های مورد مطالعه بود، تمام جایگاه­های موردمطالعه چندشکل بودند درنتیجه خطری از بابت افت هتروزیگوسیتی وجود ندارد. تفاوت اندک ‌بین مقادیر حداکثر و حداقل نشان می‌دهد که تنوع درون جمعیتی برای جمعیت‌های موردمطالعه تقریباً یکسان بود. همچنین می­توان گفت که این جایگاه­ها برای تشخیص آزمون والدینی در جمعیت گوسفندان مناسب­ هستند.

مقادیری که برای He در دو جایگاه McMA2 و McMA26 به دست آمد نزدیک به گزارش Maddox et al.  (2000) بود که به ترتیب مقادیر 9/0 و 88/0 را برای این دو جایگاه به دست آورده‌اند. در مطالعه Nanekarani et al. (2010) روی جایگاه­های McMA2، BM6444 و McMA26 برای هتروزیگوسیتی مورد انتظار به ترتیب مقادیر 91/0، 82/0 و 896/0 را گزارش کردند. مقدار He در مطالعه Sun et al.  (2010) برای جایگاه OarHH35، 66/0 بود که خیلی کمتر از مقدار به‌دست‌آمده در مطالعه حاضر بود. Esmail-Khanian et al. (2007) هتروزیگوسیتی مورد انتظار نااریب Nei را در مورد جایگاه­های McMA2 و BM6444 به ترتیب 9/0 و 93/0 گزارش کردند و نیز Mohammadi and Saberivand (2007) این پارامتر را برای جایگاه OARHH35، 86/0 گزارش کردند.

چون حد نهایی برای هتروزیگوسیتی یک است و از این‌رو مقادیر هتروزیگوسیتی به افزایش تنوع زیاد حساس نمی‌باشند، مقایسه این مقادیر به‌عنوان معیار تنوع درون جمعیتی برای نشانگرهای بسیار چندشکل همچون ریزماهواره‌ها (که در اکثر موارد هتروزیگوسیتی حدود 8/0 یا بالاتر دارند) صحیح نبوده و این تفاوت‌های میان آن‌ها اطلاعات دقیقی را بیان نمی‌کنند. لذا، شاخص اطلاعات Shannon () و محتوای اطلاعات چندشکلی (PIC) در جدول (4) خلاصه‌ شده است.

 

برخلاف هتروزیگوسیتی، که برای هر تعداد آلل، حد نهایی یک را دارد، حداکثر مقدار شاخص شانون () برابر با ln(n) است. باوجودآنکه تفسیر بیولوژیکی این معیار مشخص نیست، ممکن است برای اندازه‌گیری تنوع جایگاه‌های بسیار متغیر همچون ریزماهواره­­ها مفید باشد (Hedrick, 1999; Shannon and Weaver, 1949). نتایج بدست آمده برای شاخص شانون نیز مشابه با آنچه از هتروزیگوسیتی­ها حاصل‌شده بود (هم ازلحاظ روند و هم ازلحاظ بیشترین و کمترین مقدار)، حاصل گردید (جدول 4). بیشترین مقدار شاخص شانون از آن جایگاه McMA26 در جمعیت عربی (396/2) و کمترین مقدار نیز از آن جایگاه OarHH35 در جمعیت دالاق (926/1) بود. در اینجا نیز همچون هتروزیگوسیتی‌ها تفاوت بین جمعیت‌ها به ازای تمامی جایگاه‌ها زیاد نیست که بار دیگر بر یکسان بودن تنوع در جمعیت‌های موردمطالعه دلالت دارد. برای سه جایگاه McMA2، McMA26 و BM6444، Nanekarani et al. (2010) به ترتیب مقادیر 41/2، 33/2 و 89/1 را برای شاخص شانون گزارش کردند. در پژوهشی Kumar et al. (2006) شاخص شانون را 409/1 برای جایگاه OarHH35 گزارش نموده­اند. برای جایگاه OarHH35 در پژوهش Mohammadi and Saberivand (2007)، شاخص شانون به ترتیب 23/2 و 078/2 گزارش شده است. هنگامی‌که  در سطح جایگاه‌ها به ازای تمامی جمعیت‌ها ارزیابی می‌گردد، معیار مناسبی را نیز برای ارزیابی چندشکلی و میزان تغییرپذیری جایگاه‌های موردمطالعه فراهم می‌نماید.

همان‌طور که در جدول 4 مشاهده می­شود، نتایج حاصل از PIC مشابه با آنچه از هتروزیگوسیتی‌ها و شاخص شانون حاصل‌شده بود نیز حاصل گردید. بیشترین مقدار PIC در بین 20 ترکیب جایگاه-جمعیت مربوط به جایگاه­های McMA26 در جمعیت عربی و کمترین مقدار نیز از آن جایگاه OarHH35 در جمعیت­های دالاق و لری (82/0) بود. جایگاه­های چندشکلی همچون نشانگرهای ریزماهواره­ای را در صورت دارا بودن مقادیر PIC بالاتر از 5/0 جزو نشانگرهای دارای محتوای اطلاعاتی زیاد در نظر می­گیرند (Botstein et al. 1980). بررسی معیارهای مختلف چندشکلی نیز حاکی از چندشکلی تمامی جایگاه­های مورد مطالعه بود به‌ طوری‌ که بر اساس تعریف چندشکلی می­توان تمام نشانگرهای تکثیر یافته را 100 درصد چندشکل در نظر گرفت. مقدار PIC به ‌دست ‌آمده برای جایگاه­های McMA2، McMA26 تقریباً مشابه با نتایج Nanekarani et al. (2010) و Esmail-Khanian et al. (2007) بود (9/0). در پژوهشی Esmail-Khanian et al. (2007) برای جایگاه BM6444 مقدار 92/0 را برای PIC  گزارش کردند. در مورد اطلاعات چندشکلی Kumar et al. (2006)، Bhatia and Arora  (2007)و Arora and Bhatia (2006) به ترتیب مقادیر 84/0، 78/0 و 82/0 را برای جایگاه OarHH35 گزارش نموده­اند که مشابه با نتایج به‌دست‌آمده بود. مقادیر زیاد به دست آمده برای ارزش­های PIC در جمعیت­ها و جایگاه­های مورد مطالعه در حدی بود که از جایگاه­های بسیار چندشکل همچون ریزماهواره­ها انتظار می­رود. به‌هرحال، استفاده از آغازگرهای مورد مطالعه در این پژوهش که در سایر مطالعات نیز چندشکلی بالایی را نشان داده­اند، احتمال دستیابی به چندشکلی را در مطالعات جدید افزایش می­دهد.

نتایج به ‌دست‌آمده از برآورد معیارهای تنوع درون­ جمعیتی مانند شاخص شانون، محتوای اطلاعات چندشکلی، تعداد آلل مشاهده ‌شده، تعداد آلل مؤثر و تنوع ژنی Nei نشان‌دهنده تنوع مناسب در جمعیت‌های مورد بررسی بود که می­توان آن‌ها را به ­عنوان یک ذخیره ژنتیکی مناسب برای اهداف مختلف پرورشی و اصلاح نژادی در کشور به ‌حساب آورد و نیز با توجه به انطباق نتایج حاصله با پیش‌فرض‌های معلوم می‌توان جمع‌بندی کرد که نشانگرهای ریزماهواره مورد مطالعه در این پژوهش ابزار بسیار توانمندی برای بررسی و ارزیابی تنوع و روابط ژنتیکی درون و بین جمعیتی و روابط تکاملی میان جمعیت­ها و به ‌طورکلی مطالعات ژنتیک جمعیتی هستند.

 

سپاسگزاری

این پژوهش در آزمایشگاه ژنتیک مولکولی بخش علوم دامی دانشگاه شهید باهنر انجام شد، لذا از همراهی و همکاری تمامی اساتید، دانشجویان و کارکنان گروه علوم دامی این دانشگاه و نیز تمامی مراکز  اصلاح نژاد دام و جهاد کشاورزی سراسر کشور که برای انجام این طرح همکاری کردند، تشکر و قدر دانی می­شود.

 

 

منابع

 

Abadi MRM, Askari N, Baghizadeh A, Esmailizadeh AK (2009). A directed search around caprine candidate loci provided evidence for microsatellites linkage to growth and cashmere yield in Rayini goats. Small Ruminant Research 81: 146-151.

Arora R, Bhatia S (2006). Genetic diversity of Magra sheep from India using microsatellite analysis. Asian-Australasian Journal of Animal Sciences (AJAS) 19: 938 – 942.

Askari N, Abadi MM, Baghizadeh A (2011). ISSR markers for assessing DNA polymorphism and genetic characterization of cattle, goat and sheep populations. Iranian Journal of Biotechnology 9: 222-9.

Banabazi MH, Esmaeil-khanian S, Miraei-Ashtiani SR, Moradi Shahrbabak M (2007). Genetic variation within and between five Iranian sheep populations using microsatellite markers. Iranian Journal of Science Technology Agriculture Nature Research. 10: 4 [In Persian].

Bhatia S, Arora R (2007). Genetic diversity in Kheri-A pastoralists developed Indian sheep using microsatellite markers. Indian Journal of Biotechnology 108-112.

Botstein D, White RL, Skolnick M, Davis RW (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum.Genet 32: 314-331

Buchanan FC, Thue TD (1998). Interbreed polymorphic information content of microsatellites in cattle and sheep. Canadian Journal of Animal Science 78: 425-428.

 

Buchanan FC, Adams LJ, Littlejohn RP, Maddox JF, Crawford AM (1994). Determination of evolutionary relationships among sheep breeds using microsatellites. Genomics 22: 397-403.

Chauhan T, Lal KK, Mohindra V, Singh R, Punia O, Gopalakrishnan A, Sharma PC, Lakra WS (2007). Evaluation genetic differentiation in wild populations of the Indian major carp. Aquaculture 269: 135-149.

Crispim BDA, Seno LDO, Egito AADE, Junior FMDV, Grisolia AB (2014). Application of microsatellite markers for breeding and genetic conservation of herds of Pantaneiro sheep. Electronic Journal of Biotechnology 17: 317–321.

Davis GH, Galloway SM, Ross IK, Gregan SM, Ward J, Nimbkar BV, Ghalsasi PM, Nimbkar C, Subandriyog GD, Inounu I, Tiesnamurti B, Martyniuk E, Eythorsdottir E, Mulsant P, Lecerf F, Hanrahan JP, Bradford GE, Wilson T (2002). DNA test in prolific sheep from eight countries provide new evidence on origin of the booroola FecB mutation. Biology of Reproduction 66:1869-1874.

Esmail-Khanian S, Nejati JA, Afraz F, Daneshyar P, Ghanbari S (2007). Genetic variation amongbaluchi sheep population using microsatellitemarkers. Iranian Journal of Science Technology Agriculture Nature Research 11: 41[In Persian].

Frankham R (1994). Conservation of genetic diversity for animals. 5th world congress on genetics Applied to livestock production, University of Guelph, Ontario, Canada 21: 385-392.

Hedrick PW (1999). Genetic of populations, Second edition. Jones and Bartlett Publishers. Sudbury, MA, USA.

Henry HM, Penty JM, Pierson CA, Crawford AM (1993). Ovine microsatellites at the OARHH 35, OARHH 41, OARHH44, OARHH47and OARHH64 loci. Animal Genetics 24: 222.

Jakaria Zein MSA, Sulandari S, Subandriyo, Muladno (2012). The use of microsatellite markers to study genetic diversity in Indonesian sheep. Journal of the Indonesian Tropical Animal Agriculture (J Indonesian Trop Anim Agric) 371.

Kumar D, Sharma R, Pandey AK, Gour DS, Malik G, Ahlawat SPS, Jain A (2006). Genetic diversity and bottleneck analysis of Indian Bellary sheep by microsatellite markers. Russian Journal of Genetics 43, 9: 996–1005.

Maddox JF, Riffkin CD, Beh KJ (2000). Dinucleotide repeat polymorphism at ovin MCMA1, MCMA2, MCMA5, MCMA8, MCMA9, MCMA11, MCMA14, MCMA20, MCMA24, MCMA26  loci. Animal Genetics 31: 148-149.

Mohammadi G, Saberivand A (2007). Molecular study of polymorphisms of some microsatellite markers linked to the FecB genes in Ghezel Iranian sheep. Iranian Veterinary Journal 2: 34-43 [In Persian].

Mohammadifar A, Mohammadabadi MR (2011). Application of Microsatellite Markers for a Study of Kermani Sheep Genome. Iranian Journal of Animal Science 42, 4: 337-344 [In Persian].

Musthafa Muneeb M, Aljummah RS, Alshaik MA (2012). Genetic diversity of Najdi sheep based on microsatellite analysis. African Journal of Biotechnology 1183: 14868-14876.

Nanekarani Sh, Amirinia C, Mozafari NA, Vaez Torshizi R, Gharahdaghi AA (2010). Genetic Variation among Zandi Sheep Population Using Microsatellite Markers. Veterinary Journal of Islamic Azad University 11: 79-85 [In Persian].

Promega C (1993). DNA Silver Staining System, product DQ7050. Madison, WI, USA.

Qanbari S, Osfoori R, Eskandari Nasab MP (2007). Introgression of major genes into elite breeding flock of afshari sheep: a preliminary evaluation of polymorphism content and applicability of marker data. Iranian Journal of Animal Science 1: 39-47 [In Persian].

Saadat-Noori M, Siah-Mansoor S (1990). Sheep husbandry and management. 4th ed. Ashrafi Publ. Co, Tehran, Iran [In Persian].

Shannon CE, Weaver W (1949). Mathematical theory of communication. 1st Edn. University of Illinois Press, Urbana, IL10: 0252725484.

Shojaei M, Mohammad Abadi MR, Asadi Fozi M, Dayani O, Khezri A, Akhondi M (2010). Association of growth trait and Leptin gene polymorphism in Kermani sheep. Journal of Cell and Molecular Research 2: 67-73.

Solimani B, Chaharaein B, Rahimi Mianji GH (2011). Assessing polymorphism in BM6444, INRA13 and Oarhh35 microsatellite markers associated with inhibin gene in sanjabi sheep. Journal of Studies of Animal Sciences Iran 1: 85-90 [In Persian].     

Sun W, Chang H, Hussein Musa H, Chu M (2010). Study on relationship between microsatellite polymorphism and producing ability on litter size of Hu sheep in China. African Journal of Biotechnology 8704-8711.

Tomasco I, Wlasiuk G, Lessa EP (2002). Evaluation of polymorphism in ten microsatellite loci in Uruguayan sheep flocks. Genetics and Molecular Biology 25: 37-41.

Wang QG, Zhong FG, Li H, Wamg XH, Liu SR, Chen XJ, Gan SQ (2005). Detection of major gene on litter size in sheep. Yi Chuan 271:80-84.

Yeh FC, Yang R, Boyle T (1999). PopGene, Version 1.31. Microsoft window–based freeware for population genetic analysis, University of Alberta. Edmonton, AB, Canada.

Zamani P, Akhondi M, Mohammadabadi MR, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2013). Genetic variation of Mehraban sheep using two intersimple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812-1817.

 


Analysis of genetic diversity in five Iranian sheep population using microsatellites markers

 

Vajed Ebrahimi M.T.*1, Mohammad Abadi M.R.2, Esmailizadeh A.K.3

 

1MSc Student, Department of Animal Science, Shahid Bahonar University of Kerman, Iran.

2Associate Professor, Department of Animal Science, Shahid Bahonar University of Kerman, Iran.

3Professor, Department of Animal Science, Shahid Bahonar University of Kerman, Iran.

 

 

Abstract

Indigenous races in each country are as a national capital and strategic product in the economy and prosperity of the country that maintenance of these races are very valuable. Due to the importance of maintaining diversity and identification of genetic resources in indigenous breeds, five breeds of Iranian sheep population (Arabi, Arman, Dalagh, Karakul and Lory) using four microsatellite markers (McMA2, McMA26, OarHH35 and BM6444 ) was evaluated in terms of genetic diversity. Genomic DNA was extracted from 225 blood samples by optimized salting-out DNA extraction procedure. The PCR reactions were successfully performed with all primers. The results showed that all loci were quite polymorph. Hardy-Weinberg equilibrium test using chi-square test showed that some various combinations of loci-population were in a state of Hardy-Weinberg disequilibrium (P<0.05). The highest number of alleles was observed in loci of McMA2, McMA26 and OarHH35 of Dalagh, Arabic, Lory and Karakul populations respectively (12 alleles) and the lowest number of alleles for OarHH35 in Dalagh population. The maximum expected heterozygosity in combination locus-population of loci McMA2، McMA26، OarHH35، BM6444 were 0.885, 0.9, 0.88 and 0.87 respectively. With respect to the studied loci, it is concluded that all studied sheep populations have wide genetic diversity.

 

Keywords: Sheep, Microsatellite Markers, Genetic Variation, Heterozygosity.

 



*  نویسنده مسئول: محمدتقی واجد ابراهیمی                       تلفن: 03431322689                              Email: m.vajed@yahoo.com

[1] International Society for Animal Genetics (ISAG)

[2] Najdi

[3] Accession Number

* Corresponding Author: Vajed Ebrahimi M.T.                 Tel: 03431322689         Email: m.vajed@yahoo.com

Abadi MRM, Askari N, Baghizadeh A, Esmailizadeh AK (2009). A directed search around caprine candidate loci provided evidence for microsatellites linkage to growth and cashmere yield in Rayini goats. Small Ruminant Research 81: 146-151.
Arora R, Bhatia S (2006). Genetic diversity of Magra sheep from India using microsatellite analysis. Asian-Australasian Journal of Animal Sciences (AJAS) 19: 938 – 942.
Askari N, Abadi MM, Baghizadeh A (2011). ISSR markers for assessing DNA polymorphism and genetic characterization of cattle, goat and sheep populations. Iranian Journal of Biotechnology 9: 222-9.
Banabazi MH, Esmaeil-khanian S, Miraei-Ashtiani SR, Moradi Shahrbabak M (2007). Genetic variation within and between five Iranian sheep populations using microsatellite markers. Iranian Journal of Science Technology Agriculture Nature Research. 10: 4 [In Persian].
Bhatia S, Arora R (2007). Genetic diversity in Kheri-A pastoralists developed Indian sheep using microsatellite markers. Indian Journal of Biotechnology 108-112.
Botstein D, White RL, Skolnick M, Davis RW (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum.Genet 32: 314-331
Buchanan FC, Thue TD (1998). Interbreed polymorphic information content of microsatellites in cattle and sheep. Canadian Journal of Animal Science 78: 425-428.
 
Buchanan FC, Adams LJ, Littlejohn RP, Maddox JF, Crawford AM (1994). Determination of evolutionary relationships among sheep breeds using microsatellites. Genomics 22: 397-403.
Chauhan T, Lal KK, Mohindra V, Singh R, Punia O, Gopalakrishnan A, Sharma PC, Lakra WS (2007). Evaluation genetic differentiation in wild populations of the Indian major carp. Aquaculture 269: 135-149.
Crispim BDA, Seno LDO, Egito AADE, Junior FMDV, Grisolia AB (2014). Application of microsatellite markers for breeding and genetic conservation of herds of Pantaneiro sheep. Electronic Journal of Biotechnology 17: 317–321.
Davis GH, Galloway SM, Ross IK, Gregan SM, Ward J, Nimbkar BV, Ghalsasi PM, Nimbkar C, Subandriyog GD, Inounu I, Tiesnamurti B, Martyniuk E, Eythorsdottir E, Mulsant P, Lecerf F, Hanrahan JP, Bradford GE, Wilson T (2002). DNA test in prolific sheep from eight countries provide new evidence on origin of the booroola FecB mutation. Biology of Reproduction 66:1869-1874.
Esmail-Khanian S, Nejati JA, Afraz F, Daneshyar P, Ghanbari S (2007). Genetic variation amongbaluchi sheep population using microsatellitemarkers. Iranian Journal of Science Technology Agriculture Nature Research 11: 41[In Persian].
Frankham R (1994). Conservation of genetic diversity for animals. 5th world congress on genetics Applied to livestock production, University of Guelph, Ontario, Canada 21: 385-392.
Hedrick PW (1999). Genetic of populations, Second edition. Jones and Bartlett Publishers. Sudbury, MA, USA.
Henry HM, Penty JM, Pierson CA, Crawford AM (1993). Ovine microsatellites at the OARHH 35, OARHH 41, OARHH44, OARHH47and OARHH64 loci. Animal Genetics 24: 222.
Jakaria Zein MSA, Sulandari S, Subandriyo, Muladno (2012). The use of microsatellite markers to study genetic diversity in Indonesian sheep. Journal of the Indonesian Tropical Animal Agriculture (J Indonesian Trop Anim Agric) 371.
Kumar D, Sharma R, Pandey AK, Gour DS, Malik G, Ahlawat SPS, Jain A (2006). Genetic diversity and bottleneck analysis of Indian Bellary sheep by microsatellite markers. Russian Journal of Genetics 43, 9: 996–1005.
Maddox JF, Riffkin CD, Beh KJ (2000). Dinucleotide repeat polymorphism at ovin MCMA1, MCMA2, MCMA5, MCMA8, MCMA9, MCMA11, MCMA14, MCMA20, MCMA24, MCMA26  loci. Animal Genetics 31: 148-149.
Mohammadi G, Saberivand A (2007). Molecular study of polymorphisms of some microsatellite markers linked to the FecB genes in Ghezel Iranian sheep. Iranian Veterinary Journal 2: 34-43 [In Persian].
Mohammadifar A, Mohammadabadi MR (2011). Application of Microsatellite Markers for a Study of Kermani Sheep Genome. Iranian Journal of Animal Science 42, 4: 337-344 [In Persian].
Musthafa Muneeb M, Aljummah RS, Alshaik MA (2012). Genetic diversity of Najdi sheep based on microsatellite analysis. African Journal of Biotechnology 1183: 14868-14876.
Nanekarani Sh, Amirinia C, Mozafari NA, Vaez Torshizi R, Gharahdaghi AA (2010). Genetic Variation among Zandi Sheep Population Using Microsatellite Markers. Veterinary Journal of Islamic Azad University 11: 79-85 [In Persian].
Promega C (1993). DNA Silver Staining System, product DQ7050. Madison, WI, USA.
Qanbari S, Osfoori R, Eskandari Nasab MP (2007). Introgression of major genes into elite breeding flock of afshari sheep: a preliminary evaluation of polymorphism content and applicability of marker data. Iranian Journal of Animal Science 1: 39-47 [In Persian].
Saadat-Noori M, Siah-Mansoor S (1990). Sheep husbandry and management. 4th ed. Ashrafi Publ. Co, Tehran, Iran [In Persian].
Shannon CE, Weaver W (1949). Mathematical theory of communication. 1st Edn. University of Illinois Press, Urbana, IL10: 0252725484.
Shojaei M, Mohammad Abadi MR, Asadi Fozi M, Dayani O, Khezri A, Akhondi M (2010). Association of growth trait and Leptin gene polymorphism in Kermani sheep. Journal of Cell and Molecular Research 2: 67-73.
Solimani B, Chaharaein B, Rahimi Mianji GH (2011). Assessing polymorphism in BM6444, INRA13 and Oarhh35 microsatellite markers associated with inhibin gene in sanjabi sheep. Journal of Studies of Animal Sciences Iran 1: 85-90 [In Persian].     
Sun W, Chang H, Hussein Musa H, Chu M (2010). Study on relationship between microsatellite polymorphism and producing ability on litter size of Hu sheep in China. African Journal of Biotechnology 8704-8711.
Tomasco I, Wlasiuk G, Lessa EP (2002). Evaluation of polymorphism in ten microsatellite loci in Uruguayan sheep flocks. Genetics and Molecular Biology 25: 37-41.
Wang QG, Zhong FG, Li H, Wamg XH, Liu SR, Chen XJ, Gan SQ (2005). Detection of major gene on litter size in sheep. Yi Chuan 271:80-84.
Yeh FC, Yang R, Boyle T (1999). PopGene, Version 1.31. Microsoft window–based freeware for population genetic analysis, University of Alberta. Edmonton, AB, Canada.
Zamani P, Akhondi M, Mohammadabadi MR, Saki AA, Ershadi A, Banabazi MH, Abdolmohammadi AR (2013). Genetic variation of Mehraban sheep using two intersimple sequence repeat (ISSR) markers. African Journal of Biotechnology 10: 1812-1817.