〜 おへその乱読1000本ノック #0043 ～ 細胞の一生を追う：多層オミクスで捉えるヒト細胞老化の縦断データ

　第43回「おへその乱読1000本ノック」へようこそ。
　老化は、ゲノムやエピゲノム、転写、タンパク質、さらには細胞機能に至るまで、多層的な変化が積み重なるプロセスです。近年、オミクス解析技術の発展によりこれらを包括的に捉えることが可能になりましたが、時間軸に沿って統合的に理解することは依然として大きな課題です。ヒト個体でこれを追跡することは、技術的・倫理的に容易ではありません。
　こうした制約を補う手法として、培養細胞を用いた老化研究が用いられています。ヒト初代線維芽細胞は分裂を繰り返すことで老化様変化を示し、その進行は生体内より大幅に速く、数百日の培養がヒトの数十年に相当する時間スケールとして扱えるとされています。この系により、老化のダイナミクスを現実的な時間で追跡することが可能になります。
　今回取り上げるSturm et al (2022)は、この複製寿命モデルにおいて多層オミクスデータを縦断的に取得した研究です。DNAメチル化、遺伝子発現、ミトコンドリア機能などを同一細胞系で時系列的に測定し、さらに代謝経路への介入を組み合わせることで、老化の分子基盤に迫る設計となっています。
　老化の進行を規定する要因はどこにあるのか。そして、細胞内エネルギー代謝はその時間軸にどのように関与するのか。本研究を手がかりに、老化を動的なシステムとして捉える視点を見ていきます。

1. 背景と目的

　老化は、ゲノム、エピゲノム、トランスクリプトーム、プロテオーム、さらには臓器・個体レベルに至るまで、多層的な分子変化の蓄積として特徴づけられる現象である。近年、老化の分子基盤に関する研究の進展により、老化のホールマークと呼ばれる共通原理が提唱され、細胞および個体の機能低下や疾患リスク増大の根底にある原因として整理されてきた。
　一方で、オミクス解析技術と機械学習手法の発展により、DNAメチル化に基づくエピジェネティック時計など、暦年齢や死亡リスクを高精度に推定する分子指標が確立されている。しかし、これらの指標は主として記述的な性質を持ち、その背後にある老化の分子機構、とりわけ老化のホールマークとの因果的なつながりは十分に解明されていない。実際、本研究においても、細胞増殖に関連する遺伝子の低下や老化関連指標の変化など、細胞老化に特徴的な分子変化が観察されている（Fig.5）。

　老化の分子機構を理解するためには、多様な分子指標を長期間にわたり高頻度で縦断的に測定する研究が求められるが、ヒト個体においてこれを実現することは困難である。この課題に対する代替アプローチとして、本研究ではヒト初代線維芽細胞の複製寿命モデルを用いた細胞レベルの縦断解析を採用している（Fig.1A）。この系では、培養下での細胞老化が生体内の老化関連変化の一部を再現しつつ、老化プロセスが in vivo よりも加速して進行するため、数百日スケールでヒトの数十年に相当する変化を追跡できる。
　さらに、ミトコンドリアはエネルギー産生に加え、代謝シグナルを介して核およびエピゲノムに影響を及ぼすことが知られている（Fig.1B）。このような細胞内シグナル伝達を背景として、DNAメチル化などのエピゲノム指標はミトコンドリア機能と関連する形で変化しうる（Fig.1C）。
　本研究では、ヒト初代線維芽細胞の複製寿命に沿った多層オミクスの縦断解析に加え、ミトコンドリア代謝およびストレス応答経路への介入を組み合わせて、老化に伴う分子変化と代謝機構との関係を検討する（Fig.1D）。

2．研究戦略

　老化の分子機構を明らかにするためには、長期間にわたる縦断的な観察が理想的であるが、ヒト個体を対象とした研究には時間的・技術的制約が伴う。この課題に対する代替アプローチとして、本研究ではヒト初代線維芽細胞の複製寿命モデルを採用した。
　初代線維芽細胞は、培養下で分裂を繰り返す過程において増殖停止に至るまでの一連の変化を示し、細胞レベルでの老化様現象を再現するモデルとして用いられてきた。本研究では、この複製寿命にわたる変化を追跡することで、同一個体由来細胞における老化過程を縦断的に評価した。
　さらに、この培養系では老化関連の分子変化が加速した時間スケールで進行することが報告されている。本研究では、約200〜300日の培養期間がヒトにおける数十年の時間経過に相当する系として扱い、老化に伴う分子変化を現実的な時間枠内で観測する戦略を採用した。

3．実験デザイン

　本研究では、ヒト初代線維芽細胞を用いた複製寿命モデルにおいて、細胞老化の進行を縦断的に追跡する実験系を構築した。細胞は健常ドナー6名およびミトコンドリア疾患（SURF1変異）を有するドナー3名から取得し、年齢、性別、組織由来などの情報をTable 1に示した。
　各細胞系は培養下で継代し、約270日間にわたり維持した。細胞は一定間隔で継代およびサンプリングを行い、分子解析および細胞機能評価のための試料を継時的に取得した。継代はおおよそ5〜7日ごとに実施し、細胞が増殖停止に至るまで追跡した。
　細胞寿命の終点は、一定期間における増殖率の低下、すなわちpopulation doublingの停止により定義した。また、実験は複数のフェーズに分けて実施し、異なる条件および実験者における再現性を検証した。

4．multi-omicsにおける測定項目

　本研究では、細胞老化の進行を多面的に評価するため、複数の生物学的指標を縦断的に測定した。これらの測定項目は、細胞レベルの指標から分子レベルのオミクスデータに至るまで多岐にわたり、その概要をTable 2およびFig. 2に示した。
　細胞レベルでは、細胞数、細胞サイズ、細胞死の割合などの細胞動態に加え、顕微鏡画像による形態評価を行った。また、細胞のエネルギー代謝状態を評価するため、酸素消費速度（OCR）および細胞外酸性化速度（ECAR）を指標としたバイオエネルギー測定を実施した。
　分子レベルでは、RNAシーケンスによりトランスクリプトーム解析を行い、遺伝子発現の変化を網羅的に評価した。エピゲノムについては、Illumina EPICアレイを用いてDNAメチル化を測定し、多数のCpGサイトにおけるメチル化状態を定量した。さらに、これらのデータから複数のエピジェネティック時計を算出した。
　加えて、テロメア長の測定、ミトコンドリアDNAコピー数および変異解析、全ゲノムシーケンスを実施した。さらに、培養上清中のサイトカインおよびcell-free DNAを定量し、分泌因子の変化を評価した。
　これらの測定は、各細胞系および条件において複数の時点で繰り返し取得し、高い時間分解能をもつ縦断データとして統合した。

5．介入条件（実験操作）

　本研究では、細胞老化に関連する分子過程を検討するため、複数の代謝経路およびシグナル伝達経路に対する介入を組み合わせた実験系を構築した。各細胞系および処理条件の対応関係をFig. 3に示し、実験フェーズの構成をTable 3に示した。
　介入は主として、ミトコンドリア機能および細胞代謝に関連する経路を対象として設計した。具体的には、グルココルチコイド受容体作動薬であるデキサメタゾン（DEX）を用いて慢性的なストレス応答を誘導した。また、ATP合成酵素を阻害するoligomycinにより酸化的リン酸化（OxPhos）を抑制し、ミトコンドリア由来ATP産生の低下を引き起こした。さらに、ミトコンドリアへの基質供給を阻害するmitoNUITs（pyruvate、脂肪酸、グルタミンの取り込み阻害）を用い、TCA回路への炭素供給を制限した。
　加えて、酸素濃度を低下させる低酸素条件（3% O₂）により細胞のエネルギー代謝状態を変化させたほか、galactose培地や2-deoxy-D-glucose（2-DG）を用いて解糖系の制御を行った。また、β-hydroxybutyrateの添加によりケトン体代謝の影響を評価した。
　酸化ストレスに関連する介入として、ミトコンドリア特異的抗酸化物質であるMitoQおよびグルタチオン前駆体であるN-acetyl-L-cysteine（NAC）を用いた。また、DNAメチル化状態への直接的な介入として、5-azacytidineを用いた脱メチル化処理を実施した。
　これらの介入は単独または組み合わせて適用され、各処理条件における細胞老化の進行および分子指標の変化が縦断的に評価された。さらに、これらの実験は複数のフェーズ（Phase I–V）に分けて実施され、異なる条件および実験者間での再現性が検証された。
　また、各測定項目および条件におけるサンプル数および時間分布はFig. 4に示されており、本研究のデータセットが多様な条件および時間軸を包含することが示されている。

6．データの特徴と利用可能性

　本研究で構築したデータセットは、細胞老化の進行を多層的に捉えることを目的として、複数の生物学的指標を縦断的に統合したものである。細胞動態、バイオエネルギー、トランスクリプトーム、DNAメチル化、ゲノム情報、分泌因子など、13種類の主要なアウトカムを取得した。
　各測定項目は、細胞系および処理条件に応じて異なる頻度で取得されており、細胞寿命に沿った複数時点のデータが蓄積されている。これにより、時間分解能の異なる情報が統合された縦断的な多層オミクスデータセットが構成されている。
　トランスクリプトーム解析では数百サンプルにわたり遺伝子発現データが取得され、DNAメチル化についても多数のCpGサイトを対象とした大規模データが得られている。さらに、バイオエネルギー指標や分泌因子などの機能的指標も並行して測定されており、分子変化と細胞機能の関係を統合的に解析可能な構成となっている。
　本データセットは公開されており、専用のウェブツールを通じて可視化およびダウンロードが可能である（Fig.6）。ユーザーはドナーや処理条件を選択することで、細胞寿命に沿った時系列データを閲覧でき、各種データはCSV形式で取得できる。また、RNA-seqおよびDNAメチル化データは公共データベースを通じて提供されている。全ゲノムシーケンスデータについては共有に制限があるものの、変異情報などの一部は可視化ツールを通じて利用可能であり、追加データは著者への申請により取得できる。
　以上のように、本研究で構築したデータセットは、多様な分子層および機能指標を統合した高時間分解能の縦断データとして提供されており、細胞老化の分子基盤を多角的に解析するための基盤となる。

7．技術的妥当性（Technical Validation）

　本研究では、取得したデータの信頼性および再現性を担保するため、細胞培養および分子測定の各段階において技術的検証を行った。継代間隔は全期間にわたり記録されており、細胞寿命の進行に伴い、増殖能の低下に応じて継代間隔のばらつきが増加する傾向が確認された（Fig.7A）。
　また、健常ドナー由来の細胞系については、異なる実験フェーズおよび実験者により複数回の独立実験を実施し、増殖曲線の再現性を検証した。その結果、独立した実験間で一貫した成長パターンが確認された（Fig.7B）。
　分子測定に関しては、各種測定系において技術的ばらつきの評価を行った。テロメア長測定では複数の技術的リプリケートを用いて変動係数を算出し、測定精度を確認した（Fig.7C）。また、ミトコンドリアDNAコピー数およびcell-free DNAの定量についても、複数プレートにわたる測定により再現性を評価した（Fig.7D）。
　さらに、DNAメチル化およびRNAシーケンスデータについては、品質管理および正規化処理を実施し、プローブ特性やプレート間差異に由来するバッチ効果を補正した。これにより、データ間の比較可能性および解析の信頼性を確保した。
　以上のように、本研究では細胞培養から分子解析に至るまでの各段階において体系的な技術的検証を行い、データの品質および再現性を担保した。これにより、本データセットは細胞老化の分子基盤を解析するための信頼性の高い基盤として提供される。

8．最後に

　本研究は、ヒト初代線維芽細胞の複製寿命をモデルとして、多層オミクスデータを高時間分解能で縦断的に取得したデータセットを提示したものである。さらに、ミトコンドリア代謝やストレス応答経路への介入を組み合わせることで、細胞老化に伴う分子変化を多面的に評価可能な実験系を構築した。本データセットは、老化に伴う分子変化を統合的に解析するための基盤として活用されることが期待される。