CN114242240A - 一种分化型甲状腺癌患者疾病持续复发预测模型 - Google Patents
一种分化型甲状腺癌患者疾病持续复发预测模型 Download PDFInfo
- Publication number
- CN114242240A CN114242240A CN202111182720.0A CN202111182720A CN114242240A CN 114242240 A CN114242240 A CN 114242240A CN 202111182720 A CN202111182720 A CN 202111182720A CN 114242240 A CN114242240 A CN 114242240A
- Authority
- CN
- China
- Prior art keywords
- scale
- patient
- thyroid cancer
- corresponding score
- prediction
- 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.)
- Pending
Links
- 201000010099 disease Diseases 0.000 title claims abstract description 64
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 64
- 208000015799 differentiated thyroid carcinoma Diseases 0.000 title claims abstract description 37
- 210000001165 lymph node Anatomy 0.000 claims abstract description 52
- 230000002688 persistence Effects 0.000 claims abstract description 41
- 230000001394 metastastic effect Effects 0.000 claims abstract description 31
- 206010061289 metastatic neoplasm Diseases 0.000 claims abstract description 31
- 206010028980 Neoplasm Diseases 0.000 claims abstract description 25
- 230000009545 invasion Effects 0.000 claims abstract description 18
- 210000001685 thyroid gland Anatomy 0.000 claims abstract description 18
- 238000004393 prognosis Methods 0.000 claims abstract description 14
- 238000002271 resection Methods 0.000 claims abstract description 11
- 208000007433 Lymphatic Metastasis Diseases 0.000 claims abstract description 10
- 238000010276 construction Methods 0.000 claims abstract description 5
- 208000024770 Thyroid neoplasm Diseases 0.000 claims description 26
- 201000002510 thyroid cancer Diseases 0.000 claims description 23
- 201000010198 papillary carcinoma Diseases 0.000 claims description 6
- 201000009030 Carcinoma Diseases 0.000 claims description 3
- 230000002327 eosinophilic effect Effects 0.000 claims description 3
- 230000003325 follicular Effects 0.000 claims description 3
- 208000030901 thyroid gland follicular carcinoma Diseases 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims 1
- 206010061819 Disease recurrence Diseases 0.000 abstract description 4
- 238000011337 individualized treatment Methods 0.000 abstract description 2
- 238000013517 stratification Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 238000012549 training Methods 0.000 description 7
- 206010027476 Metastases Diseases 0.000 description 6
- 238000011088 calibration curve Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 6
- 230000004083 survival effect Effects 0.000 description 6
- 238000000034 method Methods 0.000 description 5
- 238000012216 screening Methods 0.000 description 5
- 238000011282 treatment Methods 0.000 description 5
- 108010034949 Thyroglobulin Proteins 0.000 description 4
- 102000009843 Thyroglobulin Human genes 0.000 description 4
- 230000003902 lesion Effects 0.000 description 4
- 230000007170 pathology Effects 0.000 description 4
- 230000002085 persistent effect Effects 0.000 description 4
- 238000001356 surgical procedure Methods 0.000 description 4
- 229960002175 thyroglobulin Drugs 0.000 description 4
- 102000011923 Thyrotropin Human genes 0.000 description 3
- 108010061174 Thyrotropin Proteins 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 238000002790 cross-validation Methods 0.000 description 3
- 238000002224 dissection Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000005764 inhibitory process Effects 0.000 description 3
- 230000002980 postoperative effect Effects 0.000 description 3
- 230000000306 recurrent effect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 230000009401 metastasis Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 208000030836 Hashimoto thyroiditis Diseases 0.000 description 1
- 101000984753 Homo sapiens Serine/threonine-protein kinase B-raf Proteins 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 206010020850 Hyperthyroidism Diseases 0.000 description 1
- 210000005131 Hürthle cell Anatomy 0.000 description 1
- 102100027103 Serine/threonine-protein kinase B-raf Human genes 0.000 description 1
- 208000009453 Thyroid Nodule Diseases 0.000 description 1
- 230000001093 anti-cancer Effects 0.000 description 1
- 238000007475 c-index Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 210000000750 endocrine system Anatomy 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 208000017333 follicular variant thyroid gland papillary carcinoma Diseases 0.000 description 1
- 208000014829 head and neck neoplasm Diseases 0.000 description 1
- 238000011221 initial treatment Methods 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 229940036571 iodine therapy Drugs 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 238000004223 overdiagnosis Methods 0.000 description 1
- 210000002990 parathyroid gland Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000391 smoking effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002459 sustained effect Effects 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
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- 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/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- 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/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- 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
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供了一种分化型甲状腺癌患者疾病持续复发预测模型，属于预测模型领域。本发明基于患者的肿瘤大小、淋巴结转移比例、转移淋巴结最大径、淋巴结结外侵犯、组织学类型多指标构建了一种预测分化型甲状腺癌患者预后的预测模型。利用该预测模型能够准确预测行甲状腺全切后的分化型甲状腺癌患者的预后情况，准确预测患者出现疾病持续或复发的概率，有助于针对每个患者做出个体化的预后预测。本发明提供的预测系统构建方法简单，预测准确度和区分度高，对临床上辅助判断行甲状腺全切后的分化型甲状腺癌患者的疾病持续或复发状况具有重要意义，有助于指导医师制定个体化治疗决策。
Description
技术领域
本发明属于预测模型领域，具体涉及一种分化型甲状腺癌患者疾病持续复发预测模型。
背景技术
甲状腺癌(Thyroid cancer,TC)是内分泌系统最常见的恶性肿瘤，其发病率在全球均逐渐上升，我国甲状腺癌的发病率每年增长20.1％。甲状腺癌发病率增长的一个重要原因是诊断手段的提高，尤其是高分辨超声的使用，使许多甲状腺癌被早期发现并诊断。然而，TC病死率却并未伴随发病率增长而有明显变化，这就导致了对于甲状腺癌是否存在过度诊断与过度治疗的广泛争议，因而开展甲状腺癌相关问题研究具有非常重要的意义。
分化型甲状腺癌(Differentiated thyroid cancer,DTC)是甲状腺癌中最常见的类型，占甲状腺癌的90％。近年来分化型甲状腺癌由于发病率逐年增高而引发了人们的关注。DTC的治疗策略主要包括手术、术后选择性放射性碘-131(radioiodine-131，RAI-131)治疗及促甲状腺激素(thyroid-stimulating hormone，TSH)抑制治疗。经过上述规范治疗的DTC患者10年生存率达90％以上。但是，仍有约30％的患者会出现复发，提示DTC患者的复发风险不容小觑。为了评估患者的术后复发风险，2009年美国甲状腺协会(AmericanThyroid Association，ATA)指南及2012年我国中华医学会内分泌学分会、外科分会内分泌学组、核医学分会和中国抗癌协会头颈肿瘤专业委员会联合制定的《甲状腺结节和分化型甲状腺癌诊治指南》将DTC患者术后进行高、中、低危分层。但是，研究发现，该复发风险分层系统对行甲状腺全切术后的分化型甲状腺癌患者的疾病持续或复发情况的预测效果有待进一步提高。
为了指导医师对分化型甲状腺癌患者制定个体化治疗决策，开发出一种能够准确预测分化型甲状腺癌患者行甲状腺全切后的未发生疾病、疾病持续或疾病复发概率的预测系统具有重要意义。
发明内容
本发明的目的在于提供一种能够准确预测分化型甲状腺癌患者行甲状腺全切后疾病持续或疾病复发概率的预测系统。
本发明提供了一种预测甲状腺癌患者预后的预测系统，所述预测系统是以患者的肿瘤大小、淋巴结转移比例、转移淋巴结最大直径、淋巴结结外侵犯情况、组织学类型为预测指标构建的。
进一步地，所述预测系统是列线图，该列线图包含第1～12条直线，第1～12条直线由上至下依次排列，且相互平行；每条直线代表一个标尺，标尺上有刻度；
第1个标尺表示第2～6个标尺上刻度对应的分值的标尺；
第2个标尺表示肿瘤大小；
第3个标尺表示转移淋巴结比例；
第4个标尺表示转移淋巴结最大直径；
第5个标尺表示淋巴结结外侵犯情况；
第6个标尺表示组织学类型；
第7个标尺表示第2～6个标尺上刻度对应的分值的总和的标尺；
第8个标尺表示患者6个月内未出现疾病持续或复发的概率；
第9个标尺表示患者1年内未出现疾病持续或复发的概率；
第10个标尺表示患者2年内未出现疾病持续或复发的概率；
第11个标尺表示患者3年内未出现疾病持续或复发的概率；
第12个标尺表示患者5年内未出现疾病持续或复发的概率。
进一步地，第1个标尺的刻度范围为1～100，1在最左端，100在最右端；
第2个标尺中，若患者肿瘤大小≤10mm，在第1个标尺上对应的分值为0；若患者肿瘤大小>10且≤20mm，在第1个标尺上对应的分值为21；若患者肿瘤大小>20且≤40mm，在第1个标尺上对应的分值为42；若患者肿瘤大小>40mm，在第1个标尺上对应的分值为63；
第3个标尺中，若患者转移淋巴结比例≤30％，在第1个标尺上对应的分值为0；若患者转移淋巴结比例>30且≤50％，在第1个标尺上对应的分值为17；若患者转移淋巴结比例>50％且≤80％，在第1个标尺上对应的分值为32；若患者转移淋巴结比例>80％，在第1个标尺上对应的分值为48；
第4个标尺中，若患者无淋巴结转移，在第1个标尺上对应的分值为0；若患者转移淋巴结最大直径<30mm，在第1个标尺上对应的分值为9；若患者转移淋巴结最大直径≥30mm，在第1个标尺上对应的分值为18；
第5个标尺中，若患者无淋巴结结外侵犯，在第1个标尺上对应的分值为0；若患者有淋巴结结外侵犯，在第1个标尺上对应的分值为42；
第6个标尺中，若患者为乳头状癌，在第1个标尺上对应的分值为0；
若患者为滤泡亚型乳头状癌，在第1个标尺上对应的分值为32；若患者为滤泡癌，在第1个标尺上对应的分值为66；若患者为嗜酸细胞癌，在第1个标尺上对应的分值为100；
第7个标尺的刻度范围为0～180，0在最左端，180在最右端；
第8个标尺中，患者6个月内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第8个标尺上对应的刻度值，范围为0.2～0.95；
第9个标尺中，患者1年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第9个标尺上对应的刻度值，范围为0.05～0.95；
第10个标尺中，患者2年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第10个标尺上对应的刻度值，范围为0.05～0.95；
第11个标尺中，患者3年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第11个标尺上对应的刻度值，范围为0.05～0.95；
第12个标尺中，患者5年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第12个标尺上对应的刻度值，范围为0.05～0.95。
进一步地，所述甲状腺癌患者为分化型甲状腺癌患者。
进一步地，所述分化型甲状腺癌患者为行甲状腺全切术后的患者。
进一步地，所述列线图的构建方法包括以下步骤：
(1)收集甲状腺癌患者的预测指标，录入输入模块；
(2)利用输入模块中的预测指标构建多因素Cox回归模型，绘制列线图。
进一步地，所述列线图如图2所示。
本发明还提供了一种预测甲状腺癌患者预后的设备，所述设备包括上述的预测系统。
本发明还提供了上述的预测系统在制备预测甲状腺癌患者预后的设备中的用途。
进一步地，所述甲状腺癌患者为分化型甲状腺癌患者；优选地，所述分化型甲状腺癌患者为行甲状腺全切术后的患者。
本发明中，原发灶大小即肿瘤大小，指病理报告的肿瘤原发灶最大径；淋巴结转移比例指病理报告的转移阳性淋巴结数占手术清扫淋巴结数的比例；转移淋巴结最大直径指病理报告的转移阳性淋巴结的最大径，淋巴结结外侵犯指病理报告的转移阳性淋巴结外的侵犯，组织学类型指病理报告的肿瘤组织学类型。
对于接受了碘治疗的分化型甲状腺癌患者，其无疾病、疾病持续、疾病复发的定义如下：(1)无疾病状态：没有疾病存在的证据，影像学检查阴性，甲状腺球蛋白抗体(TgAb)阴性时抑制状态下的甲状腺球蛋白(Tg)<0.2ug/L或刺激状态下Tg<1ug/L。(2)疾病持续：指生化持续或结构持续，生化持续指无分化型甲状腺癌病灶存在的影像学及病理学证据，但TgAb阴性的抑制状态下的Tg≥0.2ug/L或刺激状态下Tg≥1ug/L，或TgAb阳性；结构持续指碘治疗后6-12月随访时，结构或功能影像学检查发现病灶，并通过病理学证实。(3)疾病复发：指分化型甲状腺癌患者达到无疾病状态半年后出现疾病持续状态，包括生化复发和结构复发。
对于只接受甲状腺全切术，未接受碘治疗的分化型甲状腺癌患者，其无疾病、疾病持续、疾病复发的定义如下：(1)无疾病状态：没有疾病存在的证据，影像学检查阴性，甲状腺球蛋白抗体(TgAb)阴性时甲状腺球蛋白(Tg)<0.2ug/L。(2)疾病持续：指生化持续或结构持续，生化持续指无分化型甲状腺癌病灶存在的影像学及病理学证据，但TgAb阴性的抑制状态下的Tg≥5ug/L，或Tg在TSH水平可比情况下持续升高，TgAb不明原因持续升高；结构持续指结构或功能影像学检查发现病灶，并通过病理学证实。(3)疾病复发：指分化型甲状腺癌患者达到无疾病状态半年后出现疾病持续状态，包括生化复发和结构复发。疾病持续或复发是观察分化型甲状腺癌患者术后疗效的一个重要指标，本发明利用区域性医疗中心的真实数据，基于患者的肿瘤大小、淋巴结转移比例、转移淋巴结最大径、淋巴结结外侵犯、组织学类型多指标构建了一种预测分化型甲状腺癌患者预后的预测模型。利用该预测模型能够准确预测行甲状腺全切后的分化型甲状腺癌患者的预后情况，准确预测患者出现疾病持续或复发的概率，有助于针对每个患者做出个体化的预后预测。
本发明提供的预测系统构建方法简单，预测准确度和区分度高，对临床上辅助判断行甲状腺全切后的分化型甲状腺癌患者的疾病持续或复发状况具有重要意义，有助于指导医师制定个体化治疗决策。
显然，根据本发明的上述内容，按照本领域的普通技术知识和惯用手段，在不脱离本发明上述基本技术思想前提下，还可以做出其它多种形式的修改、替换或变更。
以下通过实施例形式的具体实施方式，对本发明的上述内容再作进一步的详细说明。但不应将此理解为本发明上述主题的范围仅限于以下的实例。凡基于本发明上述内容所实现的技术均属于本发明的范围。
附图说明
图1.LASSO筛选结果。上图：特征系数惩罚图；下图：LASSO逻辑回归模型进行特征选择，选择在距离均方误差最小时一个标准误的λ值(λ＝0.05)。图2.预测分化型甲状腺癌患者无疾病持续/复发生存期的列线图。
图3.模型列线图的校准曲线(从左至右依次为：预测1年、3年和5年的无疾病持续/复发生存期的校准曲线)。
具体实施方式
本发明所用原料与设备均为已知产品，通过购买市售产品所得。
实施例1建立预测分化型甲状腺癌患者在行甲状腺全切术后疾病持续或复发的预测系统
本发明的预测系统为预测分化型甲状腺癌患者在初始治疗后的未出现疾病持续或复发的概率的列线图，构建方法包括以下步骤。
(1)从四川大学华西医院的甲状腺恶性肿瘤数据库中筛选出2009-2019年间于四川大学华西医院就诊并确诊为分化型甲状腺癌的患者3471名，所有患者在确诊后接受甲状腺全切术(伴或不伴颈淋巴结清扫术)。经过中位时间为44.23个月的随访后，最终有418名患者发生疾病的持续或复发。
(2)结合文献复习与专家咨询，形成包括3类特征的特征池，包括(a)人口学特征：年龄、性别、民族、饮酒史、吸烟史、体质量、合并症_糖尿病、合并症_高血压、合并症_桥本氏甲状腺炎、合并症_甲状腺功能亢进。(b)手术相关特征结果：甲状腺切除术、中央区淋巴结清扫术、侧区淋巴结清扫术、淋巴结清扫术、甲状旁腺移植术。(c)术前检查特征结果：结节大小、结节位置、结节数量、淋巴结位置、弹性成像评分。(d)肿瘤相关特征结果：基因检测_BRAF、基因检测_TERT、原发灶_数量、原发灶_部位、原发灶_大小、淋巴结_中央区_转移数量、淋巴结_中央区_清扫数量、淋巴结_侧区_转移数量、淋巴结_侧区_清扫数量、淋巴结_转移数量、淋巴结_清扫数量、淋巴结_转移比例、淋巴结_大小、淋巴结_结外侵犯、组织学类型、AJCC8版分期_T分期、AJCC8版分期_N分期、AJCC8版分期_M分期、AJCC8版分期。
(3)将患者按照8：2随机拆分为训练集(2776)和测试集(695)。将连续型变量进行标准化。训练集用于建立预后模型，验证集用于对预后模型进行验证和评价。
(4)训练集中采用最小绝对收缩和选择算子方法(the least absoluteshrinkage and selection operator，LASSO)进行特征筛选，同时应用K折交叉验证法(Kfold cross validation)进行调整参数λ的优化和选择，并选择交叉验证误差最小时的λ作为模型最优值，将具有显著性的协变系数降为零，选择其余的非零系数。预实验中LASSO筛选(λ＝0.01)得到5个特征包括：原发灶大小、淋巴结转移比例、转移淋巴结最大直径、淋巴结结外侵犯、组织学类型。LASSO筛选过程如图1所示。
(5)将训练集中筛选所得5个特征输入多因素Cox回归模型，5个特征均是无疾病持续/复发生存期(RFS)的独立预测因子(P<0.05)；利用上述5个特征绘制用于预测无疾病持续/复发生存时间的列线图。
无疾病持续/复发生存期(RFS)指从手术日期到发现疾病持续/复发的日期为止，或到末次随访日期为止。
如图2所示，该列线图包含第1～12条直线，第1～12条直线由上至下依次排列，且相互平行；每条直线代表一个标尺，标尺上有刻度。第1个标尺表示第2～6个标尺上刻度对应的分值的标尺；分值范围为0～100，1在最左端，100在最右端；
第2个标尺表示肿瘤大小；若患者肿瘤大小≤10mm，在第1个标尺上对应的分值为0；若患者肿瘤大小>10且≤20mm，在第1个标尺上对应的分值为21；若患者肿瘤大小>20且≤40mm，在第1个标尺上对应的分值为42；若患者肿瘤大小>40mm，在第1个标尺上对应的分值为63；
第3个标尺表示转移淋巴结比例；若患者转移淋巴结比例≤30％，在第1个标尺上对应的分值为0；若患者转移淋巴结比例>30且≤50％，在第1个标尺上对应的分值为17；若患者转移淋巴结比例>50％且≤80％，在第1个标尺上对应的分值为32；若患者转移淋巴结比例>80％，在第1个标尺上对应的分值为48；
第4个标尺表示转移淋巴结最大直径；若患者无淋巴结转移，在第1个标尺上对应的分值为0；若患者转移淋巴结最大直径<30mm，在第1个标尺上对应的分值为9；若患者转移淋巴结最大直径≥30mm，在第1个标尺上对应的分值为18；
第5个标尺表示淋巴结结外侵犯情况；若患者无淋巴结结外侵犯，在第1个标尺上对应的分值为0；若患者有淋巴结结外侵犯，在第1个标尺上对应的分值为42；
第6个标尺表示组织学类型；若患者为乳头状癌(PTC)，在第1个标尺上对应的分值为0；若患者为滤泡亚型乳头状癌(FVPTC)，在第1个标尺上对应的分值为32；若患者为滤泡癌(FTC)，在第1个标尺上对应的分值为66；若患者为嗜酸细胞癌(Hurthle cell cancer)，在第1个标尺上对应的分值为100；
第7个标尺表示第2～6个标尺上刻度对应的分值的总和的标尺；刻度范围为0～180，0在最左端，180在最右端；
第8个标尺表示患者接受甲状腺全切术后6个月内未出现疾病持续或复发的概率(即第7个标尺上分值的总和在第8个标尺上对应的刻度值)，范围为0.2～0.95；
第9个标尺表示患者接受甲状腺全切术后1年内未出现疾病持续或复发的概率(即第7个标尺上分值的总和在第9个标尺上对应的刻度值)，范围为0.05～0.95；
第10个标尺表示患者接受甲状腺全切术后2年内未出现疾病持续或复发的概率(即第7个标尺上分值的总和在第10个标尺上对应的刻度值)，范围为0.05～0.95；
第11个标尺表示患者接受甲状腺全切术后3年内未出现疾病持续或复发的概率(即第7个标尺上分值的总和在第11个标尺上对应的刻度值)，范围为0.05～0.95；
第12个标尺表示患者接受甲状腺全切术后5年内未出现疾病持续或复发的概率(即第7个标尺上分值的总和在第12个标尺上对应的刻度值)，范围为0.05～0.95。
对于行甲状腺全切术后的分化型甲状腺癌患者，可以通过上述列线图第2～6个标尺中各变量的实际情况进行打分，各变量具体的分值刻度即为与第1列对应位置的分值。第2～6个标尺中5个变量的分值求和后可确定该患者的总分，并在第7个标尺找到对应的分值位置，继而可以根据总分在第8～12个标尺的对应位置得出该患者在接受甲状腺全切术后第6个月、第1年、第2年、第3年、第5年内未出现疾病持续或复发的预测概率。
(6)将测试集的以下5个特征输入步骤(5)构建的列线图中：原发灶大小、淋巴结转移比例、转移淋巴结最大径、淋巴结结外侵犯、组织学类型，利用列线图预测出测试集的患者在接受甲状腺全切术后第6个月、第1年、第2年、第3年、第5年内未出现疾病持续或复发的概率。
以下通过实验例证明本发明预测系统的效果。
实验例1本发明预测系统区分度和校准度验证
1、实验方法
(1)验证模型区分度：在测试集使用一致性指数(Harrell concordance index，C-Index)评估实施例1构建的模型的区分度。C-index的值介于0.5-1.0之间，0.5表示没有区分能力，1.0表示有较强区分能力，采用自助抽样分析200次来获得校正的C-Index。
(2)验证模型校准度：在测试集采用校准曲线评价实施例1构建的模型的校准度，校准曲线中对角的虚线即拟合曲线，代表实际复发情况，实线即标准曲线，表示模型预测复发的效能，越接近虚线表示预测效能越好。
2、实验结果
(1)模型区分度
结果显示，本发明的预测系统的C-index为0.74，显示出良好的区分度。
(2)模型校准度
模型预测的校准曲线如图3所示，校准曲线接近45°，说明模型预测的概率和实际概率有较高的一致性。
实验例2本发明预测系统与美国甲状腺协会的复发风险分层系统的预测效能进行比较
将训练集人群按照美国甲状腺协会的复发风险分层系统分为低危、中危、高危，并在训练集中利用Cox单因素回归模型分析该分层系统的预测能力，Cox单因素回归模型结果显示该复发风险分层系统能够在训练集中预测RFS(P＝0.001)。进一步在测试集中比较复发风险分层系统和本发明实施例1的列线图模型的预测效能，结果显示在测试集中复发风险分层系统的C-index为0.67，低于本发明中的列线图模型(C-index为0.74)。
由上述结果可知，与美国甲状腺协会的复发风险分层系统相比，本发明构建的预测系统具有更好的预测效果。
Claims (10)
1.一种预测甲状腺癌患者预后的预测系统，其特征在于：所述预测系统是以患者的肿瘤大小、淋巴结转移比例、转移淋巴结最大直径、淋巴结结外侵犯情况、组织学类型为预测指标构建的。
2.根据权利要求1所述的预测系统，其特征在于：所述预测系统是列线图，该列线图包含第1～12条直线，第1～12条直线由上至下依次排列，且相互平行；每条直线代表一个标尺，标尺上有刻度；
第1个标尺表示第2～6个标尺上刻度对应的分值的标尺；
第2个标尺表示肿瘤大小；
第3个标尺表示转移淋巴结比例；
第4个标尺表示转移淋巴结最大直径；
第5个标尺表示淋巴结结外侵犯情况；
第6个标尺表示组织学类型；
第7个标尺表示第2～6个标尺上刻度对应的分值的总和的标尺；
第8个标尺表示患者6个月内未出现疾病持续或复发的概率；
第9个标尺表示患者1年内未出现疾病持续或复发的概率；
第10个标尺表示患者2年内未出现疾病持续或复发的概率；
第11个标尺表示患者3年内未出现疾病持续或复发的概率；
第12个标尺表示患者5年内未出现疾病持续或复发的概率。
3.根据权利要求2所述的预测系统，其特征在于：第1个标尺的刻度范围为1～100，1在最左端，100在最右端；
第2个标尺中，若患者肿瘤大小≤10mm，在第1个标尺上对应的分值为0；若患者肿瘤大小>10且≤20mm，在第1个标尺上对应的分值为21；若患者肿瘤大小>20且≤40mm，在第1个标尺上对应的分值为42；若患者肿瘤大小>40mm，在第1个标尺上对应的分值为63；
第3个标尺中，若患者转移淋巴结比例≤30％，在第1个标尺上对应的分值为0；若患者转移淋巴结比例>30且≤50％，在第1个标尺上对应的分值为17；若患者转移淋巴结比例>50％且≤80％，在第1个标尺上对应的分值为32；若患者转移淋巴结比例>80％，在第1个标尺上对应的分值为48；
第4个标尺中，若患者无淋巴结转移，在第1个标尺上对应的分值为0；若患者转移淋巴结最大直径<30mm，在第1个标尺上对应的分值为9；若患者转移淋巴结最大直径≥30mm，在第1个标尺上对应的分值为18；
第5个标尺中，若患者无淋巴结结外侵犯，在第1个标尺上对应的分值为0；若患者有淋巴结结外侵犯，在第1个标尺上对应的分值为42；
第6个标尺中，若患者为乳头状癌，在第1个标尺上对应的分值为0；若患者为滤泡亚型乳头状癌，在第1个标尺上对应的分值为32；若患者为滤泡癌，在第1个标尺上对应的分值为66；若患者为嗜酸细胞癌，在第1个标尺上对应的分值为100；
第7个标尺的刻度范围为0～180，0在最左端，180在最右端；
第8个标尺中，患者6个月内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第8个标尺上对应的刻度值，范围为0.2～0.95；
第9个标尺中，患者1年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第9个标尺上对应的刻度值，范围为0.05～0.95；
第10个标尺中，患者2年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第10个标尺上对应的刻度值，范围为0.05～0.95；
第11个标尺中，患者3年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第11个标尺上对应的刻度值，范围为0.05～0.95；
第12个标尺中，患者5年内未出现疾病持续或复发的概率＝第7个标尺上分值的总和在第12个标尺上对应的刻度值，范围为0.05～0.95。
4.根据权利要求1～3任一项所述的预测系统，其特征在于：所述甲状腺癌患者为分化型甲状腺癌患者。
5.根据权利要求4所述的预测系统，其特征在于：所述分化型甲状腺癌患者为行甲状腺全切术后的患者。
6.根据权利要求2所述的预测系统，其特征在于：所述列线图的构建方法包括以下步骤：
(1)收集甲状腺癌患者的预测指标，录入输入模块；
(2)利用输入模块中的预测指标构建多因素Cox回归模型，绘制列线图。
7.根据权利要求2所述的预测系统，其特征在于：所述列线图如图2所示。
8.一种预测甲状腺癌患者预后的设备，其特征在于：所述设备包括权利要求1～7任一项所述的预测系统。
9.权利要求1～7任一项所述的预测系统在制备预测甲状腺癌患者预后的设备中的用途。
10.根据权利要求9所述的用途，其特征在于：所述甲状腺癌患者为分化型甲状腺癌患者；优选地，所述分化型甲状腺癌患者为行甲状腺全切术后的患者。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111182720.0A CN114242240A (zh) | 2021-10-11 | 2021-10-11 | 一种分化型甲状腺癌患者疾病持续复发预测模型 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111182720.0A CN114242240A (zh) | 2021-10-11 | 2021-10-11 | 一种分化型甲状腺癌患者疾病持续复发预测模型 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114242240A true CN114242240A (zh) | 2022-03-25 |
Family
ID=80743068
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111182720.0A Pending CN114242240A (zh) | 2021-10-11 | 2021-10-11 | 一种分化型甲状腺癌患者疾病持续复发预测模型 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114242240A (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106319062A (zh) * | 2016-08-23 | 2017-01-11 | 中南大学 | 一种用于甲状腺癌辅助诊断或疗效预测的微创试剂盒 |
| CN111354462A (zh) * | 2020-04-14 | 2020-06-30 | 中山大学孙逸仙纪念医院 | 晚期乳腺癌生存概率预测列线图、生存概率预测方法及患者分类方法 |
| CN112185569A (zh) * | 2020-09-11 | 2021-01-05 | 中山大学孙逸仙纪念医院 | 一种乳腺癌患者无病生存期预测模型及其构建方法 |
| CN112466464A (zh) * | 2020-12-17 | 2021-03-09 | 四川大学华西医院 | 初诊转移性前列腺癌患者的预后预测模型及其建立方法和应用 |
| CN112687394A (zh) * | 2021-01-05 | 2021-04-20 | 四川大学华西医院 | 转移性去势抵抗性前列腺癌患者在阿比特龙治疗中的预后预测模型及其建立方法和应用 |
| CN113284619A (zh) * | 2021-04-20 | 2021-08-20 | 四川大学华西医院 | 预测前列腺癌患者术后出现生化复发的方法及设备 |
-
2021
- 2021-10-11 CN CN202111182720.0A patent/CN114242240A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106319062A (zh) * | 2016-08-23 | 2017-01-11 | 中南大学 | 一种用于甲状腺癌辅助诊断或疗效预测的微创试剂盒 |
| CN111354462A (zh) * | 2020-04-14 | 2020-06-30 | 中山大学孙逸仙纪念医院 | 晚期乳腺癌生存概率预测列线图、生存概率预测方法及患者分类方法 |
| CN112185569A (zh) * | 2020-09-11 | 2021-01-05 | 中山大学孙逸仙纪念医院 | 一种乳腺癌患者无病生存期预测模型及其构建方法 |
| CN112466464A (zh) * | 2020-12-17 | 2021-03-09 | 四川大学华西医院 | 初诊转移性前列腺癌患者的预后预测模型及其建立方法和应用 |
| CN112687394A (zh) * | 2021-01-05 | 2021-04-20 | 四川大学华西医院 | 转移性去势抵抗性前列腺癌患者在阿比特龙治疗中的预后预测模型及其建立方法和应用 |
| CN113284619A (zh) * | 2021-04-20 | 2021-08-20 | 四川大学华西医院 | 预测前列腺癌患者术后出现生化复发的方法及设备 |
Non-Patent Citations (2)
| Title |
|---|
| 李春明: "89例甲状腺癌临床病理分析", 《医学信息(中旬刊)》 * |
| 隋鑫等: "术前甲状腺癌超声检查特征与颈部淋巴结转移的多因素分析", 《河北医药》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8510245B2 (en) | Bayesian clinical decision model for determining probability of transplant glomerulopathy | |
| Berlin et al. | Genomic classifier for guiding treatment of intermediate-risk prostate cancers to dose-escalated image guided radiation therapy without hormone therapy | |
| Bobdey et al. | A nomogram based prognostic score that is superior to conventional TNM staging in predicting outcome of surgically treated T4 buccal mucosa cancer: time to think beyond TNM | |
| Bhardwaj et al. | Machine learning for endometrial cancer prediction and prognostication | |
| JP2017520775A (ja) | 結腸直腸腫瘍を検出するためのタンパク質バイオマーカープロファイル | |
| CN109055562A (zh) | 一种生物标志物、预测肾细胞癌的复发和死亡风险的方法 | |
| WO2018071896A1 (en) | Methods and systems for assessing fertility based on subclinical genetic factors | |
| Zheng et al. | A three‐gene panel that distinguishes benign from malignant thyroid nodules | |
| JP2017523437A (ja) | 体軸性脊椎関節炎の疾患活動性を測定およびモニタリングするためのバイオマーカーおよび方法 | |
| CN108531597A (zh) | 一种用于口腔鳞癌早期诊断的检测试剂盒 | |
| JP2024512454A (ja) | 手術リスクスコアを生成するシステムおよび方法ならびにそれらの使用 | |
| AU730373B2 (en) | Prognostic protocol | |
| Hei et al. | A nomogram predicting contralateral central neck lymph node metastasis for papillary thyroid carcinoma | |
| US20240127955A1 (en) | System and methods for mapping cancer classifications | |
| US20140236621A1 (en) | Method for determining a predictive function for discriminating patients according to their disease activity status | |
| CN114242240A (zh) | 一种分化型甲状腺癌患者疾病持续复发预测模型 | |
| Koskas et al. | Independent external validation of radiotherapy and its impact on the accuracy of a nomogram for predicting survival of women with endometrial cancer | |
| EP4343775A1 (en) | A method for determining a medical outcome for an individual, related electronic system and computer program | |
| Fang et al. | Deep learning predicts biomarker status and discovers related histomorphology characteristics for low-grade glioma | |
| Sharma et al. | Quantitative assessment of tumor-associated tissue eosinophilia and nuclear organizing region activity to validate the significance of the pattern of invasion in oral squamous cell carcinoma: A retrospective study | |
| Farolfi et al. | frontiers Frontiers in Oncology 2022 published: 27 July REVIEW | |
| Cui et al. | Novel alternative tools for metastatic pheochromocytomas/paragangliomas prediction | |
| Reguly et al. | 3.4 kb mitochondrial genome deletion serves as a surrogate predictive biomarker for prostate cancer in histopathologically benign biopsy cores | |
| Ou et al. | Competing Risk Models for Predicting Prognosis of Cervical Cancer Patients based on Surveillance, Epidemiology and End Results (SEER) Database | |
| Wu et al. | Construction of prediction models for determining the risk of lateral lymph node metastasis in patients with thyroid papillary carcinoma based on gender stratification |
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 |