Abstract
Cucurbita genus belongs to Cucurbitaceae family. These plants are attacked with different pathogenic agents. Currently the effective strategy to control many disease of Cucurbitaceae is the usage of resistant plants. In this study, the NBS domain encoded by resistance genes was studied in Iranian native varieties of cantaloupe. For this purpose, the seeds of Iranian native cantaloupe cultivars were provided and cultivated in greenhouse. DNA was extracted from leaves of different cultivars using CTAB method. Degenerate primers were designed from conserved motives of resistance analogues genes and amplification of these genes was performed using PCR method. The results showed the isolation of an analogous for resistance gene to tomato mosaic virus (ToMV) in cantaloupe cultivars TN-92-99 and TN-92-80 that were cloned in pGEM-T plasmid and then sequenced. Blast analysis indicates that NBS sequence in Irnian cantaloupe share 92% similarity with ToMV resistance gene in watermelon. Phylogenetic tree drawing using resistance gene analogues from cantaloupe and watermelon existing in NCBI gene bank created tow main groups and seven subgroups.Bioinformatic analysis revealed that these genes sequences are conserved among some native cantaloupe caltivars of Iran. Second structure study of protein of ToMV resistance gene analogus using PSIpred software indicates that these proteins only have α helix structure. This protein is secreted in plasma membrane and its extracellular secretion is very little.
Keywords
جداسازی ، همسانهسازی، توالی یابی و بررسی بیوانفوماتیکی آنالوگ ژن مقاومت به ویروس موزاییک گوجه فرنگی در دو رقم طالبی بومیایران
فاطمه قرائی1و2 و مریم غایب زمهریر1*
1 ایران،تهران، سازمان تحقیقات، ترویج و آموزش کشاورزی، مؤسسه تحقیقات گیاهپزشکی کشور، بخش تحقیقات بیماریهای گیاهی
2 ایران، تهران، دانشگاه آزاد اسلامی واحد علوم و تحقیقات، دانشکده کشاورزی ومنابع طبیعی، گروه باغبانی
تاریخ دریافت: 23/5/95 تاریخ پذیرش: 13/8/96
چکیده
جنس Cucurbita از خانواده کدوئیان (Cucurbitacea) میباشد. این گیاهان، مورد حمله طیف وسیعی از عوامل بیماریزا قرار میگیرند. در حال حاضر تدبیر کارآمد برای کنترل بسیاری از بیماریهای کدوئیان، استفاده از گیاهان مقاوم میباشد. در این مطالعه، قلمرو NBS کد شونده با ژنهای مقاومت در ارقام طالبی بومی ایران بررسی شد. به این منظور بذور ارقام مختلف طالبی ایرانی تهیه و در گلخانه کشت داده شدند. استخراج DNA به روش CTAB از برگ ارقام مختلف انجام شد. آغازگرهای دجنره از نواحی حفاظت شده ژنهای آنالوگ مقاومت طراحی شدند و تکثیر این ژنها از روی DNA به روش PCR انجام شد. نتایج این بررسی منجر به جداسازی یک ژن آنالوگ مقاومت به ویروس موزاییک گوجه فرنگی (ToMV) از طالبی رقم TN-92-99 و TN-92-80 شد که در پلاسمید pGEM-T همسانه سازی و توالی یابی شد. آنالیز بلاست نشان داد که توالی NBS در طالبی ایرانی 99-92 درصد با ژن مقاومت به ویروس موزاییک گوجه فرنگیدر سایر کدوئیان شباهت دارد. ترسیم درخت فیلوژنی با استفاده از آنالوگهای جداسازی شده از طالبی و هندوانه موجود در بانک ژن NCBI دو گروه اصلی و هفت زیر گروه را ایجاد کرد. بررسی ساختار دوم پروتئین آنالوگ ژن مقاومت به ToMV با استفاده از نرم افزار PSIpred نشان داد که این پروتئین فقط از مارپیچ α تشکیل شده است. ترشح این پروتئین در رقم طالبی TN-92-99 در غشاء پلاسمایی است و ترشح خارج سلولی آن ناچیز است. این اولین مطالعه در رابطه باکلون، بیان و تعیین فعالیت آنزیمی ژن مقاومت به ویروس موزاییک گوجه فرنگی در طالبی بومی ایران است. تظاهر این ژن میتواند نقش مهمی در مقاومت به بیماریهای ویروسی از جمله بیماری موازییک گوجه فرنگی در کدوئیان داشته باشد که از آن در برنامه های اصلاحی می توان استفاده نمود.
واژههایکلیدی: همسانهسازی، طالبی، مقاومت، ویروس موزاییک گوجه فرنگی، PCR
* نویسنده مسئول، تلفن: 09123708128 ، پست الکترونیکی: zamharir2005@yahoo.com
مقدمه
طالبی با نام علمی Cucumis melo var reticulatus گیاهی علفی از تیره کدوئیان میباشد .عقیده بر این است که طالبی از هندوستان منشأ گرفته است (1)،اما برخی منشأ آن را ایران و هند می دانند (14).
طالبی دارای واریته های زیادی است که از نظر شکل و طعم متفاوتند و شامل C. m. var. reticulatus و
C. m. var. indorus و dudaim C. m. var. میباشند (14). کدوئیان از جمله طالبی مانند سایر گیاهان مورد حمله آفات و بیماریهای مختلفی قرار میگیرند (10).
در حال حاضر تدبیر کارآمد برای مدیریت بیماریهای مختلف گیاهان از جمله کدوئیان استفاده از ارقام مقاوم است (8). در مورد بعضی بیماریها مانند بیماریهای فایتوپلاسمایی یافتن مقاومتهای طبیعی نادر است و بیشتر تمرکز بر استفاده از روشهای زیست فنآوری برای تولید گیاهان مقاوم میباشد (11 و 12). تودههای بومی گونههای کدوئیان و خویشاوندان وحشی آنها غنیترین ذخایر ژنتیکی در خصوص تحمل یا مقاومت به تنشهای زیستی هستند. مهندسی ژنتیک در دهه های اخیر، راهکارهای جدیدی برای معرفی مقاومت به بیمارگرها در گیاهان مختلف فراهم کرده است. گیاهان سازوکارهای مختلفی را برای دفاع در مقابل بیمارگرها به کار میبرند. این سازوکارها با ژنهای مقاومت گیاه (R gene) فعال میشوند (8).
با توجه به نقش ژنهای مقاومت بر علیه عوامل بیماریزای گیاهی و اهمیت هرمبندی ژنهای مقاومت در تولید ارقام مقاوم، محققان روشهایی جهت شناسایی و نشانمندکردن ژنهای مقاومت ابداع کرده اند. با استفاده از موتیفهای حفاظت شده ژنهای مقاومت، یک راهکار مبتنی بر واکنش زنجیرهای پایمراز (PCR) جهت جداسازی ژنهای مقاومت جدید و ایجاد نشانگرهای شدیداً پیوسته با ژنهای مقاومت در گونههای مختلف گیاهان زراعی به کار برده شده است (3و 6 ). PCR با استفاده از آغازگرهای الیگونوکلئوتیدی دجنره در جداسازی توالیهای حفاظت شده بسیار حساس بوده و از این رو در جداسازی ژنهای مقاومت کارآیی بالایی دارد. اگر ساختارهای حفاظت شده در بین تعداد زیادی از ژنهای مقاومت مشترک باشند، امکان زیادی برای جداسازی ژنهای مقاومت جدید مبتنی بر همولوژی توالی وجود دارد (4).
از آنجایی که اطلاعات کمی درمورد ژنهای مقاومت و آنالوگهای آنها و همچنین ساختار و عملکرد آنها در ارقام طالبی بومی ایران وجود دارد، هدف این تحقیق جداسازی و تعیین ویژگیهای توالی آنالوگ ژن مقاومت به بیماری در طالبی، ترسیم ساختارهای پروتئینی، پیشبینی عملکرد و بررسیهای فیلوژنتیکی با سایر ژنهای این خانواده با استفاده از ابزار های بیوانفورماتیکی بود.
مواد و روشها
تهیه بذور و کشت ارقام طالبی بومی و استخراج DNA: بذور ارقام بومی طالبی در ایران شامل TN-92-99، TN-92-80، TN-92-107، TN-92-120 و TN-92-115 از بانک ژن (موسسۀ تحقیقات اصلاح نهال و بذر، کرج) تهیه و در گلخانه موسسه تحقیقات گیاهپزشکی کشور کشت شدند. به منظور کشت بذور، گلدانها کاملاً شسته و با خاک استریل پر شدند و 7-5 بذر در هر گلدان قرار داده شد و سطح بذور با خاک نرم پوشانده شد. گلدانها در دمای 25-18 درجه قرار داده شدند. یکماه بعد از کشت، DNA از برگ گیاهان طالبی به روش CTAB با اندکی تغییر (5) استخراج شد.
گزینش آغازگر و انجام PCR: در این مطالعه از ترکیب آغازگری F1/R1 و F2/R1 که در مطالعات قبلی از روی موتیفهای حفاظت شده ژنهای کد کننده پروتئین NBS–LRR طراحی شده بودند برای واکنش PCR استفاده شد (جدول 1).
جدول1- توالی وترکیب آغازگرهای مورد استفاده در مطالعه آنالوگ ژنهای مقاومت در ارقام بومی طالبی ایران
Primer |
Conserved motif |
Primer sequence |
Reference |
F1 |
P-loop |
TGSSRGGHWYRGGBAAAACTAC |
Zhang et al. (2008) Zhang et al. (2008) |
R1 |
GLPL |
HRCWARAGGVARCCCTYBACA |
|
F2 |
P-loop |
GGDGTDGGNAARACWAC |
Deng et al. (2000) |
هر واکنش PCR در حجم نهایی 25 میکرو لیتری شامل 2-1 میکرو لیترDNA، ) dNTP10 میلی مولار) با حجم 1میکرولیتر، آغازگرها (10-20 میلی مولار) با حجم نهایی 2 میکرولیتر برای هر آغازگر، 5/2 میکرولیتر بافر(X10) ، MgCl2 با حجم نهایی1 میکرولیتر، آنزیم Taq پلیمرازU)1)2/0 میکرولیتر واحد (سیناژن ،ایران) و آب دو بار تقطیر شده استریل صورت گرفت. لولهها در دستگاه PCR قرار داده شده و 35 چرخه دمایی شامل 1دقیقه در 94 درجه سانتی گراد برای واسرشتگی DNA (اولین چرخه 5 دقیقه در 94 درجه سانتی گراد)، 1دقیقه در 55 درجه سانتی گراد برای اتصال آغازگرها،5/1 دقیقه در 72 درجه سانتی گراد برای ساختن DNA انجام گرفت. در آخر به منظور توسعه طول رشته DNA به مدت 5 دقیقه در 72 درجه باقی ماندند. محصول PCR در آگارز 2/1 درصد الکتروفورز شده و بعد از رنگ آمیزی با اتیدیوم برماید، باندها با اشعه UV مرئی شدند و از آنها عکسبرداری شد.
همسانه سازی در ناقل (Vector)pGEM-T: بعد از تکثیر ژن توسط PCR و خالص سازی آن به روش لئونارد (1998) (9)، عمل اتصال بین ژن و پلاسمید pGEM-T easy vector system I، محصول شرکت پرومگا (آمریکا) صورت گرفت. پلاسمیدpGEM-T شامل راهاندازهای T7 و SP6 در طرفین ناحیه برشی (MCS) است و ژن انتخابگر آن برای گزینش E. coli ژن مقاومت به آمپیسیلین میباشد. همچنین کاست ژنی حاوی β-گلوکورونیداز تحت راه اندازهای SP6 وT7 در طرفین سایت برشی اختصاصی MCS)) جهت گزینش به روش کلونی سفید و آبی وجود دارد. پس از الحاق قطعه به پلاسمید و ایجاد پلاسمید نوترکیب، انتقال به وسیله شوک حرارتی به باکتری E. coli نژاد Top10 انجام شد. در مرحله بعد باکتریها روی محیط LBجامد حاوی آنتیبیوتیک آمپیسیلین، IPTG و X-GAL کشت داده شدند. پس از 16 ساعت تیمار در دمای 37 درجه سانتی گراد، تعدادی کلونی سفید به عنوان کلونیهای نوترکیب انتخاب و در نهایت استخراج پلاسمیدهای نوترکیب با استفاده از کیت GF-1 Plasmid DNA Extraction Kit (Vivantis, Malasia) انجام گرفت. کلونیهای نوترکیب با استفاده از تکنیک PCR مورد بررسی قرار گرفتند. پلاسمیدها با دو آغازگر T7 و SP6 توالییابی شدند.
آنالیرهای تعیین توالی و رسم درخت فیلوژنی: با استفاده از بانک اطلاعاتی NCBI و BLAST، توالی نوکلئوتیدی و اسید آمینهای آنالوگهای مختلف مقاومت در ارقام طالبی بومی ایران و چندین گونه دیگر به همراه سایر ویژگیها (ویژگیهای توالیهای مورد بررسی از جمله طول mRNA ، توالی cDNA، تعداد اگزون و اینترون، طول پروتئین و شماره دسترسی توالیها) مورد بررسی قرار گرفتند (2). همردیفی مقایسهای ژن کلون شده در وکتور pGEM-T با سایر گیاهان موجود در NCBI با استفاده از نرم افزار T-Coffee به روش ClustalW انجام شد (16). بررسی فیلوژنتیکی با استفاده از آنالوگهای مختلف مقاومت در ارقام طالبی بومی ایران و انواع ژنهای NBS با استفاده از نرمافزار MEGA6 به روش UPGMA انجام شد (13). به منظور تأیید صحت و اعتبار درختهای حاصل، تست BootStrap با 1000 تکرار صورت گرفت (7).
تعیین مشخصات و شناسایی قلمروهای آنالوگهای مختلف مقاومت در ارقام طالبی بومی ایران: توالی پروتئینی با استفاده از پایگاههای اطلاعاتی پروتئینی و UniProtKB تعیین شد. تعیین قسمتهای درون، خارج سلولی و غشایی پروتئین با استفاده از برنامه TMHMM صورت گرفت.
تعیین ساختار و محل عملکرد پروتئین: برای تعیین ساختارهای اول و دوم پروتئین از برنامههای Uniprot B،PSIpredو برای تعیین محل عملکرد آن از بانک PSORT PROTEIN استفاده شد.
نتایج
نتایج تکثیر ژنهای آنالوگ مقاومت با استفاده از ترکیبهای آغازگری F1/R1 و F2/R1 نشان داد ترکیب آغازگری F2/R1 قادر به تکثیر آنالوگهای مقاومت به بیماری در ارقام طالبی TN-92-99 و TN-92-80 بومی ایران است (شکل 1) ولی ترکیب آغازگری F1/R1 این قابلیت را ندارد.
شکل 1- تکثیر آنالوگهای مختلف ژن مقاومت در ارقام طالبی بومی ایران: چاهک 5-1 به ترتیب ارقام طالبی TN-92-107، TN-92-120، TN-92-115، TN-92-99 و TN-92-80، چاهک Tom مربوط به گوجه فرنگی به عنوان کنترل منفی و چاهک 6 واکنش انجام شده با آب است و M نشانگر اندازهای Kb 1 (فرمنتاز).
آنالیز بلاست با استفاده از برنامه NCBI نشان داد که توالی NBS در رقم طالبی TN-92-99 ایرانی 92 درصد با ژن مقاومت به ویروس موزاییک گوجه فرنگی در طالبی و در رقم طالبی TN-92-80 ایرانی 99 درصد شباهت را با ژن مقاومت NBS-LRR به ترتیب در هندوانه دارد. آنالوگهای مقاومت جداسازی شده از سایر ارقام طالبی بومی ایران شباهتی با ژنهای مقاومت ثبت شده در بانک ژن نداشتند. در بعضی از ارقام مانند TN-92-107، TN-92-115 و TN-92-120 با هیچیک از ترکیبهای آغازگری مورد استفاده در این مطالعه، ژنی تکثیر نشد (شکل 1).
نتایج همردیفی با استفاده از بانک اطلاعاتی NCBI و نرم افزار T-Coffee به روش ClustalW نشان داد که این آنالوگهای مقاومت به بیماری در ارقام طالبی TN-92-99 و TN-92-80 بومی ایران، حفاظت شدگی بالایی به طول حدود 60 جفت باز را در دو رقم نشان میدهند (شکل 2).
شکل2- همردیفی مقایسهای آنالوگهای مقاومت کلون شده در ناقل pGEM-T با آنالوگهای مشابه در هندوانه و موجود در NCBI (JN238699, JN238688, JN238689, JN238672, JN238671) با استفاده از برنامه T-Coffe.
به منظور بررسی روابط تکاملی بین ارقام طالبی بومی ایران از نظر توالی آنالوگهای مقاومت به بیماری، رسم درخت فیلوژنی انجام شد، تا میزان شباهت آنها با سایر توالیهای همساخت در طالبی و هندوانه موجود در پایگاه بررسی شود. با توجه به اینکه ژنوم جانداران حجم گستردهای از اطلاعات را در خود دارد و بررسی تمام این توالیها بسیار دشوار است، معمولا از توالیهایی استفاده میشود که در موجودات همساخت باشند. ترسیم درخت فیلوژنی با استفاده از آنالوگهای جداسازی شده از خربزه و هندوانه موجود در بانک ژن NCBI دو گروه اصلی و هفت زیر گروه را ایجاد کرد (شکل 3). گروه اول شامل آنالوگهای JN230699 ،JN230690, JN230689 و JN230671است و گروه دوم شامل آنالوگهای مقاومت 80 و 99 در ارقام طالبی بومی ایران به همراه XM-008465356 و JN230688 است که بیشترین شباهت این آنالوگها در ارقام طالبی بومی ایران با ژن مذکور مربوط به گیاه خربزه است (XM-008465356). آنالوگهای مقاومت گروه اول به صورت جداگانه از گروه دوم قرار گرفته اند.
شکل 3- درخت فیلوژنتیکی توالی همسانه شده با سایر ژنهای همتای خود در بانک اطلاعاتی NCBI با استفاده از نرمافزار MEGA6
بررسی ساختار دوم پروتئین آنالوگ ژن مقاومت همسانه شده 99 و 80، به ترتیب جداسازی شده از ارقام طالبی بومی ایران TN-92-99 و TN-92-80، با استفاده از PSIpred نشان داد که این پروتئین فقط از مارپیچ α تشکیل شده است (شکل 4). نتایج بیوانفورماتیکی مربوط به هدفگیری پروتئین، با استفاده از نرم افزار PSORT PROTEINT، نشان داد که بیشترین میزان ترشح پروتئین درغشاء پلاسمایی بوده است و ترشح خارج سلولی آن ناچیز است.
شکل4- ساختار دوم پروتئین ژن آنالوگ مقاومت 80 (بالا) و 99 (پایین) جدا شده از ارقام بومی طالبی با استفاده ازPSIpred
بحث
ژنهای خانواده NBS-LRR نقش مهمی در پاسخ مقاومت گیاهان به بیمارگرهای مختلف دارند (15). با در نظر گرفتن خسارت بیماریهای غیرقابل کنترل ویروسی و فایتوپلاسمایی در کدوئیان، اهمیت شناسایی و آنالیز ژنهای مقاوم به بیمارگرهای مذکور در کشاورزی حائز اهمیت است. با توجه به این خسارتها در باغبانی، تولید گیاهان تراریخته متحمل از طریق انتقال ژن میتواند منجر به افزایش تولیدات، کاهش مصرف سم و رسیدن به کارآیی مطلوب تولید آن در واحد سطح شود. مطالعاتی که تا کنون روی ژنهای خانواده NBS-LRR ارقام تجاری کدوئیان انجام شده نشان می دهد این ژنها در این خانواده دارای نقاط محاظت شده هستند که بر اساس این توالیها پرایمرهای مختلفی طراحی شده است (15 و 16). در این مطالعه آغازگرهای مختلفی بررسی شد و نتایج آن ما نشان داد که آغازگرF2/R1 در تمامی ارقام طالبی بومی ایران قادر به تکثیر آنالوگ ژنهای مقاومت به بیماریها می باشد.
دادههای به دست آمده در مطالعه اخیر نشان داد، آنالوگهای مختف مقاومت در ارقام طالبی بومی ایران وجود دارد که تظاهر آنها میتواند نقش مهمی در مقاومت به بیماریهای ویروسی از جمله بیماری موازییک گوجه فرنگی در کدوئیان داشته باشد. زیرا خانواده کدوئیان در مقایسه با سایر گیاهان دارای تعداد کمی ژن مقاومت است (17). به عنوان مثال تعداد ژنهای متعلق به خانواده NBS-LRR در خربزه و هندوانه 8 عدد است که روی کروموزوم 9قرار دارند و در خیار 12 عدد است (17). این ژنها برای اولین بار از ارقام ایرانی طالبی TN-92-99 و TN-92-80 از طریق PCR جداسازی و در پلاسمید pGEM کلون شدند.
درخت فیلوژنی براساس توالی آنالوگهای مقاومت 90 و 88 جداسازی شده از ارقام طالبی بومی ایران و آنالوگهای جداسازی شده از طالبی و هندوانه موجود در بانک ژن NCBI به منظور بررسی روابط تکاملی نشان داد که آنالوگهای طالبی ایرانی قرابت زیادی با آنالوگهای XM-008465356 و JN230688 جداسازی شده از هندوانه و خربزه دارند. همچنین آنالوگهای جداسازی شده از خربزه و هندوانه نیز در درخت فیلوژنی در هر دو شاخه اصلی کنار هم قرار گرفته اند که بیانگر شباهت بالای این آنالوگها در سطح مولکولی بین گونههای خربزه و هندوانه است که در مطالعات دیگر نیز این نتایج حاصل شده است (16 و 17) و به عبارت دیگر مفهوم ژن آنالوگ مقاومت را بیان می کند و اینکه این ژنها نقش مهمی در مقاومت در گونه های مختلف گیاهی دارد. نکته حائز اهمیت در این درخت وجود آنالوگهای مختلف مقاومت در خربزه و هندوانه با عملکرد یکسان است که نشان میدهد این ژن میتواند نقش کلیدی در گیاه داشته باشد و همین باعث شده تا به صورت حفاظت شده در گونههای خربزه و هندوانه در طول تکامل حفظ شود.
مطالعه ساختار دوم پروتئین نشان داد، پروتئین 99 و 80، به ترتیب جداسازی شده از ارقام طالبی بومی ایران TN-92-99 و TN-92-80، فقط از مارپیچ α تشکیل شده است. پس دارای انعطاف و حلالیت بیشتر است. انعطافپذیری بیشتر باعث پایداری کمتر پروتئین میشود، در نتیجه میتوان گفت که این پروتئین دارای پایداری نمیباشد. و این مسئله در پروتئینهای مقاومتی که در زمان تنش و واکنش به یک پاتوژن بیان می شوند مشاهده می شود (17).
براساس آنالیزهای مختلف بیوانفورماتیکی آنالوگهای مختلف مقاومت، شباهت زیادی با توالی این ژن ها در سایر کدوئیان ثبت شده در NCBI دارند. از آنجایی که خانواده NBS-LRR ها در مقاومت عمودی نقش دارند، این مطالعه راهی برای تولید گیاهان تراریخته طالبی متحمل به بیماریهای با مدیریت دشوار، نظیر بیماریهای ویروسی و فایتوپلاسمایی کدوئیان از طریق انتقال ژن در آینده است.