Simulator Hanyutan Genetika Populasi

Model Wright–Fisher dalam Genetika Populasi

Simulator Hanyutan Genetika Populasi menerapkan model diploid netral Wright–Fisher untuk mensimulasikan dan memvisualisasikan bagaimana hanyutan genetik (genetic drift) memengaruhi frekuensi alel dari generasi ke generasi [i, j]. Hanyutan genetik adalah proses stokastik di mana frekuensi alel berfluktuasi secara acak karena kesalahan pengambilan sampel (sampling error) dari satu generasi ke generasi berikutnya.

Dalam model diploid Wright–Fisher, setiap generasi baru dibentuk melalui pengambilan sampel acak dengan pengembalian dari kumpulan gen tetua [i, j]. Jika suatu populasi memiliki ukuran populasi efektif sebesar Ne individu diploid, maka terdapat total 2Ne salinan gen dalam kumpulan gen tersebut. Proses transisi dari generasi t ke generasi t+1 dimodelkan menggunakan distribusi binomial.

Pada setiap generasi, transisi Wright–Fisher mengambil sampel jumlah alel berikutnya dan membaginya dengan salinan gen yang tersedia: Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne).

Sebelum pengambilan sampel, momen-momen ini menggambarkan rentang kemungkinan generasi berikutnya: E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne).

Melalui pengambilan sampel acak ini, frekuensi alel mengalami "random walk". Tanpa adanya kekuatan evolusioner lain seperti mutasi, migrasi, atau seleksi, fluktuasi acak ini pada akhirnya akan menyebabkan salah satu alel hilang (p = 0) atau terfiksasi (p = 1) [i, j]. Setelah frekuensi mencapai batas absorpsi ini, frekuensi alel tidak dapat berubah lagi tanpa adanya mutasi atau aliran gen dari luar.

---

Ukuran Populasi Efektif vs. Ukuran Sensus

Dalam genetika populasi, terdapat perbedaan penting antara ukuran populasi sensus (jumlah total individu nyata) dan ukuran populasi efektif, Ne. Parameter yang dimasukkan dalam simulator ini adalah ukuran populasi efektif, Ne, yang diukur dalam individu pembiak.

Ukuran populasi efektif didefinisikan sebagai ukuran populasi ideal yang mengalami tingkat hanyutan genetik (atau laju kehilangan heterozigositas) yang sama dengan populasi nyata yang sedang diamati. Dalam populasi alami, Ne hampir selalu lebih kecil daripada ukuran sensus nyata karena beberapa faktor biologis:

• Rasio Jenis Kelamin yang Tidak Seimbang: Jika jumlah jantan dan betina yang bereproduksi tidak sama, kontribusi genetik ke generasi berikutnya menjadi tidak seimbang, yang mempercepat hanyutan genetik.
• Varians dalam Keberhasilan Reproduksi: Jika beberapa individu menghasilkan keturunan jauh lebih banyak daripada yang lain, ukuran populasi efektif akan menurun.
• Fluktuasi Ukuran Populasi: Peristiwa penyusutan populasi (bottleneck) secara dramatis menurunkan Ne, karena ukuran populasi efektif jangka panjang sangat dipengaruhi oleh rata-rata harmonik dari ukuran populasi dari waktu ke waktu.

---

Dinamika Fiksasi, Lenyap, dan Heterozigositas

Hanyutan genetik secara inheren mengurangi variasi genetik di dalam populasi seiring berjalannya waktu [i, j]. Salah satu ukuran variasi genetik ini adalah heterozigositas (H), yaitu proporsi individu yang membawa dua alel berbeda pada lokus yang diamati.

Proporsi pasangan alel berbeda yang diharapkan menurun selama t generasi sebagai berikut: H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ [i, j].

Formula ini menunjukkan bahwa laju penurunan heterozigositas berbanding terbalik dengan ukuran populasi efektif. Pada populasi kecil, nilai 1 ÷ (2Ne) menjadi besar, sehingga heterozigositas menurun dengan sangat cepat. Sebaliknya, pada populasi besar, penurunan variasi genetik terjadi jauh lebih lambat.

Secara teoritis, probabilitas akhir bahwa suatu alel netral akan terfiksasi (mencapai p = 1) sama dengan frekuensi awalnya, p₀ [j]. Sebaliknya, probabilitas bahwa alel tersebut akan lenyap (mencapai p = 0) adalah 1 − p₀ [j]. Namun, dalam simulasi dengan batas generasi yang terbatas (T), tidak semua populasi replikasi akan mencapai batas absorpsi ini sebelum simulasi berakhir.

Meskipun masing-masing populasi independen akan menyimpang dan akhirnya mengalami fiksasi atau lenyap, rata-rata frekuensi alel di banyak populasi replikasi akan tetap mendekati frekuensi awal p₀ [j]. Fenomena ini menunjukkan bahwa meskipun variasi genetik hilang di dalam setiap populasi individu, variasi tersebut dipertahankan ketika melihat seluruh metapopulasi secara kolektif.

---

Metode, Asumsi, dan Batasan

Halaman ini menerapkan model diploid netral Wright–Fisher yang dijelaskan oleh Nature Education dan teks Stanford Human Genetics [j]. Setiap generasi baru adalah sampel acak dari 2Ne salinan alel dari generasi sebelumnya. Sumber diperiksa pada 6 Agustus 2026 [i, j].

Asumsi Model

Model Wright–Fisher yang diidealkan ini bergantung pada asumsi-asumsi ketat berikut [i, j]:

1. Ukuran populasi efektif konstan dari generasi ke generasi.
2. Generasi diskret dan tidak tumpang tindih.
3. Pengambilan sampel acak independen (perkawinan acak sempurna).
4. Netralitas selektif (semua genotipe memiliki kebugaran atau fitness yang sama).
5. Satu lokus dengan dua alel yang tidak terpaut dengan lokus lain.
6. Tidak ada mutasi baru yang memunculkan alel baru.
7. Tidak ada migrasi atau aliran gen dari populasi lain.

Batasan Penggunaan

Jangan gunakan simulator ini sebagai prediksi jika seleksi, perubahan populasi, ketidakseimbangan jenis kelamin atau keberhasilan reproduksi, generasi tumpang tindih, perkawinan sekerabat, pautan, aliran gen, atau struktur populasi berpengaruh. Faktor keselamatan teknik tidak berlaku untuk model biologis stokastik ini.

Setiap pengambilan binomial menggunakan pengambilan sampel Bernoulli yang tepat; tidak ada pendekatan normal yang digunakan. Perhitungan mempertahankan presisi ganda JavaScript. Frekuensi berada pada kelipatan 1 ÷ (2Ne); tampilan dibulatkan hingga maksimal enam tempat desimal.

Untuk mencegah kelebihan beban komputasi pada perangkat pengguna, browser menerapkan batas beban kerja maksimum sebesar 5.000,000 pengambilan salinan alel yang direncanakan, yang dihitung sebagai 2Ne × generasi × replikasi. Jika parameter yang dimasukkan melebihi batas ini, simulasi akan diblokir demi menjaga stabilitas kinerja.

---

Panduan Parameter dan Pengaturan Populasi

Untuk menjalankan simulasi, pengguna dapat memilih skenario bawaan (Presets) atau memasukkan parameter secara manual pada antarmuka Pengaturan populasi:

• Skenario: Menyediakan pilihan cepat seperti Populasi kecil, Populasi besar, atau Alel langka.
• Ukuran populasi efektif, Ne: Dimasukkan dalam satuan individu pembiak. Nilai harus berupa bilangan bulat dari 1 hingga 100,000.
• Frekuensi awal alel, p₀: Dimasukkan dalam satuan proporsi, 0–1. Nilai harus berada di antara 0 dan 1, termasuk kedua batas tersebut.
• Generasi: Dimasukkan dalam satuan generasi. Nilai harus berupa bilangan bulat dari 1 hingga 2,000.
• Populasi independen: Dimasukkan dalam satuan replikasi. Nilai harus berupa bilangan bulat dari 1 hingga 200.
• Seed acak: String teks (maksimal 100 karakter) yang digunakan untuk menginisialisasi generator angka acak semu. Menggunakan seed yang sama dengan parameter input yang sama akan menghasilkan jalur simulasi yang persis sama, sehingga eksperimen dapat direproduksi secara konsisten. Tombol Seed baru dapat digunakan untuk mengacak nilai ini.

---

Pemrosesan Data dan Privasi

Keamanan data pengguna adalah prioritas utama dalam penggunaan alat analisis ini. Parameter dan simulasi populasi Anda tetap berada di browser ini. Seluruh proses komputasi, mulai dari penarikan sampel binomial hingga pembuatan grafik interaktif, dilakukan secara lokal pada perangkat Anda. Tidak ada data parameter, seed, atau hasil simulasi yang diunggah ke server luar.

---

Pertanyaan yang Sering Diajukan (FAQ)

Haruskah saya memasukkan ukuran sensus atau ukuran populasi efektif? Masukkan ukuran populasi efektif, Ne: ukuran populasi pembiak ideal yang akan mengalami hanyutan pada tingkat yang sama dengan populasi yang ingin Anda pelajari. Ukuran sensus bisa jauh lebih besar ketika rasio jenis kelamin, keberhasilan reproduksi, fluktuasi populasi, atau struktur populasi tidak merata.

Apakah probabilitas netral untuk fiksasi akhir sama dengan frekuensi awal? Ya, dalam model netral tanpa mutasi atau migrasi ini, alel yang dimulai pada frekuensi p₀ memiliki probabilitas fiksasi akhir p₀ dan probabilitas lenyap 1 − p₀ [j]. Tingkat yang ditunjukkan di atas hanya mencakup jumlah generasi dan pengulangan yang dimasukkan, sehingga nilainya akan bervariasi tergantung pada seed.

Mengapa frekuensi rata-rata alel dapat tetap mendekati p₀ sementara heterozigositas menurun? Hanyutan netral tidak memiliki arah pilihan di banyak populasi independen, sehingga frekuensi rata-rata tetap p₀ dalam ekspektasi [j]. Namun, populasi individu menyebar mendekati 0 atau 1, dan kedua batas tersebut memiliki heterozigositas nol [i, j]. Oleh karena itu, variasi akan hilang meskipun rata-rata lintas populasi hampir tidak berubah [j].

Apakah ini dapat memprediksi apa yang akan terjadi pada populasi nyata? Tidak. Ini menunjukkan hasil di bawah model netral Wright–Fisher yang sengaja diidealkan [i, j]. Prediksi nyata memerlukan estimasi yang dapat dipertanggungjawabkan mengenai ukuran efektif serta seleksi, mutasi, migrasi, perubahan ukuran populasi, generasi yang tumpang tindih, pautan, dan struktur populasi yang memengaruhi kasus tersebut.