Inheritance Entropy: A Model-Independent Method to Probe the Hereditary Structure of Cell Lineage Trees

本論文は、「継承エントロピー（Inheritance Entropy）」と呼ばれるモデルに依存しない手法を導入することで、ヒト骨髄間葉系幹細胞コロニーにおける不均一性が、細胞周期停止を制御する遺伝子に依存しない継承可能な要因に由来することを実証しており、これは観察されたクローン系統樹の値が非継承モデルと比較して有意に低いエントロピー値を示したことによって裏付けられている。

要約

家族の系図を想像してみてください。ただし、それは人間ではなく、シャーレの中で成長する幹細胞のコロニーです。細胞が分裂するたびに、2つの新しい枝が生まれます。通常、これらの枝は成長し続け、樹形図を埋め尽くしていきます。しかし、時として、ある細胞が分裂をやめる決断をすることがあります。その細胞は「不活性状態」（眠りについたり、引退したりするような状態）に入ります。これが起こると、その家族の系図の枝全体が成長を停止します。

科学者たちが投げかけた大きな問いは、**「なぜ、成長を止める枝もあれば、成長し続ける枝もあるのか？」**ということでした。

それは単なるランダムな偶然なのでしょうか？ それとも、止まるべき時を伝える「家族の秘密」が、親から子へと受け継がれているのでしょうか？

問題：機械の中の「ゴースト」

これを研究することの難しい点は、細胞が分裂を停止すると、その細胞は将来の樹形図から事実上消滅してしまうことです。子供を持たなかった細胞の子供を研究することはできません。それは、存在しない葉を見て、なぜ特定の枝が枯れたのかを突き止めようとするようなものです。

通常、もし変異が起こったなら、それは波紋のように広がっていくのが見えます。しかし、「停止」の場合、その波紋はむしろ「穴」となります。科学者たちは、もし停止の決定が純粋にランダムなものだとしたら、樹形図の中の「穴」は、フロントガラスに当たる雨粒のようにいたるところに散らばるはずだと気づきました。しかし、もしその決定が遺伝的なもの（親から子へ受け継がれるもの）であれば、穴は特定の枝の中に固まって現れるはずです。

解決策：「遺伝エントロピー」の測定

この問題を解決するために、チームは「遺伝エントロピー（Inheritance Entropy）」と呼ばれる新しい数学的ツールを考案しました。

家族の系図を、散らかった部屋だと考えてみてください。

• 高いエントロピー（乱雑／ランダム）： もし「停止した」細胞が樹形図のいたるところにランダムに散らばっているなら、その部屋は混沌として予測不能に見えます。これは、停止が単なる不運やランダムなノイズによるものであることを示唆しています。
• 低いエントロピー（組織的／遺伝的）： もし「停止した」細胞が、特定の1つまたは2つの枝に固まっているなら、その部屋は奇妙に整然として見えます。これは、特定の枝に対して停止を命じる「ルール」が受け継がれたことを示唆しています。

科学者たちは、この「散らかり具合のスコア（エントロピー）」を、32種類のヒト幹細胞コロニーに対して算出しました。

発見：それはランダムではない

彼らは、実際の細胞コロニーを「シャッフルされた」バージョンと比較しました。想像してみてください、実際の家族の系図を取り出し、すべての「停止した」細胞を切り取り、それらをランダムに別の枝へと貼り直す様子を。これにより、遺伝が不可能な「偽の」樹形図が作成されます。

結果： 実際のコロニーの75%において、その「散らかり具合のスコア」は、シャッフルされた偽の樹形図よりも有意に低くなっていました。

• 言い換えれば： 「停止した」細胞はランダムに散らばっていたのではなく、特定の枝に固まっていました。これは、停止するという決定が遺伝していることを証明しています。それは単なるランダムな偶然ではなく、親細胞からその子供へと受け継がれる特性なのです。

「ラグ」効果：遅延反応

科学者たちはまた、これが「いつ」起こるかについても興味深い発見をしました。「変異ラグ」を見つけたのです。

ある細胞が「停止信号（変異）」を受け取ったと想像してください。その細胞はすぐに停止するわけではありません。数世代にわたって分裂を続け、その信号を子供たちへと受け渡していきます。数世代が経過して初めて、細胞はブレーキを踏み、停止するのです。

• 発見： 「停止信号」が生成されてから、細胞が実際に分裂を停止するまでには、平均して2世代の遅延があります。
• なぜ重要か： この遅延は極めて重要です。もし細胞が即座に停止していたら、科学者たちはその信号を源へと遡ることはできなかったでしょう。遅延があったからこそ、彼らは「空っぽの」枝を見て、信号が最初に生成された正確な瞬間まで経路を辿ることができたのです。

「家族の秘密」（エピジェネティクス）

論文は、この遺伝は恐らく**エピジェネティック（後成的なもの）**であると結論付けています。

• 比喩： DNAをコンピュータのハードウェア（決して変わることのないコード）だと考えてください。エピジェネティクスはソフトウェアの設定（音量、明るさ、あるいは「おやすみモード」のようなもの）です。
• 細胞はDNAのコード自体を変えているのではありません。代わりに、いつ作業を止めるかを指示する「設定」を受け継いでいるのです。これが、同じ種類の細胞から始まっていたとしても、なぜ幹細胞コロニー同士で見た目が大きく異なる場合があるのかを説明しています。

まとめ

1. 手法： 彼らは、細胞の死がランダムなのか、あるいは組織的なものなのかを見極めるために、「散らかり具合のスコア（エントロピー）」を用いました。
2. 発見： ほとんどのコロニーにおいて、死は組織的であり、「停止」の信号が遺伝していることが証明されました。
3. タイミング： 信号が送られてから細胞が実際に停止するまでには、2世代の遅延があります。
4. 原因： これは遺伝的な変異（ハードウェアの変化）ではなく、エピジェネティックな要因（ソフトウェアの設定）によるものである可能性が高いです。

この論文は、これが直ちに疾患を治療したり移植を改善したりすることを主張するものではありませんが、幹細胞コロニーには、その成長と停止のパターンを決定する強力で遺伝的な構造があることを証明する、新しいモデルフリー（モデルに依存しない）な手法を提供しています。

技術概要

技術要約：細胞系統樹における継承エントロピー

1. 問題提起

細胞系統樹における遺伝的伝達の影響を定量化することは、生物学における根本的な課題であり続けている。シングルセル・プロファイリングは細胞状態の特性評価を可能にするが、これを古典的な系統追跡（lineage tracing）と組み合わせ、遺伝的な影響を追跡することは困難である。特有の課題は、不活性化（老化、分化、または静止期を含む、分裂を停止した細胞）において生じる。活発な増殖においては、遺伝的な変異が子孫の表現型を増幅させるが、不活性化は自己限定的である。すなわち、細胞が不活性になるとその枝は終了し、下流の情報は消失してしまう。このため、不活性化の出現がランダムな確率的事象によるものか、あるいは細胞周期の退出を制御する遺伝的（エピジェネティック）な要因によるものかを判断することが困難となる。さらに、遺伝的な変異とエピジェネティックな継承を区別することは、治療の可能性に影響を与えるコロニー間の不均一性を理解する上で極めて重要である。

2. 手法

著者らは、系統樹のトポロジー（位相構造）にシャノン・エントロピーを適用した、モデルに依存しない手法を導入している。これは特に、不活性細胞の分布に焦点を当てたものである。

A. データ取得

本研究では、32個のヒト骨髄由来間葉系幹細胞（BMSC）クローナル・コロニーの高解像度シングルセル系統追跡を利用している。実験は、すべての複製細胞が第7世代（$k_{max}=7$）に達するまで、位相差タイムラプス顕微鏡を用いて行われた。細胞は以下のように分類された：

• 活性（Active）: 分裂（mitosis）が観察されたもの。
• 不活性（Inactive, G0/Dead）: 84時間分裂が見られなかった、または細胞死に至ったもの。
• リーフ（Leaves）: $k=7$ において誕生のみが追跡された細胞。

B. 不活性不均衡（Inactivity Imbalance）

系統樹における不活性のトポロジー的影響を定量化するために、著者らはすべての分裂 $m$ における不活性不均衡（$I_m$）を次のように定義する：
$$I_m = |M_{m}^{left} - M_{m}^{right}|$$
ここで、$M$ は左側および右側の枝における欠損した $k_{max}$ リーフの数を表す。特定のノードにおける大きな不均衡は、二つの娘細胞の枝の間で不活性の確率に乖離があることを示唆しており、これは上流での変異イベントの可能性を示している。

C. 継承エントロピー（Inheritance Entropy）

核となる指標は、系統内のすべての分裂における不均衡の正規化された重み（$w_m$）から計算される継承エントロピー（$S$）である：
$$w_m = \frac{I_m}{\sum_n I_n}, \quad S = -\sum_m w_m \log w_m$$

• 高エントロピー: 不均衡が樹状構造全体に均質に散らばっており、ランダムで非遺伝的な変動であることを示唆する。
• 低エントロピー: 不均衡が少数のノードに集中しており、特定の系統（サブブランチ）の不活性確率を変化させた特定の遺伝的イベント（変異）を示唆する。

D. ヌル・アンサンブル（Null Ensembles）と検証

生物学的系統樹のエントロピーを、以下の3つの非遺伝的なヌル・アンサンブルと比較して統計的有意性を判定する：

1. ガルトン・ワトソン（GW）モデル: すべての細胞に対して一定の不活性確率（$q$）を持つ。
2. 非均質GW（NHGW）モデル: 世代依存的な不活性確率（$q_k$）を持つ。
3. ランダム・リシャッフルリング（RR）法: 推奨される手法。各世代内の細胞の位置を入れ替え、各世代における不活性細胞の総数は維持しつつ、細胞の配置をかき混ぜる。これにより、母細胞と娘細胞の遺伝的リンクを断ち切りつつ、樹状構造の制約を維持する。

P値は、ヌル・アンサンブルの樹の中で $S \leq S_{biological}$ となるエントロピーを持つ樹の割合として定義される。低いP値は、継承の強い証拠を示す。

E. モデル適合（SKH）

著者らは、データを記述するために**ストカスティック・カークウッド・ホリデー（SKH）**モデルを提案している。このモデルでは：

• 野生型細胞は、確率 $p$ で「変異」状態へと変化する。
• 変異細胞（およびその子孫）は、不活性化する確率 $q$ を持つ。
• 決定的な点として、このモデルは変異ラグ（$l$）、すなわち変異イベントから表現型の発現（不活性化）までの世代数を組み込んでいる。

3. 主な貢献

1. 継承エントロピー: 遮断された（不活性化した）枝のトポロジーのみに基づいて、系統樹における遺伝的構造を検出する、新しいモデルに依存しない指標。
2. ヌルモデルの検証: GWおよびNHGWモデルは、非均質な非遺伝的ツリーと遺伝的ツリーを区別できず（偽陽性を生む）、**ランダム・リシャッフルリング（RR）**法が堅牢で感度が高く、モデルに依存しないことを示す厳密なベンチマーク研究。
3. 変異ラグの定量化: 不活性の確率が変化した特定の分裂箇所を特定し、変異ラグ（変異と表現型の発現の間の遅延）を測定する能力。
4. ストカスティック・カークウッド・ホリデー・モデル: 不活性化への変異ラグを含むBMSC系統の遺伝的特徴を捉えることに成功した、一般化された確率論的モデル。

4. 結果

• 遺伝的構造: 32個のBMSCコロニーの実験データセットにおいて、**21個のコロニー（75%）**がRRテストを用いて有意な遺伝的特性（$P < 0.05$）を示した。残りのコロニーは、不活性細胞が少なすぎて統計的に有意な結果を得られなかったものであり、必ずしも継承の欠如を意味するものではない。
• エントロピー分析: 生物学的な樹は、ランダムにリシャッフルされた対照群よりも一貫して有意に低いエントロピーを示した。これは、不活性細胞の分布がランダムではなく、遺伝的なパターンに従っていることを証明している。
• 変異ラグ: 平均的な変異ラグは約2世代であることが判明した。これは、不活性化につながるエピジェネティックな変化が、実際の表現型の発現よりも上流で起こっていることを裏付けている。
• モデル適合: データに適合させたSKHモデルは、$p \approx 0.1$（変異確率）および $q \approx 0.46$（変異細胞の不活性確率）というパラメータを与えた。これらのパラメータを用いて生成された合成樹は、実験的なBMSC系統のトポロジー構造と酷似していた。
• 空間的 vs トポロジー的: 不活性細胞のクラスター化は、培養皿上の物理的な近接性（空間的相関）によるアーティファクトではなく、系統のトポロジー的特徴であることが分析により確認された。

5. 意義と主張

本論文は、観察されたヒトBMSCコロニーの不均一性は、純粋に確率的なメカニズムや遺伝的なメカニズムではなく、細胞周期の退出を制御するエピジェネティックな遺伝的要因によって駆動されていると主張している。

• 堅牢性: エントロピーに基づくテストは、あらゆる系統樹における差異のある成長の遺伝的起源を評価するための、単純で堅牢かつモデルに依存しない方法を提供する。
• 臨床的関連性: これらの知見は、BMSC移植で見られるコロニー間の不均一性（骨再生治療における大きな障壁）が、エピジェネティックな遺伝的起源を持っていることを示唆している。これらの遺伝的メカニズムを理解し制御することは、幹細胞治療の一貫性を向上させる可能性がある。
• 手法の進歩: 「変異ラグ」を確立することで、表現型（不活性）が自己消去的な性質を持つ場合でも、研究者は表現型の変化の起源を遡ることが可能になる。
• 今後の方向性: 著者らは、現在の手法は不活性化に焦点を当てているが、系統のトポロジーと細胞周期時間の相関関係は、細胞周期の持続時間もエピジェネティックな遺伝的メカニズムによって制御されている可能性を示唆しているものの、これにはさらなる研究が必要であると述べている。

著者らは、SKHモデルは簡略化されたものであること、また、より複雑なモデルが変異ラグの分散をより良く捉えられる可能性があることを認め、控えめな姿勢を保っている。彼らは、この手法は確率の変化を検出するものであることを強調している。つまり、系統に不活性の発現に関する変動がない場合、たとえ遺伝的なメカニズムが存在したとしても、このテストは継承を検出できない。