مطالعه چندشکلی اگزون 4 ژن -IIGF و ارتباط آن با صفات رشد و تولید شیر در بز نژاد مهابادی به روش PCR-SSCP

مطالعه چندشکلی اگزون 4 ژن -IIGF و ارتباط آن با صفات رشد و تولید شیر در بز نژاد مهابادی به روش PCR-SSCP

لیلا قره داغی*¹، حسین مرادی شهربابک²، مصطفی صادقی²، مهدی گنج خانلو²

^1و2 به ترتیب دانشجوی کارشناسی ارشد ژنتیک و اصلاح دام و استادیار گروه علوم دامی، پردیس کشاورزی و منابع طبیعی کرج دانشگاه تهران،کرج، ایران.

تاریخ دریافت: 27/08/1391، تاریخ پذیرش: 26/04/1392

چکیده

پژوهش حاضر با هدف مطالعه ارتباط چند شکلی ژن IGF-I با صفات رشد و تولید شیر در بزهای نژاد مهابادی به روش PCR-SSCP انجام شد. در این مطالعه از 89 رأس بز و 51 رأس بزغاله نژاد مهابادی موجود در مزرعه پردیس کشاورزی و منابع طبیعی دانشگاه تهران، واقع در کرج با استفاده از لوله­های خلاء حاوی EDTA از سیاهرگ وداج خونگیری شد. استخراج DNA از خون با روش نمکی بهینه یافته انجام شد و واکنش زنجیره­ای پلی­مراز (PCR) برای تکثیر قطعه 365 جفت بازی اگزون 4 ژن IGF-I صورت گرفت. برای تعیین چندشکلی نمونه­ها الکتروفورز محصولات PCR پس از تک رشته­ای شدن قطعات (روش SSCP) روی ژل پلی­اکریل­آمید و رنگ­آمیزی ژل به روش نیترات نقره انجام شد. نتایج بیانگر وجود دو الگوی باندی متفاوت بود که فراوانی الگوهای ژنوتیپی 1 و 2 به ترتیب 7/65 % و 3/34% بود. آنالیز آماری نشان داد که صفت تولید شیر ارتباط معنی­داری با الگوهای باندی دارد P<0.05))، ولی اثر ژنوتیپ­های IGF-I بر درصد چربی، پروتئین و شمار سلول­های بدنی معنی­دار نبود P<0.05))، همچنین آنالیز داده­های مربوط به صفات رشد، نشان داد که ژنوتیپ­های مختلف IGF-I، اثر معنی­داری بر وزن تولد، وزن از شیرگیری، افزایش وزن روزانه، مقدار خوراک مصرفی80 روز و وزن نهایی بزغاله­ها نداشت.    

کلمات کلیدی: ژن IGF-I، چندشکلی، PCR-SSCP ، بز نژاد مهابادی.

مقدمه

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

هورمون­های رشد شبه­انسولینی (IGF) مواد پروتئینی با وزن مولکولی کم هستند که در پاسخ به هورمون رشد آزاد شده وهورمون رشد اثرات خود را از طریق این عوامل القا می­کند. به فاکتور­های رشد شبه انسولینی سوماتومدین نیز گفته می­شود و به دلیل اینکه این مواد واسطه عمل هورمون رشد هستند به نام هورمون رشد سوماتوتروپ نیز شناخته می­شوند. این فاکتورها اثراتی شبیه انسولین بر سلول­ها دارند ولی برخلاف انسولین در مایعات بیولوژیکی تقریبا 98% آن به یکی از 6 پروتئین باند شونده (IGFBP) وصل می شود (Zhou yinli, 2007). سیستم سیگنال دهنده­ی [1]IGFs  شامل IGF-I، IGF-II، گیرنده­ی IIGF-، گیرنده­ی IGF-II و شش پروتئین باند شونده (IGFBP-6~IIGFBP) نقش مهمی در توسعه، رشد و تولیدمثل پستانداران ایفا می­کنند (Duan et al., 2005). IGF-I پلی پپتید تک زنجیره­ای با سه پل دی­سولفید دارای وزن مولکولی 7649 دالتون، و شامل 70 اسیدآمینه می باشد که هورمون رشد مترشحه از غده هیپوفیز عامل اصلی تحریک آزادسازی آن است (Renderknecht & Humbel, 1978). هپاتوسیت­ها یا سلول­های کبدی محل اولیه و اصلی ترشح IGF-I به عنوان یک هورمون اندوکرین می­باشد(Combes et al., 1997)، ولی در بعضی بافت­های محیطی نظیر استخوان نیز ممکن است آزادسازی سوماتومدین­ها به صورت اوتوکرین/ پاراکرین تحت تاثیر عوامل یا فاکتورهای موضعی صورت گیرد (Yee, 1994). نقش این هورمون در رشد بافت های مختلف بدن (سلول های ماهیچه، غضروف و استخوان)، تحریک سنتز پروتئین، افزایش متابولیسم قندها و چربی ها در بدن، محرک تقسیم میتوز، تحریک ازدیاد سلول های اپیتلیال (مخاطی) و استرومال غده پستانی، سنتز و بیان ژن کازئین و حمل گلوکز شناخته شده است (Zhou Yinli, 2007). علاوه بر این­ IGF-I سنتز شیر و تکامل غده پستانی را در طول دوره­ی شیردهی تسهیل می­کند. مطالعات نشان دادند که IGF-I یک محرک قوی برای ماندگاری سلول­ها، جلوگیری از مرگ و میر برنامه ریزی شده­ی سلول[2]­ها و بیان ژن­های مرتبط با فرآیند مرگ و میر برنامه ریزی شده­ی سلول­های غده پستانی شناخته شده است، همچنین بیان ژن­هایی که در انتقال و متابولیسم آمینواسید­ها و همچنین بیوسنتز و پایداری پروتئین­ها درگیر هستند را تنظیم می­کند.(Yinli, 2007 Zhou). ژن IGF-I دارای 6 اگزون و 5 اینترون است و در5/73 سانتی مورگانی کروموزوم 5 نقشه یابی شده است 2008)  Wibowo et al.,).

در مطالعه چند شکلی ناحیه 5´-Flanking ژن IGF-I با استفاده از روش PCR-SSCP در ده نژاد بز چینی، نتایج حاصل از توالی یابی این ناحیه حاکی از حذف دو نوکلئوتیدCA)) و همچنین وجود دو SNP به ترتیب در نوکلئوتید­های 1637(A به G)،1640(T به C ) بودکه فقط در یکی از نژادها ارتباط معنی داری با چندقلوزایی مشاهده شد (P<0.05)، (.( Ping-Qiong et al., 2011 در مطالعه دیگری روی همان ناحیه ژن IGF-I در گوسفند نژاد بلوچی 3 الگوی باندی متفاوت مشاهده شد و تفاوت معنی­داری در متوسط افزایش وزن روزانه از تولد تا شیرگیری بین افراد با ژنویپ AB وBB مشخص شد P<0.05)) و افراد با ژنوتیپ­AB  در مقایسه با افراد ژنوتیپ­های دیگر بالاترین متوسط افزایش وزن روزانه از تولد تا شیرگیری را نشان دادند P<0.05)) (Tahmoorespur et al., 2009). در بررسی اگزون 4 ژن IGF-I در سه نژاد بومی بز در چین با استفاده از روشPCR-SSCP  و توالی­یابی DNA، دو جهش در نوکلئوتیدهای 1617 (G>A) و1620(C>T) گزارش شد که از نوع جهش­های خاموش بودند. برای اولین بار این SNP ها در بز گزارش شد که به طور معنی­داری با صفت تولید کشمیر مرتبط بود، افراد دارای 4 نوع ژنوتیپ بودند که ژنوتیپ AA به طور معنی­داری تولید کشمیر پایین­تری نسبت به ژنوتیپ AB داشتند P<0.05)). وزن بدن ژنوتیپ AC به طور معنی­داری بالاتر از ژنوتیپ BB بودp<0.05) ) (Quing et al., 2011). در مطالعه چندشکلی اگزون 4 ژن IGF-I با روش PCR-RFLP در دو نژاد بز چینی، ارتباط معنی­داری بین تنوع ژنتیکی در این ناحیه با تولید کشمیر، قطرالیاف و وزن بدن بز مشاهده نشد (Wu-Jun et al., 2010).

هدف از انجام پژوهش حاضر، شناسایی چندشکلی­های موجود در اگزون 4 ژن IGF-I و برآورد میزان فراوانی ­الگوهای ژنوتیپی مختلف این جایگاه و بررسی ارتباط بین چند­شکلی­های این جایگاه ژنی با صفت تولید شیر در بزهای نژاد مهابادی و صفات رشد در بزغاله­ها  با استفاده از تکنیک PCR-SSCP بود.

مواد و روش­ها

در این مطالعه از 89 رأس بز و 51 رأس بزغاله نژاد مهابادی موجود در مزرعه پردیس کشاورزی و منابع طبیعی دانشگاه تهران واقع در کرج، با استفاده از لوله­های خلاء حاوی EDTA از سیاهرگ وداجی خونگیری انجام شد. استخراج DNA با روش نمکی تغییر یافته صورت گرفت. برای تکثیر قطعه 326 جفت بازی ژن IGF-I از آغازگرهای رفت و برگشت که توسط نرم افزارVector NTI  طراحی شد، استفاده گردید که این آغازگرها توسط شرکت متابیون سنتز شد. توالی آغازگرها به شرح زیر بود: توالی آغازگر رفت: 5`-ATT TGA ACA GAC AAG CCC AC-3` و توالی آغازگر برگشت5`- TTT GAC ACT ATG AGC CAG AAG-3` جهت بهینه سازی واکنش PCR برنامه­های حرارتی مختلفی استفاده شد. ولی در نهایت برنامه­ حرارتی زیر با 35 چرخه، ایده آل ترین شرایط برای تکثیر ژن IGF-I تشخیص داده شد. واسرشت سازی اولیه DNA به مدت 5 دقیقه در دمای 95 درجه سانتی گراد، واسرشت سازی DNA به مدت 30 ثانیه در دمای 95  درجه سانتی گراد، اتصال آغازگر به DNA به مدت 30 ثانیه در دمای 4/61 درجه سانتی گراد، بسط DNA به مدت 30 ثانیه در دمای 72 درجه سانتی گراد و بسط نهایی به مدت 10 دقیقه در دمای 72 درجه سانتی گراد انجام شد. واکنش PCR بوسیله مواد مخصوص PCR شرکت GENNET به شرح زیر با غلظت نهایی واکنش 25 میکرولیتر انجام شد: کلرید منیزیم (MgCl2) 5/1 میکرولیتر ( mM/µL5/1)، بافر (PCR buffer) 5/2 میکرولیتر (1X)، بازهای آزاد دزوکسی نوکلئوتید تری فسفات (dNTPs) 2 میکرولیتر (mM/µL 8/0)، آغاز گر رفت 1 میکرولیتر (pm/­µL 5/0)، آغازگر برگشت 1 میکرولیتر (pm/­µL 5/0)، آنزیم DNA پلی مراز Taq 3/0 میکرولیتر (unit/­µL 25/0)، DNA 2 میکرولیتر، آب دیونیزه (dH2O) 7/14 میکرولیتر.

الکتروفورز محصولات PCR روی ژل آگارز 2 درصد، با ولتاژ 92 ولت، به مدت نیم ساعت صورت گرفت. برای رنگ آمیزی ژل آگارز ازDNA Safe Stain  استفاده شد. پس از انجام الکتروفورز از دستگاه اشعه ماورای بنفش برای بررسی قطعه تکثیر شده استفاده شد. به منظور انجام واکنش SSCP، محصول PCR تکثیر شده، با استفاده از بافر SSCP ( فرمامید[3] 99% (1960 میکرولیتر)، 0.5m اتیلن دی­آمید تترا استیک اسید[4] (40 میکرولیتر)، بروموفنول­بلو[5] (0.025%)، زینول­سیانید[6]  (0.25)) به نسبت 4 به 16 به مدت 10 دقیقه در دمای 96 درجه تک رشته ای شدند که بلافاصله پس از این مرحله در داخل یخ قرار گرفتند. الکتروفورز محصولات تک رشته­ای، روی ژل پلی­اکریل­آمید 10 درصد ، با ولتاژ300 به مدت 24 ساعت انجام و رنگ آمیزی ژل با استفاده از نیترات نقره طی سه مرحله تثبیت، لکه گذاری و ظهور صورت گرفت.                                                        

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

در این تحقیق پس از تعیین ژنوتیپ دام­ها برای جایگاه مورد بررسی این اطلاعات به همراه داده های مربوط به تولید شیر، درصد پروتئین شیر، درصد چربی شیر و شمارش سلول­های بدنی وارد نرم­افزار Excel گردید و بعد از نرمال­کردن داده­ها (روش لگاریتم­گیری) از نرم­افزار9.1  SAS و رویه Mixed (به دلیل رکوردگیری تکراری برای هر حیوان) برای آنالیز داده­ها استفاده شد. در این آزمایش اثر ژنوتیپ­های ژن IGF-I به عنوان عامل ثابت تاثیرگذار بر صفات مورد مطالعه در مدل قرار داده شد به این صورت که الگوهای ژنوتیپی مشاهده شده­ی این ژن (2 الگو)، اثر سن حیوان هنگام زایش، سال زایش، و ماه رکوردگیری به عنوان عوامل ثابت و اثر صفت تعداد روزهای شیردهی به عنوان عامل کواریت، همچنین اثر حیوان به عنوان عامل تصادفی در مدل قرار داده شد. معادله مدل به شکل زیر بود.

Y_{ijklmn  =}  + A_i  + G_j  + Y_k + M_l + b (D_ijklmn -  (+Animal_n(G_i)+ e_ijklmn

Y_ijklmn: هریک از مشاهدات مربوط به مقدار شیر، درصد چربی، درصد پروتئین، شمار سلول­های بدنی، :m میانگین، G_j: اثر j امین ژنوتیپ (IGF-1)، A_i: اثر i امین سن حیوان هنگام زایش، Y_k : اثر k امین سال زایش، M_l: اثر 1 امین ماه رکوردگیری، Animal_n(G_i) _n: اثر تصادفی حیوان، b: ضریب تابعیت صفت مورد نظر از تعداد روزهای شیردهی، D_ijklmn: تعداد روزهای شیردهی، : میانگین تعداد روزهای شیردهی، e_ijklmn: اثر عوامل باقی­مانده.                                        

برای تجزیه واریانس اثر چند­شکلی ژن IGF-I بر صفات رشد، داده­های مربوط به وزن تولد، وزن از شیرگیری، افزایش وزن روزانه،  میانگین مقدار خوراک مصرفی 80 روز و تاریخ تولد مربوط به 51 بزغاله ی مورد مطالعه، پس از تعیین ژنوتیپ دام­ها و ثبت در نرم افزار Excel، با استفاده از نرم افزار SAS با رویه GLM انجام شد. در این مطالعه اثر ژنوتیپ، سال تولد، ماه تولد و سن دام هنگام زایش به عنوان عوامل ثابت، همچنین وزن مادر هنگام زایش و وزن تولد بزغاله به عنوان عامل کواریت در معادله مدل قرار داده شد.

دو معادله مدل برای تجزیه واریانس مربوط به صفات رشد به کار گرفته شد که به صورت زیر می­باشند.

Y1_ijklm=μ+G_i+Y_j+M_k+Al+b₁ (W_ijklm-) + e_ijklm

  Y2ijklmn=μ+Gi+Yj+Mk+Al+b1 (W_ijklm-) + b2 (Wijklmn -_n ) + eijklmn

Y1_ijklmn= هر یک از مشاهدات مربوط به وزن تولد و میانگین خوراک مصرفی 80 روز، هریک از مشاهدات مربوط به وزن از شیرگیری، =Y_2ijklm وزن نهایی و میانگین خوراک مصرفی 80 روز، m: میانگین صفت در جامعه، Gi=اثر i امین ژنوتیپ (IGF-1)، Yj= اثر j امین سال تولد، M_k= اثر  kامین ماه تولد، A_l= اثرl امین سن مادر هنگام زایش، W_ijklm=وزن ijklm امین مادر هنگام زایش، = میانگین وزن مادر هنگام زایش،  b₁=ضریب تابعیت Y روی W (وزن دام هنگام زایش)، W_ijklmn= اثر ijklmn امین وزن تولد، b₂=ضریب تابعیت Y روی W (وزن تولد بزغاله­ها)، = میانگین وزن تولد بزغاله ها، e_ijklmn= اثر عوامل باقیمانده.

نتایج و بحث

با استفاده از روش اسپکتوفتومتری و ژل آگارز 1 درصد، کمیّت و کیفیت DNA استخراج شده مورد تایید قرار گرفت (شکل 1 و 2).

شکل1-کیفیت نمونه های DNA بر روی ژل آگارز.

 Figure 1- The quality of DNA samples on agarose gel.

شکل 2- تعیین کمیت و کیفیت DNA با دستگاه نانودراپ.

              Figure 2- Determination of DNA quality by Nonodrop. 

قطعه 326 جفت بازی ژن IGF-I، طی واکنش PCR  تکثیر (شکل 3) و در مرحله بعد، از روش SSCP برای شناسایی تنوع در قطعه تکثیر شده استفاده شد. نتایج حاصل از SSCP و رنگ­آمیزی ژل پلی­اکریل آمید بیانگر دو الگوی باندی متفاوت بود که در شکل4 مشاهده می­شود. و فراوانی الگوهای 1 و 2 برابر با 71/65% و 29/34% بود. همانطور که Wu-Jun et al., (2010) وQuing et al.,  (2011) در مطالعه اگزون 4 ژن IGF-I در نژادهای مختلف بز چینی، چند­شکلی در این ناحیه را گزارش کرده بودند، طبق انتظار مطالعه ناحیه مذکور در بز نژاد مهابادی نیز دارای چند­شکلی می­باشد. جدول 1 آماره­های توصیفی مربوط به رکورد تولید شیر، درصد پروتئین، درصد چربی و شمار سلول­های بدنی را نشان می­دهد. بالا بودن ضریب تنوع در صفات مورد مطالعه می­تواند به دلیل رکوردگیری در دوره­های مختلف شیردهی هر کدام از حیوانات باشد.

شکل3- تکثیر اختصاصی قطعه  bp326  ژن IGF-1 بر روی ژل آگارز.

Figure 3- PCR product analyzed by electrophoresis on 2% agarose gel (326bp).

شکل 4- الگوهای حاصل از الکتروفورز محصولاتSSCP  اگزون 4  ژن IGF-1 در بزهای مهابادی.

Figure 4- SSCP pattern of IGF-1 gene in Mahabadi goats.

جدول 1- آماره­های توصیفی صفات تولید شیر، درصد چربی، درصد پروتئین و شمار ش سلول­های بدنی شیر.

Table 1- Descriptive statistical of milk production, fatty percent, protein percent and somatic cell count traits.

در مطالعه حاضر الگوهای ژنوتیپی IGF-I با مقدار تولید شیر ارتباط معنی­داری نشان دادندP<0.05))، ولی اثر این الگوها بر ترکیبات شیر شامل: درصد چربی، پروتئین و شمار سلول­های بدنی معنی­دار نبود (جدول 2). در میان سایر عوامل ثابت موجود در مدل، اثر سال زایش و ماه رکوردگیری بر تولید شیر معنی­دار بود p<0.01))، (جدول 3). ولی، سن حیوان هنگام زایش اثر معنی­داری روی این صفت نداشت. در این پژوهش رکوردگیری از بزها در ماه­های اسفند، فروردین، اردیبهشت و خرداد سال­های 89، 90 و 91 صورت گرفته، که آنالیز داده­ها نشان داد، سال 90 بیشترین تاثیر و سال 91 کمترین تاثیر را بر روی تولید شیر داشت تاثیر بیشتر سال 90 را می­توان به تغذیه بهتر در آن سال نسبت داد. همچنین اثر ماه رکورد گیری بر مقدار تولید شیر معنی­دار بود و ماه دوم بیشترین تولید شیر و ماه دوازدهم کمترین تولید شیر را داشت. می­توان این اثر را به میزان علوفه و مواد مغذی بهتر در ماه دوم نسبت داد. تعداد روزهای شیردهی که به عنوان عامل کواریت در مدل قرار داده شد تاثیر معنی­داری بر روی مقدار تولید شیر داشت p<0.05)).

جدول 2-  اثر ژنوتیپ­های ژن IGF-I بر میانگین حداقل مربعات تبدیل یافته[7] تولید شیر.

Table 2- The effect of different genotype of IGF-1 gene on least square means milk production.

 *اختلاف معنی­دار در سسطح احتمال 5 درصد.   Significant difference in 5% probability   

جدول 3- اثر عوامل ثابت مؤثر بر میانگین حداقل مربعات تبدیل یافته[8] تولید شیر.

Table 3- The effect of fixed effects on least square means milk production.

** معنی­داری در سطح احتمال 1 درصد      Significant difference in 1% probability

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

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

Table 4- Least square means and standard deviation of birth weight, weaning weight, average daily gain, feed intake and final weight traits.

در مطالعه ارتباط چند شکلی ژن IGF-I با صفت تولید شیر روی بز بومی چین، ارتباط معنی­دار این ژن با مقدار تولید شیر گزارش شد Deng et al., 2010))، که با نتایج  مطالعه حاضر مطابقت دارد. همچنین در مطالعه چند­ ­شکلی ژن فاکتور رشد شبه انسولین (IGF-I) و رابطه آن با تولید شیر روی گاوهای شیری، یک جایگزینی تک نوکلئوتیدی، گزارش شد که با تولید شیر، چربی و پروتئین در ارتباط بود (Mullen et al., 2011). در یک بررسی دیگر روی گاوهای هلشتاین ارتباط معنی­دار این ژن با مقدار تولید شیر، چربی و پروتئین گزارش شد ( Szewczuk et al., 2011).  همچنین یک جهش (C/T) در ناحیه 5´-Flanking این ژن در گاو شناسایی شده است که مرتبط با درصد چربی و پروتئین در شیر بود (Ge et al., 2001)، ولی et al. Hines (1998) در مطالعه­ی ناحیه مذکور در گاوهای هلشتاین هیچ ارتباطی بین جهش رخ داده با صفات تولید شیر پیدا نکردند. نتایج پژوهش حاضر بیانگر ارتباط معنی­دار الگوهای متفاوت ژنتیکی ژن IGF-I با صفت تولید شیر است. ولی اثر چند شکلی این ژن بر روی ترکیبات شیر شامل درصد چربی، پروتئین و شمارش سلول­های بدنی معنی­دار نبود. نقش هورمون رشد شبه انسولین (IGF-I) در رشد بافت های مختلف بدن (سلول های ماهیچه، غضروف و استخوان)، تحریک سنتز پروتئین، افزایش متابولیسم قندها و چربی ها در بدن، محرک تقسیم میتوز، تحریک ازدیاد سلول های مخاطی و استرومال غده پستانی، سنتز و بیان ژن کازئین و حمل گلوکز شناخته شده است و همچنین IGF-I سنتز شیر وتکامل غده پستانی را در طول دوره­ی شیردهی تسهیل می­بخشد، علاوه بر این در طولانی شدن دوره شیردهی نیز نقش دارد (Yinli, 2007 Zhou). با توجه به نتایج مطالعه حاضر و مطالعات گذشته که اثر پلی­مورفیسم ژن  IGF-Iبر صفات تولید شیر را تأیید می­کند، این ژن می­تواند به عنوان نشانگر ژنتیکی برای صفت تولید و ترکیبات شیر استفاده شود.