Reconstructing whole-organism cell phylogenies with resolved ancestral transcriptional states

技術要約：祖先的転写状態を解明した全個体細胞系統樹の再構築

問題提起
細胞系統追跡と単一細胞RNAシーケンシング（scRNA-seq）の統合は、発生の軌跡をマッピングするための系統樹の再構築を可能にしてきた。しかし、決定的な限界が依然として存在する。それは、内部ノード（祖先細胞）はサンプリングされないため、転写的に「沈黙」しているという点である。末端ノード（サンプリングされた細胞）のトランスクリプトームは利用可能であるが、内部ノードの祖先的転写データが欠如していることは、真の意味での縦断的な細胞状態の動態追跡を妨げ、特定の細胞状態が時間とともにどのように進化するかを観察することを困難にしている。この問題は、不変の系統を持つC. elegansにおいては解決されているが、これをより複雑な非線虫生物へとスケールアップすることは、重大な技術的および計算的な課題となっている。

手法
著者らは、ゼブラフィッシュ（Danio rerio）を用いた実験的革新と計算フレームワークの開発を組み合わせて、この課題に取り組んでいる。

1. 実験系 (scSMALT): 本研究では、SMALT（Single-cell Mutational Lineage Tracing）システムを利用し、これを適応させている。このシステムには、普遍的に発現するGFPの3'-UTR内に2つの1-kb DNAバーコード（HMF1およびHMF2）を持つ、トランスジェニック・ゼブラフィッシュ系統が含まれる。1細胞段階で塩基編集装置（AID、iSceI、およびUGIドメインからなる）を注入し、細胞分裂中にこれらのバーコード内で継続的なC·G-to-T·A変異を誘導する。

   • 統合: このシステムはscRNA-seq（10x Genomics）と結合されている。フルレングスのcDNAライブラリから1-kbのバーコードを標的として増幅する特定のワークフローにより、同一の単一細胞から内因性のトランスクリプトームと系統バーコードを同時に捕捉することが可能となる。
   • データ生成: 著者らは、受精後7日目（7 dpf）の15個体のゼブラフィッシュ幼魚から高解像度のデータを生成し、9,000個以上の細胞からバーコードを回収し、数千の変異を特定して系統樹を再構築した。

2. 計算フレームワーク (LEAP): 未観測の内部ノードにおける転写状態を推論するために、著者らは**LUGエンコード祖先投影（LEAP: LUG-encoded Ancestral Projection）**を開発した。

   • LUGの特定: この手法は、**系統コミット型アップレギュレート遺伝子（LUGs: Lineage-committed Upregulated Genes）**を特定する。これらは、特定のクレード内で一度アップレギュレートされると、その後のすべての末裔系統においてもアップレギュレートされ続ける遺伝子である。
   • インピュテーション（補完）: 系統樹全体にわたってこれらのLUGの発現を二値化することにより、著者らは内部ノードのための「LUGエンコード転写プロファイル」を構築する。これらのプロファイルは、時間分解能を持つリファレンスscRNA-seqアトラス（3 hpfから120 hpfまで）に投影され、内部ノードの細胞アイデンティティおよび発生段階をインピュテーションする。
   • 検証: フレームワークは、まず系統が解明されたC. elegansのグラウンドトゥルース・アトラスを用いてベンチマークされ、創始者系統および細胞型を内部ノードに割り当てる高い精度が示された。

主な結果

• 高解像度系統樹: 本研究は、30個の受精卵から幼魚に至る発生細胞系統樹（2種類のバーコード × 15個体）の再構築に成功した。これらの樹形図は、平均764個の末端ノードを含み、約62.1%の祖先細胞分裂イベントを解明し、高い統計的支持（中央値の超高速ブートストラップ値96%）を得た。
• 祖先的インピュテーションの成功: ゼブラフィッシュのデータにLEAPを適用することで、7,708個の内部ノードの転写状態をインピュテーションすることができた。インピュテーションされた状態は、リファレンスアトラスと強い一致を示し（iLISIスコア ~1.8）、発生の時系列的論理（深いノードほど初期の発生段階に対応する）に従っていた。
• 機能的検証: 推論された祖先状態を検証するため、著者らはインピュテーションされた「ニューロン」集団に特異的な5つのLUGに対してCRISPR-Cas9ノックアウトを行った。その結果得られたF0クリスパント（crispants）のscRNA-seqにより、神経細胞の存在量が有意に減少していることが確認され、LEAPによる細胞アイデンティティの正確性が支持された。
• 発生の波の発見: 系統情報を活用した縦断的解析により、孵化から幼魚への移行期（>72 hpf）に関連する、これまで認識されていなかった主要な運命特異化の波が明らかになった。この波は、よく知られた原腸形成のイベントとは別物であり、多能な「前駆細胞」の収縮と、運命が制限された「表皮（periderm）」集団の拡大を伴う組成変化によって駆動されていた。
• 潜在的な運命決定: 解析により、「隠れた層」としての発生の不均一性が明らかになった。具体的には、インピュテーションされた「前駆細胞」の集団（転写的には類似している）が、系統学的な状態においては不均一であることが判明した。これらのサブ集団は、最小限の転写的差異にもかかわらず、特定の末端運命への明確なバイアスを示しており、これは運命決定が顕著な転写的分岐に先行して起こり得ることを示唆している。

意義と主張
本論文は、再構築された細胞系統樹において、非線虫生物の祖先的転写状態を解明した最初の研究であると主張している。内部ノードの「転写的沈黙」を克服することで、本研究は細胞状態のダイナミクスに関する系統情報を伴った縦断的解析を可能にした。

著者らは、このアプローチが複雑な生物における系統解明された細胞アトラスの構築への道を開くと考えている。本研究の知見は、細胞系統が、状態ベースのアトラスに対する補完的なフレームワークとして機能することを示しており、標準的な転写学的クラスタリングだけでは不可視である、細胞の多様化における歴史的および機能的な次元を明らかにできることを示している。具体的には、系統解析が、転写的な差異が最小限である状況下でも、潜在的な運命決定を検出できることを実証しており、発生の軌跡に関するより包括的な視点を提供している。

認められた限界
著者らは、以下の限界についても謙虚に述べている。

• 現在の個別の系統樹の規模は、全系統的状態を解明するための統計的パワーを制限している。
• 現在のデータでは、特定された潜在的な運命バイアスが、固有の分子プログラム、空間的パターン形成、または環境的要因のいずれかに由来するかを区別することはできない（空間情報が測定されていないため）。
• scSMALTを他のモダリティ（例：scATAC-seq）と統合することで、系統的状態がエピジェネティックなプライミングと結合しているかどうかを判断する必要がある。
• バーコードの発現キネティクスや細胞数の影響により、極めて初期の段階（例：原腸形成期）における祖先状態の解明は依然として困難である。