《Nature Communications》：Vascular smooth muscle cell state trajectories mediate molecular mechanisms of coronary disease risk
编辑推荐：
本研究旨在探究血管平滑肌细胞(VSMC)在动脉粥样硬化（AS）中的复杂状态转换，及其介导遗传性冠状动脉疾病(CAD)风险的分子机制。研究人员通过建立密集时间点的单细胞转录组与表观基因组图谱，结合空间转录组学，识别了驱动细胞状态转变的关键转录因子(如Tcf21)，并揭示了其在人类CAD风险中的因果机制，为理解AS病理进程提供了高分辨率的视角。
冠状动脉疾病是全球范围内导致死亡的主要原因之一，其病理基础动脉粥样硬化是一个复杂的动态过程。在这个过程中，血管壁中的平滑肌细胞（vascular smooth muscle cell, VSMC）扮演着关键但充满谜团的角色。传统观点认为，VSMC主要是构成血管壁的“稳定”结构成分，然而越来越多的证据显示，在疾病刺激下，它们会表现出惊人的“变身”能力，可以转变成多种与疾病相关的表型，如合成型、泡沫细胞样或成骨细胞样细胞，从而参与斑块的形成、稳定与破裂。这种细胞状态的复杂转变轨迹究竟由哪些遗传因素控制？它们又是如何影响个体罹患冠心病的风险？长期以来，由于缺乏在疾病发展全程中对VSMC进行高分辨率、多组学的动态追踪，这些问题一直悬而未决。
为了深入解析动脉粥样硬化的细胞演变图谱并追溯其遗传根源，一个研究团队在《Nature Communications》上发表了一项系统性研究。他们利用小鼠疾病模型，构建了一个覆盖疾病进程密集时间点的单细胞转录组和表观基因组全景图谱，并结合了高plex空间原位数据，如同一幅描绘细胞命运变迁的“高清动态地图”。通过结合时序数据与概率性命运建模，研究人员像侦探一样，追踪了VSMC在动脉粥样硬化斑块发展过程中的状态演变路径。他们识别出驱动这些细胞状态转变的关键“指挥家”——转录因子，其发现得益于基于网络的优先级排序和计算机模拟的转录因子扰动分析。其中，一个已知的冠状动脉疾病风险基因Tcf21被重点研究。通过平行开展的基因敲除实验，团队揭示了Tcf21在调控VSMC，特别是源自第二心区（secondary heart field）的平滑肌细胞表型转换中的具体分子机制。最终，通过将小鼠图谱与人类冠状动脉疾病遗传数据进行整合分析，研究精准定位了那些真正介导疾病风险的VSMC表型，从而阐明了从遗传变异到细胞功能异常，再到疾病表型的因果生物学机制。这项研究不仅以前所未有的单细胞分辨率解析了动脉粥样硬化的发展轨迹，更重要的是，它搭建了一座连接遗传风险与细胞分子机制的桥梁，为未来开发针对特定细胞状态或通路的精准治疗策略提供了新的靶点和思路。
研究人员运用了多项关键技术构建和分析数据：1. 密集时间点取样的小鼠动脉粥样硬化模型，用于捕获疾病动态过程；2. 单细胞RNA测序（scRNA-seq）和单细胞表观基因组测序（如scATAC-seq），以在单个细胞水平解析基因表达和染色质可及性；3. 高plex多重荧光原位杂交或类似的空间转录组技术，在组织原位确定不同细胞类型的空间分布与关系；4. 计算生物学方法，包括细胞命运轨迹推断、基因共表达网络分析和基于网络的转录因子调控优先级排序；5. 计算机模拟的转录因子扰动预测，以评估转录因子对细胞状态的影响；6. 体内基因功能验证，如对候选基因Tcf21进行条件性敲除，并结合谱系追踪技术，在动物模型中验证其功能。
研究结果
构建动脉粥样硬化的单细胞多组学动态图谱
研究人员首先建立了一个载脂蛋白E敲除（ApoE-/-）的小鼠动脉粥样硬化模型，并在疾病发展的多个关键时间点采集主动脉组织进行单细胞RNA测序和单细胞染色质可及性测序。这产生了一个高分辨率的细胞图谱，涵盖了包括VSMC、内皮细胞、免疫细胞在内的多种细胞类型。通过分析，他们绘制了VSMC在疾病进程中从收缩表型向多种病变相关表型（如纤维肌原细胞样、炎症样、泡沫细胞样）分化的详细轨迹。
识别驱动VSMC状态转变的关键转录调控网络
利用时序单细胞数据，团队应用了细胞命运概率建模算法，预测了VSMC向不同终点状态分化的倾向。通过整合基因共表达网络与染色质可及性数据，他们构建了状态特异的基因调控网络。在此基础上，采用计算模拟扰动的方法，系统性地筛选了对细胞状态转换影响最大的转录因子。分析锁定了一批核心转录调控因子，其中就包括已知的冠心病风险基因Tcf21。
揭示Tcf21在VSMC表型转换中的机制
为了验证计算预测，研究对Tcf21进行了功能缺失研究。在动脉粥样硬化模型中，VSMC特异性敲除Tcf21显著改变了斑块的细胞组成，减少了某些病变相关的VSMC亚群，同时增加了其他亚型。进一步的机制探究发现，Tcf21部分通过调控源自第二心区的VSMC的特定转换程序来发挥作用。这表明，遗传风险可能通过影响特定发育起源的VSMC亚群在疾病中的重编程过程来介导疾病易感性。
整合小鼠图谱与人类遗传学定位致病性VSMC表型
最后，研究将小鼠模型中鉴定出的VSMC状态特异性基因表达特征与大规模人类冠状动脉疾病全基因组关联研究（GWAS）数据相结合。通过遗传共定位和表达数量性状位点（eQTL）分析，他们发现，人类CAD风险遗传信号显著富集于小鼠模型中特定的疾病相关VSMC状态（如炎症样状态）的调控区域。这直接证明了这些特定的细胞表型转换是连接人类遗传变异与冠心病发病的关键中间环节，从而阐明了其因果生物学机制。
结论与意义
该研究通过构建动脉粥样硬化进展的密集时间点单细胞多组学图谱，系统性地解析了血管平滑肌细胞复杂的状态转变轨迹。研究不仅描绘了细胞命运的动态变化图景，更重要的是，通过计算生物学与实验验证相结合，成功识别出驱动这些转变的关键转录因子（如Tcf21），并深入揭示了其分子作用机制。尤为关键的是，通过跨物种整合分析，研究将小鼠模型中的细胞表型与人类冠心病遗传风险直接关联，精准定位了那些真正介导疾病风险的VSMC表型，从而在单细胞分辨率上阐释了遗传风险影响冠心病发病的因果链条。这项工作的意义在于，它将传统的遗传关联信号“落地”到了具体的细胞生物学过程，为理解动脉粥样硬化的发病机理提供了前所未有的高分辨率动态视角，并为未来开发针对特定致病细胞状态或通路的精准治疗策略奠定了坚实的理论基础。
