免疫分子分型发一区5.3分杂志？这个思路值得借鉴！

2019 年 12 月 30 日
筆記

今天分享的是发表于老牌杂志ONCOIMMUNOLOGY的一篇文章，该杂志影响因子5.3，中科院分区一区。文章标题为Two immune-enhanced molecular subtypes differ in inflammation, checkpoint signaling and outcome of advanced head and neck squamous cell carcinoma，这里用到的免疫方面的分析值得借鉴。

术语

•HNSCC （head and neck squamous cell carcinoma）头颈部鳞状细胞癌

•Consensus Cluster 一致性聚类

• nearest-centroid classifier 最近质心分类

研究思路

• 数据来源：TCGA-HNSCC和GSE65858 (LHNG cohort)

• 基于1703个免疫相关基因对HNSCC进行亚型的划分

• 对分类亚型进行免疫概况（免疫细胞和免疫检查点）的展示

• 在验证集（LHNG cohort）用最近质心分类分成同样的4类后进行相同的分析

• 在测试集和验证集进行非免疫基质细胞和炎症趋化因子表达的分析

• 对4个亚型进行预后的分析（生存分析）

结果

2.1

划分亚型

这里采用（Consensus Cluster）一致性聚类的方法，利用1703个免疫相关基因将TCGA的HSNCC样本划分为4个亚型；Consensus Cluster Plus包，利用无监督聚类的方法，确定数据集中可能的聚类数目(即k，这里用的是4)以及聚类中的成员（即C1、C2、C3和C4包括哪些病人）（更详细的可以看Consensus Cluster Plus的说明书）；PCA、维恩图、热图（top 100 genes）均是另一种形式的聚类结果展示；

2.2

亚型的免疫情况

这里是热图展示不同亚型的特定免疫细胞的表达和检查点通路分子的表达；同时用箱线图展示了6种免疫细胞和8种免疫检查点分子的表达情况，这里主要看p值了解是否为显著差异；

24种免疫细胞的数据集来源于2013年发表的Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer，主要思想是利用数个marker的表达来表征细胞的情况，即从基因层面上升到了基因集的层面；

2.3

在LHNG cohort进行亚型的验证

LHNG cohort在文章多次提到，为GSE65858 数据；即将测试集（TCGA-HNSCC）换成验证集（GSE65858 ），进行上述的相同操作（热图和箱线图）；

TCGA和GEO这样作为测试集和验证集的情况用得比较多；这里对GEO数据集的分类采用了最近质心分类的方法，即取训练集每个亚型的平均值为中心，待分类的样本（验证集中的样本）归属为离得最近的中心点的亚型；

2.4

非免疫基质细胞和免疫趋化因子的表达

在TCGA和LHNG数据集，利用MCP-counter对endothelial cells和 fibroblasts的表达进行计算，此外，还进行了趋化因子的计算，并用箱线图进行表示；

MCP-counter源于2016年发表的文章Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression，R包可方便使用，能够依据表达数据对异质组织中8种免疫细胞和两个基质细胞群体进行定量；

2.5

生存曲线的绘制

为TCGA和LHNG cohort的4 个亚型作生存曲线的绘制，比较4个亚型是否会有预后差异；

结语

文章依然是TCGA作测试集，GEO作验证集展开分析；这里主要是免疫相关的一些分析，利用表达数据进行免疫细胞和免疫检查点的表征，最后再落到不同亚型的预后差异（一般分型后，生存分析总是绕不开的）上；文中提到的分析，可以通过学习相应的R包或使用网页工具来做，总的来说，难度不高，可以学起来。老规矩，原文获取方式，后台回复关键词：5.3

如果您想按照这个套路针对其他癌症做相关生信分析，或您有其他生信分析需求，也可以联系小编哦