背景

がんの進行は、遺伝的に異なるがん細胞集団が選択される過程であり、同時に低酸素・免疫抑制・線維化が蓄積される生態系の変化でもある。しかし、単一細胞RNAシーケンシングは組織を解離するため、細胞の元の位置が失われ、マルチゾーンシーケンシングは限られたポイントのみを採取するため、細胞状態と系譜の連続的な空間構造を失いやすい。がん細胞の亜克隆がどこで成長し、周囲の細胞をどのように変化させるのか、転移能力を得た細胞が原発腫瘍のどの領域から出現するのかは、まだ十分に解明されていない。

研究チームはこのギャップを埋めるために、高解像度空間トランスクリプトームと時間に応じて遺伝的マーカーが蓄積される系譜追跡を組み合わせた。対象はKras活性化とTrp53欠損により肺腺癌を発症するKP-Tracerマウスである。8〜12週齢マウスの肺胞2型細胞にCreウイルスを投与し、腫瘍発生とCRISPR–Cas9による系譜記録を同時に開始した。がん細胞の現在の転写状態だけでなく、共通祖先と増殖歴を組織座標上で復元することを目的とした。

主な発見

研究チームは、広い組織面積を準細胞レベルで読み取るSlide-seqと、単一核レベルの感度を提供するSlide-tagsを併用した。合計49個の空間トランスクリプトームアレイから100個以上の腫瘍を分析し、系譜マーカーが欠落した位置には周囲のがん細胞の空間・遺伝情報を利用した補正アルゴリズムを適用した。実際のデータから隠されたマーカーを復元した中央精度は90%であり、繰り返し補正によりデータセットごとの欠損情報の4〜58%、平均31%を回復した。ランダム予測の中央精度は67%であった。

腫瘍が進行するにつれて、肺胞様状態を維持したがん細胞は主に縁に、上皮間葉転換（EMT）と低酸素遺伝子プログラムを持つ細胞は内部に位置した。急速に拡大した高適応度亜克隆の周囲には、Arg1陽性の免疫抑制性腫瘍関連マクロファージとPostn陽性の筋線維芽細胞様がん関連線維芽細胞が集まった。系譜樹から計算した適応度は、従来のトランスクリプトームベースの適応度指標とピアソン相関係数0.4を示し、低酸素・EMT空間コミュニティが最も高い適応度と関連していた。これは亜克隆増殖、酸素欠乏、免疫・基質再編、転移促進性細胞状態が順次関与している説明を裏付ける。空間解析は、Piezo2・Robo1・Pecam1を発現しながらNkx2-1を維持しVimは発現しない未報告のがん細胞状態も発見した。

転移が広範囲に起こったマウスでは、肺を200〜500マイクロメートルの間隔で切断し、原発腫瘍を三次元的に追跡した。中隔リンパ節・肋骨・横隔膜の転移病巣はすべて、原発腫瘍T2内の空間的に限定された亜克隆と系譜的に結びついていた。転移病巣のTGF-βプログラムは該当の原発亜克隆とほぼ同じであったが（log2倍数変化−0.14、P=1.0）、コラーゲン関連発現は大幅に増加した（log2倍数変化3.81、P<10⁻⁵）。COL3A1タンパク質染色も転移部位の線維化を裏付けた。Nature Genetics論文

意義と展望

今回の研究は、攻撃的ながん細胞状態が細胞内変異だけで固定されるのではなく、亜克隆拡大によって生じた低酸素・免疫抑制・線維化環境と共に形成される空間的根拠を提示した。転移の種が原発腫瘍全体に均等に広がっているのではなく、特定の局所生態系に集中しているという結果も注目される。転移後には、既存のTGF-β・EMTプログラムが維持される一方でコラーゲン沈着が加わることで、転移定着に適した新たな環境が形成される。

ただし、結果はKras;Trp53ベースのマウス肺腺癌モデルに主に依存している。薄い組織切片は腫瘍の全体的な三次元構造を代表できない可能性があり、空間データの系譜マーカー多様性も解離型単一細胞解析より低い。欠落マーカーを近隣細胞で補正する方法は、移動性の高いがん細胞で誤差が生じる可能性もある。人間の肺腺癌コホートと治療前後の検体で同じ空間構造が再現されるかどうか、低酸素や線維化シグナルを遮断した場合に亜克隆の転移能力が実際に低下するかどうかを検証する必要がある。