《Cells》：Physiological Reference Mapping Reveals Context-Dependent Dysregulation of Regenerative Programs in Duchenne Muscular Dystrophy
编辑推荐：
杜氏肌营养不良（Duchenne muscular dystrophy，DMD）涉及慢性变性与再生，导致营养不良性肌肉干细胞（dystrophic muscle stem cells，MuSCs）与未损伤野生型组织存在差异。常规的疾病-静息对照比较可能将失调与
杜氏肌营养不良（Duchenne muscular dystrophy，DMD）涉及慢性变性与再生，导致营养不良性肌肉干细胞（dystrophic muscle stem cells，MuSCs）与未损伤野生型组织存在差异。常规的疾病-静息对照比较可能将失调与再生混为一谈。研究人员开发了一种生理参考映射框架（physiological reference-mapping framework），整合正常再生轨迹与可重复的营养不良数据集，以估计相对于再生状态和空间背景的疾病相关偏差。研究人员将该框架应用于单细胞/单核、批量及空间RNA测序数据集。在一个可重复的基因型×状态MuSC数据集中，包含野生型（wild-type，WT）和肌营养不良蛋白缺失的mdx小鼠样本，覆盖四个谱系状态，每组三个生物学重复，激活的mdx MuSC表现出过度的细胞外基质（extracellular matrix，ECM）诱导（标准化mdx-WT差异=1.431；FDR=3.69×10?9）以及减弱的MYC相关应答（?0.731；FDR=2.57×10?4），此为五个预设程序中的两项。空间分解表明，在改变的生态位组成之外，还存在显著的生态位内重塑。生理校准将ECM诱导确定为一种在营养不良中变得过度的正常短暂程序，并提示相对于静息肌肉的氧化磷酸化转录降低可能代表相对于正常再生的抑制不足。一个独立的人类肌营养不良症空间队列支持整个活检组织中ECM转录的增加。因此，生理参考映射有助于区分再生适应与DMD中背景依赖性疾病的失调。
论文解读：生理参考映射框架下的杜氏肌营养不良转录失调再解析
一、研究背景与问题
杜氏肌营养不良（Duchenne muscular dystrophy，DMD）由功能性肌营养不良蛋白缺失引起，以反复肌纤维变性、炎症、纤维化和不完全再生为特征。肌肉干细胞（muscle stem cells，MuSCs）的静息-激活-增殖-分化转变与周围微环境变化紧密协调。然而，营养不良肌肉中测得的转录组不仅反映原始遗传缺陷，还包含慢性再生活动、细胞组成变化及组织微环境重塑。此前研究多将DMD组织与未损伤野生型（wild-type，WT）组织直接比较，但这种比较可能混淆两类截然不同的生物学过程：沿正常再生轨迹的位置移动，以及真正与疾病相关的转录偏离。由于正常再生本身是一个动态、异质的过程，不同阶段具有不同的预期分子特征，因此某一通路相对于静息WT肌肉表达降低可能恰是激活状态下的正常表现，而某通路在正常修复中升高则仅在幅度、时间、持续性或定位超出生理范围时才具有病理意义。此外，空间转录组数据中组织切片间的差异可能源于生态位丰度改变或匹配生态位内表达改变，而将单个细胞或空间点视为独立生物学重复会进一步高估证据强度。因此，如何将疾病相关失调与正常再生状态位移分离，成为解读DMD转录组研究的核心问题。
二、研究目标与总体设计
研究人员开发了一种生理参考映射框架，将观察到的营养不良转录程序分解为匹配再生状态下的预期正常程序加上背景依赖性（包括细胞状态、再生时间、空间生态位、组织与遗传背景）的疾病相关偏差。该框架应用于多个独立小鼠单细胞和单核数据集、一个可重复的基因型×状态MuSC实验、正常损伤时间进程、空间转录组切片及人类肌营养不良症队列。统计推断均在动物、文库或供体层面进行，而非细胞或空间点层面。研究目的在于确定哪些转录异常在考虑正常再生状态和空间背景后仍然存在，而非提名某单一通路作为DMD的决定性机制。
三、关键技术方法
研究人员使用的主要方法包括：第一，单细胞/单核RNA测序数据处理，包括质量过滤、双细胞去除、无批次整合的图聚类及基于经典标记的细胞注释；第二，基于伪bulk的差异表达分析，使用edgeR拟似然负二项模型，以生物学重复（小鼠或供体）为统计单位；第三，正常生理参考构建与跨平台状态映射，利用正常急性激活参考时间序列，采用样本内秩次而非原始表达量进行比对；第四，预定义通路评分与CAMERA基因集检验，分析五个固定基因程序；第五，空间转录组生态位匹配与混合效应模型，将空间点按验证过的生态位归类，以小鼠×生态位汇总为推断单位；第六，敏感性与稳健性分析，包括leave-one-out重拟合、异方差一致性标准误及精确置换检验。样本队列来源包括：小鼠单核数据集PRJNA771932（3只WT、4只mdx）、可重复基因型×状态数据集GSE268313（24个文库，每个基因型-状态组3个生物学重复）、空间数据集GSE225766（3只WT DBA/2J、5只D2-mdx）、人类空间队列GSE288958（3个对照供体、8个肌营养不良症供体），以及基因编辑数据集GSE180771。
四、研究结果
1. 营养不良性MuSC沿生理激活轴偏离静息状态。对PRJNA771932的MuSC/肌原性区室分析识别出稳态静息MuSC、营养不良Pax7阳性MuSC、Runx1阳性激活祖细胞、循环肌原性祖细胞和分化中肌原性细胞。生物学重复水平上，3只WT小鼠以稳态静息MuSC为主，而4只mdx小鼠含有不同比例的激活、循环和分化状态。与独立正常激活参考（T0静息、T2早期激活、T4激活/定型）比对显示，WT样本最相似于T0，而多个mdx样本更相似于T2或T4。这表明传统mdx对WT比较会混入沿生理激活轴的位移。
2. 状态匹配分析揭示过度ECM诱导和减弱的MYC相关转录应答。在GSE268313数据集中，主成分分析显示谱系状态解释了主要转录变异（PC1 52.2%，PC2 25.9%）。在五个预设程序中，激活的mdx MuSC表现出显著过度的ECM程序诱导（标准化差异1.431，FDR=3.69×10
?9），以及减弱的MYC相关激活应答（差异?0.731，FDR=2.57×10
?4）。留一文库敏感性分析在所有24次拟合中均保留效应方向和显著性，证明结果不依赖单个文库。
3. 疾病相关转录重塑在匹配的生态位内持续存在。空间转录组分析识别出收缩、ECM/FAP、血管周、再生和免疫富集生态位。营养不良组织表现出明显的组成重塑，包括再生和免疫区域扩大、收缩区室减少。对组织水平差异的分解显示，生态位内重塑占ECM对比的84%、mTORC1的87%、MYC的94%、TNF/NF-κB的86%和OXPHOS的62%。匹配生态位分析中，D2-mdx肌肉在收缩生态位（差异0.938，FDR=0.00105）、ECM/FAP生态位（0.803，FDR=0.00273）和血管周生态位（1.182，FDR=0.000389）均显著升高ECM评分。留一鼠分析保持阳性结果。
4. 生理损伤轨迹区分再生反应与疾病相关偏差。在正常损伤后第1、3、5天，再生富集生态位中的ECM评分逐渐增加，说明ECM诱导本身是正常早期再生特征。然而，5只D2-mdx小鼠的ECM评分超过正常损伤采样区间范围。五个胶原基因（Col1a1、Col1a2、Col3a1、Col5a1、Col5a2）均在D2-mdx中一致升高（FDR=0.037）。OXPHOS在正常损伤中呈现早期抑制（第1-3天）后恢复（第5天）的模式，且在三个独立正常损伤空间数据集中均观察到类似恢复趋势。这改变了静态疾病对照对比的解释：相对于正常再生预期，营养不良肌肉中OXPHOS降低可能代表抑制不足而非功能缺陷。
5. 人类全活检ECM关联在供体校正后保留。在人类肌营养不良症空间队列中，以供体为统计单位，五个预设程序中仅ECM通过校正（差异1.718，精确p=0.00606，BH q=0.0303）。五个胶原基因（COL1A1、COL1A2、COL3A1、COL5A1、COL5A2）均升高（q=0.00758–0.0121）。在MuSC富集和再生富集区域中ECM也更高，但未达显著（均q=0.0848）。小鼠单核全图谱显示基质细胞占捕获胶原转录本的88.3–92.5%，但肌原性区室中也可检出胶原转录。
6. 额外分析。激活的mdx样本显示Foxo1和Foxo3表达升高及推断靶程序评分升高，FOXO1结果在排除共享靶标后仍稳健。38个免疫-卫星细胞互作对在多重校正后无基因型差异。e23AON治疗未达显著；而AAV9-CRISPR/Cas9介導的肌营养不良蛋白校正后，4周时OXPHOS升高（q=0.0208）、TNF/NF-κB降低（q=0.0208），18个月时ECM降低（q=0.00533）、TNF/NF-κB降低（q=0.00533），方向与疾病相关差异相反。
7. 已发表发现的结构化生理再解释。对7篇营养不良肌肉研究的79条声明进行结构化审查，其中48条有合格生理参考并得以分类：19条（40%）被确认但需细化生物学解释，16条（33%）受空间或组成条件限制，9条（19%）无需实质修改即被确认，4条（8%）被生理性重构。其余31条因缺乏适当参考、重复或状态匹配而保留为未解决。
五、讨论与结论
讨论部分指出，生理参考匹配揭示多种表观营养不良相关变化实为沿正常再生程序的位移，而另一些则作为状态或生态位匹配后的持续疾病偏差存在。ECM是生理程序在营养不良中变得定量和空间异常的清晰示例：正常再生中短暂出现的程序在营养不良中表现为幅度过大或持续存在，并在小鼠跨生态位及人类活检中一致再现。空间匹配证明这些效应不能仅由组织组成变化解释，生态位内重塑贡献了主要部分。MYC应答减弱表明激活程序的不完整执行，OXPHOS结果则提示静态比较可能误导，需以匹配再生状态的预期为参照。FOXO因子表达升高可能反映激活与干细胞池维持之间的协调异常。作者强调研究局限性：队列样本量小、年龄和肌肉差异、正常损伤参考采样稀疏、人类数据为横断面等，因此功能验证仍属必要。结论认为，正常再生参考可细化对营养不良肌肉转录变化的解读，该框架可区分生理状态位移与疾病偏差、分离组成重塑与生态位内失调，并保留未解决或阴性结果。应用于DMD，过度ECM重塑是跨平台、跨物种最可重复的残留异常，而其他表观异常则依赖于细胞状态、空间背景或所选生理对照。
