فراتحلیل دادههای بیانی RNA-Seq و Microarray برای شناسایی ژنهای مؤثر در رشد و نمو عضله گوسفند
نویسندگان
1 گروه علوم دامی، دانشکده کشاورزی، دانشگاه فردوسی مشهد، مشهد، ایران.
2 گروه علوم دامی دانشکده کشاورزی، دانشگاه فردوسی مشهد، مشهد، ایران.
3 گروه علوم دامی، دانشکده کشاورزی، دانشگاه فردوسی مشهد، مشهد، ایران.
doi
10.22067/ijasr.2023.79909.1111چکیده
امروزه میتوان با روشهای بیوانفورماتیکی دادههای حاصل از مطالعات و پلتفرمهای مختلف را ادغام و از آنها بهره برد. در این مطالعه، با ادغام دادههای ریزآرایه و RNA-Seq بافت عضله گوسفند نژاد تکسل موجود در پایگاه داده به مقایسه پروفایل ترنسکریپتومی عضله در دو مقطع سنی جنینی و بلوغ پرداخته شد. برای محاسبه مقادیر بیانی عضله ریزآرایه مربوط به دوران جنینی از بستههای نرمافزاری Limma، Biobase و GEOquery در محیط R و برای محاسبه مقادیر بیانی عضله RNA-Seq از پروتکل Tuxedo و بسته نرمافزاری HTSeq در محیط لینوکس و بستهنرمافزاری DESeq2 در محیط R استفاده شد. سپس دو نوع مقادیر بیانی ادغام شدند. نتایج نشان داد، در بافت عضله بین مقطع سنی بلوغ و جنینی بیان 62 ژن (37 ژن افزایش و 25 ژن کاهش بیان) اختلاف معناداری داشتند. با رسم شبکه ژنی بین ژنهای افتراقی، 15 ژن منتخب MYH1، ACTN3، CASQ1، TMOD4، FBP2، SLC2A4، MX1، COX4I1، SOD2، MFN2،UQCRB، UCP3، PRKAB2، PHKG2، PPP1R3C شناسایی شدند. عملکرد این ژنها در تکثیر سلولی، تشکیل میوفیبریلها و متابولیسمهای چربیزایی ثابت شدهاست. آنالیز هستیشناسی ژنهای افتراقی نقش برخی از این ژنها مثل ACTN3 وCASQ1 را در فرآیندهای زیستی مثل توسعه سلول عضلانی مخطط، مسیرهای علامتدهی Calcineurin-NFAT و JAK-STAT آشکار کرد. این مطالعه علاوهبر تأیید روش ادغامی دادههای ناهمگن، دیدی کلی از تفاوتهای ترنسکریپتومی بافت عضله گوسفند تکسل در دو مقطع مهم سنی را فراهم آورد تا ژنهای منتخب معرفیشده منبع مفیدی برای بررسیهای زیستی ژنهای مربوط به رشد و نمو عضله باشد.