آنالیز شبکه هم‌بیانی ترانسکریپتوم لوبیا به منظور شناسایی ماژول‌ها و ژن‌های هاب درگیر در مقاومت به کنه تارتن دولکه‌ای

کنه تارتن دولکه‌ای (tetranychus urticae koch) از مهم‌ترین آفات تخریب‌کننده لوبیا (phaseolus vulgaris l.) است. به‌طور کلی شبکه‌های ژنی پیچیده‌ای در ایجاد حساسیت یا مقاومت در مقابل کنه تارتن دولکه‌ای دخیل هستند، بنابراین در این تحقیق از روش سامانه‌های زیستی به‌کار برده شد. برای این منظور از داده‌های rna-seq مربوط به تنش کنه تارتن روی گیاه لوبیا استفاده شد. پس از فراهم‌کردن ماتریس بیان ژن‌ها، شبکه‌های مولکولی با استفاده از آنالیز شبکه هم‌بیان وزن‌دار (wgcna) مورد تجزیه و تحلیل قرار گرفتند. پس از ایجاد ماژول‌ها، عملکرد ژن‌ها در هر ماژول مورد بررسی و تجزیه و تحلیل قرار گرفت. بر اساس نتایج، مجموع 699 ژن با بیان افتراقی در پاسخ به تنش کنه تارتن شناسایی شدند که در 7 ماژول هم‌بیان از طریق خوشه‌بندی سلسله‌مراتبی جای گرفتند. بررسی هستی‌شناسی ژن‌ها و آنالیز برهمکنش ژن‌های کلیدی با استفاده از بانک اطلاعاتی string نشان داد پاسخ ترانسکریپتوم لوبیا به آلودگی با کنه تارتن بیشتر شامل ژن‌های کدکننده پروتئین کینازها، کاتالیزورها، فاکتورهای رونویسی، ساخت متابولیت‌ها و مسیرهای انتقال پیام هورمونی بود. ماژول فیروز‌ه‌ای بیشترین ژن‌های درگیر در مقاومت را دارا بود که این ماژول و ماژول زرد، بیشترین همبستگی را با رقم مقاوم به ترتیب پس از پنج و یک روز آلودگی داشتند. همچنین، ماژول مشکی بیشترین همبستگی با رقم حساس پس از پنج روز آلودگی را داشت. این مطالعه دانش ما را از مکانیسم‌های مولکولی دخیل در مقاومت به کنه تارتن افزایش می‌دهد. همچنین ژن‌های بررسی شده در این تحقیق می‌توانند به‌عنوان اهداف اصلاحی برای ایجاد مقاومت معرفی شوند.

co-expression network analysis of common bean transcriptome in order to identify modules and hub genes involved in resistance to tetranychus urticae (acari, tetranychidae)

two-spotted spider mite (tetranychus urticae koch) is one of the most important pests of beans (phaseolus vulgaris l.). since, complex gene networks are involved in creating sensitivity or resistance against the two-spotted spider mite; therefore, in this research we used biological system methods to identify key networks. for this purpose, we used the rna-seq data related to the two-spotted spider mite stress on common bean plant. after providing the gene expression matrix, molecular networks were analyzed using weighted co-expression network analysis (wgcna). after the modules identification, the gene functions in each module were investigated and analyzed. according to the results, a total of 699 genes were identified with differential expression in response to two-spotted spider mite stress, which were placed in 7 co-expression modules through hierarchical clustering. gene ontology and interaction analysis of key genes using the string database showed that the response of common bean transcriptome to two-spotted spider mite infestation includes genes encoding protein kinases, catalysts, transcription factors, and metabolite production and pathways of hormonal message transmission. it is notable that among the most important genes that showed co-expression, wrky and lipoxygenase were highlighted. the turquoise module had the higher number of genes involved in resistance, and this module and the yellow module had the highest correlation with the resistant variety after five and one day of contamination, respectively. also, the black module had the highest correlation with the sensitive variety after five days of contamination. in conclusion, this study increases our knowledge of the molecular mechanisms involved in resistance to the two-spotted spider mite. also, the genes examined in this research can be introduced as breeding targets to create resistance.