アレル頻度計算機 - ハーディー・ワインベルグ

遺伝子型カウントまたは直接のアレル入力から、任意の集団のアレル頻度、遺伝子型頻度、ハーディー・ワインベルグ平衡を計算します。

遺伝子型カウント(AA、Aa、aa)または直接のアレル数を入力して、p、q、期待遺伝子型頻度、HWE のカイ二乗検定を計算します。

アレル頻度計算機 - ハーディー・ワインベルグ
遺伝子型カウントまたは直接のアレル入力から、任意の集団のアレル頻度、遺伝子型頻度、ハーディー・ワインベルグ平衡を計算します。

アレル頻度計算機について

アレル頻度計算機は、集団内の遺伝的変異を特徴づける基本量、すなわち各アレルの頻度、各遺伝子型の頻度、そしてその集団がハーディー・ワインベルグ平衡(HWE)にあるかどうかを計算する集団遺伝学ツールです。これらの計算は、遺伝学や進化生物学の授業、臨床遺伝学、保全生物学、集団レベルの DNA データを扱うあらゆる研究で不可欠です。 アレル頻度は、集団遺伝学で最も単純かつ重要な要約統計量です。優性アレル A と劣性アレル a からなる二アレル系では、アレル頻度は p(A の頻度)と q(a の頻度)で表され、p + q = 1 を満たす必要があります。遺伝子型カウントから計算する場合、p = (2 × AA + Aa) / (2 × N) で、N は二倍体個体の総数です。各二倍体個体は 2 コピーを持つため、分母 2N は集団内の総アレル数を表します。直接のアレルカウントからは、p = A / (A + a)、q = a / (A + a) となります。 遺伝子型頻度は、集団のうち 3 つの二倍体遺伝子型クラス、AA(優性ホモ接合)、Aa(ヘテロ接合)、aa(劣性ホモ接合)を持つ割合を表します。観測遺伝子型頻度は、単に各カウントを N で割ったものです。ハーディー・ワインベルグ平衡下の期待遺伝子型頻度は、ハーディー・ワインベルグの原理により、それぞれ p²、2pq、q² と予測されます。 ハーディー・ワインベルグの原理は、大きくランダム交配する集団で、突然変異、移住、自然選択、遺伝的浮動がない場合、アレル頻度と遺伝子型頻度は世代を超えて一定に保たれると述べています。HWE 下の期待遺伝子型数は、AA が p² × N、Aa が 2pq × N、aa が q² × N です。観測遺伝子型数がこれらの期待値と有意に異なる場合、その集団は HWE から外れているとされ、近親交配、集団構造化、その座位に働く自然選択、または実験室データのジェノタイピングエラーを示している可能性があります。 カイ二乗検定は、HWE からの逸脱が偶然に期待される範囲を超えているかどうかを判断するための正式な統計的枠組みを提供します。検定統計量は、3 つの遺伝子型クラスについて (観測値 − 期待値)² / 期待値 を合計したものです。自由度 1(2 アレル、3 遺伝子型、1 制約)では、p = 0.05 の臨界値は 3.841 です。カイ二乗値が 3.841 未満なら HWE と矛盾せず、それを超える値は統計的に有意な逸脱を示します。 実用的な応用は多岐にわたります。臨床遺伝学では、SNP が HWE に従うかどうかを調べることは標準的な品質管理手順です。対照群のジェノタイピングデータで体系的な HWE 逸脱がある場合、バッチ効果やジェノタイピングエラーが疑われます。保全生物学では、小規模または分断化された集団で HWE から逸脱すると、近交弱勢や遺伝的ボトルネックを示すことがあります。法医学遺伝学では、DNA プロファイルのランダム一致確率を計算する積の法則の前提に HWE が使われます。疾患関連アレルの p と q を理解することで、疫学者は 2pq 式を用いて保因者頻度を推定できます。

アレル頻度の例

遺伝子型カウント、アレル頻度、ハーディー・ワインベルグ期待値を示す実際の集団シナリオです。

集団p / qHWE 評価
50 AA、30 Aa、20 aa (N=100)p = 0.6500、q = 0.3500期待値: AA 42.25、Aa 45.50、aa 12.25。カイ二乗 ≈ 7.14 — HWE から逸脱。
10 AA、80 Aa、10 aa (N=100)p = 0.5000、q = 0.5000ほとんどがヘテロ接合体。期待値: AA 25、Aa 50、aa 25。カイ二乗 ≈ 36 — HWE から大きく逸脱(ヘテロ接合体過剰)。
A アレル 120、a アレル 80p = 0.6000、q = 0.4000直接アレル入力モード。N は 100 と推定。期待値: AA 36、Aa 48、aa 16。
3 AA、2 Aa、5 aa (N=10)p = 0.4000、q = 0.6000小集団。期待値: AA 1.6、Aa 4.8、aa 3.6。サンプルサイズが小さいと HWE 検定は信頼性が低くなります。

アレル頻度計算機の使い方

  1. 入力方法を選びます。AA、Aa、aa のカウントがある場合は「遺伝子型カウント」、各アレルの直接カウントがある場合は「アレルカウント」を選択します。
  2. 遺伝子型カウントの場合: 優性ホモ接合(AA)、ヘテロ接合(Aa)、劣性ホモ接合(aa)の個体数を入力します。
  3. アレルカウントの場合: A アレルと a アレルの数を入力します。スケールされた HWE 期待カウントを得るには、任意で総個体数も入力します。
  4. 「計算」をクリックします。ツールは p と q、総 N、HWE 下の期待遺伝子型数、そして(遺伝子型入力の場合)HWE 逸脱のカイ二乗検定を表示します。
  5. 「リセット」をクリックすると、すべての項目を消去して新しい計算を開始できます。

アレル頻度計算機 FAQ

遺伝学で p と q は何を表しますか?
二アレル系では、p は優性アレル(A)の頻度、q は劣性アレル(a)の頻度です。座位上のすべてのアレルは A または a のいずれかであるため、常に合計は 1(p + q = 1)になります。アレル頻度は 0(アレルが存在しない)から 1(アレルが集団内で固定されている)までの範囲を取ります。
ハーディー・ワインベルグ平衡とは何ですか?
ハーディー・ワインベルグ平衡(HWE)は、集団内のアレル頻度と遺伝子型頻度が世代を超えて一定に保たれる理論的状態を表します。これには、大きな集団サイズ、ランダム交配、突然変異なし、移住なし、自然選択なしという 5 つの条件が必要です。HWE 下では、期待遺伝子型頻度は p²(AA)、2pq(Aa)、q²(aa)です。
なぜ集団が HWE から逸脱することがありますか?
HWE 逸脱の一般的な原因には、近親交配(ヘテロ接合性を低下させる)、集団構造化(遺伝的に異なる亜集団の混合)、その座位への強い自然選択、最近の移住、または実験データのジェノタイピングエラーがあります。ホモ接合体過剰は近親交配を示唆し、ヘテロ接合体過剰は集団構造化またはヘテロ接合体優位を示すことがあります。
HWE のカイ二乗検定はどのように解釈しますか?
カイ二乗統計量は、観測遺伝子型数と期待遺伝子型数の差を測定します。自由度 1 では、3.841 未満の値は 5% 水準で統計的に有意ではなく(HWE と一致)、3.841 を超える値は HWE からの有意な逸脱を示します。なお、サンプルサイズが非常に小さい場合(N < 20)、カイ二乗近似は信頼性が低く、正確検定が望まれます。
この計算機は X 連鎖座位に使えますか?
ここに示す標準的なハーディー・ワインベルグ式は常染色体座位に適用されます。X 連鎖座位では、男性は半接合(アレルを 1 つだけ持つ)であるため、男性における劣性表現型の頻度は q² ではなく q に等しくなります。別個の X 連鎖 HWE 解析が必要ですが、女性の遺伝子型カウントからアレル頻度(p と q)を計算する場合は同じ式を使います。
一般的な疾患変異の典型的なアレル頻度はどれくらいですか?
ゲノムワイド関連解析(GWAS)で見つかる一般的な疾患関連 SNP は、しばしば 0.05 から 0.50 のマイナーアレル頻度(MAF)を持ちます。まれなメンデル性疾患変異では、頻度が 0.01(1%)を大きく下回ることがあります。劣性疾患の保因者頻度は 2pq に等しいため、劣性アレル q = 0.01 の場合、保因者頻度はおよそ 2 × 0.99 × 0.01 ≈ 2% です。