Analisis Filogenetik dengan MEGA11: Dari Alignment hingga Tree Building

• ⏱️ 14 minit bacaan • 🌳 Bioinformatik • Filogenetik

Pokok filogenetik ialah “peta jalan evolusi” — visualisasi yang menunjukkan bagaimana organisma, gen, atau protein berkait antara satu sama lain sepanjang masa. Daripada pengenalan patogen baharu hingga penelusuran asal-usul spesies endemik, pokok ini menjadi tulang belakang hampir setiap penerbitan biologi moden. Dan untuk membuatnya, MEGA11 (Molecular Evolutionary Genetics Analysis versi 11) ialah perisian paling popular dalam kalangan ahli biologi — percuma, GUI mesra, dan cukup berkuasa untuk analisis tahap jurnal Q1. Masalahnya, ramai pelajar terperangkap: keliru memilih antara Neighbor-Joining vs Maximum Likelihood, salah baca bootstrap value, atau tidak tahu model penggantian mana yang tepat. Artikel ini membuka semua rahsianya dari sifar. Dan jika anda ada dataset jujukan besar dan memerlukan pokok yang sedia untuk penerbitan, pasukan Corvexum sedia melaksanakannya secara profesional.

Mengapa MEGA11 Menjadi Standard De Facto?

MEGA dibangunkan sejak 1993 oleh Sudhir Kumar dan pasukan (kini di Temple University). Versi 11 (keluaran 2021) membawa peningkatan besar:

  • GUI moden dengan mod gelap dan paparan lebih intuitif
  • Aligner bersepadu: MUSCLE, ClustalW, MAFFT terus dalam perisian
  • Pemilihan model automatik: ciri “Find Best DNA/Protein Models” yang memilih model penggantian terbaik berdasarkan AIC/BIC
  • 4 kaedah pembinaan pokok: NJ, MP, ML, dan UPGMA
  • Bootstrap & jackknife terbina dalam untuk ujian ketahanan topologi
  • Wizard TimeTree untuk pokok dengan skala masa (molekuler clock)

Tidak seperti perisian command-line (RAxML, IQ-TREE) yang memerlukan skrip, MEGA11 boleh dikuasai sesiapa sahaja dalam beberapa jam.

Langkah 1: Pemasangan & Penyediaan Data

Muat turun MEGA11 dari mega.com — percuma untuk akademik. Pasang dalam Windows/Mac/Linux. Sebelum mula, sediakan:

  • Jujukan pertanyaan: jujukan sasaran anda (format FASTA)
  • Jujukan homolog: 10-100 jujukan homolog dari NCBI/UniProt/Ensembl
  • Outgroup: sekurang-kurangnya satu jujukan dari spesies lebih jauh — penting untuk mengakar pokok

Tip praktikal: Ambil homolog melalui BLAST (lihat artikel sebelumnya), kemudian muat turun FASTA pukal menggunakan NCBI Batch Entrez atau E-utilities. Sasarkan 20-50 jujukan untuk pokok yang bermaklumat tetapi tidak overwhelming.

Langkah 2: Multiple Sequence Alignment

Buka MEGA11 → Align → Edit/Build Alignment → pilih Create a new alignment → DNA atau Protein. Import fail FASTA anda. Kemudian pilih aligner:

  • MUSCLE (default): pantas dan tepat untuk kebanyakan kes
  • ClustalW: klasik, baik untuk dataset kecil
  • MAFFT: terbaik untuk jujukan panjang atau bercapah

Selepas alignment selesai, periksa secara manual — alignment automatik tidak sempurna. Buang kawasan yang sangat bercapah atau penuh gap menggunakan ciri Select → Delete. Kawasan yang tidak ter-align dengan baik akan menghasilkan pokok yang mengelirukan.

Langkah 3: Pemilihan Model — Langkah Paling Kritikal

Model penggantian menentukan bagaimana bes/asid amino diandaikan berevolusi. Pemilihan model yang salah = pokok yang salah, walaupun kaedah anda betul.

Dalam MEGA11: Models → Find Best DNA/Protein Models. Perisian akan menilai 24 model (DNA) atau lebih (protein) berdasarkan kriteria:

  • AIC (Akaike Information Criterion) — paling biasa digunakan
  • BIC (Bayesian Information Criterion) — lebih konservatif
  • Log-likelihood — ukuran fit terhadap data

Model popular: T92+G+I, GTR+G+I untuk DNA; WAG+G+I, JTT+G+I untuk protein. Notasi “+G” = taburan Gamma (rate heterogeneity), “+I” = perkadaran tapak tidak berubah.

Langkah 4: Pembinaan Pokok — 4 Kaedah Utama

Kaedah Kelajuan Ketepatan Kes Penggunaan
NJ (Neighbor-Joining) Sangat pantas Sederhana Penerokaan awal, dataset besar
MP (Maximum Parsimony) Sederhana Sederhana Data morfologi, jujukan pendek
ML (Maximum Likelihood) Perlahan Sangat tinggi Penerbitan Q1, standard emas
UPGMA Sangat pantas Rendah Hanya untuk clustering kasar

Cadangan: Untuk penerbitan, sentiasa gunakan Maximum Likelihood dengan model terbaik dari langkah 3. NJ boleh jadi semakan sanity awal.

Langkah 5: Analisis Bootstrap — Ujian Kekuatan Pokok

Bootstrap menguji sejauh mana teguhnya topologi pokok anda. Pilih Phylogeny → Construct/Test Maximum Likelihood Tree → tetapkan Bootstrap method dengan 1000 ulangan (standard penerbitan).

Peraturan tafsiran:

  • Bootstrap ≥ 70%: nod dianggap kukuh / well-supported
  • Bootstrap 50-69%: nod lemah, perlu ditafsir dengan berhati-hati
  • Bootstrap < 50%: nod tidak boleh dipercayai, biasanya dibuang dari paparan

Dalam MEGA11, nilai bootstrap muncul sebagai nombor di setiap nod pokok. Untuk paper Q1, laporkan hanya nod dengan bootstrap ≥ 70%.

Langkah 6: Tafsiran Pokok

Pokok yang dihasilkan MEGA11 boleh dipaparkan dalam beberapa gaya:

  • Rectangular: klasik, baik untuk perbandingan panjang cabang
  • Radial/Circular: padat untuk pokok besar
  • Curved: estetik untuk pembentangan

Perkara yang perlu dibaca:

  • Topologi: siapa berkait rapat dengan siapa
  • Panjang cabang: berkadar dengan bilangan penggantian evolusi
  • Clade: kumpulan monofiletik (semua keturunan dari satu nenek moyang)
  • Root: akar pokok (ditentukan oleh outgroup)

Amaran: jika outgroup muncul di tengah clade ingroup, alignment atau rooting anda bermasalah — jangan teruskan tafsiran.

Langkah 7: Eksport & Penerbitan

MEGA11 mengeksport pokok dalam format:

  • PNG/TIFF: terus untuk rajah
  • SVG/PDF: vektor, boleh disunting dalam Adobe Illustrator
  • Newick: format standard untuk diimport ke FigTree atau iTOL

Untuk penerbitan, sentiasa gunakan SVG atau PDF dan sunting dalam Illustrator/Inkscape untuk menambah:

  • Scale bar (0.05 substitutions/site adalah biasa)
  • Label nilai bootstrap di setiap nod penting
  • Sorot clade yang ingin ditekankan
  • Pengkodan warna taksonomi untuk kejelasan visual
💡 Tip Penerbitan Pro: Buat pokok dengan dua kaedah berbeza (misalnya ML + NJ atau ML + Bayesian melalui MrBayes). Jika topologi terasnya konsisten, anda mempunyai bukti kukuh bahawa kesimpulan evolusi anda bukan artifak metodologi — ini sangat disukai reviewer Q1.

❓ FAQ

S: MEGA11 percuma?
A: 100% percuma untuk penggunaan akademik. Versi komersial tersedia untuk industri.

S: Berapa bilangan ulangan bootstrap yang cukup?
A: Minimum 1000 untuk penerbitan. Beberapa jurnal Q1 meminta 10000 untuk tuntutan kukuh.

S: MEGA11 boleh untuk genom penuh?
A: Kurang ideal. Untuk filogeni seluruh genom, pertimbangkan OrthoFinder + IQ-TREE. MEGA11 optimum untuk puluhan hingga ratusan jujukan.

🚀 Dataset Jujukan Banyak Tetapi Pokok Tidak Meyakinkan?

MEGA11 mudah dijalankan, tetapi memilih model, mentafsir topologi, dan menulis naratif evolusi yang betul memerlukan pengalaman bertahun-tahun. Salah pilih model, salah baca bootstrap, atau salah menentukan outgroup boleh menyebabkan kesimpulan evolusi anda ditolak reviewer Q1. Corvexum mengambil alih beban analisis filogenetik anda:

  • ✅ Kurasi jujukan homolog + alignment berkualiti (pemeriksaan manual)
  • ✅ Pemilihan model automatik dengan AIC/BIC
  • ✅ Pokok ML + NJ + bootstrap 1000 ulangan
  • ✅ Visualisasi figure-ready (SVG, PNG resolusi tinggi) + anotasi profesional
  • ✅ Naratif evolusi sedia tampal di bab hasil + konsultasi jurnal sasaran

Emel alternatif: corvexum@gmail.com


Tinggalkan Balasan

Alamat email Anda tidak akan dipublikasikan. Ruas yang wajib ditandai *

/* ============================================ CORVEXUM UNIVERSAL POST JS v2.2 ES5 + TANPA KARAKTER KHUSUS ( kurang-dari dan ampersand ) CTA Tracking + Smooth Scroll + PRISMA Animation ============================================ */ (function () { 'use strict'; function onReady(fn) { if (document.readyState !== 'loading') { fn(); } else { document.addEventListener('DOMContentLoaded', fn); } } onReady(function () { /* 1. CTA Click Tracking (siap untuk GA4/GTM) */ var btns = document.querySelectorAll('.cvx-btn'); var i; for (i = 0; i !== btns.length; i = i + 1) { (function (btn) { btn.addEventListener('click', function () { if (window.console) { if (window.console.log) { window.console.log('CTA Clicked: ' + btn.href); } } }); }(btns[i])); } /* 2. Smooth Scroll untuk anchor internal */ var anchors = document.querySelectorAll('a[href^="#"]'); var j; for (j = 0; j !== anchors.length; j = j + 1) { (function (link) { link.addEventListener('click', function (e) { var href = link.getAttribute('href'); var target = null; if (href) { if (href.length !== 1) { target = document.querySelector(href); } } if (target) { e.preventDefault(); target.scrollIntoView({ behavior: 'smooth', block: 'start' }); } }); }(anchors[j])); } /* 3. Animasi PRISMA Flow Diagram saat scroll into view */ var flows = document.querySelectorAll('.cvx-prisma-flow'); var m; var n; if (flows.length !== 0) { if (typeof window.IntersectionObserver === 'function') { var observer = new window.IntersectionObserver(function (entries) { for (n = 0; n !== entries.length; n = n + 1) { if (entries[n].isIntersecting) { entries[n].target.style.opacity = '1'; entries[n].target.style.transform = 'translateY(0)'; observer.unobserve(entries[n].target); } } }, { threshold: 0.1 }); for (m = 0; m !== flows.length; m = m + 1) { flows[m].style.opacity = '0'; flows[m].style.transform = 'translateY(20px)'; flows[m].style.transition = 'opacity 0.6s ease, transform 0.6s ease'; observer.observe(flows[m]); } } } /* 4. Auto-update tahun di footer (opsional) */ var yearSpans = document.querySelectorAll('.cvx-year'); var currentYear = new Date().getFullYear(); var p; for (p = 0; p !== yearSpans.length; p = p + 1) { yearSpans[p].textContent = String(currentYear); } }); }());