- 著者: H. Lalchungnunga, Hande Atasoy, Edward C. Schwalbe, Chris M. Bacon, Gordon Strathdee
- Corresponding author: Gordon Strathdee (Newcastle University Centre for Cancer, Biosciences Institute, Newcastle University, Newcastle, UK)
- 雑誌: British Journal of Cancer
- 発行年: 2026
- Epub日: 2025-02-01
- Article種別: Original Article
- PMID: 41174288
背景
がんにおけるエピゲノムプロファイルの再構成、特にDNAメチル化異常は、腫瘍の発生や進行に関与する普遍的な特徴として広く認識されている。これまでのゲノムワイドなメチル化解析により、B細胞悪性腫瘍をはじめとする血液がんにおいて、正常なB細胞前駆細胞と比較して数千もの差次的メチル化領域 (DMR: differentially methylated region) が同定されてきた。しかし、これらの膨大なメチル化変化のうち、実際にがん化のドライバーとして機能している領域と、単なる随伴的な変化 (パッセンジャー変化) を区別することは極めて困難であった。
先行研究において、Kondo et al. (2010) はがんにおけるDNAメチル化プロファイリングの重要性を指摘し、Kulis et al. (2010) はDNAメチル化と発がんの密接な関連を報告している。さらに、Schwalbe et al. (2021) は急性リンパ性白血病 (ALL: acute lymphoblastic leukaemia) におけるメチル化変化が、正常な記憶B細胞 (memory B-cell) で観察されるメチル化変化と高度に重複していることを示し、両者の間に極めて高い相関が存在することを明らかにした。記憶B細胞は長命かつ高増殖性の細胞であり、正常細胞プロセスとしての増殖が、がん細胞と類似したメチル化パターンを誘導する主要な要因である可能性が示唆されていた。
従来の解析手法では、がん細胞と正常な前駆細胞を単純に比較するアプローチが主流であったため、正常な細胞増殖や分化に伴う生理的なメチル化変化を「がん特異的な異常」として誤同定するリスクが常に存在していた。このため、真の疾患特異的なメチル化変化を正確に分離・同定する手法が確立されておらず、血液がんにおけるDNAメチル化異常の機能的意義の解明を阻む大きな知識ギャップ (knowledge gap) となっていた。また、液体生検 (liquid biopsy) などの臨床応用においても、増殖細胞由来の背景ノイズを排除し、真の疾患特異的マーカーを同定するための情報が不足しており、生理的変化を排除した高精度なマッピング手法の開発が課題として残されていた。本研究は、この課題を解決するために、正常細胞の生理的変化 (増殖・分化) を参照点として組み込んだ新しいバイオインフォマティクス手法を開発し、真の疾患関連メチル化変化を抽出することを試みた。
目的
本研究の目的は、B細胞悪性腫瘍における全DNAメチル化変化を、増殖由来、分化関連、および真の疾患特異的な変化の3つの主要な成分に統合的に分解する手法「統合的メチル化マッピング (integrative methylation mapping)」を開発・確立することである。この新規アプローチを用いて、5種類の主要なB細胞悪性腫瘍における疾患特異的DMRの正確な割合を定量化し、その分子生物学的およびエピゲノム的な特徴を解明することを目指す。
さらに、このマッピング手法によって絞り込まれた真の疾患特異的DMRの中から、腫瘍抑制遺伝子 (TSG: tumour suppressor gene) として機能する重要な候補遺伝子を同定し、複数のALL細胞株を用いた体外実験を通じてその生物学的機能を検証する。特に、新規腫瘍抑制遺伝子候補であるSLC22A15 (solute carrier family 22 member 15) に焦点を当て、ALL細胞の生存、アポトーシス誘導、および増殖速度に与える影響を機能的に実証することで、将来的な診断バイオマーカーや新規治療標的としての有用性を評価することを目的とする。
結果
B細胞悪性腫瘍におけるメチル化変化の大部分が増殖シグナルに起因: B細胞前駆細胞をベースラインとしたDMR解析 (beta値差異>0.2カットオフ) の結果、全B細胞悪性腫瘍を統合した解析で同定された総DMR数 5,692 のうち、増殖関連DMRが 5,004 (87.9%) と圧倒的多数を占めることが判明した (Table 1)。これに対し、分化関連DMRは 514 (9.0%)、真の疾患特異的DMRはわずか 156 (2.7%)、疾患不在DMRは 18 (0.3%) に過ぎなかった。beta値差異のカットオフを0.1に緩和した解析でも、総DMR数は大幅に増加したものの、疾患特異的DMRの割合は 1.7% (224 DMR) とさらに低下し、増殖関連DMRが 66.4% (8,628 DMR) と依然として支配的であった (Table 1)。全B細胞悪性腫瘍におけるメチル化変化の全体像は、正常な記憶B細胞における変化と極めて高い相関を示し、Pearson相関係数は r=0.94 (p<0.0001) に達した (Figure 1b, c)。この結果は、従来の単純な正常-がん比較で同定されてきたメチル化異常の大部分が、がん化そのものではなく、正常細胞の生理的な増殖プロセスを反映しているに過ぎないことを示している。
DMRグループごとの異なるエピゲノムおよび分子特徴: 各DMRグループの配列コンテキストおよびクロマチン状態をSeSAMe解析により評価したところ、明確な分子特徴の差異が認められた (Figure 2)。増殖関連DMRは、平均CpGサイト数が 7.29/DMR と最も高く、その 34.4% が転写開始点 (TSS) から 1,000 bp 以内のプロモーター領域に位置していた (Table 2)。また、増殖関連DMRでは高メチル化 (53.4%) と低メチル化 (46.6%) が同程度に観察された。一方、分化関連DMRは著しく低メチル化に偏っており、高メチル化はわずか 5.1% であった (Table 2)。疾患特異的DMRは中間の特性を示し、19.9% が高メチル化、平均CpGサイト数は 4.35/DMR であった。ヒストン修飾の解析において、増殖関連および分化関連の高メチル化DMRは、ポリコーム複合体 (PRC2) 標的領域やH3K27me3およびH3K9me3修飾領域と強い相関を示した。これに対し、疾患特異的DMRはH3K27me3領域とは相関したものの、H3K9me3領域との相関を欠いており、メチル化が導入されるエピゲノム制御機構に違いがあることが示唆された (Figure 2)。
個別疾患における疾患特異的DMRの著しい少なさ: 5種類の悪性腫瘍を個別に解析したところ、疾患特異的DMRの割合はさらに極少であることが明らかになった (Table 3)。CLLでは総DMR 6,326 のうち疾患特異的DMRはわずか 0.10% (6 DMR)、MCLでは総DMR 6,659 のうち 0.15% (10 DMR)、DLBCLでは総DMR 11,701 のうち 0.15% (17 DMR) であり、これら3つの成熟B細胞腫瘍においては、真の疾患特異的メチル化変化は事実上ほとんど存在しなかった。一方、ALLでは 4.80% (240/4,988 DMR)、PCNSLでは 5.40% (1,104/20,583 DMR) と比較的高い割合を示した (Table 3)。しかし、他疾患での変化が対象疾患の4分の1以下というより厳格な特異性フィルターを適用すると、疾患特異的DMRの割合はPCNSLで 1.5%、ALLで 3.4% にまで低下した。PCNSLにおけるDMRの絶対数の多さは、真の疾患特異性よりも、この疾患におけるメチル化変化の絶対的な振幅の大きさを反映していると考えられる。
ZAP70の汎B細胞悪性腫瘍マーカーとしての同定: 疾患特異的DMRのうち、近傍遺伝子の発現と有意な負の相関を示す領域を探索したところ、174の候補領域から10遺伝子座が同定された。このうち、すべてのB細胞悪性腫瘍においてbeta値差異>0.2を満たし、かつTSS近傍のプロモーター領域に位置する真の共通標的はZAP70 (zeta chain of T-cell receptor associated protein kinase 70) のみであった (Figure 1d)。ZAP70はCLLにおける予後因子として知られているが、本解析結果は、ZAP70のメチル化抑制がすべてのB細胞悪性腫瘍に共通する重要な疾患特異的イベントであることを示している。
SLC22A15の新規ALL腫瘍抑制遺伝子としての同定と機能検証: ALLに特異的な疾患関連DMRから、遺伝子発現と負の相関を示す腫瘍抑制遺伝子候補として THEM4、MAP9、SLC22A15、TTC12 の4遺伝子を同定した。機能検証実験において、TTC12の再発現は細胞生存に影響を与えず、MAP9は導入後の発現が確認できなかった。THEM4は t(1;19) サブタイプであるPreB697細胞でのみ選択的に細胞消失を誘導した。
これに対し、有機溶質トランスポーターをコードするSLC22A15をレンチウイルスを用いて5種のALL細胞株 (REH, PreB697, RCH-ACV, NALM-6, KOPN-8) に再発現させたところ、極めて顕著な抗腫瘍効果が観察された。遺伝子導入の検証にあたり、n=3 cells (独立した細胞培養バッチ) のスケールで再現性を確保し、導入後の細胞生存率を経時的に測定した。遺伝子導入後8日目 (day 8) にはウェスタンブロッティングによりSLC22A15タンパク質の明らかな発現が確認されたが、SLC22A15陽性 (eGFP陽性) 細胞の割合はすべての細胞株において経時的に著しく減少した (Figure 3a, b)。day 29までには、培養集団中からSLC22A15陽性細胞がほぼ完全に消失した。この細胞消失のメカニズムを解明するため、Caspase 3/7活性を測定したところ、SLC22A15の再発現後3日目 (day 3) から顕著な活性化が認められ、対照群 (empty vector) と比較して3.5-fold increase以上の有意な上昇を示し、day 4からday 5にかけてピークに達した (p<0.001, Figure 3c)。一方で、eFluor 450を用いた細胞増殖アッセイでは、対照群と比較して細胞分裂速度に有意な変化は認められず、log2FC 0.05 未満の差にとどまった (Figure 3d)。以上の結果から、SLC22A15の再発現はALL細胞の増殖を停止させるのではなく、強力なアポトーシス (細胞死) を誘導することで、細胞生存を阻害する腫瘍抑制遺伝子として機能することが実証された。
考察/結論
先行研究との違い: 本研究の成果は、B細胞悪性腫瘍におけるDNAメチル化変化の 87.9% が正常細胞の増殖プロセスに由来し、真の疾患特異的な変化はわずか 2.7% に過ぎないことを定量的に示した点にある。この知見は、従来の単純ながんと正常前駆細胞の比較によって同定された数千のDMRを「がん特異的変化」とみなしてきた多くの先行研究の解釈と対照的である。例えば、Duran-Ferrer et al. (2020) はCLLにおける増殖歴がメチル化パターンを規定することを示したが、本研究はこれを5種類のB細胞悪性腫瘍に拡張し、生理的な増殖シグナルががんのエピゲノム景観の大部分を支配していることを明確に証明した。
新規性: 本研究で初めて開発された「統合的メチル化マッピング」は、正常記憶B細胞の生理的メチル化プロファイルを明示的な参照点として用いることで、パッセンジャー変化を排除する極めてユニークなアプローチである。本研究で初めて、これまで血液がんにおける関与が報告されていない新規の腫瘍抑制遺伝子としてSLC22A15を同定し、その機能を複数のALL細胞遺伝学的サブタイプにおいて実験的に実証することに初めて成功した。SLC22A15はカルノシンなどのツビッターイオンを輸送するトランスポーターであり、その再発現がALL細胞に特異的にアポトーシスを誘導するという発見は、本研究が初めて提示する新規の知見である。
臨床応用: 本研究の知見は、無細胞DNA (cfDNA: cell-free DNA) メチル化プロファイルを利用した液体生検や微小残存病変 (MRD: minimal residual disease) モニタリングの臨床応用に直結する。従来のマーカー選定では、増殖細胞由来の生理的なメチル化変化が背景ノイズとなり、偽陽性を生む原因となっていた。本研究で絞り込まれた 156 の真の疾患特異的DMR (SLC22A15プロモーターやZAP70など) に標的を限定することで、診断アッセイの特異度と感度を大幅に向上させることが可能となる。また、SLC22A15の基質であるカルノシンは抗腫瘍活性を持つ代謝物として知られており、カルノシン代謝経路を標的とした新たな治療戦略の確立など、臨床現場への translational な貢献が期待される。
残された課題: 今後の検討課題として、いくつかの制限事項 (limitation) が挙げられる。第一に、本研究で使用したデータはすべてIllumina 450kアレイに基づくものであり、ゲノム全体のCpGサイトの一部しかカバーしていないため、エンハンサー領域などの疾患特異的変化を見落としている可能性がある。第二に、従来の重亜硫酸塩 (ビスルファイト) 処理では5-メチルシトシンと5-ヒドロキシメチルシトシンを区別できないという技術的限界が存在する。今後は、全ゲノムメチル化シーケンス (WGBS) を用いた高解像度な検証や、SLC22A15の具体的な基質輸送メカニズムの解明、さらにはT細胞腫瘍や骨髄系腫瘍など他の血液悪性腫瘍への本手法の適用拡大が必要である。
方法
データ収集と前処理: 本研究では、公開データベースであるGEO (Gene Expression Omnibus) およびEGA (European Genome-phenome Archive) から、合計n=995検体のIllumina Human Methylation 450kアレイデータを取得した。対象としたB細胞悪性腫瘍は、ALL (n=517)、慢性リンパ性白血病 (CLL: chronic lymphocytic leukaemia, n=187)、マントル細胞リンパ腫 (MCL: mantle cell lymphoma, n=86)、びまん性大細胞型B細胞リンパ腫 (DLBCL: diffuse large B-cell lymphoma, n=79)、原発性中枢神経系リンパ腫 (PCNSL: primary central nervous system lymphoma, n=95) の5種類である。正常対照群として、B細胞前駆細胞 (n=22) および成熟B細胞集団 (naive B細胞、クラススイッチ記憶B細胞など、n=9) のデータを使用した。また、対応する各集団の転写物データ (ALL n=326、CLL n=143、MCL n=122、DLBCL n=203、PCNSL n=34、正常前駆細胞 n=31、正常成熟B細胞 n=20) も取得し、統合解析に用いた。
統合的メチル化マッピングの構築: Rソフトウェア (v3.4.0) を用いてバイオインフォマティクス解析を実施した。DMRの同定にはDMRcate (differentially methylated region detection tool) パッケージを使用し、平均beta値差異のカットオフを0.2 (一部の解析では0.1) 以上、DMR内の最小CpGサイト数を2以上、有意確率をp<0.0001と定義した。B細胞前駆細胞を基準の「ベースライン」として設定し、(1) 全B細胞悪性腫瘍 vs B細胞前駆細胞、(2) 正常記憶B細胞 vs B細胞前駆細胞の2つの比較から得られたDMRセットを比較した。これにより、DMRを「増殖関連」(悪性腫瘍と記憶B細胞で共通して変化)、「分化関連」(成熟B細胞悪性腫瘍と記憶B細胞で変化し、ALLでは不変)、「疾患特異的」(悪性腫瘍でのみ変化し、記憶B細胞では不変)、「疾患不在」(記憶B細胞でのみ変化し、悪性腫瘍では不変) の4つのカテゴリに分類した。
エピゲノム特徴解析: 同定された各DMRグループのゲノムコンテキスト、クロマチン状態、ヒストン修飾、および転写因子結合部位 (TFBS: transcription factor binding site) の濃縮度を評価するため、SeSAMe (sensible step-by-step analysis of DNA methylation data) パッケージを使用した。Illumina 450kアレイの全プローブを背景対照とし、Fisher’s exact (フィッシャー正確検定) を用いて統計的有意性を算出した。また、メチル化レベルと遺伝子発現の相関分析には Pearson correlation (ピアソン相関係数) を用いた。
機能的検証実験: 同定された腫瘍抑制遺伝子候補である THEM4 (thioesterase superfamily member 4)、MAP9 (microtubule-associated protein 9)、SLC22A15、TTC12 (tetratricopeptide repeat domain 12) を、ゲートウェイテクノロジーを用いてpSINE-SIEWレンチウイルスベクターにサブクローニングした。このベクターは、同一転写産物から内部リボソーム進入部位 (IRES: internal ribosome entry site) を介してeGFP (enhanced green fluorescent protein) を同時に発現する。構築したレンチウイルスを用いて、異なる細胞遺伝学的サブタイプを代表する5種のヒトALL細胞株 (REH, PreB697, RCH-ACV, NALM-6, KOPN-8) に遺伝子を導入した。また、パッケージング細胞として HEK293T 細胞株を使用し、トランスフェクションを実施した。遺伝子導入後の細胞生存率は、eGFP陽性細胞の割合をフローサイトメトリー (FACS Fortessa X-20) で経時的に追跡することで評価した。アポトーシスは、Annexin V/PI染色アッセイおよびCaspase-Glo 3/7 Assayを用いて検証した。細胞増殖速度は、eBioscience Cell Proliferation Dye eFluor 450を用いて測定した。タンパク質の発現確認は、抗SLC22A15抗体および抗THEM4抗体を用いたウェスタンブロッティング法により実施し、β-actinを内部標準とした。統計解析には Student t-test (学生のt検定) を用いた。