CN112614546B - 一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 - Google Patents
一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 Download PDFInfo
- Publication number
- CN112614546B CN112614546B CN202011573822.0A CN202011573822A CN112614546B CN 112614546 B CN112614546 B CN 112614546B CN 202011573822 A CN202011573822 A CN 202011573822A CN 112614546 B CN112614546 B CN 112614546B
- Authority
- CN
- China
- Prior art keywords
- model
- genes
- gene
- expression level
- tcga
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 206010073071 hepatocellular carcinoma Diseases 0.000 title claims abstract description 46
- 231100000844 hepatocellular carcinoma Toxicity 0.000 title claims abstract description 41
- 238000009169 immunotherapy Methods 0.000 title claims abstract description 26
- 238000010276 construction Methods 0.000 title abstract description 12
- 230000000694 effects Effects 0.000 title abstract description 11
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 56
- 101150096316 5 gene Proteins 0.000 claims abstract description 24
- 108010088350 Lactate Dehydrogenase 5 Proteins 0.000 claims abstract description 18
- 102100035831 Filensin Human genes 0.000 claims abstract description 14
- 101000873715 Homo sapiens Filensin Proteins 0.000 claims abstract description 14
- 239000000835 fiber Substances 0.000 claims abstract description 11
- 108010014790 DAX-1 Orphan Nuclear Receptor Proteins 0.000 claims abstract description 9
- 102100039019 Nuclear receptor subfamily 0 group B member 1 Human genes 0.000 claims abstract description 9
- 108020005497 Nuclear hormone receptor Proteins 0.000 claims abstract description 6
- 239000011324 bead Substances 0.000 claims abstract description 6
- 108020004017 nuclear receptors Proteins 0.000 claims abstract description 6
- 101710085955 Structural protein 1 Proteins 0.000 claims abstract description 5
- 101150041530 ldha gene Proteins 0.000 claims abstract 2
- 230000014509 gene expression Effects 0.000 claims description 47
- 238000004393 prognosis Methods 0.000 claims description 22
- 108020004999 messenger RNA Proteins 0.000 claims description 15
- 102100034671 L-lactate dehydrogenase A chain Human genes 0.000 claims description 14
- 238000000611 regression analysis Methods 0.000 claims description 13
- 230000004083 survival effect Effects 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 13
- 238000010200 validation analysis Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 9
- 102100022627 Fructose-2,6-bisphosphatase Human genes 0.000 claims description 8
- 101000823456 Homo sapiens Fructose-2,6-bisphosphatase Proteins 0.000 claims description 8
- 238000000034 method Methods 0.000 claims description 8
- 102000006255 nuclear receptors Human genes 0.000 claims description 5
- 102000004195 Isomerases Human genes 0.000 claims description 2
- 108090000769 Isomerases Proteins 0.000 claims description 2
- 239000000090 biomarker Substances 0.000 abstract description 3
- 239000003814 drug Substances 0.000 abstract description 3
- 102100037458 Dephospho-CoA kinase Human genes 0.000 abstract 1
- 102000007399 Nuclear hormone receptor Human genes 0.000 abstract 1
- 108010049285 dephospho-CoA kinase Proteins 0.000 abstract 1
- 102000004169 proteins and genes Human genes 0.000 description 9
- 210000001519 tissue Anatomy 0.000 description 8
- 238000003364 immunohistochemistry Methods 0.000 description 7
- 206010028980 Neoplasm Diseases 0.000 description 6
- 201000007270 liver cancer Diseases 0.000 description 6
- 208000014018 liver neoplasm Diseases 0.000 description 6
- 102000008096 B7-H1 Antigen Human genes 0.000 description 4
- 108010074708 B7-H1 Antigen Proteins 0.000 description 4
- 238000010824 Kaplan-Meier survival analysis Methods 0.000 description 4
- 201000011510 cancer Diseases 0.000 description 4
- 229940076838 Immune checkpoint inhibitor Drugs 0.000 description 3
- 102000037984 Inhibitory immune checkpoint proteins Human genes 0.000 description 3
- 108091008026 Inhibitory immune checkpoint proteins Proteins 0.000 description 3
- -1 PPAT Proteins 0.000 description 3
- 210000001744 T-lymphocyte Anatomy 0.000 description 3
- 239000012274 immune-checkpoint protein inhibitor Substances 0.000 description 3
- 238000002493 microarray Methods 0.000 description 3
- 101150048714 75 gene Proteins 0.000 description 2
- 238000003766 bioinformatics method Methods 0.000 description 2
- 231100000504 carcinogenesis Toxicity 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012308 immunohistochemistry method Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000036962 time dependent Effects 0.000 description 2
- 210000004881 tumor cell Anatomy 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 101150084750 1 gene Proteins 0.000 description 1
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- 208000008964 Chemical and Drug Induced Liver Injury Diseases 0.000 description 1
- 108091007741 Chimeric antigen receptor T cells Proteins 0.000 description 1
- 206010008909 Chronic Hepatitis Diseases 0.000 description 1
- 102000004127 Cytokines Human genes 0.000 description 1
- 108090000695 Cytokines Proteins 0.000 description 1
- 101150027068 DEGS1 gene Proteins 0.000 description 1
- 206010072268 Drug-induced liver injury Diseases 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 206010016654 Fibrosis Diseases 0.000 description 1
- 108010033040 Histones Proteins 0.000 description 1
- 206010027476 Metastases Diseases 0.000 description 1
- 206010061309 Neoplasm progression Diseases 0.000 description 1
- 238000000692 Student's t-test Methods 0.000 description 1
- 108091046869 Telomeric non-coding RNA Proteins 0.000 description 1
- 239000000427 antigen Substances 0.000 description 1
- 102000036639 antigens Human genes 0.000 description 1
- 108091007433 antigens Proteins 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 229940022399 cancer vaccine Drugs 0.000 description 1
- 238000009566 cancer vaccine Methods 0.000 description 1
- 230000005773 cancer-related death Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000007882 cirrhosis Effects 0.000 description 1
- 208000019425 cirrhosis of liver Diseases 0.000 description 1
- 238000002648 combination therapy Methods 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 201000008865 drug-induced hepatitis Diseases 0.000 description 1
- 239000012636 effector Substances 0.000 description 1
- 230000004049 epigenetic modification Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012268 genome sequencing Methods 0.000 description 1
- 208000006454 hepatitis Diseases 0.000 description 1
- 210000002865 immune cell Anatomy 0.000 description 1
- 230000028993 immune response Effects 0.000 description 1
- 230000002055 immunohistochemical effect Effects 0.000 description 1
- 238000012151 immunohistochemical method Methods 0.000 description 1
- 230000001506 immunosuppresive effect Effects 0.000 description 1
- 230000001024 immunotherapeutic effect Effects 0.000 description 1
- 238000011293 immunotherapeutic strategy Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 208000019423 liver disease Diseases 0.000 description 1
- 210000001165 lymph node Anatomy 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 208000008338 non-alcoholic fatty liver disease Diseases 0.000 description 1
- 108091027963 non-coding RNA Proteins 0.000 description 1
- 102000042567 non-coding RNA Human genes 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008506 pathogenesis Effects 0.000 description 1
- 239000013610 patient sample Substances 0.000 description 1
- 238000010837 poor prognosis Methods 0.000 description 1
- 230000037452 priming Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 210000003289 regulatory T cell Anatomy 0.000 description 1
- 238000002271 resection Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000009097 single-agent therapy Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 238000012353 t test Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
- 230000005751 tumor progression Effects 0.000 description 1
- 230000009385 viral infection Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/53—Immunoassay; Biospecific binding assay; Materials therefor
- G01N33/574—Immunoassay; Biospecific binding assay; Materials therefor for cancer
- G01N33/57407—Specifically defined cancers
- G01N33/57438—Specifically defined cancers of liver, pancreas or kidney
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/20—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for electronic clinical trials or questionnaires
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Epidemiology (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Biotechnology (AREA)
- Pathology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Immunology (AREA)
- Urology & Nephrology (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Hematology (AREA)
- Biophysics (AREA)
- Oncology (AREA)
- Microbiology (AREA)
- Gastroenterology & Hepatology (AREA)
- Genetics & Genomics (AREA)
- Hospice & Palliative Care (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cell Biology (AREA)
- Bioethics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
Abstract
本发明属临床医学领域，涉及一种基于5‑基因的预测肝细胞癌免疫治疗疗效的模型及其构建方法。用于预测肝细胞癌免疫治疗疗效的模型，包括5个基因乳酸脱氢酶A(LDHA)、磷脂磷酸磷酰胺异质酶(PPAT)、珠子纤维结构蛋白1(BFSP1)、核受体亚纤维0组B成员1(NR0B1)和6‑磷脂‑2‑激酶/果糖酶‑2，6‑双磷酸酶4(PFPFB4)。本发明还提供了基于5基因的预测肝细胞癌免疫治疗疗效的模型的构建方法。本发明建立了一个新的基于5基因的免疫微环境模型，可以有效预测患者对HCC免疫治疗的效果，这5个基因可以作为潜在的生物标志物在临床应用。
Description
技术领域
本发明属临床医学领域，涉及一种基于5-基因的预测肝细胞癌免疫治疗疗效的模型及其构建方法。
背景技术
肝细胞癌(HCC)是一种高度恶性的癌症，是全球癌症相关死亡的第三大原因。5年生存期和总体存活率低于12％。大多数HCC患者是由肝硬化、慢性肝炎病毒感染、酒精相关肝病、非酒精性脂肪性肝病和药物引起的肝炎进展而来。由于HCC通常不能早期诊断，治疗效果较差。因此，研发新的诊断分子标志物，对改善患者预后十分重要。
免疫微环境在肿瘤发生中起着关键作用，与肿瘤进展和治疗效果相关。全身免疫疗法已显示出对HCC 的疗效，特别是对于没有机会进行切除或肝移植的患者。常见的免疫治疗策略包括嵌合抗原受体T细胞 (CAR-T细胞)、癌症疫苗、细胞因子疗法和免疫检查点抑制剂(ICIs)。目前，ICIs治疗是免疫治疗中最成功的一类，包括单一疗法和联合疗法。例如，临床使用最广泛的PD-1抗体和PD-L1抗体。由T细胞表达的 PD-L1通过PD-1调节淋巴结启动阶段和肿瘤细胞效应期阶段的免疫反应。利用针对PD-1和PD-L1的单克隆抗体，“耗竭”T细胞功能的恢复和免疫抑制调节性T淋巴细胞的耗竭为恶性肿瘤的治疗开辟了新的途径。然而，对HCC中PD-1抗体和PD-L1抗体的疗效进行了研究发现，只有大约25％的高表达PD-1的T细胞的 HCC患者对ICIs有应答，因此，识别对ICIs有良好应答的患者具有重要意义。
使用TNM(tumor-node-metastasis)分类预测HCC预后的传统策略有助于指导HCC临床治疗的决策。然而，他们的预测功效并不令人满意。随着测序技术和生物信息学的进展，结合生物信息学分析的基因组测序技术的使用提高了肿瘤诊断和预后的预测能力。建立基于基因的预后模型可以鉴别癌症和正常组织之间的mRNA表达模式差异。TCGA和NCBI基因表达综合数据库(GEO)中用于疾病预后的长链非编码RNA (lncRNAs)、调控表观遗传修饰的基因以及免疫相关基因的表达谱已得到越来越多的研究。然而，目前还没有成熟的模型可以稳定地预测患者对HCC免疫治疗的反应和预后。
发明内容
本发明提供了一种基于5基因的预测肝细胞癌免疫治疗疗效的模型。
一种用于预测肝细胞癌免疫治疗疗效的模型，包括5个基因乳酸脱氢酶A(LDHA)、磷脂磷酸磷酰胺异质酶(PPAT)、珠子纤维结构蛋白1(BFSP1)、核受体亚纤维0组B成员1(NR0B1)和6-磷脂-2-激酶/ 果糖酶-2，6-双磷酸酶4(PFPFB4)。
另外，本发明还提供了基于5基因的预测肝细胞癌免疫治疗疗效的模型的构建方法。
模型的构建方法，包括如下步骤：
1)筛选与HCC预后相关的免疫相关基因；
2)在TCGA-LIHC训练集中，免疫相关基因和生存数据使用R包“coxph”进行单变量Cox回归分析，鉴定与HCC预后相关的免疫相关性基因(DEGs)；
3)采用LASSO-Cox回归分析进一步缩小目标基因，基于LASSO-Cox回归模型系数乘以mRNA表达水平建立预后基因模型。
具体为：
1)筛选了与HCC预后相关的免疫相关基因4227个；
2)利用TCGA-LIHC训练集，对4227个免疫相关基因的进行单变量Cox回归分析，确定245个基因为HCC OS的潜在预后指标；
3)然后执行LASSO-Cox回归分析，以建立预后模型；确定5个基因：包括乳酸脱氢酶A(LDHA)、磷脂磷酸磷酰胺异质酶(PPAT)、珠子纤维结构蛋白1(BFSP1)、核受体亚纤维0组B成员1(NR0B1) 和6-磷脂-2-激酶/果糖酶-2，6-双磷酸酶4(PFPFB4)；风险评分是五个基因的加权基因表达水平的总和乘以各自的LASSO系数：风险评分＝[0.307×LDHA表达水平]+[0.268×PPAT表达水平+ [0.455×BFSP1表达水平]+[0.234×NROB1表达水平]+[0.109×PFKFB4表达水平]；结果分析显示，5个基因的LASSO系数均为阳性。
进一步的，该构建方法还包括：利用TCGA-LIHC训练集评价基于五基因模型的预测效率。
进一步的，该构建方法还包括：使用TCGA-LIHC验证集、TCGA-LIHC全集和GSE14520数据集验证基于5基因的HCC预后模型。
本发明建立了一个新的基于5基因的免疫微环境模型，并进行生物信息学分析，以评估该模型预测HCC 免疫治疗结果的能力。我们进行了单变量Cox回归分析和(LASSO)-Cox回归分析来构建风险模型。采用风险评分、时间依赖性ROC值和Kaplan-Meier生存分析来评估模型的预后能力。结果表明，我们的模型可以有效预测患者对HCC免疫治疗的效果，这5个基因可以作为潜在的生物标志物在临床应用。
附图说明
图1基因预测模型构建和验证流程图
图2基于TCGA-LIHC训练集，5-基因模型的预测效能
图3基于TCGA-LIHC验证集，5-基因模型的预测效能
图4基于TCGA-LIHC全集，5-基因模型的预测效能
图5基于GSE14520数据集，5-基因模型的预测效能
图6免疫组化法，检测LDHA、PPAT、BFSP1、NR0B1和PFKFB4在肝癌组织中表达
图7 5-基因模型蛋白表达水平和肝癌患者预后关系
图8基于mvigor210数据集，5-基因模型对肝癌的免疫治疗效果的预测效能
具体实施方式
(一)方法和材料
1、数据采集
共从TCGA-LIHC数据集中检索365例HCC患者样本，用于分析预后基因表达特征。对365个样本(包括训练集[n＝219]和验证集[n＝146])进行随机抽样100次，并进行置换。TNM分期、分级、OS、性别、年龄在训练集和验证集之间无显著差异(p>0.05)。GSE14520数据集(n＝221)的临床数据和mRNA表达数据来自 NCBI GEO数据库(https://www.ncbi.nlm.nih.gov/geo/)，免疫治疗数据集(Imvigor210)来自于http://research-pub.gene.com/IMvigor210CoreBiologie.。本研究中使用的所有患者数据均具有完整的临床信息，包括TNM分期、分级、生存时间、性别、年龄和免疫相关基因表达水平。
2、建立基于免疫相关基因的5基因模型
我们筛选了与HCC预后相关的免疫相关基因。HCC发病机制中的免疫相关基因从已发表文献中收集。
在TCGA-LIHC训练集中，免疫相关基因和生存数据使用R包“coxph”进行单变量Cox回归分析。我们鉴定了与HCC预后相关的免疫基因。其次，采用LASSO-Cox回归分析进一步缩小目标基因。该方法能够同时进行变量选择和参数估计，并能较好地解决回归分析中的多重共线性问题。基于LASSO-Cox回归模型系数乘以mRNA表达水平建立预后基因模型。最终获得风险评分＝[0.307×LDHA mRNA表达水平]+[0.268 ×PPAT mRNA表达水平+[0.455×mRNA BFSP1表达水平]+[0.234×mRNA NROB1表达水平]+[0.109 ×mRNA PFKFB4表达水平]。样本按中位数分为低风险组和高风险组。
3、验证五基因模型的性能和预测能力
使用TCGA-LIHC的验证集、TCGA-LIHC全集、GSE14520数据集和免疫治疗数据集(Imvigor210)，采用时间依赖性ROC分析和Kaplan-Meier log-rank分析来评估模型的性能和预后能力。
4、统计分析
使用SPSS v25(IBM,Chicago,IL,USA)、GraphPad Prism 7.0(GraphPadSoftware,La Jolla,CA,USA)和R软件 (version 3.5.1)进行统计分析。采用t检验进行统计学比较，采用Kaplan-Meier分析法估计OS,p<0.05认为有统计学意义。
(二)结果
1、筛选HCC预后基因表达特征构建HCC预后模型
图1显示了分析工作流程图(5-基因预测模型构建和验证流程图)。利用TCGA-LIHC训练集，对4227 个免疫相关基因的进行单变量Cox回归分析，筛选出245个免疫相关基因为HCC OS的潜在预后指标。
2、基于五基因的HCC预后模型的构建
然后执行LASSO-Cox回归分析，以建立预测模型。确定了5个基因：包括乳酸脱氢酶A(LDHA)、磷脂磷酸磷酰胺异质酶(PPAT)、珠子纤维结构蛋白1(BFSP1)、核受体亚纤维0组B成员1(NR0B1)和6- 磷脂-2-激酶/果糖酶-2，6-双磷酸酶4(PFPFB4)。风险评分是五个基因的加权基因表达水平的总和乘以各自的LASSO系数：风险评分＝[0.307×mRNA LDHA表达水平]+[0.268×mRNA PPAT表达水平+[0.455 ×mRNA BFSP1表达水平]+[0.234×mRNA NROB1表达水平]+[0.109×mRNA PFKFB4表达水平]。结果分析显示，5个基因的LASSO系数均为阳性。
3、利用TCGA-LIHC训练集评价基于五基因模型的预测效率
为了确定这5个基因的基因表达特征与HCC患者生存结果之间的关系，分别计算每个样本基于5个基因的模型风险评分(AUC值)(图2.基于TCGA-LIHC训练集，5-基因模型的预测效能)，并使用TCGA-LIHC训练集定义风险评分的最佳截止点(图2A)。AUC值越高，表明基于5基因的HCC预后模型的分类性能越好。在1年、3年和5年的生存率中，AUC值分别为0.80、0.77和0.73(图2B)。Kaplan-Meier生存分析显示，高危组患者预后较低危组差(p<0.0001)；图2C)。综上所述，这些结果表明建立的基于5基因的模型在预测HCC预后方面表现良好。
4、使用TCGA-LIHC验证集、TCGA-LIHC全集和GSE14520数据集验证基于5-基因的HCC预后模型。
为了验证基于5基因模型的稳定性，我们对3个数据集(TCGA-LIHC验证集、TCGA-LIHC全集和GSE14520数据集)进行了类似的流程分析。通过计算各自ROC曲线的AUC值，得到基于五基因的模型的风险评分。 TCGA-LIHC验证集1年、3年和5年生存率的AUC值分别为0.70、0.69和0.60。结果显示，高危组患者的生存率明显低于低危组，(p<0.001)；(图3基于TCGA-LIHC验证集，5-基因模型的预测效能)。与TCGA-LIHC 验证集的结果一致，高危组患者的OS较低(p均<0.001)，整个TCGA-LIHC数据集和GSE14520数据集的AUC 均大于0.6(图4基于TCGA-LIHC全集，5-基因模型的预测效能，和图5基于GSE14520数据集，5-基因模型的预测效能)。结果显示，基于5个基因的模型可以从基因表达水平预测患者的生存时间。
实施例2
组织微阵列(TMA)结构五基因模型评估
为了评估5基因模型的预后能力，我们使用了一个90例HCC患者癌组织的组织微阵列(上海芯超生物科技有限公司)进行免疫组织化学(IHC)，并测量集成光密度分析。所使用的一抗见表1。IHC评分是由三位资深病理学家在不了解患者特征的情况下独立评估得出的。IHC得分计算:包含IHC score＝[0.307×protein LDHA表达水平]+[0.268×proteinPPAT表达水平+[0.455×protein BFSP1表达水平]+[0.234×protein NROB1表达水平]+[0.109×protein PFKFB4表达水平]。采用Kaplan-Meier log-rank分析评价高IHC评分组和低IHC评分组的生存率差异。
表1.抗体信息表
用免疫组化法检测组织微阵列芯片(TMA)5个基因LDHA、PPAT、BFSP1、NR0B1和PFKFB4的蛋白表达水平。结果发现这5个基因的蛋白水平在HCC组织中存在显著差异(图6免疫组化法检测LDHA、PPAT、BFSP1、 NR0B1和PFKFB4在肝癌组织中表达)。免疫组化评分较高的组蛋白表达水平较高，疾病预后较差(图7 5-基因模型蛋白表达水平和肝癌患者预后关系)。结果表明，五基因模型、分级及TNM分期均有统计学意义(p< 0.05)。这5个基因的高蛋白表达水平预示着较差的预后，低蛋白表达水平预示着相对较好的预后。
目前缺乏能够有效预测免疫治疗药物疗效的生物标志物。因此，识别新的预测5-标记对进一步提高免疫治疗的精确性是必要的。检索抗PD-L1免疫治疗患者治疗反应数据的转录组数据(Imvigor210)，评估基于五基因模型的免疫治疗预测效果(图8基于mvigor210数据集，5-基因模型对肝癌的免疫治疗效果的预测效能)。Kaplan-Meier分析显示，风险评分值越高，生存率越低(图8a)。ROC曲线分析显示，合并肿瘤新生抗原(NEO)组、肿瘤突变负荷(tumor mutational burden,TMB)组和风险评分预后因素组的AUC值(AUC＝0.91)高于单独NEO组(AUC＝0.7)、TMB组(AUC＝0.65)和风险评分组(AUC＝0.54)(图8b)。风险评分与免疫治疗疗效无相关性(p>0.05；图8c)，g免疫细胞和肿瘤细胞亚群相关(p<0.05)；图8d、8e)。这些结果表明，将基于5基因的模型与NEO和TMB相结合，可以增强免疫治疗疗效的评估，识别对免疫治疗有反应的患者。
Claims (3)
1.用于预测肝细胞癌免疫治疗疗效的模型的构建方法，所述模型 包括5个基因乳酸脱氢酶A（LDHA）、磷脂磷酸磷酰胺异质酶（PPAT）、珠子纤维结构蛋白1（BFSP1）、核受体亚纤维0组B成员1（NR0B1）和6-磷脂-2-激酶/果糖酶-2，6-双磷酸酶4（PFPFB4），包括如下步骤：
1）筛选与HCC预后相关的免疫相关基因4227 个；
2）在TCGA-LIHC训练集中，对 4227 个免疫相关基因和生存数据使用R包“coxph”进行单变量Cox回归分析，鉴定与HCC预后相关的免疫相关性基因，确定245个基因为HCC OS的潜在预后指标；
3）采用LASSO-Cox回归分析进一步缩小目标基因，基于LASSO-Cox回归模型系数乘以mRNA表达水平建立预后基因模型；确定5个基因：包括乳酸脱氢酶A（LDHA）、磷脂磷酸磷酰胺异质酶（PPAT）、珠子纤维结构蛋白1（BFSP1）、核受体亚纤维0组B成员1（NR0B1）和6-磷脂-2-激酶/果糖酶-2，6-双磷酸酶4（PFPFB4）；风险评分是五个基因的加权基因表达水平的总和乘以各自的 LASSO 系数：风险评分 = [0.307 × LDHA 表达水平]+ [0.268 ×PPAT表达水平] + [0.455 × BFSP1表达水平] + [0.234 × NROB1 表达水平] + [0.109 ×PFKFB4表达水平] ；结果分析显示，5个基因的LASSO系数均为阳性。
2.根据权利要求1所述构建方法，该方法还包括：利用TCGA-LIHC训练集评价基于五基因模型的预测效率。
3.根据权利要求1所述构建方法，该方法还包括：使用TCGA-LIHC验证集、TCGA-LIHC全集和GSE14520数据集验证基于5基因的HCC预后模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011573822.0A CN112614546B (zh) | 2020-12-25 | 2020-12-25 | 一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011573822.0A CN112614546B (zh) | 2020-12-25 | 2020-12-25 | 一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112614546A CN112614546A (zh) | 2021-04-06 |
| CN112614546B true CN112614546B (zh) | 2022-09-02 |
Family
ID=75248056
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011573822.0A Expired - Fee Related CN112614546B (zh) | 2020-12-25 | 2020-12-25 | 一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112614546B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113161000B (zh) * | 2021-05-06 | 2024-05-28 | 复旦大学附属中山医院 | 混合细胞型肝癌的预后评分模型及其构建方法 |
| CN113930506B (zh) * | 2021-09-23 | 2022-10-18 | 江苏大学附属医院 | 一种预测肝细胞癌预后和治疗抵抗的谷氨酰胺代谢基因标签评分系统 |
| CN115019880B (zh) * | 2022-05-05 | 2024-01-09 | 中山大学附属第一医院 | 一种肝癌预后模型及其构建方法和应用 |
| CN116844685B (zh) * | 2023-07-03 | 2024-04-12 | 广州默锐医药科技有限公司 | 一种免疫治疗效果评估方法、装置、电子设备及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017060319A1 (en) * | 2015-10-09 | 2017-04-13 | F. Hoffmann-La Roche Ag | Gene expression biomarkers for personalized cancer care to epigenetic modifying agents |
| CN110580956A (zh) * | 2019-09-19 | 2019-12-17 | 青岛市市立医院 | 一组肝癌预后标志物及其应用 |
| CN111402949A (zh) * | 2020-04-17 | 2020-07-10 | 北京恩瑞尼生物科技股份有限公司 | 一种肝细胞肝癌患者诊断、预后和复发统一模型的构建方法 |
| CN112011616A (zh) * | 2020-09-02 | 2020-12-01 | 复旦大学附属中山医院 | 预测肝细胞癌肿瘤免疫浸润和术后生存时间的免疫基因预后模型 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2854665A1 (en) * | 2011-11-10 | 2013-05-16 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Gene expression signatures of neoplasm responsiveness to therapy |
| WO2016126883A1 (en) * | 2015-02-03 | 2016-08-11 | Cedars-Sinai Medical Center | Biomarker based prognostic model for predicting overall survival in patients with metastatic clear cell kidney cancer |
-
2020
- 2020-12-25 CN CN202011573822.0A patent/CN112614546B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017060319A1 (en) * | 2015-10-09 | 2017-04-13 | F. Hoffmann-La Roche Ag | Gene expression biomarkers for personalized cancer care to epigenetic modifying agents |
| CN110580956A (zh) * | 2019-09-19 | 2019-12-17 | 青岛市市立医院 | 一组肝癌预后标志物及其应用 |
| CN111402949A (zh) * | 2020-04-17 | 2020-07-10 | 北京恩瑞尼生物科技股份有限公司 | 一种肝细胞肝癌患者诊断、预后和复发统一模型的构建方法 |
| CN112011616A (zh) * | 2020-09-02 | 2020-12-01 | 复旦大学附属中山医院 | 预测肝细胞癌肿瘤免疫浸润和术后生存时间的免疫基因预后模型 |
Non-Patent Citations (3)
| Title |
|---|
| Exploration and Verification of a Six-RNA Binding Proteins-based Prognosis Evaluating Model for Hepatocellular Carcinoma;Miao Chen;《Research Square》;20201203;第1-30页 * |
| Identification and Validation of a Novel Metabolism‑Related Prognostic Signature in Hepatocellular Carcinoma;Zhihao wang;《Research Square》;20200819;第1-32页 * |
| 基于生物信息学分析PKM2在肺腺癌中的表达及临床意义;张轶等;《肿瘤预防与治疗》;20200925(第09期);43-49 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112614546A (zh) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112614546B (zh) | 一种用于预测肝细胞癌免疫治疗疗效的模型及其构建方法 | |
| AU2015249113B2 (en) | Lung cancer biomarkers and uses thereof | |
| CN113450873B (zh) | 一种预测胃癌预后和免疫治疗适用性的标志物及其应用 | |
| KR101921945B1 (ko) | 폐암 바이오마커 및 그것의 용도 | |
| CN111394456B (zh) | 早期肺腺癌患者预后评估系统及其应用 | |
| US20120101002A1 (en) | Lung Cancer Biomarkers and Uses Thereof | |
| CN112133365A (zh) | 评估肿瘤微环境的基因集、评分模型及其应用 | |
| JP2016500512A (ja) | 肝臓サンプルの分類ならびに限局性結節性異形成、肝細胞腺腫および肝細胞癌の診断のための新規方法 | |
| CN113462776B (zh) | m6A修饰相关联合基因组在预测肾透明细胞癌患者免疫治疗疗效中的应用 | |
| CN115410713A (zh) | 一种基于免疫相关基因的肝细胞癌预后风险预测模型构建 | |
| CA2885518A1 (en) | A method for prognosis of global survival and survival without relapse in hepatocellular carcinoma | |
| WO2010063121A1 (en) | Methods for biomarker identification and biomarker for non-small cell lung cancer | |
| CN116042832B (zh) | 一种预测非小细胞肺癌免疫治疗获益程度及预后的生物标志物及应用 | |
| JP2022524484A (ja) | がん患者の生存率を予測する方法 | |
| CN111690740A (zh) | 肝细胞癌早筛试剂盒及其制备方法和用途 | |
| CN113430266A (zh) | G6pc及其所在基因组在制备肾透明细胞癌诊断或预后评估系统中的应用 | |
| CN112961920A (zh) | 一种预测肝胆肿瘤患者免疫检查点抑制剂联合靶向治疗疗效的方法及产品 | |
| CN113528509A (zh) | 一种用于预测胃癌免疫治疗的肿瘤微环境评分系统的构建方法以及分子探针 | |
| CN117038067A (zh) | 神经内分泌型前列腺癌风险预测方法及其应用 | |
| US20220065872A1 (en) | Lung Cancer Biomarkers and Uses Thereof | |
| CN112481380A (zh) | 一种评估晚期膀胱癌抗肿瘤免疫治疗反应性和预后生存的标志物及其应用 | |
| Feng et al. | Aggrephagy-related LncRNAs index: A predictor for HCC prognosis, immunotherapy efficacy, and chemosensitivity | |
| CN117187386A (zh) | 一种非小细胞肺癌免疫治疗预后相关基因、预后评分模型及其应用 | |
| CN116403648B (zh) | 一种基于多维分析建立的小细胞肺癌免疫新分型方法 | |
| US20240117441A1 (en) | Tumor microenvironment by liquid biopsy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220902 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |