4+非肿瘤纯生信。氧化应激+WGCNA+药物预测筛序关键基因



今天给同学们分享一篇非肿瘤+氧化应激+WGCNA的生信文章“Identification of oxidative stress-related biomarkers associated with the development of acute-on-chronic liver failure using bioinformatics”,这篇文章于2023年10月10日发表在Scientific Reports期刊上,影响因子为4.6。

0ab25fcbb8ea5958c314184cb2d04e4f.jpeg

急性慢性肝衰竭(ACLF)是一种独特的综合征,其特征是慢性肝病急性失代偿、多器官功能衰竭和高短期死亡率。ACLF的病理生理机制主要包括强烈的全身性炎症、免疫功能障碍、线粒体功能障碍、代谢变化和氧化应激。迄今为止,除了肝移植外,对于ACLF还没有明确有效的治疗策略。ACLF管理的主要原则是治疗诱发事件或相关并发症,并提供器官支持治疗。迫切需要进一步探索ACLF发展和进展的病因学,以制定更有效的治疗策略,降低该疾病的高死亡率。

fc9fc520b684b09bb9c513b33c71e9d1.jpeg

 图1 本研究的流程图


1. 鉴定差异表达基因(DEGs)和差异表达有机物基因(DEOSGs)

从数据集GSE142255中筛选出了1228个不同表达的基因(DEGs),包括432个上调基因和796个下调基因(图2A、B)。为了进一步确定这些DEGs的来源和通路,作者使用BTMs作为参考基因集进行了GSEA分析。结果显示,上调的DEGs在与TLR和炎症信号传导、单核细胞、细胞周期和转录、单核细胞表面标志、免疫激活通用簇和细胞外基质相关的细胞系或通路中显著富集,而下调的DEGs在与T细胞激活、信号传导和分化、T细胞表面标志、NK细胞和MHC-TLR7-TLR8簇相关的细胞系或通路中显著富集(图2C)。在1228个DEGs和1093个氧化应激相关基因(OSGs)之间进行交集分析后,得到了79个差异表达的氧化应激相关基因(DEOSGs)(图2D)。

89f0ed3f5b1e45d6e6191bc3d5878ca0.jpeg

图2 ACLF中的DEGs和DEOSGs


2. DEOSGS的功能富集分析

GO富集分析显示,上调的差异表达氧化应激基因(DEOSGS)在与炎症反应调节、对肽的反应、对细菌源分子的反应、对脂多糖的反应、对氧化应激的反应等生物过程中显著富集(图3A),而下调的DEOSGS在与肽酰丝氨酸修饰、肽酰丝氨酸磷酸化、炎症反应调节、金属离子转运调节和对氧化应激的反应等生物过程中显著富集(图3B)。Reactome富集通路分析显示,上调的DEOSGS在与先天免疫反应、中性粒细胞颗粒溶解、白细胞介素信号传导、细胞因子信号传导和白细胞介素4/13信号传导等信号通路中显著富集(图3C),而下调的DEOSGS在与细胞因子信号传导、白细胞介素信号传导、白细胞介素4/13信号传导、细胞对热应激的反应和白细胞介素12刺激后的JAK STAT信号传导等信号通路中显著富集(图3D)。

697b182cc1ba6cd2c9e0b850fe96fbae.jpeg

图3 不同表达的氧化应激基因(DEOSGS)的功能富集分析


3. 不同表达的氧化应激基因(DEOSGS)的功能富集分析

在基因集GSE142255中,使用氧化应激相关基因(OSGs)的表达信息进行基因集富集分析(GSEA)。使用c5.go.bp.v2023.1.Hs.symbols.gmt作为基因集的GSEA显示,ACLF与与细菌防御反应、真菌防御反应、创伤炎症反应、线粒体基因表达、线粒体翻译和肽代谢过程相关的生物过程显著相关(图4A)。使用c2.cp.reactome.v2023.1.Hs.symbols.gmt作为基因集的GSEA显示,ACLF与与先天免疫系统、中性粒细胞颗粒溶解、抗微生物肽、DNA修复和翻译相关的途径显著相关(图4B)。

1fd39bcfbecb832e5ac3fa9b247a69fa.jpeg

图4 基因集富集分析(GSEA)


4. 免疫浸润分析

如功能富集分析所示,DEOSGs和OSGs在生物过程和免疫应答相关的通路中得到了显著富集。因此,作者进一步使用CIBERSORT算法研究了GSE142255中基因表达与免疫细胞浸润之间的关系,并将结果可视化为热图(图5A、B)。在ACLF患者和健康受试者样本中,六种免疫细胞的丰度存在显著差异(P < 0.05),包括B细胞记忆型、CD4记忆型静止T细胞、CD4记忆型活化T细胞、CD8 T细胞、M0巨噬细胞和中性粒细胞(图5C)。在这六种不同浸润的免疫细胞中,中性粒细胞与CD8 T细胞呈显著负相关(P < 0.05)(图5D)。

a9b30e9ea7a7e5ed43fd7b78b4975551.jpeg

图5 免疫浸润分析


5.&nbsp;加权基因共表达网络分析(WGCNA)

为了确定与ACLF高度相关的免疫相关基因,作者进行了WGCNA分析。作者发现了一个异常样本并将其删除(图6A)。根据尺度独立性和平均连接性,作者将β = 10设置为软阈值(图6B)。基于这个指数,生成了15个基因共表达模块(GCMs)(图6C)。其中,棕色、蓝色、青色、棕黄色和绿宝石模块与T细胞CD8、CD4记忆静止细胞和中性粒细胞呈高度相关(相关系数> 0.7,P < 0.001)(图6D)。因此,这些模块被视为后续分析的关键模块。随后,将包含在关键模块中的3668个基因命名为免疫相关模块基因。功能富集分析表明,这些模块基因主要富集在与免疫反应、激活或调节相关的通路中(图6E,F)。

ee92572822e192e3ffb2dccd84f09747.jpeg

图6 基于加权基因共表达网络分析基因集


6.&nbsp;构建PPI网络和筛选中心基因

通过79个差异表达的氧化应激相关基因(DEOSGs)和3668个免疫相关模块基因的交集,共鉴定出67个免疫相关氧化应激差异表达基因(OSDEGS)(图7A)。将DEOSGs上传至STRING数据库后,构建了蛋白质相互作用(PPI)网络,并使用Cytoscape软件进行可视化(图7B)。在PPI网络中,通过Cytoscape的MCODE插件鉴定出最密集连接的簇(18个节点和117条边),即簇1(图7C)。通过MCC算法在簇1中筛选出的前10个中心基因包括C-C趋化因子4(CCL4)、C-C趋化因子5(CCL5)、整合素α-M(ITGAM)、干扰素γ(IFNG)、髓过氧化物酶(MPO)、基质金属蛋白酶-9(MMP9)、NF-κB抑制因子α(NFKBIA)、Toll样受体2(TLR2)、Toll样受体4(TLR4)和组织抑制金属蛋白酶-1(TIMP1)(图7D)。

3d9b31282344e99428a2a345f600984b.jpeg

图7 蛋白质-蛋白质相互作用(PPI)网络和关键基因


7.&nbsp;关键基因的鉴定

为了确定ACL肝衰竭的诊断和预后价值,作者进行了受试者工作特征曲线(ROC)分析。在GSE142255数据集中,所有前10个中心基因的曲线下面积(AUC)值均大于0.7(图8A-C)。在GSE180014数据集中,AUC大于0.7的中心基因包括IFNG、MMP9、MPO、CCL5、ITGAM、TLR2、TLR4和TIMP1(图8D-F)。在GSE139602数据集中,AUC大于0.7的中心基因包括IFNG、MPO、CCL5、ITGAM、NFKBIA、TLR2、TLR4和TIMP1(图8G-I)。在GSE168049数据集中,AUC大于0.7的中心基因包括MPO、CCL5、MMP9、ITGAM、TLR2、TLR4和TIMP1(图8J-L)。总之,共有6个中心基因(MPO、CCL5、ITGAM、TLR2、TLR4和TIMP1)在四个数据集中一致表现出良好的诊断或预后价值,并因此被命名为关键基因。

fa907c77058744943579b337171387f8.jpeg

图8


8.&nbsp;关键基因的基因表达水平

为了确定ACLF和HC之间关键基因表达的差异,作者对这两个组进行了独立样本t检验。在数据集GSE142255和GSE180014中,与HC组的基因表达水平相比,ACLF组的CCL5显著下调,而其他5个关键基因显著上调(P < 0.001)(图9A,B)。在数据集GSE139602和GSE168049中,发现3个关键基因(CCL5,MPO和TLR2)没有显著调节(图9C,D)。

2cbd61a8515333fbb7094b451164419a.jpeg

图9


9.&nbsp;关键基因的功能富集分析

GO富集分析显示,关键基因主要富集在与对细菌来源分子的响应、对脂多糖的响应、髓系白细胞激活、巨噬细胞激活、脂多糖介导的信号通路、细胞对干扰素-γ的响应、基质金属蛋白酶分泌的调节以及对脂肽多糖酸的细胞响应等生物过程相关的领域(图10A)。Reactome富集分析显示,关键基因主要富集在与先天免疫系统、Toll样受体级联、白细胞介素信号传导、IRAK4缺陷TLR2/4和TLR内源性配体调节、免疫系统疾病、白细胞介素10信号传导、Toll样受体TLR1/TLR2级联以及抗原处理交叉呈递等通路相关的领域(图10B)。

0fabe6f998a5dd67d51f3cca12504ed6.jpeg

图10 关键基因的功能富集分析


10.&nbsp;关键基因的基因集富集分析

为了探索关键基因的潜在功能,作者对单基因进行了基因集富集分析(GSEA),参考基因集为“c2.cp.kegg.v7.4.symbols.gmt”。结果显示,MPO、CCL5、ITGAM、TLR2、TLR4和TIMP1在高表达的队列中的基因在免疫和代谢相关途径中高度富集,例如抗原处理和呈递、自身免疫性甲状腺疾病、肠道免疫网络IGA产生、原发性免疫缺陷、鞘脂代谢、果糖和甘露醇代谢、半胱氨酸和甲硫氨酸代谢、丁酸代谢、嘌呤代谢、糖胺聚糖生物合成硫酸肝素、PPAR信号通路、糖胺聚糖降解、糖脂生物合成乳糖和新乳糖系列、泛酸和辅酶A生物合成以及氨酰核糖酸生物合成(图11A-F),表明这些关键基因可能参与免疫应答和生物代谢过程。

bc86b466b2af700585b6eacebbe9dcb7.jpeg

图11 TLR2(A),TLR4(B),CCL5(C),ITGAM(D),MPO(E)和TIMP1(F)的基因集富集分析(GSEA)


11.&nbsp;关键基因与免疫细胞的相关性

为了确定关键基因与免疫细胞浸润的相关性,作者进行了Pearson相关分析。结果显示,CCL5与其他5个关键基因呈负相关,而这些关键基因之间呈正相关(图12A);除了CCL5(正相关)和MPO(无显著性相关)外,所有关键基因与CD8 T细胞呈显著负相关;MPO与TIMP1分别与静止记忆CD4 T细胞和活化记忆CD4 T细胞呈显著负相关和正相关。CCL5与TLR4分别与中性粒细胞呈负相关和正相关(图12B)。

54ef0c354594e019ba7d3e656983722d.jpeg

图12 关键基因之间的相关性


12.&nbsp;识别调控关键基因的转录因子和miRNA

为了确定miRNA和转录因子对关键基因在转录水平上的调控影响,作者获取了miRNA和转录因子与关键基因的相互作用网络,并使用Cytoscape软件进行可视化。总共有23个、22个、19个、7个、4个和3个miRNA与TLR4、TIMP1、CCL5、TLR2、MPO和ITGAM相互作用(图13A)。总共有8个、7个和5个转录因子与ITGAM、TLR2/TLR4/TIMP1和MPO/CCL5相互作用(图13B)。两个miRNA(hsa-mir-129-2-3p和hsa-mir-210-3p)与3个关键基因相互作用,两个转录因子(FOXC1和HNF4A)与4个关键基因相互作用,表明这些miRNA和转录因子可能与关键基因之间有更密切的相互作用。

07174323dbebfe66e44fda3fae499cf1.jpeg

图13 关键基因的调控网络


13.&nbsp;潜在有效药物的鉴定

为了探索针对关键基因的潜在有效药物,作者从药物-基因相互作用数据库(DGIdb)获取了基因-化学物质相互作用网络,并使用Cytoscape软件进行可视化(图14)。在6个关键基因中,TLR2、ITGAM、TLR4和MPO分别与4个、10个、11个和19个药物/化学物质有交集。在41种药物中,芬太尼、二甲基亚砜和茶碱与两个关键基因(MPO和ITGAM)有交集。

f4a0604261009782412cb05b40a32e57.jpeg

图14 药物-关键基因相互作用网络


总结

总之,本研究利用综合生物信息学方法探讨了氧化应激在ACLF发生和进展中的作用。作者确定了与ACLF的诊断和预后密切相关的六个免疫和氧化应激相关基因。此外,作者构建了一个关键基因的调控网络,并探索了一些针对这些关键基因的药物。本研究加深了作者对氧化应激对ACLF发病机制和预后的影响的理解,并为有效的分子靶向治疗策略提供了更多的见解。作者仍需要进行更多的前瞻性研究和进一步的实验证实作者的结果。




本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/121277.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

数据库数据恢复—Oracle数据库报错ORA-01110错误的数据恢复案例

Oracle数据库故障&#xff1a; 北京某公司一台运行oracle数据库的服务器&#xff0c;机房意外断电导致该服务器重启&#xff0c;重启后发现oracle数据库报错。该Oracle数据库没有备份。 Oracle数据库数据恢复过程&#xff1a; 1、北亚企安数据恢复工程师检查该oracle数据库的数…

input改造文件上传,el-table的改造,点击上传,拖拽上传,多选上传

第一个input标签效果 第二个input标签的效果 el-table的改造效果 <template><div class"outerBox"><div class"analyze" v-if"status"><div class"unFile"><div class"mainBox"><img clas…

黔院长 | 黄帝内经:人有四经十二从!

"人有四经十二从"这句话出自《黄帝内经素问》&#xff0c;“四经”指的是与四时相应的正常脉象&#xff0c;也是指四个主要经络&#xff1a;太阳经、少阳经、太阴经和少阴经。在中医理论当中这些经络被认为是人体气血运行的通道。 而“十二从”则表示人体的十二个经脉…

计算机毕设 基于CNN实现谣言检测 - python 深度学习 机器学习

文章目录 1 前言1.1 背景 2 数据集3 实现过程4 CNN网络实现5 模型训练部分6 模型评估7 预测结果8 最后 1 前言 Hi&#xff0c;大家好&#xff0c;这里是丹成学长&#xff0c;今天向大家介绍 一个深度学习项目 基于CNN实现谣言检测 1.1 背景 社交媒体的发展在加速信息传播的…

ChatGLM3 本地部署的解决方案

大家好,我是herosunly。985院校硕士毕业,现担任算法研究员一职,热衷于机器学习算法研究与应用。曾获得阿里云天池比赛第一名,CCF比赛第二名,科大讯飞比赛第三名。拥有多项发明专利。对机器学习和深度学习拥有自己独到的见解。曾经辅导过若干个非计算机专业的学生进入到算法…

Java 将list集合的字符串格式转为Map

Java 将list集合的字符串格式转为Map List<Object> list new ArrayList<>(); Map<String,String> map1 new HashMap<>(); map1.put("fileName","测试1"); map1.put("level","1"); list.add(map1);Map<S…

[Java]前中后序遍历二叉树/递归与非递归

一、递归方法 首先&#xff0c;树形结构都是由递归方式定义的。那么递归是怎么用的&#xff1f; 1、终止条件&#xff1b;2、调用自身 分析 1、什么时候停止&#xff1f; 当结点值为空的时候&#xff0c;返回null&#xff1b; 2、如何调用自身&#xff1f; 以前序遍历为例&…

php之 角色的权限管理(RBAC)详解

RBAC&#xff08;Role-based access control&#xff09;是一种常见的权限管理模型&#xff0c;通过将用户分配至特定的角色&#xff0c;以及为角色分配访问权限&#xff0c;实现了权限管理的目的。以下是关于RBAC的详细解释&#xff1a; 角色&#xff1a;RBAC模型的核心是角色…

Ubuntu编译 PCL 1.13.1 详细流程

Ubuntu编译 PCL 1.13. 详细流程 一、编译环境二、虚拟机准备1. 虚拟机扩容2. 配置交换分区 三、Cmake - gui 生成 MakeFile1. 解决 flann 依赖问题2. 配置 Cmake 四、编译安装1.编译&#xff1a;2. 安装 一、编译环境 Ubuntu&#xff1a;Ubuntu 20.04 VMware&#xff1a;VMwar…

如何学好C++?学习C和C++的技巧是什么?

如何学好C?学习C和C的技巧是什么&#xff1f; 你这三个问题&#xff0c;前两个都是意思是差不多的&#xff0c;那么怎么怎么学习C/C我来问答一下&#xff1a;最近很多小伙伴找我&#xff0c;说想要一些C资料&#xff0c;然后我根据自己从业十年经验&#xff0c;熬夜搞了几个通…

【嵌入式开源库】timeslice的使用,完全解耦的时间片轮询框架构

完全解耦的时间片轮询框架构 简介项目代码timeslice.htimeslice.clist.hlist.c 创建工程移植代码实验函数说明timeslice_task_inittimeslice_task_addtimeslice_tak_deltimeslice_get_task_num 结尾 简介 timeslice是一个时间片轮询框架&#xff0c;他是一个完全解耦的时间片轮…

三十七、【进阶】验证索引的效率

1、准备工作&#xff1a; 创建一张表&#xff0c;该表中有一千万条数据&#xff0c;名为tb_sku&#xff1b; 2、使用主键查询&#xff1a; select * from tb_stu where id1\G; 3、使用非索引查询&#xff1a; 4、给sn字段创建索引&#xff1a; 在创建过程中&#xff0c;发现…

如何使用 nvm-windows 这个工具来管理你电脑上的Node.js版本

nvm-windows 是一个用于管理在 Windows 上安装的多个 Node.js 版本的工具。以下是安装和使用 nvm-windows 的步骤&#xff1a; 第1步&#xff1a;下载 nvm-windows 访问 nvm-windows 的 GitHub发布页面.下载最新版本的 nvm-setup.zip 文件。 第2步&#xff1a;安装 nvm-wind…

conda虚拟环境笔记收录

1、安装conda 增加执行权限&#xff1a; chmod x Anaconda3-2023.03-1-Linux-x86_64.sh 开始执行&#xff1a;./Anaconda3-2023.03-1-Linux-x86_64.sh2、查看版本 conda --version3、查看当前虚拟环境 虚拟环境和全局环境有前缀可见 如果不进行设置&#xff0c;重新启动就变成…

MySQL的基础(一)

MySQL的基础&#xff08;一&#xff09; SQLSQL的语法特点主要包括以下几点&#xff1a;一、 SQL - DDL -- 数据定义语言1.1 数据库操作1.1 显示现有的数据库1.2 创建数据库1.3 删除数据库1.4 使用 1.2 数据表操作1.2.1 表查询1.2.2 表创建1.2.3 修改表 1.2.4 小结 二、SQL - D…

calcite 校验层总结

1、校验的作用 1&#xff09;完善语义信息 例如在SQL语句中&#xff0c;如果碰到select * 这样的指令&#xff0c;在SQL的语义当中&#xff0c;“*” 指的是取出对应数据源中所有字段的信息&#xff0c;因此就需要根据元数据信息来展开。 2&#xff09;结合元数据信息来纠偏…

特殊类设计[下] --- 单例模式

文章目录 5.只能创建一个对象的类5.1设计模式[2.5 万字详解&#xff1a;23 种设计模式](https://zhuanlan.zhihu.com/p/433152245)5.2单例模式1.饿汉模式1.懒汉模式 6.饿汉模式7.懒汉模式7.1饿汉模式优缺点:7.2懒汉模式1.线程安全问题2.单例对象的析构问题 8.整体代码9.C11后可…

UVa129 Krypton Factor(困难的串)

1、题目 2、题意 如果一个字符串包含两个相邻的重复子串&#xff0c;则称它是“容易的串”&#xff0c;其他串称为“困难的串”。例如&#xff0c;BB、ABCDACABCAB、ABCDABCD都是容易的的串&#xff0c;而D、DC、ABDAB、CBABCBA 都是困难的串。 输入正整数 k k k 和 L L L&a…

Linux虚拟机部署与发布项目(Windows版本)

目录 前言 一、虚拟机部署项目的流程 二、单机项目 1. 本机测试 2.虚拟机部署项目 三、前后端分离项目 前言 在软件开发过程中&#xff0c;部署和发布项目是非常重要的一环。使用虚拟机技术可以方便、灵活且可重复使用地部署和发布项目。本篇博客将介绍如何在 Windows 环…

【密评】商用密码应用安全性评估从业人员考核题库(十四)

商用密码应用安全性评估从业人员考核题库&#xff08;十四&#xff09; 国密局给的参考题库5000道只是基础题&#xff0c;后续更新完5000还会继续更其他高质量题库&#xff0c;持续学习&#xff0c;共同进步。 3251 单项选择题 根据GB/T 39786《信息安全技术 信息系统密码应用…