From Genes to Brains: Molecular Evolution and Mammalian Brain Cellular Diversity

技術要約：遺伝子から脳へ：分子進化と哺乳類の脳細胞多様性

問題提起
哺乳類の脳は、サイズ、細胞組成、および組織化において広範な多様性を示すが、その表現型の多様化を駆動する分子進化的な変化については、まだ完全には理解されていない。これまでの研究では、特に霊長類において、遺伝子発現やタンパク質進化における系統特異的な変化が記録されているが、ゲノムワイドな評価によって、分子的な分岐が哺乳類における脳細胞の表現型の定量的差異を系統的に追跡しているかどうかを判断する研究は不足している。このような種間解析における中心的な課題は、共通の進化史を考慮することであり、これにより種は統計的に独立した観察値ではなくなる。本研究では、系統比較フレームワークを利用して、特定の遺伝子と表現型の関連を共通の祖先由来のアーティファクト（偽の相関）から区別することにより、タンパク質配列の進化が哺乳類における脳細胞組織化の定量的な差異と共変するかどうかを調査する。

方法論
本研究では、18種の哺乳類（Homo sapiens、Mus musculus、Loxodonta africanaを含む）における比較神経解剖学的データと、オーソログなタンパク質配列を統合した。

1. データの統合と形質の選択:

   • 表現型: 著者らは、Herculano-Houzelら [2015] による29の脳および身体関連形態形質（例：ニューロン数、細胞密度、構造質量）を用いた。
   • 次元削減: 形質間に強い相関があるため、徹底的なサブセット選択分析が行われた。標準化された形質の分散の88.12%を説明する5つの代表的な形質が選択された：皮質ニューロン数（cx_N）、小脳の単位質量あたりのニューロン数（cb_NpM）、体質量（body_M）、皮質質量の割合（cx_percM）、および小脳ニューロンの割合（cb_percN）。
   • ゲノムデータ: 18種すべての10,395個のヒト遺伝子に対するオーソログタンパク質配列をOrthoDB v12.2から取得した。

2. 進化速度の解析:

   • アライメントと樹状図の構築: アミノ酸配列はMUSCLE v5.3を用いてアライメントされ、遺伝子特異的な枝の長さは、VertLife/MamPhyに由来する共通の18種哺乳類トポロジーに制約された上で、IQ-TREE v3.0.1を用いて推定された。
   • RERconverge: 著者らは、遺伝子特異的な相対進化速度（RER）と5つの代表的形質の進化的な変化との間の関連をテストするためにRERconvergeを用いた。RERは、平方根変換された分子経路長と加重回帰を用いて算出された。
   • 統計検定: 相関はスペアマンの順位相関係数を用いてテストされた。解析は、対数変換された表現型（主要解析）および未変換の表現型（感度分析）の両方に対して行われた。有意性はBenjamini-Hochberg補正（FDR < 0.05）によって決定された。

3. ダウンストリームの特性評価:

   • 発現プロファイリング: 両方の解析で有意であった遺伝子は、13種類の成体ヒト脳組織におけるGTEx v11データを用いて特性評価された。
   • 高発現比較: 上位500個の脳高発現遺伝子のRER相関プロファイルを、残りの遺伝子と比較した。これにはKolmogorov–Smirnov検定が用いられた。
   • 機能的濃縮: オーソロジー・フィルタリングによって導入される機能的なバイアスを評価するために、10,395個の保持された遺伝子に対してKEGGパスウェイ濃縮解析を行った。

主な結果

• 遺伝子と形質の関連: 主要な解析（対数変換された表現型）では、23個のユニークな遺伝子を含む25個の有意な遺伝子–形質関連が特定された。未変換の解析では、28個の遺伝子を含む29個の関連が特定された。決定的なことに、19個の遺伝子が両方の解析で回収されており、これは、表現型の変換に依存しない堅牢なサブセットであることを示している。
• 形質の不均一性: 有意な関連は不均一に分布していた。最も多くの関連が見られたのは、小脳ニューロンの割合（cb_percN）および皮質質量の割合（cx_percM）であり、体質量や小脳の単位質量あたりのニューロン数については、より少ない数の関連が検出された。
• 発現パターン: 19個の共有候補遺伝子は、成体ヒト脳組織において不均一な発現を示した。一部の遺伝子（例：ABHD12, NKTR）は広範かつ高度に発現していたが、他の遺伝子（例：IL22RA2, DYTN, CEP55）は成体脳組織における発現が低く、進化速度の関連が、一様に高い成体脳発現を持つ遺伝子に限定されないことを示唆している。
• 高発現遺伝子のプロファイル: 高度に脳で発現している遺伝子（上位500個）は、残りのゲノムとは異なるRER関連プロファイルを示した。3つの形質（cx_N, cx_percM, cb_percN）について、相関係数の分布は高発現グループにおいて有意に正の値へとシフトしていた。
• 機能的濃縮: 保持された遺伝子セット（10,395個）は、代謝経路（例：クエン酸回路、解糖系）、細胞恒常性、糖鎖バイオシンセシス、軸索ガイダンス、および神経活性リガンド–受容体相互作用に濃縮されていた。対照的に、除外された遺伝子は、嗅覚伝達、免疫関連パスウェイ、およびリボソーム/スプライソソームの機能に濃縮されており、これはオーソログの厳格な要件を反映していると考えられる。

意義と主張
著者らは、本研究がタンパク質コード遺伝子の相対的進化速度を哺乳類の脳細胞構造の定量的変異に結びつける比較フレームワークを提供するものであると主張している。主な貢献は以下の通りである：

1. 特定の相関因子の特定: 本研究は、進化速度が量的脳表現型と共変する特定の遺伝子セットを特定しており、単一遺伝子または系統特異的なケーススタディを超えた、ゲノムワイドな評価を実現している。
2. 変換に対する堅牢性: 対数変換および未変換の両方の解析において相当数の共通遺伝子が回収されたことは、多くの検出された関連が表現型の表現スケールに対して堅牢であることを示唆している。ただし、一部（特に体質量に関して）はスケールに敏感である。
3. 進化シグナルの文脈化: 進化速度の有意な関連が、成体における高い脳発現を持つ遺伝子と低い発現を持つ遺伝子の両方に存在することを実証することで、本研究は、脳の多様化の分子的相関が、成体ヒト脳における構成的な高発現に一律に限定されないことを示唆している。
4. 形質依存的な進化: 関連の不均一な分布は、哺乳類の脳の多様化が、脳全体のサイズを規定する単一の分子軸ではなく、形質依存的な分子相関を伴うという見解を支持している。

著者らは、これらの関連は、直接的な因果的遺伝子–表現型の関係を示す証拠ではなく、分子的変化と表現型的変化の間の「進化的な結合」の候補として解釈されるべきであると明示している。また、種セットが比較的少ないこと（18種）、オーソログ・フィルタリングによる系統特異的な遺伝子や調節進化の除外、および広範な哺乳類の発現コンテキストのプロキシとしての成体ヒト発現データの依存性といった限界についても認めている。本研究は、哺乳類の脳の多様性の分子基盤を完全に解明するためには、これらのフレームワークをより大きな種セットへと拡張し、調節配列の進化を取り入れることが必要であると結論付けている。