《PLOS Computational Biology》：Cd47 deficiency in mouse brain under excitotoxic stress
编辑推荐：
CD47是一种跨膜蛋白，在免疫和神经发育中具有公认的作用，近期被发现是人类、非人灵长类和小鼠中丰富的突触前蛋白。为研究Cd47如何影响神经元应激与损伤，研究人员对暴露于惊厥剂量（KA25；25 mg/kg）或亚惊厥剂量（KA3；3 mg/kg）海人酸（kain
CD47是一种跨膜蛋白，在免疫和神经发育中具有公认的作用，近期被发现是人类、非人灵长类和小鼠中丰富的突触前蛋白。为研究Cd47如何影响神经元应激与损伤，研究人员对暴露于惊厥剂量（KA25；25 mg/kg）或亚惊厥剂量（KA3；3 mg/kg）海人酸（kainic acid, KA）的野生型（wild-type, WT）和Cd47敲除（knockout, KO）小鼠的海马进行了单核RNA测序（single-nucleus RNA sequencing, snRNA-seq）。将细胞类型注释与加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA）相结合，揭示了不同的神经元和胶质细胞基因模块，其共表达结构因基因型和KA暴露而异。值得注意的是，某些模块仅在Cd47缺陷细胞中表现出KA相关的共表达结构变化。对这些模块的富集分析突出了翻译、生物能量学和突触相关基因集。使用Wasserstein距离比较转录组-wide分布表明，在所有十五个两两比较中，载体处理的KO细胞与暴露于KA3的WT细胞最为相似。这些发现提示Cd47缺失将海马转录组推向一种类似于轻度兴奋性毒性应激的状态，暗示Cd47是神经元应激适应的候选调节因子，并为检验该作用提供了假说。
论文解读
研究背景与科学问题
神经元应激反应是神经元抵御多种损伤（包括过度兴奋应激、氧化应激、免疫应答和低氧）的基础适应性机制。当这种反应过度或失调时，可导致兴奋性毒性等神经元损伤，并与阿尔茨海默病、帕金森病和肌萎缩侧索硬化等神经退行性疾病相关。因此，阐明调控神经元应激反应的分子调节因子是转化神经科学的核心目标之一。
CD47，又称整合素相关蛋白（integrin-associated protein, IAP），是一种广泛表达的跨膜糖蛋白，在先天免疫系统中作为“别吃我”信号发挥作用：通过与其受体信号调节蛋白α（signal regulatory protein-alpha, SIRPα）结合，抑制巨噬细胞和小胶质细胞对健康宿主细胞的吞噬。越来越多的证据表明CD47在中枢神经系统（central nervous system, CNS）功能中不可或缺：它丰富表达于突触前元件，参与神经突生长和神经元发育，并通过与微胶质SIRPα受体的相互作用保护突触免受补体介导的吞噬清除。在脑缺血实验模型中，Cd47缺陷可减少梗死体积并抑制缺血后炎症。
然而，Cd47在调控神经元对应激的内在反应中的细胞自主作用仍不完全清楚，特别是Cd47缺失如何改变神经元的转录景观及其对直接神经元应激源的易感性尚属未知。回答这一知识空白对于确定CD47信号是否可作为以神经元应激和损伤为特征的疾病的可行治疗靶点至关重要。研究人员因此提出假说：Cd47是神经元对过度兴奋应激和兴奋性毒性反应的关键调节因子。
研究设计与方法
研究人员使用亚惊厥剂量（KA3，3 mg/kg）和惊厥剂量（KA25，25 mg/kg）的海人酸（kainic acid, KA）对C57BL/6野生型（WT）小鼠和Cd47敲除（KO）小鼠（B6.129S7-Cd47tm1Fpl/J，购自Jackson Laboratory）进行腹腔注射，以模拟过度兴奋应激和兴奋性毒性。实验队列共15只小鼠（7只WT和8只KO，8雄7雌，约3月龄，体重20–25 g），分为载体（Veh；每组3只）、KA3（每组3只）和KA25（WT 1只，KO 2只）三组。注射后12小时处死动物，收集海马组织进行核分离，使用10X Genomics Chromium Single Cell 3' v3试剂盒建库并在HiSeq 4000平台测序，经Cell Ranger 7.0.1比对至mm10基因组，随后通过Scanpy进行质量控制（去除空液滴、低质量核和计算双联体），最终64,569个核被分配到七个细胞群体。
研究采用三层次多尺度分析框架：第一层次为单基因差异表达分析（Approach 1），在单核水平和伪bulk水平识别WT与KO之间的差异表达基因；第二层次为基于WGCNA的基因模块分析（Approach 2），构建细胞类型特异的基因共表达网络，通过模块-性状相关性、Jaccard相似性、模块中心性分析和功能富集识别核心模块和枢纽基因；第三层次为全转录组分布相似性分析（Approach 3），使用核密度估计（kernel density estimation, KDE）估计各基因型-处理组合的转录组概率分布，并以归一化1-Wasserstein距离量化两两之间的分布差异。
研究结果
Cd47 KO以背景依赖方式改变兴奋性毒性应激的转录景观（Approach 1）
通过无偏Leiden聚类和标记基因比对，研究人员注释了七个主要细胞类型：CA1、CA3、齿状回（dentate gyrus, DG）、抑制性神经元，以及星形胶质细胞、少突胶质细胞和小胶质细胞。CA3锥体神经元表现出最多的核水平差异，这与CA3海马锥体神经元对KA最易感的形态学数据一致。在基线（载体处理）条件下，Cd47 KO的CA3核较WT更高表达囊泡谷氨酸转运体Slc17a7及大量胞质核糖体蛋白和氧化磷酸化亚基基因（如Rps28和Rpl27），研究人员将其描述为翻译和生物能量学转录本的协调性转变。突触相关基因在两个方向上均有改变：KO中更高表达Kcnj3（GIRK1钾通道亚基）、Negr1（IgLON家族粘附分子）和Tmeff2；更低表达Neurod6（bHLH转录因子）、Mdga1（GPI锚定粘附分子）、Kctd4、Wscd2、Chgb（嗜铬粒蛋白B）和Ttr（转甲状腺素蛋白）。在KA暴露下，基因型对比转向突触传递、可塑性和应激反应相关基因，WT较KO更高表达Tsc1（mTOR通路负调控因子）、Prrt1（AMPA受体转运）、Tmem108（突触囊泡运输）、Tango2和Ttr。在KA3条件下，小鼠水平分析中任何细胞类型均无基因达到统计显著性，这与全转录组排名中KO-KA3和WT-KA3高度相似一致。KA25基因型对比仅作描述性呈现，因为仅含1只WT和2只KO动物，最小可达置换p值为0.333，无法支持推断性结论。在载体条件下，小鼠水平伪bulk检验恢复了核心基线观察：Apoe、Rplp1和Fth1——共表达分析中Cd47相关神经元模块连接度最高的三个基因——在检测到它们的每个细胞类型中KO均高于WT，且无细胞类型显示不一致方向。Rplp1（倍数变化=4.22，校正p=0.046）和Fth1（倍数变化=3.92，校正p=0.046）在星形胶质细胞中达到FDR显著性，其余对比在同一方向上名义显著。
使用WGCNA鉴定和表征细胞类型特异基因模块（Approach 2）
研究人员跨所有基因型-处理组和细胞类型应用WGCNA，共获得110个细胞类型特异模块。模块-性状相关性分析揭示了基因型和处理相关的共表达结构变化。在基因型差异模块中，13个模块在载体处理的KO与WT之间表现出显著不同的模块-性状相关性。神经元模块CA1–1和CA3–1——均以Rplp1、Apoe和Fth1为最高连接基因，且均富集胞质翻译和氧化磷酸化——表现出最大的基因型差异；胶质模块中Astro-1效应最大。在KA响应模块中，亚惊厥KA暴露下两种基因型均显示保守的CA3–10模块-性状相关性增加和DG-12模块-性状相关性降低。然而Cd47状态明显调节了这一模式：Astro-1特征基因在KO而非WT中改变其性状关联，小胶质模块在KO中的变化幅度小于WT，而CA1–1、CA3–1和Astro-1中显著的负向变化仅见于KO。通过Jaccard相似性分析，七个神经元模块对既高度重叠又显著富集。CA1–1和CA3–1表现出最强且最显著的重叠（Jaccard=0.47），两者独立地富集胞质核糖体和氧化磷酸化条目，构成两个锥体群体共有的翻译-生物能量学基因集。功能富集分析揭示了三个主题：第一，翻译-生物能量学基因集将蛋白质合成机制转录本与氧化磷酸化组分转录本联系起来，在CA1–1、CA3–1和Microglia-9中突出；第二，CA3–2模块富集“非运动纤毛组装”基因，其特征基因关联在WT中随处理改变而在KO中不改变，提示CD47可能是纤毛转录反应所必需的；第三，富集突触和神经递质释放蛋白转录本的DG-1和DG-12模块在KA25处理下表现出最大的相关性变化，表明这些转录本在惊厥剂量下共表达最不稳定。小鼠水平线性模型（n=12，限于Veh和KA3条件，85个保留模块）中，仅三个模块达到FDR显著性，均在抑制性神经元中且均为处理项（KA3对载体）：Inhibitory-6（FDR=0.015）、Inhibitory-18（FDR=0.023）和Inhibitory-22（FDR=0.023），各显示KA3下特征基因低于载体。无基因型主效应或基因型-处理交互项在任何细胞类型中通过校正。
共表达枢纽分析鉴定Cd47敏感模块内高度连接基因
通过度中心性量化模块内基因连接度，CA1–1和CA3–1模块的中心性得分比其他神经元模块高出近十倍，表明其内部共表达异常强。两模块在metacell重评估中均高度保留（CA1–1的Zsummary=35.85，CA3–1的Zsummary=25.58）。跨多个神经元群体的泛神经元枢纽基因排名中，Fth1（铁稳态和氧化防御）、Rplp1（翻译控制和应激适应性蛋白质合成）和Apoe（脂质运输和神经元-胶质修复信号）位居前列。研究人员强调，高共表达中心性不建立这些基因的功能重要性或调控作用，仅将其作为功能验证的候选基因。Apoe和Rplp1在Approach 1和Approach 2中同时被识别，两种方法测量同一表达矩阵的不同属性（平均表达差异对共变异结构），其共现具有信息量但并非独立重复。
跨条件变异和分布相似性分析（Approach 3）
在十五个两两比较中，KO-Veh与WT-KA3组的分布最为相似。这一发现表明，Cd47缺失在基础条件下产生的转录组谱与野生型细胞暴露于低剂量KA诱导的轻度兴奋性毒性应激的谱高度接近。与此一致，KA3条件下的基因型对比在每个细胞类型中均为零结果：一旦施加KA3，WT和KO海马转录组在动物水平上不可区分。Wasserstein距离是确定性度量而非样本级检验统计量，每个分布通过扩散流形上的核密度估计整合数千个核的信息。
讨论与结论
研究人员提出，Cd47缺失在基线条件下使海马细胞转录组类似于野生型细胞在轻度兴奋性毒性应激下的状态，提示Cd47可能设定海马神经元启动应激反应的阈值。研究人员类比CD47已确立的“别吃我”功能，提出Cd47依赖性“别应激我”假说作为组织性假说而非已证实的机制。基线状态下，KO动物表现为跨细胞类型的趋同转录差异：神经元和胶质细胞更高表达抗氧化防御和脂质代谢相关转录本（Apoe、Cst3、Mt1、Ptgds），更低表达突触和剪接相关转录本（Ttr、Pla2g3、Doc2b、Srsf7）；小胶质细胞更高表达Apoe、Clu、Mt1、Camk1d，更低表达补体通路转录本（C1qa、Selplg）；少突胶质细胞在髓鞘和细胞骨架转录本（Mylk、Gpm6a、Sept4）上存在差异。在兴奋性毒性挑战下，基因型对比特征改变：KA3下KO神经元更高表达AMPA受体转运（Prrt1）和细胞信号抑制（Sag）相关基因；KA25下描述性图景由Ttr、Mt1、Fkbp5和Apoe主导，伴随Actb和Camk1d等关键细胞骨架和突触组分的下调。KO小鼠的胶质炎症特征相对减弱，与CD47在调节巨噬细胞和小胶质细胞参与中的已知作用一致。
三个枢纽基因的功能假说分别指向：Fth1可能通过铁稳态调节影响KA驱动的氧化应激和铁死亡样损伤的易感性；Rplp1所在的模块可能反映与翻译反应相关的协调调控并贡献于细胞韧性；Apoe相关的脂质运输、修复和清理过程在Cd47缺失背景下因“别吃我”信号缺失而可能改变，脂质组学、吞噬通量测定和细胞因子测量可区分其协调作用与相关性作用。APOE是晚发性阿尔茨海默病最强的全基因组关联研究（Genome-Wide Association Study, GWAS）风险因子，为网络分析与人类疾病遗传学提供了直接汇合点。
研究结论：Cd47缺失在基础条件下将海马转录组推向类似于轻度兴奋性毒性应激的状态，Fth1、Rplp1和Apoe是受Cd47缺失影响最大的基因网络中最高度连接的基因，可作为未来功能验证的候选。研究局限性包括：KA25组样本量不足以支持推断性结论；核来自同一动物不构成独立重复；三种分析方法应用于同一数据集不构成独立生物学证据；缺乏蛋白质水平、形态学或电生理学验证；两种剂量无法解析剂量-反应曲线形状；性别虽在基因型水平大致平衡但无法作为生物学变量正式检验。
