有很多同学看见肿瘤方向做生信数据挖掘发了很多SCI文章,就觉得非肿瘤方向发不了。非肿瘤生信热度虽然不如肿瘤生信,也没有它好发表,原因一则数据集少,二则临床信息不够。
但是!非肿瘤生信也有它的优势,相同的研究套路在肿瘤领域迁移到非肿瘤领域,不失为发SCI的“捷径”,或许可以发到更好的分数。今天小编就和大家分享1篇非肿瘤生信文章。
牙周炎是在牙周支持组织中发生的一种慢性炎症,是由对牙菌斑中细菌的炎症反应引起的。此外,牙周组织的炎症是宿主免疫防御反应与微生物相互作用的结果,牙周炎的发病机制与免疫密切相关。因此联想到牙周炎 免疫 m6A。
文章题为:“m6A调节因子介导的RNA甲基化修饰模式参与牙周炎的免疫微环境调控。”,发表于J Cell Mol Med杂志,影响因子5.310。
Journal of cellular and molecular medicine, 25(7), 3634–3645.
https://doi.org/10.1111/jcmm.16469
一、期刊信息
科研猫公众号后台对话框发送期刊名称(支持缩写、简写),即可查询期刊最新影响因子、近5年IF和发文量、期刊预警、分区等~
二.材料和方法
- 从GEO数据库中下在基因表达数据GSE16134,其中包括310例(69个健康样本和241个牙周炎样本);
- 从STRING数据库研究了疾病中23个m6A RNA甲基化调控因子的蛋白相互作用、基于lasso模型筛选特征基因,并利用ROC评估模型;
- ssGSEA评估免疫浸润细胞数量和免疫反应活性, 免疫反应基因集来自 ImmPort 数据库;
- 无监督聚类分析和功能富集分析;
三.结果
1.健康和牙周炎样本之间的m6A调节因子
- 23个m6A调节因子之间存在密切相关性;
- CBLL1和YTHDC2是所有样本和牙周炎样本中表达最相关的m6A调节因子;
- 正常样本和疾病样本中大部分m6A调控因子在转录水平上存在显著差异。
2.m6A调节因子参与牙周炎进程
- 单变量逻辑回归分析揭示15种与牙周炎相关的m6A调节因子;
- 多因素logistic回归分析建立15种m6A RNA甲基化调节因子的识别特征,并计算牙周炎的风险评分;
- 健康和疾病之间的风险分布,其中牙周炎的风险评分远高于健康样本;
- 通过ROC曲线分析m6A调节因子对健康和牙周炎样本的辨别能力,并通过AUC值进行评估
3.m6A调节因子与牙周炎的免疫特征相关
- 对失调的m6A调剂因子与免疫浸润性细胞、免疫反应基因集、HLA基因集,进行相关性分析
- A.免疫细胞层面:最正相关:ALKBH5-单核细胞,小提琴图表明牙周炎中ALKBH5的表达较高,单核细胞较多;最负相关:FMR1-单核细胞,小提琴图表明牙周炎中FMR1的表达较低。
- B.免疫反应层面:最正相关:ELAVL1-TNF家族成员受体,小提琴图显示,表明在牙周炎中ELAVL1的表达更高,TNF家族成员受体反应更活跃。最负相关:CBLL1-细胞因子,小提琴图显示,表明在牙周炎中CBLL1表达较高,细胞因子反应较不活跃。
4.23种调节因子介导的m6A RNA甲基化修饰模式在牙周炎中的作用
- 基于23种m6A调节因子的表达对牙周炎样本进行无监督一致性聚类分析;
- 根据m6A调节因子的表达水平,患者被分为3个亚型;
- 除 YTHDC1 外,其余 m6A 调节因子在 3 个 m6A 修饰模式之间的表达存在显着差异;
5.不同m6A修饰模式免疫微环境特征差异
- 免疫浸润性细胞、免疫反应基因集、HLA家族基因表达在在三种m6A修饰模式中存在差异;
- 亚型-2和亚型-3 m6A修饰介导主动免疫反应,pattern-1在牙周炎中导致轻度免疫反应,而pattern-2和pattern-3介导的免疫反应不同。再次证明m6A甲基化修饰对牙周炎不同免疫微环境的塑造具有重要的调控作用。
6.3个m6A修饰模式之间潜在的生物学功能特征多样性
- GSVA富集分析三种m6A模式亚型HALLMARKS通路和KEGG通路
- ①图AB:亚型-1和亚型-2;②图CD:亚型-1和亚型-3;③图EF:亚型-2和亚型-3;
7.牙周炎m6A表型相关基因的鉴定及功能分析
- 确定三种m6A修饰模式中共有的DEGs,1631个常见基因被视为m6A表型相关基因;
- GO-BP功能富集分析:参与免疫调剂的m6A表型相关基因:①主要参与细胞周期和蛋白质修饰的调节;②与抗原过程、Fc信号转导、免疫细胞增殖调控和细胞因子介导的通路显着相关
- 构建 m6A 修饰模式相关的综合基因图谱,WGCNA 鉴定不同 m6A 修饰相关的基因-基因模块,WGCNA采用前25%的变异基因进行分析,排除异常样本;
- G图绿色模块的hub基因与m6A修饰模式-1高度相关。
ps:本文进行了三种亚型分析,因此在进行数据收集的时候注意样本量哦。