〜 おへその乱読1000本ノック #0058 ～ゲノムワイドDNAメチル化プロファイルからみるヒト老化速度の定量化

　第58回『おへその乱読1000本ノック』へようこそ。
　前回は、エピジェネティック老化指標の背後にある遺伝的基盤に注目し、複数のエピジェネティッククロックや関連バイオマーカーがどのような分子ネットワークと結びついているのかを俯瞰しました。今回からは、前回整理した6種類のエピジェネティックバイオマーカーを、それぞれ原著論文に立ち返って順に読み解いていきます。
　すでに第21回において、複数組織に適用可能なHorvathクロックを取り上げました。そこで今回は、それに続く2つめとして、Hannumら（2013）によって報告された、いわゆるHannumクロックを取り上げます。ここでは、全血のDNAメチル化データを基盤として構築された、血液ベースのエピジェネティッククロックに目を向けます。こうした対比は、エピジェネティック年齢という指標が、どの組織をもとに、何を反映しているのかを考えるうえでも重要です。
　Hannumクロックは、エピジェネティッククロック研究の初期を支えた代表的なモデルのひとつであり、血液中のDNAメチル化パターンから年齢を推定する枠組みを示しました。Horvathクロックと並んで、その後のさまざまな老化指標の基盤として参照される重要な研究です。
　本稿では、このHannumクロックを手がかりに、血液ベースのエピジェネティッククロックとはどのようなものか、そしてそれが老化研究のなかでどのような位置を占めるのかを見ていきます。

1. イントロダクション

1.1 老化の個人差と老化バイオマーカー探索

　ヒトの老化は一様に進行するものではなく、性差や生活習慣の違いによってその進み方には個人差が生じる。このため、加齢に伴う生理機能の低下や疾患リスクを予測・理解するための老化バイオマーカーが求められてきた。これまでにも、テロメア長や遺伝子発現が年齢と関連する指標として注目されてきたが、近年ではエピゲノム、特にDNAメチル化の状態も、加齢と深く関係することが明らかにされている。

1.2 DNAメチル化と加齢、エピジェネティックドリフト

　なかでもDNAメチル化は、長期的な加齢に伴って変化し、代謝疾患やがんなどの加齢関連疾患とも結びつくことが報告されてきた。また、一卵性双生児において加齢とともにメチル化状態の差が広がるエピジェネティックドリフトも知られており、エピゲノムは固定的なものではなく、時間とともに変動する動的な景観として捉えられるようになってきた。

1.3 老化メチローム研究における課題

　しかし、こうした加齢に伴うメチル化変化を体系的に記述し、個人ごとの老化速度の違いを定量的に比較できるかどうかは、当時なお大きな課題として残されていた。著者らは、老化メチロームの変化には、環境曝露に応答した比較的一貫した変化と、複製エラーや化学的攪乱などに由来するより偶発的な変化の両方が関与すると考え、こうした変化を網羅的に測定することで、加速あるいは遅延した老化に関わる要因を捉えられる可能性を提示している。

1.4 本研究の目的

　そこで本研究では、幅広い年齢層の大規模集団を対象にゲノムワイドなDNAメチル化解析を行い、ヒトの老化速度を定量化できる予測モデルの構築を目指した。そのうえで、得られた老化モデルが性差や遺伝的背景とどのように関係するのか、さらにエピジェネティックドリフトや転写変化とどのようにつながるのかを検討することが、本研究の基本的な目的として提示されている。

2. 加齢に伴うメチローム変化の全体像

2.1 大規模コホートにおけるDNAメチル化プロファイリング

　本研究では、19歳から101歳までの幅広い年齢層を含むヒト集団（計656名）を対象に、全血由来のDNAを用いたゲノムワイドなメチル化解析が実施された。解析には、約48万CpG部位のメチル化状態を測定可能なアレイが用いられ、各CpGについてメチル化率（0〜1）として定量化されている 。
　品質管理ののち、性染色体を除いたCpGを対象に解析が行われ、本データセットは当時として最大規模かつ高解像度の老化メチローム解析基盤となっている。

2.2 年齢と関連するCpG部位の同定

　各CpG部位について年齢との関連を統計的に評価した結果、全体の約15%にあたる70,387箇所で、メチル化率と年齢との有意な相関が認められた。これらの多くは独立データセットにおいても再現されており、加齢に伴うDNAメチル化変化が広範かつ再現性の高い現象であることが示されている 。
　これらの結果は、特定の遺伝子座に限定された変化ではなく、ゲノム全体にわたる体系的な変動として、加齢に伴うメチローム変化が生じていることを示唆する（Figure 1）。

2.3 加齢に伴うメチル化パターンの特徴

　年齢関連CpGの分布を俯瞰すると、メチル化率が上昇する部位と低下する部位の双方が存在し、加齢は単一方向の変化ではなく、複雑な再編成として進行することが示される。また、これらの変化はCpGアイランドやその周辺領域（shore, shelf）、エンハンサーやプロモーターなど、ゲノム機能領域ごとに偏りをもって分布しており、加齢に伴うメチル化変化が特定のゲノム文脈と結びついていることが示される（Table S2） 。

　さらに、これらの年齢関連CpGは、細胞間コミュニケーション、細胞増殖、発生・成長といった生理機能に関連する遺伝子近傍に多く見られ、機能的にも偏りを有していることが示されている（Table S1） 。

2.4 老化メチロームという視点

　以上の結果から、DNAメチル化状態は固定的なものではなく、加齢に伴って体系的に変化する動的な特徴量として捉えられる。このような全体的変化は、単一のバイオマーカーではなく、複数のCpGの組み合わせとして、老化メチロームとして記述されるべきものである。ここで示されたゲノムワイドな変化の存在は、後続の解析で行われる年齢予測モデルの構築や、老化の進行速度の定量化に向けた基盤となる。

3. 老化予測モデルの構築

3.1 モデル構築の基本戦略

　前章で示されたように、加齢に伴うDNAメチル化変化はゲノム全体に広く分布している。本研究では、これらの情報を統合することで、メチロームから年齢を定量的に予測するモデルの構築が試みられた。解析には、変数選択と回帰を同時に行うElastic Net法が用いられ、多数のCpGマーカーの中から年齢予測に寄与する特徴量が抽出された。さらに、ブートストラップ法を併用することで、モデルの安定性と再現性が担保されている 。

3.2 年齢予測モデルの精度と構成

　最終的に、71のCpGマーカーから構成される年齢予測モデルが構築された。このモデルにより、メチル化データから推定された年齢と実際の年齢との間には極めて高い相関（R ≈ 0.96）が認められ、平均誤差は約3.9年と報告されている 。さらに、独立した検証コホートにおいても同様に高い予測精度が確認されており、本モデルの汎用性と再現性が支持されている（Figure 2B, C）。

3.3 モデルに選択されたCpGの生物学的特徴

　モデルに含まれるCpGマーカーの多くは、老化や加齢関連疾患に関与する遺伝子近傍に位置していた。具体的には、神経内分泌系の調節に関与する遺伝子や、代謝・肥満に関連する転写因子などが含まれており、老化と代謝機能や疾患との結びつきを反映していると考えられる 。このことは、エピジェネティッククロックが単なる統計的モデルではなく、生理学的変化を背景とした分子指標として機能している可能性を示唆する。

3.4 モデルの位置づけ

　本研究で構築された年齢予測モデルは、DNAメチル化情報を統合することで暦年齢を高精度に再現できることを示した点で、エピジェネティッククロックの基本的枠組みを確立したものといえる。このモデルは、後に提案される多様なエピジェネティッククロックの原型の一つとして位置づけられ、単なる年齢推定にとどまらず、個体ごとの老化の進行速度を評価する基盤となる。

4. 老化の進行速度と個体差

4.1 老化の進行速度（AMAR）の定義

　前章で構築された年齢予測モデルは、単に暦年齢を推定するだけでなく、個体間の老化の進み方の違いを評価するための基盤としても利用できる。本研究では、メチル化データから予測された年齢を実年齢で割った値として、apparent methylomic aging rate（AMAR）が定義された。この指標は、メチロームの変化が暦年齢に対してどの程度進んでいるかを示すものであり、1を基準として、これより大きい場合は加速した老化、小さい場合は遅延した老化を意味する 。

4.2 老化速度の個体差と外れ値

　年齢予測モデルの結果を俯瞰すると、実年齢と予測年齢が大きく乖離する個体が存在することが確認される。これらの外れ値は単なる測定誤差ではなく、個体ごとの生物学的状態の違いを反映している可能性がある。AMARを用いることで、こうした個体差を定量的に評価することが可能となり、メチロームに基づく老化のばらつきが明確に可視化される。

4.3 性差と老化速度

　AMARと臨床因子との関連を解析した結果、性別が老化速度に有意に影響することが示された。具体的には、男性のメチロームは女性と比較して約4%速く老化する傾向が認められた（Figure 2D）。一方で、BMIについては有意な関連は認められず、少なくとも本解析条件においては、体格よりも性差の影響が強く現れていることが示唆される。

4.4 老化速度という概念の意義

　AMARの導入により、本研究は年齢を当てるモデルから、老化の進み方を測るモデルへと視点を拡張している。この点は、エピジェネティッククロック研究における重要な転換点である。すなわち、DNAメチル化は単なる時間の指標ではなく、個体ごとの生理的状態や老化の進行度を反映する指標として機能しうる。この概念は、後に提案されるフェノタイプ関連クロックや死亡リスク予測モデルへと発展していく基盤となる。

5. 遺伝的要因とメチローム老化

5.1 老化速度の個体差を規定する要因へのアプローチ

　前章で示したように、メチロームに基づく老化速度（AMAR）には個体差が存在する。このばらつきがどのような要因によって規定されているのかを明らかにするため、本研究では遺伝的変異とDNAメチル化との関連に着目した。特に、加齢関連CpGのメチル化状態に影響を与える遺伝的変異、すなわちmethylation quantitative trait loci（meQTL）の同定が試みられている 。

5.2 meQTLの同定とその特徴

　エクソームシーケンスにより同定された遺伝的変異と、加齢関連CpGのメチル化状態との関連解析の結果、303のmeQTLが同定された。これらの多くは、対応するCpG近傍に位置するcis作用型の変異であり、特定の遺伝子領域においてメチル化状態を調節していることが示唆される（Figure 3A）。さらに、これらの一部は独立コホートにおいても再現されており、遺伝的背景がメチロームの加齢変化に影響を与えることが支持されている。

5.3 老化関連遺伝子とメチル化制御

　同定されたmeQTLの中には、DNA修復や細胞周期制御、代謝など、老化や加齢関連疾患に関与する遺伝子が含まれていた。例えば、MBD4はメチル化DNAに結合するタンパク質をコードし、ゲノム安定性の維持に関与することが知られている。また、NEK4やGTPBP10といった遺伝子も、それぞれ細胞周期制御やミトコンドリア機能と関連している。これらの結果は、加齢に伴うメチル化変化が単なる受動的な変動ではなく、遺伝的に制御されたプロセスの一部である可能性を示唆する（Figure 3B, C）。

5.4 遺伝的変異と老化速度の関連

　さらに興味深いことに、一部のmeQTLは単にメチル化状態と関連するだけでなく、老化速度（AMAR）そのものとも有意に関連していた。すなわち、特定の遺伝的変異が、メチロームの加齢変化の速度に影響を与える可能性が示された。この結果は、老化速度の個体差が、少なくとも一部は遺伝的要因によって規定されていることを示すものである。

6. 多組織における適用可能性と腫瘍における老化加速

6.1 血液ベースモデルの他組織への適用

　本研究で構築された老化予測モデルは全血由来のデータに基づくものであるが、その適用可能性を検証するため、他組織におけるDNAメチル化データへの適用が試みられた。解析には、The Cancer Genome Atlas（TCGA）に含まれる正常組織（乳腺、腎臓、肺、皮膚）のデータが用いられている 。
　その結果、いずれの組織においても、予測年齢と実年齢の間には有意な相関が認められた（Figure 4A）。一方で、組織ごとに回帰直線の傾きや切片に系統的なずれが存在しており、血液ベースモデルは複数組織に共通する加齢シグナルを一定程度捉える一方で、組織特異的な差も反映していることが示された。

6.2 組織特異的な補正と予測精度

　こうした組織ごとの差に対して、線形補正（傾きおよび切片の補正）を行うことで、各組織における年齢予測誤差は血液と同程度まで改善された（Figure 4B）。この結果は、血液由来のメチロームに基づく老化モデルが、そのままでは組織差の影響を受けるものの、適切な補正を加えることで他組織にも適用しうることを示している。

6.3 腫瘍組織における老化の加速

　さらに、正常組織と対応する腫瘍組織の比較が行われた結果、腫瘍では予測年齢が有意に高く、見かけ上の老化速度（AMAR）が約40%増加していることが示された（Figure 4C, D）。この傾向は特定の組織に限らず、複数の組織で一貫して観察されており、腫瘍においてはメチロームの加齢様変化が顕著に進行していることが示唆される。

7. エピジェネティックドリフトとメチロームの変動性

7.1 加齢に伴うメチル化変動の増大

　本研究では、加齢に伴うDNAメチル化の平均値の変化だけでなく、個体間におけるメチル化のばらつき（deviance）にも着目している。各CpGについて、期待される平均メチル化値からの偏差（二乗距離）を指標として解析した結果、27,800箇所において、このばらつきが年齢と有意に関連することが示された 。
　特に注目すべき点として、そのほとんど（99.8%）が加齢に伴うばらつきの増加として観察されており、メチロームは年齢とともにより不均一な状態へと変化していくことが示されている（Figure 5）。

7.2 エピジェネティックドリフトの定量化

　このようなメチル化のばらつきの増大は、「エピジェネティックドリフト」として知られる現象と対応するものである。本研究では、このドリフトを定量的に評価する指標として、メチル化の分散や偏差に基づく解析が導入されている。さらに、個体ごとのメチローム全体のばらつきの大きさは、老化速度（AMAR）とも正の相関を示しており、メチル化の不均一性の増大が老化の進行と結びついている可能性が示唆される（Figure 5） 。

7.3 メチロームのエントロピー増加

　加齢に伴うメチル化変化は、情報理論的な観点からも評価されている。すなわち、各CpGのメチル化率が極端な値（0または1）から中間値（約0.5）へと近づくことで、メチローム全体のエントロピー（不確実性）が増大することが示された。実際に、加齢に伴ってメチル化のエントロピーが有意に増加し、このエントロピー指標もまた老化速度（AMAR）と強く関連することが確認されている（Figure 6） 。

8. トランスクリプトームとの関連

8.1 メチローム変化と遺伝子発現の対応

　本研究では、DNAメチル化の加齢変化が遺伝子発現にどのように対応するかを検討するため、同一個体におけるトランスクリプトーム解析が行われた。加齢に伴ってメチル化が変化するCpG近傍の遺伝子について、その発現量との関連が評価されている 。その結果、メチル化変化と遺伝子発現との間には有意な関連が認められ、エピゲノムの変化が転写レベルの変動と結びついていることが示された（Figure 7A）。

8.2 発現変化に基づく老化シグナル

　さらに、遺伝子発現データのみを用いた解析においても、加齢に伴う体系的な変化が観察された。すなわち、トランスクリプトームにおいても、年齢に関連する発現変動が広範に存在し、これらの情報から老化に関する特徴量を抽出できることが示されている。この結果は、DNAメチル化に基づく老化指標が、単なるエピジェネティックな変化ではなく、転写レベルの変化と整合的な生物学的現象を反映していることを支持する（Figure 7B）。

8.3 メチロームとトランスクリプトームの統合的理解

　以上の結果から、加齢に伴うDNAメチル化変化は、遺伝子発現の変動と連動しており、機能的な側面を伴う変化であることが示された。すなわち、メチロームは単なる時間の記録ではなく、細胞機能の変化を反映する分子基盤の一部として位置づけられる。

9. 結論

　本研究は、ヒト全血におけるDNAメチル化のゲノムワイド解析を通じて、加齢に伴うメチローム変化の体系的な特徴を明らかにした。メチル化状態は年齢とともに広範に変化し、その情報を統合することで暦年齢を高精度に予測可能であることが示された。
　さらに、メチロームに基づく年齢推定は、単なる年齢の再現にとどまらず、個体ごとの老化の進行速度（AMAR）を定量化する指標として機能することが示された。この老化速度には個体差が存在し、その一部は遺伝的変異によって規定されていることが示唆される。
　また、本モデルは血液に限らず複数の組織に適用可能であり、腫瘍組織においては老化の加速が観察されるなど、メチローム変化が生理的状態や疾患と関連することが示された。
　加えて、加齢に伴うメチル化変化は平均値の変動だけでなく、ばらつきの増加やエントロピーの上昇としても現れ、メチロームの情報的な秩序が徐々に失われていく過程が示された。これらの変化は老化速度とも関連しており、老化が単なる時間依存的な変化ではなく、系のゆらぎの増大を伴う現象であることを示唆する。
　最後に、これらのメチローム変化は遺伝子発現の変動とも対応しており、エピジェネティックな変化が機能的な側面と結びついていることが確認された。