CN110197701A - 一种新型多发性骨髓瘤诺模图构建方法 - Google Patents
一种新型多发性骨髓瘤诺模图构建方法 Download PDFInfo
- Publication number
- CN110197701A CN110197701A CN201910325039.3A CN201910325039A CN110197701A CN 110197701 A CN110197701 A CN 110197701A CN 201910325039 A CN201910325039 A CN 201910325039A CN 110197701 A CN110197701 A CN 110197701A
- Authority
- CN
- China
- Prior art keywords
- gene
- expression
- module
- genes
- construction method
- 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.)
- Granted
Links
- 201000010099 disease Diseases 0.000 title claims abstract description 14
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 14
- 238000010276 construction Methods 0.000 title claims abstract description 12
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 60
- 230000014509 gene expression Effects 0.000 claims abstract description 51
- 230000004083 survival effect Effects 0.000 claims abstract description 13
- 102000015736 beta 2-Microglobulin Human genes 0.000 claims abstract description 8
- 108010081355 beta 2-Microglobulin Proteins 0.000 claims abstract description 8
- 238000003012 network analysis Methods 0.000 claims abstract description 7
- 210000002966 serum Anatomy 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- DDRJAANPRJIHGJ-UHFFFAOYSA-N creatinine Chemical compound CN1CC(=O)NC1=N DDRJAANPRJIHGJ-UHFFFAOYSA-N 0.000 claims description 14
- 102000001554 Hemoglobins Human genes 0.000 claims description 8
- 108010054147 Hemoglobins Proteins 0.000 claims description 8
- 229940109239 creatinine Drugs 0.000 claims description 7
- 108010074051 C-Reactive Protein Proteins 0.000 claims description 6
- 102100032752 C-reactive protein Human genes 0.000 claims description 6
- 101001027634 Homo sapiens Kinesin-like protein KIF21B Proteins 0.000 claims description 5
- 101001002191 Homo sapiens Postmeiotic segregation increased 2-like protein 5 Proteins 0.000 claims description 5
- 102100020952 Postmeiotic segregation increased 2-like protein 5 Human genes 0.000 claims description 5
- 210000001185 bone marrow Anatomy 0.000 claims description 5
- 210000004180 plasmocyte Anatomy 0.000 claims description 5
- 101710088194 Dehydrogenase Proteins 0.000 claims description 4
- 102100024175 Centriole, cilia and spindle-associated protein Human genes 0.000 claims description 3
- 102100025832 Centromere-associated protein E Human genes 0.000 claims description 3
- 108010009392 Cyclin-Dependent Kinase Inhibitor p16 Proteins 0.000 claims description 3
- 102100030421 Fatty acid-binding protein 5 Human genes 0.000 claims description 3
- 101000980796 Homo sapiens Centriole, cilia and spindle-associated protein Proteins 0.000 claims description 3
- 101000914247 Homo sapiens Centromere-associated protein E Proteins 0.000 claims description 3
- 101000802406 Homo sapiens E3 ubiquitin-protein ligase ZNRF3 Proteins 0.000 claims description 3
- 101001062855 Homo sapiens Fatty acid-binding protein 5 Proteins 0.000 claims description 3
- 101000837008 Homo sapiens Sigma intracellular receptor 2 Proteins 0.000 claims description 3
- 101000837581 Homo sapiens Ubiquitin-conjugating enzyme E2 T Proteins 0.000 claims description 3
- 102100037690 Kinesin-like protein KIF21B Human genes 0.000 claims description 3
- 108091006941 SLC39A10 Proteins 0.000 claims description 3
- 102100028662 Sigma intracellular receptor 2 Human genes 0.000 claims description 3
- 102100028705 Ubiquitin-conjugating enzyme E2 T Human genes 0.000 claims description 3
- 102000006083 ZNRF3 Human genes 0.000 claims description 3
- 102100035243 Zinc transporter ZIP10 Human genes 0.000 claims description 3
- 238000000034 method Methods 0.000 claims description 3
- 235000012571 Ficus glomerata Nutrition 0.000 claims description 2
- 244000153665 Ficus glomerata Species 0.000 claims description 2
- 102100034613 Annexin A2 Human genes 0.000 claims 1
- 108090000668 Annexin A2 Proteins 0.000 claims 1
- 102100024458 Cyclin-dependent kinase inhibitor 2A Human genes 0.000 claims 1
- 230000004077 genetic alteration Effects 0.000 claims 1
- 231100000118 genetic alteration Toxicity 0.000 claims 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 239000000523 sample Substances 0.000 description 25
- 208000034578 Multiple myelomas Diseases 0.000 description 19
- 206010035226 Plasma cell myeloma Diseases 0.000 description 19
- 238000012360 testing method Methods 0.000 description 19
- 238000012549 training Methods 0.000 description 12
- 238000004393 prognosis Methods 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 238000011088 calibration curve Methods 0.000 description 4
- 210000004027 cell Anatomy 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000000611 regression analysis Methods 0.000 description 4
- 101150041972 CDKN2A gene Proteins 0.000 description 2
- 101150012801 Ccsap gene Proteins 0.000 description 2
- 101150050296 Cenpe gene Proteins 0.000 description 2
- 102000009508 Cyclin-Dependent Kinase Inhibitor p16 Human genes 0.000 description 2
- 101100202313 Danio rerio slc39a10 gene Proteins 0.000 description 2
- 101150077503 Fabp5 gene Proteins 0.000 description 2
- 101100202314 Homo sapiens SLC39A10 gene Proteins 0.000 description 2
- 101100202315 Mus musculus Slc39a10 gene Proteins 0.000 description 2
- 101150020537 Tmem97 gene Proteins 0.000 description 2
- 101150041027 ZNRF3 gene Proteins 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 108700042657 p16 Genes Proteins 0.000 description 2
- 101150117250 slc39a10 gene Proteins 0.000 description 2
- 208000024891 symptom Diseases 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 101150071951 ube2t gene Proteins 0.000 description 2
- 208000031404 Chromosome Aberrations Diseases 0.000 description 1
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 108060003951 Immunoglobulin Proteins 0.000 description 1
- 238000010824 Kaplan-Meier survival analysis Methods 0.000 description 1
- 208000002774 Paraproteinemias Diseases 0.000 description 1
- 150000008065 acid anhydrides Chemical class 0.000 description 1
- 238000003705 background correction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002559 cytogenic effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 102000018358 immunoglobulin Human genes 0.000 description 1
- 210000004698 lymphocyte Anatomy 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002493 microarray Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 208000010626 plasma cell neoplasm Diseases 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
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
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
-
- 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
- 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
- G16B45/00—ICT specially adapted for bioinformatics-related data visualisation, e.g. displaying of maps or networks
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biotechnology (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Biology (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Artificial Intelligence (AREA)
- Primary Health Care (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Bioethics (AREA)
- Pathology (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Genetics & Genomics (AREA)
- Molecular Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供了一种基于十基因特征、血清β2‑微球蛋白和LDH的多发性骨髓瘤诺模图构建方法，包括以下步骤：S1，从GEO数据库中获取MM患者的基因表达谱GSE24080，并对所述基因表达谱GSE24080中的基因进行预处理，获得表达值方差最大的前25％的5413个基因；S2，对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块；S3，通过Pearson相关检验将所述功能模块与临床信息之间的相关性进行评估，以确定最显著的模块；S4，使用Cox比例危险度模型对最显著的模块中的基因进行单变量生存分析，并通过LASSO回归筛选出由10个最佳基因组成的评分模型；S5，建立基于所述评分模型、血清β2M和高LDH的新型诺模图来预测新诊断MM患者的3年OS和5年OS。
Description
技术领域
本发明涉及一种基于十基因特征、血清β2-微球蛋白和LDH的多发性骨髓瘤诺模图构建方法。
背景技术
多发性骨髓瘤(MM)是一种恶性浆细胞病，其肿瘤细胞起源于骨髓中的浆细胞，而浆细胞是B淋巴细胞发育到最终功能阶段的细胞。目前，多发性骨髓瘤的异质性日益受到重视，然而，世界上还没有建立有效、简便的MM预后模型。
发明内容
本发明提供了一种基于十基因特征、血清β2-微球蛋白和LDH的多发性骨髓瘤诺模图构建方法，可以有效解决上述问题。
本发明是这样实现的：
一种基于十基因特征、血清β2-微球蛋白和LDH的多发性骨髓瘤诺模图构建方法，包括以下步骤：
S1，从GEO数据库中获取MM患者的基因表达谱GSE24080，并对所述基因表达谱GSE24080中的基因进行预处理，获得表达值方差最大的前 25％的5413个基因；
S2，对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块；
S3，通过Pearson相关检验将所述功能模块与临床信息之间的相关性进行评估，以确定最显著的模块；
S4，使用Cox比例危险度模型对最显著的模块中的基因进行单变量生存分析，并通过LASSO回归筛选出由10个最佳基因组成的评分模型如下：风险评分＝0.239*ZNRF3的表达水平+0.219*UBE2T的表达水平 +0.164*CCSAP的表达水平+0.161*CENPE的表达水平+0.152*PMS2P5的表达水平+0.147*TMEM97的表达水平+0.131*CDKN2A的表达水平 +0.123*SLC39A10的表达水平+0.107*KIF21B的表达水平+0.002*FABP5的表达水平；
S5，建立基于所述评分模型、血清β2M和高LDH的新型诺模图来预测新诊断MM患者的3年OS和5年OS。
本发明的有益效果是：本发明建立一个十基因风险评分、血清β2M和 LDH来构建诺模图。该诺模图能较好地预测个体患者的预后。该模型的C 指数(0.729)优于ISS模型(0.613)。该诺模图可能比R-ISS更有助于定量预测患者的生存率。
附图说明
为了更清楚地说明本发明实施方式的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1A为样品聚类树和临床特征热图。
图1B为分析不同软阈值加权系数的网络拓扑结构图。
图1C为将基于拓扑重叠的不同的基因树状图与分配的模块颜色进行聚类图。
图1D为模块-样本特征关联性热图。
图2为LASSO回归法绘制回归系数图。
图3A为GSE24080训练集多发性骨髓瘤患者的Kaplan-Meier曲线图。
图3B为GSE24080训练集多发性骨髓瘤患者的受试者工作特征曲线 (ROC)图。
图3C为GSE24080训练集多发性骨髓瘤患者的10个预后相关基因的表达谱。
图3D为GSE24080测试集多发性骨髓瘤患者的Kaplan-Meier曲线图。
图3E为GSE24080测试集多发性骨髓瘤患者的受试者工作特征曲线 (ROC)图。
图3F为GSE24080测试集多发性骨髓瘤患者的10个预后相关基因的表达谱。
图3G为GSE24080全集多发性骨髓瘤患者的Kaplan-Meier曲线图。
图3H为GSE24080全集多发性骨髓瘤患者的受试者工作特征曲线 (ROC)图。
图3I为GSE24080全集多发性骨髓瘤患者的10个预后相关基因的表达谱。
图3J为GSE57317数据集多发性骨髓瘤患者的Kaplan-Meier曲线图。
图3K为GSE57317数据集多发性骨髓瘤患者的受试者工作特征曲线 (ROC)图。
图3L为GSE57317数据多发性骨髓瘤患者的10个预后相关基因的表达谱。
图4A为通过综合风险评分系统(IRSS)获得的本发明训练集中多发性骨髓瘤患者的Kaplan-Meier曲线图。
图4B为通过国际分期系统(ISS)获得的训练集中多发性骨髓瘤患者的Kaplan-Meier曲线图。
图4C为通过综合风险评分系统(IRSS)获得的本发明测试集中多发性骨髓瘤患者的Kaplan-Meier曲线图。
图4D为通过国际分期系统(ISS)获得的测试集中多发性骨髓瘤患者的Kaplan-Meier曲线图。
图4E为通过综合风险评分系统(IRSS)获得的本发明全集中多发性骨髓瘤患者的Kaplan-Meier曲线图。
图4F为通过国际分期系统(ISS)获得的全集中多发性骨髓瘤患者的 Kaplan-Meier曲线图。
图5A为本发明提供的多发性骨髓瘤MM预后诺模图。
图5B为本发明用于预测GSE24080全集患者3年OS的校准曲线。
图5C为本发明用于预测GSE24080全集患者5年OS的校准曲线。
图5D为本发明用于预测GSE24080验证集患者3年OS的校准曲线。
图5E为本发明用于预测GSE24080测试集患者3年OS的校准曲线。
具体实施方式
为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。因此，以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
本发明实施例提供一种基于十基因特征、血清β2-微球蛋白(β2M)和乳酸脱氢酶(LDH)的多发性骨髓瘤诺模图构建方法，其包括以下步骤：
S1，从GEO数据库中获取MM患者的基因表达谱GSE24080，并对所述基因表达谱GSE24080中的基因进行预处理，获得表达值方差最大的前 25％的5413个基因；
S2，对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块；
S3，通过Pearson相关检验将所述功能模块与临床信息之间的相关性进行评估，以确定最显著的模块；
S4，使用Cox比例危险度模型对最显著的模块中的基因进行单变量生存分析，并通过LASSO回归筛选出由10个最佳基因组成的十基因评分模型：风险评分＝0.239*ZNRF3的表达水平+0.219*UBE2T的表达水平 +0.164*CCSAP的表达水平+0.161*CENPE的表达水平+0.152*PMS2P5的表达水平+0.147*TMEM97的表达水平+0.131*CDKN2A的表达水平+0.123*SLC39A10的表达水平+0.107*KIF21B的表达水平+0.002*FABP5的表达水平；
S5，将所述十基因评分模型或血清β2M或LDH高于截止值时各因子得分为1分，否则为0分为规则建立了综合风险评分系统。
在步骤S1中，在GSE24080数据集中，总共包括了549个来自新诊断 MM患者的样本。本发明中以340份纳入整体治疗2(TT2)试验中的样本作为训练集；其余219份样本中包括纳入整体治疗3(TT3)试验的214份样本和纳入微阵列质量控制的5份样本为测试集。训练集人群的病人特征 (包括年龄、性别、IgA、β2M、C反应蛋白、肌酐、LDH、血红蛋白、骨髓浆细胞和细胞抗体)与测试集人群相似。
为了WGCNA基因共表达网络分析，对来自GSE24080的原始数据进行相同的预处理，以便进行背景校正和归一化，从GSE24080的原始数据共鉴定出21653个基因。进一步的，对所述基因表达谱GSE24080中的21653 个基因进行预处理，获得表达值方差最大的前25％的5413个基因，并以5413 个基因为输入进行聚类分析，如图1A所示。图1A中，颜色强度与年龄、高β2微球蛋白(β2M)、C反应蛋白(CRP)、肌酐(CREAT)、乳酸脱氢酶 (LDH)、血红蛋白(HGB)、骨髓浆细胞(BMPC)及较长生存月成正比。白色代表女性、IgA类型和没有细胞遗传学异常样本。在生存状态下，白色意味着病人活着，红色意味着病人死亡。
在步骤S2中，所述对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块的步骤包括：
S21，选择无尺度拓扑指数达到0.85的最低加权系数8来生成分层聚类树，如图1B所示。图1B左侧图显示软阈值加权系数(x轴)对无尺度拟合指数(y轴)的影响。图1B右侧图显示了软阈值加权系数(x轴)对平均连通度(度，y轴)的影响。
S22，结合拓扑重叠矩阵，采用层次平均连锁聚类方法检测每个基因网络的基因模块。
S23，使用Dynamic Tree Cut(deepSp.＝2)识别20个共表达的功能模块，如图1C所示。其中，图1C每个着色分支指示高度连接基因的模块。
在步骤S3中，所述临床信息包括年龄、性别、免疫球蛋白IgA类型、血清β2-微球蛋白、C反应蛋白、肌酐、乳酸脱氢酶、血红蛋白、骨髓浆细胞、细胞遗传学异常、无事件生存时间和状态、生存期时间和状态。另外，请参照图1D，每行对应一个模块特征基因，每列对应一个临床特征。每个单元格包含第一行中的对应相关性和第二行中的P值。所述最显著的模块的包括：聚集于黑色模块中的240个基因，其分别与EFS时间和OS时间性状的相关性最强。根据相关系数，我们发现，聚集于黑色模块中的基因 (240个基因)与EFS时间和OS时间性状的相关性最强(Pearson的 R2＝-0.25，P-value＝4e-6和Pearson的R2＝-0.25，P-value＝3e-6；)。根据显著性检验的P值，黑色模块与MM生存率、HGB呈负相关，与MM相关预后指标β2M、CRP、CREAT、LDH及CytoAbn呈正相关。因此，黑色模块被认为是一个有代表性的生存相关模块。
在步骤S4中，请参照图2及表1，对182个与MM患者OS不良显著相关的基因进行LASSO回归分析
表1.LASSO回归分析与训练集340名患者总体生存有关的基因
计算训练集各样本的风险评分，以中位数(9.426)为阈值。340例患者分为高危(n＝170)和低危(n＝170)。Kaplan-Meier生存分析显示高危和低危患者的生存率有显著性差异(HR＝3.068,95％CI为2.089-4.505,Log-rank检验 P<0.001,图3A所示)。图3D及3G分别表示测试集以及全集。高危患者中位OS为69.0个月，低危患者未达到中位OS。此外，根据时间相关ROC 分析得出的预后模型3年生存率预测的AUC为0.749，明显高于已发表的基因模型EMC92(AUC＝0.71)、UAMS-70(AUC＝0.737)和UAMS-17 (AUC＝0.717)，表明该预测模型可很好预测MM患者OS的能力(图3B 所示，以及图3E及3H)。图3C显示了训练集中10个预后基因的表达谱，结果表明高危组中10个基因的表达水平高于低危组。图3F及3I显示了测试集以及全集10个预后基因的表达谱，结果表明高危组中10个基因的表达水平高于低危组。
在GSE57317外部验证数据集中，十基因模型可将患者分为OS明显不同的高危组(n＝34)的和低危组(n＝21)。高危组的OS显著短于低危组 (HR＝8.445，95％CI 1.088-65.581，P＝0.041，图3J)。ROC曲线预测3年OS 的AUC达到0.859(图3K)，与EMC92、UAMS-70(AUC＝0.737)和UAMS-17 相比也具有可比性。图3L显示了GSE57317数据集十个预后基因的表达谱。
在步骤S5中，通过单因素Cox回归分析，十基因评分模型和年龄、β2M、肌酐、LDH、HGB、BMPC和CytoAbn的临床协变量对预后有一定的预测价值(如表2所示)。我们发现，由十基因评分模型计算的HR高于任何一个临床协变量，表明其较高的预测效率。将年龄、β2M、肌酐、LDH、HGB、 BMPC、CytoAbn与十基因评分模型进行多元Cox回归分析，结果表明十基因评分模型、β2M和CytoAbn是OS的独立预后因子(如表2所示)。
表2.每个数据集多发性骨髓瘤患者总体生存率的单因素和多因素分析
在新的综合风险评分系统(IRSS)中，十基因评分模型或β2M或LDH 高于截止值得分为1分，否则为0。患者被分成三组：低危，评分为0(无因素)；中危，评分为1-2(三分之一或两个因素)；高危，评分为3(所有三个因素)。在训练集中分出以下三组(n＝340)：101例(30％)患者为低危组；190例(56％)患者为中危组；49例(14％)患者为高危组。如图 4A所示，5年期OS分别为85.73％、64.42％和34.0％。低危组和中危组未达到中位OS，而高危组中位OS为43个月。中危组与低危组相比，死亡的风险更高(HR＝2.852，95％CI为1.683-4.833，P<0.001)，高危组与中危组相比死亡风险也更高(HR＝2.349，95％CI为1.562-3.531，P<0.001)。如图 4B所示，ISSⅠ、Ⅱ和Ⅲ期患者的5年OS率分别为76.04％、59.87％和 47.36％。ISSⅠ期和ISSⅡ期未达到中位OS，ISSⅢ期中位OS为43.5个月。因此，与ISS相比，这种新的IRSS可以更清楚地将患者分为三个不同的风险组。
然后将这种新的IRSS应用于GSE24080的测试集和全集。测试集或全集的患者被分为三个危险组，也可以观察到类似的模式(图4C-4F)。正如所预期的，对于中危组与低危组相比以及高危组与中危组相比，其死亡的风险均增加。然而，测试集的结果表明，ISS II期与ISS I期相比，死亡的风险没有显著增加(HR＝1.449，95％CI为0.664-3.163，P＝0.352)。显然， IRSS提供了比ISS分类更高的预测精度。
在步骤S4之后，可以进一步包括：
S6，建立基于所述十基因评分模型、血清β2M和高LDH的诺模图来预测新诊断MM患者的3年OS和5年OS。
图5A显示了诺模图中MM的3年和5年OS的预测。结合β2M和LDH 的十基因评分模型的C指数(0.729；95％CI 0.649～0.809；P<0.001)优于ISS 模型(0.613；95％CI 0.537～0.701；P<0.01)，也优于EMC92模型(0.653； 95％CI,0.557to 0.749；P<0.01),UAMS-70模型(0.665；95％CI,0.573to 0.757；P<0.01)和UAMS-17模型(0.666；95％CI,0.576to 0.756；P<0.01)。请参照图5B-E，与训练集一致，测试集和整个集合的诺模图的C指数分别为0.72和0.754。3或5年生存率的校准图与实际观测值也有较好的相关性。
本发明实施例进一步提供一种用于检测多发性骨髓瘤的基因探针组合物，包括ZNRF3基因探针、UBE2T基因探针、CCSAP基因探针、CENPE 基因探针、PMS2P5、TMEM97基因探针、CDKN2A基因探针、SLC39A10 基因探针、KIF21B基因探针、以及FABP5基因探针。
本发明另一实施例进一步提供一种用于检测多发性骨髓瘤的基因探针组合物，由ZNRF3基因探针、UBE2T基因探针、CCSAP基因探针、CENPE 基因探针、PMS2P5、TMEM97基因探针、CDKN2A基因探针、SLC39A10 基因探针、KIF21B基因探针、以及FABP5基因探针组成。
本发明实施例还进一步提供一种用于检测多发性骨髓瘤的试剂盒，含有上述的基因探针组合物。
以上所述仅为本发明的优选实施方式而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于十基因特征、血清β2-微球蛋白和LDH的多发性骨髓瘤诺模图构建方法，其特征在于，包括以下步骤：
S1，从GEO数据库中获取MM患者的基因表达谱GSE24080，并对所述基因表达谱GSE24080中的基因进行预处理，获得表达值方差最大的前25％的5413个基因；
S2，对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块；
S3，通过Pearson相关检验将所述功能模块与临床信息之间的相关性进行评估，以确定最显著的模块；
S4，使用Cox比例危险度模型对最显著的模块中的基因进行单变量生存分析，并通过LASSO回归筛选出由10个最佳基因组成的评分模型如下：风险评分＝0.239*ZNRF3的表达水平+0.219*UBE2T的表达水平+0.164*CCSAP的表达水平+0.161*CENPE的表达水平+0.152*PMS2P5的表达水平+0.147*TMEM97的表达水平+0.131*CDKN2A的表达水平+0.123*SLC39A10的表达水平+0.107*KIF21B的表达水平+0.002*FABP5的表达水平；
S5，建立基于所述评分模型、血清β2M和高LDH的新型诺模图来预测新诊断MM患者的3年OS和5年OS。
2.根据权利要求1所述的构建方法，其特征在于，在步骤S2中，所述对所述5413个基因进行WGCNA基因共表达网络分析，以识别共表达的功能模块的步骤包括：
S21，选择无尺度拓扑指数达到0.85的最低加权系数8来生成分层聚类树；
S22，结合拓扑重叠矩，采用层次平均连锁聚类方法检测每个基因网络的基因模块；
S23，使用Dynamic Tree Cut(deepSp.＝2)识别20个共表达的功能模块。
3.根据权利要求1所述的构建方法，其特征在于，所述临床信息包括年龄、性别、免疫球蛋白IgA类型、血清β2-微球蛋白、C反应蛋白、肌酐、乳酸脱氢酶、血红蛋白、骨髓浆细胞、细胞遗传学异常、无事件生存时间和状态、生存期时间和状态。
4.根据权利要求1所述的构建方法，其特征在于，在步骤S3中，所述最显著的模块的包括：聚集于黑色模块中的240个基因，其分别与EFS时间和OS时间性状的相关性最强。
5.根据权利要求1所述的构建方法，其特征在于，所述新型诺模图的C指数为：0.729；95％CI0.649～0.809；P<0.001。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910325039.3A CN110197701B (zh) | 2019-04-22 | 2019-04-22 | 一种新型多发性骨髓瘤诺模图构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910325039.3A CN110197701B (zh) | 2019-04-22 | 2019-04-22 | 一种新型多发性骨髓瘤诺模图构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110197701A true CN110197701A (zh) | 2019-09-03 |
| CN110197701B CN110197701B (zh) | 2021-08-10 |
Family
ID=67751980
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910325039.3A Active CN110197701B (zh) | 2019-04-22 | 2019-04-22 | 一种新型多发性骨髓瘤诺模图构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110197701B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101392285A (zh) * | 2006-09-21 | 2009-03-25 | 富士胶片株式会社 | 多发性骨髓瘤的检测方法及抑制方法 |
| WO2010111208A1 (en) * | 2009-03-23 | 2010-09-30 | University Of Miami | Mitochondrial inhibitors and uses thereof |
| CN101932724A (zh) * | 2007-10-05 | 2010-12-29 | 环太平洋生物技术有限公司 | 胃肠癌的增生标签及预后 |
| CN102186987A (zh) * | 2008-04-24 | 2011-09-14 | 阿肯色大学托管委员会 | 一种基于基因表达谱的用来标示高风险多发性骨髓瘤的基因组指纹的方法和用途 |
| CN107463796A (zh) * | 2017-07-12 | 2017-12-12 | 北京航空航天大学 | 基于基因共表达网络传播分析的早期致病因子探测方法 |
| CN108470111A (zh) * | 2018-05-09 | 2018-08-31 | 中国科学院昆明动物研究所 | 一种基于多基因表达特征谱的胃癌个性化预后评估方法 |
| CN108733980A (zh) * | 2018-05-09 | 2018-11-02 | 中国科学院昆明动物研究所 | 一种基于多基因表达特征谱的胶质母细胞瘤个性化预后评估方法 |
-
2019
- 2019-04-22 CN CN201910325039.3A patent/CN110197701B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101392285A (zh) * | 2006-09-21 | 2009-03-25 | 富士胶片株式会社 | 多发性骨髓瘤的检测方法及抑制方法 |
| CN101932724A (zh) * | 2007-10-05 | 2010-12-29 | 环太平洋生物技术有限公司 | 胃肠癌的增生标签及预后 |
| CN102186987A (zh) * | 2008-04-24 | 2011-09-14 | 阿肯色大学托管委员会 | 一种基于基因表达谱的用来标示高风险多发性骨髓瘤的基因组指纹的方法和用途 |
| WO2010111208A1 (en) * | 2009-03-23 | 2010-09-30 | University Of Miami | Mitochondrial inhibitors and uses thereof |
| CN107463796A (zh) * | 2017-07-12 | 2017-12-12 | 北京航空航天大学 | 基于基因共表达网络传播分析的早期致病因子探测方法 |
| CN108470111A (zh) * | 2018-05-09 | 2018-08-31 | 中国科学院昆明动物研究所 | 一种基于多基因表达特征谱的胃癌个性化预后评估方法 |
| CN108733980A (zh) * | 2018-05-09 | 2018-11-02 | 中国科学院昆明动物研究所 | 一种基于多基因表达特征谱的胶质母细胞瘤个性化预后评估方法 |
Non-Patent Citations (4)
| Title |
|---|
| PETER LANGFELDER 等: "WGCNA: an R package for weighted correlation network analysis", 《BMC BIOINFORMATICS》 * |
| XIAO-PING LIU 等: "Development and Validation of a 9-Gene Prognostic Signature in Patients With Multiple Myeloma", 《ORIGINAL RESEARCH》 * |
| 徐声鸣 等: "多发性骨髓瘤患者血清lncRNA HOTAIR表达量与血清β2微球蛋白及预后的相关性", 《中国实验诊断学》 * |
| 曲莹: "血清microRNA、长链非编码RNA及CRBN蛋白表达在多发性骨髓瘤中的临床意义", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110197701B (zh) | 2021-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110232974A (zh) | 一种新型多发性骨髓瘤综合风险评分方法 | |
| Adhikari et al. | A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia | |
| Holden et al. | Rapid radiation, borrowing and dialect continua in the Bantu languages | |
| CN109478231A (zh) | 帮助区别良性和恶性放射线照相明显肺结节的方法和组合物 | |
| CN105849279A (zh) | 用于识别疾病诱导的突变的方法和系统 | |
| CN107422126A (zh) | 心血管危险事件预测及其用途 | |
| JP2014521926A (ja) | 細胞におけるバイオマーカーの発現の積率による解析 | |
| CN111524594A (zh) | 目标人群血液系统恶性肿瘤筛查系统 | |
| CN113327679A (zh) | 一种肺栓塞临床风险及预后评分方法与系统 | |
| CN113223714B (zh) | 一种用于预测子痫前期风险的基因组合、子痫前期风险预测模型及其构建方法 | |
| JP2018517192A (ja) | 生物学的経路の関連性の評価を行うコンピュータ実装方法、コンピュータ・システム、およびコンピュータ・プログラム | |
| CN116287204A (zh) | 检测特征基因的突变情况在制备静脉血栓栓塞症风险检测产品中的应用 | |
| CN109929934A (zh) | 免疫相关基因在结直肠癌预后的试剂盒和系统中的应用 | |
| CN107169264B (zh) | 一种复杂疾病诊断系统 | |
| Li et al. | Current progress and future prospects in phylofloristics | |
| CN110223733A (zh) | 一种新型多发性骨髓瘤预后基因的筛查方法 | |
| KR102467999B1 (ko) | 위암의 다층 다요인 패널과 Computational biological network modeling을 통한 위암 발암에 대한 에티옴 모형 | |
| CN106446601B (zh) | 一种大规模标注lncRNA功能的方法 | |
| CN110197701A (zh) | 一种新型多发性骨髓瘤诺模图构建方法 | |
| CN110218789A (zh) | 一种新型多发性骨髓瘤的基因探针检测组合物及试剂盒 | |
| CN109811057A (zh) | 缺氧相关基因在结直肠癌预测系统中的应用 | |
| CN114678062B (zh) | 基于多组学特征的肝细胞癌预后预测系统及其预测方法 | |
| CN106650284B (zh) | 一种疾病康复评价系统 | |
| CN115033758A (zh) | 肾透明细胞癌预后标记基因的应用及筛选、预后预测方法 | |
| Avila Santos et al. | BioDeepfuse: a hybrid deep learning approach with integrated feature extraction techniques for enhanced non-coding RNA classification |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220608 Address after: 1st and 2nd floors, No.21 Lane 908, Ziping Road, Pudong New Area, Shanghai, 200120 Patentee after: SHANGHAI TISSUEBANK MEDICAL LABORATORY Co.,Ltd. Address before: No.20, Chazhong Road, Taijiang District, Fuzhou City, Fujian Province Patentee before: THE FIRST AFFILIATED HOSPITAL OF FUJIAN MEDICAL University |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A new method for constructing nomographs of Multiple myeloma Effective date of registration: 20230719 Granted publication date: 20210810 Pledgee: Industrial Bank Co.,Ltd. Shanghai Nanhui Branch Pledgor: SHANGHAI TISSUEBANK MEDICAL LABORATORY Co.,Ltd.|SHENZHEN TISSUEBANK PRECISION MEDICINE CO.,LTD.|SHANGHAI TISSUEBANK BIOTECHNOLOGY Co.,Ltd.|Shanghai dishuobeiken Gene Technology Co.,Ltd.|SHANGHAI TISSUEBANK BIOTECHNOLOGY CO.,LTD. Registration number: Y2023310000384 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Granted publication date: 20210810 Pledgee: Industrial Bank Co.,Ltd. Shanghai Nanhui Branch Pledgor: SHANGHAI TISSUEBANK MEDICAL LABORATORY Co.,Ltd.|SHENZHEN TISSUEBANK PRECISION MEDICINE CO.,LTD.|SHANGHAI TISSUEBANK BIOTECHNOLOGY Co.,Ltd.|Shanghai dishuobeiken Gene Technology Co.,Ltd.|SHANGHAI TISSUEBANK BIOTECHNOLOGY CO.,LTD. Registration number: Y2023310000384 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A novel nomogram construction method for multiple myeloma Granted publication date: 20210810 Pledgee: Industrial Bank Co.,Ltd. Shanghai Nanhui Branch Pledgor: SHANGHAI TISSUEBANK MEDICAL LABORATORY Co.,Ltd.|SHENZHEN TISSUEBANK PRECISION MEDICINE CO.,LTD.|SHANGHAI TISSUEBANK BIOTECHNOLOGY Co.,Ltd.|Shanghai dishuobeiken Gene Technology Co.,Ltd.|SHANGHAI TISSUEBANK BIOTECHNOLOGY CO.,LTD. Registration number: Y2024310000533 |