The Analysis of IGF-I Gene Polymorphism and its Association with Breeding Values of Production and Reproduction Traits in Raieni Cashmere Goats

Document Type : Research Paper

Authors

1 M. Sc. Student, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran.

2 2Associate professor, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran.

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

4 4Ph.D Student, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran.

Abstract

The Insulin-like growth factor I (IGF-I) has a similar structure with Insulin hormone. This hormone which is transcribed by (IGF-I) gene has an important functional in growth and development of different body tissues. The aim of this study was to investigate the polymorphism of this gene and its relationship with the estimated breeding values of the birth type, birth weight and fleece weight in Raieni Cashmere goats. Breeding values were estimated using records on 13020 Raieni cashmere goats originated from 336 sire and 3617 dams. Univariate animal model was applied for the genetic analysis of the investigated traits using ASReml. To study IGF-1 gene polymorphism, 94 Raieni cashmere goats were selected based on their estimated breeding values for these traits. PCR was used for amplification of 363 bp fragment of exon 4 of IGF-I gene. Thereafter the animal’s genotype was determined using PCR-RFLP. Three genotypes including AA, AB and BB were observed. The frequency of these genotypes was 0.22, 0.65 and 0.13, respectively. No significant associations were found between these genotypes and the investigated production and reproduction traits. The results of this study show that, there is important genetic diversity for the IGF-I in the investigated population but to study the association between the polymorphism of this gene and the studied traits, more data should be used as currently few animals were genotyped.  

Keywords


ارتباط چند شکلی بخشی از ژن IGF-I با ارزش اصلاحی صفات تولیدی و تولید مثلی در بز کرکی راینی

محمد حسین ملایی1، مسعود اسدی فوزی*2، محمدرضا محمدآبادی3، مهدیه منتظری4 و5

1 دانش آموخته کارشناسی ارشد، گروه علوم دامی، دانشکده کشاورزی، دانشگاه شهید باهنر کرمان، کرمان، ایران

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

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

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

5 انجمن پژوهشگران جوان، دانشگاه شهید باهنر کرمان، کرمان، ایران

 

چکیده

فاکتور رشد شبه انسولین یک (IGF-I) ساختار مشابهی با هورمون انسولین دارد. این هورمون که به وسیله ژن IGF-I کد می­شود نقش مهمی در رشد و نمو بافت­های مختلف بدن دارد. این تحقیق به منظور بررسی چند شکلی این ژن و ارتباط آن با ارزش اصلاحی صفات تیپ تولد، وزن تولد و وزن الیاف در بز کرکی راینی انجام شد. بدین منظور از رکوردهای 13020 حیوان مربوط به 336 پدر و 3617 مادر استفاده گردید. جهت آنالیز ژنتیکی صفات مورد بررسی از مدل حیوانی یک متغیره و نرم­افزار ASReml استفاده شد. جهت بررسی اثر چند شکلی ژن IGF-I بر روی ارزش اصلاحی برآورد شده برای صفات مورد بررسی تعداد 94 راس از گله ایستگاه پرورش و اصلاح نژاد بز کرکی راینی انتخاب و خون­گیری بعمل آمد. واکنش زنجیره­ای پلیمراز جهت تکثیر قطعه 363 جفت بازی اگزون 4 ژن IGF1 انجام گرفت. سپس ژنوتیپ حیوانات با استفاده از روش PCR-RFLP تعیین شد. در این تحقیق سه ژنوتیپ AA، AB و BB با فراوانی­های ژنوتیپی 22/0، 65/0 و 13/0 مشاهده گردید، اما بین این ژنوتیپ­ها و صفات تولیدی و تولید مثلی مورد بررسی رابطه معنی­داری مشاهده نشد. نتایج حاصل از این تحقیق نشان می­دهد که در جمعیت مورد بررسی از نظر جایگاه ژنI  IGF- تنوع قابل ملاحظه­ای وجود دارد اما با توجه به تعداد کم نمونه­های مورد استفاده در این تحقیق، پیشنهاد می­گردد جهت بررسی ارتباط بین این تنوع ژنتیکی و صفات مورد بررسی از نمونه­های بیشتری استفاده گردد. 

واژه­های کلیدی: ژن IGF-I، ارزش اصلاحی، PCR-RFLP، بز کرکی راینی.



مقدمه

امروزه با بکارگیری برنامه­های اصلاح نژادی کیفیت و کمیت برخی صفات مهم اقتصادی دام­ها از نسلی به نسل بعد بهبود یافته است. در طراحی این برنامه­های اصلاح نژادی پارامترهای ژنتیکی صفات هدف و هم­چنین صفاتی که به عنوان ملاک انتخاب دام­های برتر مورد استفاده قرار می­گیرند، نیاز می­باشد. دقت برآورد این پارامترهای ژنتیکی در میزان کارایی برنامه­های اصلاح نژادی موثر می­باشند (Asadi Fozi, 2005).

یکی از راه­های افزایش تولیدات دام، شناخت پتانسیل­های ژنتیکی آنها و انتخاب دام­های با پتانسیل ژنتیکی برتر می­باشد (Jack & Felix, 1996). در حوزه ژنتیک و اصلاح دام، اطلاع از ساختار ژنتیکی جمعیت­ها می­تواند کمک بزرگی جهت برنامه­ریزی برای طرح­های اصلاح نژادی و از همه مهمتر، حفظ ذخایر ژنتیکی باشد. روش­های مولکولی و استفاده از نشانگرهای مولکولی در این زمینه یکی از بهترین گزینه­ها به حساب می­آید، زیرا با توجه به اطلاعات زیادی که به دست می­دهد می­تواند نتایجی که از تجزیه و تحلیل رکوردها با روش­های آماری به دست آمده است را تأیید و تکمیل نموده و حتی ممکن است که آنها را رد کند (Alinaghizadeh et al., 2010). به علاوه، استفاده از ژنتیک ملکولی فواید زیادی دارد که یکی از این فواید تعیین ژنوتیپ افراد برای یک جایگاه خاص می­باشد (Mousavizadeh et al., 2009) هم­چنین استفاده از نشانگرهای ژنتیکی در انتخاب و اصلاح نژاد حیوانات ممکن است به طور مهیجی پیشرفت ژنتیکی را تسریع کند (Javanmard et al., 2008). از طرفی مطالعه تنوع ژنتیکی نژادهای بومی برای حفاظت از منابع ژنتیکی ذخایر بومی لازم و ضروری است (Mohammadi et al., 2009). حفاظت باید بر اساس دانش عمیقی از منابع ژنتیکی نژادهای خاص باشد، لذا تلاش برای شناسایی و تعیین خصوصیات ژنتیکی نژادهای بومی و محلی بسیار اهمیت دارد (Shojaei et al., 2010; Zamani et al., 2013).

ژنIGF-I  روی کروموزوم شماره 4 بز قرار دارد و دارای 6 اگزون و 5 اینترون می­باشد. فاکتور رشد شبه انسولین (IGF-I) در بسیاری از فرآیندهای فیزیولوژیکی بدن مثل رشد، شیردهی، تولید مثل و سیستم ایمنی نقش دارد (Walkden-brown et al., 2008). هورمون رشد یکی از عواملی است که ممکن است سبب افزایش سنتز IGF-1 در کبد شود و تولید موضعی این هورمون را در دیگر بافت­ها تحریک کند (Abdolmohammadi, 2008). ژنIGF-I  عامل افزایش دهنده جذب گلوکز در بافت­های محیطی بوده که باعث تحریک ساخت گلیکوژن شده و دارای اثری شبیه انسولین می­باشد که با افزایش جذب اسید آمینه منجر به ساخت پروتئین می­شود (Roite et al., 2001). شواهد نشان می­دهد که جهش در IGF-I با غلظت هورمون IGF-I، صفات لاشه، صفات رشد و تولید مثل در برخی گونه­ها نظیر گوسفند، گاو و بز ارتباط دارد (Thomas et al., 2007; Velazquez et al., 2005; Kim et al., 2009; Karimi Kurdistani et al., 2012). در گاوهای گوشتی آنگوس و سمینتال مشخص گردید که اثر ژن IGF-1 بر وزن بدن و ضخامت چربی پشت معنی­دار است (Curia et al., 2005). جهش در اینترون شماره 4 این ژن در بز با وزن شش ماهگی، یک سالگی و افزایش وزن روزانه از سه تا دوازده ماهگی همبستگی معنی­داری دارد (Karimi Kurdistani et al., 2012).

حدود 30 میلیون راس بز کرکی در سراسر جهان وجود دارد که 5/4 تا 5 میلیون راس از آن­ها (حدود 20 درصد بزهای جهان) در ایران پرورش داده می­شوند (Baghizadeh et al., 2009). اگر چه پژوهش­های زیادی روی بزهای کرکی راینی انجام شده است (Alinaghizadeh et al., 2010; Askari et al., 2011; Askari et al., 2009; Askari et al., 2010; Moghadaszadeh et al., 2015; Mohammadabadi et al., 2009; Shamsaddini et al., 2016; Askari et al., 2008; Hassani et al., 2010; Mohammadabadi, 2012; Tohidi nezhad et al., 2015)، ولی تاکنون همبستگی ژن IGF-I و ارزش اصلاحی صفات مهم مطالعه نشده است. لذا، هدف از تحقیق حاضر، بررسی چندشکلی ژن IGF-I و ارتباط آن با ارزش اصلاحی صفات تیپ تولد، وزن تولد و وزن الیاف در بز راینی بود.


مواد و روش­ها

به منظور برآورد ارزش اصلاحی برخی صفات تولیدی و تولید مثلی بز کرکی راینی در مجموع تعداد 10986 رکورد مربوط به صفت تیپ تولد از 1270 بز دو قلو و 4228 بز تک قلو، 13020 رکورد مربوط به صفت وزن تولد و 7651 رکورد مربوط به صفت وزن الیاف بین سال­های 1370 تا 1390 مورد استفاده قرار گرفت. لازم به ذکر است که داده­های مربوط به صفت وزن الیاف مربوط به 5 سال متفاوت می­باشند. ساختار داده­ها در جدول 1 موجود می­باشد.

برای آنالیز ژنتیکی صفات مورد بررسی شامل تیپ تولد، وزن تولد و وزن الیاف و برآورد ارزش اصلاحی آنها از مدل حیوانی یک متغیره استفاده شد (مدل 1). لازم به توضیح است که در مدل یک متغیره­ای که برای هر کدام از صفات استفاده گردید اثرات ثابت و تصادفی مهم که مختص همان صفت بود، گنجانده و مابقی حذف گردیدند. برای سهولت انجام آنالیزها از مدل یک متغیره استفاده گردید. برای اجرای این مدل­ها نرم­افزارv3.1  ASReml بکار گرفته شد (Glimur et al., 2009).

[1]   y = Xb + Z1u + Z2 m + Z3mpe + Z4 pe + e

 


جدول 1- ساختار داده­های استفاده شده در این مطالعه.

Table 1- Data structure used in this study.

صفات (Traits)

تعداد حیوان

Number of animal

تعداد رکورد

Number of record

تعداد پدر

Number of sire

تعداد مادر

Number of dam

حداقل

Minimum

حداکثر

Maximum

انحراف معیار

Standard deviation

وزن تولد

(Birth weight)

13020

13020

336

6317

0.5

3.8

0.39

تیپ تولد

 (Birth type)

5498

10986

307

2017

1

2

0.33

وزن الیاف

(Fleece weight)

7651

19948

380

3296

0.35

0.97

0.13

 

 

در این مدل، y بردار مشاهدات؛ b بردار اثرات ثابت؛  uبردار اثرات تصادفی ژنتیک مستقیم حیوان؛ m بردار اثرات ژنتیکی مادری، mpe بردار اثرات محیط دائمی مادر، pe بردار اثرات محیط دائمی حیوان و e بردار اثرات تصادفی اشتباه می­باشند. هم­چنین X، 1Z، Z2، Z3 و Z4 به ترتیب ماتریس ضرایب برای اثرات ثابت، اثرات ژنتیکی حیوان، اثرات ژنتیکی مادری، اثرات محیط دائمی مادری و اثرات محیط دائمی حیوان می­باشند. با توجه به اینکه برای صفات تیپ تولد و وزن الیاف رکورد­ سال­های مختلف (رکورد­های تکرار شده) وجود داشت برای آنالیز ژنتیکی این صفات از مدل تکرار­پذیری استفاده شد. به­ عبارت دیگر، در این مدل از اثرات محیط دائمی حیوان نیز استفاده شد اما برای آنالیر ژنتیکی وزن تولد این اثر از مدل ]1 [حذف گردید.

در مدل آنالیز ژنتیکی صفات مورد بررسی مدل ]1 [اهمیت هرکدام از اثرات تصادفی اضافی شامل اثرات ژنتیکی مادری و محیط دائمی مادری با استفاده از آزمون لگاریتم حداکثر درستنمایی[1] مورد بررسی قرار گرفت. بر این اساس اثرات مهم در مدل گنجانده و اثراتی که معنی­دار نبودند از مدل آنالیز ژنتیکی حذف شدند. در مدل آنالیز ژنتیکی تیپ تولد، اثر گله، سال زایش مادر، شکم زایش مادر و فصل زایش مادر به عنوان اثرات ثابت در نظر گرفته شدند. هم­چنین در مدل آنالیز ژنتیکی وزن الیاف اثرات گله، جنس، سن مادر، تیپ تولد و سال و ماه کرک چینی به عنوان اثرات ثابت لحاظ شدند. اثرات ثابت مورد استفاده در مدل آنالیز ژنتیکی وزن تولد عبارت از اثر گله، جنس، تیپ تولد، سال و ماه تولد و سن مادر بر اساس ماه بودند. لازم به ذکر است که برای تمام اثرات ثابت مورد استفاده در مدل تصحیح صورت گرفته است.

به منظور بررسی ارتباط چندشکلی ژن IGF-I با ارزش اصلاحی هر کدام از صفات تیپ تولد، وزن تولد و وزن الیاف از بین بزهای کرکی راینی موجود در ایستگاه، تعداد 94 راس بز انتخاب و خون­گیری شدند. به این منظور ابتدا دام­های موجود بر اساس ارزش اصلاحی تیپ تولدشان به سه گروه مختلف شامل حیوانات دارای ارزش اصلاحی بالا، متوسط و پایین تقسیم گردیدند. در مرحله بعد از هر کدام از این گروه­ها تعداد حدود 30 نمونه انتخاب شد. حیواناتی ­که برای تیپ تولد در گروه ارزش اصلاحی بالا قرار داشتند برای وزن الیاف و وزن بدن در گروه ارزش اصلاحی پائین قرار گرفتند و بالعکس. به دلیل اینکه بین ارزش­های اصلاحی وزن بدن و وزن الیاف همبستگی مثبت و بین ارزش اصلاحی این دو صفت و ارزش اصلاحی تیپ تولد همبستگی منفی وجود داشت.  

DNA ژنومی با استفاده از کیتDiatum DNA  استخراج گردید. از آنجائی­که استخراج DNA با کیفیت و خلوص مطلوب، شرط لازم برای به دست آوردن تکرارپذیری بالا برای بیشتر نشانگرهای مولکولی است (Bechmann & Soller, 1987)، پس از استخراج، کیفیت و کمیت DNA با دو روش اسپکتروفتومتری و الکتروفورز بر روی ژل آگارز 1 درصد تعیین گردید. سپس با استفاده از آغازگرهای طراحی شده توسط شرکت ژن فن آوران، یک قطعه 363 جفت بازی از اگزون 4 ژن فاکتور رشد شبه انسولین به وسیله واکنش زنجیره پلیمراز تکثیر شد. انتخاب توالی آغازگرها بر اساس مطالعه Liu et al. (2014) صورت گرفت و توالی آنها به شرح زیر می­باشد:

F: 5'-CACAGCGTATTATCCCAC-3'

R: 5'-GACACTATGAGCCAGAAG-3'

 در مرحله بعد، واکنش زنجیره­ای پلیمراز در حجم 25 میکرولیتر طی 35 سیکل انجام شد. سیکل حرارتی PCR به صورت مرحله واسرشته سازی اولیه در °C95 به مدت 5 دقیقه، واسرشته سازی در °C94 به مدت 40 ثانیه، اتصال در °C 55 به مدت 30 ثانیه، تکثیر در °C72 به مدت 35 ثانیه و تکثیر نهایی در °C72 به مدت 10 دقیقه انجام شد. در مرحله بعد، نمونه­های سنتز شده با آنزیم برشی HaeIIIساخت شرکت Fermentase برش و روی ژل آگاروز 2 درصد جهت تعیین ژنوتیپ افراد الکتروفورز شدند. برای خواندن آلل­ها از نرم­افزار OneD Scan نسخه یک استفاده شد (Zamani et al., 2015). سپس به منظور تجزیه و تحلیل معیارهایی از قبیل فراوانی آللی، تعداد آلل مشاهده شده، تعداد آلل موثر، هتروزیگوسیتی (He) و شاخص شانون (I) از نرم­افزار PopGene (Yeh et al., 1999) استفاده گردید.

در پایان نیز به منظور بررسی ارتباط ژنوتیپ­ها با ارزش اصلاحی صفات وزن تولد، تیپ تولد و وزن الیاف مدل آماری [2] مورد استفاده قرار گرفت. این مدل به وسیله نرم­افزارv3.1  ASReml آنالیز گردید (Glimur et al., 2009).

[2] 

در مدل [2] برابر ارزش اصلاحی (برای هر کدام از صفات) می­باشد. µ برابر با میانگین، برابر با اثر i امین ژنوتیپ و  برابر با مقادیر باقیمانده می­باشد.

 

نتایج و بحث

در تحقیق حاضر ارتباط بین چند شکلی ژن IGF-I با ارزش اصلاحی صفات تیپ تولد، وزن تولد و وزن الیاف مورد بررسی قرار گرفت. برای برآورد ارزش اصلاحی این صفات از مدل حیوانی یک متغیره استفاده گردید. در مدل آنالیز ژنتیکی ابتدا اهمیت اثرات ثابت مورد بررسی قرار گرفت. بر این اساس گله و سال تولد بر هر سه صفت مورد مطالعه اثر معنی­دار داشتند (001/0 p<). اثر سن مادر (شکم زایش) بر تیپ تولد و وزن تولد معنی­دار بود (01/0p<) که این می­تواند به دلیل افزایش جثه حیوان در نتیجه بهبود شرایط رحم در سنین بالاتر باشد (Chungyan et al., 2009). هم­چنین جنسیت بر وزن تولد و وزن الیاف اثر معنی­داری نشان داد (001/0 p<). به گونه­ای که حیوانات نر دارای وزن تولد و وزن الیاف بیشتری نسبت به ماده­ها بودند که می­تواند ناشی از تفاوت واکنش­های فیزیولوژی و ژنتیکی در جنس نر و ماده، از قبیل هورمون­ها باشد (Hermiz et al., 2009).

محققین در سال 2008 بیان کردند که وزن تولد بزهای نژاد Draa برای حیوانات نر 19 درصد بیشتر از حیوانات ماده می­باشد (Boujenane & Hazzab, 2008). نتایج حاصل از تحقیق حاضر با نتایج گزارش شده برای نژادهای بز Boer، Barbari و بومی عراقی مطابقت دارد (Chungyan et al., 2008; Bharathidhasan et al., 2009; Hermiz et al., 2009). اختلاف وزن الیاف نرها و ماده­ها در بزهای کرکی یکساله استرالیا 51 گرم گزارش شده است (Restall & Pattie, 1989). در مطالعه­ای دیگر میانگین حداقل مربعات وزن الیاف در بزهای نر آنقوره استرالیا در سنین مختلف بیشتر از ماده­ها گزارش گردید (Walkden-brown et al., 2008). در مطالعه­ای دیگر بیان شده است که بزهای نر آنقوره فرانسه نسبت به ماده­ها 250 گرم الیاف بیشتری تولید می­کنند (Allain & Roguet, 2003). در تحقیقی که امامی میبدی بر روی بزهای کرکی راینی انجام داد، اختلاف بین وزن بیده را در دو جنس را 62/66 گرم گزارش نموده است (Emami Meybodi, 1993). در تحقیقی دیگر متوسط وزن الیاف در بزهای نر نژاد راینی حدود 40 گرم بیشتر از ماده­ها گزارش شده است (Mohebi Nejad & Asadi Fozi, 2012).

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

به منظور بررسی چند شکلی ژن IGF-I، طی واکنش PCR، قطعه 363 جفت بازی اگزون 4 ژن فاکتور رشد شبه انسولین (IGF-I) تکثیر شد. سپس الگوهای ژنوتیپی حاصل از PCR-RFLP بیانگر 3 الگوی تنوع این جایگاه در بز کرکی راینی بود که فراوانی آنها محاسبه گردید (جدول 2). در این تحقیق سه ژنوتیپ AA (bp 363)، AB (bp 363، bp 264) و BB (bp 264) مشاهده گردید (شکل 1). سودمندی یک نشانگر ژنتیکی به تعداد و فراوانی آللی آن مربوط می­شود که این فاکتورها در میزان هتروزیگوسیتی و چندشکلی آن تاثیر دارند. با توجه به نتایج جدول 2، بیشترین فراوانی ژنوتیپی و آللی به ترتیب مربوط به ژنوتیپ AB و آلل A می­باشد. این الگوهای ژنوتیپی متفاوت بین افراد درون یک نژاد، بیانگر وجود جهش­های تک نوکلئوتیدی است که منجر به بروز تفاوت­های ژنتیکی بین افراد درون یک نژاد شده است. فراوانی سه الگوی باندی در گوسفند نژاد ماکوئی به ترتیب 06/0، 42/0 و 52/0 گزارش شده است (Haji Hossenloo et al., 2011). هم­چنین ژنوتیپ­های مربوط به ژن IGF1 در گوسفندان زل به ترتیب با فراوانی­های 67/56، 67/36 و 67/6 درصد مشاهده شدند (Ghazikhani et al., 2012). با استفاده از تکنیک SSCP یک موتاسیون تک نوکلئوتیدی در ژن IGF-I گزارش شده است (Li et al., 2006). هم­چنین دو جهش تک نوکلئوتیدی در پروموتور ژن IGF-I در گوسفندان خالص Polypay گزارش گردیده است (Yilmaz et al., 2005).

 

 

 

شکل 1- هضم آنزیمی قطعه 363 جفت بازی با استفاده از آنزیم HaeIII

Figure 1- Enzyme digestion of 363 bp fragment using HaeIII.

 


جدول 2- فراوانی ژنوتیپی و آللی مرتبط با ژن IGF-I در بز کرکی راینی.

Table 2- Genotype and allele frequency of IGF-I gene in Rayeni Cashmere goats.

ژنوتیپ

Genotype

فراوانی ژنوتیپی

Genotypic frequency

آلل

Allele

فراوانی آللی

Allelic frequency

BB

0.13

A

0.54

AB

0.65

B

0.46

AA

0.22

 

 

 

 

نتایج حاصل از این مطالعه نشان داد که جمعیت مورد مطالعه از تعادل هاردی-وینبرگ انحراف دارد (05/0 p<). دلایل انحراف از تعادل هاردی-وینبرگ را می­توان به علت انتخاب حیوانات در ایستگاه بر اساس شکل ظاهری و خطای نمونه­گیری دانست. در جمعیت مورد مطالعه دو معیار شاخص شانون و هتروزیگوسیتی مورد انتظار به ترتیب 69/0 و 5/0 برآورد گردید که بیانگر تنوع قابل قبول در این جمعیت می­باشد (جدول 3). با توجه به این نکته که تنوع اساس کارهای اصلاح نژادی است، این جمعیت می­تواند منبع بسیار خوبی برای انجام کارهای اصلاح نژادی و طراحی تلاقی­های مناسب جهت حفظ این نژاد باشد.

 

 

جدول 3- معیارهای تنوع ژنتیکی حیوانات مورد بررسی.

Table 3- The criteria diversity of all animals

هتروزیگوسیتی مشاهده شده

Observation Heterozygosity

 

هتروزیگوسیتی مورد انتظار

Expected Heterozygosity

تعداد آلل مشاهده شده

Observed number of alleles

تعداد آلل موثر

Effective number of alleles

شاخص شانون

Shannon Index

شاخص نئی

Nei Index

0.67

0.50

2

1.9

0.69

0.50

 

 

در این تحقیق ارتباط بین چند شکلی ژن IGF-I و ارزش­های اصلاحی برآورد شده مورد بررسی قرار گرفت. نتایج حاصل نشان می­دهد اگر چه میانگین عددی ارزش اصلاحی تخمین زده شده برای صفات تیپ تولد، وزن الیاف و وزن تولد در ژنوتیپ­های مختلف یکسان نمی­باشد با این حال از نظر آماری معنی­دار نیستند (p>0.05) (جدول 4). با توجه به مقادیر انحراف معیار (SE) این مساله می­تواند به دلیل کوچک بودن جمعیت مورد مطالعه باشد. میانگین حداقل مربعات ارزش اصلاحی تیپ تولد، وزن تولد و وزن الیاف در هرکدام از سه ژنوتیپ ذکر شده در جدول 4 آورده شده است.

 

ارتباط ارزش اصلاحی صفات با انواع ژنوتیپ­ها

بین ژنوتیپ­های مشاهده شده و ارزش اصلاحی برآورد شده حیوانات برای صفت تیپ تولد رابطه معنی­داری مشاهده نشد. اگر چه بین ژنوتیپ مشاهده شده و ارزش اصلاحی صفات ارتباط معنی­داری تشخیص داده نشد با این حال مشاهده شد که افراد دارای آلل A دارای میانگین ارزش اصلاحی بالاتری نسبت به میانگین ارزش اصلاحی افراد دارای آلل B برای صفت تیپ تولد می­باشند.

بین ارزش اصلاحی صفت وزن تولد و ژنوتیپ­های مشاهده شده رابطه معنی­داری وجود نداشت. با این حال مشاهده شد که افراد دارای آلل B دارای میانگین ارزش اصلاحی بالاتری برای صفت وزن تولد نسبت به میانگین ارزش اصلاحی افراد دارای ژنوتیپ A بودند.

رابطه معنی­داری بین ژنوتیپ­های مشاهده شده و ارزش اصلاحی صفت وزن الیاف وجود ندارد. با این حال مشاهده شد که حیواناتی که تنها دارای آلل A می­باشند دارای میانگین ارزش اصلاحی منفی برای صفت وزن الیاف بودند و افرادی که تنها دارای آلل B می­باشند دارای میانگین ارزش اصلاحی مثبتی برای صفت وزن الیاف می­باشند. اما، افرادی که هم دارای آلل A و هم آلل B در جایگاه مورد نظر می­باشند دارای میانگین ارزش اصلاحی بالاتری نسبت به دیگر ژنوتیپ­ها می­باشند که می­تواند به دلیل اثر هتروزیس در این ژنوتیپ باشد.

 

 

جدول 4- میانگین حداقل مربعات ارزش اصلاحی ژنوتیپ­های مورد بررسی همراه با خطای معیار برای صفات مورد بررسی.

Table 4- Least-square mean and standard error of estimated breeding values of the different genotypes for the studied traits.

 

ژنوتیپ

Genotype

AA

AB

BB

تیپ تولد (birth type)

0.0036±0.0001

0.0079±-0.0022

0.0097±-0.0007

وزن تولد (birth weight)

0.0307±0.0111

0.0663±0.0200

0.0817±0.0842

وزن الیاف (fleece weight)

0.0003±-0.0003

0.0190±0.0107

0.0013±0.0013

 

 

در پژوهشی  Thomas et al.(2007)، اثر معنی­دار چند شکلی در ناحیه پروموتور ژن IGF-I را بر وزن بدن گوسفند نشان دادند. تحقیقات گذشته نشان می­دهد بین میزان غلظت IGF-I در سرم خون و رشد الیاف در بز کرکی ارتباط معنی­داری وجود دارد (Ma et al., 2010). جهش در ژن IGF-I با صفات تولید مثلی در بز از قبیل تیپ تولد و نرخ آبستنی در اولین لقاح همبسته می­باشد (Patton et al., 2007). در تحقیقی دیگر گزارش شده است که جهش در ژن IGF-I در بز نژاد Guliuma تیپ تولد را تحت تاثیر قرار می­دهد (Ping-qing et al., 2011).   Velazquez et al. (2008) بیان کردند می­توان از IGF-I به عنوان یک نشانگر برای انتخاب گاوهای با پتانسیل باروری بالا بخصوص در گله­هایی که دارای شایستگی بالا برای تولید شیر می­باشند، استفاده نمود.

 نتایج حاصل جدول 4 نشان می­دهد که در هر کدام از این صفات بین میانگین­ها تفاوت معنی­داری وجود ندارد. بعبارت دیگر بین چند شکلی این ژن و صفات مورد بررسی ارتباط معنی­داری وجود ندارد. بالا بودن خطای معیار این میانگین­ها که خود نتیجه تعداد کم نمونه­های مورد بررسی است می­تواند در نتیجه حاصل موثر باشد. با مقایسه میانگین ارزش اصلاحی هر صفت در ژنوتیپ­های مختلف نیز یک روند خاص حاصل می­شود بطوریکه به نظر می­رسد آلل A ژن IGF-I یک آلل مناسب برای صفت تیپ تولد و آلل B این ژن یک آلل مناسب برای صفات وزن تولد و وزن الیاف در جمعیت مورد مطالعه باشد. بعبارت دیگر ممکن است با افزایش تعداد نمونه­ها نتیجه متفاوت و دقیق­تری بدست آورد.

 

نتیجه گیری

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


منابع

Abdolmohammadi A (2008). Study of genetic variation for four candidate genes using PCR-RFLP and HRM and their association with reproduction and production traits in Holstein cows of Iran. Thesis of Ph.D. University of Tehran, Iran  (In Farsi).

Alinaghizadeh H, Mohammadabadi MR, Zakizadeh S (2010). Exon 2 of BMP15 gene polymorphismin Jabal Barez red goat. Journal of Agricultural Biotechnology 2: 69-80 (In Farsi).

Allain D, Roguet J (2003). Genetic and nongenetic factors influencing mohair production traits within the national selection scheme of Angora goats in France. Livestock Production Science 82: 129-137.

Asadi Fozi M (2005). Appropriate models for genetic analysis and efficient selection of Merino sheep. Ph.D. Thesis. University of New England, Armadale, Australia.

Askari N, Baghizadeh A, Mohammadabadi MR (2008). Analysis of the genetic structure of Iranian indigenous Raieni Cashmere goat populations using microsatellite markers.Biotechnology 2: 1-4 (In Farsi).

Askari N, Baghizadeh A, Mohammadabadi MR (2010). Study of genetic diversity in four populations of Raeini Cashmere goat using ISSR markers. Modern Genetics Journal 5: 49-56 (In Farsi).

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

Askari N, Mohammadabadi MR, Beygi Nassiry MT, Baghizadeh A, Fayazi J (2009). Study of genetic diversity of Raeini Cashmere goat based on microsatellite markers. Journal of Agricultural Science 18: 155-161 (In Farsi).

Baghizadeh A, Bahaaddini M, Mohamadabadi MR, Askari N (2009). Allelic variations in exon 2 of caprine MHC class II DRB3 gene in Raeini Cashmere goat. American-Eurasian Journal of Agricultural and Environmental Science 6: 454-459.

Bechmann JS, Soller M (1987). Molecular markers in the genetic improvement of farm animals. Biotechnology 5: 573-576.

Bharathidhasan A, Narayanan R, Gopu P, Subramanian A, Prabakaran R, Rajendran R (2009). Effect of nongenetic factors on birth weight, weaning weight and pre weaning gain of Barbari goat.  Journal of Veterinary & Animal Science 5: 99-103.

Boujenane I, Hazzab AE (2008). Genetic parameter for direct and mater effect on body weight of Draa goat. Small Ruminant Research 69: 149-156.

Chungyan Z, Ligu Y, Zhang S (2008). Variance components and genetic parameters for weight and size at birth in Boer goat. Livestock Production Science 115: 73-79.

Chungyan Z, Chen SL, Li X, Xu DQ, Ying Z, Yang LG (2009). Genetic and phenotypic parameter estimates for reproduction traits in the Boer dam. Livestock Science 125: 60-65.

Curia RA, De HN, Oliveirab A, Silveirab C, Lopesa CR (2005). Association between IGF-I, IGF-IR and GHRH gene polymorphisms and growth and carcass traits in beef cattle. Livestock Production Science 94: 159-167.

Emami Meybodi MA, Eftekhar Shahrodi F, Nikkhah A (1993). The estimate genetic parameters for some economic traits in Raieni Cashmere goats. Ms.c. Thesis. Ferdosi University, Mashhad, Iran (In Farsi).

Ghazikhani A, Moradi Shahr_e Babak H, Sadeghi M, Faraji R (2012). The genetic analysis exon4 of IGF-1 gene and its relationship with carcass, fatty acid and blood metabolic traits in Zel sheep. Journal of Animal Science and Research 12: 51-61. (In Farsi). 

Hermiz HN, Alkass JE, Hobi AA, Asofi MK (2009). Genetic and phenotypic parameters of body weights in Iraqi local goat and their crosses with Damascus. Journal of University of Duhok 12: 189-194.

Glimur AR, Gogel Bj, Cullis BR, Thompson R (2009). ASReml user guide. NSW Department of Industry and Investment, VSN International Ltd, Hemel Hempstead, HP1 1ES, UK.

Haji Hossenloo A, Pirani N, Hashemi A, Philkosh Moghadam F, Jaefari Sh (2011). The study of polymorphism IGF-1 gene using PCR-SSCP in Makoee sheep. The First Conference on the Role of Biotechnology in Animal Science, Isfahan University of Technology. 

Hassani MN, Asadi Fozi M, Esmailizadeh AK, Mohammadabadi MR (2010). A genetic analysis of growth traits in Raieni Cashmere goat using multivariate animal model. Iranian Journal of Animal Science 41: 323-329 (In Farsi).

Jack ES, Felix CS (1996). Genetic marker map construction and their application in plant breeding. Hort Science 31: 729-741.

Javanmard A, Mohammadabadi MR, Zarrigabayi GE, Gharahedaghi AA, Nassiry MR, Javadmansh A, Asadzadeh N (2008). Polymorphism within the intron region of the bovine leptin gene in Iranian Sarabi cattle (Iranian Bos taurus). Russian Journal of Genetics 44: 495-497.

Karimi Kurdistania Z, Rostamzadeha J, Rashidia A, Davisb ME (2012). Evaluation of insulin-like growth factor-I gene polymorphism on growth traits and yearling fleece weight in goats. Small Ruminant Research 111: 10-15.

Kim ES, Shi X, Cobanoglu O, Weigel K, Berger PJ, Kirkpatrick BW (2009). Refined mapping of twinning-rate quantitative trait loci on Bovine chromosome 5 and analysis of insulin like growth factor-1 as a positional candidate gene. Journal of Animal Science 87: 835-843.

Li MH, Adamowicz T, Switonski M, Ammosov I, Ivanova Z,  Kiselyova T, Popov R, Kantanen J (2006). Analysis of population differentiation in North Eurasian cattle (Bos taurus) using single nucleotide polymorphisms in three genes associated with production traits. Animal Genetics 37: 390-392.

Liu WJ, Fang GX, Fang Y, Tian KC, Huang XX, Yao XK, Wang M, Yu H, Huang YZ, Xin JJ, Xin YP, Yu SG, Chen H (2010). The polymorphism of IGF-1 gene on two goat breeds in China. Journal of Animal and Veterinary Advances 9: 790-794.

Ma RN, Deng CJ, Zhang XM (2010). S-novel SNP of alpha-lactalbumin gene in Chinese dairy goats. Molecular Biology 44: 608-612.

Moghadaszadeh M, Mohammadabadi MR, Esmailizadeh Koshkoieh A (2015). Association of exon 2 of BMP15 gene with the litter size in the Raieni Cashmere goat. Genetics in the 3rd Millennium 13: 4062-4067 (In Farsi).

Mohammadabadi MR (2012). Relationships of IGFBP-3 gene polymorphism with cashmere traits in Raieni Cashmere goat. Modern Genetics Journal 7: 115-120 (In Farsi).

Mohammadabadi MR, Askari N, Baghizadeh A, Esmailizadeh A (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.

Mohammadi A, Nassiry MR, Mosafer J, Mohammadabadi MR, Sulimova GE (2009). Distribution of BoLA-DRB3 allelic frequencies and identification of a new allele in the Iranian cattle breed Sistani (Bos indicus). Russian Journal of Genetics 45: 198-202.

Mohebi Nejad E, Asadi Fozi M (2012). Effect of environmental factors on fleece weight of Raeini Cashmere goats. Journal of Animal Science Research 22: 1-9 (In Farsi).

Mousavizadeh A, Mohammadabadi MR, Torabi A, Nassiry MR, Ghiasi H, Esmailizadeh AK (2009). Genetic polymorphism at the growth hormone locus in Iranian Talli goats by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Iranian Journal of Biotechnology 7: 51-53 (In Farsi).

Patton J, Kenny DA, McNamara S, Mee JF, O'Mara FP, Diskin MG, Murphy JJ (2007). Relationships among milk production, energy balance, plasma analyzes and reproduction in Holstein-Friesian cows. Journal of Dairy Science 90: 649-658.

Ping-qing W, Ying T, Bao-Yun Z, Ming-Xing C, La-Mei1 D, Qi1 F, Chong-Xu (2011). DNA polymorphisms of 5´-flanking region of Insulin-like growth factor 1 gene and their association with reproduction traits in goats. Agricultural Sciences in China 10: 1609-1617.

Restall BJ, Pattie WA (1989). The inheritance of cashmere in Australian goats. 2. Genetic parameters and breeding value. Livestock Production Science 21: 251-261.

Roite D, Bondy C, Yakar S, Liu JL, Butler A (2001). The somatomedin hypothesis endocrine. Journal of Clinical Endocrinology 22: 53-74.

Shamsalddini S, Mohammadabadi MR, Esmailizadeh AK (2016). Polymorphism of the prolactin gene and its effect on fiber traits in goat. Russian Journal of Genetics. 52: 405-408.

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.

Thomas MG, Enns RM, Shirley KL, Garcia MD, Garrett AJ, Silver GA (2007). Associations of DNA polymorphisms in growth hormone and its transcriptional regulators with growth and carcass traits in two populations of Brangus bulls. Genetics and Molecular Research 6: 222-237.

Tohidi nezhad F, Mohammadabadi MR, Esmailizadeh AK, Najmi Noori A (2015). Comparison of different levels of Rheb gene expression in different tissues of Raieni Cashmir goat. Journal of Agricultural Biotechnology 6: 35-50. (In Farsi).

Velazquez MA, Newman M, Christie MF, Cripps P, Crowe MA, Smith RF, Dobson H (2005). The usefulness of a single measurement of insulin-like growth factor-1 as a predictor of embryo yield and pregnancy rates in a Bovine MOET program. Theriogenology 64: 1977-1994.

Velazquez MA, Spicer LJ, Wathes DC (2008). The role of endocrine insulin-like growth factor-I (IGF-I) in female Bovine reproduction. Domestic Animal Endocrinology 35: 325-342.

Walkden-brown S, Sunduimijid B, Olayemi B, van der Werf J, Ruvinsky A (2008). Breeding fiber goat for resistance to worm infections. Rural Industries Research and Development Corporation, University of New England, Armidale, Australia, pp. 70-87.

Yeh FC, Yang R, Boyle T (1999). PopGene. Version 1.31. Microsoft window-based free software for population genetic analysis. University of Alberta, Canada.

Yilmaz A, Davis ME, Hines H, Chung H (2005). Detection of two nucleotide substitutions and putative promoters in the 5´ flanking region of the ovine IGF-I gene. Journal of Applied Genetics 46: 307-309.

 

Zamani P, Akhondi M, Mohammadabadi MR (2015). Associations of Inter-Simple Sequence Repeat loci with predicted breeding values of body weight in Sheep. Small Ruminant Research 132:123–127

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.

 

 

 

 

 

 


The Analysis of IGF-I Gene Polymorphism and its Association with Breeding Values of Production and Reproduction Traits in Raieni Cashmere Goats

 

Molaei M.H. 1, Asadi Fozi M.*2, Mohammadabadi M.R.3, Montazeri M.4,5

 

1M. Sc. Student, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, 76169-133, Iran.

2Associate professor, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, 76169-133, Iran.

3Professor, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, 76169-133, Iran.

4Ph.D Student, Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, 76169-133, Iran.

5Young Researchers Society, Shahid Bahonar University of Kerman, Kerman, Iran.

 

Abstract

The Insulin-like growth factor I (IGF-I) has a similar structure with Insulin hormone. This hormone which is transcribed by (IGF-I) gene has an important functional in growth and development of different body tissues. The aim of this study was to investigate the polymorphism of this gene and its relationship with the estimated breeding values of the birth type, birth weight and fleece weight in Raieni Cashmere goats. Breeding values were estimated using records on 13020 Raieni cashmere goats originated from 336 sire and 3617 dams. Univariate animal model was applied for the genetic analysis of the investigated traits using ASReml. To study IGF-1 gene polymorphism, 94 Raieni cashmere goats were selected based on their estimated breeding values for these traits. PCR was used for amplification of 363 bp fragment of exon 4 of IGF-I gene. Thereafter the animal’s genotype was determined using PCR-RFLP. Three genotypes including AA, AB and BB were observed. The frequency of these genotypes was 0.22, 0.65 and 0.13, respectively. No significant associations were found between these genotypes and the investigated production and reproduction traits. The results of this study show that, there is important genetic diversity for the IGF-I in the investigated population but to study the association between the polymorphism of this gene and the studied traits, more data should be used as currently few animals were genotyped.  

Keywords: Raeini Cashmere goat, IGF-I gene, Breeding value, PCR-RFLP.

 



* نویسنده مسئول: مسعود اسدی فوزی                             تلفن: 09132410774                                   Email: masadi@uk.ac.ir       

[1] Log Likelihood Ratio Test

* Corresponding Author: Asadi Fozi M.                          Tel: 09132410774              Email: masadi@uk.ac.ir

Abdolmohammadi A (2008). Study of genetic variation for four candidate genes using PCR-RFLP and HRM and their association with reproduction and production traits in Holstein cows of Iran. Thesis of Ph.D. University of Tehran, Iran  (In Farsi).
Alinaghizadeh H, Mohammadabadi MR, Zakizadeh S (2010). Exon 2 of BMP15 gene polymorphismin Jabal Barez red goat. Journal of Agricultural Biotechnology 2: 69-80 (In Farsi).
Allain D, Roguet J (2003). Genetic and nongenetic factors influencing mohair production traits within the national selection scheme of Angora goats in France. Livestock Production Science 82: 129-137.
Asadi Fozi M (2005). Appropriate models for genetic analysis and efficient selection of Merino sheep. Ph.D. Thesis. University of New England, Armadale, Australia.
Askari N, Baghizadeh A, Mohammadabadi MR (2008). Analysis of the genetic structure of Iranian indigenous Raieni Cashmere goat populations using microsatellite markers.Biotechnology 2: 1-4 (In Farsi).
Askari N, Baghizadeh A, Mohammadabadi MR (2010). Study of genetic diversity in four populations of Raeini Cashmere goat using ISSR markers. Modern Genetics Journal 5: 49-56 (In Farsi).
Askari N, Mohammadabadi MR, Baghizadeh A (2011). ISSR markers for assessing DNA polymorphism and genetic characterization of cattle, goat and sheep populations. Iranian Journal of Biotechnology 9: 222-229 (In Farsi).
Askari N, Mohammadabadi MR, Beygi Nassiry MT, Baghizadeh A, Fayazi J (2009). Study of genetic diversity of Raeini Cashmere goat based on microsatellite markers. Journal of Agricultural Science 18: 155-161 (In Farsi).
Baghizadeh A, Bahaaddini M, Mohamadabadi MR, Askari N (2009). Allelic variations in exon 2 of caprine MHC class II DRB3 gene in Raeini Cashmere goat. American-Eurasian Journal of Agricultural and Environmental Science 6: 454-459.
Bechmann JS, Soller M (1987). Molecular markers in the genetic improvement of farm animals. Biotechnology 5: 573-576.
Bharathidhasan A, Narayanan R, Gopu P, Subramanian A, Prabakaran R, Rajendran R (2009). Effect of nongenetic factors on birth weight, weaning weight and pre weaning gain of Barbari goat.  Journal of Veterinary & Animal Science 5: 99-103.
Boujenane I, Hazzab AE (2008). Genetic parameter for direct and mater effect on body weight of Draa goat. Small Ruminant Research 69: 149-156.
Chungyan Z, Ligu Y, Zhang S (2008). Variance components and genetic parameters for weight and size at birth in Boer goat. Livestock Production Science 115: 73-79.
Chungyan Z, Chen SL, Li X, Xu DQ, Ying Z, Yang LG (2009). Genetic and phenotypic parameter estimates for reproduction traits in the Boer dam. Livestock Science 125: 60-65.
Curia RA, De HN, Oliveirab A, Silveirab C, Lopesa CR (2005). Association between IGF-I, IGF-IR and GHRH gene polymorphisms and growth and carcass traits in beef cattle. Livestock Production Science 94: 159-167.
Emami Meybodi MA, Eftekhar Shahrodi F, Nikkhah A (1993). The estimate genetic parameters for some economic traits in Raieni Cashmere goats. Ms.c. Thesis. Ferdosi University, Mashhad, Iran (In Farsi).
Ghazikhani A, Moradi Shahr_e Babak H, Sadeghi M, Faraji R (2012). The genetic analysis exon4 of IGF-1 gene and its relationship with carcass, fatty acid and blood metabolic traits in Zel sheep. Journal of Animal Science and Research 12: 51-61. (In Farsi). 
Hermiz HN, Alkass JE, Hobi AA, Asofi MK (2009). Genetic and phenotypic parameters of body weights in Iraqi local goat and their crosses with Damascus. Journal of University of Duhok 12: 189-194.
Glimur AR, Gogel Bj, Cullis BR, Thompson R (2009). ASReml user guide. NSW Department of Industry and Investment, VSN International Ltd, Hemel Hempstead, HP1 1ES, UK.
Haji Hossenloo A, Pirani N, Hashemi A, Philkosh Moghadam F, Jaefari Sh (2011). The study of polymorphism IGF-1 gene using PCR-SSCP in Makoee sheep. The First Conference on the Role of Biotechnology in Animal Science, Isfahan University of Technology. 
Hassani MN, Asadi Fozi M, Esmailizadeh AK, Mohammadabadi MR (2010). A genetic analysis of growth traits in Raieni Cashmere goat using multivariate animal model. Iranian Journal of Animal Science 41: 323-329 (In Farsi).
Jack ES, Felix CS (1996). Genetic marker map construction and their application in plant breeding. Hort Science 31: 729-741.
Javanmard A, Mohammadabadi MR, Zarrigabayi GE, Gharahedaghi AA, Nassiry MR, Javadmansh A, Asadzadeh N (2008). Polymorphism within the intron region of the bovine leptin gene in Iranian Sarabi cattle (Iranian Bos taurus). Russian Journal of Genetics 44: 495-497.
Karimi Kurdistania Z, Rostamzadeha J, Rashidia A, Davisb ME (2012). Evaluation of insulin-like growth factor-I gene polymorphism on growth traits and yearling fleece weight in goats. Small Ruminant Research 111: 10-15.
Kim ES, Shi X, Cobanoglu O, Weigel K, Berger PJ, Kirkpatrick BW (2009). Refined mapping of twinning-rate quantitative trait loci on Bovine chromosome 5 and analysis of insulin like growth factor-1 as a positional candidate gene. Journal of Animal Science 87: 835-843.
Li MH, Adamowicz T, Switonski M, Ammosov I, Ivanova Z,  Kiselyova T, Popov R, Kantanen J (2006). Analysis of population differentiation in North Eurasian cattle (Bos taurus) using single nucleotide polymorphisms in three genes associated with production traits. Animal Genetics 37: 390-392.
Liu WJ, Fang GX, Fang Y, Tian KC, Huang XX, Yao XK, Wang M, Yu H, Huang YZ, Xin JJ, Xin YP, Yu SG, Chen H (2010). The polymorphism of IGF-1 gene on two goat breeds in China. Journal of Animal and Veterinary Advances 9: 790-794.
Ma RN, Deng CJ, Zhang XM (2010). S-novel SNP of alpha-lactalbumin gene in Chinese dairy goats. Molecular Biology 44: 608-612.
Moghadaszadeh M, Mohammadabadi MR, Esmailizadeh Koshkoieh A (2015). Association of exon 2 of BMP15 gene with the litter size in the Raieni Cashmere goat. Genetics in the 3rd Millennium 13: 4062-4067 (In Farsi).
Mohammadabadi MR (2012). Relationships of IGFBP-3 gene polymorphism with cashmere traits in Raieni Cashmere goat. Modern Genetics Journal 7: 115-120 (In Farsi).
Mohammadabadi MR, Askari N, Baghizadeh A, Esmailizadeh A (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.
Mohammadi A, Nassiry MR, Mosafer J, Mohammadabadi MR, Sulimova GE (2009). Distribution of BoLA-DRB3 allelic frequencies and identification of a new allele in the Iranian cattle breed Sistani (Bos indicus). Russian Journal of Genetics 45: 198-202.
Mohebi Nejad E, Asadi Fozi M (2012). Effect of environmental factors on fleece weight of Raeini Cashmere goats. Journal of Animal Science Research 22: 1-9 (In Farsi).
Mousavizadeh A, Mohammadabadi MR, Torabi A, Nassiry MR, Ghiasi H, Esmailizadeh AK (2009). Genetic polymorphism at the growth hormone locus in Iranian Talli goats by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Iranian Journal of Biotechnology 7: 51-53 (In Farsi).
Patton J, Kenny DA, McNamara S, Mee JF, O'Mara FP, Diskin MG, Murphy JJ (2007). Relationships among milk production, energy balance, plasma analyzes and reproduction in Holstein-Friesian cows. Journal of Dairy Science 90: 649-658.
Ping-qing W, Ying T, Bao-Yun Z, Ming-Xing C, La-Mei1 D, Qi1 F, Chong-Xu (2011). DNA polymorphisms of 5´-flanking region of Insulin-like growth factor 1 gene and their association with reproduction traits in goats. Agricultural Sciences in China 10: 1609-1617.
Restall BJ, Pattie WA (1989). The inheritance of cashmere in Australian goats. 2. Genetic parameters and breeding value. Livestock Production Science 21: 251-261.
Roite D, Bondy C, Yakar S, Liu JL, Butler A (2001). The somatomedin hypothesis endocrine. Journal of Clinical Endocrinology 22: 53-74.
Shamsalddini S, Mohammadabadi MR, Esmailizadeh AK (2016). Polymorphism of the prolactin gene and its effect on fiber traits in goat. Russian Journal of Genetics. 52: 405-408.
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.
Thomas MG, Enns RM, Shirley KL, Garcia MD, Garrett AJ, Silver GA (2007). Associations of DNA polymorphisms in growth hormone and its transcriptional regulators with growth and carcass traits in two populations of Brangus bulls. Genetics and Molecular Research 6: 222-237.
Tohidi nezhad F, Mohammadabadi MR, Esmailizadeh AK, Najmi Noori A (2015). Comparison of different levels of Rheb gene expression in different tissues of Raieni Cashmir goat. Journal of Agricultural Biotechnology 6: 35-50. (In Farsi).
Velazquez MA, Newman M, Christie MF, Cripps P, Crowe MA, Smith RF, Dobson H (2005). The usefulness of a single measurement of insulin-like growth factor-1 as a predictor of embryo yield and pregnancy rates in a Bovine MOET program. Theriogenology 64: 1977-1994.
Velazquez MA, Spicer LJ, Wathes DC (2008). The role of endocrine insulin-like growth factor-I (IGF-I) in female Bovine reproduction. Domestic Animal Endocrinology 35: 325-342.
Walkden-brown S, Sunduimijid B, Olayemi B, van der Werf J, Ruvinsky A (2008). Breeding fiber goat for resistance to worm infections. Rural Industries Research and Development Corporation, University of New England, Armidale, Australia, pp. 70-87.
Yeh FC, Yang R, Boyle T (1999). PopGene. Version 1.31. Microsoft window-based free software for population genetic analysis. University of Alberta, Canada.
Yilmaz A, Davis ME, Hines H, Chung H (2005). Detection of two nucleotide substitutions and putative promoters in the 5´ flanking region of the ovine IGF-I gene. Journal of Applied Genetics 46: 307-309.
 
Zamani P, Akhondi M, Mohammadabadi MR (2015). Associations of Inter-Simple Sequence Repeat loci with predicted breeding values of body weight in Sheep. Small Ruminant Research 132:123–127
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.