Popülasyon Genetiği Sürüklenme Simülatörü

Wright–Fisher Modeli ve Genetik Sürüklenmenin Temelleri

Popülasyon Genetiği Sürüklenme Simülatörü, diploid popülasyonlarda nötr genetik sürüklenmenin etkilerini görselleştirmek ve analiz etmek için tasarlanmış bir araçtır. Simülatör, popülasyon genetiğinin temel taşlarından biri olan diploid Wright–Fisher modelini temel alır. Bu model, her yeni neslin bir önceki nesildeki alel havuzundan rastgele örnekleme yoluyla oluşturulduğu idealize edilmiş bir üreme sistemini simüle eder.

Nötr sürüklenme sürecinde, seçilim baskısı olmaksızın alel frekanslarının nesiller boyunca nasıl değiştiği gözlemlenebilir. Simülatör, kullanıcıların belirlediği parametreler doğrultusunda bağımsız popülasyonların izlediği yolları grafik üzerinde göstererek fiksasyon (alelin popülasyonda sabitlenmesi, p = 1) ve kayıp (alelin tamamen yok olması, p = 0) oranlarını hesaplar.

Popülasyon Kurulumu ve Parametre Girişleri

Simülasyonu başlatmak için kullanıcıların "Popülasyon kurulumu" arayüzünde belirli parametreleri tanımlaması gerekir. Bu parametreler, simülasyonun biyolojik sınırlarını ve hesaplama yükünü doğrudan belirler:

• Senaryolar: Hızlı başlangıç için hazır şablonlar sunar. Bu şablonlar arasında "Küçük popülasyon", "Büyük popülasyon" ve "Nadir alel" seçenekleri yer alır.
• Etkin popülasyon büyüklüğü, Ne: "üreyen birey" cinsinden girilir. Bu değer 1 ile 100,000 arasında bir tam sayı olmalıdır.
• Başlangıç alel frekansı, p₀: "oran, 0–1" cinsinden girilir. Her iki sınır dahil olmak üzere 0 ile 1 arasında bir değer olmalıdır.
• Nesil Sayısı: "nesil" cinsinden girilir. 1 ile 2,000 arasında bir tam sayı olmalıdır.
• Bağımsız Popülasyonlar: "tekrar" sayısıdır. 1 ile 200 arasında bir tam sayı olmalıdır.
• Rastgele başlangıç değeri (seed): Simülasyon yollarının birebir yeniden üretilebilmesini sağlayan, en fazla 100 karakterden oluşan bir metindir. Kullanıcılar dilerse "Yeni başlangıç değeri" butonuna tıklayarak bu değeri rastgele değiştirebilirler.

İş Yükü Sınırları ve Hata Yönetimi

Tarayıcı tabanlı hesaplama performansını korumak amacıyla, simülatör maksimum 5,000,000 planlanan alel kopyası çekimine izin verir. Bu sınır 2Ne × nesil sayısı × bağımsız popülasyon sayısı formülüyle hesaplanır. Bu sınır aşıldığında simülasyon bloke edilir ve kullanıcıya şu hata mesajı gösterilir:

"Bu kurulum, tarayıcının ‹limit› sınırının üzerinde, ‹draws› alel kopyası çekimi planlıyor. Popülasyon büyüklüğünü, nesil sayısını veya tekrar sayısını azaltın."

Diğer geçersiz girişlerde ise şu hata mesajları tetiklenir:

• Sayısal olmayan girişlerde: ‹field›: “‹token›” bir sayı değil.
• Ondalık sayı girildiğinde: ‹field›: bir tam sayı girin.
• Sınır dışı değerlerde:
  • "Etkin popülasyon büyüklüğü 1 ile 100,000 arasında bir tam sayı olmalıdır."
  • "Başlangıç alel frekansı, her iki sınır dahil olmak üzere 0 ile 1 arasında olmalıdır."
  • "Nesil sayısı 1 ile 2,000 arasında bir tam sayı olmalıdır."
  • "Bağımsız popülasyon sayısı 1 ile 200 arasında bir tam sayı olmalıdır."
• Eksik veya uzun başlangıç değerlerinde:
  • "Simülasyonun yeniden üretilebilmesi için rastgele bir başlangıç değeri (seed) girin."
  • "Rastgele başlangıç değerini (seed) 100 karakter veya daha az tutun."

Yöntem, Varsayımlar ve Sınırlar

Bu simülatör, Nature Education ve Stanford Human Genetics kaynaklarında açıklanan nötr diploid Wright–Fisher modelini uygulamaktadır. Kaynaklar 6 Ağustos 2026'da kontrol edilmiştir.

Matematiksel Model ve Örnekleme

Her nesilde, Wright–Fisher geçişi bir sonraki alel sayısını örnekler ve bunu mevcut gen kopyalarına böler:

Xₜ₊₁ ~ Binomial(2Ne, pₜ); pₜ₊₁ = Xₜ₊₁ ÷ (2Ne)

Örneklemeden önce, bu momentler bir sonraki olası nesillerin aralığını tanımlar:

E[Xₜ₊₁] = 2Ne·pₜ; Var(Xₜ₊₁) = 2Ne·pₜ(1 − pₜ); Var(pₜ₊₁) = pₜ(1 − pₜ) ÷ (2Ne)

Farklı alel çiftlerinin beklenen payı (heterozigotluk), t nesil boyunca şu şekilde azalır:

H₀ = 2p₀(1 − p₀); Hₜ = H₀(1 − 1 ÷ (2Ne))ᵗ

Her binomial çekim kesin Bernoulli örneklemesi kullanır; normal yaklaşım kullanılmaz. Hesaplamalar JavaScript çift duyarlılığını korur. Frekanslar 1 ÷ (2Ne) artışlarında yer alır; ekran en fazla altı ondalık basamağa yuvarlar.

Modelin Varsayımları

Simülasyonun arkasındaki matematiksel model belirli basitleştirici varsayımlara dayanmaktadır:

1. Sabit etkin popülasyon büyüklüğü.
2. Kesikli ve örtüşmeyen nesiller.
3. Rastgele bağımsız örnekleme.
4. Seçilimsel nötrlük (tüm aleller eşit uyum başarısına sahiptir).
5. Bağlantısız tek bir iki alelli lokus.
6. Mutasyon veya göç olmaması.

Seçilim, popülasyon değişimi, cinsiyet veya üreme başarısı dengesizliği, örtüşen nesiller, akraba evliliği, bağlantı, gen akışı veya popülasyon yapısı önem taşıdığında bu araç bir tahmin aracı olarak kullanılmamalıdır. Bu stokastik biyolojik model için bir mühendislik güvenlik faktörü geçerli değildir.

Simülasyon Sonuçları ve Çıktılar

Simülasyon çalıştırıldığında, arayüzde durum mesajları, grafiksel yollar ve özet istatistikler güncellenir.

Durum Mesajları

Simülasyonun farklı aşamalarında kullanıcıya şu mesajlar gösterilir:

• Hazır durumdayken: "Popülasyon senaryosunu değiştirin veya bu başlangıç değerli örneği tekrar çalıştırın."
• Hesaplama sırasında: "Sonraki nesiller örnekleniyor…"
• Tamamlandığında: "“‹seed›” başlangıç değeri kullanılarak ‹replicates› popülasyon ‹generations› nesil boyunca simüle edildi."

Grafiksel Gösterim

Alel frekansı yolları grafiği, "Tekrarlanan simüle edilmiş popülasyonlar için nesle göre alel frekansı" ARIA etiketine sahiptir. Grafiğin X ekseni "Nesil", Y ekseni ise "Alel frekansı p" olarak etiketlenmiştir. Grafik üzerinde "Beklenen p = p₀" referans çizgisi yer alır. Grafik altında ise kaç yolun çizildiği şu formatta belirtilir:

"‹total› yoldan ‹shown› tanesi gösteriliyor; her bir tekrar özete dahil edilmiştir."

Özet İstatistikler ve Değer Değişimleri

Simülasyon tamamlandıktan sonra aşağıdaki özet veriler sunulur:

• p = 1 değerinde fikse olmuş: Fiksasyona ulaşan popülasyonların oranı.
• p = 0 değerinde kaybolmuş: Alelin tamamen yok olduğu popülasyonların oranı.
• T anına kadar absorbe olmamış: Belirlenen nesil süresince fiksasyona veya kayba uğramamış popülasyonlar.
• Ortalama nihai p: Tüm tekrarların son nesildeki ortalama alel frekansı.
• T anında beklenen H: Teorik olarak beklenen heterozigotluk değeri.
• T anında ortalama simüle edilen H: Simülasyonda gözlemlenen ortalama heterozigotluk.
• Ortalama absorpsiyon nesli: Popülasyonların fiksasyon veya kayba ulaşma nesillerinin ortalaması (eğer hiçbir popülasyonda absorpsiyon gerçekleşmediyse "ulaşılmadı" yazar).

Ayrıca ilk nesle ait teorik ve gözlemlenen değerlerin karşılaştırması şu şekilde sunulur:

• "İlk nesil 2 × ‹population› = ‹copies› alel kopyası çekimi kullanır."
• "E[X₁] = ‹copies› × ‹frequency› = ‹mean› A kopyası."
• "İlk nesil için Var(X₁) = ‹copies› × ‹frequency› × (1 − ‹frequency›) = ‹variance›; SD(p₁) = ‹sd›."
• "Grafiğe dökülen ilk tekrarda X₁ = ‹sampled› A kopyası örneklendi, dolayısıyla p₁ = ‹sampled› ÷ ‹copies› = ‹frequency›."

Veri Gizliliği ve İşlem Güvenliği

Kullanıcıların simülatör üzerinde yaptığı tüm işlemler ve girdikleri parametreler tamamen yerel düzeyde kalır. Parametreleriniz ve simüle edilen popülasyonlar bu tarayıcıda kalır. Hiçbir veri dış sunuculara yüklenmez veya üçüncü taraflarla paylaşılmaz.

---

Sıkça Sorulan Sorular (FAQ)

Sayım büyüklüğünü mü yoksa etkin popülasyon büyüklüğünü mü girmeliyim?

Etkin popülasyon büyüklüğünü (Ne) girin: incelemek istediğiniz popülasyonla aynı oranda sürüklenecek idealize edilmiş bir üreme popülasyonunun büyüklüğü. Cinsiyet oranları, üreme başarısı, popülasyon dalgalanmaları veya popülasyon yapısı dengesiz olduğunda sayım büyüklüğü çok daha büyük olabilir.

Nihai fiksasyonun nötr olasılığı başlangıç frekansına eşit midir?

Evet, mutasyon veya göçün olmadığı bu nötr modelde, p₀ başlangıç frekansına sahip bir alelin nihai fiksasyon olasılığı p₀, kaybolma olasılığı ise 1 − p₀'dır. Yukarıda gösterilen oranlar yalnızca girilen nesil ve tekrar sayısını kapsar, bu nedenle başlangıç değerine (seed) göre değişiklik gösterecektir.

Heterozigotluk düşerken ortalama alel frekansı neden p₀ civarında kalabilir?

Nötr sürüklenmenin birçok bağımsız popülasyon genelinde tercih edilen bir yönü yoktur, bu nedenle ortalama frekansın beklenen değeri p₀ olarak kalır. Ancak bireysel popülasyonlar 0 veya 1'e doğru yayılır ve her iki sınırda da heterozigotluk sıfırdır. Bu nedenle, popülasyonlar arası ortalama neredeyse hiç değişmeden kalırken varyasyon kaybedilir.

Bu simülatör gerçek bir popülasyonun ne yapacağını tahmin edebilir mi?

Hayır. Bu simülatör, kasıtlı olarak idealize edilmiş nötr bir Wright–Fisher modeli altındaki sonuçları gösterir. Gerçek tahminler; etkin büyüklüğün ve durum için önem taşıyan seçilim, mutasyon, göç, değişen popülasyon büyüklüğü, örtüşen nesiller, bağlantı ve popülasyon yapısı gibi faktörlerin savunulabilir tahminlerini gerektirir.