- 著者: Daniel Price, Grigory Zemlyanskiy, Watanya Trakarnphornsombat, Ignasi Forne, Nadezda Volkova, Louis Dubusse, Alexandra J. Cooper, Axel Imhof, Aliaksandra Radzisheuskaya, et al.
- Corresponding author: Aliaksandra Radzisheuskaya (Division of Cell and Molecular Biology, The Institute of Cancer Research, London, UK)
- 雑誌: Nature Genetics
- 発行年: 2026
- Epub日: N/A
- Article種別: Original Article
- DOI: 10.1038/s41588-026-02675-y
背景
ヒストンのPTM (post-translational modification; 翻訳後修飾) はクロマチン構造の調節とゲノム機能の制御に中心的な役割を果たす。しかし、哺乳類における個別のPTMの機能的寄与を解明することは、ゲノム全体に分散した多数のヒストン遺伝子コピーの存在により著しく困難であった。これまでの機能研究の多くは、特定のPTMを付加する酵素の阻害や遺伝子欠損に依存してきたが、こうした手法は酵素の基質多様性、触媒活性に依存しない機能、そして機能的冗長性によって解釈が複雑化するという根本的な限界を抱えていた (Ciabrelli et al. Sci. Adv. 2023; Hunt et al. Mol. Cell 2022)。
酵母 (Saccharomyces cerevisiae) では2コピーのヒストン遺伝子という単純な構成を活かした系統的なヒストン変異ライブラリー解析が実施されてきたが (Nakanishi et al. Nat. Struct. Mol. Biol. 2008; Dai et al. Cell 2008)、哺乳類では正準型ヒストンH3.1/H3.2を9+3の合計12遺伝子が担い、非正準型バリアントH3.3が別途H3f3a (histone H3.3 gene A) とH3f3b (histone H3.3 gene B) の2遺伝子から発現するという複雑な構成が障壁となっていた。ショウジョウバエでは単一ゲノム遺伝子座に集中するヒストン遺伝子の利点を活かした網羅的変異解析が行われ (McKay et al. Dev. Cell 2015; Zhang et al. DevCell 2019)、H3K4・H3K9・H3K14・H3K79など複数のリジン残基の発育致死性が確認されているが、哺乳類での同等の機能マップは未解明のまま残されていた。哺乳類において全ヒストンH3リジン残基の適応度への寄与を定量評価できる系統的高スループット手法が不足していたことがこの空白の根本原因である。
哺乳類への応用として、CRISPR塩基編集を用いたH3K27R変異マウス胚性幹細胞 (embryonic stem cell; ESC) 作製が試みられたが (Sankar et al. Nat. Genet. 2022、Sankar et al. NatGenet 2022)、塩基編集は編集窓の狭さ・バイスタンダー変異・アミノ酸置換の選択肢の制限という問題点を有していた。CRISPR prime editingはCas9ニッカーゼと逆転写酵素の融合タンパク質、およびprime editing guide RNA (pegRNA) を用いて二本鎖DNA切断や鋳型DNAを必要とせず正確な塩基置換・小挿入・小欠失を導入できる次世代技術であるが (Anzalone et al. Nature 2019; Nelson et al. NatBiotechnol 2022)、哺乳類における多コピーヒストン遺伝子への系統的適用はこれまで実証されていなかった。本研究はprime editingを活用した高スループットプラットフォームの開発とH3全リジン残基の機能マッピングを目的とした。
目的
マウスESCにおいてCRISPR prime editingを用いて正準型・非正準型ヒストンH3の全リジン残基を系統的にK→R (lysine-to-arginine) 変異させるプラットフォームを確立し、各リジン残基のESC適応度への寄与を定量的に評価することで、哺乳類ヒストンH3リジン機能マップを構築する。
結果
PE5bシステムの最適化とH3K23R完全変異株の確立：PE5b (prime editing version 5b) を最適化するため、PEmax (prime editing maximized) リバーストランスクリプターゼ融合型Cas9ニッカーゼ、MLH1dn (dominant-negative MLH1 mutant)、およびBFP (blue fluorescent protein) レポーターを含むTPMB (TetO promoter, PEmax, MLH1dn, BFP reporter) コンストラクトのドキシサイクリン誘導型ESC株を樹立した。H3K23を標的とした検証では、PE4およびPE5bが最高の編集効率と最低のインデル頻度を示した (n=3 biological replicates; Fig. 1c)。K23R変異クローンの35%以上が80%以上の編集効率を達成し、そのうち2クローンが18全H3.1/H3.2アレルで完全変異を確認された (Fig. 1i)。Western blotによりH3K23acシグナルはほぼ完全に消失 (q<0.0001, n=3 independent experiments, Fig. 1j) しており、prime editingが正準型ヒストンH3全コピーへの正確な変異導入に成功したことが示された。H3R23K逆変換epegRNAを導入することでH3K23ac発現が完全に回復し、可逆的ヒストン変異導入も実証された (Fig. 3d)。H3K14K/H3K18K/H3K23Rの3残基同時変異epegRNAを用いた組合せ編集では、K→R三重変異は単一変異と比較して約4-fold to 5-fold lower編集効率にとどまり、多重変異の困難さが明らかとなった (Fig. 3g)。
高スループットスクリーンによる適応度必須リジン残基の同定：13のリジン残基すべてについてK→RおよびシノニマスK→K対照epegRNAを設計し、3生物学的反復 (n=3 independent biological replicates) でTMPB ESCをトランスフェクションして11日間培養した後、各リジン47クローンをqPCRでプレスクリーニングし、最上位12クローンをアンプリコンNGSで解析した (Fig. 6b)。K→Rの編集効率が対照K→Kと比較して著明に低下し、かつPAM (protospacer adjacent motif) only変異 (編集不完全産物) の蓄積が認められた残基を「適応度への負の効果あり」と判定した。その結果、H3K4R (最大編集率11〜12%、2/36クローンのみ; control K→Kと比較して約8-fold lower) ・H3K9R (最大39%) ・H3K14R (最大67%) ・H3K18R (最大50%) ・H3K79R (最大46%) の5残基で強いfitness defectが確認された (Fig. 6c–f)。一方、H3K23R・H3K36R・H3K37R・H3K56R・H3K64R・H3K115R・H3K122Rでは高編集クローンが多数回収され、ESC増殖において許容性を示した。H3K27Rは軽度の編集効率低下のみで、マイルドな増殖障害と合致する表現型を示した (Fig. 6c,d)。H3K64Rはむしろ対照K→Kよりも編集効率が高く、わずかな適応度ゲインを示した。
H3.3バリアントによるH3K18ac代償機構の解明：部分編集H3K18R変異株においてH3K18acシグナルがほぼ正常レベルに維持されていることから、H3.3バリアントによる代償を仮説した。H3.3K18R epegRNA/nickRNAをWT ESCおよび部分H3K18R ESCに導入したところ、H3f3a (H3.3をコードする遺伝子1)・H3f3b (H3.3をコードする遺伝子2) 両遺伝子で完全変異が達成された (Fig. 4e,f)。部分H3K18R + H3.3K18R二重変異ではH3K18acシグナルが有意に低下した (q=0.0003、n=3、Fig. 4g)。CUT&Tag (C&T) プロファイリングにより、H3K18Rにおいては増加したH3K18acがH3.3高局在領域 (エンハンサー・プロモーター) と強く相関し、H3.3K18R追加によりこの再分布が解消されることが示された (Fig. 4h–k)。さらに、液体クロマトグラフィー-タンデム質量分析 (LC-MS/MS) によりH3.1/H3.2K18R細胞ではH3.3/H3.1・H3.2比が有意に増加 (q=0.0104、n=3、Fig. 4l) しており、H3.3がクロマチン中への取込みを増加させることでH3K18ac維持に寄与するという補償機構が独立に支持された。
H3K56のゲノム安定性保全機能の哺乳類での保存：酵母・ショウジョウバエにおいてH3K56アセチル化 (H3K56ac) がDNA複製・修復後のクロマチン回復とゲノム安定性に関与することが知られていた。本研究では完全変異H3K56R ESCにトポイソメラーゼII阻害剤エトポシドを投与したところ、H3K56K対照と比較して約2-fold higher感受性が示された (p<0.05、n=3、Fig. 7j)。この結果はH3K56のゲノム安定性制御機能が酵母・ショウジョウバエから哺乳類まで進化的に保存されていることを初めて直接的に実証するものである。DNA損傷応答においてH3K56acが新規クロマチン取込みの合図として機能することは原核生物から真核生物まで広く保存された原理と一致する。
組合せ変異が明らかにするリジン間のfunctional crosstalk：単一変異体のうち適応度に負の影響を与えないH3K23R・H3K27R・H3K36R・H3K37Rの組合せ二重変異を逐次的に導入した。H3K27R+H3K36RおよびH3K27R+H3K37R二重変異では、それぞれの単一変異体と比較してコロニー形成能が著明に低下した (q<0.0001、q=0.0001、n=3、Fig. 8e)。RNA-seq解析によりH3K27R+H3K36R二重変異では5,967遺伝子に有意な遺伝的相互作用効果が検出され (q<0.01)、k平均クラスタリングにより epistatic・buffering・emergent応答が識別された (Fig. 8f,g)。Emergentグループでは多能性維持因子 Klf4・Chd4・Sin3a の下方制御と分化ドライバー Fgfr2・Dkk の上方制御が認められ、自己複製能低下と整合した転写変化パターンが示された。H3K36me2/me3とH3K27me2/me3がPRC2-EED相互作用を介して競合することが知られており (Strahl and Briggs Mol. Cell 2024)、本結果はこの拮抗関係の機能的帰結を哺乳類ESCで直接示したものである。
考察/結論
① 先行研究との違い： 従来、哺乳類における個別ヒストンリジンの機能解析はヒストン修飾酵素の遺伝子欠損・阻害剤を通じた間接的アプローチに依存しており、酵素の多基質性・非酵素機能・代償機能による解釈の困難が指摘されてきた。ショウジョウバエではヒストン遺伝子が単一ゲノム遺伝子座に集中するため系統的変異解析が可能であったが、哺乳類の遺伝子構成と異なり直接的な知見の外挿には限界があった。これまでの哺乳類ヒストン変異研究はCRISPR塩基編集によるH3K27R変異体作製（Sankar et al. 2022）に限定されており、全リジン残基を網羅する高スループット解析は存在しなかった。本研究のPE5b + TPMB ESCシステムは、対照K→K epegRNAとの比較によってpriming biasを除いた純粋な生物学的効果の定量を可能にし、従来の遺伝子KOアプローチとは異なるリジン特異的解析を実現した点で根本的に異なる。
② 新規性： 本研究は哺乳類ヒストンH3の全13リジン残基を標的とするCRISPR prime editing高スループットプラットフォームを新規に開発し、哺乳類系での初のリジン別適応度機能マップを確立した。特に、H3K18acの維持においてH3.3バリアントが代償的機能を担い、H3.3の当該ゲノム領域への優先的取込みがこの代償の分子機構であることを新規に示した。さらに、H3K56のゲノム安定性機能が哺乳類で保存されること、H3K27とH3K36の組合せ変異がemergent転写応答を誘導して自己複製能を障害することは、これまで報告されていない知見である。
③ 臨床応用・トランスレーショナル意義： 本研究で確立されたプラットフォームはがんのオンコヒストン (H3K27M・H3G34R/V等のヒストン変異型ドライバー) の機能解析に直接応用可能であり、小児神経膠腫や骨肉腫などヒストン変異を持つ腫瘍の病態理解に貢献できる臨床応用上の意義がある。H3K18Q (アセチル化模倣) 変異がオンコジェニックとして報告されており、CBP/p300阻害剤A485 (5〜10 μM) によるH3K18ac低下がESC増殖を障害した知見は、H3K18acを標的としたエピジェネティクス療法の有望性を示唆する。また、組合せ変異による転写制御の emergent応答の解析は、複数のエピジェネティック異常を持つがん細胞の新規治療標的同定に向けた方法論的基盤を提供する。
④ 残された課題： K→R変異は特定リジン残基の全PTMを同時に除去するため、個別の修飾 (メチル化・アセチル化等) への寄与を特定することが困難である。特に複数の修飾が共存するH3K27 (メチル化・アセチル化) のような残基では修飾特異的な結論を導くためにK→A/K→Q等の他の置換変異との組合せ検討が今後の課題として残されている。さらにH3K64・H3K115・H3K122のK→R変異はESC増殖への障害を示さなかったが、ショウジョウバエでは発生過程での致死性・生殖障害が報告されており、分化条件やストレス条件下での哺乳類表現型解析は未解明な部分が多い。今後の研究としてH3.3K9R変異体の解析によるバリアント間の機能的冗長性・特異性の解明、pan-H3 K9R三重変異体の表現型確認なども重要な方向性である。
方法
本研究はmouse embryonic stem cell line E14Tg2a (マウスESC、Jens Brickman博士より提供) を用いた基礎実験研究である。piggyBacトランスポゾンベクターを用いてPEmax-2A-MLH1dn-2A-BFP (TPMB) を安定導入したドキシサイクリン誘導型ESC株を構築した。各リジン残基に対してPrimeDesign (v.0.2) またはDeepPrime (v.3.6) でepegRNA・nickRNAを設計し、Lipofectamine 3000を用いて一過性トランスフェクション後、ドキシサイクリン (1 μg/ml) 添加下でPE5b prime editingを実施した。編集効率の評価にはOpentrons Flexロボットを用いたアンプリコンNGSを採用し (CRISPResso2 v.2.3.4で解析)、18全H3.1アレルを単一プライマーペアで網羅的にジェノタイピングした。47クローンをqPCRでプレスクリーニング後、上位12クローンをNGS定量した。ヒストンPTM変化はWestern blot (n=3 independent experiments; 一方向または二方向ANOVA + Dunnett/Tukey多重比較補正) およびSNAP-ChIP (Sequential Native Antibody-Panel ChIP, KAcyl-Stat panel) で評価した。ゲノムワイドなPTM分布はCUT&Tag (C&T, E. coli スパイクイン正規化) およびCUT&RUN (C&R, Drosophila スパイクイン正規化) で解析し、大規模ゲノム変化は浅層全ゲノムシーケンシング (shallow WGS、mm10 C57BL/6J reference genome、Bowtie 2 v.2.5.4、ACEパッケージv1.30.0) で検出した。H3.3/H3.1・H3.2比は質量分析 (LC-MS/MS、Orbitrap Exploris 480) でS31含有ペプチドの相対定量により評価した。RNA-seqはpoly(A)ライブラリー (n=3) を用いてNEBNext Ultra II kitで作製し、HiSat2でmm10にアラインメント後、DESeq2 (v.1.52.0) で因子設計モデルによる遺伝子相互作用解析を実施した (Benjamini-Hochberg補正、FDR <0.01)。コロニー形成アッセイは200 cells/wellを12ウェルプレートに播種し6日後にアルカリホスファターゼ染色 (Leukocyte Alkaline Phosphatase Kit) で定量した。統計：一方向/二方向ANOVA + 多重比較補正 (Tukey/Dunnett/Sidak)、Fisher exact test (Benjamini-Hochberg補正)、両側t検定 (Benjamini-Krieger-Yekutieli補正)。データはSRA (PRJNA1268213)、GEO (GSE319420/GSE319506/GSE319059) に登録。