نوع مقاله : انگل شناسی و بیماری های انگلی
نویسندگان
1 دانش آموخته علوم درمانگاهی، دانشکده دامپزشکی، دانشگاه فردوسی مشهد، مشهد، ایران
2 گروه پاتوبیولوژی، دانشکده دامپزشکی، دانشگاه فردوسی مشهد، مشهد، ایران
چکیده
کلیدواژهها
عنوان مقاله [English]
نویسندگان [English]
BACKGROUND: Cystic echinococosis (CE) is known to be one of the most important zoonotic diseases in different parts of Iran. Even though it causes major health problems, there is limited information regarding its transmission cycles and strain of this infection in eastern Iran.
OBJECTIVES: The present study aimed to characterize the strain of Echinococcus granulosus cysts in the slaughtered sheep and goats in Birjand area using morphological and molecular criteria.
METHODS: Isolates of E. granulosus were collected from sheep (30) and goats (30) from Birjand slaughterhouse and characterized employing both DNA (PCR-RFLP of ITS1) and morphological criteria (metacestode rostellar hook dimensions). In addition, the fragments of the genes coding for ITS-1 were sequenced.
RESULTS: Among the two different identified strains/genotypes (sheep and camel), the sheep strain appeared to be the most common genotype of E. granulosus affecting sheep and goats. All of the 30 sheep samples and 20 out of 30 goat samples were infected with sheep strain. However, the camel genotype was only observed in the goats and 10 out of 30 goat isolates were infected with the camel genotype. The camel genotypes had RFLP patterns, which were different from the RFLP patterns of the rest of isolates (sheep strain).
CONCLUSIONS: The results of this study revealed that the ‘sheep’ strain was the most prevalent strain in sheep and goats in this area. Moreover, the camel genotype (G6) was confirmed to trigger infection in the slaughtered goats of Birjand area.
کلیدواژهها [English]
کیست هیداتید بیماری زئونوزی است که توسط مرحله نوزادی سستود سگ متعلق به جنس اکینوکوکوس (خانواده تنیده) ایجاد میشود. هیداتیدوز از بیماریهای بومی کشور هست و حیوانات اهلی زیادی ازجمله گوسفند، بز، شتر و گاو میزبانان واسط آن میباشند (5،12،14). بیماری در انسان موجب ضایعات و آسیبهای بافتی شده و گاهی منجر به مرگ میگردد. هیداتیدوزیس همچنین در حیوانات بهخصوص دامها موجب زیانهای اقتصادی فراوان میشود. سگ، میزبان اصلی بیماری، با دفع تخم انگل همراه با مدفوع در مزارع، سبب آلودگی آب، گیاهان و سبزیها میشود. انسان در اثر تماس نزدیک با سگ و مصرف سبزیها و آب آلوده و نشخوارکنندگان با تغذیه از علوفه آلوده به تخم، به بیماری مبتلا میگردند. کیست هیداتید یک مشکل اقتصادی و بهداشتی برای کشورهایی است که بیماری بهصورت بومی در آنجا وجود دارد. کیست هیداتید در بسیاری از مناطق آسیا، اروپا، آمریکای جنوبی، خاورمیانه، استرالیا و نیوزلند بهصورت بومی وجود دارد (21).
اکینوکوکوس گرانولوزوس واجد سویههای فراوانی است. این سویهها از نظر ریختشناسی، همهگیریشناسی، درمان و کنترل باهم اختلاف دارند بهطوریکه این سویهها ویژگیهای آسیبشناسی، اختصاصات میزبانی، بیولوژی تکامل، فیزیولوژی، ژنتیک، بیماریزایی انگل برای میزبان، الگوی چرخه زندگی انگل، آنتی ژنیسیتی، روند انتقال آلودگی و حساسیت آنها نسبت به عوامل دارویی را در مواردی تحت تأثیر قرار میدهند. تمامی موارد فوق در تهیه واکسنها، تشخیص، پاسخ درمانی، همهگیریشناسی و کنترل بیماری تأثیر گذارند (20). اکینوکوکوس گرانولوزوس دارای 10 واریته (G1-10) است. رایجترین سویه زئونوز انسانی در اکثر کشورها شامل جنوب و شرق آفریقا، قسمتهایی از آسیا، استرالیا و جنوب آمریکا (آرژانتین) سویه G1 (Sheep strain) است. این سویه علاوه بر انسان در گوسفند، گاو، بز، بوفالو، شتر و خوک نیز گزارش شده است. سویهG2 (Tasmanian sheep strain) در گوسفند و انسان، سویه buffalo strain)) G3 در بوفالو، سویهی (horse strain =Echinococcus equinus) G4 در اسب و الاغ، سویه G5 (E. ortolepi) در گوسفند، بز، گاو، بوفالو، خوک و انسان، سویه G6 (E. intermedius) در شتر، گاو، انسان، گوسفند و بز، سویه G7 (E. intermedius) در خوک، گاو، بز و انسان، سویه (G8 (cervid strain در گوزن و انسان، سویه G9 در انسان و سویه (G10 (fennoscandian strain در گوزن گزارش شده است (19،21).
امروزه از تکنیکهای مختلفی از جمله روشهای ریختشناسی و مولکولی بهمنظور تشخیص سویههای اکینوکوکوس گرانولوزوس استفاده میشود. در روش ریختشناسی، سویههای اکینوکوکوس گرانولوزوس بر اساس اندازه و تعداد قلابهای روستلوم پروتواسکولکسها مشخص میشوند. پروتواسکولکسهای داخل هر کیست حاوی 2 ردیف قلاب و هر قلاب دارای یک تیغه و یک دسته و گارد است که تعداد و طول و عرض قلابها در شناسایی سویههای انگل اهمیت دارد (4). از جمله تکنیکهای مولکولی استفاده شده برای تعیین سویه شامل: تعیین ترادف ژنهای میتوکندریال cox-1 وnad-1 و آنالیز مناطق DNA ریبوزومال شاملITS1 و ITS2 با تکنیکهای PCR-RFLP و PCR-RAPD و SSCP هست (19).
گوسفند و بز نقش مهمی در چرخه انتقال هیداتیدوزیس دارند. مطالعات انجام شده در نواحی مختلف ایران بیانگر این است که سویه G1 شایعترین سویه در دامهای ایران میباشد (6،7،8،12). در مطالعات انجام شده اخیر در بزهای نواحی مرکزی و شمال شرق کشور بیانگر حضور سویه G6 , G7 در این دام میباشد (6،14). اما در خصوص پراکندگی سویههای اکینوکوکوس گرانولوزوس در گوسفندان و بزهای استان خراسان جنوبی تنها یک مطالعه روی تعداد معدودی کیست هیداتید جدا شده از گوسفندان این منطقه انجام شده است ولی روی بز مطالعهای انجام نشده است (9). لذا هدف از انجام مطالعه حاضر بررسی سویههای کیست هیداتید جدا شده از گوسفندان و بزهای کشتار شده در کشتارگاه بیرجند واقع در استان خراسان جنوبی با دو روش ریختشناسی و بررسی شاخصهای مورفومتریک قلابهای پروتواسکولکس و روش مولکولی با روشPCR-RFLP ناحیهی rDNA-ITS1 هست.
نمونهگیری: منطقه مورد مطالعه شهرستان بیرجند مرکز استان خراسان جنوبی در ارتفاع ۱۴۷۰ متری از سطح دریا قرار گرفته است و جمعیت آن 020/178 نفر میباشد. بیرجند، شهرستانی کوهستانی است و در آن کوهها و درههای عمیق و حاصلخیزی وجود دارد. آب و هوای این شهرستان بیابانی و نیمه بیابانی است و در این شهرستان رودخانه عمدهای وجود ندارد و متوسط بارندگی سالیانه به 150 میلیمتر میرسد. این استان ازنظر جمعیت دامی دارای 1 میلیون و ۷۵۰ هزار رأس گوسفند و بز است.
جهت انجام مطالعه حاضر طی دوره زمانی بهمن 94 تا تیر 96، تعداد 30 عدد کبد و ریه گوسفند و تعداد 30 عدد کبد و ریه بز آلوده به کیست هیداتید از کشتارگاه بیرجند جمعآوری گردید و به آزمایشگاه انگلشناسی دانشکده دامپزشکی دانشگاه فردوسی مشهد انتقال یافت.
بررسیهای مورفولوژیکی: در شرایط آسپتیک با استفاده از سرنگ استریل 10 میلی لیتر مایع درون کیستها را کشیده و در لولههای آزمایشگاهی خالی گردید. لولههای حاوی کیست هیداتید جمعآوری شده از هر عضو (کبد و ریه) با دور 800 به مدت 2 دقیقه سانتریفوژ گردید. سپس مایع رویی خارج گردید و پروتواسکولکسهای موجود در ته لولهها جمعآوری شد و با استفاده از پیپت پاستوراستریل تعدادی از پروتواسکولکسها روی لام قرار داده شد و پس از اضافه نمودن لاکتوفنل لامل قرار داده شد. با استفاده از عدسی چشمی حاوی لام میکرومتری شاخصهای ریختشناسی قلابهای روستلوم هر پروتواسکولکس شفاف شده شامل: طول و عرض قلابها، طول و عرض تیغه و طول و عرض دسته هر قلاب اندازهگیری گردید. پروتواسکولکسهای باقیمانده در ته لوله پس از چند بار شستشو با آب مقطر استریل در میکرو تیوبهای 5/1 میلی لیتری به همراه مقداری اتانول 70 درصد تا زمان استخراج DNA در20- درجه نگهداری شد (تصویر 1).
استخراج DNA: برای استخراج DNA از مخلوط پروتواسکولکسهای جدا شده از هر کیست (در حدود 300 پرتواسکولکس) استفاده گردید.
برای استخراج DNA با استفاده از کیت MBST (ساخت کشور ایران) طبق دستورالعمل شرکت سازنده صورت گرفت. پس از استخراج DNA غلظت DNA تخلیص شده توسط نانو دراپ تعیین شد که غلظت آن 100 نانوگرم در لیتر بود.
واکنش زنجیرهای پلی مرازPCR : بعد از انجام استخراج DNA، واکنش PCR بر روی ژن ITS1 انجام شد. برای انجام واکنش PCR در حجم نهایی 25 میکرولیتر به میزان 5 میکرولیتر از الگوی DNA استخراجشده، 5 /12 میکرولیتر مستر میکس (شرکت سینا ژن)، 1 میکرولیتر پرایمر Forward ،1 میکرولیتر پرایمر Reverse و 5/5 میکرولیتر آب مقطر دیونیزه اضافه گردید. در هر واکنش از یک کنترل مثبت و یک کنترل منفی (فاقد الگوی DNA) استفاده شد. پرایمر مورد استفاده از ژن ITS1 بر اساس مطالعه انجام شده توسط Moghaddas و همکاران در سال 2014 (13) استفاده گردید.
برنامه PCR جهت تکثیر ژن ITS1 شامل واسرشت اولیه 95 درجه سانتیگراد به مدت 5 دقیقه، 35 سیکل شامل واسرشت 94 درجه سانتیگراد مدت 45 ثانیه، اتصال پرایمر 50 درجه سانتیگراد به مدت 45 ثانیه، تکثیر 72 درجه سانتیگراد به مدت 45 ثانیه و سیکل تکثیر نهایی 72 درجه سانتیگراد به مدت 10 دقیقه بود. بعد از انجام PCR، محصول واکنش بر روی ژل 7/1 درصد آگارز مورد بررسی قرار گرفت.
PCR- RFLP بر روی ژن ITSI: برای تعیین سویههای کیست هیداتید از روش هضم آنزیمی با استفاده از آنزیم Bsh 1236I استفاده گردید (13). لازم به ذکر است که آنزیم برش دهندهیBsh1236I(BST UI) توالی CGCG را شناسایی کرده و محل بین GC را در توالی′ 5′…CG⁄ CG…3 برش میدهد.
مواد تشکیلدهنده واکنش PCR- RFLP بر روی ژن ITS1 در حجم کلی واکنش 25 میکرولیتر شامل 10 میکرولیتر محصول واکنش PCR، 2 میکرولیتر آنزیمBsh 1236I ، 2 میکرولیتر بافر 10x Buffer و 16 میکرولیتر آب مقطر دیونیزه است. سپس نمونهها به مدت 16 ساعت در دمای 37 درجه سانتیگراد قرار داده شد. بعد از 16 ساعت برای غیرفعال شدن آنزیم نمونهها در دمای 65 درجه سانتیگراد به مدت 20 دقیقه قرار گرفت. سپس جهت مشاهده چگونگی برشها و وزن باندها و تشخیص سویهها بر اساس آن، نمونهها بر روی ژل آگارز 3 درصد قرار داده شد. سپس با استفاده از دستگاه UV-Eluminator باندهای DNA (محصولاتPCR) بر روی ژل بررسی و از آنها عکس گرفته شد.
تعیین توالی: تعداد 5 نمونه از 60 نمونه PCR شده برای تأیید نتایج PCR-RFLP بهمنظور تعیین توالی به شرکت Bioneer کره فرستاده شد. آنالیز توالیهای سکانس شده توسط نرمافزار CLC Main Workbench صورت گرفت.
نتایج بررسی ریختشناسی کیست هیداتید: نتایج بررسی شاخصهای ریختشناسی و مورفومتریک قلابهای روستلوم پروتواسکولکسها بیانگر حضور یکسویه در گوسفندان بود و مشابه با شاخصهای مورفولوژی سویه گوسفندی بود. ولی در بزهای کشتار شده دو نوع شاخص قلاب شناسایی گردید که از لحاظ شاخصهای مورفولوژی و مورفومتریک تفاوت معنیداری با یکدیگر داشتند (جدول1). شاخصهای مورفولوژی در قلابهای نوع اول مشابه سویه گوسفندی بود درحالیکه شاخصهای مورفولوژی قلابهای نوع دوم مشابه سویه شتری بود.
نتایج مولکولی (نتایج بررسی مولکولی پروتواسکولکسهای استخراج شده از کیستهای هیداتید): اندازه قطعه تکثیر شده بعد از انجام واکنش PCR بهوسیله پرایمرهای اختصاصی (Eg-F, Eg-R) و پس از الکتروفورز 462 جفت باز است (تصویر 2).
نتیجه واکنش PCR-RFLP: در30 نمونه بررسی شده از کیستهای هیداتید گوسفندان در مطالعه حاضر بر اساس الگوی برشی با آنزیم BSH1236I ، الگوی 4 باند مشاهده گردید. با توجه به مطالعات قبلی در این زمینه سویه تعیین شده (بر اساس تعداد و وزن باندها) سویه G1-G3 است (تصویر 3).
در 20 نمونه بررسی شده از کیستهای هیداتید بزهای کشتار شده در مطالعه حاضر بر اساس الگوی برشی با آنزیم BSH1236I، الگوی 4 باند مشاهده گردید که بر اساس تعداد و وزن باندها سویه G1-G3 است. در 10 نمونه از کیستهای هیداتید بز الگوی 3 باند مشاهده گردید که بر اساس تعداد و وزن باندها سویه G6/G7 است (تصویر 3).
نتیجه تعیین توالی: پس از اصلاح کردن توالیها در مرحله بعد توالیها در NCBI بلاست گردید. مقایسهی توالیها با توالیهای ژن ITS1 نشاندهنده وجود سویه G1 در کل نمونههای جدا شده از کیستهای هیداتید گوسفندان و 20 عدد از بزهای کشتار شده و سویه G6 در تعداد 10 عدد از بزهای کشتارگاه بیرجند است.
با توجه به تنوعهای درونگونهای اکینوکوکوس گرانولوزوس اکتفا نمودن به اطلاعات به دست آمده از شاخصهای ریختشناسی جهت تمیز سویه و گونه، قابل اعتماد نیست. بنابراین تلفیقی از روشهای ریختشناسی و مولکولی با هم میتواند اطلاعات مقبولتری را در رابطه با شناسایی تنوع در سویههای اکینوکوکوس گرانولوزوس فراهم کند.
ویژگیهای ریختشناسی اگر چه دارای محدودیت هستند، اما هنوز هم بهعنوان ابزاری مهم در شناسایی و تاکسونومی انگل اکینوکوس گرانولوزوس شناخته میشوند. بررسیهای ریختشناسی یک روش سریع و آسان برای شناسایی سویهها در زمینة مطالعات اپیدمیولوژیکی میباشند. در مطالعه حاضر 6 شاخص ریختی قلابهای روستلوم از قبیل تعداد کلی قلابهای کوچک و بزرگ، طول قلاب کوچک و بزرگ، طول تیغه قلاب کوچک و بزرگ، طول دسته قلاب کوچک و بزرگ و عرض کلی قلاب کوچک و بزرگ مورد بررسی قرار گرفت. نتایج مورفولوژی و مورفومتریک در کل کیستهای هیداتید جدا شده از گوسفندان دارای الگوی مشابهی بودند و بیانگر این بود که همه این کیستها ازنظر سویه یکسان هستند. از طرف دیگر در کیستهای هیداتید جدا شده از بزها پارامترهای ریختشناسی یکسان نبودند و تعداد 20 کیست دارای الگوی مشابه بودند که از نظر این پارامترها با کیستهای هیداتید گوسفندان تفاوت معنیداری نداشتند و بیانگر سویه مشابه با گوسفندان بودند ولی تعداد 10 کیست باقی مانده الگوی متفاوتی با بقیه داشتند که بیانگر حضور سویه متفاوت با قبلیها بود. بر اساس نتایج ریخت شناسی میتوان به این نتیجه رسید که در گوسفندان کشتار شده در ناحیه مورد مطالعه یک سویه ولی در بزهای این منطقه دو سویه حضور دارد.
نتایج روش مولکولی مورد استفاده در این مطالعه با روش PCR-RFLP روی قطعه ITS1 و با استفاده از آنزیم Bsh1236I با نتایج مورفولوژی همخوانی داشتند و الگوی هضم آنزیمی در کل گوسفندان این منطقه یکسان بود و بیانگر حضور سویه G1-G3 بود. از طرف دیگر در کیستهای هیداتید جدا شده از بزهای این منطقه دو نوع الگوی برش مشاهده گردید. الگوی برش تعداد 20 کیست جدا شده از بزها با الگوی برش در کیستهای هیداتید گوسفندان یکسان بود که با توجه به الگوی وزن باندها سویه G1-G3 بود. در حالی که الگوی برش در تعداد 10 کیست هیداتید جدا شده از بزها با بقیه کیستهای جدا شده از گوسفندان و بزها یکسان نبودند که بر اساس الگوی وزن باندها سویه G6/G7 تشخیص داده شد. بر اساس تعیین توالی انجام شده از این الگوها در کل کیستهای هیداتید جدا شده از گوسفندان و تعداد 20 عدد از کیستهای هیداتید جدا شده از بزهای این ناحیه سویه G1 و در تعداد 10 کیست هیداتید جدا شده از بزها سویه G6 مورد تأیید قرار گرفت.
نتایج حاصل از مطالعه حاضر با مطالعه انجام شده قبلی در این منطقه که روی تعداد 3 کیست هیداتید جدا شده از گوسفندان انجام شده بود مشابهت دارد که بیانگر حضور سویهG1 در گوسفند است (9) ولی روی کیستهای هیداتید جدا شده از بز در این منطقه، مطالعه حاضر اولین مطالعه است. البته در مطالعه قبلی سویه G6 از تعداد 5 شتر و 8 انسان گزارش شده است که بیانگر حضور چرخه آلودگی با سویه G6 در انسان و دامهای این منطقه است و با توجه به حضور سویهG6 در انسانهای این منطقه نقش بزها در انتقال آلودگی توجیه پذیر است.
بیشترین مطالعات انجام شده در ایران روی کیستهای هیداتید جدا شده از گوسفندان انجام شده است و همه این مطالعات بیانگر این است که گوسفند مهمترین میزبان واسط برای سویه G1 است. ازجمله در مطالعه Fadakar و همکاران در سال 2015 (6) به روش ITS1-RFLP در 30 نمونه گوسفند در استان خراسان رضوی و شمالی، مطالعه Khademvatan و همکاران در سال 2013 (11) بر روی 141 نمونه گوسفند در استان خوزستان با تکنیک ITS1-RFLP، Shahnazi و همکاران در سال 2011 (17) در منطقه اصفهان، Rostami Nejad و همکاران در سال 2008 (16) در گوسفندان نواحی شمالی و غرب کشور Karimi and Dianatpour در سال 2008 (10) با روش nad1-RFLP در گوسفندان منطقه استان فارس و در مطالعه Zhang و همکاران در سال 1998 (22) از مناطق مختلف ایران، سویه G1 را گزارش کردند. اما در بعضی مطالعات ازجمله، مطالعه Hajialilo و همکاران در سال 2012 (8) در منطقه جنوب شرقی ایران و مطالعه Abedi و همکاران در سال 2019 (1) در استان مرکزی در 86-61 درصد گوسفندان سویه G1 و در بقیه گوسفندان سویه G3 (37 درصد) گزارش شده است. البته در مطالعهFasihi Harandi و همکاران در سال 2002 (7) در استان کرمان به روش ITS1-RFLP در 40 نمونه گوسفند سویه G1 و در 3 نمونه گوسفند سویه G6 پیدا کردند. در مطالعه دیگر انجام شده توسط Ahmadi and Dalimi در سال 2006 (2) با استفاده از 4 آنزیم به روش ITS1-RFLP دوسویه G1 و G6 را در نمونههای گوسفند و شتر گزارش کردند.
از طرف دیگر مطالعات کمتری روی کیستهای هیداتید بز در نواحی مختلف ایران انجام شده است و اگرچه مشابه گوسفند سویه G1 بیشترین فراوانی را در بزها دارد ولی نکته قابل توجه در این مطالعات گزارش وجود سویه G6 و G7 از بزها میباشد. ازجمله این مطالعات میتوان به مطالعه Fadakar و همکاران در سال 2015 (6) که در 24 نمونه بز در استان خراسان رضوی و شمالی سویه G1 و در 4 نمونه بز از خراسان شمالی سویه G7 پیدا کردند. Rajablo و همکاران در سال 2012 (14) در مطالعه روی سویههای کیست هیداتید در بزهای نواحی مختلف ایران تعداد 17 کیست از نوع سویه G1 و تعداد 3 کیست سویه G6 بودند که با نتایج حاصل از این طرح همخوانی دارد. با توجه به اینکه شتر میزبان اصلی سویه G6 است این سویه در مناطق مرکزی و نواحی مرزی و بیابانی که عمدتاً محل پرورش و نگهداری از شتر است دیده میشود. نظر به پرورش شتر در منطقه استان خراسان جنوبی حضور این سویه در بزهای منطقه قابل توجیه است و لذا امکان دارد که بزها در سیر تکاملی و پراکندگی این سویه هم نقش داشته باشند.
اکینوکوکوس گرانولوزوس دارای 10 سویه هست، که 7 سویه انسان را درگیر میکند. رایجترین سویهی زئونوز انسانی در اکثر کشورها شامل جنوب و شرق اروپا، شمال و شرق افریقا، قسمتهایی از آسیا، استرالیا و جنوب امریکا (آرژانتین) سویه G1 هست (5). مهمترین میزبان این سویه گوسفند است ولی این سویه در گاو، بز، شتر، گاومیش و خوک نیز گزارش شده است. در ایران هم مطالعاتی مبنی بر وجود هیداتیدوز در سه سیکل مجزا وجود دارد: 1- سیکل اهلی بین سگها و چهارپایان2- سیکل بیابانی بین سگها و شتر 3- سیکل وحشی بین گوشتخواران وحشی و نشخوارکنندگان وحشی. در سیکل اهلی مهمترین میزبان واسط گوسفند هست، که در بین سایر دامها نیز از نظر میزان باروری کیستها بیشترین درصد باروری کیست (88 درصد) را به خود اختصاص میدهد (4). با توجه به نتایج مطالعه حاضر، گوسفند نقش مهمی در بقای چرخه انتقال هیداتیدوز بین انسان و دام دارد. با توجه به نتایج بهدستآمده به نظر میرسد که سویه G1 اکینوکوکوس گرانولوزوس موجود در گوسفندان استان به عنوان مهمترین سویه آلودهکننده انسان مطرح است. از طرف دیگر با توجه به مطالعه آلودگی انسان به سویه G6 از منطقه بیرجند و حضور این سویه در بزهای منطقه، نقش این دام در انتقال آلودگی و تداوم چرخه آلودگی زندگی انگل در این منطقه قابلتوجه است (9). با در نظر داشتن چرخه زندگی و دوره کوتاهتر بلوغ سویه G6 در مقایسه با سویهG1 در میزبان نهایی لازم است در اقدامات بهداشتی جهت کنترل و پیشگیری بیماری در انسان به این نکته توجه گردد.
این پایاننامه با حمایت مالی معاونت پژوهشی دانشگاه فردوسی مشهد انجام گرفت و نویسندگان مطالعه حاضر از معاونت پژوهشی دانشگاه فردوسی مشهد به دلیل حمایت مالی این پایاننامه تشــکر و قدردانی مینمایند.
بین نویسندگان تعارض در منافع گزارش نشده است.
References