نوع مقاله : میکروبشناسی و ایمنی شناسی
نویسندگان
1 دانش آموخته دانشکده دامپزشکی دانشگاه تهران، تهران، ایران
2 گروه میکروبیولوژی و ایمونولوژی، دانشکده دامپزشکی دانشگاه تهران، تهران، ایران
3 گروه پاتوبیولوژی، دانشکده دامپزشکی دانشگاه شیراز، شیراز، ایران
چکیده
کلیدواژهها
عنوان مقاله [English]
نویسندگان [English]
BACKGROUND: Major histocompatibility complex (MHC) encodes for highly variable molecules, most of which are responsible for foreign antigen recognition and activation of immune responses in the host. LEI0258 microsatellite, located in the poultry MHC region, is a suitable genetic marker for determining MHC haplotypes and genetic diversity in poultry.
OBJECTIVES: Considering the fact that there is no report on the frequency and types of MHC alleles and population genetic analysis in Ross 308 poultry in Iran, the present study aimed to investigate the diversity of MHC haplotypes of Ross 308 broilers by LEI0258 microsatellite.
METHODS: A total of 216 blood samples were collected from two productive herds of Ross 308 broilers. After extracting DNA of the blood samples and amplifying LEI0258 microsatellite alleles, genotyping of MHC haplotypes was performed using agarose gel electrophoresis and fragment analysis techniques.
RESULTS: A total of seven alleles and 21 genotypes were identified for LEI0258 microsatellite in these two groups. the highest and the lowest frequencies belonged respectively to allele 385 bp (42.86 %) and allele 300 bp (4.33 %). Heterozygous 207/385 was found to be the dominant genotype in both populations. According to the similarity matrix analysis, there was an 84.56 % similarity between the two groups.
CONCLUSIONS: The results obtained in this study revealed a high level of heterozygosity (85.71 % and 91.35 %) and deviation from Hardy–Weinberg equilibrium (P<0.0001) in these two Ross populations. Ross 308 broiler chickens had lower allelic diversity and higher genetic similarity compared to the native ones. These findings provided additional information on the use of MHC as a candidate gene marker in genetic improvement and resource conservation in broiler populations.
کلیدواژهها [English]
تمام مهرهداران به جز ماهیان فاقد آرواره دارای ژنهای Major histocompatibility complex (MHC) یا مجتمع عمده پذیرش بافتی هستند. این ناحیه ژنتیکی توسط Peter Gorer در سال 1936 که بر روی آنتیژنهای گروههای خونی موش کار میکرد شناسایی شد (1). جایگاه ژنهای MHC در گونههای مختلف متفاوت است و در انسان درکروموزوم ٦، موش کروموزوم 17 و در ماکیان در میکروکروموزوم 16 قرار دارد. نامگذاری کلی پروتئینهای MHC نیز در بین گونههای مختلف متفاوت است چنانچه در انسان Human leukocyte antigen (HLA)، در موش H-2 و در ماکیان مجتمع B نامیده میشود (2). MHC دارای 3 دسته جایگاه ژنتیکی است که با توجه به ماهیت و نقشی که ایفا میکنند به سه کلاس II ،I و III طبقه بندی میشوند. مهمترین نقش ژنهای MHC عرضه پادگن به سلولهای حساس به پادگن است. اعمال دیگر نظیر شناسایی و آمادهسازی پادگن نیز بر عهده محصولات ژنی MHC است. مولکولهای MHC کلاس I و II به پپتیدهایی با منشأهای متفاوت متصل میشوند. مولکولهای MHC کلاس I به پپتیدهای مشتق از پروتئینهای سیتوزولی یا داخل سلولی (آنتیژنهای درونزاد) متصل میشوند و آنتیژن را به سلولهایT سیتوتوکسیک عرضه میکنند. مولکولهای MHC کلاس II به پپتیدهای مشتق از پروتئینهای خارج سلولی که توسط فاگوسیتوز یا اندوسیتوز به داخل سلول آورده شدهاند (آنتیژنهای برونزاد) متصل میشوند و آنتیژن را به سلولهای T کمکی عرضه میکنند (1،3). برخلاف پستانداران که ژنهای MHC در آنها در یک منطقه ژنومی بزرگ قرار دارند و از اعضای مشابه چندین خانواده ژنی تشکیل شدهاند، این منطقه در طیور کوچکتر و فشردهتر است. اندازه جایگاه Bدر طیور 92 کیلو باز، شامل 19 ژن و در حدود یک بیستم MHC انسان است. جایگاه B شامل ژنهای کلاس (B-F) I ، کلاس II (B-L) و کلاسIV ( (B-G است (2،4). این اندازه کوچک و سادگی MHC طیور و ارتباط قوی آللهای خاصی در لکوسB با مقاومت و حساسیت به بسیاری از بیماریهای عفونی مهم، پاسخ به واکسن و صفات تولیدی سبب افزایش توجهات به این جایگاه ژنی برای مطالعه روند تکامل در سطح مولکولی و درک عملکرد تغییرات MHC در پدیدههای ایمنی و غیر ایمنی شده است (14-5).
پلیمورفیسم در ناحیه MHC طیور توسط روشهای مختلفی شناسایی میشود. یکی از این روشها استفاده از ریزماهواره LEI0258 است که توسط McConnell و همکاران در سال 1999 شناسایی شد (15). این ریزماهواره تترانوکلئوتیدی است و در جایگاه ژنی B ماکیان بین ژنهای BF و BG واقع شده است (تصویر 1) (16). ریزماهواره LEI0258 اولین بار توسط Fulton و همکاران در سال 2006 برای تعیین هاپلوتیپهای طیور مورد استفاده قرار گرفت (17). ریزماهواره LEI0258 به دلیل داشتن آللهای زیاد و همچنین دامنه وسیع اندازه آللهای آن بسیار جالب توجه است و تفاوت چشمگیر اندازه آللهای این ریزماهواره به آسانی توسط یک الکتروفورز ساده قابل جداسازی و تشخیص است. همچنین این ریزماهواره در عدم تعادل پیوستگی با آللهای MHC است. عدم تعادل پیوستگی معرف یک ارتباط غیرتصادفی از آللها در دو یا چند جایگاه ژنتیکی پیوسته است. به عبارت سادهتر برخی از آللها بیشتر تمایل دارند که همراه یکدیگر در یک هاپلوتیپ حضور یابند. از آنجایی که در تمام دستگاه B این عدم تعادل پیوستگی وجود دارد، نشانگری که در این ناحیه یا یک مکان دیگر از کروموزوم 16 واقع شده باشد، بسیار مفید و کارآمد خواهد بود و امکان تحقیقات بیشتری در رابطه با MHC و کروموزوم 16 را فراهم خواهد ساخت (17). در نتیجه با توجه به اینکه در غالب موارد ارتباط مستقیمی بین یک آلل ریزماهواره LEI0258 و یک هاپلوتیپ خاص از MHC وجود دارد پس این ریزماهواره برای تایپینگ هاپلوتیپهای MHC در طیور مفید است (17). مناسب بودن این مارکر برای تجزیه و تحلیل تنوع ژنتیکی ناحیه MHC در جمعیتهای مختلف طیور در سراسر جهان نشان داده شده است (23-16).
مطالعات نشان میدهند که علاوه بر نقش بسیار ارزشمند ریزماهواره LEI0258 در تعیین ژنوتیپهای MHC، تعدادی از آللهای این ریزماهواره با برخی پارامترهای مربوط به سلامتی از جمله پاسخ پادتن علیه واکسن نیوکاسل، انگلهای دستگاه گوارش، وزن بدن و صفات تولیدی در ارتباط هستند (11،21،24-27). در نتیجه میتوان از این ریزماهواره جهت حفظ و انتخاب این صفات ارزشمند در برنامههای اصلاحنژادی استفاده کرد. هدف از مطالعه حاضر شناسایی هاپلوتیپها و تنوع MHC با استفاده از نشانگر LEI0258 در دو جمعیت مختلف جوجه گوشتی راس 308 است.
نمونهگیری و استخراج DNA: در مطالعه حاضر از تعداد 216 قطعه جوجه گوشتی سویه تجاری راس 308 استفاده شد. 112 جوجه از یک گله مولد و 104 جوجه از گله مولد دیگری تهیه شدند. 1 الی 5/1 میلیلیتر خون از ورید بالی هر جوجه اخذ شد. نمونههای خون در تیوبهای واجد ماده ضد انعقاد EDTA جمعآوری و تا زمان استخراج DNA در دمای 20- درجه سانتیگراد نگهداری شدند. استخراج DNA با استفاده از کیت تجاری iNtRON Biotechnology Genomic DNA Extraction kit و طبق پروتکل کیت صورت گرفت.
سنجش کیفیت DNA استخراج شده: پس از استخراج DNA نمونههای خون، کیفیت نمونهها توسط دستگاه اسپکتروفوتومتر T70+UV/VISمورد بررسی قرار گرفت(PG Instruments, Lutterworth, United Kingdom) .
غلظت DNA و جذب نوری 280/260 نانومتر جهت بررسی خلوص DNA مورد اندازهگیری قرار گرفت. غلظت DNA در تمام نمونهها بالای 50 نانوگرم در میکرولیتر و جذب نوری نزدیک به 9/1 بود.
افزودهسازی آللهای ریزماهواره LEI0258: واکنش PCR در حجم نهایی 25 میکرولیتر شامل 20 نانوگرم DNA ژنومی، 10 پیکومول از هر کدام از آغازگرهای LEI0258، 5/1 میلی مولار MgCl2 و یک واحد آنزیم تک پلیمراز (Fermentas، آلمان) انجام شد. چرخههای دمایی شامل یک مرحله واسرشتگی اولیه 95 درجه سانتیگراد به مدت 2 دقیقه، 30 چرخه سه مرحلهای، شامل مرحله واسرشتهسازی 95 درجه سانتیگراد به مدت 30 ثانیه، مرحله اتصال 63 درجه سانتیگراد به مدت 30 ثانیه، مرحله بسط 72 درجه سانتیگراد به مدت 40 ثانیه و در انتها یک مرحله بسط انتهایی 72 درجه سانتیگراد به مدت 5 دقیقه بود. آغازگرهای مورد استفاده برای افزودهسازی طبق گزارش McConnell و همکاران در سال 1999 (به
شماره دسترسی Z83781 بانک ژن) مورد استفاده قرار گرفتند (15). توالی آغازگرهای رفت و برگشت در جدول 1 نشان داده شده است.
الکتروفورز و آنالیز محصولاتPCR : 5 میکرولیتر از محصول PCR با 2 میکرولیتر از x6 لودینگ بافر مخلوط و روی ژل آگاروز 3 درصد با اختلاف پتانسیل 60 ولت و 25 میلی آمپر به مدت 3 ساعت الکتروفورز شد. اندازه باندهای ایجاد شده در نتایج الکتروفورز بر اساس DNA مارکر 50 جفت باز سنجیده شد. رنگ آمیزی با DNA Safe Stain (شرکت سیناکلون) انجام گرفت و تصاویر ژلها با دستگاه ژل داک EBOX مدل VX5 (Vilber Lourmat, France) ثبت و اندازه هر کدام از محصولات PCR (آللهای ریزماهواره) توسط نرم افزار دستگاه بر اساس جفت باز تعیین گردید. تمام نمونهها با هم مقایسه و نمونههای مشابه مشخص گردیدند.
آنالیز قطعهای آللهای ریزماهواره LEI0258: از آنجایی که ممکن است تفاوت اندازه برخی از آللهای ریزماهواره در حدی باشد که این تفاوت با الکتروفورز بر روی ژل آگارز 3 درصد قابل تشخیص نباشد. لذا برای تأیید نتایج بهدست آمده از الکتروفورز و برای تفکیک دقیقتر آللها از یکدیگر از روش آنالیز قطعهای نیز استفاده شد. جهت آمادهسازی نمونهها ابتدا پرایمر رفت ریزماهواره LEI0258 با رنگ فلورسنت 6-FAM نشاندار شده (رنگ به قسمت' 5 پرایمر متصل شد) و سپس واکنش PCR با شرایط مشابه مرحله قبل تکرار شد. به منظور انجام آنالیز قطعهای، 5/0 میکرولیتر از محصول PCR با 5/0 میکرولیتر از نشانگر وزنی (Gene Scan 500 LIZ Size standadrd) GS500LIZ (Applied Biosystems) و 9 میکرولیتر فرمامید(HiDi Formamide) ترکیب شده و 3 دقیقه در 95 درجه سانتیگراد قرار گرفتند. قطعات DNA تک رشتهای شده بلافاصله بر روی یخ انتقال داده شده و سپس با استفاده از دستگاه تعیین توالی ABI 3130 genetic analyzer تفکیک شدند (Applied Biosystems, Foster City, CA, USA). نتایج حاصل با استفاده از دو نرم افزار Peack scanner نسخه 0/1 و Gene marker نسخه 6/1 آنالیز شدند.
روشهای تحلیل: تعداد آللها و ژنوتیپهای مشاهده شده، فراوانی آللها و ژنوتیپها، محتوای اطلاعات پلیمرفیسم (PIC)، میزان هتروزیگوسیتی و هوموزیگوسیتی مورد مشاهده و مورد انتظار برای هر ژن بهطور جداگانه توسط نرم افزار POP GENE نسخه 32/1 مورد ارزیابی قرار گرفت. میزان تنوع ژنتیکی در جمعیت توسط دو پارامتر تعداد آلل (Na) و میزان هتروزیگوسیتی (He) در هر جایگاه، بر اساس فرمول ارائه شده توسط Nei و همکاران در سال 1973 تعیین شد. انحراف از تعادل هاردی-وینبرگ (HWE) در هر جایگاه با استفاده از آزمون دقیق فیشر بررسی شد.
در بررسی پلیمورفیسم جایگاه MHC جوجه گوشتی راس 308 در گروه اول (112 نمونه) در مجموع 6 آلل و در گروه دوم (104 نمونه) در مجموع 7 آلل متنوع از ریزماهواره LEI0258 در محدوده 195 تا 448 جفت باز شناسایی شد (جدول 2). در گروه 1 بیشترین فراوانی متعلق به آلل 385 جفت باز (86/42 درصد) و آللهای 195 و 300 جفت باز کمترین فراوانی (14/7 درصد) را دارا بودند. با توجه به هم بارز بودن آللهای ریزماهواره از مجموع این 6 آلل در گروه 1 تعداد 13 ژنوتیپ (2 ژنوتیپ هوموزیگوت و 11 ژنوتیپ هتروزیگوت) مشاهده شد که ژنوتیپ 385/207 بیشترین فراوانی (1/26 درصد) و ژنوتیپ 362/362 کمترین فراوانی (89/0 درصد) را داشتند (جدول 3). در گروه 2 بیشترین فراوانی متعلق به آلل 207 جفت باز (21/32 درصد) و آللهای 263 و 300 جفت باز به ترتیب با 37/3 و 33/4 درصد کمترین فراوانی را دارا بودند. از مجموع این 7 آلل در این جمعیت تعداد 17 ژنوتیپ (2 ژنوتیپ هوموزیگوت و 15 ژنوتیپ هتروزیگوت) مشاهده شد که ژنوتیپ 385/207 بیشترین فراوانی (84/28 درصد) و ژنوتیپهای 263/195، 263/263 و 385/263 کمترین فراوانی (96/0 درصد) را داشتند (جدول 3). میزان هتروزیگوسیتی و هوموزیگوسیتی مورد انتظار و مشاهده شده در گروههای 1و 2 در جدول 4 آورده شده است. بهطور کلی نتایج حاصل از مطالعه حاضر سطح بالایی از هتروزیگوسیتی (گروه اول 17/85 درصد و گروه دوم آنچه از یک ریزماهواره معمول انتظار میرود
35/91 درصد) و انحراف از تعادل هاردی-وینبرگ را در این دو جمعیت نشان داد (0001/0P<). تشابه ژنتیکی محاسبه شده بین دو جمعیت نیز 56/84 درصد بود. تصویر2 الکتروفورز آللهای مختلف ریزماهواره LEI0258 بر روی ژل آگارز 3 درصد را نشان میدهد.
بهطور کلی نتایج حاصل از تکنیکهای آگارز ژل الکتروفورز و آنالیز قطعهای با یکدیگر همخوانی داشتند.
تاکنون از روشهای مختلفی جهت مشخص کردن هاپلوتیپهای MHC در طیور استفاده شده است. روش سرولوژی برای اولین بار توسط Briles و همکاران در سال 1982 استفاده شد (28) ولی به علت مشکلاتی از جمله واکنشهای پیچیده و احتمال خطای بالا، زمانبر بودن آزمایش و نبود آنتیسرمهای مربوط به هاپلوتیپهای جدید، روشهای مولکولی از جمله PCR-SSP و RFLP به مرور جایگزین این روش شدند؛ از مزایای روشهای مولکولی، عدم نیاز به آنتیسرمهای مخصوص است اما برای این روشها نیز حضور پروبهای خاص مربوط به ناحیه مورد مطالعه لازم و ضروری است و از سوی دیگر تنوع در سطح جمعیت مربوط به توالیهای نوکلئوتیدی است که آنزیمهای به کار رفته توانایی شناسایی و برش این توالیها را داشته باشند. این تکنیکها همیشه برای آزمایشگاهها عملی نیستند و در تعداد زیاد نمونه کارایی لازم را نداشته و مقرون به صرفه نیز نیستند (17،20).
یکی از روشهای تقریباً نوین در این زمینه استفاده از ریزماهواره LEI0258 است. این ریزماهواره تنوع بسیار بیشتری از را نشان میدهد و علت این تنوع بالا ممکن است مربوط به این باشد که این ریزماهواره در ناحیهای که ژنهای MHC قرار دارند واقع شده است و میتواند تحت تأثیر بهگزینی تنوع دهنده (Diversifying selection) قرار گیرد. این نوع از انتخاب، تنوع را نه تنها در ژنهای تحت تأثیر آن بلکه در نواحی مجاور هم ایجاد میکند (17،23).
سهولت استفاده، تنوع قابل تشخیص آللهای ریزماهواره LEI0258 و عدم تعادل پیوستگیاش با سایر آللهای MHC، امروزه این ریزماهواره را به مرکز توجهات جهت تایپینگ MHC طیور تبدیل کرده و تاکنون در بررسی تنوع MHC-B در نژادهای بومی و تجاری کشورهای مختلف از جمله چین، هند، ایران، کره، آرژانتین، تانزانیا، کامرون و نیجریه استفاده شده است (39-16،22،23،37). در ایران تاکنون گزارشی مبنی بر تحلیل تنوع ریزماهواره LEI0258 برای سویه راس صورت نگرفته است. در مطالعه حاضر جهت شناسایی آللهای ریزماهواره LEI0258 از تکنیک آگارز ژل الکتروفورز و آنالیز قطعهای استفاده شد. در مطالعات مشابهی در طیور بومی ایران، نیجریه و آرژانتین از آگارز ژل الکتروفورز و همچنین در مطالعات متعددی از روش آنالیز قطعهای برای شناسایی آللهای ریزماهواره LEI0258 با موفقیت استفاده شده است (21،23،30،35).
Haunshi و همکاران در سال 2020 تنوع ژنتیکی MHC را در دو جمعیت طیور بومی هند (Ghagus و Nicobari) در مقایسه با نژاد White Leghorn (WLH) با استفاده از نشانگر LEI0258 مورد بررسی قرار دادند. در نژادهای Ghagus، Nicobari و WLH به ترتیب 23، 14 و 10 آلل از این ریزماهواره شناسایی شد که از این میان، آللهای 195 و 385 جفت باز که به ترتیب در Ghagusو WLH یافت شدند در مطالعه حاضر نیز شناسایی شدند (21).
Duan و همکاران در سال 2019 توزیع آللی LEI0258 را در 9 جمعیت مختلف jungle fowl بهدست آوردند. در مجموع کل آللهای شناسایی شده در این جمعیتها، سه آلل 207، 263 و 448 جفت باز در مطالعه حاضر نیز شناسایی شدند (19).
مطالعه دیگری توسط Ncube و همکاران در سال 2014 بر روی چهار جمعیت طیور روستایی آفریقای جنوبی، زیمباوه و مالاوی انجام شد. در مجموع 22 آلل در محدوده 448-188 جفت باز مشاهده شد و به طور متوسط 12 آلل در هر جمعیت تشخیص داده شد. در دو جمعیت مورد مطالعه از آفریقای جنوبی، آلل 448 جفت باز با درصد فروانی کمتر از 5 درصد شناسایی شد (34). همچنین Touko و همکاران در سال 2015 نیز به بررسی تنوع ژنتیکی MHC-B طیور بومی کامرون از طریق تایپینگ مولکولی با نشانگرهای ریزماهواره LEI0258 و MCW0371 پرداختند، آلل 448 جفت باز برای ریزماهواره LEI0258 یافت شد (37). در مطالعه حاضر نیز آلل 448 جفت باز شناسایی شد.
Mwambene و همکاران در سال 2019 مطالعهای بر روی 10 اکوتیپ بومی تانزانیا انجام دادند. از مجموع 30 آلل شناسایی شده، آلل 300 جفت باز در تمام جمعیتها و آلل 385 جفت باز در هشت جمعیت مشاهده شدند و در مطالعه حاضر نیز این آللها شناسایی شدند (22).
اگرچه تعداد آللهای مشاهده شده در مطالعه حاضر در محدوده تعداد آللهای گزارش شده در سایر جمعیتهای تجاری است، اما کمتر از تعداد آللهایی است که در جمعیتهای طیور بومی کشورهای مختلف گزارش شده است. در مطالعهای که توسط Nikbakht و همکاران در سال 2013 بر روی سه جمعیت مرندی، خراسان و آرین از ایران انجام شد، 23 آلل در جمعیت مرندی، 15 آلل در جمعیت خراسان و 11 آلل در جمعیت گوشتی آرین تشخیص داده شد (23). همچنین در مجموع 42، 37 و 30 آلل از این ریزماهواره به ترتیب در طیور بومی کامرون، هند و تانزانیا گزارش شده است (21،22،37). تنوع ژنتیکی بالاتر طیور بومی در مقایسه با طیور تجاری میتواند به این دلیل باشد که طیور بومی اغلب در شرایط آزاد نگهداری میشوند و تحت کنترل شدید برنامههای مدیریتی اصلاح نژاد نیستند. طیور بومی با پاتوژنهای متعددی مواجه شده و شیوع بیماریها در آنها بیشتر است و در نتیجه افزایش تنوع به آنها این امکان را میدهد که به پاتوژنهای محیطی بیشتری پاسخ داده و بخت زندهمانی بیشتری داشته باشند. در حالی که طیور تجاری اغلب در شرایط بهداشتی و پرورشی دقیق به منظور بهبود صفات تولیدی قرار داشته و در این بهگزینی، مقاومت در برابر بیماریها انتخاب دوم است. جمعیتهای تجاری تنوع ژنتیکی بسیار کمی داشته و آللهای محدودتری از ریزماهواره LEI0258 در مقایسه با طیور بومی را نشان میدهند (38). در این مطالعه در بررسی دو جمعیت طیور راس، در یک گروه آلل 263 جفت باز شناسایی نشد و در گروه دیگر نیز این آلل کمترین فراوانی را داشت که این امر میتواند نشانهای از کاهش تنوع در این جمعیت به علت از بین رفتن این آلل نایاب باشد لذا ضروری است در مطالعات آتی این امر مورد توجه قرار گیرد.
نتیجهگیری نهایی: مطالعه حاضر با هدف ارزیابی تنوع آللی و ژنتیکی در دو جمعیت راس 308 با استفاده از نشانگر LEI0258 انجام شده است. از ریزماهواره LEI0258 میتوان به عنوان یک نشانگر DNA در مطالعات ایمنیشناسی، ژنتیکی و همچنین به عنوان یک روش مناسب برای تایپینگ MHCدر طیور استفاده کرد. استفاده از این نشانگر برای طیور میتواند اطلاعات ارزشمندی در رابطه با تنوع و فراوانی آللهای موجود در ذخایر ژنتیکی کشور به دست دهد و امکان مقایسه آن با سایر جمعیتهای جهان را برای ما فراهم سازد. جمعیتهای راس 308 در مقایسه با جمعیتهای بومی ایران از تنوع آللی کم و تشابه ژنتیکی بالایی برخوردارند. تنوع کم در آللها نشاندهنده اصلاح نژاد و انتخاب پرنده به منظور تولید گوشت است که ممکن است به حساسیت به بیماریها منجر شود. از سوی دیگر تشابه فراوانی آللی بین دو جمعیت حاکی از اجداد مشابه است که میتواند به دلیل انحصار تولید سویه راس رخ دهد. یافتههای مطالعه حاضر اطلاعات بیشتری جهت معرفی آللهای سویه راس و استفاده از MHC به عنوان یک نشانگر ژنی برای حفظ منابع و بهبود ژنتیکی در جمعیت جوجههای گوشتی فراهم میکند.
نویسندگان از معاونت پژوهشی دانشگاه تهران و دانشکده دامپزشکی جهت تأمین اعتبار این مطالعه در راستای طرح پژوهشی مصوب به شماره ۷۵۰۲۰۱۵/۶/39 قدردانی میکنند.
بین نویسندگان تعارض در منافع گزارش نشده است.