- 著者: Chu N, Patel A, Wildes F, Navarro A, Conforte A, Zecha J, et al.
- Corresponding author: 未記載 (多施設共同 CPTAC)
- 雑誌: Nature Cancer
- 発行年: 2026
- Epub日: 2026-06-12
- Article種別: Resource
- DOI: 10.1038/s43018-026-01175-6
背景
急性骨髄性白血病 (AML) は遺伝的・表現型的・臨床的に高度に不均一な疾患であり、WHO 2022 分類や ELN (European LeukemiaNet; 欧州白血病ネット) 勧告に基づくリスク層別化が進んでいるものの、分子異常が多彩なため統一的な治療戦略の確立は困難なまま未解明の部分が多い。DiNardo et al. (2018) はイボシデニブが IDH1 変異 AML に対して持続的寛解を誘導することを示したが (DiNardo et al. NEnglJMed 2018)、サブタイプ横断的な分子基盤と、対応する治療的脆弱性の体系的な記述は不足していた。Rowles et al. (2026) はがんにおけるメタボロミクスの活用可能性を概説したが (Rowles et al. NatRevCancer 2026)、実際の AML コホートへの多層オミクス統合は限定的であった。さらに Mandal et al. (2023) は AML 特異的サーフェオミクス標的を同定したものの (Mandal et al. NatCancer 2023)、PTM (post-translational modification; タンパク質翻訳後修飾)・代謝・脂質を同時に統合した包括的アトラスは存在しなかった。Ruggero et al. (2026) は翻訳制御を介したタンパク質発現リモデリングが治療標的となりうることを示し (Ruggero et al. JClinInvest 2026)、AML においても翻訳因子 EIF4F (eukaryotic initiation factor 4F; 真核生物翻訳開始因子 4F) 複合体の役割が注目されていた。このように、AML の分子的多様性を「遺伝子型」と「細胞分化」の双方から定量的に記述し、対応する薬物標的を実証する統合フレームワークの手薄な状態が大きな課題であった。
目的
治療歴のない AML 患者 n=173 例を対象に、RNA・グローバルタンパク質・リン酸化・アセチル化・糖鎖修飾・代謝産物・脂質を含む 13 種のオミクスモダリティを統合し、タンパク質中心の AML サブタイプ AML-8 (acute myeloid leukemia 8-subtype classification) を定義するとともに、サブタイプ別の治療的脆弱性を外部コホートで検証する。
結果
AML-8 タンパク質中心サブタイプの確立と MYC/mTOR 二極対立による AML 生物学の再定義： SNF (similarity network fusion; 類似性ネットワーク融合) アルゴリズムを RNA と全タンパク質データに適用し、n=163 例から 8 つのタンパク質中心サブタイプ (AML-8) を同定した。C1 は NPM1 変異を持つ成熟単球性 AML、C2 は NPM1/FLT3-ITD (FLT3 fms-like tyrosine kinase, internal tandem duplication)/DNMT3A 三重変異の未分化 Mito-AML (既報 Mito-AML の精緻化)、C3 は FOXC1・HOXB8/9・SOX4・POU2F1 の外れ値発現を示す新規 NPM1 変異サブタイプ、C4 は CEBPA 変異 AML、C5-C7 は MDS 関連 sAML、C8 は RUNX1-RUNX1T1 転座 t(8;21) の GMP (granulocyte-monocyte progenitor; 顆粒球・単球前駆細胞) 様 AML と特徴づけられた。C3 では CIITA (class II immune transcription activating factor) プロモーターの高メチル化によって MHC (major histocompatibility complex; 主要組織適合複合体) クラス II が抑制されており、免疫回避機構が示唆される (Fig. 2)。連続軸として MYC 活性と mTOR 活性は鋭く対立し、Pearson r=0.65 の強い逆相関 (p<0.001) を示し、未分化 (MYC 高/mTOR 低) と分化 (mTOR 高/MYC 低) という AML の二極を定量化した。111,993 の分子フィーチャーを包含する本データセットは既存のいかなる単一オミクス分類よりも AML の分子的多様性を高解像度で描写する。
MAP1A が未分化 AML の普遍的タンパク質マーカーとして同定され、ベネトクラックス感受性と相関する： 多変量解析により、微小管会合タンパク質 MAP1A (microtubule-associated protein 1A) が NPM1 変異 AML から MDS 関連 AML にわたる全サブタイプの未分化 AML に共通するタンパク質マーカーとして浮上した (Fig. 3)。MAP1A は外部 BeatAML-proteomics コホートにおいてベネトクラックス薬剤応答 AUC と有意な逆相関を示し、未分化・ベネトクラックス感受性 AML の診断マーカーとして機能する。さらに MAP1A は正常造血幹前駆細胞と白血病幹細胞を区別するシグネチャーとも一致した。EIF4F 複合体 (EIF4A1/EIF4E/EIF4G1 から同定した 40-55 の翻訳標的) は GMP 様 C8 AML で顕著に低下し、同サブタイプのアミノ酸・ジペプチド蓄積やプロテアーゼ (ELANE, PRSS57) 高発現と呼応した (Fig. 3)。未分化 AML では逆に OXPHOS (oxidative phosphorylation; 酸化的リン酸化) が亢進し、ミトコンドリア・mtDNA の増加や ATP 産生系酵素 (HK1/2/3, FBP1/2) の発現増大が観察された。
CEBPA 変異 (C4) のミトコンドリア超アセチル化と GMP 様 C8 の全細胞低アセチル化という対称的な PTM 異常が判明する： 非メジアン補正アセチル化データの解析により、C4 サブタイプにおけるミトコンドリア限局性の超アセチル化と、C8 における全細胞コンパートメント横断的な低アセチル化という対照的なパターンを発見した (Fig. 5b, c)。C4 の超アセチル化は酵素非依存性のミトコンドリア反応であり、アセチル-CoA 蓄積によって駆動される。PTM-MOFA (post-translational modification multi-omics factor analysis; 翻訳後修飾多変量因子解析) の因子 F1 がアセチル化データに支配されかつ C4 に強く関連することを確認した。C4 では FA (fatty acid; 脂肪酸) β 酸化と lysine 分解経路の活性が有意に上昇し、TCA (tricarboxylic acid; トリカルボン酸) 回路遺伝子の同時上方制御が観察された。注目すべき例外は SUCLG2 の低下であり、TCA 回路フラックスの逆転とアセチル-CoA 蓄積が誘導されると考えられる (Fig. 5g, h)。一方 C8 では HAT (histone acetyltransferase; ヒストンアセチル化酵素) および HDAC (histone deacetylase; ヒストン脱アセチル化酵素) のタンパク量は変化しないが、EP300 K1794 の脱アセチル化 (EP300 HAT 活性に重要な部位) や CREBBP K1216 の変化を含む PTM レベルの規制異常が生じており、PTM 駆動型の制御機構を示唆した (Fig. 5j, k)。これらアセチル化異常は分化ベースの AML 分類では完全に説明されず、AML-8 固有の生物学を反映する。
FunMap ネットワーク解析が MTA1 (NuRD 複合体) をパノビノスタット耐性の主要ドライバーとして同定し CRISPR KO が感受性を回復させる： AML コホートと外部 6 データセットを統合した FunMap コファンクションネットワーク (13,729 遺伝子, 134,100 エッジ) を構築し、ICE (iterative clique enumeration; 反復クリーク列挙) アルゴリズムで 274 の完全連結クリークを抽出した (Fig. 6c)。このうち未分化 AML で高発現かつ不良な生存に関連する 7 クリークを同定し、なかでも C38 (HDAC1・HDAC2 を含む 11 タンパク質、NuRD (nucleosome remodeling and deacetylase; ヌクレオソームリモデリング脱アセチル化複合体)・Sin3 複合体関連) がパノビノスタット感受性と最も強い負の相関を示した。C38 の中で MTA1 (NuRD 複合体スキャフォールドタンパク) が最強の負相関を持ち (Fig. 6f)、かつ正常骨髄対照と比べて AML で有意に高発現であることを確認した。CRISPR-Cas9 による MTA1 ノックアウトは MOLM-14 および MONO-MAC-6 の 2 細胞株においてパノビノスタット感受性を回復させた。さらに MTA1 過剰発現は既知のパノビノスタット感受性 MOLM-14 quizartinib 耐性株に耐性を付与し (Fig. 6g, h, i)、MTA1 が HDAC 阻害剤耐性の必要十分因子であることを実証した。
代謝・脂質プロファイルの AML 表現型別分類と mTOR・ATP1B3・FGR を軸とする治療的脆弱性の体系化： メタボロミクス解析 (n=91 検体, 1,780 フィーチャー, 433 同定) では ICA (independent component analysis; 独立成分分析) により 8 因子を抽出し、F1 (mTOR/分化型) と F4 (MYC/未分化型) が AML の二極を定量した (Fig. 4)。脂質解析 (n=96 検体, 1,365 脂質) では、フォスファチジルセリン (PS) とスフィンゴミエリン (SM) が分化型で高く、フォスファチジルコリン (PC) が未分化型で高く、ヘキソシルセラミド (HexCer) および Hex2Cer が GMP 様 C8 で特異的に蓄積した。分化型 AML では PUFA (polyunsaturated fatty acid; 多価不飽和脂肪酸) 22:6 含有エーテルリン脂質が増加しフェロトーシスシグネチャーが豊富で、AIFM2 および GPX4 (フェロトーシス抑制タンパク) が発現上昇していた (Fig. 4j, k)。これらは分化型 AML でのフェロトーシス誘導療法の可能性を示唆する。治療的脆弱性については、XGBoost (gradient boosting classifier; 勾配ブースティング分類器) で BeatAML-proteomics コホートのサブタイプを予測し FDA 承認薬への応答を解析した結果、mTOR はラパマイシン応答が C1・C8 で確認された。ATP1B3 (膜タンパク) が C1・C8 のトップ 10 標的に挙がり (Fig. 7d)、合成致死パートナー ATP1B1 の低発現が依存性を裏付けた (Fig. 7f)。FGR (Feline sarcoma proto-oncogene kinase; Src ファミリーキナーゼ) は単球性分化スコアとの最強正相関を示し (Fig. 7h)、BeatAML コホートにおいてダサチニブ+ベネトクラックス併用が分化型 AML で改善効果を示す n=281 例は非改善 n=106 例より FGR タンパク発現が有意に高く、FGR が治療応答バイオマーカーとなりうることが示された (Fig. 7j)。
考察/結論
① 先行研究との違い: これまでの AML 分類は WHO/ELN の遺伝子型ベース (Severens et al.) または Jayavelu et al.・Pino et al. による細胞分化ベースの単一軸分類に留まり、両者を同時に統合した枠組みは存在しなかった。本研究が提示する AML-8 は遺伝子型 (NPM1, CEBPA, RUNX1-RUNX1T1) と分化状態 (未分化/分化型/GMP 様) を統合した点でこれまでと異なる次元のサブタイプ体系であり、既報の Mito-AML が NPM1/FLT3-ITD/DNMT3A 三重変異の未分化サブタイプ (C2) に対応することを精緻化した。
② 新規性: MAP1A が NPM1 変異 AML から MDS 関連 AML まで未分化表現型の普遍的タンパク質マーカーとして機能し、ベネトクラックス感受性と連動するという知見は新規に得られた発見である。また MTA1 が NuRD 複合体スキャフォールドとしてパノビノスタット耐性を主導するという機構も、FunMap コファンクションネットワークという新規なアプローチによって初めて同定された。CEBPA 変異 C4 サブタイプにおけるミトコンドリア非酵素的超アセチル化が SUCLG2 欠失 + FA β 酸化亢進 → アセチル-CoA 蓄積という代謝経路に起因することも、PTM と代謝の統合解析がなければ見えなかった新規な知見である。
③ 臨床応用: MAP1A の発現が未分化 AML とベネトクラックス感受性を同定するバイオマーカーとして機能するならば、臨床応用として治療前の MAP1A 測定がベネトクラックス+アザシチジン療法の患者選択に役立つ可能性がある。FGR タンパク発現はダサチニブ+ベネトクラックス併用の恩恵を受ける分化型 AML 患者を前向きに選別するバイオマーカーとして臨床的意義をもつ。mTOR 阻害薬ラパマイシンの C1・C8 サブタイプへの感受性は独立した ex vivo モデルで検証されており、サブタイプ別 mTOR 阻害戦略の臨床での展開を支持する。MTA1 を標的とした NuRD 複合体阻害が HDAC 阻害剤耐性を克服しうるという知見も、将来の臨床応用に向けた開発基盤となる。
④ 残された課題: 本コホートは欧州系を主体とし健常対照が限られるため、全地域・全 AML サブタイプを網羅するには至っておらず、残された課題として国際的な多人種コホートによる検証が不可欠である。機械学習によるサブタイプ転移の精度向上、AML-8 と小児 AML・正常造血との比較、そして MAP1A・FGR・MTA1 の機構的実証と前向き臨床試験での予測妥当性の評価が今後の課題として残る。技術的限界として一括バルクプロファイリングは細胞集団変化と真の制御変化を分離できず、2-HG (イソクエン酸) 等の一部代謝物は MS で定量できなかった点も課題である。
方法
研究デザイン: 多施設共同前向き観察研究。CPTAC (Clinical Proteomic Tumor Analysis Consortium; 臨床プロテオミクス腫瘍解析コンソーシアム) プロトコルに準拠し、全例で書面同意取得・倫理委員会承認済み。過去 12 か月以内の他悪性腫瘍歴または前治療歴を有する患者は除外。
対象: 治療歴のない AML 患者 n=173 例。主要変異の内訳: Fusion 19% (n=33), NPM1 31% (n=53), DNMT3A 33% (n=57), FLT3 36% (n=63), CEBPA 12% (n=20), TP53 12% (n=21), TET2 17% (n=30), sAML 26% (n=45)。腫瘍純度中央値 79% (IQR 73-92%)。サブタイプ分類は WHO 2022 基準に基づく階層的アルゴリズムを使用。
オミクスプラットフォーム: RNA-seq, WGS, global proteome (TMT 標識), リン酸化, アセチル化, 糖鎖修飾 (site-level), ATAC-seq, metabolomics (RP + HILIC, n=91), lipidomics (正・負イオンモード, n=96/97) の計 13 モダリティを統合。PTM 正規化は TMT-Integrator ratio と親タンパク質発現の線形回帰残差を使用 (非メジアン補正アセチル化を別途解析)。RBC 汚染は 14 マーカー遺伝子スコアで推定し limma removeBatchEffect で除去。
サブタイプ解析: SNF (SNFtool R パッケージ) を RNA + タンパク質データで適用、Euclidean 距離・α=0.5・20 近傍・10 回拡散反復を使用。クラスター数 3-10 および特徴量 10-100% の安定性を評価。転写因子活性は SCENIC および decoupleR (DoRothEA/CollecTRI regulon) で推定。代謝解析は MOFA2 による三種 MOFA (multi-omics factor analysis; RNA/protein, metabolomics/lipidomics, PTM)。
FunMap & 治療的脆弱性: AML 特異的コファンクションネットワークを 6 外部 AML データセットと統合し LLR (log-likelihood ratio; 対数尤度比) > 3.912 で構築 (13,729 遺伝子, 134,100 エッジ)。ICE アルゴリズム (最小サイズ 5) で 274 全連結クリーク抽出。治療標的は RNA ベース XGBoost 分類器 (500 遺伝子) で DepMap 細胞株を AML-8 へマッピング後、タンパク過剰発現 + CRISPR 依存スコアで 5 ティア評価。
検証コホート: BeatAML-proteomics コホート (外部 AML データセット) で薬剤応答 AUC を比較 (単剤ベネトクラックス/ダサチニブ, 最良応答; 分化型 n=32 vs 未分化型 n=31)。FGR 改善群 n=281 vs 非改善群 n=106。CRISPR-Cas9 MTA1 KO 実験: MOLM-14 および MONO-MAC-6 細胞株でパノビノスタット感受性を評価 (Supplementary Table 33)。
統計: 差次発現解析に limma + eBayes 関数 (log2 scale)。RNA-seq は voom 正規化。PTM-SEA でキナーゼ活性推定。ssGSEA および fGSEA で経路エンリッチメント (FDR < 0.05)。生存解析は Kaplan-Meier 法。ピアソン相関係数使用。