《Frontiers in Immunology》：Integrative single-cell and bulk transcriptomic analysis identify R-loop-associated gene signatures driving immune-inflammatory networks in acute pancreatitis
编辑推荐：
摘要专业翻译
背景： R-loop失调已被证实可通过DNA损伤和细胞质单链DNA（ssDNA）释放参与慢性炎症过程，然而其在急性胰腺炎（AP）中的作用尚未被系统研究。本研究旨在通过整合生物信息学分析和体外实验验证，鉴定AP中R-loop相关生物标志物，并阐明其
摘要专业翻译
背景： R-loop失调已被证实可通过DNA损伤和细胞质单链DNA（ssDNA）释放参与慢性炎症过程，然而其在急性胰腺炎（AP）中的作用尚未被系统研究。本研究旨在通过整合生物信息学分析和体外实验验证，鉴定AP中R-loop相关生物标志物，并阐明其在免疫炎症网络中的潜在作用。
方法： 整合来自GEO数据库的单细胞和批量转录组数据集，采用多种评分算法计算各细胞类型的R-loop相关基因（RLRG）评分。应用高维加权基因共表达网络分析（hdWGCNA）、差异表达分析及多种机器学习算法（LASSO、SVM-RFE和Boruta）筛选候选基因特征。开展免疫浸润和功能富集分析以探究所鉴定基因特征的调控机制。建立雨蛙素诱导的AR42J细胞模型进行实验验证，并通过siRNA介导的基因敲低评估其功能作用。
结果： 成纤维细胞被鉴定为RLRG评分升高且在AP组与对照组间差异显著的主要细胞类型。拟时序分析揭示AP促进成纤维细胞向终末分化方向进展，高RLRG评分细胞表现出更高的干性和细胞间通讯活性。研究鉴定出三个候选基因特征——Anp32b、Myl12a和Nras，其在AP队列中表达显著升高。这些基因与免疫细胞群体显著相关，并富集于免疫炎症、凋亡、细胞骨架和蛋白酶体通路。在雨蛙素诱导的AR42J细胞模型中证实了Anp32b和Myl12a的上调表达。siRNA介导的Myl12a敲低可减弱胰蛋白酶活化并减少炎症细胞因子释放，同时S9.6抗体免疫荧光染色显示Myl12a敲低显著减少核内R-loop积累。
结论： 本研究结果揭示了R-loop相关基因特征，并通过功能验证确定Myl12a为解析AP发病机制的潜在候选分子。本研究为AP相关免疫炎症提供了机制线索，有待进一步研究验证。
论文解读
一、研究背景与意义
急性胰腺炎（acute pancreatitis, AP）是一种以剧烈腹痛、淀粉酶和脂肪酶升高及/或器官功能障碍为特征的胰腺炎症性疾病，全球总体发病率约为每10万人33.74例，且存在显著的区域差异。胆结石、酒精滥用和高甘油三酯血症是AP的主要病因，其中胆结石诱导的AP在全球最为常见。约20%–30%的AP患者可进展为中重度胰腺炎并常伴随并发症，在感染性胰腺坏死病例中死亡率可高达30%。因此，深入探究AP的分子机制、明确调控炎症激活和组织损伤的关键通路具有重要的临床意义。
R-loop是哺乳动物基因组中普遍存在的动态三链核酸结构，参与DNA修复、细胞增殖分化、RNA转录加工、基因表达调控和DNA甲基化等多种生物学过程。由DNA损伤引发的R-loop异常积累可导致细胞质单链DNA释放并激活cGAS-STING通路，从而促进慢性炎症和组织损伤。R-loop失调已被证实与多种无菌性炎症疾病相关，包括神经炎症和线粒体依赖性上皮坏死驱动的自发性肠道炎症。值得注意的是，AP以腺泡细胞DNA损伤、大量胞质核酸释放和强烈的cGAS-STING介导的无菌性炎症为特征，这些病理事件与病理性R-loop积累的下游后果高度重叠。然而，尽管R-loop与胰腺炎之间的关联日益受到关注，但R-loop相关基因特征及调控AP发病的关键介质仍不明确，尚无系统性研究表征AP中R-loop相关生物标志物及其功能作用。本研究旨在填补这一知识空白。
二、研究内容与主要技术方法
研究人员整合了来自GEO数据库的多个AP相关数据集，其中GSE109227和GSE65146合并为训练队列，GSE121038和GSE3644合并为验证队列，GSE181276和GSE198183合并为单细胞队列。研究采用Seurat包进行单细胞数据处理和细胞类型注释，通过irGSEA包中的AUCell、UCell、singscore、ssGSEA和AddModuleScore五种算法计算各细胞类型的RLRG评分。利用hdWGCNA进行高维加权基因共表达网络分析，结合limma包差异表达分析和LASSO、SVM-RFE、Boruta三种机器学习算法筛选候选基因特征。研究还运用了Monocle拟时序分析、CytoTRACE干性评分、Velocyto RNA速率分析和CellChat细胞通讯分析。免疫浸润分析采用ImmuCellAImouse包评估33种免疫细胞的浸润丰度。此外，研究构建了转录因子调控网络和lncRNA-miRNA-mRNA调控网络，并通过DsigDB数据库进行药物预测和AutoDock Vina分子对接。实验验证部分采用雨蛙素诱导的AR42J大鼠胰腺外分泌细胞建立AP体外模型，通过siRNA介导的Myl12a基因敲低进行功能研究。
三、研究结果
3.1 单细胞转录组分析鉴定成纤维细胞为关键细胞
经过质量控制后，聚类分析获得18个细胞簇并注释为11种细胞类型：内分泌细胞、红细胞、星状细胞、成纤维细胞、内皮细胞、间皮细胞、导管细胞、腺泡细胞、巨噬细胞、B细胞和T细胞。在对照组中腺泡细胞最为丰富，而AP组中成纤维细胞最为丰富。在各类细胞中，成纤维细胞表现出最高的RLRG评分且AP组与对照组间差异显著，因此被定义为关键细胞。拟时序分析显示AP促进成纤维细胞沿无分支轨迹向终末分化方向进展，高RLRG评分细胞在早期阶段表现出更低的化程度和更高的干性特征。细胞通讯分析表明高评分成纤维细胞与其他细胞类型表现出更多的相互作用数量和更强的互作强度，关键配体-受体对包括Mif-(Cd74+Cxcr4)（成纤维细胞→B细胞）、Mif-(Cd74+Cd44)（成纤维细胞→巨噬细胞）和Csf1-Csf1r（成纤维细胞→巨噬细胞），MIF、SPP1、PTN、TGFb和MK信号通路活性显著增强。
3.2 机器学习鉴定候选基因特征
hdWGCNA鉴定出四个与成纤维细胞RLRG评分显著相关的关键模块（绿色、红色、grey60和品红色），共获得1010个关键模块相关基因。差异表达分析鉴定出1006个DEGs，其中上调903个、下调103个。两者交集获得79个基因，显著富集于506个GO条目和26条KEGG通路，包括肌动蛋白丝组织、黏附连接、肌动蛋白结合和肌动蛋白细胞骨架调控等。其中38个基因在训练和验证队列中均表现出显著差异表达且趋势一致。LASSO、SVM-RFE和Boruta三种算法分别鉴定出4、15和31个特征基因，三者交集获得三个共享基因：Anp32b、Myl12a和Nras，在AP组中均显著上调。这三个基因在成纤维细胞、内分泌细胞和腺泡细胞中表现出显著的组间差异表达，且在成纤维细胞分化轨迹起点处表达较高。
3.3 免疫微环境与候选基因特征的关联
免疫浸润分析显示33种免疫细胞中有28种在AP组与对照组间表现出显著浸润差异，包括B细胞、树突状细胞和巨噬细胞。相关性分析表明三个候选基因特征与这些免疫细胞群体显著相关，其中Myl12a与巨噬细胞呈最强正相关（Cor = 0.912），与自然杀伤细胞呈最强负相关（Cor = -0.819）。
3.4 候选基因特征的功能富集
GSEA分析显示Anp32b、Myl12a和Nras分别与71、70和70条KEGG通路显著相关，与11、10和10条HALLMARK通路显著相关。三个候选基因特征共享功能相似性，均显著关联于免疫炎症相关通路（如B细胞受体信号通路、NOD样受体信号通路、IL-6-JAK-STAT3信号）、凋亡和细胞存活相关通路（如凋亡、p53信号通路、PI3K-AKT-mTOR信号）、细胞骨架和细胞黏附相关通路（如黏附连接、ECM-受体互作）以及蛋白酶体和蛋白分泌相关通路。
3.5 GSVA揭示广泛的通路失调
GSVA分析显示AP组中免疫和炎症相关通路显著激活，包括B细胞受体信号通路、Toll样受体信号、FcγR介导的吞噬作用、NOD样受体信号和IL-6-JAK-STAT3信号。细胞骨架和细胞黏附相关通路（如肌动蛋白细胞骨架调控、黏附连接、紧密连接）和蛋白酶体通路同样被激活，而缬氨酸、亮氨酸和异亮氨酸降解等代谢通路则被抑制。
3.6 调控网络构建、药物预测与分子对接
预测获得9个转录因子，其中MAX、USF1、SRF和REST与所有候选基因特征均有互作。lncRNA-miRNA-mRNA调控网络包含3个mRNA、36个miRNA和8个lncRNA，代表性调控对包括Oip5os1-mmu-let-7c-5p-Nras、Meg3-mmu-miR-361-5p-Myl12a和Mir17hg-mmu-miR-673-5p-Anp32b。共预测到29种药物与候选基因特征存在互作，其中维生素E与Anp32b和Nras均有互作。分子对接结果显示Clofarabine-Anp32b（结合能-5.493 kcal/mol）和temozolomide-Myl12a（结合能-5.498 kcal/mol）的结合相对稳定。
3.7 Myl12a的表达与功能验证
在雨蛙素诱导的AR42J细胞AP模型中，α-淀粉酶、IL-6和TNF-α水平升高证实模型诱导成功。Anp32b和Myl12a在AP组中显著上调，与生物信息学结果一致，而Nras无显著差异。选择Myl12a进行后续功能分析，siRNA敲低后胰蛋白酶活性较AP模型组显著降低，提示胰腺酶活化减弱。S9.6抗体免疫荧光染色显示Myl12a敲低后核内R-loop水平显著降低。
四、讨论与结论
成纤维细胞在AP中扮演"双刃剑"角色：一方面通过NF-κB-HIF-1α信号通路增强糖酵解依赖的CXCL1产生，招募CXCR2阳性中性粒细胞浸润并加剧胰腺炎症；另一方面通过旁分泌胰岛素样生长因子I（IGF-I）增强腺泡细胞增殖，发挥组织修复功能。本研究发现AP胰腺组织中成纤维细胞比例显著增加，证实成纤维细胞是参与AP发生发展的关键细胞群体，且与R-loop稳态失衡密切相关。
三个候选基因特征Anp32b、Myl12a和Nras均与成纤维细胞功能存在内在生物学联系。ANP32B与成纤维细胞中I型胶原生物合成相关，MYL12A介导成纤维细胞的收缩表型，NRAS参与成纤维细胞活化、增殖和促纤维化特征获得。体外AP模型中Myl12a敲低显著降低胰蛋白酶活性并减少R-loop积累，提示Myl12a可能通过调控R-loop稳态参与AP的腺泡细胞损伤和炎症反应。鉴于Myl12a作为肌球蛋白轻链家族蛋白和RhoA/ROCK信号轴的下游靶点，其可能通过调控非肌球蛋白II收缩性参与腺泡细胞形态重塑、内皮和肠上皮连接失稳及白细胞跨内皮迁移，进而促进AP相关的肠道屏障功能障碍。
本研究存在以下局限性：生物信息学分析结果主要来源于公共数据库，样本异质性可能引入偏倚；未纳入不同病因和不同严重程度的AP患者进行分层分析；仅采用雨蛙素刺激的AR42J细胞作为体外模型，未在成纤维细胞中进行直接功能验证；三个候选基因特征中仅Myl12a完成了功能验证，Anp32b和Nras的生物学功能及其与R-loop稳态的具体调控关系尚未通过系统的过表达或敲低及临床实验确认；机制推断基于关联分析，未直接测量R-loop水平，也未阐明这些基因特征调控R-loop稳态的具体分子机制。
研究结论： 本研究通过整合转录组数据和单细胞测序数据，系统筛选并鉴定了AP中与R-loop稳态失衡相关的关键细胞群体和候选基因特征，阐明了这些分子通过多维免疫炎症通路协同调控AP发生发展的潜在机制。这些发现增强了对AP免疫炎症网络和转录失调的理解，为AP的早期检测和靶向治疗建立了新的理论框架和分子靶点。
