多组学分析揭示溃疡性结肠炎缓解期脾虚湿困证与活动期湿热阻滞证间的潜在分子机制及靶向中药预测

目的 综合分析溃疡性结肠炎（ulcerative colitis,UC）相关scRNA-seq和RNA-seq数据集,探讨UC缓解期脾虚湿困证与活动期湿热阻滞证间的潜在分子机制,并挖掘潜在干预中药。方法 运用差异基因表达分析和韦恩图鉴定出UC缓解期脾虚湿困证和活动期湿热阻滞证相关基因及两者间的对话基因。通过京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）、基因集变异分析（gene set variation analysis,GSVA）及Spearman相关性分析深入探讨对话基因参与的生物过程及潜在功能。通过单细胞转录组分析探讨UC缓解期与活动期的差异情况及对话基因在其中的作用。基于循环算法构建UC缓解期脾虚湿困证向活动期湿热阻滞证发生发展的风险预测模型。运用外部数据集、动物实验与实时定量聚合酶链反应（real-time quantitative polymerase chain reaction,RT-PCR）进一步验证对话基因的表达情况。借助SoFDA数据库构建“基因-中医症状/现代医学症状”网络;借助TCMIP和COREMINE数据库预测对话基因的潜在靶向中药。结果 鉴定出31个UC缓解期脾虚湿困证相关基因,其主要富集在代谢相关途径;鉴定出160个UC活动期湿热阻滞证相关基因,其主要富集在炎症免疫相关途径。鉴定出22个UC缓解期脾虚湿困证与活动期湿热阻滞证之间的对话基因,上调对话基因与炎症免疫相关途径呈显著相关性,下调对话基因与代谢途径呈显著相关性。单细胞转录组分析显示,相比于正常对照及UC缓解期,UC活动期中B细胞占比上调;B细胞在UC活动期湿热阻滞证中作用贡献较大,而在UC缓解期脾虚湿困证中的作用贡献较小;UC活动期B细胞的受体信号通路活性明显高于UC缓解期;上调的对话基因CD55在B细胞中表达水平较高,CD55⁺B细胞与CD55⁻B细胞之间存在差异的细胞通讯及细胞代谢水平。6个风险预测基因（CD55、PCSK1、SERPING1、ACAT1、SLC26A2、HMGCS2）可以较为准确预测UC缓解期脾虚湿困证向UC活动期湿热阻滞证发展,动物实验及RT-PCR验证了其表达水平。通过数据库检索得到人参、黄芪、枸杞子等45味潜在靶向中药。结论 鉴定出22个关键对话基因,它们可能通过影响炎症免疫和代谢相关途径诱导UC缓解期脾虚湿困证向UC活动期湿热阻滞证发展。B细胞在UC缓解期脾虚湿困证向活动期湿热阻滞证发生发展过程中发挥重要作用。CD55可能通过影响炎症免疫和代谢相关途径,促进B细胞活性,从而促进UC缓解期脾虚湿困证向活动期湿热阻滞证发展。22个关键对话基因与UC相关中医症状和现代医学症状关联密切。

Revelation of potential molecular mechanism between spleen deficiency and dampness-stagnation syndrome in remission stage and damp-heat stagnation syndrome in active stage of ulcerative colitis based on multiomics analysis and prediction of targeted traditional Chinese medicine