〜 おへその乱読1000本ノック #0038 ～ エイジングクロック、エントロピー、そして年齢若返りの可能性と限界

　老化研究の分野では近年、エイジングクロックと呼ばれる指標が大きな注目を集めています。DNAメチル化や血液データ、生活習慣に関する情報などをもとに、生物学的な年齢を推定しようとする試みです。これらの指標は、健康状態や寿命の予測にも関係することが示され、老化研究や医療の現場でも活用が広がりつつあります。
　しかし一方で、こうしたエイジングクロックが老化の本質をどこまで捉えているのかについては、まだ十分に理解されているとは言えません。老化とは、単に時間が経過することなのでしょうか。それとも、生体の中で起こる無数の変化が積み重なった結果なのでしょうか。
　本論文は、この問題を少し異なる視点から捉えます。老化を、生体の中で起こる多数の小さな状態変化が確率的に積み重なっていく過程として理解し、それをエントロピーという物理学の概念と結びつけて説明しようとするものです。エイジングクロックが示す数値の背後には、こうした統計的な変化の蓄積があるのかもしれません。
　ここでは、エイジングクロック、エントロピー、そして若返りの可能性と限界という視点から、この論文の考え方を整理していきます。

1. イントロダクション

1.1 生物学的年齢という指標の登場

　老化は多くの生物学的階層に現れる複雑な過程であり、慢性疾患の発症率や死亡率が加齢とともに指数関数的に増加することが知られている。このような多様な表現型変化を統合的に捉える試みとして、生物学的年齢（biological age）という概念が提案されてきた。生物学的年齢モデルは、DNAメチル化、身体活動データ、血液検査などの生物医学データから、個体の年齢や死亡リスクを予測するように構築される。これらの指標は、全死亡率や慢性疾患の発症率・重症度と関連することが示されており、近年では臨床試験においても利用されつつある。

1.2 生物学的年齢と老化過程の関係

　しかし、生物学的年齢のダイナミクスや、それが老化そのものとどのように関係しているのかについては、依然として十分には理解されていない。
　例えば、DNAメチル化年齢は、老化がほとんど進行しないとされる種においても増加する場合がある。また、人間においても、生物学的年齢は概日リズムやストレス、喫煙などの生活習慣の影響を受け、一時的に変動することが知られている。
　さらに、加齢に伴って生体状態が恒常性へ回復するまでの時間や、生物学的年齢の揺らぎの大きさは増加する。回復が遅い個体の割合は約8年ごとに倍増することが報告されており、この時間スケールは人間の死亡率倍増時間と近い値である。このことは、生体状態の変動と死亡率との関係を理解する必要性を示唆している。
　今後、生物学的年齢指標が捉える変動と、死亡率・疾患・介入効果との因果的な関係を理解することが、老化研究と医療応用の両面で重要な課題となる。

1.3 大規模データ解析から見える老化の共通パターン

　本研究では、生物医学データに共通して現れる老化シグネチャーの特徴を明らかにすることを目的とした。具体的には、白血球DNAメチル化データおよびUK Biobankに収録された電子医療記録の大規模データセットに対して主成分解析を行った。その結果、老化に伴う変化は、多数のメタ安定状態の間で起こる希少かつ独立な確率的遷移として説明できることが示された。また、データの分散の大部分は、年齢とともに線形に増加する単一の要因によって説明され、この要因はDNAメチル化年齢や慢性疾患数と強く相関していた。

1.4 複雑な調節ネットワークとしての老化

　これらの観察結果を説明するために、本研究では複雑な調節ネットワークにおける老化の半定量的モデルを提案する。生体は多数の機能単位から構成され、それらは発生の終了時にメタ安定状態に初期化されると仮定する。老化は、これらの機能単位が確率的な状態遷移を繰り返すことで、生体全体が熱力学的平衡へと緩和していく過程として理解される。すなわち、老化は多数の微視的状態変化の蓄積として現れる現象である。

1.5 熱力学的生物学的年齢（tBA）という概念

　このような微視的状態変化の累積効果は、平均値と分散が時間とともに増加する確率変数として記述できる。本研究では、この量を熱力学的生物学的年齢（thermodynamic biological age, tBA）と定義する。
　tBAは年齢とともにほぼ線形に増加し、老化過程で生成されるエントロピー、すなわち情報の喪失量に比例すると考えられる。この視点から見ると、老化は生体状態の不可逆的ドリフトを引き起こし、生理的回復力の低下や慢性疾患発症率および死亡率の指数関数的増加につながる現象として理解できる。さらに、このような老化のエントロピー的性質は、完全な年齢逆転の可能性に制約を与える可能性がある。

2. 理論モデル

　　  老化を多数の機能単位の確率的状態遷移として捉える枠組み

　本研究では、老化を複雑な生体調節ネットワークのダイナミクスとして理解するための理論モデルを提案する。著者らは、生体を多数の機能単位（functional units; FUs）から構成されるシステムとして捉える。機能単位とは、DNAメチル化状態、疾病状態、分子構造変化など、生体内のさまざまな生理的状態を表す最小単位である。
　各機能単位は複数の状態をとり得るが、発生の終了時には特定の状態に安定化していると仮定される。この状態は完全な平衡状態ではなく、メタ安定状態である。したがって、時間の経過とともに確率的な揺らぎによって状態遷移が生じる可能性がある。

　この考え方を模式的に示したのが図3Aである。機能単位はエネルギーランドスケープのポテンシャル井戸に対応し、二つの状態（σ=0またはσ=1）として表現される。通常は活性化エネルギー障壁によって保護されているため状態遷移は稀であるが、熱的ゆらぎや環境ノイズによって時折遷移が生じる。
　生体には膨大な数の機能単位が存在するため、個々の遷移は稀であっても、系全体では多数の状態変化が累積していく。図3Bは、このような機能単位の遷移頻度の分布を概念的に示している。大部分の機能単位では遷移は非常に稀であり、個体の一生の間に一度も起こらない場合も多い。しかし、機能単位の総数が極めて多いため、全体としては遷移の総数は徐々に増加する。
　このとき重要な量となるのが、時間tまでに生じた状態遷移の総数 Z(t) である。個々の機能単位の遷移確率は小さいが、機能単位の数が非常に多いため、Z(t)は時間とともにほぼ線形に増加する。この累積変化が、生体の老化状態を表す基本的な変数であると著者らは考える。
　さらに、生体内の機能単位は独立して存在するわけではなく、相互作用によって結びついたネットワークを形成している。これらの相互作用は、調節場（regulatory fields）と呼ばれる変数によって表現される。調節場は、他の機能単位の状態や外部要因（生活習慣や環境ストレスなど）によって変化し、生体状態のダイナミクスを決定する。
　このネットワークにおいて、確率的状態遷移の累積は生体の調節構造に徐々に影響を与え、生理的回復力（resilience）の低下を引き起こす。図3Cはこの過程を概念的に示している。若い個体では安定したポテンシャル井戸に対応する状態にあるが、状態遷移が蓄積するとエネルギーランドスケープが徐々に変形し、生体状態の揺らぎが増大する。
　このように、本モデルでは老化を多数の微視的状態変化の確率的蓄積として理解する。個々の変化は小さく、短期的には生体機能に大きな影響を与えないが、長期的にはそれらの累積が生体状態を徐々に変化させる。この累積量 Z(t) が後に導入される熱力学的生物学的年齢（tBA）の基礎となる。

3. DNAメチル化データの解析

　　  老化シグネチャーの主要成分は年齢とともに線形に増加する

　老化に伴う分子変化の統計的特徴を調べるため、本研究では白血球のDNAメチル化データセットを解析した。DNAメチル化は加齢に伴って系統的に変化することが知られており、生物学的年齢を推定するエピジェネティッククロックの基盤として広く利用されている。

　本研究では、年齢と有意に相関するCpGサイトを選択したうえで主成分解析（PCA）を行い、老化に伴う主要な変動パターンを抽出した。その結果、最も大きな分散を説明する第一主成分（DNAm-PC1）が、年齢とともにほぼ線形に増加することが示された（図1A）。また、この第一主成分の分散も年齢とともに線形に増加することが確認された（図1B）。
　このような平均値と分散の双方が時間とともに増加する挙動は、確率的な拡散過程に特徴的であり、DNAメチル化変化が多数の独立したランダムな状態変化の蓄積として説明できる可能性を示唆する。さらに、この主成分はHorvathのDNAメチル化年齢と強い相関を示し、エピジェネティッククロックが捉えている老化シグネチャーと密接に関係していることが確認された。
　一方、第三主成分（DNAm-PC3）も年齢と相関する変動として同定されたが、その増加は線形ではなく、年齢とともに加速する非線形的な挙動を示した（図1C）。また、この主成分の分散は年齢とともに急速に増大し、その逆数を外挿すると120〜150歳付近でゼロに近づくことが示された（図1D）。これは、生体状態の揺らぎが高齢期に著しく増大する可能性を示唆している。
　さらに、主成分の負荷量の分布を調べたところ、DNAm-PC1は互いに逆方向に変化する二つの大きなCpGサイト群から構成されていることが明らかとなった（図1E）。すなわち、加齢に伴って一部のCpGサイトはメチル化が増加し、別のCpGサイトではメチル化が減少する。このことから、第一主成分はDNAメチル化状態の「分極」の変化、すなわち多数の状態遷移の総数を反映していると解釈される。
　また、DNAm-PC3に強く寄与するCpGサイトの遺伝子領域を解析した結果、自然免疫応答やがん関連経路といった生物学的プロセスが富化していることが示された（図1G, H）。この結果は、老化に伴う特定の生理機能の変化が、第二の老化シグネチャーとして現れている可能性を示している。
　以上の結果から、DNAメチル化データにおける主要な老化シグネチャーは、多数の独立した状態変化の累積として現れる統計的パターンを示すことが明らかとなった。この観察は、老化を確率的な状態遷移の蓄積として理解する理論モデルと整合的である。

4. 医療記録データの解析

　　  慢性疾患の累積も同様の統計構造を持つ

　DNAメチル化データにおいて観察された老化の統計的特徴が一般的な現象であるかを検証するため、本研究ではUK Biobankに収録された電子医療記録（electronic medical records; EMR）を解析した。解析対象には、慢性疾患として分類される111種類の診断項目が用いられた。
　各個体は、これらの疾患の有無を示す二値ベクトルによって表現される。若年期の多くの個体では慢性疾患がほとんど存在しないため、これらの状態は初期状態として強く偏極していると考えられる。しかし加齢とともに、特定の疾患が発症することで状態が変化していく。

　この疾患状態ベクトルに対してPCAを行った結果、DNAメチル化データと同様に、年齢とともに変化する主要な成分が存在することが明らかとなった。第一主成分（EMR-PC1）は年齢とともにほぼ線形に増加し（図2A）、慢性疾患の総数と強く相関することが確認された（図2B）。すなわち、この主成分は個体が経験した疾患状態の変化の累積を反映していると解釈される。
　また、この主成分の分散は年齢とともに増加し（図2C）、DNAメチル化データで観察されたものと同様に確率的拡散過程の特徴を示した。さらに、縦断データを用いて自己相関関数を解析した結果、主成分の変動は時間遅れに対して線形に増加することが示された（図2D）。この挙動は、ドリフトを伴う確率過程に典型的な性質であり、老化に伴う健康状態の変化が多数の独立した状態変化の蓄積として説明できることを示唆している。
　これらの結果は、老化に伴う統計的パターンがDNAメチル化のような分子レベルの指標だけでなく、慢性疾患の発症という臨床的指標にも共通して現れることを示している。したがって、生体における老化の主要なシグナルは、さまざまな生物医学データに共通する普遍的な統計構造を持つ可能性がある。

5. 熱力学的生物学的年齢（tBA）

　　  老化の累積変化を表す単一変数として定義

　DNAメチル化データおよび電子医療記録データの解析から、老化に伴う主要な変動は、年齢とともに増加する単一の統計的要因によって記述できることが示された。この結果を説明するため、本研究では、機能単位の状態遷移の累積を表す基本変数 Z(t) を導入する。
　各機能単位の状態遷移は個別には稀であるが、生体には膨大な数の機能単位が存在するため、全体としては遷移の総数が時間とともに増加する。この遷移の累積数 Z(t) は、主成分解析で得られた第一主成分と対応し、DNAメチル化変化や慢性疾患の発症数の増加として観測される。
　このような状態遷移の累積は確率過程として振る舞い、その平均値と分散は時間とともに増加する。したがって、Z(t) は老化の進行を表す自然な指標とみなすことができる。本研究では、この量を 熱力学的生物学的年齢（thermodynamic biological age; tBA） と定義する。

　さらに、状態遷移の統計分布からシャノンエントロピーを計算すると、DNAメチル化データおよび医療記録データの双方において、エントロピーが年齢とともにほぼ線形に増加することが示された（図4）。これは、老化の進行が生体状態のエントロピー増加、すなわち情報の喪失として記述できることを意味する。
　この観点から見ると、老化は生体状態が時間とともに不可逆的に変化していく過程であり、その進行は機能単位の状態遷移の累積量、すなわち tBA によって特徴づけられる。エイジングクロックは、この tBA を異なる生物医学データから推定する指標として理解することができる。
　また、tBA の増加は生体の調節ネットワークにおける平均場として作用し、生理機能のドリフトや回復力の低下を引き起こす。さらに、この累積的変化は慢性疾患の発症率や死亡率の指数関数的増加を説明する要因となる可能性がある。

6. ディスカッション

老化のエントロピー的性質と若返りの可能性と限界

　本研究では、DNAメチル化データおよび電子医療記録データの解析を通じて、老化に伴う変化が多数の確率的状態遷移の累積として説明できることを示した。この結果は、老化が単一の分子機構によって決定される過程ではなく、複雑な調節ネットワークにおける微視的状態変化の統計的な蓄積として理解できることを示唆している。
　提案されたモデルでは、生体の老化は機能単位の状態遷移の累積によって特徴づけられ、その累積量は熱力学的生物学的年齢（tBA）として表される。tBAは老化に伴うエントロピー生成量に対応し、生体状態の情報喪失を表す指標とみなすことができる。このような視点から見ると、老化は生体状態の不可逆的ドリフトとして理解される。
　この累積的変化は生体の回復力（resilience）の低下を引き起こす。モデルによれば、生体の調節ネットワークには回復速度によって特徴づけられる動的安定性が存在するが、状態遷移が蓄積するにつれて回復速度は徐々に低下する。理論的には、ある臨界点において回復速度がゼロに近づき、生体は恒常性を維持できなくなる。この臨界点はおよそ100〜150歳付近に位置する可能性が示唆されている。
　この観点から、人間の老化には強いエントロピー的成分が存在すると考えられる。エントロピー増加は本質的に不可逆的な過程であるため、老化を完全に巻き戻すことは理論的に困難である可能性がある。したがって、特定の老化関連経路や生理機能を一時的に改善することは可能であっても、生体全体の老化状態を完全に若返らせることは容易ではないと考えられる。
　一方で、このモデルは老化速度を低下させる可能性についても示唆している。機能単位の状態遷移は調節ネットワークにおけるノイズの強さ、すなわち有効温度に依存すると考えられる。この有効温度を低下させることができれば、状態遷移の速度が減少し、老化の進行を大きく遅らせる可能性がある。
　例えば、ハダカデバネズミのような長寿哺乳類では、DNA修復能力やストレス耐性が高く、調節ネットワークにおけるノイズが低い可能性が指摘されている。このような生理学的特性は、老化速度を低下させる要因として機能している可能性がある。
　以上のように、本研究は老化をエントロピー生成過程として理解する理論的枠組みを提示した。この視点は、生体の老化ダイナミクスを理解する新たな手がかりを提供するとともに、老化速度を制御する生物学的要因の探索に重要な示唆を与えるものである。