Somatic mutations and single-nucleus transcriptomics reveal uniquely human properties of neuronal aging

問題提起

ヒトの脳の老化は、自然な加齢に伴う神経変性疾患をほとんど示さない哺乳類を用いてモデル化されることが多い。体細胞変異（単一塩基バリアント［SNV］および短鎖挿入・欠失［Indel］）はすべての細胞に蓄積するが、ゲノムの安定性、変異蓄積率、および種間の寿命の関係については、依然として理解が進んでいない。先行研究における腸管クリプト（隠窩）では、長寿の種は体細胞変異をより緩やかな速度で蓄積するため、哺乳類間で終末期の変明負担が比較的一定になることが示唆されていた。しかし、この進化的な「天井」が、細胞の補充が不可能なポストミトーシス（分裂後）状態のニューロンにも適用されるのかどうかは不明である。ニューロンはDNA損傷の不可逆的な蓄積に直面している。ヒトの神経変性疾患に対する脆弱性と、近縁種と比較して数十年に及ぶ寿命の可能性を考慮すると、種を超えたニューロンのゲノム変化を体系的に比較することが極めて重要である。

手法

著者らは、多種解析によるシングルセルアプローチを用いて、背外側前頭前野（PFC）における体細胞変異の景観とトランスクリプトームの変化を特徴付けた。

1. SCAN2-Zooの開発: 非ヒト種における体細胞変異を解析するために、著者らは既存のヒト特異的なアルゴリズムであるSCAN2を改良し、SCAN2-Zooを開発した。この新しいパイプラインは、非ヒトのデータセットでは希少なフェーズ化されたヘテロ接合型SNPの必要性を、最大変異アレル頻度（maxVAF）アプローチに置き換えることで、多様な二倍体種において高い忠実度で体細胞変異のコーリングを維持している。
2. サンプル収集とソーティング: 本研究では、ヒト、チンパンジー、アカゲザル、マーモセット、マウス、フェレットの6種類の哺乳類から得られた死後脳組織を利用した。細胞型の不均一性を最小限に抑えるため、著者らはNeuN、SOX6、SOX10に対する抗体を用いた蛍光活性化核ソーティング（FANS）戦略を開発し、純粋なグルタミン酸作動性ニューロンおよびMGE-GABA作動性ニューロン（パルブアルブミン/ソマトスタチン発現）の集団を単離した。
3. シングルセル全ゲノムシーケンシング（scWGS）: ソーティングされた単一核は、プライマリーテンプレート増幅（PTA）による全ゲノム増幅を経て、25–30Xのカバー率でシーケンシングされた。変異はSCAN2-Zooを用いてコールされた。
4. META-CSによる検証: 知見を検証するため、ヒトとアカゲザルのサブセットに対し、真の体細胞変異とアーティファクトを区別するためにストランドコンセンサスに依存する技術である**マルチプレックス・エンドタグ増幅補完鎖（META-CS）**を用いた。
5. 変異シグネチャー解析: 著者らは、3塩基配列スペクトラムに対して非負値行列因子分解（NMF）を用いたde novo変異シグネチャー抽出を行い、COSMICデータベースと比較した。また、神経変性疾患（アルツハイマー病、CTE、FTD、ALS）における既存の単一ニューロンデータも解析した。
6. シングル核トランスクリプトミクス（snRNA-seq）: 機能的影響を評価するため、成体および高齢のヒト、チンパンジー、アカゲザルのソーティングされた核に対してsnRNA-seqを実施した。環境RNAのアーティファクトを軽減するため、異なる種の細胞を同一の反応内で処理し、種のアサインメントは一意分子識別子（UMI）のミスマッチに基づきバイオインフォマティクス的に行った。差分的遺伝子発現（DGE）解析は擬似バルクカウントに対して行われた。

主な結果

1. 年間変異率の保存性と、終末期負担の乖離
腸管クリプトで観察されたパターンとは対照的に、皮質ニューロンにおける体細胞SNVおよびIndelの年間蓄積率は、6種すべてにおいて驚くほど類似していた（~18–24 SNVs/年）。

• ヒト: ~18 SNVs/年。
• マウス/フェレット: ~18–24 SNVs/年。
• 帰結: ヒトは著しく長寿であるため、ヒトのニューロンにおける終末期の変異負担は、マウスよりも遥かに高く（12倍以上のSNVと数百のIndel）、チンパンジーやアカゲザルよりも有意に高い。

2. 種特異的な変異シグネチャー
年間率は同様であるが、基礎となる変異メカニズムは異なる：

• SBS-C (ヒト特異的): C>T置換（特にCpGおよび非CpG部位）が支配的なシグネチャーであり、ヒトでは非ヒト霊長類（ほぼゼロ）と比較して有意に速い速度で蓄積する（総変異の~15%）。このシグネチャーは、塩基除去修復（BER）不全および酸化ストレスに関連するCOSMIC SBS30に類似している。
• SBS-A (非ヒト霊長類特異的): T>AおよびT>C置換が支配的なシグネチャーであり、アカゲザルやその他の非ヒト霊長類で高い割合で蓄積するが、ヒトでは低い。
• 疾患との相関: ヒストリカルなコントロールと比較して、アルツハイマー病、CTE、FTD、およびALSの患者のニューロンにおいて、ヒト特異的なSBS-Cシグネチャーが有意に上昇している。

3. ヒトの老化におけるトランスクリプトームの調節不全
シングル核トランスクリプトミクスにより、ヒトの老化ニューロンは、チンパンジーやアカゲザルでは同程度には観察されない広範な転写変化を起こすことが明らかになった。

• 遺伝子ダウンレギュレーション: 加齢したヒトのニューロンは、数百の遺伝子の広範なダウンレギュレーションを示すが、加齢した非ヒト霊長類では、ダウンレギュレーションされる遺伝子は著しく少ない。
• 機能経路: ヒトでダウンレギュレートされた遺伝子は、プロテオスタシス（シャペロン、リボソームタンパク質、プロテオソーム成分）およびミトコンドリア機能（呼吸複合体）に高度に濃縮されている。
• 特定のマーカー: 長寿種であるホッキョククジラに関連する遺伝子CIRBPは、高齢のヒトにおいて特異的にダウンレギュレート（>2倍）されていた。

意義と主張

本論文は、ヒトのニューロンの老化が、分子レベルにおいて、その近縁種や標準的な動物モデルと比較して、一意に深刻であることを主張している。

• 進化的制約: 著者らは、ニューロンにおける体細胞変異蓄積の年間率には（シナプス可塑性の維持や即時的な毒性の回避のために）進化的制約があるかもしれないが、終末期の負担にはそのような制約はないと提案している。これにより、長寿のヒトにおいて潜在的に有害となるレベルまで負担が蓄積する一方で、短命な種ではそのレベルに達しないという現象が生じる。
• 動物モデルの限界: これらの知見は、マウスやアカゲザルといった動物モデルが、「一意のヒト」の晩期状態（高い体細胞変異負担、特定の変異シグネチャー［SBS-C］、および広範なプロテオスタシス／ミトコンドリア崩壊を特徴とする状態）を再現できていないことを示唆している。
• 神経変性との関連: 酸化ストレス修復経路に関連するヒト特異的なSBS-Cシグネチャーの蓄積と、神経変性疾患におけるその上昇は、ヒトの老化に内在するゲノム不安定性が神経変性を誘発する可能性があることを示唆している。
• 衰退のメカニズム: 広範な遺伝子ダウンレギュレーションは、コーディング領域における高い変異負担のみによって駆動されているわけではない（ダウンレギュレートされた遺伝子に過剰な負担は見られなかったため）、これは、老化するゲノムによって引き起こされる他の形態の遺伝的不安定性やストレス経路が転写の低下を駆動していることを示唆している。

要約すると、本研究は、ヒトのニューロンは他の哺乳類と同様の速度で体細胞変異を蓄積するが、より長いタイムフレームを経て、結果として一意に高い変異負荷と、酸化ストレスに関連する特定の変異シグネチャーをもたらすことを確立した。このゲノム不安定性は、プロテオスタシスとミトコンドリア機能の深刻な、種特異的な崩壊と一致しており、現在の非ヒトモデルでは捉えきれない老化の分子的景観を定義している。