Low-heteroplasmy mitochondrial DNA mutations improve clonal reconstruction of human cells

技術要約：低ヘテロプラスミーのミトコンドリアDNA変異によるヒト細胞のクローン再構成の向上

問題提起
ヒト組織におけるクローンの関係性を再構成することは、発生、老化、および疾患を理解する上で不可欠である。ソマティックなミトコンドリアDNA（mtDNA）変異は、内因性のシングルセル・バーコードとして機能するが、従来の系統追跡は高ヘテロプラスミー変異（通常10%超）への依存により制限されてきた。これらの高頻度変異は検出が容易である一方、数値的には乏しく、選択によって形成されている可能性があるため、利用可能な系統情報のわずかな一部しか捉えられない。対照的に、より豊富な低ヘテロプラスミー変異（10%未満）は、シングルセル・シーケンシングにおける技術的ノイズと区別することが困難であり、そのために系統解析から除外されてきた。低ヘテロプラスミーの観察結果が、真の系統履歴をコードしているのか、それとも単なるアーティファクトであるのかは、これまで検証されてこなかった。

手法
著者らは、低ヘテロプラスミーのmtDNA変異の有用性を評価するために、グラウンドトゥルース（正解）を用いたベンチマーク・アプローチを利用した。

• 実験系： 健康なドナー由来の一次ヒト造血幹・前駆細胞（HSPC）を用い、レンチウイルスによる系統およびRNA回収バーコード（LARRY v2）ライブラリを導入した。これらの細胞をex vivoで増殖させ、mtDNAおよびレンチウイルス・バーコードの標的キャプチャを用いたシングルセル・マルチオーム・シーケンシングを用いてプロファイリングした。レンチウイルス・バーコードは、クローンの同一性を定義する独立したグラウンドトゥルース（lenti-clones）を提供した。
• 変異コーリングとフィルタリング： 本研究では、mtDNA変異プロファイリングのためにReDeeMパイプラインを採用した。著者らは、オリジナルの戦略（filter1）に基づいた更新されたポストコンセンサス・フィルタリング・ワークフロー（filter2と呼称）であるredeemR2.0を導入した。filter2の主要な構成要素は以下の通りである：
  • エッジに偏ったアーティファクトを除去するための9-bpフラグメント端トリミング。
  • 細胞間の分布がランダムな技術的ノイズから逸脱している変異を保持するための、変異ごとの二項適合度検定（FDR 0.05）。
  • 集団頻度（MITOMAP）に基づく生殖細胞系列様の変異の除去、およびホモプラスミーまたはブラックリスト領域の除外。
  • 統計的フィルターを通過する限り、1つの細胞あたりわずか1つの内因性ユニーク分子識別子（eUMI）によって支持される変異の保持。
• 系統推論： 著者らは、ミトコンドリアの遺伝とヘテロプラスミーのドリフトをモデル化する信頼度を考慮したフレームワークであるMitoDriftを用い、フィルタリングされたmtDNA変異プロファイルから系統樹を再構成した。
• ベンチマーク指標： パフォーマンスは、推論されたクレイド（系統群）をグラウンドトゥルースであるlenti-clonesと比較することで評価された。指標には以下が含まれる：
  • 再現率（Recall）： 推論されたクレイドによって回収されたグラウンドトゥルース・クローンの割合。
  • 適合率（Precision）： グラウンドトゥルースのクローンと一致した推論されたクレイドの割合。
  • ジャカード類似度（Jaccard Similarity）： 推論されたクレイドとグラウンドトゥルースのクローンの間のマッチングを定義するための閾値（≥0.5および≥0.75）。
  • 精度・再現率曲線の下端面積（AUPRC）： 様々なMitoDriftサポート閾値にわたって算出。

主な結果

• 低ヘテロプラスミー・コールの検証： 厳格なfilter2ワークフローを適用した後、単一のeUMIによって支持される、あるいは10%未満のヘテロプラスミーに該当する変異コールは、真のmtDNA変異の分子学的特徴を保持していた。これらには、遷移優位のシグネチャー（トランスバージョン比 ~0.03）とリード断片に沿った一様な分布が含まれ、背景ノイズと区別された。
• 改善されたクローン回収： LARRYのグラウンドトゥルースを用いた結果、本研究は、全ヘテロプラスミー・スペクトラムを保持することが、制限的なカットオフを用いるよりも大幅に優れていることを示した。
  • 10%超のヘテロプラスミー・カットオフと比較して、全filter2セットを使用した場合、Jaccard閾値0.5においてAUPRCは0.322から0.640へ（2倍の向上）、Jaccard閾値0.75においてAUPRCは0.203から0.433へ（2.1倍の向上）増加した。
  • 1-eUMI変異を含めることでさらに性能が向上し、0.5の閾値においてAUPRCは0.573から0.640へと上昇した。
• 系統情報の分布： 個々の変異コールに対する分析により、低ヘテロプラスミー変異は無視できないノイズではなく、実質的なシグナルの源であることが明らかになった。
  • 系統情報を持つコール（特定のクローン内で濃縮されているもの）の中で、52.9%が細胞あたりのヘテロプラスミー10%未満であり、42.0%が細胞あたりわずか1〜3つのeUMIによって支持されていた。
  • 1-eUMIグループおよび最低ヘテロプラスミーグループは、高サポートグループよりもノイズが多いものの、ランダムな背景と比較して系統情報を持つコールが有意に濃縮されており、クローン回収における集団的な価値が確認された。
• フィルタリング戦略： redeemR2.0/filter2戦略は、UMIベースのエラー補正を欠くmgatkのような設定やfilter1を含む他のすべてのテストされたアプローチを、精度・再現性能において上回った。

意義および主張
本論文は、低ヘテロプラスミーのmtDNA変異が、豊富かつ真正なクローン履歴の記録を構成していることを確立している。著者らは以下の結論を得ている：

1. シグナルの拡大： 制限的なヘテロプラスミー閾値（例：10%超）は、大幅な感度のコストを伴う。この閾値を下回る変異、たとえ分子的な支持が低いものであっても、それらを含めることは、利用可能な系統マーカーのプールを大幅に拡大させる。
2. 集団的有用性： 個々の低分子コールは注意深く解釈されるべきであるが、それらの集団的な包含は、真のクローン関係の回収を大幅に改善する。
3. ベストプラクティス： 著者らは、redeemR2.0/filter2を、mtDNAベースのシングルセル系統追跡におけるデフォルトのプリプロセッシング戦略として推奨する。このアプローチは、厳格な分子レベルのフィルタリングと全ヘテロプラスミー・スペクトラムの包含を組み合わせることで、遺伝子工学を必要とせずにヒト組織における高解像度のクローン再構成を可能にする。
4. フレームワーク： 本研究は、グラウンドトゥルースにアンカリングされたベストプラクティスのフレームワークを提供しており、適切なエラー補正（eUMIコンセンサスおよび統計的フィルタリング）を用いることで、「ノイズ」としての低ヘテロプラスミー変異を高価値な系統データへと変換できることを示している。

著者らは、グラウンドトゥルース（レンチウイルス・バーコーディング）がバーコードのサイレンシングやドロップアウトにより不完全である可能性があることに留意しているが、観察された系統回収の利得は堅牢であるとしている。彼らは、特定のパラメータは細胞型によって異なる可能性があるものの、厳格なフィルタリングを伴う低ヘテロプラスミー変異を活用するという原則は、ヒト細胞の歴史における自然なミトコンドリア変異の解像度を高めるために広く適用可能であることを強調している。