- 著者: Saori Goto, Vikram Deshpande, Ömer H. Yilmaz, Norihiro Goto
- Corresponding author: Ömer H. Yilmaz (Koch Institute, MIT / Broad Institute / MGH; ohyilmaz@gmail.com); Norihiro Goto (Weill Cornell Medicine, Cornell University; nog4004@med.cornell.edu)
- 雑誌: Cell Stem Cell (Vol. 33, in press)
- 発行年: 2026
- Epub日: 2026-06-05
- Article種別: Original Article
- PMID: 42330962
背景
大腸癌 (CRC; colorectal cancer) は癌関連死の主要因の一つであり、転移が CRC 関連死の大半を占め、肝が遠隔転移の優勢部位である。化学療法や免疫チェックポイント阻害の進歩で一部の CRC 患者の予後は改善したが、肝転移を有する患者が現行治療で完全寛解に至ることは稀である。先行研究では、matched 原発腫瘍と肝転移の whole-exome sequencing 比較が pro-metastatic driver 変異を同定できなかった一方 (先行研究)、single-cell RNA sequencing (scRNA-seq) 研究は原発と肝転移で明瞭に異なる遺伝子発現プロファイルを示し、エピジェネティック変化が肝転移促進に critical な役割を果たす可能性を示唆してきた。しかし腫瘍遺伝学・微小環境・患者背景の著しい異質性が「driver」エピジェネティック変化を偶発的「passenger」変化から区別することを困難にしていた。CRC には健常結腸と同様に LGR5+ stem-like 細胞と LGR5- 細胞の階層が存在し、LGR5+ 細胞が腫瘍増殖維持に必須である一方、Xu et al. 2022 や Ganesh et al. 2020 の系譜追跡研究 (先行研究) が LGR5- 細胞の肝転移における critical な役割を示した。EMP1 で標識される LGR5- 細胞の亜集団は原発切除後の転移再発を駆動し、basal-like/squamous 転写状態を示す。fetal-like reversion は PROX1・SOX17・YAP/TAZ など複数経路で誘導され、ATRX と HNF4A の喪失が結腸系譜同一性を破壊して basal-like/squamous 可塑性を促すことも報告された。既存研究では、これら fetal-like・basal-like/squamous という multilineage program を束ねる統一的制御因子が不足しており、転移性 CRC を理解・標的化するうえでの gap in knowledge であった。著者らは先行して SOX17 が早期免疫逃避を起動することを示しており (Goto et al. Nature 2024)、本研究はこの gap を、肝転移由来 organoid の連続同所移植による時系列エピジェネティック解析で埋める。
目的
本研究は、① 肝転移由来 organoid の連続同所移植により pro-metastatic CRC organoid ライブラリを樹立し、肝転移進展中に生じるエピジェネティック・転写変化の時系列軌跡を捉えること、② RNA-seq と ATAC-seq の統合により「driver」となる pro-metastatic signature を同定し、これを束ねる中心的転写因子を motif/TF 活性解析で特定すること、③ その因子 (GATA6) の喪失が肝転移を促進する機構を loss/gain-of-function 実験と CUT&RUN で解明し、既報の SOX17 経路との関係を明らかにすること、を目的とする。
結果
連続移植による pro-metastatic CRC organoid ライブラリの樹立：Apc・Kras[G12D]・P53 変異を導入したマウス結腸癌 organoid (naive AKP organoids、AKP = Apc/Kras/P53 変異型) を免疫正常マウス結腸に同所移植すると、マウスは原発腫瘍による腸閉塞で生存中央値 78.5 日となり、約 30% のみが肝転移を発症し (1 匹あたり平均 0.8 個)、肝転移の低頻度が解析の障害であった (Fig 1B-E、n=27 匹)。肝転移由来 organoid を再移植する操作を 3 ラウンド反復し、原発腫瘍由来 organoid を対照とした。各ラウンドで生存時間 (= 原発腫瘍増殖の指標) は両群で差がなかったが、肝転移由来 organoid 移植群は肝転移マウスの割合と 1 匹あたり転移数が有意に高かった (Fig 1F-H、n=5-10 匹/群)。TGF-β 処理で 3° 肝転移 organoid の増殖が抑制されたことから TGF-β シグナルは保存され、転移能亢進が TGF-β 経路の遺伝子変化に依らないことを確認した (Fig S1B)。
RNA-seq/ATAC-seq 統合による pro-metastatic signature の同定：3° 肝転移 organoid を naive AKP と比較した differential gene expression (DGE) 解析では、有意に 1,624 遺伝子が上昇し 1,258 遺伝子が低下した (fold change >2 または <0.5、adjusted p<0.05、Fig 2A)。Gene Ontology では “tissue development” が濃縮した。GSEA は fetal-like intestine program (Clu, Spp1, Aqp5) の上昇と、LGR5+ 幹細胞シグネチャ (Lgr5, Ascl2, Smoc2) の低下を示した。さらに regenerative cell signature・EpiHR (epithelial-specific high-risk)・basal-like/squamous pancreatic cancer signature も濃縮した (Fig 2D、NES 最大 1.86、FDR 0.0008)。ATAC-seq promoter accessibility では 577 peak が増加し 417 peak が減少し (Fig S1H)、移植ラウンドと ATAC シグナルは正相関した。次に RNA-seq と ATAC-seq を統合し、247 の共通上昇遺伝子 (Tff1, Chl1 等) と 213 の共通低下遺伝子 (Pctp, Gata6 等) を “pro-metastatic signature” と定義した。この signature は ssGSEA で移植ラウンドと正相関し (Spearman)、ヒト CRC 肝転移 scRNA-seq でも保存された (Fig 2G-I、Fig S2C,D)。
pro-metastatic エピジェネティック変化は GATA6 喪失に収束する：HOMER による motif enrichment で 3° 肝転移 organoid は AP-1 family motif の accessibility 増加と GATA family (Gata1/Gata2/Gata6)・HNF・TCF/LEF motif の減少を示し、TOBIAS footprinting でも GATA family (Gata2/3/5/6) の occupancy 低下が確認された (Fig 3A、Fig S2G)。motif・promoter accessibility・遺伝子発現の 3 層を統合すると、GATA6 が全解析に共通する唯一の転写因子として同定された (Fig 3B)。Gata6 promoter の active mark H3K4me3 accessibility は連続移植で漸減し Gata6 発現も低下した (Fig 3C,D)。臨床的関連として、ヒト CRC の免疫染色で GATA6 は原発腫瘍の大半で high/mid (≥30% 陽性、82.1%、32/39 例) だが肝転移では大半が low/negative (95%、34/36 例) であった (Fig 3F,G、n=39 原発 / n=36 肝転移)。paired scRNA-seq でも肝転移上皮で GATA6 が有意低下し (paired two-tailed t 検定、Fig 3H、n=24)、原発の GATA6 発現は matched 肝転移と正相関した (Spearman、Fig 3I)。TCGA COAD/READ で stage IV 原発 CRC は GATA6 低発現で、低 GATA6 は全生存・無病生存不良と関連した (log-rank 検定、Fig 3J,K)。
GATA6 欠失は原発腫瘍に影響せず肝転移を促進する：naive AKP organoid で CRISPR-Cas9 により GATA6 を欠失させると (immunoblot で確認、Fig 4B)、対照の spheroid 構造に対し crypt-like fold と緩徐な増殖を示したが、GATA6-null adenoma と異なり CRC では自己複製能は保たれ 12 継代後も生存した (Fig 4C,D、Fig S3B,C)。同所移植で対照と GATA6-null はともに中分化浸潤腺癌の原発腫瘍を形成し、原発腫瘍サイズ・生存時間に差はなかったが、肝転移数は GATA6-null で有意に増加した (unpaired two-tailed t 検定、Fig 4J-M、n=6-9 匹/群)。脾臓内注入 (spleno-portal 循環経由) でも GATA6-null 群で 6 週後の肝転移数が有意増加し (Fig 4N-P、n=11 匹)、逆に 3° 肝転移 organoid への GATA6 再発現は肝転移を有意に減少させた (gain-of-function、Fig S3J,K)。GATA6-null organoid の RNA-seq は 841 遺伝子上昇・566 遺伝子低下を示し (Fig 5B)、pro-metastatic signature の強い濃縮、Hnf4a 低下と fetal-like marker (Psca, Glis3, Il33) 上昇を認め、flow cytometry で LY6D+ basal-like/squamous 細胞が対照 ~0.6% から GATA6-null ~5% へ拡大した (Fig 5I,J)。In situ hybridization と RT-qPCR で GATA6-null organoid は Lgr5 mRNA が約 200-fold 低下し LGR5- 状態へ移行した (Fig 5K,L)。
HNF4A 直接抑制と H3K27ac 再構築を介した SOX17 非依存の系譜可塑性：GATA6 の CUT&RUN で 895 の consensus GATA6 peak を同定した。GATA6-bound かつ発現変動する 155 遺伝子のうち 93 が GATA6 loss で低下し (GATA6 は activator/repressor 両機能)、その中に basal-like/squamous を駆動する HNF4A が含まれた (Fig 6E,F)。実際、GATA6-null organoid への HNF4A 再発現は basal-like/squamous marker を低下させ、HNF4A 抑制が basal-like/squamous 再プログラム化に必要であることが示された (Fig 6G-I)。次に、早期 CRC で fetal-like reversion を誘導する SOX17 との関係を検討した。GATA6 は SOX17 制御領域に直接結合せず、SOX17 発現に有意な影響を与えなかった (Fig 6J)。一方で SOX17-bound 遺伝子の約 1/4 が GATA6-bound と重複し、co-immunoprecipitation では SOX17 と GATA6 が同一複合体に存在した (Fig 7C-E)。GATA6-null での SOX17 過剰発現は Lgr5 低下・Prox1 上昇を相加的に増強しなかった。すなわち GATA6 喪失は SOX17 過剰発現 (SOX17-OE) の誘導プログラムを largely recapitulate し、追加の SOX17 活性を不要にする (SOX17 非依存)。最後に H3K27ac CUT&RUN で、GATA6-null は 573 peak の H3K27ac 有意増加を示した。Prox1 等 fetal-like program entry 遺伝子近傍で H3K27ac 亢進を認め、GATA6 喪失が直接標的を超えた広範な H3K27ac 制御領域の再構築を誘導することが示された (Fig 7F-K)。
考察/結論
本研究は、肝転移由来 organoid の連続同所移植という新規モデルで、CRC 肝転移を駆動する統一的エピジェネティック制御因子として GATA6 喪失を同定した。先行研究との違いとして、matched 原発・肝転移の whole-exome sequencing が pro-metastatic driver 変異を同定できなかったこれまでのアプローチと対照的に、本研究は driver エピジェネティック変化を passenger 変化から区別するために「共通の初発 oncogenic driver を持ち、ラウンドごとに転移能が漸増する isogenic organoid ライブラリ」を用いる設計で、背景異質性の交絡を排除した点が異なる。新規性として、fetal-like reversion と basal-like/squamous reprogramming という 2 つの multilineage program を束ねる中心因子が GATA6 喪失であることを本研究で初めて示し、GATA6 が HNF4A を直接抑制し H3K27ac を再構築して LGR5 陰性化を伴う可塑性を駆動するという機構を提示した。さらに、GATA6 喪失による fetal-like 誘導が既報の SOX17 経路とは独立 (SOX17-independent) に作動し、GATA6 が SOX17-OE のプログラムを largely recapitulate するという上流・下流関係を明らかにした点は、著者ら自身の SOX17 研究 (Goto et al. Nature 2024) を発展させる新規知見である。この GATA6-LGR5 軸は、LGR5- 細胞が転移を担うという知見や、早期免疫逃避と CSC 可塑性を統合した総説の枠組みとも整合する (Goto et al. TrendsCancer 2025)。TGF-β 活性化 CAF による免疫逃避との連関も CRC 転移の理解に重要である (Xie et al. NatCommun 2022)。臨床応用の観点では、GATA6 発現がヒト CRC 原発腫瘍で肝転移リスクと予後を反映する (低 GATA6 = stage IV・予後不良) ことから、GATA6 が原発腫瘍生検段階での転移リスク層別化バイオマーカーとなりうる。また GATA6 喪失下流の HNF4A 抑制や H3K27ac 再構築、basal-like/squamous 状態は、bench-to-bedside な治療標的 (エピジェネティック修飾酵素・系譜可塑性阻害) の候補を提示する。残された課題として、GATA6 喪失を治療的に回復・代償する方法、GATA6-null 転移細胞の免疫微小環境との相互作用、ヒト CRC における GATA6 低発現の予測的意義の前向き検証が今後の検討として求められ、GATA6 依存的な可塑性制御ネットワークの全容解明が転移予防戦略の鍵となる。
方法
マウス正常結腸 organoid に Apc (sgApc + Cas9)・Kras[G12D]・P53f/f・Rosa-LSL-tdTomato を導入した naive AKP organoid を作製し、C57BL/6 系統の免疫正常マウス (immunocompetent syngeneic mice) の結腸に同所移植 (orthotopic transplantation) した。肝転移・原発腫瘍から organoid を樹立し 3 ラウンド連続移植して pro-metastatic organoid ライブラリを構築した。転写・クロマチン解析には RNA sequencing (RNA-seq)、assay for transposase-accessible chromatin using sequencing (ATAC-seq)、H3K4me3・H3K27ac・GATA6・SOX17 の CUT&RUN (cleavage under targets and release using nuclease) を用い、DGE は fold change >2 または <0.5・adjusted p<0.05 を有意とした。転写因子解析は HOMER (hypergeometric optimization of motif enrichment) による motif enrichment と TOBIAS (transcription factor occupancy prediction by investigation of ATAC-seq signal) による footprinting、機能濃縮は Gene Ontology・GSEA・ssGSEA・TissueEnrich を用いた。GATA6 の機能検証は CRISPR-Cas9 による GATA6-KO、doxycycline 誘導性 shGata6 knockdown、GATA6/HNF4A/SOX17 の過剰発現 (再発現) organoid を作製し、免疫正常マウスへの同所移植および脾臓内注入 (intrasplenic injection、spleno-portal 循環経由、6 週後評価) で肝転移を評価した。検証法は immunoblot・RT-qPCR・免疫染色・in situ hybridization・flow cytometry・co-immunoprecipitation。ヒト検体は CRC 原発腫瘍 (n=39) と肝転移 (n=36) の GATA6 免疫染色、公開 scRNA-seq (Moorman et al. 2025、matched primary-liver、n=24)、TCGA COAD/READ RNA-seq を再解析した。統計は log-rank 検定 (Benjamini-Hochberg 補正)・chi-squared 検定・Mann-Whitney 検定・one-way ANOVA・unpaired/paired two-tailed t 検定・Spearman rank correlation を用い、データは mean ± SD、p<0.05 を有意とした。