کاربرد یادگیری ماشین در داده‌های پروتئومیک
کشف تنوع پروتئین‌ها، یک رویکرد مهم در تشخیص بیماری در علوم زیستی است. معیار فعلی برای بررسی و توضیح اطلاعات حاصل از متغیرهای بیولوژیکی متعدد، موضوعی است که در علوم زیست‌شناسی، امیکس (omics) نامیده می‌شود. با استفاده از تکنیک‌های داده کاوی چند متغیره در دو رویکرد اصلی، یادگیری ماشین و مبتنی بر آمار، این تنوع مورد بررسی قرار می‌گیرد. به‌طور معمول، مطالعات پروتئومیک، بسته به بستر تحلیلی یا روش استفاده شده برای تولید داده، می‌توانند صدها یا هزاران متغیر (p) در هر مشاهده (n) تولید کنند. بسیاری از روش‌های طبقه‌بندی به دلیل n≪p محدود می‌شوند و به همین ترتیب، پیش از طبقه‌بندی، نیاز به پیش‌پردازشی برای کاهش ابعاد دارند. یکی از محدودیت‌های چنین روش‌هایی، عدم وجود یک مدل است که امکان تفسیر معنی دار ویژگی‌های به کار رفته در طبقه‌بندی را فراهم می‌سازد. این مشکل، ممکن است با استفاده از یک رویکرد مبتنی بر مدل آماری حل شود که در آن، نه تنها اهمیت پروتئین مشخص است، بلکه می‌توانند، به راحتی با یک قانون طبقه‌بندی قابل تفسیر، ترکیب شوند. بدین خاطر، اخیراً تکنیک‌های یادگیری ماشین به دلیل توانایی طبقه‌بندی نمونه‌های ناشناخته، محبوبیت زیادی در این زمینه کسب کرده‌اند.
__فهرست__
در طی یک دهه گذشته، مطالعات پروتئینی با استفاده از روش‌های سنتی پیشرفته، از تجزیه و تحلیل ژن‌های واحد تا آزمایش پروتئین‌های متعدد، رشد چشمگیری داشته‌است.
پروتئومیکس، یک رشته تحقیقاتی است که هدف آن، ارزیابی دقیق و توصیف پروتئین‌های یک ارگانیسم، نمونه یا بافت است. فناوری‌های فعلی این امکان را به دانشمندان می‌دهد تا در یک بستر واحد، اطلاعات صدها تا هزاران پروتئین یا پپتید را به‌طور همزمان جمع‌آوری کنند. چنین روش‌هایی، به دلیل توان بالا و اندازه‌گیری همزمان تعداد زیادی متغیر تولید شده در یک آزمایش، در پروژه‌های کشف نشانگرهای زیستی استفاده می‌شوند. با وجود آن که ارزیابی پروتئومیک، خروجی تحقیقات در رشته‌های مختلف را بهبود بخشیده‌است، اما مشکلاتی در تفسیر و تجزیه تحلیل متغیرهای همزمان اندازه‌گیری شده، ایجاد کرده‌است. این مشکلات، مشابه مشکلاتی است که در تحقیقات مربوط به بیان ژن، دیده می‌شود. هنگام تصمیم‌گیری در مورد تکنیک تحلیلی مناسب، باید تعدادی از خصوصیات یک مجموعه داده پروتئومیک در نظر گرفته شوند. از مهم‌ترین خصوصیات که در اصطلاح نفرین ابعاد نامیده می‌شود، زمانی است که n≪p باشد. یعنی، تعداد مشاهدات (n) به مراتب از تعداد متغیرها (p) که همان پروتئین یا پپتیدها هستند، کوچکتر باشد. این مسئله، می‌تواند منجر به مشکلاتی شود که قابلیت تعمیم و در نتیجه کاربرد بالینی ابزارهای تشخیص را محدود می‌کند. زمانی‌که متغیرها بسیار همبسته هستند، این رشته، روش‌های تجزیه و تحلیلی را ارائه می‌دهد که واحدهای جرمی را نامناسب دانسته و جزئیاتی که در مطالعات قبلی نادیده گرفته می‌شده‌اند را بررسی می‌کند. همبستگی زیاد، در نتیجه تنظیم شدن پروتئین‌های خاصی است که می‌توانند در تنظیم بالا/پایین دیگر پروتئین‌ها، تأثیر داشته باشند. به عنوان مثال یک پروتئین، ممکن است بین مردها و زن‌ها متفاوت باشد اما به شدت، با وضعیت پاتولوژی یک بیماری همبستگی داشته باشد.
خواص داده‌های پروتئومیک که در بالا مورد بحث قرار گرفت، می‌تواند مانع ایجاد یک طبقه‌بندی مناسب، برای تمایز بین گروه‌های مختلف شود (به عنوان مثال، تمایز بین افراد بیمار و غیر بیمار). یکی از راه‌حل‌های این مشکل که اغلب، قبل از تحلیل چند متغیره استفاده می‌شد، کاهش ابعاد داده‌های خام است. متداول‌ترین روش فیلتر کردن است؛ به عنوان مثال می‌توان به انتخاب متغیرهای آماری مرتبط (پروتئین یا پپتیدها)، پیش از ایجاد مدل اشاره کرد. این روش امکان حذف متغیرهای اضافی و زائد را فراهم می‌سازد.
از طرف دیگر، روش رایج‌تر در تحقیقات مربوط به کشف نشانگرهای زیستی، حذف همه متغیرهای مجموعه داده و سپس به‌کارگیری رویکردهای کاهش ابعاد است. چنین تلاش‌هایی اجازه می‌دهد تا واحدهای جرمی که در طبقه‌بندی تأثیرگذار هستند، آشکار شده و بقیه در پس‌زمینه باقی بمانند. همچنین کاهش ابعاد، برای اثرات متغیرهای بسیار همبسته، به عنوان یک ویژگی کلیدی داده‌های پروتئومی و ژنومی به‌حساب می‌آید. علاوه بر این، تحلیلگر ممکن است از یک روش ترکیبی برای انتخاب متغیر و کاهش ابعاد استفاده کند تا بتواند مجموعه‌ای از نشانگرهای زیستی را ارائه دهد که منجر به یک طبقه‌بندی با نتایج مناسب شود. با این حال، این تصمیم اغلب وابسته به رویکرد انتخابی تحلیلگران برای طبقه‌بندی است.
در گذشته، از تکنیک‌های مختلفی برای تجزیه و تحلیل داده‌های پروتئومیک استفاده شده‌است. روش‌های محاسباتی از جمله ماشین‌های بردار پشتیبانی (SVM)، شبکه‌های عصبی مصنوعی (ANN) و جنگل‌های تصادفی (PLS-RF) و همچنین روش‌های مبتنی بر مدل، مانند حداقل مربعات جزئی و واکنش زنجیره‌ای پلیمراز در آنالیز افتراقی خطی (PLS-LDA و PCR-LDA) هستند.
Willingale و همکاران در سال ۲۰۰۶، از SVM , ANN، الگوریتم‌های ژنتیکی و درخت تصمیم‌گیری برای طبقه‌بندی داده‌های تولید شده از بیماران مبتلا به نارسایی قلبی استفاده کردند. آنها برای طبقه‌بندی خود، مدلی بر اساس یک مجموعه داده آموزشی متشکل از ۱۰۰ نمونه نارسایی قلبی و ۱۰۰ نمونه تحت کنترل ایجاد کرده و آن را با ۳۲ نمونه نارسایی قلبی و ۲۰ نمونه تحت کنترل آزمون کردند. مدل نهایی توانست ۸۸٫۵ درصد از داده‌های آزمون‌شده را به درستی طبقه‌بندی کند.
اسمیت و همکاران در سال ۲۰۰۷، از SVM برای طبقه‌بندی پاسخ فاز اول روش‌های درمانی جدید برای بیماران مبتلا به تومور رکتال استفاده کردند. حساسیت این مدل بین ۲۵ تا ۸۷٫۵ درصد بود. با این حال، مشکل اصلی این مطالعه، تعداد کم مشاهداتی بود (۲۰ نفر) که طبقه‌بندی بر اساس آن انجام شد.
Purohit و Rocke در سال ۲۰۰۳ از تکنیک‌های طبقه‌بندی نظارت شده و بدون نظارت و همچنین کاهش ابعاد به روش PCA استفاده کردند که در ادامه منجر به تجزیه و تحلیل خوشه ای برای طبقه‌بندی داده‌های پروتئومیک بین نمونه‌های سالم و بیمار شد. علاوه بر این، آنها با استفاده از روش‌های رگرسیون لجستیک و آنالیز افتراقی، ترکیبی از PLS و PCR برای طبقه‌بندی ارائه دادند.
لیو و همکاران در سال ۲۰۰۸ از روشهای مبتنی بر PLS برای طبقه‌بندی داده‌های مربوط به سرطان تخمدان استفاده کردند. آنها PLS در آنالیز افتراقی خطی، PLS در کی نزدیکترین همسایه، PLS در رگرسیون لجستیک و PLS در شبکه‌های عصبی مصنوعی را با طیف گسترده‌ای از روشهای طبقه‌بندی مبتنی بر PCA مقایسه کردند. یافته‌های آنها حاکی از آن بود که کاهش بعد PLS حاصل از طبقه‌بندی رگرسیون لجستیک، نتایج بهتری را از روشهای مبتنی بر PCA و سایر رویکردهای PLS ایجاد می‌کند.
رجالتی و همکاران در سال ۲۰۰۹ از PLS استفاده کردند تا با کاهش ابعاد حاصل از آنالیز افتراقی، نمونه‌های مایع مغزی نخاعی (CSF) و نمونه‌های CSF جابه‌جاشده با استانداردهای پپتید را طبقه‌بندی کند.
یک ماتریس داده پروتئومیک (X) شامل متغیرهای پاسخ به شکل پروتئین یا پپتید است. هر سطر آن نشان‌دهنده یک مشاهده یا یک نمونه بوده و هر ستون آن مربوط به یک پروتئین، پپتید یا نسبت جرم به وزن است. توجه کنید که اصطلاح واحد جرمی، برای نشان‌دادن پروتئین‌ها، پپتیدها یا نسبت جرم به وزن استفاده می‌شود. در مسائل طبقه‌بندی، یک بردار از متغیرهای خام (y) برای شناسایی کلاس نمونه مشاهده‌شده، کدگذاری شده‌است. در موارد چند کلاسی، به جای بردار y، از ماتریس Y استفاده می‌شود. در ادامه برخی از روش‌های مختلف یادگیری ماشین برای طبقه‌بندی این داده‌ها شرح داده می‌شود:
formula_1 متمرکز سازی داده‌ها، شامل جابه‌جایی داده‌های اصلی است به گونه ای که بردار میانگین کل، برابر با صفر باشد. این در حالی است که کاهش داده‌ها به واریانس واحد، باعث می‌شود همه متغیرها، به‌طور مساوی و صرف نظر از واریانس هریک، در نحوه ارائه مشاهدات در فضای ابعادی کاهش یافته، سهیم باشند. این ویژگی، در صورتی مفید است که بزرگی واریانس هر واحد جرمی به اهمیت نسبی، آن مربوط نباشد؛ بنابراین اگر فقط داده‌های متمرکز (ماتریس کواریانس) وارد الگوریتم PCA شوند، نسبت جرم به وزن‌های هر داده، تأثیر بیشتری در چگونگی مشاهده‌ها در فضای ابعادی کمتر خواهد داشت. در حالی که متمرکز سازی و مقیاس‌بندی (همبستگی)، اثرات ناشی از نسبت جرم به وزن با واریانس‌های بزرگ را کاهش می‌دهد. اگر همه متغیرها، واحد و نوع یکسان داشته باشند، از ماتریس کواریانس استفاده می‌شود. به‌طور خلاصه، PCA سعی در ایجاد ترکیبات خطی در متغیرهای اصلی که به‌طور خطی مستقل (متعامد) از یکدیگر هستند، دارد. این کار به گونه‌ای انجام می‌شود که فاصله اقلیدسی در بین مشاهدات حفظ شود.
formula_14