遺伝医学におけるLODスコアとは？連鎖解析の基礎から次世代シークエンサー（NGS）解析への応用まで臨床遺伝専門医が徹底解説

LODスコア（ロッドスコア）とは、2つの遺伝子座が同じ染色体上で近くにあり「一緒に遺伝するかどうか（＝連鎖）」を統計的に評価するための指標です。スコアが3.0以上なら連鎖あり（偶然の一致より1,000倍も確からしい）、反対に−2.0以下ならその領域から原因遺伝子を除外、と判断します。一見すると統計の専門的な話に見えますが、ハンチントン病やBRCA1/BRCA2など、数多くの病気の原因遺伝子は、このLODスコアを「羅針盤」にして発見されてきました。そして現在の全エクソーム解析（WES）・全ゲノム解析（WGS）による診断も、その統計学的な土台の上に成り立っています。

1. LODスコアとは：定義と遺伝医療における意義

LOD（ロッド）スコアは、「Logarithm of the Odds（オッズの対数）」の略称です。2つの遺伝子座が同じ染色体上で近くに存在し、一緒に次世代へ受け継がれる確率の相対的な高さを評価する統計学的な指標です。

ヒトのDNAをすべて取り出して一直線に伸ばすと、その長さは約2メートルにもなります。この広大な「体の設計図」の中から、たった一つの病的な変化（バリアント）を見つけ出す作業は、地球全体の地図の中から一つの番地を探し当てるような挑戦です。LODスコアは、この遺伝子の位置探索において歴史的・統計学的な基盤となってきました。なぜLODスコアを学ぶことが大切かといえば、それは「原因がわからない病気の遺伝子を特定する」という遺伝診療の出発点そのものだからです。

連鎖解析では、病気と「目印（遺伝子マーカー）」が家系の中でどのように伝わっていくかを観察します。そして、観察されたデータが「連鎖している」という仮説を支持するのか、それとも「偶然の一致（連鎖していない）」に過ぎないのかを、数学的に判定するためにLODスコアを用います。LODスコアが正の値なら連鎖の存在を示唆し、負の値なら連鎖していない可能性が高いことを示します。

2. 遺伝的連鎖の基礎：組換えがすべてのカギ

LODスコアを理解するには、まず「なぜ近い遺伝子は一緒に遺伝するのか」を知る必要があります。そのカギを握るのが、精子や卵子が作られる減数分裂のときに起こる組換え（乗換え）という現象です。

減数分裂では、父親由来と母親由来の相同染色体どうしが部分的にDNAを物理的に交換します。これが組換えです。2つの遺伝子座の距離が遠ければ、その間で組換えが起こりやすく、別々に遺伝します（メンデルの独立の法則）。逆に距離が近ければ、その間で組換えが起こる確率は著しく低くなり、2つはセットで一緒に遺伝します。この「組換えの起こりにくさ」を数値にしたものが、次の章で出てくる組換え価θ（シータ）です。

3. LODスコアの数学的構造と組換え価θ

LODスコアは、尤度比（ゆうどひ）の常用対数（底が10の対数）をとることで計算されます。連鎖の度合いを表す指標として、組換え価θ（シータ）を使います。

具体的には、ある組換え価θを仮定したときのデータの尤度 L(θ) と、連鎖がない（θ＝0.5）と仮定したときの尤度 L(0.5) の比の常用対数として定義されます。

計算例：最も確からしい距離を探す（最尤推定）

ある家系で、子どもが親から遺伝子を受け継ぐ過程を観察し、組換えを起こした子（R）が1人、起こさなかった子（NR）が5人（合計6回の減数分裂）だったとします。連鎖がない場合の尤度は 0.5の6乗＝0.015625 です。任意のθに対するLODスコアは次のように計算できます。

θの値をいろいろ変えてZを計算すると、下のグラフのような曲線が描けます。LODスコアが最大になるθの値が、そのデータから推測される最も確からしい遺伝的距離（最尤推定量）になります。この例では、θ＝約0.17（＝R÷(R＋NR)＝1/6）のときにLODスコアが最大（約0.63）となります。

組換えが1回・非組換えが5回観察された家系モデルにおけるLODスコアの推移。θが約0.17の地点でLODスコアが最大（約0.63）となり、これがこのデータから推測される最も可能性の高い遺伝的距離となる。

なぜ「対数」を使うのか：加法性という強み

LODスコアが対数を使う最大の理由は「加法性」にあります。遺伝性疾患の研究では、ひとつの家系だけで十分な証拠（LOD＞3.0）を得られることはまれで、複数の家系のデータを統合する必要があります。通常、確率（尤度）を統合するには掛け算が必要ですが、対数をとっていれば単なる「足し算」で統合できます。これにより、世界中の研究グループが集めた家系ごとのLODスコアを簡単に合算し、原因遺伝子のマッピング精度を高めることができるのです。

4. 「LODスコア＝3.0」という古典的な閾値の意味

連鎖解析では、長らく「LODスコアが3.0以上なら有意な連鎖がある」という基準が用いられてきました。なぜ「3.0」なのでしょうか。

LODスコア3.0は、尤度比でいうと10の3乗＝1,000倍を意味します。つまり「連鎖している確率が、偶然の一致である確率の1,000倍高い」ということです。しかし、ヒトゲノムの中からランダムに2つの遺伝子座を選んだとき、それらが偶然に同じ染色体上で連鎖している確率は約50分の1にすぎません。裏を返せば、「連鎖していない確率の方が約50倍高い」という事前確率が存在するのです。

したがって、観測データから1,000倍の連鎖の証拠（LOD 3.0）が得られても、この事前確率（1/50）を考慮すると、最終的な信頼性のオッズは 1,000 × (1/50) ＝ 20 となり、20対1（約95%の信頼性、偽陽性率約5%）に落ち着きます。LOD 3.0という一見高い閾値は、ゲノムの広大さに由来する低い事前確率を補正するために設定された、極めて合理的な数値なのです。なお、LODスコアが−2.0を下回る場合は、その染色体領域から原因遺伝子を明確に除外する強力な証拠となります。

LODスコアは、標準的な統計のp値（有意確率）にも変換できます。一般に、LODスコア3.0は、おおよそp値0.0001（1万分の1）に相当し、極めて強力な統計的証拠とみなされます。

5. 遺伝医学の歴史を切り拓いた統計学者たち

現代の精密な遺伝子検査は、数十年にわたるアルゴリズム構築の歴史の上に成り立っています。LODスコアにまつわる3つの大きな節目を紹介します。

① ニュートン・モートンの創唱（1955年）

LODスコアという概念をヒトの連鎖解析に初めて導入したのが、集団遺伝学者のニュートン・モートン（Newton E. Morton）です。1955年、彼は「連鎖検出のための逐次検定」という画期的な論文を発表し、常染色体上の連鎖を評価する数式を確立しました。これが、複雑な家系から情報を取り出して原因遺伝子をマッピングする出発点になりました。

② エルストン・スチュワート・アルゴリズム（1971年）

世代をまたぐ巨大な家系図で全員の遺伝確率を計算するのは、かつては天文学的な計算量が必要でした。1971年、ロバート・エルストンとジョン・スチュワートは、家系図の末端から順に確率を計算して情報を折りたたんでいく再帰的アルゴリズムを発表しました。これにより計算時間が家系の人数に対して直線的にしか増えなくなり、コンピューターによる大規模家系のLODスコア計算が現実のものとなりました。一方で、扱うマーカーの数に対しては計算量が急増するという制約があり、後のランダー・グリーン・アルゴリズムと相互補完的に使われています。

③ ランダーとクルーグリャクの提言（1995年）

1990年代、ゲノム全体を一度に調べるゲノムワイド・スキャンが可能になると、新たな問題が生じました。何千ものマーカーを同時に調べると、偶然にLOD 3.0を超える「偽陽性」が多発する（多重検定の問題）のです。1995年、エリック・ランダーとレオニード・クルーグリャクは、ゲノム全体での有意水準を約5%に保つために、古典的なLOD 3.0ではなくLOD 3.3を新たな基準とすべきだと提唱しました。また、LOD 2.2を「連鎖の示唆」と定義しました。

これらの統計基盤によって、嚢胞性線維症（CFTR遺伝子）、ハンチントン病（HTT遺伝子）、家族性アルツハイマー病、遺伝性乳がん・卵巣がん症候群（BRCA1／BRCA2遺伝子）など、数多くの画期的な原因遺伝子の特定が実現しました。

6. 連鎖解析の2つのアプローチ

LODスコアを用いた連鎖解析には、対象となる病気の性質によって2つのアプローチがあります。

パラメトリック連鎖解析（モデル依存型）

従来のLODスコア分析がこれにあたります。原因遺伝子が常染色体顕性（優性）か常染色体潜性（劣性）かといった遺伝形式、アレル（対立遺伝子）の頻度、そして浸透率などのパラメーターが既知であると仮定して計算します。遺伝モデルが正確に設定されていれば、単一遺伝子疾患（メンデル遺伝病）の原因遺伝子探索において絶大な威力を発揮します。一方で、遺伝モデルが不明確だったり浸透率が不完全だったりすると、誤って連鎖を否定してしまう（偽陰性の）リスクが高まります。

ノンパラメトリック連鎖解析（モデルフリー型）

遺伝形式が不明確な多因子疾患には、パラメーターの設定を必要としない手法が使われます。代表的なのが、病気を持つ兄弟姉妹のペア（罹患同胞対：ASP）を用いる方法です。病気の兄弟姉妹が、偶然（期待値50%）よりも高い確率で親から「同じDNA領域（IBD）」を受け継いでいるかどうかを検証します。

7. 次世代シークエンサー（NGS）とLODスコアの強力なシナジー

かつての連鎖解析は、原因遺伝子が潜む「広い候補領域」を特定するのが限界で、そこから真の原因変異を見つけるには膨大な時間がかかりました。しかし、DNAを高速・低コストで解読できる次世代シークエンサー（NGS）の登場で、状況は一変しました。

WESやWGSを行うと、一人から数万〜数百万個ものバリアントが検出されます。この膨大なノイズの中から、真の原因変異を絞り込む「強力なフィルター」として、古典的なLODスコア（連鎖解析）が再び脚光を浴びています。大家系のデータからLOD＞3.0で候補領域を絞り込み、その領域内のバリアントだけをNGSで精査することで、かつては診断不可能だった希少疾患の原因究明が飛躍的に進んでいます。実際に、精子鞭毛多発形態異常症（DNAH1遺伝子）や急性ポルフィリン症（HMBS遺伝子）などで、新規の原因遺伝子が特定された報告が相次いでいます。

「複数の臓器に症状がある」「発達に遅れがある」のに、従来の染色体検査やマイクロアレイ検査では原因が特定できないお子さんがいます。こうした原因不明の症状に対して、ミネルバクリニックでは全エクソーム検査（WES）や全ゲノム検査（WGS）を、必要に応じてトリオ解析として提供しています。これらの土台にあるのが、本記事で解説してきたLODスコアと連鎖解析の考え方です。

8. よくある誤解と、LODスコアの限界

限界として最も重要なのは、パラメトリック解析が遺伝モデルの設定に依存する点です。遺伝形式や浸透率の仮定が実態とずれると、本当は連鎖しているのに見逃す（偽陰性）危険があります。また、多くのマーカーを同時に調べると偽陽性が増えるため、ランダーとクルーグリャクが提唱した厳格な基準（ゲノムワイドでLOD 3.3）が用いられます。数値の解釈には、必ず専門的な判断が必要になるのです。

9. 臨床遺伝専門医からのメッセージ

よくある質問（FAQ）

関連記事

参考文献

• [1] National Human Genome Research Institute. LOD Score (Genetics Glossary). [genome.gov]
• [2] Morton NE. Sequential tests for the detection of linkage. Am J Hum Genet. 1955;7(3):277-318. [PubMed]
• [3] Elston RC, Stewart J. A general model for the genetic analysis of pedigree data. Hum Hered. 1971;21(6):523-542. [PubMed]
• [4] Lander E, Kruglyak L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995;11(3):241-247. [Nature Genetics]
• [5] Ott J, et al. Linkage analysis in the next-generation sequencing era. Hum Hered. 2011. [PMC3267991]
• [6] Hodge SE, et al. All LODs Are Not Created Equal. Am J Hum Genet. [PMC1287176]
• [7] Statistics for Human Genetics. Chapter 17: LOD scores (NYU). [NYU Biostatistics]