一种甲状腺癌早期辅助诊断标志物的制作方法

1.本发明涉及医学领域，特别涉及一种甲状腺癌早期辅助诊断标志物。


背景技术：

2.甲状腺癌(thyroid cancer)是最常见的甲状腺恶性肿瘤，约占全身恶性肿瘤的1％。包括甲状腺乳头状癌、滤泡性甲状腺癌、未分化甲状腺癌和髓样癌。除髓样癌外，绝大部分甲状腺癌起源于滤泡上皮细胞。其中乳头状癌(papillary thyroid cancer,ptc)最为常见，占所有甲状腺恶性肿瘤的90％以上[xing,mingzhao；haugen,bryan r；schlumberger,martin(2013).progress in molecular-based management of differentiated thyroid cancer.the lancet,381(9871),1058
–
1069.]。据统计，成人甲状腺结节的患病率约为5-10％，其中60岁以上的人群最为严重，高达50-70％[guth s,theune u,aberle j,et al.very high prevalence of thyroid nodules detected by high frequency(13mhz)ultrasound examination.eur j clin invest 2009；39:699-706.]。影像学检查是比较普遍的甲状腺诊断方法，该方法大多依靠医生经验判断，存在一定的结果误差，且影像学对人身体有一定的辐射伤害。细针穿刺活检也是临床上常见的甲状腺癌诊断技术，该方法根据穿刺物的细胞学形态对结节的良恶性进行评估。由于甲状腺良恶性肿瘤的细胞学特征经常发生重叠，因此约有10-30％的细针穿刺诊断为不明确的细胞学结果[cibas es,ali sz.the 2017bethesda system for reporting thyroid cytopathology.thyroid.2017；27(11):1341-6.]。不确定的穿刺结果导致约60％的患者遭受过度治疗或者漏诊[stewart r,leang yj,bhatt cr,grodski s,serpell j,lee jc.quantifying the differences in surgical management of patients with definitive and indeterminate thyroid nodule cytology.eur j surg oncol.2020；46(2):252-7.]。这不仅增加了患者经济和身心负担，还占据了大量的公共卫生资源，导致医疗保健系统的巨额财务成本。因此，鉴别甲状腺良恶性肿瘤有利于临床医师采取更精准的治疗方案，具有重要的临床和公共卫生意义。
[0003]
表观遗传学是一种不涉及dna序列改变但可遗传的基因表达调控方式，并能够遗传给下一代[nicoglou a,merlin f.epigenetics:away to bridge the gap between biological fields.stud hist philos biol biomed sci.2017；66:73-82]。dna甲基化是表观遗传调控的重要方式之一，是指在dna甲基化转移酶的作用下，在基因组cpg二核苷酸的胞嘧啶5’碳位共价键结合一个甲基基团[bird a.perceptions of epigenetics.nature.2007；447:396-398]。大量研究表明，dna甲基化能引起染色质结构、dna构象、dna稳定性及dna与蛋白质相互作用方式的改变，从而控制基因表达[moore ld,le t,fan g.dnamethylation and its basic function.neuropsychopharmacology.2013；38:23-38]。
[0004]
目前，临床上关于甲状腺癌诊断标志物的灵敏度和特异性很有限，尤其是缺乏早期诊断的标志物，因此更为敏感、特异的早期分子标记物亟待发掘。


技术实现要素：

[0005]
本发明的目的是提供一种甲状腺癌早期辅助诊断标志物。
[0006]
第一方面，本发明要求保护一种试剂盒。
[0007]
本发明要求保护的试剂盒包括用于检测bcar3基因甲基化水平的物质；所述试剂盒的用途为如下中的至少一种：
[0008]
(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0009]
(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0010]
(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0011]
(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
[0012]
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。
[0013]
进一步地，(2)和(4)中所述不同亚型可为病理分型，如组织学分型。
[0014]
进一步地，(3)和(5)中所述不同分期可为临床分期。
[0015]
在本发明的具体实施方式中，(2)中所述区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤具体可为如下任一种：区分或辅助区分甲状腺良性肿瘤和甲状腺乳头状癌、区分或辅助区分甲状腺良性肿瘤和甲状腺滤泡癌、区分或辅助区分甲状腺良性肿瘤和甲状腺髓样癌、区分或辅助区分甲状腺良性肿瘤和甲状腺未分化癌。
[0016]
在本发明的具体实施方式中，(3)中所述区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤具体可为如下任一种：区分或辅助区分甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤、区分或辅助区分甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤、区分或辅助区分甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤、区分或辅助区分甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0017]
在本发明的具体实施方式中，(4)中所述区分或辅助区分甲状腺恶性肿瘤不同亚型具体可为如下任一种：区分或辅助区分甲状腺乳头状癌和甲状腺滤泡癌、区分或辅助区分甲状腺乳头状癌和甲状腺髓样癌、区分或辅助区分甲状腺乳头状癌和甲状腺未分化癌、区分或辅助区分甲状腺滤泡癌和甲状腺髓样癌、区分或辅助区分甲状腺滤泡癌和甲状腺未分化癌、区分或辅助区分甲状腺髓样癌和甲状腺未分化癌。
[0018]
在本发明的具体实施方式中，(5)中所述区分或辅助区分甲状腺恶性肿瘤不同分期具体可为如下任一种：区分或辅助区分ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤、区分或辅助区分ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、区分或辅助区分ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、区分或辅助区分ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、区分或辅助区分ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、区分或辅助区分ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0019]
进一步地，所述bcar3基因甲基化水平可为bcar3基因中如下(a1)-(a3)所示片段中的全部或部分cpg位点的甲基化水平：
[0020]
(a1)seq id no.1所示的dna片段或与其具有80％以上同一性的dna片段；
[0021]
(a2)seq id no.2所示的dna片段或与其具有80％以上同一性的dna片段；
[0022]
(a3)seq id no.3所示的dna片段或与其具有80％以上同一性的dna片段。
[0023]
更进一步地，所述全部或部分cpg位点可为如下任一所示cpg位点：
[0024]
(b1)bcar3基因中seq id no.1、seq id no.2和seq id no.3所示3个dna片段中的
任意一个或多个cpg位点；
[0025]
(b2)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点(表1)和seq id no.2所示的dna片段上的所有cpg位点(表2)；
[0026]
(b3)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点(表1)和seq id no.3所示的dna片段上的所有cpg位点(表3)；
[0027]
(b4)bcar3基因中seq id no.2所示的dna片段上的所有cpg位点(表2)和seq id no.3所示的dna片段上的所有cpg位点(表3)；
[0028]
(b5)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点(表1)、seq id no.2所示的dna片段上的所有cpg位点(表2)和seq id no.3所示的dna片段上的所有cpg位点(表3)；
[0029]
(b6)bcar3基因中seq id no.2所示的dna片段中的全部cpg位点(表2)或任意10个或任意9个或任意8个或任意7个或任意6个或任意5个或任意4个或任意3个或任意2个或任意1个cpg位点；
[0030]
(b7)bcar3基因中seq id no.2所示的dna片段上如下5项所示cpg位点的全部或任意4项或任意3项或任意2项或任意1项：
[0031]
第1项：seq id no.2所示的dna片段自5’端第117-118位所示cpg位点；
[0032]
第2项：seq id no.2所示的dna片段自5’端第329-330位所示cpg位点；
[0033]
第3项：seq id no.2所示的dna片段自5’端第353-354位和第355-356位和第357-358位和第359-360位所示cpg位点；
[0034]
第4项：seq id no.2所示的dna片段自5’端第384-385位所示cpg位点；
[0035]
第5项：seq id no.2所示的dna片段自5’端第581-582位所示cpg位点。
[0036]
在本发明的具体实施方式中，有些相邻的甲基化位点在利用飞行时间质谱进行dna甲基化分析时由于几个cpg位点位于一个甲基化片段上，峰图无法区分(无法区分的位点在表5中有记载)，因而在进行甲基化水平分析、以及构建和使用相关数学模型时将其按照一个甲基化位点进行处理。
[0037]
在所述试剂盒中，所述用于检测bcar3基因甲基化水平的物质可包含(或为)用于扩增bcar3基因全长或部分片段的引物组合。
[0038]
进一步地，所述部分片段可为如下中至少一个片段：
[0039]
(g1)seq id no.1所示的dna片段或其包含的dna片段；
[0040]
(g2)seq id no.2所示的dna片段或其包含的dna片段；
[0041]
(g3)seq id no.3所示的dna片段或其包含的dna片段；
[0042]
(g4)与seq id no.1所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段；
[0043]
(g5)与seq id no.2所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段；
[0044]
(g6)与seq id no.3所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段。
[0045]
在本发明的具体实施方式中，所述引物组合为引物对a和/或引物对b和/或引物对c。
[0046]
所述引物对a为引物a1和引物a2组成的引物对；所述引物a1为seq id no.4或seq id no.4的第11-35位核苷酸所示的单链dna；所述引物a2为seq id no.5或seq id no.5的第32-56位核苷酸所示的单链dna；
[0047]
所述引物对b为引物b1和引物b2组成的引物对；所述引物b1为seq id no.6或seq id no.6的第11-35位核苷酸所示的单链dna；所述引物b2为seq id no.7或seq id no.7的第32-56位核苷酸所示的单链dna；
[0048]
所述引物对c为引物c1和引物c2组成的引物对；所述引物c1为seq id no.8或seq id no.8的第11-35位核苷酸所示的单链dna；所述引物c2为seq id no.9或seq id no.9的第32-56位核苷酸所示的单链dna。
[0049]
根据需要，所述试剂盒中还可含有储存有数学模型和/或数学模型使用方法的介质。
[0050]
所述数学模型按照包括如下步骤的方法获得：
[0051]
(c1)分别检测n1个a类型样本和n2个b类型样本的基因甲基化水平；
[0052]
(c2)取步骤(c1)获得的所有样本的基因甲基化水平数据，按照a类型和b类型的分类方式，通过二分类逻辑回归法建立数学模型，确定分类判定的阈值。
[0053]
其中，n1和n2均可为10以上正整数。
[0054]
所述数学模型使用方法包括如下步骤：
[0055]
(d1)检测待测样本的基因甲基化水平；
[0056]
(d2)将步骤(d1)获得的所述待测样本的基因甲基化水平数据代入所述数学模型，得到检测指数；然后比较检测指数和阈值的大小，根据比较结果确定所述待测样本的类型是a类型还是b类型。
[0057]
在本发明的具体实施方式中，所述阈值设为0.5。大于0.5归为一类，小于0.5归为另外一类，等于0.5作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组是a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0058]
在实际应用中，所述阈值也可根据最大约登指数确定(具体可为最大约登指数对应的数值)。大于阈值归为一类，小于阈值归为另外一类，等于阈值作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0059]
所述a类型样本和所述b类型样本为如下任一种：
[0060]
(e1)甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0061]
(e2)甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0062]
(e3)甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0063]
(e4)不同亚型的甲状腺恶性肿瘤；
[0064]
(e5)不同分期的甲状腺恶性肿瘤。
[0065]
进一步地，(e2)和(e4)中所述不同亚型可为病理分型，如组织学分型。
[0066]
进一步地，(e3)和(e5)中所述不同分期可为临床分期。
[0067]
在本发明的具体实施方式中，(e2)中所述甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤具体可为如下任一种：甲状腺良性肿瘤和甲状腺乳头状癌、甲状腺良性肿瘤和甲状腺滤泡癌、甲状腺良性肿瘤和甲状腺髓样癌、甲状腺良性肿瘤和甲状腺未分化癌。
[0068]
在本发明的具体实施方式中，(e3)中所述甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤具体可为如下任一种：甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0069]
在本发明的具体实施方式中，(e4)中所述甲状腺恶性肿瘤不同亚型具体可为如下任一种：甲状腺乳头状癌和甲状腺滤泡癌、甲状腺乳头状癌和甲状腺髓样癌、甲状腺乳头状癌和甲状腺未分化癌、甲状腺滤泡癌和甲状腺髓样癌、甲状腺滤泡癌和甲状腺未分化癌、甲状腺髓样癌和甲状腺未分化癌。
[0070]
在本发明的具体实施方式中，(e5)中所述甲状腺恶性肿瘤不同分期具体可为如下任一种：ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0071]
第二方面，本发明要求保护一种系统。
[0072]
本发明要求保护的系统，包括：
[0073]
(f1)用于检测bcar3基因甲基化水平的试剂和/或仪器；
[0074]
在(f1)中，所述用于检测bcar3基因甲基化水平的试剂可为前文第一方面中所述的用于检测bcar3基因甲基化水平的物质(如引物对)。用于检测bcar3基因甲基化水平的仪器可为飞行时间质谱检测仪。当然所述用于检测bcar3基因甲基化水平的试剂中还可包含进行飞行时间质谱所用的其他常规试剂。
[0075]
(f2)装置，所述装置包括单元x和单元y；
[0076]
所述单元x用于建立数学模型，包括数据采集模块、数据分析处理模块和模型输出模块；
[0077]
所述数据采集模块被配置为采集(f1)检测得到的n1个a类型样本和n2个b类型样本的bcar3基因甲基化水平数据；
[0078]
所述数据分析处理模块被配置为接收来自于所述数据采集模块的所述n1个a类型样本和n2个b类型样本的bcar3基因甲基化水平数据，按照a类型和b类型的分类方式，通过二分类逻辑回归法建立数学模型，确定分类判定的阈值；
[0079]
其中，n1和n2均可为10以上正整数。
[0080]
所述模型输出模块被配置为接收来自于所述数据分析处理模块建立的所述数学模型，并进行输出；
[0081]
所述单元y用于确定待测样本类型，包括数据输入模块、数据运算模块、数据比较模块和结论输出模块；
[0082]
所述数据输入模块被配置为输入(f1)检测得到的待测者的bcar3基因甲基化水平数据；
[0083]
所述数据运算模块被配置为接收来自于所述数据输入模块的所述待测者的bcar3基因甲基化水平数据，并将所述待测者的bcar3基因甲基化水平数据代入所述单元x中的所述数据分析处理模块建立的所述数学模型，计算得到检测指数；
[0084]
所述数据比较模块被配置为接收来自于所述数据运算模块计算得到的检测指数，
并将所述检测指数与所述单元x中的所述数据分析处理模块中确定的所述阈值进行比较；
[0085]
所述结论输出模块被配置为接收来自于所述数据比较模块的比较结果，并根据所述比较结果输出所述待测样本的类型是a类型还是b类型的结论；
[0086]
所述a类型样本和所述b类型样本为如下任一种：
[0087]
(e1)甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0088]
(e2)甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0089]
(e3)甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0090]
(e4)不同亚型的甲状腺恶性肿瘤；
[0091]
(e5)不同分期的甲状腺恶性肿瘤。
[0092]
在本发明的具体实施方式中，所述阈值设为0.5。大于0.5归为一类，小于0.5归为另外一类，等于0.5作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组是a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0093]
在实际应用中，所述阈值也可根据最大约登指数确定(具体可为最大约登指数对应的数值)。大于阈值归为一类，小于阈值归为另外一类，等于阈值作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0094]
进一步地，(e2)和(e4)中所述不同亚型可为病理分型，如组织学分型。
[0095]
进一步地，(e3)和(e5)中所述不同分期可为临床分期。
[0096]
在本发明的具体实施方式中，(e2)中所述甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤具体可为如下任一种：甲状腺良性肿瘤和甲状腺乳头状癌、甲状腺良性肿瘤和甲状腺滤泡癌、甲状腺良性肿瘤和甲状腺髓样癌、甲状腺良性肿瘤和甲状腺未分化癌。
[0097]
在本发明的具体实施方式中，(e3)中所述甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤具体可为如下任一种：甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0098]
在本发明的具体实施方式中，(e4)中所述甲状腺恶性肿瘤不同亚型具体可为如下任一种：甲状腺乳头状癌和甲状腺滤泡癌、甲状腺乳头状癌和甲状腺髓样癌、甲状腺乳头状癌和甲状腺未分化癌、甲状腺滤泡癌和甲状腺髓样癌、甲状腺滤泡癌和甲状腺未分化癌、甲状腺髓样癌和甲状腺未分化癌。
[0099]
在本发明的具体实施方式中，(e5)中所述甲状腺恶性肿瘤不同分期具体可为如下任一种：ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤。
[0100]
第三方面，本发明要求保护前文第一方面中所述用于检测bcar3基因甲基化水平的物质在制备产品中的应用；所述产品的用途为如下至少一种：
[0101]
(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0102]
(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0103]
(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0104]
(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
[0105]
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。
[0106]
第四方面，本发明要求保护前文第一方面中所述用于检测bcar3基因甲基化水平的物质和所述储存有数学模型和/或数学模型使用方法的介质在制备产品中的应用；所述产品的用途为如下至少一种：
[0107]
(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0108]
(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0109]
(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0110]
(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
[0111]
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。
[0112]
第五方面，本发明要求保护前文第一方面中所述储存有数学模型和/或数学模型使用方法的介质在制备产品中的应用；所述产品的用途为如下至少一种：
[0113]
(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0114]
(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0115]
(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0116]
(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
[0117]
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。
[0118]
第六方面，本发明要求保护甲基化bcar3基因作为标志物在制备产品中的应用；所述产品的用途为如下至少一种：
[0119]
(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；
[0120]
(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0121]
(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0122]
(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
[0123]
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。
[0124]
在第三至第六方面中，(1)-(5)的具体含义同第一方面中所述。
[0125]
另外，本发明还要求保护一种区分待测样本为a类型样本还是b类型样本的方法。该方法可包括如下步骤：
[0126]
(a)可按照包括如下步骤的方法建立数学模型：
[0127]
(a1)分别检测n1个a类型样本和n2个b类型样本的bcar3基因甲基化水平(训练集)；
[0128]
(a2)取步骤(a1)获得的所有样本的bcar3基因甲基化水平数据，按照a类型和b类型的分类方式，通过二分类逻辑回归法建立数学模型，确定分类判定的阈值。
[0129]
其中，(a1)中的n1和n2均为10以上的正整数。
[0130]
(b)可按照包括如下步骤的方法确定所述待测样本为a类型样本还是b类型样本：
[0131]
(b1)检测所述待测样本的bcar3基因甲基化水平；
[0132]
(b2)将步骤(b1)获得的所述待测样本的bcar3基因甲基化水平数据代入所述数学模型，得到检测指数；然后比较检测指数和阈值的大小，根据比较结果确定所述待测样本的类型是a类型还是b类型。
[0133]
在本发明的具体实施方式中，所述阈值设为0.5。大于0.5归为一类，小于0.5归为
另外一类，等于0.5作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组是a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0134]
在实际应用中，所述阈值也可根据最大约登指数确定(具体可为最大约登指数对应的数值)。大于阈值归为一类，小于阈值归为另外一类，等于阈值作为不确定的灰区。其中a类型和b类型为相对应的两分类，二分类的分组，哪一组是a类型，哪一组是b类型，要根据具体的数学模型来确定，无需约定。
[0135]
所述a类型样本和所述b类型样本可为前文(e1)-(e5)中的任一种。
[0136]
以上任一所述数学模型在实际应用中可能会根据dna甲基化的检测方法以及拟合方式不同有所改变，要根据具体的数学模型来确定，无需约定。
[0137]
在本发明的实施例中，所述模型具体为ln(y/(1-y))＝b0+b1x1+b2x2+b3x3+
…
+bnxn，其中y为因变量即将待测样品的一个或者多个甲基化位点的甲基化值代入模型以后得出的检测指数，b0为常量，x1-xn为自变量即为该测试样品的一个或者多个甲基化位点的甲基化值(每一个值为0-1之间的数值)，b1-bn为模型赋予每一个位点甲基化值的权重。
[0138]
本发明的实施例中建立的一个具体模型为用于区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤的模型，所述模型具体为：ln(y/(1-y))＝3.793-1.197*bcar3_b_2
ꢀ‑
6.621*bcar3_b_3+0.790*bcar3_b_4.5.6.7+0.128*bcar3_b_8+2.105*bcar3_b_9。所述bcar3_b_2为seq id no.2所示的dna片段自5’端第117-118位所示cpg位点的甲基化水平；所述bcar3_b_3为seq id no.2所示的dna片段自5’端第329-330位所示cpg位点的甲基化水平；所述bcar3_b_4.5.6.7为seq id no.2所示的dna片段自5’端第353-354位和第355-356位和第357-358位和第359-360位所示cpg位点的甲基化水平；所述bcar3_b_8为seq id no.2所示的dna片段自5’端第384-385位所示cpg位点的甲基化水平；所述bcar3_b_9为seq id no.2所示的dna片段自5’端第581-582位所示cpg位点的甲基化水平。所述模型的阈值为0.5。通过模型计算的检测指数大于0.5的患者候选为甲状腺恶性肿瘤患者，小于0.5的患者候选为甲状腺良性肿瘤患者。
[0139]
在上述各方面中，所述检测bcar3基因甲基化水平为检测肿瘤组织样本中bcar3基因甲基化水平。
[0140]
在本发明中，甲状腺恶性肿瘤组织中的bcar3基因中seq id no.1、seq id no.2和seq id no.3所示的dna片段上甲基化位点的甲基化水平明显低于甲状腺良性肿瘤。
[0141]
在本发明中，不同临床特征甲状腺癌如：乳头状癌、滤泡癌、髓样癌和未分化癌肿瘤组织中的bcar3基因中seq id no.1、seq id no.2和seq id no.3所示的dna片段上甲基化位点的甲基化水平越来越低。
[0142]
在本发明中，随着甲状腺恶性肿瘤的分期的增加组织中的bcar3基因中seq id no.1、seq id no.2和seq id no.3所示的dna片段上甲基化位点的甲基化水平越来越低。
[0143]
以上任一所述bcar3基因具体可参见genbank登录号：nm_001261408.2(gi:1675178356)，转录变体1；nm_003567.4(gi:1519311780)，转录变体2；nm_001261409.1(gi:387157908)，转录变体3；nm_001261410.2(gi:1675175859)，转录变体4；nm_001308251.1(gi:815891261)，转录变体5；nm_001412043.1(gi:2302515426)，转录变体7；nm_001412044.1(nm_001412044.1gi:2303311949)，转录变体8；nm_001412045.1(gi:2303318054)，转录变体9；nm_001412046.1(gi:2303321770)，转录变体10；nm_001412048.1
(gi:2303321774)，转录变体12；nm_001412049.1(gi:2303321776)，转录变体13；nm_001412050.1(gi:2303321778)，转录变体14；nm_001412051.1(gi:2303321780)，转录变体15；nm_001412052.1(gi:2303321783)，转录变体16；nm_001412054.1(gi:2302515512)，转录变体17；nm_001412055.1(gi:2302515475)转录变体18；nm_001412056.1(gi:2302515533)，转录变体19；nm_001412057.1(gi:2302515479)转录变体20；nm_001412058.1(gi:2302515417)，转录变体21；nm_001412059.1(gi:2302515497)，转录变体22；nm_001412060.1(gi:2302515409)，转录变体23；nm_001412061.1(gi:2302515435)，转录变体24；nm_001412062.1(gi:2302515466)，转录变体25；nm_001412064.1(gi:2302515451)，转录变体27；nm_001412065.1(gi:2302515508)，转录变体28；nm_001412066.1(gi:2302515449)，转录变体29；nm_001412067.1(gi:2302515528)，转录变体30；nm_001412072.1(gi:2302515488)，转录变体31；nm_001412073.1(gi:2302515437)，转录变体32；nm_001412074.1(gi:2302515492)，转录变体33；nm_001412075.1(gi:2302515490)，转录变体34；nm_001412076.1(gi:2302515414)转录变体35；nm_001412077.1(gi:2302515526)转录变体36。
[0144]
本发明证明活检样本bcar3甲基化可作为甲状腺良性肿瘤和甲状腺恶性肿瘤、不同亚型或不同分期的甲状腺恶性肿瘤的鉴别诊断的潜在标志物。本发明对鉴别甲状腺良性肿瘤和甲状腺恶性肿瘤、不同亚型或不同分期的甲状腺恶性肿瘤，以及指导制定合理的临床治疗方案均有重要的科学意义和临床应用价值。
附图说明
[0145]
图1为数学模型示意图。
[0146]
图2为甲状腺良恶性肿瘤数学模型举例说明。
具体实施方式
[0147]
下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。以下提供的实施例可作为本技术领域普通技术人员进行进一步改进的指南，并不以任何方式构成对本发明的限制。
[0148]
下述实施例中的实验方法，如无特殊说明，均为常规方法，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
[0149]
实施例1、用于检测bcar3基因甲基化位点的引物设计
[0150]
本次检测选择了bcar3基因三个片段(bcar3_a片段、bcar3_b片段和bcar3_c片段)上的cpg位点进行甲基化水平和甲状腺恶性肿瘤的相关性分析。
[0151]
bcar3_a片段(seq id no.1)位于hg19参考基因组chr1：94241531-94242129，负义链。
[0152]
bcar3_b片段(seq id no.2)位于hg19参考基因组chr1：94240654-94241338，负义链。
[0153]
bcar3_c片段(seq id no.3)位于hg19参考基因组chr1：94239588-94240266，负义链。
[0154]
bcar3_a片段中的位点信息如表1所示。
[0155]
bcar3_b片段中的位点信息如表2所示。
[0156]
bcar3_c片段中的位点信息如表3所示。
[0157]
表1、bcar3_a片段中cpg位点信息
[0158]
cpg位点cpg位点在序列中的位置bcar3_a_1seq id no.1自5’端第26-27位bcar3_a_2seq id no.1自5’端第105-106位bcar3_a_3seq id no.1自5’端第182-183位bcar3_a_4seq id no.1自5’端第224-225位bcar3_a_5seq id no.1自5’端第321-322位bcar3_a_6seq id no.1自5’端第566-567位
[0159]
表2、bcar3_b片段中cpg位点信息
[0160]
cpg位点cpg位点在序列中的位置bcar3_b_1seq id no.2自5’端第26-27位bcar3_b_2seq id no.2自5’端第117-118位bcar3_b_3seq id no.2自5’端第329-330位bcar3_b_4seq id no.2自5’端第353-354位bcar3_b_5seq id no.2自5’端第355-356位bcar3_b_6seq id no.2自5’端第357-358位bcar3_b_7seq id no.2自5’端第359-360位bcar3_b_8seq id no.2自5’端第384-385位bcar3_b_9seq id no.2自5’端第581-582位bcar3_b_10seq id no.2自5’端第606-607位bcar3_b_11seq id no.2自5’端第659-660位
[0161]
表3、bcar3_c片段中cpg位点信息
[0162]
cpg位点cpg位点在序列中的位置bcar3_c_1seq id no.3自5’端第26-27位bcar3_c_2seq id no.3自5’端第58-59位bcar3_c_3seq id no.3自5’端第240-241位bcar3_c_4seq id no.3自5’端第272-273位bcar3_c_5seq id no.3自5’端第278-279位bcar3_c_6seq id no.3自5’端第312-313位bcar3_c_7seq id no.3自5’端第358-359位bcar3_c_8seq id no.3自5’端第366-367位bcar3_c_9seq id no.3自5’端第370-371位bcar3_c_10seq id no.3自5’端第429-430位bcar3_c_11seq id no.3自5’端第459-460位bcar3_c_12seq id no.3自5’端第566-567位
bcar3_c_13seq id no.3自5’端第653-654位
[0163]
针对四个片段(bcar3_a片段、bcar3_b片段和bcar3_c片段)设计特异性pcr引物，如表4所示。seq id no.4、seq id no.6和seq id no.8为正向引物；seq id no.5、seq id no.7和seq id no.9为反向引物。seq id no.4、seq id no.6和seq id no.8中自5’端第1至10位为非特异性标签，第11至35位为特异性引物序列；seq id no.5、seq id no.7和seq id no.9中自5’端第1至31位为非特异性标签，第32至56位为特异性引物序列。引物序列中不包含snp和cpg位点。
[0164]
表4、bcar3甲基化引物序列
[0165][0166][0167]
实施例2、bcar3基因甲基化检测及结果分析
[0168]
一、研究样本
[0169]
经患者知情同意，共收集380例甲状腺良性肿瘤组织和598例甲状腺恶性肿瘤组织。甲状腺癌分期以美国癌症联合会(ajcc)第八版分期系统为判断标准。根据病理类型，甲状腺恶性肿瘤包括甲状腺乳头状癌、甲状腺滤泡癌、髓样癌和未分化癌四大类。本次收集的598例甲状腺恶性肿瘤患者中包括甲状腺乳头状癌380例，甲状腺滤泡癌138例，甲状腺髓样癌44例，甲状腺未分化癌36例。按照病理分期划分，598例甲状腺恶性肿瘤患者中有470例ⅰ期患者，68例ⅱ期患者，24例ⅲ期患者，36例ⅳ期患者。
[0170]
二、甲基化检测
[0171]
1、提取肿瘤组织中的总dna。
[0172]
2、将步骤1制备的组织样本总dna进行重亚硫酸盐处理(参照qiagen的dna甲基化试剂盒说明书操作)。重亚硫酸盐处理后，原来cpg位点中未发生甲基化的胞嘧啶(c)被转化成尿嘧啶(u)，而发生甲基化的胞嘧啶保持不变。
[0173]
3、以步骤2经过重亚硫酸盐处理的dna为模板，采用表4中的3对特异引物对通过dna聚合酶按照常规pcr反应要求的反应体系进行pcr扩增，所有引物都采用常规的标准pcr反应体系，且都按照以下程序进行扩增。
[0174]
pcr反应程序为：95℃，4min
→
(95℃，20s
→
56℃，30s
→
72℃，2min)45个循环
→
72℃，5min
→
4℃，1h。
[0175]
4、取步骤3的扩增产物，通过飞行时间质谱进行dna甲基化分析，具体方法如下：
[0176]
(1)向5μl pcr产物中加入2μl虾碱性磷酸盐(sap)溶液(0.3ml sap[0.5u]+1.7ml h2o)然后按照以下程序在pcr仪中孵育(37℃,20min
→
85℃,5min
→
4℃,5min)；
[0177]
(2)取出2μl步骤(1)得到的sap处理后的产物，根据说明书加入5μl t-cleavage反应体系中，然后在37℃孵育3h；
[0178]
(3)取步骤(2)的产物，加入19μl去离子水，再用6μg resin在旋转摇床进行去离子化孵育1h；
[0179]
(4)2000rpm室温离心5min，将微量上清由nanodispenser机械手臂上样384spectrochip；
[0180]
(5)飞行时间质谱分析；获得的数据用spectroacquire v3.3.1.3软件收集，通过massarray epityper v1.2软件实现可视化。
[0181]
上述飞行时间质谱检测使用的试剂均来试剂盒(t-cleavage masscleave reagent auto kit，货号：10129a)；上述飞行时间质谱检测使用的检测仪器为massarray
○
r analyzer chip prep module 384，型号：41243；上述数据分析软件为检测仪器自带软件。
[0182]
5、对步骤4得到的数据进行分析。
[0183]
数据统计分析由spss statistics 23.0进行。
[0184]
非参数检验用于两组之间的比较分析。
[0185]
多个cpg位点的组合对于不同样品分组的鉴别效果通过逻辑回归和受试者曲线的统计学方法得以实现。
[0186]
所有的统计检验都是双侧的，p值《0.05被认为具有统计学意义。
[0187]
通过质谱实验，共获得26个可以区别的峰图。采用spectroacquire v3.3.1.3软件根据“甲基化水平＝甲基化片段的峰面积/(非甲基化片段的峰面积+甲基化片段的峰面积)”公式可自动通过计算峰面积得到每个样本在每个cpg位点的甲基化水平。
[0188]
三、统计分析
[0189]
采用中位数表示甲基化水平，采用非参数检验比较两组或多组之间的甲基化水平差异。通过受试者工作特征曲线(receiver operating characteristic curve，roc曲线)评价单个cpg位点和多个cpg位点组合的鉴别诊断价值。以双侧p＜0.05为差异有统计学意义，所有数据均通过spss 25.0进行统计分析。
[0190]
四、结果分析
[0191]
1、甲状腺良性肿瘤和甲状腺恶性肿瘤bcar3基因甲基化水平分析
[0192]
以380例甲状腺良性肿瘤和598例甲状腺恶性肿瘤的组织样本为研究材料分析bcar3基因中所有cpg位点的甲基化水平。甲状腺良性肿瘤、甲状腺恶性肿瘤、甲状腺乳头状癌、甲状腺滤泡癌、髓样癌、未分化癌、ⅰ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤病例bcar3基因3个片段的甲基化水平的具体结果见表5。结果表明，甲状腺良性肿瘤bcar3基因的甲基化水平中位数为0.75(iqr＝0.54-0.84)，甲状腺恶性肿瘤bcar3基因的甲基化水平中位数为0.70(iqr＝0.52-0.79)，不同亚型bcar3基因的甲基化水平中位数依次为：乳头状癌0.68(iqr＝0.51-0.77)、滤泡癌0.63(iqr＝0.48-0.76)、髓样癌0.58(iqr＝0.45-0.72)、未分化癌0.52(iqr＝0.37-0.68)，不同分期bcar3基因的甲基化水平中位数依次为：i期为0.73(iqr＝0.53-0.81)、ii期为0.66(iqr＝0.50-0.79)、ⅲ期为0.61(iqr＝0.46-0.77)、ⅳ期为0.55(iqr＝0.36-0.69)。通过比较分析几者之间的甲基化水平，结果发现甲状腺良性肿瘤bcar3基因中所有cpg位点的甲基化水平显著
高于甲状腺恶性肿瘤中bcar3基因中所有cpg位点的甲基化水平的(表5)。另外，本发明中列出的几类甲状腺癌亚型如：乳头状癌、滤泡癌、髓样癌和未分化癌组织中的bcar3基因中seq id no.1、seq id no.2和seq id no.3所示的dna片段上甲基化位点的甲基化水平依次呈下降趋势(表5)，与甲状腺良性肿瘤之间的差异性越来越明显。另外，随着甲状腺恶性肿瘤分期的增加组织中bcar3基因中seq id no.1、seq id no.2和seq id no.3所示的dna片段上甲基化位点的甲基化水平也呈现降低趋势(表5)，与甲状腺良性肿瘤之间差异性也越来越明显。
[0193]
表5、甲状腺良性肿瘤和甲状腺恶性肿瘤及其各亚型、各分期的bcar3基因甲基化水平
[0194][0195][0196]
注：表中的cpg位点均指可区分的cpg位点。
[0197]
2、肿瘤组织中bcar3基因甲基化水平可以区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤
[0198]
通过比较分析380例甲状腺良性肿瘤病例和598例甲状腺恶性肿瘤病例的bcar3甲基化水平，结果发现，甲状腺恶性肿瘤、甲状腺乳头状癌、甲状腺滤泡癌、甲状腺髓样癌和甲状腺未分化癌患者中bcar3_a片段、bcar3_b片段和bcar3_c片段甲基化水平显著低于甲状腺良性肿瘤患者中对应片段的甲基化水平(p＜0.05)。具体结果见表6。
[0199]
表6、甲状腺良性肿瘤与不同亚型甲状腺恶性肿瘤之间的bcar3基因甲基化水平差异
[0200][0201][0202]
注：表中的cpg位点均指可区分的cpg位点。
[0203]
3、肿瘤组织中bcar3基因甲基化水平可以区分不同亚型的甲状腺恶性肿瘤
[0204]
通过比较分析不同分型的甲状腺恶性肿瘤(380例甲状腺乳头状癌、138例甲状腺滤泡癌、44例甲状腺髓样癌和36例甲状腺未分化癌)病例的bcar3甲基化水平，结果发现，甲状腺乳头状癌、甲状腺滤泡癌、甲状腺髓样癌和甲状腺未分化癌患者的bcar3基因甲基化水平之间有显著差异(p＜0.05)。具体结果见表7。
[0205]
表7、甲状腺恶性肿瘤各亚型之间的bcar3基因甲基化水平差异
[0206][0207][0208]
注：表中的cpg位点均指可区分的cpg位点。
[0209]
4、肿瘤组织中bcar3基因甲基化水平可以区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤
[0210]
通过比较分析380例甲状腺良性肿瘤病例和不同分期的甲状腺恶性肿瘤患者(470例ⅰ期患者、68例ⅱ期患者、24例ⅲ期患者和36例ⅳ期患者)的bcar3甲基化水平，结果发现，bcar3_a片段、bcar3_b片段和bcar3_c片段甲基化水平不同分期与甲状腺良性肿瘤之间存在显著差异(p＜0.01)。具体结果见表8。
[0211]
表8、甲状腺良性肿瘤与不同分期甲状腺恶性肿瘤之间的bcar3基因甲基化水平差异
[0212][0213]
注：表中的cpg位点均指可区分的cpg位点。
[0214]
5、肿瘤组织中bcar3基因甲基化水平可以区分不同分期的甲状腺恶性肿瘤
[0215]
通过比较分析不同分期的甲状腺恶性肿瘤患者(470例ⅰ期患者、68例ⅱ期患者、24例ⅲ期患者和36例ⅳ期患者)的bcar3甲基化水平，结果发现，ⅰ期甲状腺恶性肿瘤、ⅱ期甲状腺恶性肿瘤、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者的bcar3基因甲基化水平之间有显著差异(p＜0.05)。具体结果见表9。
[0216]
表9、甲状腺恶性肿瘤不同分期之间的bcar3基因甲基化水平差异
[0217][0218]
注：表中的cpg位点均指可区分的cpg位点。
[0219]
6、bcar3基因甲基化用于辅助癌症诊断的数学模型的建立
[0220]
本发明建立的数学模型可以用于达到如下目的：
[0221]
(1)区分甲状腺恶性肿瘤和甲状腺良性肿瘤；
[0222]
(2)区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；
[0223]
(3)区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；
[0224]
(4)区分甲状腺恶性肿瘤不同亚型；
[0225]
(5)区分甲状腺恶性肿瘤不同分期。
[0226]
数学模型的建立方法如下：
[0227]
(a)数据来源：步骤一中列出380例甲状腺良性肿瘤和598例甲状腺恶性肿瘤(甲状腺乳头状癌380例，甲状腺滤泡癌138例，甲状腺髓样癌44例和甲状腺未分化癌36例)的组织样本的目标cpg位点(表1-表3中的一种或多种的组合)甲基化水平(检测方法同步骤二)。
[0228]
(b)模型建立
[0229]
根据需要选取任意两类不同类型患者数据即训练集(例如：甲状腺良性肿瘤和甲状腺恶性肿瘤患者、甲状腺良性肿瘤和甲状腺乳头状癌患者、甲状腺良性肿瘤和甲状腺滤泡癌患者、甲状腺良性肿瘤和甲状腺髓样癌患者、甲状腺良性肿瘤和甲状腺未分化癌患者、
甲状腺乳头状癌和甲状腺滤泡癌患者、甲状腺乳头状癌和甲状腺髓样癌患者、甲状腺乳头状癌和甲状腺未分化癌患者、甲状腺滤泡癌和甲状腺髓样癌患者、甲状腺滤泡癌和甲状腺未分化癌患者、甲状腺髓样癌和甲状腺未分化癌患者、甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者)作为用于建立模型的数据，使用sas，r，spss等统计软件使用二分类逻辑回归的统计方法通过公式建立数学模型。数学模型公式计算出的最大约登指数对应的数值为阈值或直接设定0.5为阈值，待测样品经过测试和代入模型计算后得到的检测指数大于阈值归为一类(b类)，小于阈值归为另外一类(a类)，等于阈值作为不确定的灰区。在对新的待测样品进行预测来判断属于哪一类时，首先通过dna甲基化的测定方法检测该待测样品bcar3基因上一个或者多个cpg位点的甲基化水平，然后将这些甲基化水平的数据代入上述数学模型，计算得到所述待测样本对应的检测指数，然后比较所述待测样本对应的检测指数和阈值的大小，根据比较结果确定所述待测样本属于哪一类样本。
[0230]
举例：如图1所示，将训练集中bcar3基因单个cpg位点的甲基化水平或者多个cpg位点组合的甲基化水平的数据通过sas、r、spss等统计软件使用二分类逻辑回归的公式建立用于区分a类和b类的数学模型。该数学模型在此为二类逻辑回归模型，具体为：ln(y/(1-y))＝b0+b1x1+b2x2+b3x3+
…
+bnxn，其中y为因变量即将待测样品的一个或者多个甲基化位点的甲基化水平代入模型以后得出的检测指数，b0为常量，x1-xn为自变量即为该测试样品的一个或者多个甲基化位点的甲基化水平(每一个值为0-1之间的数值)，b1-bn为模型赋予每一个位点甲基化水平的权重。具体应用时，先根据训练集中已经检测的样本的一个或者多个dna甲基化位点的甲基化水平(x1-xn)及其已知的分类情况(a类或者b类，分别对y赋值0和1)建立数学模型，由此确定该数学模型的常量b0以及各个甲基化位点的权重b1-bn，并由该数学模型计算出的以最大约登指数对应的数值为阈值或直接设定0.5为划分的阈值。待测样品经过测试和代入模型计算后得到的检测指数即y值大于阈值归为b类，小于阈值归为a类，等于阈值作为不确定的灰区。其中a类和b类为相对应的两分类(二分类的分组，哪一组a类，哪一组是b类，要根据具体的数学模型来确定，在此不做约定)。对受试者的样品进行预测来判断属于哪一类时，首先收集受试者的活检样本(即肿瘤组织)，然后从中提取dna。将提取的dna通过重亚硫酸盐转化后，用dna甲基化的测定方法对受试者bcar3基因的单个cpg位点的甲基化水平或者多个cpg位点组合的甲基化水平进行检测，然后将检测得到的甲基化数据代入上述数学模型。如果该受试者的bcar3基因一个或者多个cpg位点的甲基化水平代入上述数学模型后计算出来的值即检测指数大于阈值，则该受试者判定与训练集中检测指数大于阈值的归属一类(b类)；如果该受试者的bcar3基因一个或者多个cpg位点的甲基化水平数据代入上述数学模型后计算出来的值即检测指数小于阈值，则该受试者跟训练集中检测指数小于阈值的归属一类(a类)；如果该受试者的bcar3基因一个或者多个cpg位点的甲基化水平数据代入上述数学模型后计算出来的值即检测指数等于阈值，则不能判断该受试者是a类还是b类。
[0231]
举例：如图2所示，举例说明bcar3_b的如下5个可区分cpg位点(bcar3_b_2、bcar3_b_3、bcar3_b_4.5.6.7、bcar3_b_8、bcar3_b_9)的甲基化以及数学建模在鉴别甲状腺良性肿瘤和甲状腺恶性肿瘤组织中的应用：将甲状腺良性肿瘤和甲状腺恶性肿瘤患者训练集(在此为：380例甲状腺良性肿瘤和598例甲状腺恶性肿瘤患者)中已经检测的bcar3_b的上述5个可区分cpg位点的甲基化水平的数据通过spss软件或r软件使用二分类逻辑回归的公式建立用于鉴别甲状腺恶性肿瘤患者的数学模型。该数学模型在此为二类逻辑回归模型，由此确定该数学模型的常量b0以及各个甲基化位点的权重b1-bn，在此例中具体为：ln(y/(1-y))＝3.793-1.197*bcar3_b_2
ꢀ‑
6.621*bcar3_b_3+0.790*bcar3_b_4.5.6.7+0.128*bcar3_b_8+2.105*bcar3_b_9，其中y为因变量，即将待测样品的bcar3_b的上述5个可区分cpg位点的甲基化水平代入模型以后经换算得出的检测指数。bcar3_b_4、bcar3_b_5、bcar3_b_6和bcar3_b_7位于同一片段，故以bcar3_b_4.5.6.7代表这两个位点甲基化水平的平均值。通过设定0.5为阈值的情况下，待测样品的bcar3_b的上述5个可区分cpg位点的甲基化水平经过测试得到的值代入模型进行计算，得到的检测指数即y值小于阈值归为甲状腺良性肿瘤患者，大于阈值归为甲状腺恶性肿瘤患者，等于阈值则不确定为甲状腺良性肿瘤患者还是甲状腺恶性肿瘤患者。此模型的曲线下面积(auc)计算结果为0.81(表14)。具体受试者判断方法举例如下所示，从两位受试者(甲，乙)分别收集活检样本提取dna，将提取的dna通过重亚硫酸盐转化后，用dna甲基化的测定方法对受试者的(bcar3_b_2、bcar3_b_3、bcar3_b_4.5.6.7、bcar3_b_8、bcar3_b_9)这5个可区分cpg位点的甲基化水平进行检测。然后将检测得到的甲基化水平数据信息代入上述数学模型。甲受试者的bcar3_b的上述5个可区分cpg位点的甲基化水平数据代入上述数学模型后计算出来的值为0.81大于0.5，则甲受试者判定为甲状腺恶性肿瘤患者(与临床诊断相符)；乙受试者的bcar3_b的上述5个可区分cpg位点的甲基化水平数据代入上述数学模型后计算出来的值为0.40小于0.5，则乙受试者判定甲状腺良性肿瘤患者(与临床诊断相符)。
[0232]
(c)模型效果评价
[0233]
根据上述方法，分别建立用于发现甲状腺良性肿瘤和甲状腺恶性肿瘤患者、甲状腺良性肿瘤和甲状腺乳头状癌患者、甲状腺良性肿瘤和甲状腺滤泡癌患者、甲状腺良性肿瘤和甲状腺髓样癌患者、甲状腺良性肿瘤和甲状腺未分化癌患者、甲状腺乳头状癌和甲状腺滤泡癌患者、甲状腺乳头状癌和甲状腺髓样癌患者、甲状腺乳头状癌和甲状腺未分化癌患者、甲状腺滤泡癌和甲状腺髓样癌患者、甲状腺滤泡癌和甲状腺未分化癌患者、甲状腺髓样癌和甲状腺未分化癌患者、甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者的数学模型，并且通过受试者曲线(roc曲线)对其有效性进行评价。roc曲线得出的曲线下面积(auc)越大，说明模型的区分度越好，分子标志物越有效。采用不同cpg位点进行数学模型构建后的评价结果如表10、表11、表12和表13所示。表10、表11、表12和表13中，1个cpg位点代表bcar3_b扩增片段中任意一个cpg位点的位点，2个cpg位点代表bcar3_b扩增片段中任意2个cpg位点的组合，3个
cpg位点代表bcar3_b扩增片段中任意3个cpg位点的组合，
……
以此类推。表中的数值为不同位点组合评价结果的范围值(即任意个cpg位点组合方式的结果均在此范围内)。
[0234]
上述研究结果显示，bcar3基因甲基化对于各组的鉴别能力(甲状腺良性肿瘤和甲状腺恶性肿瘤患者、甲状腺良性肿瘤和甲状腺乳头状癌患者、甲状腺良性肿瘤和甲状腺滤泡癌患者、甲状腺良性肿瘤和甲状腺髓样癌患者、甲状腺良性肿瘤和甲状腺未分化癌患者、甲状腺乳头状癌和甲状腺滤泡癌患者、甲状腺乳头状癌和甲状腺髓样癌患者、甲状腺乳头状癌和甲状腺未分化癌患者、甲状腺滤泡癌和甲状腺髓样癌患者、甲状腺滤泡癌和甲状腺未分化癌患者、甲状腺髓样癌和甲状腺未分化癌患者、甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者)随着bcar3基因上甲基化位点个数的增加而增加。
[0235]
除此以外，在表1-表3所示的cpg位点中，还存在少数几个优选位点的组合比多个非优选位点组合的鉴别能力更好的情况。例如表14、表15、表16和表17所示的bcar3_b_2、bcar3_b_3、bcar3_b_4.5.6.7、bcar3_b_8、bcar3_b_9这5个可区分cpg位点的组合是bcar3_b中任意5个组合的优选位点。
[0236]
综上所述，bcar3基因上的cpg位点及其各种组合，bcar3_a片段上的cpg位点及其各种组合，bcar3_b片段上的cpg位点及其各种组合，bcar3_b_2、bcar3_b_3、bcar3_b_4.5.6.7、bcar3_b_8、bcar3_b_9位点及其各种组合，bcar3_c片段上的cpg位点及其各种组合，以及bcar3_a、bcar3_b和bcar3_c上的cpg位点及其各种组合的甲基化水平对甲状腺良性肿瘤和甲状腺恶性肿瘤患者、甲状腺良性肿瘤和甲状腺乳头状癌患者、甲状腺良性肿瘤和甲状腺滤泡癌患者、甲状腺良性肿瘤和甲状腺髓样癌患者、甲状腺良性肿瘤和甲状腺未分化癌患者、甲状腺乳头状癌和甲状腺滤泡癌患者、甲状腺乳头状癌和甲状腺髓样癌患者、甲状腺乳头状癌和甲状腺未分化癌患者、甲状腺滤泡癌和甲状腺髓样癌患者、甲状腺滤泡癌和甲状腺未分化癌患者、甲状腺髓样癌和甲状腺未分化癌患者、甲状腺良性肿瘤和ⅰ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅱ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅲ期甲状腺恶性肿瘤患者、甲状腺良性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅱ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅰ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅲ期甲状腺恶性肿瘤患者、ⅱ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者、ⅲ期甲状腺恶性肿瘤和ⅳ期甲状腺恶性肿瘤患者具有判别能力。
[0237]
表10、bcar3_b的cpg位点及其组合用于区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤
[0238][0239]
表11、bcar3_b的cpg位点及其组合用于区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤
[0240]
[0241][0242]
表12、bcar3_b的cpg位点及其组合用于区分不同亚型甲状腺恶性肿瘤患者
[0243][0244]
表13、bcar3_b的cpg位点及其组合用于区分不同分期甲状腺恶性肿瘤患者
[0245][0246]
表14、bcar3_b的最佳cpg位点及其组合用于区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤
[0247][0248]
注：表中的cpg位点均指可区分的cpg位点。
[0249]
表15、bcar3_b的最佳cpg位点及其组合用于区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤
[0250][0251]
注：表中的cpg位点均指可区分的cpg位点。
[0252]
表16、bcar3_b的最佳cpg位点及其组合用于甲状腺恶性肿瘤不同亚型之间的区分
[0253][0254]
注：表中的cpg位点均指可区分的cpg位点。
[0255]
表17、bcar3_b的最佳cpg位点及其组合用于甲状腺恶性肿瘤不同分期之间的区分
[0256][0257]
注：表中的cpg位点均指可区分的cpg位点。
[0258]
以上对本发明进行了详述。对于本领域技术人员来说，在不脱离本发明的宗旨和范围，以及无需进行不必要的实验情况下，可在等同参数、浓度和条件下，在较宽范围内实施本发明。虽然本发明给出了特殊的实施例，应该理解为，可以对本发明作进一步的改进。总之，按本发明的原理，本技术欲包括任何变更、用途或对本发明的改进，包括脱离了本技术中已公开范围，而用本领域已知的常规技术进行的改变。

技术特征：
1.一种试剂盒，包括用于检测bcar3基因甲基化水平的物质；所述试剂盒的用途为如下中的至少一种：(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；(5)区分或辅助区分甲状腺恶性肿瘤不同分期。2.根据权利要求1所述的试剂盒，其特征在于：所述bcar3基因甲基化水平为bcar3基因中如下(a1)-(a3)所示片段中的全部或部分cpg位点的甲基化水平：(a1)seq id no.1所示的dna片段或与其具有80％以上同一性的dna片段；(a2)seq id no.2所示的dna片段或与其具有80％以上同一性的dna片段；(a3)seq id no.3所示的dna片段或与其具有80％以上同一性的dna片段。3.根据权利要求1或2所述的试剂盒，其特征在于：所述全部或部分cpg位点为如下任一所示cpg位点：(b1)bcar3基因中seq id no.1、seq id no.2和seq id no.3所示3个dna片段中的任意一个或多个cpg位点；(b2)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点和seq id no.2所示的dna片段上的所有cpg位点；(b3)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点和seq id no.3所示的dna片段上的所有cpg位点；(b4)bcar3基因中seq id no.2所示的dna片段上的所有cpg位点和seq id no.3所示的dna片段上的所有cpg位点；(b5)bcar3基因中seq id no.1所示的dna片段上的所有cpg位点、seq id no.2所示的dna片段上的所有cpg位点和seq id no.3所示的dna片段上的所有cpg位点；(b6)bcar3基因中seq id no.2所示的dna片段中的全部cpg位点或任意10个或任意9个或任意8个或任意7个或任意6个或任意5个或任意4个或任意3个或任意2个或任意1个cpg位点；(b7)bcar3基因中seq id no.2所示的dna片段上如下5项所示cpg位点的全部或任意4项或任意3项或任意2项或任意1项：第1项：seq id no.2所示的dna片段自5’端第117-118位所示cpg位点；第2项：seq id no.2所示的dna片段自5’端第329-330位所示cpg位点；第3项：seq id no.2所示的dna片段自5’端第353-354位和355-356位和第357-358位和第359-360位所示cpg位点；第4项：seq id no.2所示的dna片段自5’端第384-385位所示cpg位点；第5项：seq id no.2所示的dna片段自5’端第581-582位所示cpg位点。4.根据权利要求1-3中任一所述的试剂盒，其特征在于：所述用于检测bcar3基因甲基化水平的物质包含用于扩增bcar3基因全长或部分片段的引物组合；进一步地，所述部分片段为如下中至少一个片段：(g1)seq id no.1所示的dna片段或其包含的dna片段；
(g2)seq id no.2所示的dna片段或其包含的dna片段；(g3)seq id no.3所示的dna片段或其包含的dna片段；(g4)与seq id no.1所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段；(g5)与seq id no.2所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段；(g6)与seq id no.3所示的dna片段或其包含的dna片段具有80％以上同一性的dna片段；更进一步地，所述引物组合为引物对a和/或引物对b和/或引物对c；所述引物对a为引物a1和引物a2组成的引物对；所述引物a1为seq id no.4或seq id no.4的第11-35位核苷酸所示的单链dna；所述引物a2为seq id no.5或seq id no.5的第32-56位核苷酸所示的单链dna；所述引物对b为引物b1和引物b2组成的引物对；所述引物b1为seq id no.6或seq id no.6的第11-35位核苷酸所示的单链dna；所述引物b2为seq id no.7或seq id no.7的第32-56位核苷酸所示的单链dna；所述引物对c为引物c1和引物c2组成的引物对；所述引物c1为seq id no.8或seq id no.8的第11-35位核苷酸所示的单链dna；所述引物c2为seq id no.9或seq id no.9的第32-56位核苷酸所示的单链dna。5.根据权利要求1-4中任一所述的试剂盒，其特征在于：所述试剂盒中还含有储存有数学模型和/或数学模型使用方法的介质；所述数学模型按照包括如下步骤的方法获得：(c1)分别检测n1个a类型样本和n2个b类型样本的基因甲基化水平；(c2)取步骤(c1)获得的所有样本的基因甲基化水平数据，按照a类型和b类型的分类方式，通过二分类逻辑回归法建立数学模型，确定分类判定的阈值；所述数学模型使用方法包括如下步骤：(d1)检测待测样本的基因甲基化水平；(d2)将步骤(d1)获得的所述待测样本的基因甲基化水平数据代入所述数学模型，得到检测指数；然后比较检测指数和阈值的大小，根据比较结果确定所述待测样本的类型是a类型还是b类型；所述a类型样本和所述b类型样本为如下任一种：(e1)甲状腺良性肿瘤和甲状腺恶性肿瘤；(e2)甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(e3)甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(e4)不同亚型的甲状腺恶性肿瘤；(e5)不同分期的甲状腺恶性肿瘤。6.一种系统，包括：(f1)用于检测bcar3基因甲基化水平的试剂和/或仪器；(f2)装置，所述装置包括单元x和单元y；所述单元x用于建立数学模型，包括数据采集模块、数据分析处理模块和模型输出模
块；所述数据采集模块被配置为采集(f1)检测得到的n1个a类型样本和n2个b类型样本的bcar3基因甲基化水平数据；所述数据分析处理模块被配置为接收来自于所述数据采集模块的所述n1个a类型样本和n2个b类型样本的bcar3基因甲基化水平数据，按照a类型和b类型的分类方式，通过二分类逻辑回归法建立数学模型，确定分类判定的阈值；所述模型输出模块被配置为接收来自于所述数据分析处理模块建立的所述数学模型，并进行输出；所述单元y用于确定待测样本类型，包括数据输入模块、数据运算模块、数据比较模块和结论输出模块；所述数据输入模块被配置为输入(f1)检测得到的待测者的bcar3基因甲基化水平数据；所述数据运算模块被配置为接收来自于所述数据输入模块的所述待测者的bcar3基因甲基化水平数据，并将所述待测者的bcar3基因甲基化水平数据代入所述单元x中的所述数据分析处理模块建立的所述数学模型，计算得到检测指数；所述数据比较模块被配置为接收来自于所述数据运算模块计算得到的检测指数，并将所述检测指数与所述单元x中的所述数据分析处理模块中确定的所述阈值进行比较；所述结论输出模块被配置为接收来自于所述数据比较模块的比较结果，并根据所述比较结果输出所述待测样本的类型是a类型还是b类型的结论；所述a类型样本和所述b类型样本为如下任一种：(e1)甲状腺良性肿瘤和甲状腺恶性肿瘤；(e2)甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(e3)甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(e4)不同亚型的甲状腺恶性肿瘤；(e5)不同分期的甲状腺恶性肿瘤。7.权利要求1-4中任一所述用于检测bcar3基因甲基化水平的物质在制备产品中的应用；所述产品的用途为如下至少一种：(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；(5)区分或辅助区分甲状腺恶性肿瘤不同分期。8.权利要求1-4中任一所述用于检测bcar3基因甲基化水平的物质和权利要求5中所述储存有数学模型和/或数学模型使用方法的介质在制备产品中的应用；所述产品的用途为如下至少一种：(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；
(5)区分或辅助区分甲状腺恶性肿瘤不同分期。9.权利要求5中所述储存有数学模型和/或数学模型使用方法的介质在制备产品中的应用；所述产品的用途为如下至少一种：(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；(5)区分或辅助区分甲状腺恶性肿瘤不同分期。10.甲基化bcar3基因作为标志物在制备产品中的应用；所述产品的用途为如下至少一种：(1)区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤；(2)区分或辅助区分甲状腺良性肿瘤和不同亚型的甲状腺恶性肿瘤；(3)区分或辅助区分甲状腺良性肿瘤和不同分期的甲状腺恶性肿瘤；(4)区分或辅助区分甲状腺恶性肿瘤不同亚型；(5)区分或辅助区分甲状腺恶性肿瘤不同分期。

技术总结
本发明公开了一种甲状腺癌早期辅助诊断标志物。本发明要求保护一种试剂盒，包括用于检测BCAR3基因甲基化水平的物质；所述试剂盒可用于区分或辅助区分甲状腺良性肿瘤和甲状腺恶性肿瘤、甲状腺良性肿瘤和不同亚型/不同分期的甲状腺恶性肿瘤；区分或辅助区分甲状腺恶性肿瘤不同亚型/不同分期。本发明揭示了相较于甲状腺良性肿瘤，甲状腺癌患者组织中BCAR3基因的低甲基化现象，本发明对提高甲状腺癌早期诊疗效果和降低甲状腺癌死亡率，以及指导制定合理的临床治疗方案均有重要的科学意义和临床应用价值。意义和临床应用价值。


技术研发人员：狄飞飞 王俊 袁红丽
受保护的技术使用者：南京腾辰生物科技有限公司
技术研发日：2023.03.29
技术公布日：2023/8/4