Cover Image

🤔 你的基因真的是个好Biomarker吗?

在筛选出几个感兴趣的基因后,我们通常会面临一连串的灵魂拷问: “它在治疗前后有差异吗?” “它在应答者和非应答者中表达不一样吗?” “它能预测生存期吗?” “它的预测效果比现有的签名(如TIDE, IFNG)更好吗?”

为了回答这些问题,你可能需要画箱线图、做生存分析、画ROC曲线、去查阅几十篇文献找别人的签名公式然后手算……

这一套下来,没个几天搞不定。而且,万一公式输错了呢?

⚡ tigeR:Biomarker评估的加速器

tigeR包把这一整套评估流程封装得行云流水。不管你是想看单个基因的表现,还是想和已发表的Signature一决高下,tigeR都能轻松搞定。

1. 综合分析:一张图看懂所有(integrate_analysis)

不想分开跑差异分析和生存分析?integrate_analysis函数让你一次看个够。

library(tigeR)
# 输入数据集和感兴趣的基因(如PD-L1/CD274)
integrate_analysis(SE = MEL_GSE91061, geneSet = "CD274")

这个函数会直接返回差异表达(Responders vs Non-Responders)和生存分析的结果。简单、粗暴、有效。

2. 差异可视化:颜值即正义(diff_biomk)

想画出可以直接放在文章里的图?diff_biomk满足你。

治疗前后对比:

diff_biomk(SE = MEL_GSE91061, gene = 'CD274', type = 'Treatment')

应答者 vs 非应答者:

diff_biomk(SE = MEL_GSE91061, gene = 'CD274', type = 'Response')

无需调整ggplot2的繁琐参数,默认配色就已经足够专业美观。

3. 生存分析:一键KM曲线(surv_biomk)

P <- surv_biomk(SE = MEL_GSE91061, gene = 'CD274')
P$plot # 显示KM曲线

自动进行最佳截断值(cutoff)选择,计算P值,绘制风险表。你需要的,它都有。

4. 🌟 杀手锏:23种内置免疫评分一键计算(score_biomk)

这绝对是tigeR最让人生畏的功能。

你是否还在为复现别人的Signature而头秃?tigeR内置了23种主流的免疫治疗相关Signature,包括:

  • TIDE相关: T cell-inflamed GEP, IFNG_Sig
  • TME相关: Angiogenesis, Teffector, Myeloid_inflammatory
  • 以及: IPS, CD8Teffector, TLS, EMT等等

只需要一个函数:

# 计算所有23种评分
sig_res <- score_biomk(MEL_GSE78220)

返回的就是一个整整齐齐的矩阵,每一列是一个Signature的得分。有了这个,你就可以轻松地把你的基因和这些”大佬”们进行相关性分析,或者比较谁的AUC值更高。

5. ROC曲线大比拼(roc_biomk)

说到AUC,tigeR也为你准备好了roc_biomk。你可以方便地评估你的基因或者内置Signature的预测性能。

# 评估 Tcell_inflamed_GEP 的预测效果
sig_roc <- roc_biomk(MEL_PRJEB23709,
                     Weighted_mean_Sigs$Tcell_inflamed_GEP,
                     method = "Weighted_mean",
                     rmBE = TRUE,
                     response_NR = TRUE)
                     
# 查看ROC曲线
sig_roc[[2]]

🎯 总结

tigeR包将Biomarker评估的复杂度降到了最低。

  • 全方位:差异、生存、ROC一网打尽。
  • 广覆盖:内置23种经典Signature,省去手动复现的痛苦。
  • 高颜值:出图直接达到发表级要求。

如果你正在寻找新的肿瘤免疫治疗Biomarker,tigeR绝对是你不可或缺的验证工具。

下一期,我们将进入更微观的世界,看看tigeR如何用10种算法解构肿瘤微环境(TME)


参考资料

  1. tigeR Github: https://github.com/YuLab-SMU/tigeR
  2. iMeta Paper: http://doi.org/10.1002/imt2.229