《Journal of Translational Autoimmunity》：Systematic dysregulation of immune-related alternative polyadenylation in systemic lupus erythematosus contributes to patient stratification
编辑推荐：
背景：选择性多聚腺苷酸化（APA）是一种常见的转录后基因调控机制。尽管APA失调与多种人类疾病相关，但其在系统性红斑狼疮（SLE）中的作用仍 largely 未被探索。研究人员系统地表征了SLE中免疫相关APA变化，并研究了其用于临床有用的患者分层的潜力。方法
背景：选择性多聚腺苷酸化（APA）是一种常见的转录后基因调控机制。尽管APA失调与多种人类疾病相关，但其在系统性红斑狼疮（SLE）中的作用仍 largely 未被探索。研究人员系统地表征了SLE中免疫相关APA变化，并研究了其用于临床有用的患者分层的潜力。方法：研究人员整合了来自多个SLE队列的RNA-seq数据（涉及外周血样本），并使用DaPars2量化APA。研究人员开发了一个评分框架，以识别关键的免疫相关APA事件（ImmAPA）及其在SLE患者分组中的效用。研究人员还使用Connectivity Map（CMap）探索了跨ImmAPA组药物响应的潜在差异，并在一个接受免疫抑制药物治疗的独立SLE队列中检验了这些关联。结果：SLE患者在各队列中表现出3′非翻译区（3′UTR）延长的复发趋势。二十个关键的APA调控因子在SLE患者中与这些变化的关联增强。ImmAPA评分可将患者分为ImmAPA-high和ImmAPA-low亚型，反映激活状态，并与疾病严重程度相关。与干扰素（interferon, IFN）特征相比，ImmAPA评分捕获了互补的转录组变异。CMap分析预测了可逆转这些APA改变的药物，并得到后续分析的支持。结论：研究人员的研究揭示了SLE中外周血转录组水平上向3′UTR延长的复发广泛转变。ImmAPA框架为SLE患者分层提供了新视角，并为SLE的精准治疗策略提供了见解。
**论文解读文章**
**研究背景**
系统性红斑狼疮（systemic lupus erythematosus, SLE）是一种影响多器官系统的复杂自身免疫性疾病，临床异质性显著，治疗疗效差异大。现有研究主要聚焦于基因表达水平和干扰素（interferon, IFN）为中心的转录程序，这些程序虽可用于患者分层，但未能充分反映SLE中免疫状态的多样性，且对转录后调控层次的系统性贡献关注不足。选择性多聚腺苷酸化（alternative polyadenylation, APA）是mRNA代谢的关键步骤，通过选择不同的多聚腺苷酸化位点（polyadenylation site, PAS）产生不同3′非翻译区（3′ untranslated region, 3′UTR）长度的异构体，影响mRNA稳定性、翻译效率和定位。APA失调与多种疾病相关，但在SLE中的研究多局限于个别基因，缺乏系统性分析。因此，有必要从APA调控角度揭示SLE的分子异质性，开发新的分层方法以改善精准诊断和治疗。
**研究内容与结论**
研究人员整合了多个SLE队列（GSE72509、GSE110685、GSE122459、GSE163073及B细胞数据集GSE175913）的外周血RNA-seq数据，使用DaPars2量化APA事件，并开发了ImmPathAPAScore框架系统识别免疫相关APA事件（ImmAPA）。研究发现SLE患者普遍存在3′UTR延长趋势，20个关键APA调控因子与这些变化关联增强，且遗传预测的部分调控因子表达与SLE风险相关。基于ImmAPA模块构建的ImmAPA评分可稳定地将患者分为ImmAPA-high和ImmAPA-low亚型，该评分与疾病活动指数（SLEDAI）呈正相关，并提供了超越传统IFN评分的互补信息。通过Connectivity Map（CMap）分析，研究人员预测了可逆转APA变化的候选药物，并在外部霉酚酸酯（mycophenolate mofetil, MMF）治疗队列中初步观察到亚型间的临床差异。该研究发表在《Journal of Translational Autoimmunity》。
**主要技术方法**
1. **数据资源**：从GEO获取5个SLE队列（全血样本：GSE72509、GSE110685；外周血单个核细胞PBMC：GSE122459、GSE163073；B细胞亚群：GSE175913）。
2. **APA定量**：使用DaPars2计算远端多聚腺苷酸化位点使用指数（percentage of distal poly(A) site usage index, PDUI），鉴定差异APA事件。
3. **ImmPathAPAScore框架**：将PDUI与ImmPort免疫通路基因集整合，通过单样本基因集富集分析（ssGSEA）量化APA事件与免疫通路的耦合，筛选可信的ImmAPA（202个事件）。
4. **患者分层**：基于ImmAPA事件进行无监督聚类，计算ImmAPA评分，划分ImmAPA-high和ImmAPA-low亚型。
5. **药物重定位分析**：利用CMap基于亚型特异性差异表达基因（DEG）进行药物转录连接分析，并在外部MMF治疗队列（GSE49454）验证。
**研究结果**
**3.1 系统性APA改变**
通过比较SLE患者与健康对照（HC）的PDUI，发现SLE患者在多个队列中整体PDUI显著升高，表明广泛3′UTR延长。延长事件在至少两个队列中一致重复，比例58.0%–90.7%；B细胞数据集则呈缩短趋势，故排除于跨队列整合。延长相关基因富集于T细胞受体信号、NF-κB信号等SLE相关免疫通路。
**3.2 关键APA调控因子驱动3′UTR延长**
分析83个已知APA调控因子与APA事件的关联，在SLE中关联强度高于HC，且以正相关为主。筛选出20个关键调控因子，孟德尔随机化（Mendelian randomization, MR）分析显示CPSF4、EIF3B、GTF2I、PABPC1的遗传预测表达与SLE风险相关。这些调控因子的靶基因富集于抗病毒天然免疫通路（如RIG-I样受体信号）；IFN响应基因（如IFI44）与多个调控因子关联，且同时呈现上调与3′UTR延长。蛋白-蛋白互作（PPI）网络显示多数调控因子高度互作，提示协同调控。
**3.3 整合免疫通路分析识别ImmAPA**
通过ImmPathAPAScore，研究人员筛选出202个在至少两个队列中重复出现且富集于≥20个免疫通路的APA事件，定义为ImmAPA。不同免疫通路的ImmPathAPAScore模式差异显著，表明APA改变优先耦合特定免疫功能模块。
**3.4 ImmAPA评分定义免疫激活的互补维度并进行患者分层**
ImmAPA评分与大多数免疫通路活性呈一致正相关，并与SLEDAI呈中等正相关（Spearman ρ=0.254, p=0.0158），而IFN评分与SLEDAI相关性不显著。嵌套线性回归显示，加入ImmAPA评分可改善IFN评分对SLEDAI的解释能力。ImmAPA评分呈双峰分布，可稳定分出ImmAPA-high和ImmAPA-low亚型。ImmAPA-high组中白介素-1（Interleukin-1, IL-1）家族信号活性显著高于HC和ImmAPA-low组。
**3.5 ImmAPA亚型与转录组重塑及治疗响应相关**
CIBERSORTx分析显示ImmAPA-high组CD14+单核细胞比例升高。差异表达分析发现ImmAPA-high组有2493个DEGs，ImmAPA-low组仅80个，重叠仅40个基因，表明ImmAPA-high与广泛免疫转录激活相关。CMap分析预测不同药物对两亚型转录扰动呈现不同连接分数，候选作用机制（mechanism of action, MOA）集中于脱氢酶抑制剂、JAK抑制剂等。在外部MMF治疗队列中，ImmAPA-high患者治疗后SLEDAI低于ImmAPA-low患者，且ImmAPA-high与较低的肾脏受累及低补体（C3、C4）发生率相关。
**讨论与结论**
讨论部分指出，APA重塑是SLE中跨队列的普遍转录组水平模式，3′UTR延长可能通过增加miRNA和RNA结合蛋白（RBP）结合位点影响mRNA稳定性。20个关键APA调控因子形成高度互作网络，提示协同调控远端PAS使用。ImmAPA评分与IFN评分弱相关，提供转录后层次补充信息，但临床验证仍需前瞻性队列。研究局限性包括依赖bulk RNA-seq数据受细胞组成影响、需实验验证代表性ImmAPA事件、队列异质性及缺乏标准化治疗指标等。结论：本研究拓展了对SLE转录后调控的认识，ImmAPA框架从APA调控角度实现患者分层，揭示了SLE异质性延伸至3′端加工和转录后调控层面，为未来机制研究和精准分层提供了可解释框架。
