单细胞多组学分析揭示MeCP2磷酸化调控作用的神经元偏向性

《Neurobiology of Disease》：Single-cell multiome profiling reveals neuronal bias of modulatory role of MeCP2 phosphorylation

MeCP2经历多种翻译后修饰，例如响应不同刺激的磷酸化。多个MeCP2磷酸化位点已被单独鉴定和研究。然而，多位点磷酸化的组合效应尚未被直接检验。基于先前对MeCP2磷酸化位点单突变的研究，研究人员创建了一个携带两种突变的小鼠模型：丝氨酸80变为丙氨酸和丝氨酸421变为谷氨酸（S80A;S421E，或A;E），其中MeCP2蛋白磷酸化被固定在活跃神经元中观察到的状态。进行了一系列行为测试以表征整体动物水平的功能结果，随后进行了两个独立的单细胞/单核水平多组学分析，以表征改变的分子通路，以及在转录和染色质水平上偏向谷氨酸能神经元的调控模式。虽然对多维数据集的整合分析确定了由MeCP2磷酸化调控的细胞类型特异性分子靶点，但空间分辨的转录分析为靶点验证提供了一个正交平台，并从空间角度提供了额外信息。总之，这项研究的结果表明，MeCP2 S80和S421磷酸化的组合效应不仅仅是协同的，而且调控偏向于谷氨酸能兴奋性神经元。更重要的是，这项研究首次探索了MeCP2磷酸化突变模型中单细胞多组学图谱的变化，并将其与功能结果联系起来。

**论文解读：单细胞多组学揭示MeCP2磷酸化对神经元调控的偏向性**

**研究背景、现存问题与研究动机**
甲基-CpG结合蛋白2（MeCP2）是一种关键的表观遗传调控因子，在中枢神经系统（CNS）中高度表达，其功能异常导致Rett综合征（RTT），一种以认知缺损、运动障碍和自闭症样特征为主要表现的神经发育障碍。MeCP2可经历多种翻译后修饰，其中磷酸化修饰最为重要，尤其是丝氨酸80（S80）和丝氨酸421（S421）位点：S421在神经元激活后被磷酸化，而S80在静息神经元中组成型磷酸化并在激活后去磷酸化。尽管已有研究分别探讨了S80A（去磷酸化突变）和S421E（模拟磷酸化突变）单突变的影响，但这两个位点同时发生修饰时的组合效应尚未被直接探究。因此，研究人员构建了携带S80A和S421E双突变的小鼠模型，以模拟活性神经元中MeCP2的固定磷酸化状态，并从整体行为、电生理、单细胞多组学及空间转录组层面系统解析其功能后果。

**研究内容、结论与意义**
该研究通过一系列行为测试、全细胞膜片钳记录、单核多组学（snMultiome）测序、多路误差校正荧光原位杂交（MERFISH）成像及计算分析，综合揭示了Mecp2^S80A;S421E（A;E）突变小鼠中谷氨酸能兴奋性神经元偏倚的转录和染色质可及性改变，并发现该双突变产生的表型并非单突变的简单叠加，而是呈现出协同、拮抗及涌现三种不同模式。研究成果发表在《Neurobiology of Disease》，首次在单细胞分辨率下描述MeCP2磷酸化突变的多组学图谱，并建立了分子改变与功能结局之间的关联，为理解MeCP2磷酸化的组合编码机制及其在神经发育障碍中的作用提供了重要资源。

**主要关键技术方法（不超过250字）**
研究人员利用以下核心技术（样本来源：雄性Mecp2^S80A;S421E/y小鼠及其野生型同窝对照，来自美国威斯康星大学麦迪逊分校）：（1）行为学测试（包括跑步轮、旷场、恐惧条件化等）；（2）急性海马脑片CA1锥体神经元全细胞膜片钳记录；（3）基于10× Genomics平台的单核多组学（snMultiome，同时检测RNA表达和ATAC-seq）分析，对3对9周龄小鼠海马进行测序；（4）使用SCENIC框架进行基因调控网络（GRN）推断；（5）MERFISH成像（300基因panel）独立验证转录变化并获取空间信息；（6）Tangram算法将snRNA-seq数据映射到MERFISH空间坐标，进行亚区域特异性分析。

**研究结果**

**3.1 Mecp2 S80A;S421E小鼠表现出活动减少、中央区域探索增加和情境记忆增强**
通过行为测试发现：A;E小鼠在跑步轮暗期活动减少（p=0.05），旷场总移动距离和垂直活动减少（p<0.001），但中央区域探索时间显著增加（p=0.032），且高架十字迷宫和明暗箱实验未发现焦虑差异；恐惧条件化测试中情境记忆显著增强（p<0.001），而线索记忆无变化，提示海马功能改变。

**3.2 电生理分析显示Mecp2 S80A;S421E海马CA1锥体神经元自发性兴奋性突触后电流（sEPSC）增强**
全细胞膜片钳记录发现，A;E小鼠CA1神经元sEPSC频率（p=0.021）和振幅（p=0.011）均显著增加，而静息膜电位、输入电阻、放电阈值等固有膜特性无显著差异，提示突触前活动增强和突触后反应增加。

**3.3 单核基因表达谱揭示Mecp2 S80A;S421E海马中谷氨酸能神经元偏倚的转录变化**
snRNA-seq分析海马60,889个核，鉴定出7种主要细胞类型。差异表达基因（DEG）高度富集于神经元，尤其是谷氨酸能神经元：DG-IMN Glut（4761个）、IT-ET Glut（3496个），而非神经元（如星形胶质细胞仅108个）。GRN分析显示，IT-ET Glut中27个调控子（regulon）在基因型间差异富集，其中Fos、JunB、Egr1等与记忆相关的调控子偏向野生型（WT），而DG-IMN Glut中49个调控子差异富集，Foxp1偏向A;E。GO富集分析表明，偏向WT的调控子富集于“细胞通讯调节”、“学习或记忆”，而偏向A;E的调控子富集于“突触组织”、“轴突导向”。细胞-细胞相互作用分析发现A;E中有更多互作对（18对 vs 14对），且涉及非神经元的互作权重更高。

**3.4 Mecp2 S80A;S421E海马染色质可及性特征与基因表达数据集一致**
snATAC-seq分析发现差异可及性峰（DAP）同样富集于谷氨酸能神经元，尤其IT-ET Glut（2589个DAP），且大多数DAP在A;E中可及性增加。峰注释显示多数DAP位于转录起始位点（TSS）附近。峰-基因连接分析鉴定了Fos和Cck两个基因：两者在A;E中基础表达降低，但染色质可及性增加，且存在相应转录因子（Egr1、Zbtb18）的结合基序，提示可能形成“预备”染色质状态以增强环境刺激响应。

**3.5 MERFISH验证转录变化并利用空间组织信息增强细胞分型**
对3对小鼠海马进行300基因MERFISH成像，获得50,409个细胞。超过半数的10×来源DEG被验证（IT-ET Glut中45%验证，DG-IMN Glut中85%验证），Cck等关键GRN靶点也得到验证。进一步将IT-ET Glut分为CA1、CA2、CA3、Mossy亚区，发现CA1区域DEG最多（132个），且10×数据通过Tangram映射后同样显示CA1区域DEG和DAP最多，说明主要分子改变集中于CA1。

**讨论与结论**
讨论部分指出：双突变的表型并非单突变效应的简单相加——在跑步轮实验中呈现协同效应（均减少活动），在旷场中涌现新表型（单突变无变化，双突变增加中央探索），在恐惧条件化中则表现为拮抗（单突变减少记忆，双突变增强记忆）。这表明多个磷酸化位点之间存在复杂的交叉对话机制，可能通过补偿性重编程或缓冲能力超过阈值而暴露潜在脆弱性。单细胞多组学揭示的谷氨酸能神经元偏倚与S80/S421磷酸化动态响应神经元活动的已知功能一致。染色质可及性增加与转录水平变化之间的解离（如Fos和Cck基础低表达但可及性增加）提示一种“预备”表观遗传状态，而非持续转录激活。MERFISH验证增强了结果可靠性，同时平台差异（测序 vs 成像、空间采样范围不同）解释了一部分不一致。
研究结论翻译：总之，本研究证明了Mecp2 S80A;S421E联合磷酸化突变产生复杂的表型结局。在分子水平上，研究人员观察到Mecp2 S80A;S421E海马中谷氨酸能神经元偏倚的转录和染色质变化。通过使用单核多组学和独立的MERFISH验证，研究提供了细胞类型分辨的MeCP2功能见解，并提供了高质量数据和有价值的跨平台比较。这些发现增进了对组合翻译后修饰如何调节MeCP2功能的理解，并强调了研究神经发育障碍中多位点调控的重要性。