Hardy-Weinberg Denge Hesaplayıcısı
Popülasyonun Hardy-Weinberg dengesinde olup olmadığını test etmek için dominant homozigot, heterozigot ve resesif homozigot bireylerin gözlenen sayılarını girin; alel frekanslarını (p ve q), beklenen genotip frekanslarını ve ki-kare istatistiğini hesaplayın.
Hardy-Weinberg Frekansları ve Ki-Kare Testi
Hardy-Weinberg Denge Formülü
Burada:
- \(CH\) — dominant homozigot (AA) birey sayısı
- \(H\) — heterozigot (Aa) birey sayısı
- \(RH\) — resesif homozigot (aa) birey sayısı
- \(N\) — toplam popülasyon büyüklüğü, \(N = CH + H + RH\)
- \(p\), \(q\) — A ve a alellerinin frekansları (\(p + q = 1\))
- \(E_{CH}\), \(E_H\), \(E_{RH}\) — denge durumunda beklenen genotip sayıları
- \(\chi^2\) — dengeden sapmayı ölçen ki-kare istatistiği
İlk iki denklem genotip sayılarından alel frekanslarını tahmin eder (heterozigotlar her alelden bir kopya taşır), ortadaki üçü Hardy-Weinberg dengesi altında beklenen genotip sayılarını verir ve sonuncusu, gözlenen sayıların bu beklentilerden ne kadar saptığını nicelleştirir.
Hardy-Weinberg Dengesi – Hesaplama Örneği
Örnek: 100 bireyden oluşan bir örneklemde 50 dominant homozigot (AA), 30 heterozigot (Aa) ve 20 resesif homozigot (aa) birey bulunmaktadır.
\[ p = \frac{50 + 0.5 \cdot 30}{100} = 0.65, \qquad q = \frac{20 + 0.5 \cdot 30}{100} = 0.35 \]
\[ E_{CH} = 0.65^2 \cdot 100 = 42.25, \quad E_H = 2(0.65)(0.35) \cdot 100 = 45.5, \quad E_{RH} = 0.35^2 \cdot 100 = 12.25 \]
\[ \chi^2 = \frac{(50 - 42.25)^2}{42.25} + \frac{(30 - 45.5)^2}{45.5} + \frac{(20 - 12.25)^2}{12.25} \approx 11.61 \]
Sonuç: \(p = 0.65\), \(q = 0.35\); beklenen sayılar 42,25 AA, 45,5 Aa, 12,25 aa; \(\chi^2 \approx 11.61\). 11,61 değeri kritik değer olan 3,84'ü (1 serbestlik derecesi, \(\alpha = 0.05\)) aştığından, popülasyon Hardy-Weinberg dengesinden anlamlı biçimde sapmaktadır — heterozigot eksikliği söz konusudur.
Hesaplayıcı Nasıl Kullanılır?
- İlk alana dominant homozigot birey (AA genotipi) sayısını girin.
- İkinci alana heterozigot birey (Aa genotipi) sayısını girin.
- Üçüncü alana resesif homozigot birey (aa genotipi) sayısını girin.
- Sonuçları okuyun: p ve q alel frekansları, her genotipin beklenen sayıları ve Ki-kare (χ²) istatistiği.
- Popülasyonun denge içinde olup olmadığını değerlendirmek için ki-kare değerini aşağıdaki kritik değerle karşılaştırın.
Hardy-Weinberg Dengesi Nedir?
Hardy-Weinberg ilkesi, evrimsel kuvvetlerin bulunmadığı durumlarda bir popülasyondaki alel ve genotip frekanslarının kuşaktan kuşağa sabit kaldığını belirtir. Frekansları \(p\) ve \(q\) olan (\(p + q = 1\)) A ve a olmak üzere iki alelli bir gen için, bir kuşak rastgele çiftleşme sonrasında genotip frekansları \(p^2\) (AA), \(2pq\) (Aa) ve \(q^2\) (aa) olur. Bu ilişki yalnızca ideal koşullar altında geçerlidir: mutasyon olmaması, göç (migrasyon) olmaması, doğal seçilim olmaması, sonsuz büyük popülasyon büyüklüğü (genetik sürüklenme olmaması) ve rastgele çiftleşme.
Gerçek popülasyonlar bu koşulların tümünü nadiren tam olarak sağladığından, Hardy-Weinberg dengesi bir boş (null) model olarak hizmet eder. Gözlenen genotip sayıları beklenen \(p^2 : 2pq : q^2\) oranlarıyla uyumluysa, o lokusta evrimsel kuvvetlerin etkin olduğuna dair kanıt yoktur. Uyum yoksa sapma, bozulan bir veya daha fazla varsayıma işaret eder — seçilim, rastgele olmayan çiftleşme (inbreeding), genetik sürüklenme, mutasyon, göç veya genotipleme hatası.
Ki-kare uyum iyiliği testi, gözlenen ve beklenen sayılar arasındaki uyuşmazlığı nicelleştirir. Alel frekansları aynı veriden tahmin edildiği için testin \(3 - 1 - 1 = 1\) serbestlik derecesi vardır. Geleneksel 0,05 anlamlılık düzeyinde kritik değer 3,84'tür: 3,84'ün altındaki bir χ² değeri dengeyle uyumludur; daha büyük bir değer ise anlamlı bir sapmaya işaret eder. Alel frekansları p ve q, denge testi reddedilse bile yararlı özet istatistikler olmaya devam eder.
Bu araç popülasyon genetiğinde, evrimsel biyolojide ve sınıf alıştırmalarında yaygın olarak kullanılır. Genetik danışmanlık bağlamında alel frekansları ayrıca taşıyıcılık riski tahmininin de temelini oluşturur — genotip kombinasyonlarının ebeveynden yavrulara nasıl geçtiğini görmek için Kan Grubu Kalıtımı Hesaplayıcısı'na bakın.
Ki-Kare Sonucunun Yorumlanması
| Serbestlik derecesi | χ² kritik değeri (α = 0,05) |
|---|---|
| 1 | 3,841 |
| 2 | 5,991 |
Bu test için, bir alel frekansı veriden tahmin edildiği için 1 serbestlik derecesi kullanın. χ² değeri 3,841'in altındaysa dengeye ilişkin sıfır hipotezini kabul edin; üzerindeyse sapma %5 düzeyinde istatistiksel olarak anlamlıdır.
Varsayımlar ve Sınırlamalar
- Hesaplayıcı, tam baskınlık gösteren (AA ve Aa, aa'dan fenotipik olarak ayırt edilebilir ve genotipler doğrudan sayılabilir, örn. genetik testle) ve tam olarak iki alele sahip tek bir otozomal lokus varsayar.
- Denge; rastgele çiftleşme, seçilim olmaması, mutasyon olmaması, göç olmaması ve büyük popülasyon gerektirir. Anlamlı bir χ² değeri, bu varsayımlardan en az birinin bozulduğunu gösterir — hangisinin olduğunu söylemez.
- Beklenen sayıların küçük olduğu durumlarda (pratik bir kural olarak herhangi bir \(E < 5\)) ki-kare yaklaşımı güvenilmezdir. Bu durumda genotip sınıflarını birleştirin veya daha büyük bir örneklem toplayın.
- Genotipleme hataları, örnekleme yanlılığı ve bir araya getirilen alt popülasyonlar (Wahlund etkisi), dengesizliği taklit eden görünür heterozigot eksikliklerine yol açabilir.
- Test yalnızca tek bir lokustaki sapmayı saptar; tüm bir popülasyonla ilgili sonuçlara varmak için birden fazla lokusun incelenmesi gerekir.
Sıkça Sorulan Sorular
Hesaplayıcı olmadan alel frekanslarını elle nasıl hesaplarım?
İlgili alelin toplam kopya sayısını bulun ve örneklemdeki tüm gen kopyalarına bölün. Her homozigot tek bir alelden iki kopya taşır, her heterozigot ise her alelden birer kopya taşır. Dolayısıyla p = (2·AA + Aa) / (2N) ve q = (2·aa + Aa) / (2N); burada N birey sayısıdır. Toplamda 2N alel bulunduğundan p + q tam olarak 1'e eşit olmalıdır — bu, aritmetiğinizin doğru olduğuna dair hızlı bir kontroldür.
Bu test neden 2 yerine 1 serbestlik derecesi kullanıyor?
Üç genotip sınıfı bulunduğundan normalde sabit toplam için bir serbestlik derecesi kaybedersiniz (3 − 1 = 2). Ancak burada p alel frekansı aynı örneklemden tahmin edilir ve q, q = 1 − p ile belirlenir; böylece bir serbestlik derecesi daha kaybedilir: 3 − 1 − 1 = 1. Bu yüzden %5'lik kritik değer 5,99 değil 3,84'tür. Yanlış serbestlik derecesi kullanmak, Hardy-Weinberg alıştırmalarında en sık yapılan hatalardan biridir.
Hardy-Weinberg dengesi için iyi veya normal bir ki-kare değeri nedir?
"İdeal" bir değer yoktur — istatistik, örneklem büyüklüğüne ve alel frekanslarının kendisine bağlıdır. Dengeyle uyumlu bir popülasyon, basitçe kritik değer olan 3,84'e (α = 0,05, 1 sd) kıyasla küçük bir χ² üretir. Çok küçük χ² değerleri, idealize edilmiş ders kitabı verilerinde yaygındır ve beklenir. Çok büyük değerler sapmaya işaret eder; ancak dikkat: büyük örneklemelerde çok küçük biyolojik farklar bile istatistiksel olarak anlamlı hale gelebilir, bu yüzden gözlenen ve beklenen sayıları ayrıca doğrudan karşılaştırın.
Alel frekansları negatif veya 1'den büyük olabilir mi?
Hayır. Alel frekansları, alel kopyalarının sayısının toplam gen kopyası sayısına (2N) bölünmesiyle elde edildiğinden hem p hem de q her zaman 0 ile 1 arasındadır ve p + q tam olarak 1'e eşittir. Bu aralığın dışında bir değer alırsanız aritmetik veya veri girişi hatası yapmışsınızdır — örneğin fenotip sayılarını (AA ile Aa'nın ayırt edilemediği durumlar) genotip sayıları sanmak gibi. Genotip frekanslarının 0 olabileceğini, ancak alel frekanslarının yalnızca bir alel tamamen yoksa 0 veya 1'e ulaşabildiğini unutmayın.
Örneklemimde bir genotip tamamen eksikse ne olur?
Eksik bir genotip sınıfı p ve q'nun tahmini için sorun değildir — örneğin 0 AA, 40 Aa, 20 aa yine geçerli frekanslar verir (p = 0,33, q = 0,67). Asıl sorun ki-kare testindedir: herhangi bir beklenen sayı yaklaşık 5'in altındaysa yaklaşım güvenilmez hale gelir. Nadir alellerde beklenen homozigot sayıları çoğu zaman çok küçüktür. Bu durumda daha büyük bir örneklem toplayın ya da küçük örneklemeleri doğru ele alan bir exact test kullanın (genetik yazılımlarında kullanılan exact Hardy-Weinberg testi gibi).
Alel frekanslarını, resesif bir hastalık için taşıyıcı sıklığını tahmin etmek amacıyla nasıl kullanabilirim?
Hastalık resesifse ve lokus dengedeyse, etkilenen bireyler (aa) q² sıklığında görülür; dolayısıyla q = √(etkilenenlerin sıklığı) şeklinde tahmin edebilirsiniz. Etkilenmemiş kişiler arasındaki taşıyıcı sıklığı yaklaşık 2pq'dir. Örneğin 2.500 yenidoğandan 1'i etkileniyorsa q = 1/50 = 0,02 olur ve popülasyonun yaklaşık %3,9'u tek kopya taşır. Bu, dengenin geçerli olduğunu varsayar; gerçek popülasyonlar sapabildiğinden, pratikte gerçek risk değerlendirmesi için bir genetik danışmana başvurun.
Ki-kare değerim anlamlı ama popülasyon iyi görünüyor — sahte bir sapmaya ne yol açabilir?
Birkaç teknik neden gerçek dengesizliği taklit edebilir: heterozigotları yanlış okuyan genotipleme hataları, genetik olarak farklı alt popülasyonlardan gelen verilerin birleştirilmesi (görünür homozigotluğu şişiren Wahlund etkisi), rastgele olmayan örnekleme ya da yalnızca küçük beklenen sayıların χ² yaklaşımını hassasiyetsiz kılması. Ayrıca genotiplerin gerçekten ayırt edilip edilmediğini de kontrol edin — AA ve Aa, fenotiplerden eksik baskınlık varsayımlarıyla çıkarıldıysa, girilen verinin kendisi yanlış olabilir. Birden fazla bağımsız lokusun test edilmesi, bir veri sorununu gerçek biyolojik bir sinyalden ayırmaya yardımcı olur.