基于色氨酸代谢基因的肝癌预后评估方法技术

本发明专利技术公开了一种基于色氨酸代谢基因的肝癌预后评估方法，基于公式（1）计算肝癌患者的风险指数；根据阈值“0”划分RiskScore高低风险组，采用Kaplan

【技术实现步骤摘要】
基于色氨酸代谢基因的肝癌预后评估方法


[0001]本专利技术涉及肝癌研究领域，尤其是涉及基于色氨酸代谢基因的肝癌预后评估方法。

技术介绍

[0002]肝癌是异质性高，死亡率最高的难治性恶性肿瘤，其死亡率居肿瘤死亡率第二位，仅次于肺癌。针对肝癌的治疗手段有很多，但治疗后的五年生存率极低，仅为12％左右，有一部分原因是目前不能通过有效手段评估肝癌治疗手段的预后情况，导致肝癌的临床治疗方案具有一定的盲目性，无法针对肝癌患者进行个性化治疗，也无法有效分析治疗数据，为优化肝癌临床治疗方案提供数据支持。

技术实现思路

[0003]本专利技术目的在于提供一种基于色氨酸代谢基因的肝癌预后评估方法。
[0004]为实现上述目的，本专利技术采取下述技术方案：本专利技术所述的基于色氨酸代谢基因的肝癌预后评估方法，包括以下内容：基于公式（1）计算肝癌患者的风险指数；根据阈值“0”划分RiskScore高低风险组，采用Kaplan
‑
Meier法绘制生存曲线，评估肝癌预后效果；RiskScore=Σβ
i
×
Exp
i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ 公式（1）；其中，i代表基于色氨酸代谢基因与肝癌预后显著相关的第i个基因，β为基因表达水平的Cox回归系数，Exp为基因的表达水平。
[0005]进一步地，所述基于色氨酸代谢基因与肝癌预后显著相关的基因包括CDK1，TROAP，G6PD，MMP1，BAIAP2L2，PTTG1，LCAT，CYP2C9，CFHR3，SLC22A10。
[0006]进一步地，评估包括预后特征、临床特征、免疫细胞的相对丰度和浸润程度、免疫逃逸的可能性、免疫检查点基因差异表达、对传统的化疗药物的响应程度、基因代谢通路差异性。
[0007]本申请所述方法的具体原理如下：专利技术人从TCGA GDC API下载TCGA
‑
LIHC的突变数据和拷贝数变异数据，以及TCGA
‑
LIHC的RNA
‑
Seq数据。针对TCGA
‑
LIHC的RNA
‑
Seq数据去掉没有临床随访信息；没有生存时间；没有Status状态的样本，共筛选出包含360例原发肿瘤样本和50例正常样本。
[0008]通过GEO数据库获取GSE14520、GSE76427、HCCDB18数据的表达谱数据，经过筛选分别得到242例、115例、389例肝癌样本。筛选原则为：下载对应芯片平台的注释信息，根据注释信息将探针映射到基因上，去除一个探针匹配到多个基因的探针。当多个探针匹配到一个基因时取中位数作为该基因表达值。
[0009]通过MSigDB数据库色氨酸代谢通路“KEGG_TRYPTOPHAN_METABOLISM”获取色氨酸代谢基因。
[0010]首先，通过三个方面验证色氨酸代谢基因与肝癌患者的生存预后的相关性。
[0011]（1）针对色氨酸代谢相关的40个基因，使用hallmark基因集进行了基因集变异分析比较色氨酸代谢基因突变组和非突变组，发现色氨酸代谢基因突变能导致功能改变，提示其影响肝癌患者的生存预后。
[0012]（2）通过原发肝癌患者的色氨酸代谢基因拷贝数变异（英文缩写为：CNV）比较，确定CNV增加的患者比CNV丢失的患者表现出更高的色氨酸代谢基因表达水平，提示色氨酸代谢基因表达水平影响肝癌患者的生存预后。
[0013]（3）通过比较原发肿瘤组织和癌旁正常组织中色氨酸代谢基因的mRNA（messenger RNA, 信使RNA）变化，确定两者色氨酸代谢基因表达水平存在显著差异，同样提示色氨酸代谢基因表达水平影响肝癌患者的生存预后。
[0014]因此，本申请专利技术人提供了基于色氨酸代谢基因的肝癌预后评估方法，用于评估肝癌治疗手段的预后情况，为有效分析治疗数据，优化肝癌临床治疗方案和制定个性化治疗方案提供数据支持。
[0015]专利技术人通过单变量Cox回归分析，确定了40个色氨酸代谢基因中有8个基因与肝癌预后相关，并在计算了与肝癌预后相关的色氨酸代谢基因在肝癌中的表达之间的成对相关后发现，这8个不同的色氨酸代谢基因之间也存在着显著相关性。于是，基于这8个与肝癌预后相关的色氨酸代谢基因表达谱的一致聚类，对患者进行分类，根据累计分布函数（英文缩写为：CDF）确定最佳聚类个数为2，记为cluster 1和cluster 2。
[0016]然后，专利技术人通过单因素分析，包括预后特征、临床特征、免疫细胞的相对丰度和浸润程度、免疫逃逸的可能性、免疫检查点基因差异表达、对传统的化疗药物的响应程度、基因代谢通路差异性，进一步确认cluster 2预后较好，cluster 1预后较差。
[0017]故，通过分析cluster 1和cluster 2的差异基因，找到739个差异基因，通过survival包的coxph函数进行单因素cox分析，共识别出189个对肝癌预后影响较大的基因，再使用lasso回归对这189个差异表达基因进一步压缩，最终选择了lambda= 0.066时的10个基因作为目标监测基因。以这10个基因构建肝癌预后评估模型。基于公式（1）计算肝癌患者的风险指数；根据阈值“0”划分RiskScore高低风险组，即，将公式（1）计算的RiskScore分值以阈值“0”为分界点，划分为高风险组和低风险组。并采用Kaplan
‑
Meier法绘制生存曲线（英文缩写为KM），可用于评估肝癌预后效果。
[0018]本专利技术优点在于根据色氨酸代谢表型的10个关键基因，构建了肝癌临床预后评估模型，该模型具有较强的鲁棒性，并且独立于临床病理学特征，能够独立评估肝癌治疗手段的预后情况，为有效分析治疗数据，优化肝癌临床治疗方案和制定个性化治疗方案提供数据支持。
附图说明
[0019]图1是本专利技术所述方法绘制的KM曲线与ROC曲线对比图。
[0020]图2是本专利技术所述方法的评估结果与临床分级对比图。
[0021]图3是本专利技术所述方法高/低风险组与免疫细胞评分的对比图。
[0022]图4是本专利技术所述方法高/低风险组部分免疫检查点基因表达差异图。
具体实施方式
[0023]下面将对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0024]本专利技术所述的混凝土大坝裂缝可视化交互展示方法，以TCGA数据集样本中10个基于色氨酸代谢基因、与肝癌预后显著相关的基因表达水平，按照公式（1）分别计算每个样本的风险得分，评估患者为高风险组或低风险组，并绘制生存曲线。
[0025]（1）验证本申请所述方法1,3,5年的预后预测评估效率使用 R 软件包 timeROC 对 RiskScore的评估结果进行 ROC 分析本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于色氨酸代谢基因的肝癌预后评估方法，其特征在于：基于公式（1）计算肝癌患者的风险指数；根据阈值“0”划分RiskScore高低风险组，采用Kaplan
‑
Meier法绘制生存曲线，评估肝癌预后效果；RiskScore=Σβ
i
×
Exp
i 公式（1）；其中，i代表基于色氨酸代谢基因与肝癌预后显著相关的第i个基因，β为基因表达水平的Cox回归系数，Exp为基因的表达水平。...

【专利技术属性】
技术研发人员：李兰娟，薛晨，顾心雨，鹿娟，
申请(专利权)人：浙江大学医学院附属第一医院，
类型：发明
国别省市：