استنباط بیزی تبارزایش

استنباط بیزی فیلوژنی از تابع درست نمایی برای ایجاد کمیتی به نام احتمال پسین درخت‌ها (درخت‌های فیلوژنی) استفاده می‌کند. محتمل‌ترین درخت فیلوژنتیک برای داده‌های داده شده با استفاده از مدل فرگشت و بر اساس احتمالات پیشین تولید خواهد شد. دلیل محبوبیت رویکرد بیزی، پیشرفت سرعت محاسبات و یکپارچه الگوریتم‌های زنجیره مارکوف مونت کارلو (MCMC) می‌باشد. استنباط بیزی در فیلوژنتیک مولکولی و سیستماتیککاربرد دارد.

ریشه و اساس استنباط بیزی فیلوژنی

استنباط بیزی روش احتمالی توسعه یافته مبتنی بر قضیه بیز می‌باشد که توسط توماس بیز ارائه شده‌است. استنباط بیزی در سال ۱۷۶۳ پس از مرگ بیز منتشر شده‌است. این اولین بیان احتمال معکوس و اساس استنباط بیزی می‌باشد. پی یر سیمون لاپلاس به‌طور مستقل و بدون آگاهی از کار بیز، قضیه بیز را در سال ۱۷۷۴ توسعه داده است.
استنباط بیزی تا سال ۱۹۰۰ به دلیل محدودیت محاسبات به‌طور گسترده‌ای استفاده می‌شد. در سال ۱۹۰۰ شاهد انتقال از استنباط بیزی به استنباط فراوانی‌گرایانه بودیم. بر اساس قضیه بیز رویکرد بیزی احتمال پیشین درخت P(A) را با احتمال درست نمایی داده P(B) ترکیب می‌کند تا توزیع احتمال پسین P(A|B) را روی درخت تولید کند. احتمال پسین درخت، احتمال درست بودن درخت را نشان می‌دهد. درخت با بالاترین احتمال پسین، بهترین گزینه برای نمایش فیلوژنی است. استنباط بیز مقدمه‌ای برای روش‌های زنجیره مارکوف مونت کارلو (MCMC) بود که توسط نیکولاس متروپلیس در سال ۱۹۵۳ عرضه شد. این روش استنباط بیزی منقلب شده‌است و phylogeneticistها تا دهه 1990s به‌طور گسترده‌ای از این روش استفاده کردند. از مزایای ای روش نسبت به روش‌های مرسوم بیشینه صرفه جویی و بیشینه درست نمایی می‌توان به امکان عدم قطعیت فیلوژنتیک، استفاده از اطلاعات پیشین و اختلاط مدل‌های پیچیده تکامل اشاره کرد که برای روش‌های مرسوم محدودیت تحلیل‌های محاسباتی دارند. با وجود غلبه بر پیچیدگی عملیات تحلیلی احتمال پسین، هنوز شامل پیچیدگی جمع روی تمام درختان است و برای هر درخت ادغام تمام ترکیبات جایگزینی مقادیر پارامتر مدل و طول انشعاب‌ها.
روش MCMC را می‌توان در سه مرحله توصیف کرد: اول با استفاده از مکانیزم تصادفی حالت جدید برای زنجیره مارکوف پیشنهاد می‌شود. دوم، احتمال درست بودن حالت جدید محاسبه می‌شود. سوم، یک متغیر تصادفی جدید در بازه (۰٬۱) پیشنهاد می‌شود. اگر این مقدار جدید کمتر از احتمال پذیرش باشد، حالت جدید پذیرفته شده‌است و حالت زنجیره آپدیت می‌شود. این فرایند هزاران یا میلیون‌ها بار اجرا می‌شود. مقدار زمانی که یک درخت مجزا در طول زنجیره‌ای مشاهده می‌شود تقریب معتبری از احتمال پسین آن می‌باشد. تعدادی از رایج‌ترین الگوریتم‌های مورد استفاده در روش‌های MCMC شامل الگوریتم متروپلیس‌هستینگز، MCMCMC (MC³) و الگوریتم محلی لارگت و سیمون می‌باشد.

بیشینه صرفه جویی و حداکثر درست نمایی

روش‌های زیادی برای بازسازی درختان فیلوژنتیک وجود دارند، هر کدام مزایا و معایب خود را دارد و پاسخ دقیقی برای این سؤال وجود ندارد: "بهترین روش کدام است؟". بیشینه صرفه جویی(MP) و حداکثر درست نمایی (ML) در روش‌های مرسومی هستند که برای تخمین فیلوژنی مورد استفاده قرار می‌گیرند. هر دو روش همانند روش‌های بیزی، از اطلاعات کاراکترها هم به‌طور مستقیم استفاده می‌کنند.
بیشینه صرفه‌جویی یکی (یا بیشتر) درخت بهینه را براساس ماتریس فاصله کاراکترها برای گروه خاصی از تاکسون و این نیازمند مدل تکاملی تغییر نیست. MP ساده‌ترین توضیح برای یک مجموعه داده می‌دهد، بازسازی یک درخت فیلوژنتیک که تا حد ممکن تغییرات کمی را در سراسر توالی داشته باشد، و این درخت کمترین تعداد مراحل تکاملی برای توضیح رابطه بین تاکسون‌ها را نشان می‌دهد. حمایت از شاخه‌های درخت توسط درصد بوت استرپ نمایش داده می‌شود. به همین دلیل سادگی آن است که به‌طور گسترده‌ای از آن استفاده شده‌است، انتقاداتی بر MP وارد شده‌است و توسط روش‌های حداکثر درست نمایی و بیزی به پس زمینه رانده شده‌است. MP ارائه مشکلات و محدودیت‌هایی را ارائه می‌کند. Joseph Felsenstein در سال (۱۹۷۸) نشان داد که MP ممکن است از نظر آماری ناسازگار باشد به این معنی که با داشتن داده‌های بیشتر و بیشتر (به عنوان مثال طول توالی) انباشته شده باشد، ممکن است نتایج به درخت نادرست همگرا شود و منجر به جذب شاخه طولانی - یک پدیده فیلوژنتیکی است که در آن تاکسون‌های با شاخه‌های بلند (تغییرات حالت کاراکتر متعدد)، مایل به این هستند که بیشتر از آن چیزی که در واقعیت هستند، در فیلوژنی ارتباط نزدیکتر داشته باشند - شود.
مثل بیشینه صرفه جویی، حداکثر درست نمایی، درختان دیگری را ارزیابی می‌کند. ضمناً احتمال اینکه هر درخت داده‌های داده شده را براساس مدل تکامل توصیف کند، در نظر می‌گیرد. در این روش درخت با بالاترین احتمال توصیف داده‌ها، انتخاب خواهد شد. به عبارت دیگر مقایسه‌ای بین درختان مختلف انجام می‌دهد، بر اساس اینکه چقدر هر درخت داده‌های مشاهده شده را پیش‌بینی می‌کنند. معرفی یک مدل تکامل در تحلیل ML، مزیتی نسبت به روش MP ارائه می‌کند: احتمال تعویض نوکلئوتید و نرخ این تعویض‌ها در نظر می‌گیرد، روابط فیلوژنتیک تاکسون را در قالب واقعی تر نشان می‌دهد. یک نکته مهم این روش، طول شاخه می‌باشد که تیغ اوکام آن را نادیده می‌گیرد، تغییرات در طول شاخه‌های بلندتر نسبت به رشته‌های کوتاهتر با احتمال بیشتر رخ می‌دهد. این رویکرد ممکن است جذب شاخه طولانی را نادیده بگیرد و سازگاری بیشتری از ML نسبت به MP توضیح دهد. اگر چه از دیدگاه تئوری توسط بسیاری به عنوان بهترین روش برای استنباط فیلوژنی در نظر گرفته شده، ML محاسبات سنگینی دارد و اکتشاف تمام درختان (تعداد بسیار زیاد درختان) تقریباً غیرممکن است. استنباط بیزی همچنین شامل یک مدل تکاملی است و از مزایای اصلی آن بر روش‌های MP و ML این است که از نظر محاسباتی کارآمد تر از روش‌های مرسوم می‌باشد، ضمناً می‌تواند منبع عدم قطعیت را تعیین نماید و قادر به این است که مدل‌های پیچیده تکامل را نیز لحاظ نماید.

نرم‌افزار MRBAYES

MrBayes یک نرم‌افزار رایگان است که استنباط بیزی در فیلوژنی را انجام می‌دهد. این نرم‌افزار توسط John P. Huelsenbeck و Frederik Ronquist در سال ۲۰۰۱ نوشته شده‌است. با افزایش محبوبیت روش‌های بیزی، MrBayes نرم‌افزار منتخب بسیاری از متخصصین فیلوژنتیک مولکولی شد. این نرم‌افزار برای سیستم عامل‌های مکینتاش، ویندوز و یونیکس عرضه شده‌است و واسط خط فرمان دارد. این برنامه از الگوریتم استاندارد MCMC و الگوریتم MCMCMC استفاده می‌کند. MrBayes ماتریس‌های هم تراز یابی توالی (DNA یا اسیدهای آمینه) را از NEXUS format استاندارد می‌خواند.
MrBayes با استفاده از MCMC احتمالات پسین درختان را تقریب می‌زند. کاربر می‌تواند مفروضات مدل جایگزینی، احتمالات پیشین و جزئیات تحلیل MCMCMV را تغییر بدهد. همچنین اجازه حذف و اضافه کردن تاکسون و کاراکترها را به کاربر می‌دهد. این برنامه از استانداردترین مدل جایگزینی DNA، 4x4 یا مدل تعویضی استقاده می‌کند. با این فرض که تغییرات در سراسر نوکلئوتید با احتمال برابر رخ می‌دهد. همچنین تعداد 20x20 مدل جایگزینی اسید آمینه و مدل‌های کدون جایگزینی DNA را پیاده‌سازی کرده‌است. این نرم‌افزار روش‌های مختلف برای ساده‌سازی فرض نرخ جایگزینی یکسان در سراسر مکان‌های نوکلئوتید ارائه می‌دهد. MrBayes همچنین می‌تواند تطبیق عدم اطمینان حالات اجدادی به درخت فیلوژنتیک و پارامترهای مدل را استنباط کند.
MrBayes 3 نسخه کاملاً بازسازی شده تجدید سازمان یافته نسخه اصلی MrBayes می‌باشد. جدیدترین مورد این نسخه، توانایی این نرم‌افزار به تطبیق ناهمگونی از مجموعه داده است. این چارچوب جدید به کاربر اجازه می‌دهد تا مدل‌ها را ترکیب کند و از مزایای کارایی الگوریتم تحلیل بیزی MCMC در هنگام برخورد با نوع داده‌های مختلف استفاده کند (به عنوان مثال پروتئین، نوکلئوتید و ریخت‌شناسی (زیست‌شناسی)). این نرم‌افزار به صورت پیش فرض از الگوریتم MCMCMC استفاده می‌کند.
MrBayes 3.2 نسخه جدید MrBayes منتشر شده در سال ۲۰۱۲ می‌باشد. این نسخه جدید به کاربران اجازه می‌دهد تحلیل‌های چندگانه را به صورت موازی اجرا کنند. همچنین محاسبات سریعتر درست نمایی را فراهم کرده‌است و می‌توان این محاسبات را به GPU محول کرد. نسخه ۳٫۲ خروجی گسترده‌تری سازگار با FigTree و دیگر نرم‌افزارهای مشاهده درخت‌ها فراهم کرده‌است.