- 著者: Elisa de Launoit, Sophie Skriabine, Etienne Doumazane, Mathilde Bizou, Remi Pomme, Gabriele Lienhard, Domitille Rajot, Alina Mikhailenko, Alba Vieites-Prado, Mathilde Lannes, Natalia Fabiano, Louise Mathe, Gael Cagnone, Laura Mouton, Mathieu D. Santin, Severine Candelier, Patricia Gaspar, Charly V. Rousseau, Alexandre Dubrac, Nicolas Renier
- Corresponding author: Alexandre Dubrac (CHU Sainte-Justine, Montreal), Nicolas Renier (ICM Paris Brain Institute)
- 雑誌: Cell
- 発行年: 2026
- Epub日: N/A
- Article種別: Original Article
- PMID: 41985458
背景
出生時の哺乳類における脳血管ネットワークは極めて未成熟であり、アストロサイトやペリサイトによる被覆も不十分である。生後の脳発達期 (postnatal period) においては、脳容積の急速な拡大と複雑な神経回路の成熟が同時に進行する。これに伴い、増大する代謝需要を支えるための活発な血管新生と再構築が必須となるが、全脳規模での定量的かつ時空間的な記述はこれまで技術的限界から困難であった。先行研究である Girouard et al. (2006) や Iadecola (2017) では、成体脳における神経血管カップリングの重要性が示されているものの、発達期における全脳レベルの動態は不明であった。また、Paredes et al. (2018) などの既報においても、血管発達を制御する分子シグナルと実際の血管構造変化との時空間的な相関関係は、皮質などの局所的な領域や網膜モデルでの解析に限定されており、脳全体を俯瞰した統合的な理解には至っていなかった。特に、どの分子シグナルがどの発達段階において、どの脳領域の血管成長を駆動しているのかという問いに対して、空間トランスクリプトミクスと3次元血管構造データを統合して解析できる発達期ブレインアトラスが存在しないことが大きな課題であった。さらに、神経活動依存的な血管形成 (activity-dependent angiogenesis) の概念は提唱されていたものの、感覚経験の開始 (開眼やひげによる探索行動の開始など) と血管網の発達が全脳規模でどのように時間的・空間的に一致しているのかを実証したデータは不足していた。脳血管発達における性差や、グリア細胞の成熟に伴う時空間的な連関についても定量的な評価が不十分であり、これらの解明を阻む大きな知識ギャップが存在していた。
目的
本研究の目的は、組織透明化技術、光シート顕微鏡、3次元血管グラフ再構築、および空間トランスクリプトミクス技術を統合した、マウス生後発達脳アトラス「LAMBADA (Light-sheet Atlas of Mouse Brain Angiogenesis During postnatal Development with mApple)」を構築することである。これにより、生後マウス脳全体における脈管形成 of 時空間的ダイナミクスを定量的に明らかにし、血管発達を駆動する分子基盤と神経回路成熟との相互作用を全脳規模で解明することを目指した。特に、生後発達期の脳組織収縮を補正するための超高磁場 MRI (magnetic resonance imaging) データの統合や、AI (artificial intelligence) 技術を用いた血管グラフの精度向上を達成し、従来の2次元的あるいは局所的な解析限界を打破して、脳血管発達の全体像を3次元で精密に記述することを目的とした。
結果
3相にわたる生後脳血管発達プログラムの定量的同定: 全脳の血管長密度、分岐点密度、断端密度、および血管の tortuosity (屈曲度) の時系列推移を定量解析した結果、生後の血管発達が明確に異なる3つの相 (Phase) を経て進行することが明らかになった (Figure 2)。Phase 1 (P3-P7、等尺性拡大相) では、脳容積が急速に拡大する一方で、血管長密度および分岐点密度はほぼ一定に維持されていた。この時期、断端血管 (sprouting を行う tip cell を有する未成熟な血管) の割合が極めて高く、P5 および P7 においては全断端の 80% 以上に tip cell が観察された (Figure S3)。これは、血管網が脳の等尺性拡大 (isometric expansion) に追従して活発に新生している状態を反映している。Phase 2 (P7-P21、領域特化相) に入ると、血管長密度が急激に上昇し、脳領域ごとの特異的な発達軌跡が顕著となった。特に Phase 2b (P12-P21) において、小脳、中脳、視床、後脳などの高密度化領域と、線条体、海馬、大脳皮質腹側部などの低密度化領域への分岐が明確に観察された (Figure 3)。例えば、聴覚皮質においては P12 から P14 にかけて、皮質第IV層特異的に水平方向の血管方向性 (anisotropy) が急速に確立されることが示された (Figure 4)。Phase 3 (P21-成体、精緻化相) では、血管長密度が安定化する一方で、コラーゲンIV (collagen IV) 陽性の安定化した末梢血管の割合が P21 以降に有意に増加し、血管網の成熟と安定化が進行した (Figure S3)。また、SM22 陽性の動脈性血管の比率は成体値 (約100本の血管枝に対して1本の割合) に収束した。
血管発達を駆動する分子機構の動的相転換: 空間トランスクリプトミクスデータと血管構造データの相関解析により、血管新生を駆動する主要な分子シグナルが発達段階に応じて劇的にシフトすることが実証された (Figure 6)。Phase 1 においては、低酸素応答シグナルである Vegfa および Wnt7b の発現パターンが、全脳および領域別の断端血管密度と強い正の相関を示した (Spearman r>0.5、n=12 planes)。これは、初期の等尺性拡大が組織の相対的低酸素状態に応答した古典的な VEGFA (vascular endothelial growth factor A) 経路によって主導されていることを示している。特に、P3 から P21 の移行期に Vegfa 陽性細胞 (n=15000 cells) の 1.5-fold decrease (log2FC -0.58) を伴う発現低下が確認された。しかし、Phase 2 に移行すると Vegfa と血管新生との相関は消失し、代わって Sema3e (前脳で正相関) や Slit2 (全脳で一貫して負相関、Spearman r←0.3) などの軸索ガイダンス分子や、Apln (apelin) および Wnt9a などのリガンドが血管リモデリングの主要な相関因子として浮上した (Figure 7)。さらに、Allen Brain ADMBA (Allen Brain Developing Mouse Brain Atlas) の公開データを用いた P4 と P14 の全脳相関解析 (n=220 genes) からも、P4 特異的遺伝子群が「HIF (hypoxia-inducible factor) による低酸素応答」や「血管新生」に富化しているのに対し、P14 特異的遺伝子群は「シナプス成熟」や「神経活動依存的プロセス」に強く富化していることが PANTHER (protein analysis through evolutionary relationships) パスウェイ解析により示された (Figure 7)。
神経回路成熟およびグリア細胞との時空間的同調: cNMF 解析により抽出された30の遺伝子プログラムの解析から、神経細胞の成熟段階と血管網の再構築が極めて精緻に同調していることが明らかになった (Figure 7)。Cortex において、Phase 1 でピークを示すプログラムは Sox4、Sox11、Dcx などの幼若ニューロンおよび軸索伸長マーカーに富化していた。Phase 2 でピークを示すプログラムは Satb1 や Cux1 などの皮質ニューロンのアイデンティティ確立およびシナプス成熟 (Lin7c、Nptx2) に関連していた。Phase 3 でピークを示すプログラムは、Fos、Egr1 などの IEG (immediate early gene) や Pvalb などの活動依存的なシナプス伝達機能に富化しており、これらは後期血管リモデリングパラメーターと空間的に有意に相関していた (p<0.05)。また、非神経細胞プログラムの解析から、アストロサイトの成熟に伴い発現するアンジオテンシノゲン (Agt) や Sema4a が血管密度と強い負の相関を示し (Figure 7)、アストロサイトの成熟が血管新生に対して「制動 (brake)」をかけ、ネットワークを安定化させる分子スイッチとして機能していることが示唆された。
血管発達における性的二型の出現: 雌雄別の全脳血管密度解析 (各群 n=5 mice) を実施した結果、P5 および P12 の段階では全脳のどの領域においても有意な性差は認められなかった。しかし、思春期前後に相当する P30 においては、雄のマウスで脳梁、視床下部の tuberal 核および supraoptic 核の血管密度が雌に比べて有意に高かった。逆に、雌のマウスでは diagonal band nucleus (dbn) および bed nucleus of the stria terminalis (bnst) において、雄よりも有意に高い血管密度が観察された (Figure S4)。この結果は、生後初期の「ミニ・パブティ (mini-puberty)」期ではなく、思春期における性ステロイドホルモンなどの分泌変化が、脳の特定の性的二型領域における血管再構築を駆動していることを示唆している。
技術的革新による血管グラフの精度向上: 本研究で開発された Video Swin Transformer ベースの3次元画像分類 AI 検出器の導入により、従来の TubeMap パイプラインで課題となっていた、血管の急峻な屈曲や輝度低下に起因する人工的な断端 (artefactual endings) を成体脳の再構築データにおいて 91% 削減することに成功した (p<0.001, Figure S3)。これにより、真の sprouting 血管 (tip cell 陽性) と技術的ノイズを極めて高い精度で分離することが可能となり、生後発達期における微細な血管新生ダイナミクスの正確な定量化が初めて実現した。この分類器の性能は、ROC AUC (receiver operating characteristic area under the curve) で 0.96 という極めて高い精度を達成した。
考察/結論
先行研究との違い: 従来の生後脳血管発達研究は、Coelho-Santos et al. (2021) などの報告に代表されるように、大脳皮質表面の局所的な観察や網膜モデルにおける 2次元的な解析に依存しており、皮質下核や脳幹、小脳を含む脳全体の3次元的な発達動態を同時に定量評価することは不可能であった。本研究は、新しく構築した LAMBADA アトラスを用いることで、全脳 62 領域における生後血管発達の定量的軌跡を初めて同一の空間基準で記述した点において、これまでの局所的研究と決定的に異なっている。
新規性: 本研究は、生後の脳血管発達が単なる等尺性の拡大 (Phase 1) から、領域特異的な高密度化 (Phase 2b) を経て、精緻化・安定化 (Phase 3) に至る3相のプログラムに従うことを全脳規模で初めて実証した。さらに、空間トランスクリプトミクスとの統合解析により、血管発達の駆動因子が、初期の Vegfa 主導の組織低酸素応答から、後期の感覚入力およびシナプス成熟に伴う神経活動依存的シグナル (Apln や Wnt9a) へと動的にシフトする「分子スイッチ」の存在を新規に解明した。特に、Phase 2b の開始時期 (P12-P14) がマウスの開眼や活発な whisking (ひげ運動) の開始といった行動学的・感覚的マイルストーンと時間的に一致することは、神経活動に伴う局所的な代謝需要の増大が血管網の領域特化を駆動するという仮説を強力に支持するものである。
臨床応用: 本研究で得られた知見は、小児神経疾患や発達障害における神経血管相互作用の病態解明に向けた臨床応用の基盤となる。例えば、早産児における脳室内出血や脳白質傷害、あるいは先天性心疾患術後の脳発達遅延において、生後初期の酸素環境の変化や感覚入力の遮断が、Phase 2 の領域特化プログラムにどのような破綻をもたらすかを定量的に評価するための標準対照データとして極めて有用である。また、自閉スペクトラム症である ASD (autism spectrum disorder) などの発達障害モデルにおいて、神経回路の成熟不全が activity-dependent な血管形成の異常を伴うかという問いに対し、本アトラスは translational な解析アプローチを提供し、脳血管発達の破綻を防ぐ治療介入の臨床的意義を提示する。
残された課題: 今後の検討課題として、第一に、本研究で同定された Apln や Wnt9a などの活動依存的リガンドが、in vivo において実際に神経活動依存的な血管形成を直接駆動しているかについての因果関係の検証 (条件的ノックアウトマウスや光遺伝学を用いた検証) が必要である。第二に、マウスで同定された3相の発達プログラムが、ヒトの胎児期から新生児期・乳幼児期におけるどの発達段階に対応するのかという種間ホモロジーの同定が挙げられる。第三に、P30 前後で出現した性的二型な血管分布について、ゴナドトロピンや性ステロイドホルモンが血管内皮細胞や周囲のグリア細胞に与える直接的な作用機序の解明が残されている。これらの limitation を克服することで、LAMBADA アトラスは発達期脳医学および神経血管カップリング研究の強力なプラットフォームとしてさらに発展すると考えられる。
方法
アトラス構築と画像取得: 生後3日目である P3 (postnatal day 3) から21日目である P21 (postnatal day 21) にかけての C57BL/6J マウス (計50個体以上、各タイムポイント n=15 から n=20 mice) の脳を対象とした。iDISCO+ (immunostaining-enabled 3D imaging of solvent-cleared organs) 法を用いて組織を透明化し、光シート顕微鏡を用いて25 μm等方解像度で autofluorescence 信号をスキャンした。得られた画像から平均テンプレートを作成し、Allen CCFv3 (Common Coordinate Framework version 3) に準拠して62の脳領域を手動でアノテーションした。さらに、生体内の脳容積変化を正確に反映させるため、11.7T 超高磁場 MRI スキャンデータを用いて、sigmoidal growth model に基づく容積補正を行い、透明化による組織収縮の影響を排除した。
3次元血管グラフの再構築とAIによる補正: 内皮細胞マーカーである CD31 および Podocalyxin、ならびに動脈壁マーカーである SM22 (smooth muscle 22 alpha、別名 transgelin) に対する共役一次抗体を用いて全脳免疫染色を施した。TubeMap パイプラインを用いて、9つのタイムポイントから計50個の半脳血管グラフ (n=50 reconstructions) を再構築した。この際、画像の局所的な輝度低下や急峻な屈曲に起因する人工的な血管断端 (artefactual interrupted vessels) を識別するため、Video Swin Transformer をベースとした3次元画像分類 AI 検出器を開発し、ClearMap 3.1 ソフトウェアに統合した。
空間トランスクリプトミクスとの統合: P5、P7、P9、P12、P21、P60 のマウス脳から3つの矢状断面を切り出し、10x Genomics 社の Visium 空間トランスクリプトミクス技術を用いて全トランスクリプトームデータを取得した (各断面 n=2 から n=3 replicates)。取得したデータを LAMBADA アトラスの3次元座標系にアラインメントし、遺伝子発現パターンと血管構造パラメーター (血管長密度、分岐点密度、断端密度など) との相関を Spearman’s rank correlation coefficient を用いて算出した。さらに、cNMF (consensus non-negative matrix factorization) 解析を適用し、2,451個の遺伝子から30種の遺伝子プログラム (gene programs) を抽出して血管リモデリングとの空間的・時間的相関を解析した。統計解析には Fisher’s exact test などを用いた。