- 著者: Shiman Hu, Jiaxian Liu, Jiaqi Zhou, Jiaohua Chen, Tiantian Qin, Yi Ching Esther Wan, Xiaoxuan Zhu, Danyi Wang, Chuting Shao, Yabin Chen, Xin Wang, Junhong Han, Hoi Leong Xavier Wong, Robert S. Weiss, Haojie Jin, Mo Chen, Qing Li, Yogen Saunthararajah, Haiyun Gan, Kui Ming Chan
- Corresponding author: Kui Ming Chan (City University of Hong Kong)
- 雑誌: Science Advances
- 発行年: 2026
- Epub日: 2026-07-10
- Article種別: Original Research Article (basic science, molecular/epigenetic oncology)
- DOI: 10.1126/sciadv.adx4982
背景
がんゲノムにおけるヒストン変異の発がん的役割は、小児脳腫瘍における H3K27M 変異の発見 (ヒストン H3.3 の第 27 番リジンのメチオニン置換) を端緒として急速に解明されてきた。こうした変異ヒストンは「オンコヒストン (oncohistone)」と総称され、ヒストン修飾パターンをゲノムワイドに改変し転写プログラムを根底から書き換える。H2B サブタイプにおいても G53D 変異 (膵臓癌)、E76K 変異 (乳癌)、E113K 変異 (複数がん種) など複数の機能獲得型変異が同定されており、Wan らの包括的解析でヒストン H2B 変異の発がん的意義が注目されるようになった。
ヌクレオソームのアシッドパッチ (nucleosome acidic patch) は H2A (E56, E61, E64, D90, E91, E92) と H2B (E105, E113) の計 8 残基から構成されるヌクレオソーム表面の酸性クラスターであり、Imitation Switch (ISWI) 型クロマチンリモデリング複合体やその他クロマチン調節因子の結合足場として機能する。特に SMARCA5 は ATP 依存的にヌクレオソームを転写開始点 (TSS: transcription start site) 付近でスライドさせることで転写因子結合部位を開くことが示されており (Hargreaves et al. NatGenet 2021)、アシッドパッチの電荷状態がこの反応の効率を制御する。H2BE113K 変異は中性リジンによって E113 の負電荷が喪失するため、アシッドパッチ依存的なクロマチンリモデリングの基盤を変容させる可能性が示唆されていた。
乳癌は 2022 年に世界で約 230 万人が新規診断され、約 67 万人が死亡した最も患者数の多い固形がんである。転移が最大の死亡原因であり、とりわけ三重陰性乳癌 (TNBC) ではサブタイプ多様性と分子標的治療の制限から予後改善が課題となっている (Yan et al. Nature 2026)。しかし、ヒストン変異が乳癌の転移能に及ぼす因果的影響と機序は従来ほとんど解明されておらず、H2BE113K がどのようなクロマチン変動を経て転移表現型を促進するかは未解明であった。また、標的治療可能な下流分子経路の同定に関する知見は不足しており (Gilan et al. Science 2020)、このギャップを埋めることが本研究の臨床的動機となった。
目的
乳癌に見出されるヒストン H2B の第 113 番グルタミン酸 (E) のリジン (K) への置換変異 (H2BE113K) が、クロマチン動態の変容を介して乳癌の転移能を促進するか否かを CRISPR-Cas9 ノックイン細胞モデルおよびマウスモデルで検討し、その分子機序と治療的介入の可能性を明らかにする。
結果
H2BE113K変異の乳癌における分布と機能的局在： The Cancer Genome Atlas (TCGA) および Memorial Sloan Kettering (MSK) データベースの横断解析において、H2BE113K 変異は乳癌全体の約 0.6% に存在し、そのうち 6 例が浸潤性乳管癌であった。同変異は肺癌、頭頸部癌、膀胱癌にも散在したが、乳癌における頻度が最も高かった。クロマチン分画実験では、FLAG タグを付与した E113K 変異体 H2B が WT H2B と同様にクロマチン結合分画に局在し、ヌクレオソームへの正常な組み込みを示した。ヌクレオソーム結晶構造に基づく三次元マッピングから、E113 残基はアシッドパッチの 8 残基中の 1 つとして ISWI 型リモデラーとの相互作用面に位置することが確認され、電荷逆転変異 (E→K) はアシッドパッチの電気的特性を根本的に改変すると予測された (Fig 1)。
H2BE113K変異によるコロニー形成能の選択的亢進： MDA-MB-231 (三重陰性乳癌細胞株) に CRISPR-Cas9 を用いて H2BE113K 変異をノックインした 2 クローン (EK38, EK41) を 2 種類の WT クローン (WT1, WT38) と比較した。コロニー形成アッセイ (10 日間培養、crystal violet 染色) では E113K 変異細胞が WT と比較して有意に高いコロニー形成能を示した (P<0.0001、mean ± SD、独立実験 n=3 回) (Fig 2A-B)。一方、CCK8 細胞増殖アッセイ (OD450nm) では WT と E113K の間に有意差は認められず、scratch assay による遊走能も変化しなかった。すなわち H2BE113K はアンカー非依存性増殖 (コロニー形成) を選択的に亢進させるが、単純な細胞増殖速度や遊走能には影響しないことが示された (Fig 2C-D)。
ゲノムワイドエピゲノム解析によるクロマチン動態変容の同定： RNA-seq 解析 (P<0.01 閾値、STAR v2.7.7a + Cufflinks v2.2.1 解析) で H2BE113K 細胞の発現変動遺伝子 (DEG: differentially expressed gene) を同定したところ、上方制御遺伝子は細胞接着、クロマチンリモデリング、細胞周期制御に関連する経路に濃縮された。トランスポゼース利用クロマチンアクセシビリティ解析 (ATAC-seq: assay for transposase-accessible chromatin with high-throughput sequencing) では上方制御遺伝子の 45.8% がクロマチンアクセシビリティの増大を示した一方、増大部位の 69% は DEG ではない遺伝子に対応しており、H2BE113K がクロマチン開口を広汎に誘発することが示された (Fig 3)。
FLAG タグを利用したクロマチン切断と溶出法 (CUT&RUN: cleavage under targets and release using nuclease) (yeast spike-in 正規化、MACS2 P<0.001) では H2BE113K 変異体 H2B の結合ピークが WT との共局在パターン (G1 サブグループ: ピーク全体の 50.8%) と E113K 特異的局在パターン (G2 サブグループ: 49.2%) に K-means クラスタリングで分離された。G2 サブグループでは H2B 結合強度の増大量が G1 の 2 倍 (10% vs 5%) であり、全 ATAC-seq ピーク (64,151 ピーク) のうち FLAG ピークとオーバーラップするものは 9,440 個 (14.7%) であった。さらにクロマチンアクセシビリティが増大したピークの 93% が G2 ピーク領域と重複しており、G2 サブグループが E113K 変異誘発性クロマチン開口の主要な場であることが示された。RNA-seq 上方制御遺伝子・G2 FLAG ピーク・ATAC-seq 増大ピークの三重重複解析から 221 候補遺伝子 (269 ゲノム領域) が絞り込まれた (Fig 4)。
G3BP2の下流エフェクターとしての同定とSMARCA5依存的発現制御機序： 三重重複候補 221 遺伝子中、GTPase 活性化タンパク結合タンパク 2 (G3BP2: GTPase-activating protein-binding protein 2) が特に注目された。H2BE113K 細胞では G3BP2 の mRNA が 1.5 倍上昇 (P=0.0029) し、タンパク質は 2.3 倍に増加した (P<0.0001)。TCGA 乳癌コホート (810 腫瘍 vs 291 正常組織) では G3BP2 の高発現が確認され、Kaplan-Meier 生存解析では高発現群の予後不良が示された。G3BP2 shRNA ノックダウン (sh20, sh21) は E113K 変異細胞のコロニー形成能を WT 非処理細胞相当まで回復させ (Fig 5B)、G3BP2 が H2BE113K 誘発性コロニー形成増強の必須エフェクターであることが証明された。
ChIP-qPCR 解析 (MNase 消化、FLAG ビーズ、Protein G Sepharose) により、E113K 変異細胞では WT と比較して G3BP2 プロモーター領域 (-1 ヌクレオソーム位置) への SMARCA5 (ISWI ファミリー クロマチンリモデリング ATPase) の濃縮が有意に増大した (P<0.05 および P<0.01) (Fig 5C)。SMARCA5 は ISWI 型クロマチンリモデリング複合体の触媒サブユニットとして TSS 付近のヌクレオソームをスライドし転写因子結合部位を露出させることが知られており、H2BE113K によるアシッドパッチ変化が SMARCA5 動員を促進し G3BP2 発現を亢進させるというモデルが支持された。
SMARCA5阻害薬ED2-AD101の変異選択的抗腫瘍効果： First-in-class SMARCA5 阻害薬 ED2-AD101 (特許番号 US9926316B2) を H2BE113K 変異細胞および WT 細胞に処置したところ、50% 増殖抑制濃度 (IC50) は両細胞間で同等であった。しかしコロニー形成抑制効果は E113K 変異細胞において WT を有意に上回った (P<0.001) (Fig 5D-E)。RNA-seq で ED2-AD101 処置 H2BE113K 細胞を未処置 E113K 細胞と比較したところ、E113K で特異的に上昇していた遺伝子群が選択的に抑制され、Wnt シグナル経路・細胞恒常性関連経路の遺伝子が濃縮されており、腫瘍微小環境との相互作用経路も含まれた (Fig 5F)。
H2BE113Kノックインマウスにおける肺転移促進と転移関連遺伝子プログラム： C57BL/6 バックグラウンドに H2BE113K 変異をノックインしたマウス (H2BE遺伝子ターゲティング、Gemparmatech 社作製) を MMTV-PyMT 乳腺腫瘍マウスモデル (マウス乳腺腫瘍ウイルス-ポリオーマウイルス中間 T 抗原; Jackson Laboratory #022974) と交配した。H&E 染色と Ki67 免疫染色を行った WT/WT (n=19) および WT/E113K (n=17) の比較では、原発腫瘍の発生潜時、腫瘍重量、生存期間に有意差は認められなかった。一方、肺転移発生率は 2 バッチ (n=8 と n=11) ともに H2BE113K マウスで有意に増加し、肺転移巣の Ki67 染色陽性率も WT を上回った (P<0.05) (Fig 6A-B)。
肺転移巣 RNA-seq (GRCm38/mm10 アライメント、DESeq2 解析、Molecular Signatures Database (MSigDB) 遺伝子セット) では主成分分析 (PCA: principal component analysis) において WT と E113K の転移巣が明確に分離し、6 クラスターの DEG が同定された。クラスター 3 (145 遺伝子) は細胞外マトリックス (ECM: extracellular matrix) 組織化・白血球接着に濃縮され、クラスター 4 (149 遺伝子上昇) とクラスター 6 (62 遺伝子低下) は H2BE113K 転移巣に特異的であった。クラスター 6 の遺伝子は EGFR シグナリング正制御・神経突起形成・管腔形態形成に濃縮されていた (Gene Ontology (GO) 解析、P<0.05) (Fig 6C-G)。細胞モデルと in vivo モデルの両方で上方制御が確認された遺伝子として SEMA3C (セマフォリン 3C) と GPX3 (グルタチオンペルオキシダーゼ 3) が抽出され、RT-qPCR で E113K 転移巣における mRNA 増加が有意に確認された (P<0.05、独立実験 n=3) (Fig 6H)。
考察/結論
① 先行研究との違い: これまで解析されてきたオンコヒストン変異の多くは H3K27M (小児神経膠腫)、H3K36 変異 (軟骨芽細胞腫・頭頸部癌)、H2BG53D (膵臓癌) などであり、原発腫瘍の増殖制御や修飾酵素活性変化に焦点が当てられてきた。これと異なり、本研究が示した H2BE113K の最大の特徴は「原発腫瘍の成長には影響を与えず、転移能を選択的に亢進させる」という表現型分離にある。MMTV-PyMT モデルでは腫瘍発生潜時・腫瘍重量・生存期間が不変のまま肺転移が有意に増加しており、この転移選択性は H2BE113K を他のオンコヒストンと区別する際立った特徴である。また、H2B アシッドパッチ残基 (E113) が SMARCA5/ISWI 型リモデラーとの相互作用を制御する機能的ホットスポットであることも、以前の H2BG53D や H2BE76K の報告では明確にされていなかった新たな知見である。
② 新規性: 本研究は初めて CRISPR-Cas9 ノックイン変異体を用いて H2BE113K がヌクレオソームアシッドパッチを介し SMARCA5 の G3BP2 プロモーター動員を促進することを ChIP-qPCR で実証した。これは ISWI 型リモデリング複合体の基質特異性がアシッドパッチの電荷変化によって「書き換えられる」という、新規な oncohistone-remodeler 相互作用モデルを初めて提示するものである。さらに、CUT&RUN による G2 特異的ピークと ATAC-seq・RNA-seq の三重重複から 221 候補遺伝子を体系的に同定した多層エピゲノムアプローチも、オンコヒストン機能解析の新規なパラダイムを示す。また、chromatin-targeting アプローチとして BET タンパク質 BD1/BD2 選択的阻害が示されてきた文脈で (Gilan et al. Science 2020)、本研究は ISWI 型リモデラー SMARCA5 の変異選択的阻害という新規な治療コンセプトを初めて in cellulo で実証した。
③ 臨床応用: H2BE113K は乳癌の約 0.6% に存在するため、2022 年の世界年間診断数 230 万人を基準とすると潜在的対象患者は約 1.4 万人と推計される。本研究は MMTV-PyMT マウスモデルで肺転移選択的な促進効果を実証しており (Wculek et al. Nature 2015)、転移抑制に特化した介入が標準治療の抗腫瘍効果を補完できる可能性がある。ED2-AD101 は SMARCA5 に対する first-in-class 阻害薬として変異細胞でのコロニー形成を選択的に抑制しており、H2BE113K 変異スクリーニングと組み合わせた precision oncology アプローチへの展開が期待される。Wnt 経路など ED2-AD101 が抑制する経路は複数のがん種で治療標的として研究されており、H2BE113K 以外の適応への拡張可能性も考察される。
④ 今後の課題: SEMA3C と GPX3 が H2BE113K 肺転移において具体的にどのような機能的役割を果たすか (細胞接着促進か免疫回避かなど) は未解明であり、今後の課題として残されている。本解析は単一の TNBC 細胞株 (MDA-MB-231) と単一のマウス系統 (C57BL/6-PyMT) に基づいており、他の乳癌サブタイプや患者由来 xenograft モデルでの検証が必要である。また、H2BE113K 変異を持つ実臨床患者のコホートにおける転移パターンや生存転帰との関連、および ED2-AD101 の前臨床毒性・PK 評価が今後の研究として残されている。さらに SMARCA5 が G3BP2 以外のどのゲノム領域で転写プログラムを変容させるかの網羅的理解は、H2BE113K 依存性転移の全貌解明に向けた重要な方向性である。
方法
研究デザイン: CRISPR-Cas9 ノックイン細胞モデル (MDA-MB-231 / H2BC17 遺伝子) + in vivo マウスモデル (C57BL/6 × MMTV-PyMT) による基礎科学研究。
細胞株: MDA-MB-231 (DMEM + 10% FBS、37°C、5% CO2)。CRISPR-Cas9 ノックイン: pBlueScript 修復鋳型 + lentiCRISPRv2 (sgRNA/Cas9 共発現) を PEI (DNA:PEI=1:3) でトランスフェクション → neomycin (2 mg/ml) 選択 → Sanger sequencing で同定。
マウス: C57BL/6 H2BE113K ノックインマウス (Gemparmatech 社) × MMTV-PyMT (#022974、Jackson Laboratory)。WT/WT (n=19) vs WT/E113K (n=17)。動物倫理: City University of Hong Kong 承認 (no. A-0581)。
RNA-seq (細胞): NEBNext rRNA 除去 → ライブラリ調製 → HiSeq 150bp PE シーケンス。STAR v2.7.7a (hg38) アライメント、Cufflinks v2.2.1 発現定量、Cuffdiff P<0.01 差次発現解析。
CUT&RUN-seq: 細胞 500 万個、FLAG 抗体 (Sigma F1804) / H2A 抗体 (Abcam ab15653)、pA-MNase 消化、yeast spike-in 正規化、HiSeq X 10 (~2000 万 PE リード)。MACS2 v2.2.5 ピーク呼び出し (P<0.001)、BEDTools v2.29.2 カウント、Cluster 3.0 K-means クラスタリング (G1/G2 分類)。
ATAC-seq: 細胞 5 万個、digitonin 溶解、Tn5 トランスポゼース挿入、Bowtie2 v2.3.4.1 (hg38) アライメント、MACS2 v2.2.5 ピーク呼び出し (P<0.001)。
ChIP-qPCR: MNase 消化クロマチン、FLAG ビーズ (Sigma F7425) / H2A 抗体 (Abcam ab156533) / H3 抗体 (Immunoway YM3038)、Protein G Sepharose、chelex DNA 抽出 → qPCR。
in vivo RNA-seq: Trim Galore v0.6.10 → Hisat2 v2.2.1 (GRCm38/mm10) → HTSeq v0.11.3 → DESeq2 差次発現解析。分子シグネチャデータベース MSigDB (Molecular Signatures Database; Reactome / KEGG / canonical / hallmark) 遺伝子セットで分子シグネチャ解析。
統計: 独立 2 群は unpaired two-tailed Student’s t-test (mean ± SD、n=3 独立実験)。P<0.05 を有意とした。
データ共有: 次世代シーケンスデータは Gene Expression Omnibus (GEO) 登録番号 GSE295336 に登録。