环状RNA-微RNA-转录因子-信使RNA竞争性内源RNA(ceRNA)调控网络的构建与分析揭示类风湿关节炎相关的核心基因

背景：环状RNA(circular RNAs, circRNAs)、微RNA(microRNAs, miRNAs)和转录因子(transcription factors, TFs)参与驱动类风湿关节炎(rheumatoid arthritis, RA)的免疫失调，但其在外周血中整合的调控架构尚未明确。因此，研究人员构建了一个circRNA-miRNA-TF-mRNA竞争性内源RNA(competing endogenous RNA, ceRNA)网络，并在胶原诱导关节炎(collagen-induced arthritis, CIA)模型中进行实验验证。 方法：研究人员从基因表达综合数据库(Gene Expression Omnibus, GEO)检索公开的RA数据集，包括GSE189338(circRNA)、GSE124373(miRNA)和GSE17755(mRNA)。使用limma包鉴定差异表达转录本。采用基因本体论(Gene Ontology, GO)和京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes, KEGG)分析对差异表达mRNA进行功能注释。TF-mRNA配对关系取自TRRUST v2，经实验验证的miRNA-靶基因互作取自miRTarBase，circRNA-miRNA互作使用CircBank预测。候选基因在比较毒理基因组数据库(Comparative Toxicogenomics Database, CTD)中进行交叉核对。选取网络组分使用CIA大鼠外周血进行实时荧光定量聚合酶链反应(quantitative real-time PCR, RT-qPCR)验证。 结果：研究人员鉴定出920个差异表达mRNAs、282个差异表达miRNAs和980个差异表达circRNAs。功能富集显示mRNA特征与白细胞黏附、凋亡相关过程、灶性黏附及RA相关的免疫炎症通路相关。整合TF、miRNA和circRNA各层得到初始网络含6个circRNAs、4个miRNAs、4个TFs和24个mRNAs。为进行生物学验证，研究人员优先选取6个circRNAs、4个miRNAs、3个TFs和8个mRNAs进行RT-qPCR。大多数候选分子表达变化与发现数据集一致，而miR-195-5p组间无显著差异。剔除不支持分支后，精炼网络由5个circRNAs、3个miRNAs、3个TFs和8个mRNAs组成。 结论：本研究定义了RA中基于血液的circRNA-miRNA-TF-mRNA调控框架，并提名了一组聚焦的候选生物标志物供后续机制与临床验证。研究结果应被视为假设生成性的，而非ceRNA因果调控的决定性证据。

本文解读基于论文《Construction and Analysis of ceRNA Regulatory Networks Reveal the Core Genes Associated with Rheumatoid Arthritis》，发表于《International Journal of General Medicine》。

研究背景：类风湿关节炎（rheumatoid arthritis, RA）是一种慢性系统性自身免疫病，以持续性滑膜炎、进行性软骨破坏和骨侵蚀为特征。尽管治疗有所进展，RA在早期诊断、患者分层和分子监测方面仍存在未满足需求。外周血可重复采样且反映全身免疫活化，是理想的转化生物标志物来源。非编码RNA（non-coding RNAs, ncRNAs）——特别是环状RNA（circular RNAs, circRNAs）和微RNA（microRNAs, miRNAs）——已被报道在RA滑膜组织、成纤维样滑细胞膜细胞（fibroblast-like synoviocytes, FLSs）及外周血单核细胞（peripheral blood mononuclear cells, PBMCs）中异常表达。竞争性内源RNA（competing endogenous RNA, ceRNA）假说认为，含有共同miRNA应答元件的RNA可通过竞争结合miRNA而相互调控。既往RA网络研究多局限于circRNA-miRNA-mRNA三元关系，未纳入转录因子（transcription factor, TF）这一可放大或约束转录后效应的关键调控层。此外，多数ceRNA互作为预测所得，缺乏实验验证。为此，研究人员整合多层级调控信息构建circRNA-miRNA-TF-mRNA网络，并在动物模型中实验验证，以筛选RA相关的候选血液生物标志物及调控关系。

主要关键技术方法：研究人员下载GEO数据库中三个RA外周血/ PBMCs数据集——GSE189338（RA与健康人PBMCs的circRNA芯片，n=4 vs 4）、GSE124373（PBMCs的miRNA芯片，n=28 vs 18）、GSE17755（外周血细胞mRNA芯片，n=112 RA vs 对照）。用limma进行差异表达分析（|log₂FC|>0.2, 名义P<0.05），对差异mRNAs做GO和KEGG富集。TF-mRNA对取自TRRUST v2，miRNA靶互作取自miRTarBase 2022，circRNA-miRNA互作用CircBank预测，取三层交集构建整合网络并用Cytoscape可视化。候选基因经比较毒理基因组数据库（Comparative Toxicogenomics Database, CTD）交叉核对RA相关性后，依网络拓扑中心性及跨层桥接作用优选靶标。选用雄性Sprague-Dawley大鼠建立胶原诱导关节炎（collagen-induced arthritis, CIA）模型，取外周血提取总RNA，以GAPDH为内参（circRNA/mRNA/TF）、U6为内参（miRNA），通过RT-qPCR验证所选分子，2^-ΔΔCt法计算相对表达量，两组间比较用非配对学生t检验。

研究结果：

概述分析流程（Overview of the Analytical Workflow）：对三个GEO数据集预处理后，差异表达分析筛出980个circRNAs、282个miRNAs和920个mRNAs。差异mRNAs进行富集分析和网络构建，部分候选进入CIA模型验证。

mRNA特征的功能学特点（Functional Characteristics of the mRNA Signature）：差异表达mRNAs显著富集于白细胞细胞间黏附（leukocyte cell-cell adhesion）、凋亡信号调控（regulation of apoptotic signaling）、灶性黏附（focal adhesion）及泛素相关结合功能等GO条目；KEGG富集提示炎症与免疫相关通路，符合RA免疫病理中白细胞招募、基质活化及存活信号的经典认识。

整合调控网络的构建（Construction of the Integrated Regulatory Network）：将差异mRNAs与TRRUST交集得46对TF-mRNA（涉及6个TFs、42个靶基因）；通过miRTarBase关联差异miRNAs得含12个miRNAs的miRNA-TF/mRNA网络；再用CircBank预测并与差异circRNAs取交得6个circRNAs锚定于该调控核心。初始网络含6个circRNAs、4个miRNAs（hsa-miR-29b-3p、hsa-miR-424-5p、hsa-miR-195-5p、hsa-miR-4646-5p）、4个TFs（SPI1、HIF1A、MYB、DNMT1等中入选4个）和24个mRNAs。

RT-qPCR靶标的优先选取（Prioritization of RT-qPCR Targets）：因动物实验验证通量有限，依网络位置、跨层连接度及CTD中RA相关性，优先选定6个circRNAs（circ0086684、circ0001605、circ0043947、circ0009723、circ0015911、circ0029987）、4个miRNAs、3个TFs（DNMT1、MYB、HIF1A）和8个mRNAs（TLR2、MIF、HDAC1、ENO1、STAT4、RB1、CXCR4、GATA3）进行RT-qPCR。

CIA模型的验证（Validation of the CIA Model）：CIA组大鼠出现渐进性爪肿胀及升高关节炎评分，对照组稳定；组织学示滑膜增生、炎细胞浸润及软骨/骨损伤，Micro-CT示骨破坏及骨小梁受损，证实CIA炎症性关节炎模型构建成功。

RT-qPCR验证与网络精炼（RT-qPCR Validation and Refinement of the Network）：CIA组外周血中，所选circRNAs、3个TFs及8个mRNAs普遍较对照组上调；miR-4646-5p、miR-29b-3p、miR-424-5p呈预期方向显著改变；miR-195-5p组间无显著差异被剔除。最终精炼网络含5个circRNAs、3个miRNAs（miR-4646-5p、miR-29b-3p、miR-424-5p）、3个TFs（DNMT1、MYB、HIF1A）和8个mRNAs（TLR2、MIF、HDAC1、ENO1、STAT4、RB1、CXCR4、GATA3）。

表达方向模式的解释（Interpretation of the Expression-Direction Pattern）：验证子网络偏向同向差异表达反映其为从异源独立数据集筛选并经动物模型测试的小子集，不代表所有ceRNA关系均功能协同，精炼网络应视为优先调控假说而非已解析的ceRNA因果机制。

讨论总结：本研究核心价值在于从宽泛计算网络缩小至经实验支持的候选分子集。保留的mRNAs（TLR2、MIF、HDAC1、STAT4、HIF1A相关通路分子CXCR4、ENO1、RB1、GATA3）均有文献支持参与RA免疫活化、炎症或滑膜成纤维细胞存活。非编码RNA层中，miR-29b-3p与RA单核细胞凋亡炎症、miR-424在RA滑膜成纤维细胞中异常有前期依据，miR-4646-5p在RA中研究极少仅为候选，miR-195-5p未验证提示网络预测需实验把关。局限含发现集来自不同队列与平台、circRNA样本量小、筛选阈值宽松用于探索、动物血而非人PBMCs致种属及细胞构成差异、缺独立人队列验证及功能缺失/获得实验，故ceRNA关系仍为假定。研究定位为优先排序（prioritization）而非机制确证论文，所提精炼网络及候选分子可供后续人源细胞、大样本临床队列及干预实验检验。

结论（Conclusion）翻译：通过整合公共转录组数据集、转录因子信息及胶原诱导关节炎模型中的实验筛选，研究人员鉴定出一个与类风湿关节炎相关的精炼circRNA-miRNA-TF-mRNA调控网络。保留的分子代表候选外周血生物标志物及可用于后续研究的简明框架。这些发现为未来的人体研究和机制验证提供了聚焦的基础。