生信小课堂
影响因子：8.44
关于非肿瘤生信，我们也解读过很多，主要有以下类型
1 单个疾病WGCNA+PPI分析筛选hub基因。
2 单个疾病结合免疫浸润，铁死亡，自噬等基因集，机器学习算法等。
3 两种相关疾病联合分析，包括非肿瘤结合非肿瘤，非肿瘤结合肿瘤或者非肿瘤结合泛癌分析
4 基于分型的非肿瘤生信分析
目前非肿瘤生信发文的门槛较低，有需要的朋友欢迎交流
研究背景：
冠状动脉疾病(CAD)仍是全球主要死亡原因之一。动脉粥样硬化病变通常导致进行性冠脉狭窄、心肌缺血、坏死，最终累积急性心血管综合征(ACS)的风险。外周生物标记物是监测冠状动脉疾病(CAD)进展的越来越重要的非侵入性方法。近年来的证据表明，免疫紊乱在CAD病理过程的启动和发展中起着决定性作用。
研究结果：
一、外周血中差异mRNA表达的鉴别
1、本研究中一共有199例CAD患者和218例对照组的外周血转录组。PCA分析GSE20680，GSE20681 和GSE42148表明，批次效应在三个基因集中已被去除（图2A，B）。
2、共鉴定出762个DEGs，其中上调基因396个，下调基因366个(图2C)。
3、热图显示CAD患者与对照组之间的前17 DEGs，提示其可能参与CAD的病理过程（图2D）。
二、DEGs的功能富集分析
1、蛋白质-蛋白质相互作用(PPI)网络说明了蛋白质水平上的交互关系（图3A）.
2、GO分析显示DEGs涉及的分子功能，生物过程和细胞成分（图3B）。
3、KEGG通路富集分析显示，DEGs集中于5种与免疫相关的信号通路，其中以细胞因子与细胞因子受体相互作用的信号通路最多（图3C）。
三、DE-IRGs测定和免疫相关细胞图谱分析
1、韦恩图（Venn）显示，上述CAD患者的DEGs与IRGs中，交集到58个DE-IRGs（图 4A）。
2、热图显示，主要有28种免疫细胞发生改变（图4B），箱式图显示CAD患者和对照组间，免疫细胞富集分数的差异（图4C）。其中，活化的树突状细胞、肥大细胞、中性粒细胞、Th17细胞和MDSC在CAD中的比例增加，活化的B细胞和CD56细胞的比例减少。
3、气泡图显示，发生改变的28中免疫细胞存在相关性（图4D）。例如，活化的CD4+ T细胞和Th2细胞具有相关性（r=0.55），中性粒细胞和血浆树突细胞具有相关性（r=0.62）。
四、CAD中14个关键外周DE-IRGs的鉴定及免疫诊断模型的构建
1、通过LASSO回归分析，最终确定14个基因，构建CAD诊断模型（图5A、B）。
2、利用ROC曲线评价14个基因的诊断效果（图5C、D），其中14个基因ROC曲线下面积：CCR9 (AUC=0.686), CER1 (AUC=0.604), CSF2 (AUC=0.688), CXCL2 (AUC=0.589), HTR3C (AUC=0.699), IL13RA1 (AUC=0.650), IL1A (AUC=0.667), INSL5 (AUC=0.624), MBL2 (AUC=0.619), MMP9 (AUC=0.649), MSR1(AUC=0.679), NR4RA2(AUC = 0.657), NTS(AUC = 0.667) and TNFRSF19 (AUC=0.686)。
3、ROC曲线分析上述14个基因的诊断模型训练集（图5E、F），当14个基因都包含在一个诊断模型中时，训练数据集和测试数据集的AUC分别为0.968和0.859。
五、CAD患者外周血免疫特性分析
1、箱式图显示，外周血中所鉴定的mRNA中，14个关键外周DE-IRGs在CAD组与对照组之间有显著差异（图6A）。其中, CCR9, CER1, CSF2, IL13RA1, INSL5, MBL2, MMP9, MSR1, NTS 和TNFRSF19基因上调表达, 而CXCL2, HTR3C, IL1A,和NR4A2基因表达下调。
2、气泡图显示，这些基因大多数具有正相关的关系（图6B）。
3、多种免疫成分在CAD诊断和病理机制中的重要作用，热图显示了14个CAD中DE-IRGs和28种免疫细胞相关性，以此分析免疫成分与冠心病之间的相互关系（图6C）。
六、在体内验证10个候选DE - IRGs与CAD的相关性
1、小鼠主动脉弓油红O染色显示，20周小鼠主动脉弓内大量脂质斑块聚集，导致管腔狭窄（图7）。
2、对于体重发现，CAD模型组平均体重明显高于对照组（图7B）。
3、CAD组外周循环低密度脂蛋白(LDL)和甘油三酯(TRIG)水平升高（图7C、D）。
4、全血总mRNA检测结果显示，CAD模型组CCR9、CSF2、IL13RA1、NTS表达明显高于对照组（图7E-H、M）。
5、但是IL1A，在模型组中也有升高，与预测数据结果相反（图7I）。
6、然而，CXCL2、INSL5、MSR1、NR4A2和TNFRSF19在CAD模型组和对照组间无差异 (图7G，J-L，N)。
七、IL13RA1激活JAK1/STAT3通路，调节CAD中巨噬细胞功能
1、细胞实验中，ox-LDL处理使BMDMs中IL13RA1 mRNA的表达增加，并且激活了JAK1/STAT3信号通路，而且使VEGF-C的产生增加（图8A，B-E）。
2、动物实验中，在HFD饮食诱导的20周龄LDL缺失小鼠(CAD小鼠模型)的主动脉中，IL13RA1 mRNA表达、p-JAK1、p-STAT3、VEGF-C和IL-6蛋白水平升高(图8C，D-F)。
3、免疫荧光图可见转染siRNA的RAW264.7细胞中有良好的转染效果(图8G)， IL13RA1 mRNA表达被显著减少 (图8H)，并且磷酸化的JAK1和STAT3、VEGFC、TGF-β和αSMA的表达水平降低 (图8I，J)。
4、IL13RA1敲除组中SCARB1和ox-LDL处理组的CD36 mRNA表达水平降低（图8M、N）。然而，IL13RA1敲除却增加了IL-6的表达量 (图8I，J)。
5、油红O染色结果显示IL13RA1敲除组巨噬细胞吸收的脂滴减少（图8K，L）。
研究总结：
在该研究中，作者利用数据库中筛选的14个关键性的差异表达免疫相关基因(DE -IRGs)的基础上，构建了晚期CAD诊断模型。然后，作者使用低密度蛋白受体缺陷小鼠(典型的CAD动物模型)来确定候选基因与CAD的相关性，其中IL13RA1可能通过调节JAK1/STAT3通路参与动脉粥样硬化的炎症、纤维化和胆固醇分泌过程。
