单基因gsea_筛到5分的核心基因以后你可以怎么做?

这一次我们从一些已经发表的文章拆解,我们来看看,你找到了一个核心基因以后,你可以怎么做呢?我们就不说那么多废话了,直接用几篇文章的解读来带着大家领会一下如何去进行下一步的分析。


Case1:预后标志物+免疫浸润

第一篇文章是2019年发表在 Front Oncol(IF=4.1437)名为:UBASH3B Is a Novel Prognostic Biomarker and Correlated With Immune Infiltrates in Prostate Cancer的文章。

3aef516d77c8db3e6ffe80e713536f03.png
c6a5390d49abdbe7665f186814e62c56.png

不得不说,随着免疫检查点抑制剂在各种肿瘤治疗中大放异彩,和免疫检查点抑制剂疗效相关的一些标志物也一起受到了广泛的关注,比如说我们的第一个case“免疫浸润”就是一个很好的例子。

免疫浸润不仅仅常用到肿瘤的免疫检查点抑制剂疗效中,还有其他的疾病的发生发展也都和“免疫浸润”息息相关,学习方法和套路;然后从这个所谓的“套路思维”中跳出来,这才是希望大家可以达到的。


  • 通过肿瘤组织和正常组织的mRNA表达:发现核心基因在肿瘤组织中显著升高;
  • 在转移的组织和恶性程度更高的组织中核心基因也更高;

其实我觉得还不如用Cox来筛选来得好,有那么多基因都满足上面的两点,作者为什么要研究UBASH3B这个基因?或者是通过高通量筛选出只有UBASH3B满足上面的两个标准,这样的方法更为严谨。

  • 通过核心基因的高低表达把人群分成2组做差异分析并得到差异基因(DEGs)
  • 通过上面的得到的DEGs进行通路富集来解释为什么会引起生存差异
  • 然后发现DEGs和显著激活差异的通路免疫浸润相关
  • 最后做了一下核心基因和免疫浸润相关基因以及通路的相关性分析

这样的思路的确是比较简单粗暴的,肯定是可以这样的做的,但是显然也可以做得更好

小结:

首先,可以建议大家可以补做CIBERSORT或者是xCell这样的分析,不建议做Timer或者是Estimate,因为信息量太少了。其次,作者是有用到自己的数据的,但是如果能有20-30个样本并且是带随访的RNAseq数据来进行筛选,感觉说服力就强很多了,当然作者为了增强说服力,从多种维度进行验证,弥补了缺陷。最后,如果大家能有现成的基因敲除鼠,或者是能够有免疫重建的PDX模型,那么能在现在的基础上补做一些机制实验,在10分左右还是很有希望的。


Case2:单基因多组学验证生存

第二篇文章是2020年发表在J Clin Med(IF=5.688)上名为Opposite Roles of BAP1 in Overall Survival of Uveal Melanoma and Cutaneous Melanoma的文章;

98efc3716e07961ff112a6baeddbc4bd.png
567b1f1a210f21304f6d06671dba79ed.png

这篇文章很简单,大致内容如下:

  • 核心基因的低表达组和缺失组有着更坏的生存;
  • 然后用Cox回归分析也同样证明了第一步的结论
  • 扩展分析了核心基因的表达量和诊断年龄显著相关

是的,就这样戛然而止了,很诧异。反思了一下,这个可能是一个“约稿”,毕竟JCM是MDPI的杂志,他们杂志就爱搞这种“特刊+约稿”的事情。当然,按照题目,作者发现了核心基因分别在uveal melanoma (UM) and cutaneous melanoma (CM)的生存意义是不一样的…这也算是一个重大发现吧…

那么我们从这篇文章的分析跳出来,看看我们还可以做哪些东西,从而做得更好呢?如果找到一个和生存相关的基因,接下来除了重复性质的验证,你还可以做什么呢?

  • 找到最合适的临界值,建议xTile或者是tROC,而不是直接的median;
  • GSEA或者是ssGSEA来解释生存预后以及和核心基因的关系;
  • 横向分析有很多,除了作者做的年龄;还有可以想到的临床分期分型肿瘤大小,一些热门的score
  • 从核心基因本身出发探寻通路和互作关系;

本来有一篇文章要分析,但是因为篇幅性质,就没有放进来。那是讲的筛选出来的基因是一个自噬相关的基因。筛选出那个基因了以后,就可以看看那个基因和自噬相关的通路那些是相关的;如果运气好还可以在GEO看看,说不定能找到一些别人调控过该基因的分析;

当然还可以通过疾病类型进行扩展分析,比如说研究NSCLC的时候,把LUAD和LUSC分开来做亚组补充分析的信息量;


Case3:基因突变+TMB+免疫浸润

第二篇文章是2020年发表在Aging(IF=5.515)上名为EP300 mutation is associated with tumor mutation burden and promotes antitumor immunity in bladder cancer patients的文章。

7c9e5e4c03733edbca42abc438a8dbf4.png
9004dae4cd0822e8c8378f776b8922ac.png

这篇文章的核心基因是EP300,但是它和上面的文章不一样的是,它是做的基因突变,当然随着多组学数据的开放,你可以做拷贝数变异,可以做蛋白芯片,可以做甲基化,可以做乙酰化等;RNA seq的数据最多并且干预调控比较容易,大部分的时候大家都还是以RNAseq的数据为核心方式进行展开自己的文章。

这篇文章为什么要研究EP300 mutation,作者一步一步筛选出来就很有逻辑,大家可以看一看:

定义Frequently Mutated Gene(FMG)→把TCGA和ICGC两个数据集中的FMGs取交集(11个基因)→TMB的差异分析(9个显著)→生存分析中EP300 mutation显著

那么全文思路大致如下:

  • 筛选了核心基因:和TMB以及生存存在关联;
  • Cox回归分析验证核心基因和生存的关系;
  • GSEA分析看核心基因突变与否和通路激活程度的关系;
  • 核心基因突变与否和CIBERSORT的免疫浸润之间的关系;

讨论也比较简单,这篇文章比较大的一个亮点是作者很有逻辑地筛选出来了EP300 mutation,至于后面的分析的确是比较少了;该作者既没有看那个基因突变是否和其他组学的关联,也没有看是否有热点的突变;扩展的分析也就止步于GSEA和CIBERSORT。

其实,我在想如果是能够从GSEA分析中找到一类通路,然后和CIBERSORT的结合起来讲一个故事可能会更好或者是找到一些通路,进行干预他们对一些表型的关系?如果是有一些临床的样本,结合上现实生活中ICIs治疗的案例,配上case report的一些描述和影像资料,我觉得文章还能再上一个档次。


三个案例分享完了,希望大家能够从里面学到很多新的东西,请点下右下角的在看哟,哈哈哈其实厉害的同学会说,这不就是所谓的单基因套路吗?那我问你套路香吗?

当然香呀!!!!!

并且套路都没学会就开始高级创新了??????

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

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

相关文章

安卓 原生okhttp使用get与post获取网络数据

网址 https://square.github.io/okhttp/ 配置 依赖 在module的build.gradle中&#xff1a; implementation com.squareup.okhttp3:okhttp:3.14.7implementation com.squareup.okio:okio:1.17.5AndroidManifest.xml <uses-permission android:name"android.permissi…

STL源码剖析 迭代器的概念和traits编程技法

迭代器&#xff1a;依序巡防某个聚合物(容器)所含的各个元素&#xff0c;但是不需要暴露这个聚合物的内部表述方式核心思想&#xff1a;将容器和算法分开&#xff0c;彼此独立设计容器和算法的泛型化&#xff0c;均可以使用模板&#xff0c;使用迭代器连接容器和算法例子 templ…

.sql文件如何执行_干货|一条SQL查询语句是如何执行的

作者&#xff1a;wanber链接&#xff1a;https://blog.nowcoder.net/n/9e120e8f1314466bb44fe706b283dc20

STL源码剖析 5中迭代器型别

最常使用的5种迭代器的型别 为 value_type、difference_type、pointer、reference、iterator_category。如果想要自己开发的容器和STL进行适配&#xff0c;就需要定义上述5种类型 iteraor_traits 必须针对传入的型别为 pointer 或者 pointer-to-const设计偏特化版本 template &…

Python学习16 正则表达式3 练习题

用户名匹配 1.用户名匹配&#xff1a;由数字、大小写字母、下划线_、中横线-组成&#xff0c;长度为6-12位&#xff0c;不能以数字开头。 import re usernameab578_-SDF resultre.search(^[a-zA-Z_-][0-9a-zA-Z_-]{5,12}$,username) print(result)邮箱 2.验证输入的邮箱&…

加载tf模型 正确率很低_深度学习模型训练全流程!

↑↑↑关注后"星标"Datawhale每日干货 & 每月组队学习&#xff0c;不错过Datawhale干货 作者&#xff1a;黄星源、奉现&#xff0c;Datawhale优秀学习者本文从构建数据验证集、模型训练、模型加载和模型调参四个部分对深度学习中模型训练的全流程进行讲解。一个成…

Python学习17 Turtle库绘图

学习网址&#xff1a;https://docs.python.org/zh-cn/3/library/turtle.html Turtle库 Turtle库是Python语言中一个很流行的绘制图像的函数库&#xff0c;一个小乌龟&#xff0c;在一个横轴为x、纵轴为y的坐标系原点&#xff08;画布中心&#xff09;&#xff0c;(0,0)位置开…

android ros 节点编写_嵌入式的我们为什么要学ROS

前言本来是要写一篇STM32移植ROS的一个小lib库&#xff0c;ROS一般都是需要跑在Linux上的&#xff0c;STM32使用就是当成一个ROS通讯的小节点&#xff0c;但是写文章时间不够&#xff0c;所以就简单做一篇ROS的介绍文章&#xff0c;分享给嵌入式的小伙伴们。ROS现在在机器人领域…

STL源码剖析 __type_traits

traits编程 弥补了C本身的不足STL只对迭代器进行规范制定出了iterator_traits&#xff0c;SGI在此基础上进一步扩展&#xff0c;产生了__type_traits双下划线的含义是这个是SGI内部使用的东西&#xff0c;不属于STL标准iterator_traits 负责萃取迭代器的特性__type_traits负责萃…

java 学生成绩

题目 对学生成绩大于60分的&#xff0c;输出“合格”。低于60分的&#xff0c;输出“不合格” 代码 使用/除法简化代码 package l1_switch_case;import java.util.Scanner;public class SwitchDemo2 {public static void main(String[] args) {Scanner scanner new Scanne…

STL源码剖析 序列式容器|Vector

容器的概观和分类 array 数组 、list 链表、tree树 、stack堆栈、queue队列、hash table散列表、set集合、map映射表根据数据在容器中的排列顺序&#xff0c;将上述数据结构分为序列式和关联式两种类型SGI STL使用内缩方式来表达基层和衍生层之间的关系衍生不是派生&#xff0…

ansible 修改文件变量_Ansible Playbook中的变量与引用

Ansible是一个系列文章&#xff0c;我会尽量以通俗易懂、诙谐幽默的总结方式给大家呈现这些枯燥的知识点&#xff0c;让学习变的有趣一些。Ansible自动化运维前言前面有说到使用playbook来搞一些复杂的功能&#xff0c;我们使用YAML来写playbook&#xff0c;就像我们用其它语言…

java 判断日期为第几天

题目1 编写程序&#xff1a;从键盘上输入2019年的“month”和“day”&#xff0c;要求通过程序 输出输入的日期为2019年的第几天。 代码1 从12月往下加日期数 package l1_switch_case; import java.util.Scanner; public class SwitchDemo4 {public static void main(Strin…

STL源码剖析 list概述

目录 list的节点(node) list迭代器 list 的构造和内存管理 list 的元素操作 list相较于vector连续的线性空间就显得很复杂&#xff0c;他的存储空间是不连续的&#xff0c;好处是每次插入和删除一个元素的时候&#xff0c;只需要配置或者释放一个元素的空间 插入和删除十分的…

vsftp不允许切换到其它目录_IntelliJ IDEA如何对project的目录进行筛选显示?

如果你的项目很庞大&#xff0c;同一个功能用到的各种文件散落在多个文件夹&#xff0c;开发时切换不便&#xff0c;可以利用scope功能&#xff0c;只显示该功能用到的文件&#xff0c;让project列表十分清爽&#xff0c;提高开发效率。本文使用的IDEA版本为2020.1。1、打开sco…

java 年份对应的中国生肖

题目 编写一个程序&#xff0c;为一个给定的年份找出其对应的中国生肖。 中国的生肖基于12年一个周期&#xff0c; 每年用一个动物代表&#xff1a; rat、ox、tiger、rabbit、dragon、snake、horse、sheep、monkey、 rooster、dog、pig。 提示&#xff1a;2019年&#xff1a;猪…

密码学专题 对称加密算法

一般来说&#xff0c;使用OpenSSL对称加密算法有两种方式&#xff0c;一种是使用API函数的方式&#xff0c;一种是使用OpenSSL提供的对称加密算法指令方式。本书将介绍对称加密算法的指令方式OpenSSL的对称加密算法指令主要用来对数据进行加密和解密处理&#xff0c;输入输出的…

网络防火墙单向和双向_单向晶闸管与双向晶闸管之间的不同之处

晶闸管是回一个可以控导点开关&#xff0c;能以弱电去控制强电的各种电路。晶闸管常用于整流&#xff0c;调压&#xff0c;交直流变化&#xff0c;开关&#xff0c;调光等控制电路中。具有提交小&#xff0c;重量轻&#xff0c;耐压高&#xff0c;容量大&#xff0c;效率高&…

java 遍历100以内的偶数,偶数的和,偶数的个数

题目 遍历100以内的偶数&#xff0c;偶数的和&#xff0c;偶数的个数 代码 package l2_for; /*遍历100以内的偶数&#xff0c;偶数的和&#xff0c;偶数的个数*/ public class ForDemo1 {public static void main(String[] args) {//方法1&#xff1a;int sum1 0,count10;f…

python版本切换_怎么切换python版本

展开全部 &#xff08;1&#xff09;分别安2113装 python-2.7.12.amd64.msi python-3.5.2-amd64.exe &#xff08;python官网下载的&#xff09; 顺序无所谓&#xff08;为5261了看着4102方便&#xff0c;我把安装路径修改统一了1653&#xff09; &#xff08;2&#xff09;配置…