- 著者: Siyao Guo, Yutong Zou, Canfeng Zhang, Hui Han, Yishan Li, Yucong Sun, Zhaoyu Wang, Wanrui Zhang, Jiang Zhang, Yuli Gan, Qiang Zhang, Jieyi Ma, Siyi Zheng, Hongshen Qiu, Yan Zhu, Wange Lu, Yong Bao, Meng Zhao, Shuibin Lin
- Corresponding author: Yong Bao, Meng Zhao, Shuibin Lin (Sun Yat-sen University)
- 雑誌: Nature Cancer
- 発行年: 2026
- Epub日: 2026-05-21
- Article種別: Original Article
- PMID: 42457944
背景
mRNAのポリアデニル化（polyadenylation）は、転写物の安定性と翻訳効率を決定する極めて重要なmRNA成熟ステップである。生理的条件下において、ポリA鎖の長さは動的に制御されており、組織特異的な機能やストレス応答に関与することが報告されている (Passmore and Coller 2022, Subtelny et al. 2014)。がん細胞においても、代替ポリアデニル化（alternative polyadenylation; APA）などの異常なRNAプロセシングが、遺伝子発現の不均衡を招き、腫瘍形成や分化停止を促進することが既報となっている (Davis et al. 2022, Lee et al. 2018)。
特に急性骨髄性白血病（acute myeloid leukemia; AML）においては、治療抵抗性を持つ白血病幹細胞（leukemia stem cells; LSC）の維持メカニズムの解明が、再発防止と根治治療の実現に向けた最重要課題である。近年の研究により、エピトランスクリプトーム修飾やRNAプロセシングの異常がLSCの自己複製能に寄与することが示唆されているが、ポリA鎖の「長さ」そのものがAMLの病態進行や代謝リプログラミングにどのように関与しているかは未解明である。これまで、ポリA鎖の長さの動的制御ががんの進展に果たす役割についての知見は不足しており、その分子メカニズムを特定するための研究は未開拓の領域であった。したがって、AMLにおける異常なポリA鎖延長（hyperadenylation）の有無とその機能的意義を明らかにすることは、新たな治療標的を同定する上で極めて重要な gap であると考えられた。
目的
本研究の目的は、AML患者および細胞モデルにおいてmRNAのポリアデニル化状態を網羅的に解析し、異常なポリA鎖延長を駆動する酵素を同定することである。具体的には、フルレングスmRNAシーケンシング（full-length mRNA sequencing; FLAM-seq）を用いて、AMLにおいて特異的にポリA鎖が延長している標的遺伝子群を同定し、その制御因子を探索した。さらに、同定した酵素がAMLのLSCの自己複製能や病態進行にどのように寄与しているかを検証し、その下流で制御される代謝リプログラミングの分子メカニズム、特にTCA（tricarboxylic acid）サイクルへの影響を明らかにすることを目的とした。最終的に、この経路を標的とした薬理学的阻害剤であるcordycepinが、正常造血への影響を抑えつつAMLの進行を抑制できるかを検証し、新たな治療戦略としての有効性を提示することを目指した。
結果
PAPOLAによるmRNAの過剰ポリアデニル化と予後不良の相関: 研究グループは、AML患者由来のCD34+芽球と正常な造血幹細胞・前駆細胞（hematopoietic stem and progenitor cells; HSPC）を用いてFLAM-seq解析を実施した。その結果、AMLサンプルでは正常対照群と比較して、細胞代謝、細胞周期進行、ストレス応答に関連する一連のmRNAにおいて顕著なポリA鎖の延長（hyperadenylation）が認められた (Fig. 1b,c)。この現象を駆動する酵素を探索したところ、ポリAポリメラーゼalpha（poly(A) polymerase alpha; PAPOLA）がAMLにおいて特異的に高発現しており、TCGA (The Cancer Genome Atlas) データセットを用いたパンキャンサー解析でも、PAPOLAの発現上昇は他の癌種よりもAMLで最も顕著であった (Fig. 1d,e)。臨床検体を用いた検証でもPAPOLAのmRNAおよびタンパク質レベルの上昇が確認され (n=3 participants per group)、TARGET (The Oncology Genome and Transcriptome) データベースを用いたカプラン・マイヤー生存分析では、PAPOLA高発現群において全生存期間（OS）が有意に短縮することが示された (Fig. 1h)。また、コックス回帰分析により、PAPOLA高発現はAMLの予後不良を予測する独立した因子であることが同定された。
PAPOLA欠損によるAML進行の抑制とLSC自己複製能の低下: PAPOLAの機能的意義を検証するため、shRNAを用いたノックダウン（KD）およびマウスでの遺伝子欠損モデルを用いた。ヒト原発AML細胞においてPAPOLAをKDしたところ、細胞増殖が有意に抑制され、コロニー形成能が低下し、アポトーシスが増加した (Fig. 1j-l)。NOD/SCIDマウスを用いた異種移植モデルでは、PAPOLA欠損AML細胞を注入した群で腫瘍負荷が大幅に減少し、肝臓および脾臓の重量が低下し、生存期間が有意に延長した (Fig. 1m-u, n=6 mice per group)。また、MLL-AF9 (Mixed Lineage Leukemia-AF9) 駆動の白血病モデルにおいてPapolaをノックアウト（KO）したところ、末梢血および骨髄中のGFP+ AML細胞の割合が著しく低下し (Fig. 2d,e)、肝脾腫が抑制された (Fig. 2g,h)。さらに、二次移植モデルを用いた限界希釈分析（limiting dilution assay; LDA）の結果、Papola KO群におけるLSC頻度は 1/8,953 であり、対照群の 1/621 と比較して著しく低下していた (Fig. 3e)。一方で、Mx1-Creを用いた条件付きKOマウスでは、正常な造血系細胞やHSPCのサブポピュレーションに有意な影響は見られず、PAPOLAが正常造血には不要であることが示された (Fig. 3n-y)。
PAPOLA-GSTM2軸による代謝リプログラミングの駆動: PAPOLAがどのようにAMLを促進するかを解明するため、RNA-seq、RIP-seq (RNA immunoprecipitation sequencing)、FLAM-seqを統合的に解析した。その結果、PAPOLAが直接結合し、ポリA鎖の長さを制御して安定化させる標的として、酸化ストレス応答因子であるGSTM2（glutathione S-transferase mu 2）を同定した (Fig. 5a)。PAPOLAをKDするとGSTM2 mRNAのポリA鎖が短縮し、mRNAの半減期（t 1/2）が短縮してタンパク質発現が低下した (Fig. 5f-j)。GSTM2のKDは、PAPOLA KDと同様にミトコンドリアの酸素消費速度（oxygen consumption rate; OCR）を低下させ、ATP産生量を減少させ、一方で乳酸産生量および細胞外酸性化速度（extracellular acidification rate; ECAR）を上昇させた (Fig. 5p-t)。in vivoでのGstm2 KOモデルにおいても、AMLの進行が抑制され、生存期間が延長し、同様の代謝リプログラミング（OCR低下、ECAR上昇）が認められた (Fig. 6d-m, n=6 mice per group)。また、PAPOLA欠損AML細胞にGSTM2を過剰発現させたところ、AMLの進行と代謝状態がレスキューされたが、この効果はポリA鎖が長いGSTM2 mRNAを導入した場合にのみ認められた (Fig. 5u-y)。
GSTM2-HNE-DLD軸を介したTCAサイクルの制御: GSTM2が代謝を制御する詳細なメカニズムを解析したところ、GSTM2の欠損により、脂質過酸化産物である4-hydroxynonenal (HNE) が蓄積することが判明した (Fig. 7c,d)。蓄積したHNEはミトコンドリアタンパク質を分解し、特にTCAサイクルの鍵酵素であるDLD（dihydrolipoamide dehydrogenase）のタンパク質レベルを著しく低下させた (Fig. 7c)。LC-MS (liquid chromatography-mass spectrometry) を用いたメタボローム解析では、PAPOLAまたはGSTM2の欠損により、NAD+/NADH比が上昇し、TCAサイクル中間体であるsuccinyl-CoA、succinate、fumarateが減少し、上流のα-ketoglutarate (α-KG) が蓄積することが示された (Fig. 7a,b)。この代謝ブロックを回避するため、DLDの下流代謝物であるmethylsuccinateを投与したところ、PAPOLA欠損AML細胞の増殖能とOCR、ATP産生が部分的に回復した (Fig. 7e-m)。これにより、PAPOLA-GSTM2-HNE-DLD軸がTCAサイクルの完全性を維持し、ミトコンドリア酸化リン酸化を駆動することでAMLの進行を促進していることが証明された。
CordycepinによるPAPOLA阻害の治療効果: PAPOLAを標的とした治療戦略として、天然のポリAポリメラーゼ阻害剤であるcordycepinの有効性を検証した。Cordycepin処理により、GSTM2 mRNAのポリA鎖長が短縮し、その発現レベルが低下することが確認された (Extended Data Fig. 8b-d)。in vitroでは、dose-dependentにAML細胞の増殖が抑制され、アポトーシスが誘導された (Extended Data Fig. 8e-j)。MLL-AF9駆動AMLマウスモデルにおいてcordycepinを投与したところ、骨髄中のGFP+ AML細胞の割合が低下し (Fig. 8b-d)、肝脾重量が減少し、生存期間が有意に延長した (Fig. 8f-h, n=6 mice per group)。また、cordycepinはOCRを抑制し、乳酸産生を増加させることで代謝リプログラミングを誘導した (Fig. 8i-o)。重要な点として、cordycepin投与後の野生型マウスでは、末梢血の白血球分画や骨髄中のHSPCサブポピュレーションに有意な変化は見られず、正常造血への影響は極めて限定的であった (Fig. 8p,q)。
考察/結論
先行研究との違い: 従来のRNAプロセシング研究では、主に代替ポリアデニル化（APA）による3’ UTRの短縮などががんの進展に関与することが報告されてきた。しかし、本研究はそれらと異なり、ポリA鎖の「長さ」そのものが動的に制御され、それがAMLの代謝状態を決定する「代謝レオスタット」として機能していることを初めて示した。単なる配列の選択ではなく、ポリA鎖の延長という定量的変化が、mRNAの安定性と翻訳効率を介して細胞の代謝表現型を劇的に変化させるという視点は、これまでのRNA制御研究とは対照的なアプローチである。
新規性: 本研究で初めて、PAPOLAによるmRNAの過剰ポリアデニル化がGSTM2の安定性を高め、HNEの除去を通じてTCAサイクル酵素DLDを保護するという、RNAプロセシングとミトコンドリア代謝を直接的に結びつける新規なシグナル軸（PAPOLA-GSTM2-HNE-DLD軸）を同定した。特に、ポリA鎖の長さがGSTM2の機能維持に必須であり、それがLSCの自己複製能を担保していることを証明した点は極めて独創的である。この知見は、RNAの末端修飾が単なる保護機能を超えて、代謝経路の完全性を維持するスイッチとして機能することを明示したものである。
臨床応用: 本知見は、AMLにおけるPAPOLAの過剰発現が予後不良因子であることから、PAPOLAを標的とした治療の臨床的意義を強く示唆している。特に、天然化合物であるcordycepinが、正常造血への影響を最小限に抑えつつ、LSCの代謝脆弱性を突いてAMLを抑制できることを示したことは、bench-to-bedsideのtranslationalな展開において極めて有望である。LSC特異的な代謝依存性を標的とすることで、正常造血幹細胞への毒性を回避した高選択的な治療法の開発が可能になると期待される。
残された課題: 今後の検討課題として、PAPOLAがGSTM2以外にどのような標的mRNAを制御し、AMLの多面的な形質を制御しているかを網羅的に解明する必要がある。また、cordycepinの臨床応用における最適な投与量や、他の化学療法剤との併用効果についての検証が limitation として残されており、さらなる臨床試験での検証が待たれる。特に、患者個々のPAPOLA発現レベルに基づいた層別化治療の可能性について、さらなる検証が必要である。
方法
本研究では、AML患者由来のCD34+芽球および正常HSPCを用いて、FLAM-seq (full-length mRNA sequencing) によるポリA鎖長の網羅的解析を行った。FLAM-seqでは、poly(U)ポリメラーゼを用いてポリA末端にG/I tailを付加し、PacBio Sequelシステムによるロングリードシーケンシングを実施した。PAPOLAの機能解析には、shRNAを用いたノックダウンおよび、Papola fl/fl マウス、Vav-iCreマウス、Mx1-Creマウスを用いた条件付きノックアウトモデルを構築した。
AMLのin vivoモデルとして、MLL-AF9またはHoxa9-Meis1を導入したHSPCを致死性照射後のC57BL/6Jマウスに移植した。LSCの自己複製能は、二次移植後の限界希釈分析（LDA）およびELDAソフトウェアを用いて定量化した。代謝解析には、Agilent XF96 extracellular flux analyzerを用い、OCR（酸素消費速度）およびECAR（細胞外酸性化速度）を測定した。また、LC-MS (liquid chromatography-mass spectrometry) を用いてTCAサイクル中間体およびNAD+/NADH比を定量した。
統計解析には、GraphPad Prism 8を使用し、2群比較には unpaired two-tailed Student’s t-test、多群比較には one-way ANOVA または two-way ANOVA に続き Dunnett’s multiple-comparison test を適用した。生存分析には Kaplan-Meier 法を用い、log-rank test で有意差を検定した。予後因子の解析には Cox regression 分析を用いた。細胞株としては HL-60, U-937, HEK293T を使用し、マウス系統には C57BL/6J, NOD/SCID を用いた。