CN113528658A - 甲状腺肿瘤良恶性鉴别模型及其应用 - Google Patents

Abstract

本发明涉及用于甲状腺肿瘤良恶性鉴别和恶性潜能未定的滤泡性肿瘤评估的甲基化标志物、试剂盒、装置、计算机可读介质和系统。

Description

甲状腺肿瘤良恶性鉴别模型及其应用

技术领域

本发明涉及分子生物学技术和计算机领域，涉及滤泡状甲状腺肿瘤良恶性鉴别模型及其构建方法和应用。

背景技术

甲状腺肿瘤是一种十分常见的肿瘤，有良性和恶性之分。良性肿瘤分为滤泡型腺瘤(follicular thyroid adenoma,FTA)和乳头型腺瘤(papillary thyroid adenoma,PTA)两种，其中FTA是最常见的良性肿瘤，约占甲状腺腺瘤的70％～80％。恶性肿瘤分为乳头状癌(papillary thyroid carcinoma,PTC)，滤泡癌(follicular thyroid carcinoma,FTC)，未分化癌(anaplastic thyroid carcinoma)和髓样癌(medullary thyroidcarcinoma)。后两类较为少见，PTC的恶性程度低，预后较好，而FTC是一种恶性程度偏高的恶性肿瘤。除此之外，临床上还具有一类良恶性之间的交界性肿瘤，即：形态上具有非浸润性和可疑恶性肿瘤细胞核特征的包裹性滤泡型肿瘤，或者具有可疑包膜浸润但不伴有细胞核特征的肿瘤。这类交界性肿瘤被命名为：恶性潜能未定的甲状腺肿瘤(UMP)。

临床上常根据病理特征对甲状腺肿瘤进行良性、恶性以及UMP的划分。然而FTA和FTC的病理特征存在交叉，即部分FTA和FTC不能完全区分；此外，病理上的UMP中存在部分患者的肿瘤可能是被漏诊的恶性甲状腺癌。常见的恶性肿瘤标志物并不具备区分良恶性甲状腺癌的功能。一些基因突变，比如RAS和BRAF的突变虽然与肿瘤恶性程度和不良预后相关，但是其在良恶性肿瘤鉴定中的灵敏性和特异性都较低；而TERT启动子区突变具有很高的恶性肿瘤特异性，但其在甲状腺癌中的突变率很低，PTC中仅11.7％，FTC中仅13.9％。因此，临床上还需要一种简便高效的模型和算法，辅助甲状腺癌的良恶性判断，以及UMP的恶性潜能评估。

发明内容

针对现有技术中滤泡型甲状腺肿瘤良恶性鉴别的诊断操作复杂、价格昂贵以及针对恶性潜能未定的甲状腺肿瘤诊断难度大等不足，本发明提供了一种适合于滤泡型甲状腺肿瘤诊断的模型和算法，所述模型涉及的标志物对于滤泡型甲状腺肿瘤良恶性鉴别以及对恶性潜能未定的甲状腺肿瘤的评估都具有较好的灵敏度和特异度，对于及时诊断、改善甲状腺肿瘤预后、降低死亡率有很重要的意义。

本发明还提供上述标志物的筛选方法，通过该方法所得的标志物对于滤泡型甲状腺肿瘤都具有很好的灵敏度和特异性，对于甲状腺肿瘤的治疗有重要的临床意义。

本发明还提供鉴别滤泡型甲状腺肿瘤良恶性的模型及其构建方法，该模型构建方法简单，方便快捷，对于良性、恶性甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤都具有较好的灵敏度和特异度。

因此，本发明提供一种甲基化标志物，所述标志物选自分离自动物的以下任意一个、多个或全部基因的基因间区、内含子、外显子、启动子或UTR，或包含其的基因组片段，或与其具有70％序列相同性的变体，所述片段的长度为5-1500bp：TBX15、LHX4、MYOG、HMX3、KRT85、GPC5、FOXA1、SIX6、GSE1、WNK4、HOXB3、GATA6-AS1、CBLN2、POU3F3、KIAA2012、LRRTM1、FAM95A、SIM2、PITX2、CRIPAK、RBM47、TMEM174、FRMD1、MYL10、MIR153-2、LINC00689、HOXA5、MYO1G、FIGNL1、SOX17、AJM1、MIR3675、KAZALD1、NEURL1、LINC00173、AQP5、KRT86、FRY、KCTD12、LINC00239、ZC3H18、RFNG、TNPO2、LYL1、EPHX3、ADAMTS10、LINC01270、OGFRP1、LINC00887、LINC00884、DDAH2、SNX10、MRPL23-AS1、MYEF2、RORA、NUDT16L1、LINC01530、PROC、PRMT2、LGALS1。在一个或多个实施方案中，所述标记物是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

本发明提供一种甲基化标志物，所述标志物选自分离自动物的以下任意一个、多个或全部：TBX15的基因间区、LHX4的内含子、MYOG的外显子、MYOG的3’UTR、HMX3的启动子、KRT85的基因间区、GPC5的启动子、FOXA1的基因间区、FOXA1的外显子、SIX6的启动子、GSE1的基因间区、WNK4的启动子、HOXB3的启动子、GATA6-AS1的启动子、CBLN2的启动子、POU3F3的启动子、KIAA2012的内含子、LRRTM1的启动子、FAM95A的基因间区、SIM2的内含子、PITX2的内含子、CRIPAK的基因间区、RBM47的启动子、TMEM174的基因间区、FRMD1的基因间区、MYL10的基因间区、MIR153-2的内含子、LINC00689的启动子、HOXA5的启动子、MYO1G的3’UTR、FIGNL1的外显子、SOX17的基因间区、AJM1的启动子、MIR3675的外显子、KAZALD1的启动子、NEURL1的启动子、LINC00173的基因间区、AQP5的启动子、KRT86的内含子、FRY的启动子、KCTD12的启动子、LINC00239的基因间区、ZC3H18的外显子、RFNG的启动子、TNPO2的启动子、LYL1的启动子、EPHX3的启动子、ADAMTS10的启动子、LINC01270的基因间区、OGFRP1的基因间区、LINC00887的基因间区、LINC00884的启动子、DDAH2的启动子、SNX10的基因间区、MRPL23-AS1的基因间区、MYEF2的启动子、RORA的启动子、NUDT16L1的外显子、LINC01530的基因间区、PROC的基因间区、PRMT2的基因间区、LGALS1的启动子，或包含其的基因组片段，或与其具有70％序列相同性的变体，所述片段的长度为5-1500bp。在一个或多个实施方案中，所述标记物是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

本发明提供一种甲基化标志物，所述标志物选自分离自动物的以下任意一个、多个或全部：TBX15的基因间区、LHX4的第1个内含子、MYOG的第34个外显子、MYOG的3’UTR、HMX3的启动子、KRT85的基因间区、GPC5的启动子、FOXA1的基因间区、FOXA1的的第2个外显子、SIX6的启动子、GSE1的基因间区、WNK4的启动子、HOXB3的启动子、GATA6-AS1的启动子、CBLN2的启动子、POU3F3的启动子、KIAA2012的第15个内含子、LRRTM1的启动子、FAM95A的基因间区、SIM2的第1个内含子、PITX2的第2个内含子、CRIPAK的基因间区、RBM47的启动子、TMEM174的基因间区、FRMD1的基因间区、MYL10的基因间区、MIR153-2的第12个内含子、LINC00689的启动子、HOXA5的启动子、MYO1G的3’UTR、FIGNL1的第8个外显子、SOX17的基因间区、AJM1的启动子、MIR3675的第1外显子、KAZALD1的启动子、NEURL1的启动子、LINC00173的基因间区、AQP5的启动子、KRT86的第1个内含子、FRY的启动子、KCTD12的启动子、LINC00239的基因间区、ZC3H18的的第10个外显子、RFNG的启动子、TNPO2的启动子、LYL1的启动子、EPHX3的启动子、ADAMTS10的启动子、LINC01270的基因间区、OGFRP1的基因间区、LINC00887的基因间区、LINC00884的启动子、DDAH2的启动子、SNX10的基因间区、MRPL23-AS1的基因间区、MYEF2的启动子、RORA的启动子、NUDT16L1的第7个外显子、LINC01530的基因间区、PROC的基因间区、PRMT2的基因间区、LGALS1的启动子，或包含其的基因组片段，或与其具有70％序列相同性的变体，所述片段的长度为5-1500bp。在一个或多个实施方案中，所述标记物是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

本发明提供一种核酸分子，所述核酸分子的长度为5-1500bp、6-1400bp、7-1300bp、7-1261bp，所述核酸分子具有(1)动物基因组的包含选自以下一种、多种或全部序列的序列，或与其具有至少70％相同性的变体，所述变体中的甲基化位点未发生突变：chr1:119543868:119543879，chr1:180202465:180202541，chr1:203044677:203044823，chr1:203045227:203045293，chr10:124896784:124897154，chr12:52745149:52745173，chr13:92051635:92051674，chr14:38071275:38071410，chr14:38091345:38091744，chr14:60977845:60977865，chr16:85620063:85620331，chr17:40937184:40937480，chr17:46631705:46632377，chr18:19747115:19747127，chr18:70209347:70209475，chr2:105474371:105474381，chr2:203036052:203036236，chr2:80530770:80530826，chr2:95401381:95401629，chr21:38077572:38077685，chr4:111554966:111554989，chr4:1400384:1400459，chr4:40632502:40632519，chr5:72594262:72594269，chr6:168502296:168502347，chr7:101241802:101241926，chr7:157485865:157486031，chr7:158798629:158798674，chr7:27182264:27183525，chr7:45002063:45002754，chr7:50467662:50468309，chr8:55366694:55366747，chr9:139740669:139740676，chr1:16862044:16862199，chr10:102822370:102822647，chr10:105253375:105253477，chr12:116946248:116946304，chr12:50355085:50355570，chr12:52652220:52652362，chr13:32605660:32605843，chr13:77459521:77459792，chr14:102172655:102172687，chr16:88599725:88600068，chr17:80009015:80009025，chr19:12831808:12832195，chr19:13213485:13213513，chr19:13213644:13213814，chr19:15344092:15344411，chr19:8674674:8674749，chr20:48902548:48902611，chr22:42710260:42710349，chr3:193987426:193987681，chr3:194208192:194208617，chr6:31696240:31696334，chr7:26415826:26415917，chr11:2000109:2000154，chr15:48470425:48470556，chr15:60883371:60883395，chr16:4715122:4715220，chr19:52104749:52104928，chr2:128158537:128158621，chr21:48087183:48088183，chr22:38071168:38071189，chr7:26415938:26416562(2)(1)的互补序列。在上述一个或多个实施方案中，所述核酸分子是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

在一个或多个实施方案中，所述核酸分子的序列选自以下中的任意一种、一种以上或全部的组合：chr1:119543868:119543879，chr1:180202465:180202541，chr1:203044677:203044823，chr1:203045227:203045293，chr10:124896784:124897154，chr12:52745149:52745173，chr13:92051635:92051674，chr14:38071275:38071410，chr14:38091345:38091744，chr14:60977845:60977865，chr16:85620063:85620331，chr17:40937184:40937480，chr17:46631705:46632377，chr18:19747115:19747127，chr18:70209347:70209475，chr2:105474371:105474381，chr2:203036052:203036236，chr2:80530770:80530826，chr2:95401381:95401629，chr21:38077572:38077685，chr4:111554966:111554989，chr4:1400384:1400459，chr4:40632502:40632519，chr5:72594262:72594269，chr6:168502296:168502347，chr7:101241802:101241926，chr7:157485865:157486031，chr7:158798629:158798674，chr7:27182264:27183525，chr7:45002063:45002754，chr7:50467662:50468309，chr8:55366694:55366747，chr9:139740669:139740676，chr1:16862044:16862199，chr10:102822370:102822647，chr10:105253375:105253477，chr12:116946248:116946304，chr12:50355085:50355570，chr12:52652220:52652362，chr13:32605660:32605843，chr13:77459521:77459792，chr14:102172655:102172687，chr16:88599725:88600068，chr17:80009015:80009025，chr19:12831808:12832195，chr19:13213485:13213513，chr19:13213644:13213814，chr19:15344092:15344411，chr19:8674674:8674749，chr20:48902548:48902611，chr22:42710260:42710349，chr3:193987426:193987681，chr3:194208192:194208617，chr6:31696240:31696334，chr7:26415826:26415917，chr11:2000109:2000154，chr15:48470425:48470556，chr15:60883371:60883395，chr16:4715122:4715220，chr19:52104749:52104928，chr2:128158537:128158621，chr21:48087183:48088183，chr22:38071168:38071189，chr7:26415938:26416562。

本发明提供一种标志物，包含下述核酸分子或与其具有至少70％相同性的变体和算法的组合中的任意一种、一种以上或全部，所述变体中的甲基化位点未发生突变：chr1:119543868:119543879-mhl3、chr1:180202465:180202541-mhl3、chr1:203044677:203044823-mhl3、chr1:203045227:203045293-mhl3、chr10:124896784:124897154-mhl3、chr12:52745149:52745173-mhl3、chr13:92051635:92051674-mhl3、chr14:38071275:38071410-mhl3、chr14:38091345:38091744-mhl3、chr14:60977845:60977865-mhl3、chr16:85620063:85620331-mhl3、chr17:40937184:40937480-mhl3、chr17:46631705:46632377-mhl3、chr18:19747115:19747127-mhl3、chr18:70209347:70209475-mhl3、chr2:105474371:105474381-mhl3、chr2:203036052:203036236-mhl3、chr2:80530770:80530826-mhl3、chr2:95401381:95401629-mhl3、chr21:38077572:38077685-mhl3、chr4:111554966:111554989-mhl3、chr4:1400384:1400459-mhl3、chr4:40632502:40632519-mhl3、chr5:72594262:72594269-mhl3、chr6:168502296:168502347-mhl3、chr7:101241802:101241926-mhl3、chr7:157485865:157486031-mhl3、chr7:158798629:158798674-mhl3、chr7:27182264:27183525-mhl3、chr7:45002063:45002754-mhl3、chr7:50467662:50468309-mhl3、chr8:55366694:55366747-mhl3、chr9:139740669:139740676-mhl3、chr1:203044677:203044823-mhl、chr13:92051635:92051674-mhl、chr14:60977845:60977865-mhl、chr18:19747115:19747127-mhl、chr7:27182264:27183525-mhl、chr1:16862044:16862199-pdr、chr10:102822370:102822647-pdr、chr10:105253375:105253477-pdr、chr12:116946248:116946304-pdr、chr12:50355085:50355570-pdr、chr12:52652220:52652362-pdr、chr13:32605660:32605843-pdr、chr13:77459521:77459792-pdr、chr14:102172655:102172687-pdr、chr16:88599725:88600068-pdr、chr17:80009015:80009025-pdr、chr19:12831808:12832195-pdr、chr19:13213485:13213513-pdr、chr19:13213644:13213814-pdr、chr19:15344092:15344411-pdr、chr19:8674674:8674749-pdr、chr20:48902548:48902611-pdr、chr22:42710260:42710349-pdr、chr3:193987426:193987681-pdr、chr3:194208192:194208617-pdr、chr6:31696240:31696334-pdr、chr7:26415826:26415917-pdr、chr10:102822370:102822647-umhl、chr11:2000109:2000154-umhl、chr15:48470425:48470556-umhl、chr15:60883371:60883395-umhl、chr16:4715122:4715220-umhl、chr19:52104749:52104928-umhl、chr2:128158537:128158621-umhl、chr21:48087183:48088183-umhl、chr22:38071168:38071189-umhl、chr7:26415938:26416562-umhl。在一个或多个实施方案中，所述标记物是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

在一个或多个实施方案中，标志物是下述核酸分子或与其具有至少70％相同性的变体和算法的组合中的任意一种、一种以上或全部，所述变体中的甲基化位点未发生突变：chr18:19747115:19747127_mhl3，chr15:48470425:48470556_umhl，chr19:12831808:12832195_pdr，chr17:80009015:80009025_pdr，chr10:105253375:105253477_pdr，chr13:92051635:92051674_mhl3，chr7:101241802:101241926_mhl3，chr14:38091345:38091744_mhl3，chr16:4715122:4715220_umhl，chr16:85620063:85620331_mhl3，chr18:19747115:19747127_mhl，chr7:45002063:45002754_mhl3，chr19:15344092:15344411_pdr，chr12:50355085:50355570_pdr，chr10:102822370:102822647_pdr，chr7:27182264:27183525_mhl，chr2:128158537:128158621_umhl，chr2:80530770:80530826_mhl3，chr2:105474371:105474381_mhl3，chr1:119543868:119543879_mhl3，chr2:203036052:203036236_mhl3，chr8:55366694:55366747_mhl3，chr10:124896784:124897154_mhl3，chr1:203044677:203044823_mhl3，chr2:95401381:95401629_mhl3，chr7:158798629:158798674_mhl3，chr4:1400384:1400459_mhl3，chr1:16862044:16862199_pdr，chr19:13213485:13213513_pdr，chr19:52104749:52104928_umhl，chr7:26415826:26415917_pdr，chr1:203045227:203045293_mhl3，chr17:40937184:40937480_mhl3，chr14:102172655:102172687_pdr，chr22:38071168:38071189_umhl，chr9:139740669:139740676_mhl3，chr21:48087183:48088183_umhl，chr5:72594262:72594269_mhl3，chr4:111554966:111554989_mhl3，chr12:116946248:116946304_pdr，chr12:52652220:52652362_pdr，chr7:26415938:26416562_umhl，chr19:8674674:8674749_pdr，chr17:46631705:46632377_mhl3，chr14:60977845:60977865_mhl3，chr13:92051635:92051674_mhl，chr18:70209347:70209475_mhl3，chr20:48902548:48902611_pdr，chr6:168502296:168502347_mhl3，chr21:38077572:38077685_mhl3。

在一个或多个实施方案中，标志物是下述核酸分子或与其具有至少70％相同性的变体和算法的组合中的任意一种、一种以上或全部，所述变体中的甲基化位点未发生突变：chr14:38091345:38091744_mhl3，chr22:42710260:42710349_pdr，chr16:85620063:85620331_mhl3，chr16:4715122:4715220_umhl，chr18:19747115:19747127_mhl3，chr12:50355085:50355570_pdr，chr2:105474371:105474381_mhl3，chr21:38077572:38077685_mhl3，chr17:40937184:40937480_mhl3，chr3:194208192:194208617_pdr，chr2:203036052:203036236_mhl3，chr1:203044677:203044823_mhl，chr10:102822370:102822647_pdr，chr1:203044677:203044823_mhl3，chr7:27182264:27183525_mhl3，chr22:38071168:38071189_umhl，chr1:16862044:16862199_pdr，chr8:55366694:55366747_mhl3，chr14:60977845:60977865_mhl3，chr21:48087183:48088183_umhl，chr9:139740669:139740676_mhl3，chr20:48902548:48902611_pdr，chr1:203045227:203045293_mhl3，chr7:45002063:45002754_mhl3，chr10:105253375:105253477_pdr，chr13:92051635:92051674_mhl，chr14:60977845:60977865_mhl，chr7:158798629:158798674_mhl3，chr18:70209347:70209475_mhl3，chr7:27182264:27183525_mhl。

在一个或多个实施方案中，标志物是下述核酸分子或与其具有至少70％相同性的变体和算法的组合中的任意一种、一种以上或全部，所述变体中的甲基化位点未发生突变：chr22:42710260:42710349_pdr，chr17:80009015:80009025_pdr，chr8:55366694:55366747_mhl3，chr2:95401381:95401629_mhl3，chr1:203044677:203044823_mhl3，chr19:13213485:13213513_pdr，chr7:27182264:27183525_mhl3，chr7:101241802:101241926_mhl3，chr18:19747115:19747127_mhl3，chr15:60883371:60883395_umhl。

在上述一个或多个实施方案中，所述核酸分子长度为5-1500bp、6-1400bp、7-1300bp、7-1261bp。

在上述一个或多个实施方案中，所述核酸分子的序列编号参考hg19。

在上述一个或多个实施方案中，所述标记物或核酸分子是与甲状腺肿瘤良恶性相关的基因组DNA甲基化标志物。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

在上述一个或多个实施方案中，所述动物是哺乳动物，优选人。

在上述一个或多个实施方案中，所述甲基化位点是连续的CG。

在上述一个或多个实施方案中，所述序列包括DNA正义链或反义链。

在上述一个或多个实施方案中，所述核酸分子用作检测样品中相应序列的DNA甲基化水平的内标或对照。

本发明第二方面提供检测DNA甲基化的试剂，所述试剂检测样品中本文第一方面所述标志物或核酸分子中的一种或多种的DNA甲基化水平。

在一个或多个实施方案中，所述样品来自哺乳动物，优选人。

在一个或多个实施方案中，所述试剂是选自以下方法的一个或多个中所用的试剂：基于重亚硫酸盐转化的PCR(例如甲基化特异性PCR)、DNA测序(如亚硫酸氢盐测序、全基因组甲基化测序、简化甲基化测序)、甲基化敏感的限制性内切酶分析法、荧光定量法、甲基化敏感性高分辨率熔解曲线法、基于芯片的甲基化图谱分析、质谱(例如飞行质谱)。

优选地，所述试剂选自以下一种或多种：重亚硫酸盐及其衍生物、PCR缓冲液、聚合酶、dNTP、引物、探针、甲基化敏感或不敏感的限制性内切酶、酶切缓冲液、荧光染料、荧光淬灭剂、荧光报告剂、外切核酸酶、碱性磷酸酶、内标、对照物。

在一个或多个实施方案中，所述试剂包括引物。所述引物检测本文所述区域或序列的甲基化水平。在一个或多个实施方案中，所述引物可以是基因组测序的引物，例如全基因组测序引物或针对部分基因组的测序引物，也可以是用于扩增本文所述区域或序列的PCR引物或用于扩增区域中一个或多个甲基化标志物的PCR引物。

在一个或多个实施方案中，所述引物是利用简化甲基化测序法检测所述标志物的甲基化水平的引物或用于扩增一个或多个标志物的PCR引物。

在一个或多个实施方案中，所述试剂包括探针。所述探针的序列的5’端标记荧光报告基团，3’端标记淬灭基团。优选地，所述探针检测本文所述区域或序列的甲基化水平。

本发明还提供鉴别甲状腺肿瘤良恶性的试剂盒，包含本文所述的标志物和/或本文所述的试剂。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

本发明还提供检测DNA甲基化的试剂和任选的本文所述的标志物中的一种或多种在制备用于鉴定样品中甲状腺肿瘤良恶性的试剂盒中的用途，所述试剂检测样品中本文所述标志物中的一种或多种的DNA甲基化水平。在一个或多个实施方案中，甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

在用途的一个或多个实施方案中，所述样品来自哺乳动物，优选人。所述样品优选来自组织、细胞或者体液，例如甲状腺组织或血液。在一个或多个实施方案中，所述样品是甲状腺结节或肿瘤的活检物，优选是细针穿刺活检物。在一个或多个实施方案中，所述样品是血浆。

在用途的一个或多个实施方案中，所述样品包括基因组DNA或cfDNA。

在一个或多个实施方案中，所述试剂是选自以下方法的一个或多个中所用的试剂：基于重亚硫酸盐转化的PCR(例如甲基化特异性PCR)、DNA测序(如亚硫酸氢盐测序、全基因组甲基化测序、简化甲基化测序)、甲基化敏感的限制性内切酶分析法、荧光定量法、甲基化敏感性高分辨率熔解曲线法、基于芯片的甲基化图谱分析、质谱(例如飞行质谱)。

优选地，所述试剂选自以下一种或多种：重亚硫酸盐及其衍生物、PCR缓冲液、聚合酶、dNTP、引物、探针、甲基化敏感或不敏感的限制性内切酶、酶切缓冲液、荧光染料、荧光淬灭剂、荧光报告剂、外切核酸酶、碱性磷酸酶、内标、对照物。

在一个或多个实施方案中，所述试剂包括引物。所述引物检测本文所述区域或序列的甲基化水平。在一个或多个实施方案中，所述引物可以是基因组测序的引物，例如全基因组测序引物或针对部分基因组的测序引物，也可以是用于扩增本文所述区域或序列的PCR引物或用于扩增区域中一个或多个标志物的PCR引物。

在一个或多个实施方案中，所述引物是利用简化甲基化测序法检测所述标志物的甲基化水平的引物或用于扩增一个或多个标志物的PCR引物。

在一个或多个实施方案中，所述试剂包括探针。所述探针的序列的5’端标记荧光报告基团，3’端标记淬灭基团。优选地，所述探针检测本文所述区域或序列的甲基化水平。

本发明第四方面提供一种鉴别滤泡状甲状腺肿瘤良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能的方法，包括：

(1)获取样品中一种或多种标志物的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理各标志物的甲基化水平，

(3)通过构建模型利用步骤(2)的经处理的甲基化水平获得评分，

(4)根据评分鉴别甲状腺肿瘤良恶性或评估甲状腺肿瘤的恶性潜能。

在一个或多个实施方案中，获取标志物的甲基化水平包括(1)检测全基因组DNA的甲基化水平或包含所述标志物的基因组DNA的甲基化水平，并选取标志物的甲基化水平；或(2)检测所述标志物的甲基化水平。

在一个或多个实施方案中，所述样品来自健康对象或甲状腺结节对象，优选滤泡状甲状腺结节对象。

在一个或多个实施方案中，所述样品是甲状腺组织样品。

在一个或多个实施方案中，所述样品是甲状腺结节或肿瘤样品。

在一个或多个实施方案中，所述样品是良性(FTA)和/或恶性甲状腺肿瘤(FTC)组织样品；

在一个或多个实施方案中，mhl算法如下：

其中l为该甲基化标志物的长度；P(MHi)是从该标志物起始位置到i位点的区域中，由NGS检测到的完全甲基化的读数数量占到该区域所有读数数量的比例；w_i是从该标志物起始位置到i位点的区域的长度的比重，其为i。

在一个或多个实施方案中，mhl3算法如下：

其中l为该甲基化标志物的长度；P(MHi)是从该标志物起始位置到i位点的区域中，由NGS检测到的完全甲基化的读数数量占到该区域所有读数数量的比例；w_i是从该标志物起始位置到i位点的区域的长度的比重，其为i³。

在一个或多个实施方案中，umhl算法如下：

其中l为该甲基化标志物的长度；P(MHi)是从该标志物起始位置到i位点的区域中，由NGS检测到的完全未甲基化的读数数量占到该区域所有读数数量的比例；w_i是从该标志物起始位置到i位点的区域的长度的比重，其为i。

在一个或多个实施方案中，pdr算法如下：

PDR＝不一致读数/总读数

即标志物中既有甲基化又有未甲基化胞嘧啶的读数数量与所有读数数量的比值。

在一个或多个实施方案中，步骤(1)的标志物是核酸分子，所述核酸分子的长度为5-1500bp、6-1400bp、7-1300bp、7-1261bp，所述核酸分子具有(1)动物基因组的包含选自以下一种、多种或全部序列的序列，或与其具有至少70％相同性的变体，所述变体中的甲基化位点未发生突变：chr1:119543868:119543879，chr1:180202465:180202541，chr1:203044677:203044823，chr1:203045227:203045293，chr10:124896784:124897154，chr12:52745149:52745173，chr13:92051635:92051674，chr14:38071275:38071410，chr14:38091345:38091744，chr14:60977845:60977865，chr16:85620063:85620331，chr17:40937184:40937480，chr17:46631705:46632377，chr18:19747115:19747127，chr18:70209347:70209475，chr2:105474371:105474381，chr2:203036052:203036236，chr2:80530770:80530826，chr2:95401381:95401629，chr21:38077572:38077685，chr4:111554966:111554989，chr4:1400384:1400459，chr4:40632502:40632519，chr5:72594262:72594269，chr6:168502296:168502347，chr7:101241802:101241926，chr7:157485865:157486031，chr7:158798629:158798674，chr7:27182264:27183525，chr7:45002063:45002754，chr7:50467662:50468309，chr8:55366694:55366747，chr9:139740669:139740676，chr1:16862044:16862199，chr10:102822370:102822647，chr10:105253375:105253477，chr12:116946248:116946304，chr12:50355085:50355570，chr12:52652220:52652362，chr13:32605660:32605843，chr13:77459521:77459792，chr14:102172655:102172687，chr16:88599725:88600068，chr17:80009015:80009025，chr19:12831808:12832195，chr19:13213485:13213513，chr19:13213644:13213814，chr19:15344092:15344411，chr19:8674674:8674749，chr20:48902548:48902611，chr22:42710260:42710349，chr3:193987426:193987681，chr3:194208192:194208617，chr6:31696240:31696334，chr7:26415826:26415917，chr11:2000109:2000154，chr15:48470425:48470556，chr15:60883371:60883395，chr16:4715122:4715220，chr19:52104749:52104928，chr2:128158537:128158621，chr21:48087183:48088183，chr22:38071168:38071189，chr7:26415938:26416562(2)(1)的互补序列。

在一个或多个实施方案中，所述核酸分子的序列选自以下任意一种、一种以上或全部的组合或与其具有至少70％相同性的变体，所述变体中的甲基化位点未突变：chr1:119543868:119543879，chr1:180202465:180202541，chr1:203044677:203044823，chr1:203045227:203045293，chr10:124896784:124897154，chr12:52745149:52745173，chr13:92051635:92051674，chr14:38071275:38071410，chr14:38091345:38091744，chr14:60977845:60977865，chr16:85620063:85620331，chr17:40937184:40937480，chr17:46631705:46632377，chr18:19747115:19747127，chr18:70209347:70209475，chr2:105474371:105474381，chr2:203036052:203036236，chr2:80530770:80530826，chr2:95401381:95401629，chr21:38077572:38077685，chr4:111554966:111554989，chr4:1400384:1400459，chr4:40632502:40632519，chr5:72594262:72594269，chr6:168502296:168502347，chr7:101241802:101241926，chr7:157485865:157486031，chr7:158798629:158798674，chr7:27182264:27183525，chr7:45002063:45002754，chr7:50467662:50468309，chr8:55366694:55366747，chr9:139740669:139740676，chr1:16862044:16862199，chr10:102822370:102822647，chr10:105253375:105253477，chr12:116946248:116946304，chr12:50355085:50355570，chr12:52652220:52652362，chr13:32605660:32605843，chr13:77459521:77459792，chr14:102172655:102172687，chr16:88599725:88600068，chr17:80009015:80009025，chr19:12831808:12832195，chr19:13213485:13213513，chr19:13213644:13213814，chr19:15344092:15344411，chr19:8674674:8674749，chr20:48902548:48902611，chr22:42710260:42710349，chr3:193987426:193987681，chr3:194208192:194208617，chr6:31696240:31696334，chr7:26415826:26415917，chr11:2000109:2000154，chr15:48470425:48470556，chr15:60883371:60883395，chr16:4715122:4715220，chr19:52104749:52104928，chr2:128158537:128158621，chr21:48087183:48088183，chr22:38071168:38071189，chr7:26415938:26416562。

在一个或多个实施方案中，步骤(2)中标志物及其对应处理算法选自下述核酸分子或与其具有至少70％相同性的变体和算法的组合中的任意一种、多种或全部，所述变体中的甲基化位点未发生突变：chr1:119543868:119543879-mhl3、chr1:180202465:180202541-mhl3、chr1:203044677:203044823-mhl3、chr1:203045227:203045293-mhl3、chr10:124896784:124897154-mhl3、chr12:52745149:52745173-mhl3、chr13:92051635:92051674-mhl3、chr14:38071275:38071410-mhl3、chr14:38091345:38091744-mhl3、chr14:60977845:60977865-mhl3、chr16:85620063:85620331-mhl3、chr17:40937184:40937480-mhl3、chr17:46631705:46632377-mhl3、chr18:19747115:19747127-mhl3、chr18:70209347:70209475-mhl3、chr2:105474371:105474381-mhl3、chr2:203036052:203036236-mhl3、chr2:80530770:80530826-mhl3、chr2:95401381:95401629-mhl3、chr21:38077572:38077685-mhl3、chr4:111554966:111554989-mhl3、chr4:1400384:1400459-mhl3、chr4:40632502:40632519-mhl3、chr5:72594262:72594269-mhl3、chr6:168502296:168502347-mhl3、chr7:101241802:101241926-mhl3、chr7:157485865:157486031-mhl3、chr7:158798629:158798674-mhl3、chr7:27182264:27183525-mhl3、chr7:45002063:45002754-mhl3、chr7:50467662:50468309-mhl3、chr8:55366694:55366747-mhl3、chr9:139740669:139740676-mhl3、chr1:203044677:203044823-mhl、chr13:92051635:92051674-mhl、chr14:60977845:60977865-mhl、chr18:19747115:19747127-mhl、chr7:27182264:27183525-mhl、chr1:16862044:16862199-pdr、chr10:102822370:102822647-pdr、chr10:105253375:105253477-pdr、chr12:116946248:116946304-pdr、chr12:50355085:50355570-pdr、chr12:52652220:52652362-pdr、chr13:32605660:32605843-pdr、chr13:77459521:77459792-pdr、chr14:102172655:102172687-pdr、chr16:88599725:88600068-pdr、chr17:80009015:80009025-pdr、chr19:12831808:12832195-pdr、chr19:13213485:13213513-pdr、chr19:13213644:13213814-pdr、chr19:15344092:15344411-pdr、chr19:8674674:8674749-pdr、chr20:48902548:48902611-pdr、chr22:42710260:42710349-pdr、chr3:193987426:193987681-pdr、chr3:194208192:194208617-pdr、chr6:31696240:31696334-pdr、chr7:26415826:26415917-pdr、chr10:102822370:102822647-umhl、chr11:2000109:2000154-umhl、chr15:48470425:48470556-umhl、chr15:60883371:60883395-umhl、chr16:4715122:4715220-umhl、chr19:52104749:52104928-umhl、chr2:128158537:128158621-umhl、chr21:48087183:48088183-umhl、chr22:38071168:38071189-umhl、chr7:26415938:26416562-umhl。

在一个或多个实施方案中，步骤(3)中的模型是随机森林模型或支持向量机模型。

在一个或多个实施方案中，步骤(3)中的模型是随机森林模型，其由R语言的caret软件包中的train函数构建。

在一个或多个实施方案中，步骤(4)包括：当评分满足阈值时，则鉴定甲状腺肿瘤为良性或恶性。

在一个或多个实施方案中，步骤(1)中的检测包括但不限于：基于重亚硫酸盐转化的PCR(例如甲基化特异性PCR)、DNA测序(如亚硫酸氢盐测序、全基因组甲基化测序、简化甲基化测序)、甲基化敏感的限制性内切酶分析法、荧光定量法、甲基化敏感性高分辨率熔解曲线法、基于芯片的甲基化图谱分析、质谱(例如飞行质谱)。优选地，步骤(1)中的检测是测序。

在一个或多个实施方案中，所述方法在步骤(1)之前还包括：样品DNA的抽提、质检、和将DNA上未甲基化的胞嘧啶转化为不与鸟嘌呤结合的碱基。在一个或多个实施方案中，所述转化使用酶促方法进行，优选脱氨酶处理，或所述转化使用非酶促方法进行，优选用亚硫酸氢盐或重硫酸盐处理，更优选使用亚硫酸氢钙、亚硫酸氢钠、亚硫酸氢钾、亚硫酸氢铵、重硫酸钠、重硫酸钾和重硫酸铵处理。

本发明还提供一种装置，其特征在于，所述装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：

(1)获取样品中本文所述标志物中的一种或多种或全部的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，

(3)通过构建模型利用步骤(2)的经处理的甲基化水平获得评分，

(4)根据评分鉴别滤泡状甲状腺肿瘤良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

在一个或多个实施方案中，mhl、mhl3、umhl、pdr如本文第四方面所述。

在一个或多个实施方案中，所述装置用于鉴别滤泡状甲状腺肿瘤的良恶性或用于评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

本发明还提供一种计算机可读存储介质，其特征在于，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现下述步骤：

(1)获取样品中本文所述标志物中的一种或多种或全部的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，

(3)通过构建模型利用步骤(2)的经处理的各标志物的甲基化水平获得评分，

(4)根据评分鉴别滤泡状甲状腺肿瘤的良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

在一个或多个实施方案中，mhl、mhl3、umhl、pdr如本文第四方面所述。

本发明还提供一种用于鉴别滤泡状甲状腺肿瘤的良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能的系统，其特征在于，包括：

采集装置，用于获取样品中本文所述标志物中的一种或多种或全部的甲基化水平，

数据处理装置，使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，并通过构建模型利用经处理的各标志物的甲基化水平获得评分，

判定装置，用于根据评分鉴别甲状腺肿瘤的良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

在一个或多个实施方案中，采集装置包括样品处理装置和测序装置。

在一个或多个实施方案中，采集装置包括输入所述甲基化水平的装置。

在一个或多个实施方案中，所述模型是随机森林模型或支持向量机模型。

在一个或多个实施方案中，mhl、mhl3、umhl、pdr如本文第四方面所述。

本发明还提供一种构建鉴别滤泡状甲状腺肿瘤良恶性的模型的方法，包括：

(1)获取肿瘤样品和对照样品的基因组DNA中候选位点或片段的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述位点或片段的甲基化水平，

(3)筛选出经处理的甲基化水平在肿瘤样品和对照样品之间存在显著差异的位点或片段和任选的其对应算法，所述位点或片段即为甲基化标志物，

(4)根据甲基化标志物的经处理的甲基化水平，使用R语言的caret软件包中的train函数构建鉴别滤泡状甲状腺肿瘤良恶性的模型，其为随机森林模型。

在一个或多个实施方案中，mhl、mhl3、umhl、pdr如本文第四方面所述。

在一个或多个实施方案中，所述甲状腺肿瘤是滤泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

在一个或多个实施方案中，所述肿瘤样品包括良性(FTA)和/或恶性甲状腺肿瘤(FTC)组织样品；

在一个或多个实施方案中，所述对照样品来自(1)无所述肿瘤的对象的正常组织，或(2)与肿瘤样品来源相同的对象的正常组织。

在一个或多个实施方案中，获取候选位点或片段的甲基化水平包括(1)检测全基因组DNA的甲基化水平或包含所述候选位点或片段的基因组DNA的甲基化水平，并选取候选位点或片段的甲基化水平；或(2)检测所述候选位点或片段的甲基化水平。

在一个或多个实施方案中，步骤(3)包括：(3.1)使用R语言软件包对各样品的经处理的甲基化水平进行预处理得到横坐标为甲基化标志物，纵坐标为各样品的二维矩阵；和(3.2)根据经处理的甲基化水平在肿瘤样品和对照样品之间是否存在显著差异筛选位点或片段和任选的其对应算法。

在一个或多个实施方案中，显著差异为p值小于0.05。

在一个或多个实施方案中，所述甲基化标志物是本文第一方面所述的标志物。

在一个或多个实施方案中，所述肿瘤样品中来自良性甲状腺肿瘤患者的样品为30个以上，来自恶性甲状腺肿瘤患者的样品为30个以上。

在一个或多个实施方案中，构建随机森林模型时的参数为：method＝‘rf’，ntree＝500，trControl＝trainControl(method＝"repeatedcv",savePredictions＝T,classProbs＝T,number＝3,repeats＝10,allowParallel＝TRUE)。

在一个或多个实施方案中，甲基化标志物及其对应算法如下所示并且所述模型的鉴别阈值为0.5：chr1:119543868:119543879-mhl3、chr1:180202465:180202541-mhl3、chr1:203044677:203044823-mhl3、chr1:203045227:203045293-mhl3、chr10:124896784:124897154-mhl3、chr12:52745149:52745173-mhl3、chr13:92051635:92051674-mhl3、chr14:38071275:38071410-mhl3、chr14:38091345:38091744-mhl3、chr14:60977845:60977865-mhl3、chr16:85620063:85620331-mhl3、chr17:40937184:40937480-mhl3、chr17:46631705:46632377-mhl3、chr18:19747115:19747127-mhl3、chr18:70209347:70209475-mhl3、chr2:105474371:105474381-mhl3、chr2:203036052:203036236-mhl3、chr2:80530770:80530826-mhl3、chr2:95401381:95401629-mhl3、chr21:38077572:38077685-mhl3、chr4:111554966:111554989-mhl3、chr4:1400384:1400459-mhl3、chr4:40632502:40632519-mhl3、chr5:72594262:72594269-mhl3、chr6:168502296:168502347-mhl3、chr7:101241802:101241926-mhl3、chr7:157485865:157486031-mhl3、chr7:158798629:158798674-mhl3、chr7:27182264:27183525-mhl3、chr7:45002063:45002754-mhl3、chr7:50467662:50468309-mhl3、chr8:55366694:55366747-mhl3、chr9:139740669:139740676-mhl3、chr1:203044677:203044823-mhl、chr13:92051635:92051674-mhl、chr14:60977845:60977865-mhl、chr18:19747115:19747127-mhl、chr7:27182264:27183525-mhl、chr1:16862044:16862199-pdr、chr10:102822370:102822647-pdr、chr10:105253375:105253477-pdr、chr12:116946248:116946304-pdr、chr12:50355085:50355570-pdr、chr12:52652220:52652362-pdr、chr13:32605660:32605843-pdr、chr13:77459521:77459792-pdr、chr14:102172655:102172687-pdr、chr16:88599725:88600068-pdr、chr17:80009015:80009025-pdr、chr19:12831808:12832195-pdr、chr19:13213485:13213513-pdr、chr19:13213644:13213814-pdr、chr19:15344092:15344411-pdr、chr19:8674674:8674749-pdr、chr20:48902548:48902611-pdr、chr22:42710260:42710349-pdr、chr3:193987426:193987681-pdr、chr3:194208192:194208617-pdr、chr6:31696240:31696334-pdr、chr7:26415826:26415917-pdr、chr10:102822370:102822647-umhl、chr11:2000109:2000154-umhl、chr15:48470425:48470556-umhl、chr15:60883371:60883395-umhl、chr16:4715122:4715220-umhl、chr19:52104749:52104928-umhl、chr2:128158537:128158621-umhl、chr21:48087183:48088183-umhl、chr22:38071168:38071189-umhl和chr7:26415938:26416562-umhl。

在一个或多个实施方案中，模型的预测值大于或等于阈值则鉴定甲状腺肿瘤为恶性；模型的预测值小于阈值则鉴定甲状腺肿瘤为良性。

本发明还提供了一种鉴定甲状腺肿瘤良恶性的模型，该模型按照上述甲状腺肿瘤模型的构建方法的实施方案构建而得。在一个或多个实施方案中，所述甲状腺肿瘤滤是泡状甲状腺肿瘤或恶性潜能未定的甲状腺肿瘤。

本发明还提供一种筛选甲基化标志物的方法，包括以下步骤：

(1)获取肿瘤样品和对照样品的基因组DNA中候选位点或片段的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述位点或片段的甲基化水平，

(3)筛选出经处理的甲基化水平在肿瘤样品和对照样品之间存在显著差异的位点或片段和任选的其对应算法，所述位点或片段即为甲基化标志物。

在一个或多个实施方案中，所述肿瘤样品是滤泡状甲状腺肿瘤样品或恶性潜能未定的甲状腺肿瘤样品。

在一个或多个实施方案中，所述肿瘤样品包括良性(FTA)和/或恶性甲状腺肿瘤(FTC)组织样品；

在一个或多个实施方案中，所述对照样品来自(1)无所述肿瘤的对象的正常组织，或(2)与肿瘤样品来源相同的对象的正常组织。

在一个或多个实施方案中，获取候选位点或片段的甲基化水平包括(1)检测全基因组DNA的甲基化水平或包含所述候选位点或片段的基因组DNA的甲基化水平，并选取候选位点或片段的甲基化水平；或(2)检测所述候选位点或片段的甲基化水平。

在一个或多个实施方案中，步骤(3)包括：(3.1)使用R语言软件包对各样品的经处理的甲基化水平进行预处理得到横坐标为甲基化标志物，纵坐标为各样品的二维矩阵；和(3.2)根据经处理的甲基化水平在肿瘤样品和对照样品之间存在显著差异筛选位点或片段和任选的其对应算法。

在一个或多个实施方案中，显著差异为p值小于0.05。

在一个或多个实施方案中，所述甲基化标志物是本文第一方面所述的甲基化标志物。

在一个或多个实施方案中，mhl、mhl3、umhl、pdr如本文第四方面所述。

附图说明

图1是26例疑似甲状腺癌患者的甲基化标志物预测ROC曲线。灰色部分为95％置信区间。

图2是36例UMP样本的甲基化标志物恶性潜能评分及其基因变异。

具体实施方式

本发明探究了DNA甲基化水平与甲状腺肿瘤良恶性的关系。旨在通过无创的方法利用DNA甲基化标志物群作为甲状腺肿瘤良恶性鉴别标志物，提高甲状腺肿瘤无创诊断的精确性。本文中，甲状腺肿瘤包括临床上的甲状腺结节和肿瘤。

对于临床上根据病理鉴定无法明确进行FTA/FTC鉴定的患者，以及鉴定为UMP的甲状腺肿瘤患者，首先要收集他们的肿瘤组织样本(手术切除组织或者穿刺组织)，进行基因组DNA的抽提。可以是直接新鲜或冻存组织抽提，也可以先制作成石蜡包埋(FFPE)组织，然后再进行DNA抽提。对抽提后的DNA样本进行甲基化测序，可以是全基因组甲基化测序(WGBS)，或者是简并甲基化测序(RRBS)，还可以是根据标志物设计引物进行定制测序。

本文所述“甲基化标志物”或“标志物”既可以是核酸分子，也可以是核酸分子与其甲基化算法的组合。所述核酸分子可以是分离自动物的核酸分子，也可以是人工合成的具有动物基因组片段序列的核酸分子。本文的甲基化标志物中涉及的核酸分子可以不是以基因为单位，而是以CpG位点有连锁效应的区段为单位，即甲基化单倍型单元(MHB)。因此，本发明中，两个不同的标志物可以来自相同基因。而且，本发明中标志物也可位于基因间区。

发明人经过多次的尝试，最终确定了70个DNA区域和算法的甲基化特征，最能表现滤泡状甲状腺良性和恶性肿瘤之间的差异。这70个DNA区域的位置和算法见表1，其在样品中的甲基化水平变化见表2。

表1

表2

本文中，所示序列以及序列表所示序列视为正义链。当正义链5’-3’方向为CpG时，反义链上对应位置在5’-3’方向也为CpG。因此，在提及甲基化位点时，即包含正义链上甲基化位点的胞嘧啶，也包含反义链上该位点处对应碱基(鸟嘌呤)上游(5’端)相邻位置的胞嘧啶。

本文中，甲基化水平代表一个或多个位点处于甲基化状态的比例。一个区域(或一组位点)的甲基化水平是该区域中所有位点(或组中所有位点)的甲基水平的综合体现。因此，区域的甲基化水平上升或下降并不表示区域中所有甲基化位点的甲基化水平都上升或下降。本领域知晓将检测DNA甲基化的方法(例如简化甲基化测序)所得结果转化为甲基化水平的过程。示例性实施例利用软件Bismark(v0.17.0)获得CpG位点的甲基化水平。

本文中，检测样品DNA甲基化的方法本领域周知，例如基于重亚硫酸盐转化的PCR(例如甲基化特异性PCR(Methylation-specific PCR,MSP))、DNA测序(如亚硫酸氢盐测序(Bisulfite sequencing,BS)、全基因组甲基化测序(Whole-genome bisulfitesequencing,WGBS)、简化甲基化测序(Reduced Representation Bisulfite Sequencing,RRBS))、根据标志物设计引物进行PCR或测序、甲基化敏感的限制性内切酶分析法(Methylation-Sensitive Dependent Restriction Enzymes)、荧光定量法、甲基化敏感性高分辨率熔解曲线法(Methylation-sensitivity High-resolution Melting,MS-HRM)、基于芯片的甲基化图谱分析、质谱(例如飞行质谱)。在一个或多个实施方案中，检测包括检测基因或位点处的任一条链。

因此，本发明涉及检测DNA甲基化的试剂。本领域周知上述检测DNA甲基化的方法中所用的试剂。示例性地，检测DNA甲基化的试剂可包含以下的一种或多种：重亚硫酸盐及其衍生物、PCR缓冲液、聚合酶、dNTP、引物、探针、甲基化敏感或不敏感的限制性内切酶、酶切缓冲液、荧光染料、荧光淬灭剂、荧光报告剂、外切核酸酶、碱性磷酸酶、内标、对照物。在涉及DNA扩增的检测方法中，检测DNA甲基化的试剂包括引物。所述引物序列为甲基化特异的或非特异的。优选地，所述引物的序列包括非甲基化特异的封闭序列(Blocker)。封闭序列可以提高甲基化检测的特异性。检测DNA甲基化的试剂还可包括探针。通常，探针的序列的5’端标记荧光报告基团，3’端标记淬灭基团。示例性地，所述探针的序列包含MGB(Minorgroove binder)或者LNA(Locked nucleic acid)。MGB和LNA用于提高Tm(meltingtemperature)值，增加分析的特异性，提高探针设计的灵活性。

在示例性的实施方式中，本发明使用简化基因组甲基化测序(RRBS)检测DNA甲基化。简化基因组甲基化测序是利用限制性内切酶对基因组进行酶切，经Bisulfite处理，对基因组CpG区域进行测序的技术。步骤包括：1、利用限制性内切酶对基因组进行酶切；2、文库构建，包括末端修复、加A尾和接头3、片段长度分选；4、重亚硫酸盐转换；5、PCR扩增；6、测序。本文中，文库用Illumina Hiseq 2500测序仪进行双端测序，测序量为每个样本35～40M。示例性地，简化基因组甲基化测序所用试剂包括：血浆核酸纯化试剂盒、连接酶、重亚硫酸盐及其衍生物、dNTP、聚合酶、引物、无核酸酶水，任选的磁珠、乙酸钠、糖原。

本文中，样品来自哺乳动物，优选人。样品可来自任何器官(例如甲状腺)、组织(例如上皮组织、结缔组织、肌肉组织和神经组织)、细胞(例如甲状腺结节活检物)或者体液(例如血液、血浆、血清、组织液、尿液)。通常，只要所述样品包含基因组DNA或cfDNA(Circulating free DNA or Cell free DNA)即可。cfDNA称为循环游离DNA或者细胞游离DNA，是释放到血浆中的降解的DNA片段。示例性地，所述样品是甲状腺结节活检物，优选是细针穿刺活检物。或者，所述样品是血浆或cfDNA。

本发明还涉及鉴定甲状腺结节性质的试剂盒，包含本文所述，特别是本文第三方面所述的试剂。所述试剂盒还可包含本文所述，特别是第一方面所述的核酸分子作为内标或阳性对照。除了所述试剂和核酸分子之外，试剂盒还包含检测DNA甲基化所需的其他试剂。示例性地，检测DNA甲基化的其他试剂可包含以下的一种或多种：重亚硫酸盐及其衍生物、PCR缓冲液、聚合酶、dNTP、引物、探针、甲基化敏感或不敏感的限制性内切酶、酶切缓冲液、荧光染料、荧光淬灭剂、荧光报告剂、外切核酸酶、碱性磷酸酶、内标、对照物。

本文所述“引物”是指在核苷酸聚合作用起始时，引导合成的一种具有特定核苷酸序列的核酸分子。引物通常是人工合成的两段寡核苷酸序列，一个引物与靶区域一端的一条DNA模板链互补，另一个引物与靶区域另一端的另一条DNA模板链互补，其功能是作为核苷酸聚合作用的起始点。体外人工设计的引物被广泛用于聚合酶链反应(PCR)、qPCR、测序和探针合成等。通常，引物被设计为扩增的产物长度为50～150bp、60-140、70-130、80-120bp。

本文试剂中所含引物可以是基因组测序的引物，例如全基因组测序引物或针对基因组某一区域的测序引物，也可以是用于扩增特定区域的PCR引物或用于扩增区域中一个或多个甲基化位点的PCR引物。

例如，检测DNA区域所用的引物可以是全基因组测序引物，所述引物可以获得很多扩增产物，这些扩增产物可以包含所述区域或在拼接后包含所述区域。根据全基因组测序结果，在测序后获得该区域中的每个甲基化位点(CpG)的甲基化状态，从而获取整个区域的甲基化水平。

又如，检测DNA区域所用的引物可以是针对包含该区域的DNA进行测序的引物，所述引物可以获得较多扩增产物，这些扩增产物可以包含所述区域或在拼接后包含所述区域。在测序后获得该区域中的每个甲基化位点(CpG)的甲基化状态，从而获取整个区域的甲基化水平。

再如，检测DNA区域所用的引物可以是扩增该区域中一个或多个甲基化位点的PCR引物。这些引物的扩增产物可以包含区域中某一个或多个或全部甲基化位点，在对扩增产物中包含的甲基化位点进行检测后，获取整个区域的甲基化水平。

因此，检测某一区域所用的引物的扩增产物可以仅包含该区域中的一个或多个甲基化位点，也可以同时包含其他区域中的一个或多个甲基化位点。而且检测某一区域所需的引物可为一对或多对，例如1对、2对、3对、4对、5对、6对、7对、8对、9对、10对等，其中任一对引物的扩增产物至少包含该区域中的至少一个甲基化位点。

上述引物的描述同样适用于本文所述的其他DNA区域。本领域知晓全基因组测序引物或针对特定区域或区域中某一位点的PCR引物的设计方法。

本文术语“变体”或“突变体”是指与参照序列相比，通过一个或多个核苷酸的插入、缺失或取代使核酸序列发生变化同时保留其与其他核酸杂交能力的多核苷酸。本文任一实施方案所述的突变体包括与参照序列具有至少70％，优选至少80％，优选至少85％，优选至少90％，优选至少95％，优选至少97％的序列相同性并保留参照序列的生物学活性的核苷酸序列。可采用例如NCBI的BLASTn计算两条比对的序列之间的序列相同性。突变体还包括在参照序列的和核苷酸序列中具有一个或多个突变(插入、缺失或取代)、同时仍保留参照序列生物学活性的核苷酸序列。所述多个突变通常指1－10个以内，例如1－8个、1－5个或1－3个。取代可以是嘌呤核苷酸与嘧啶核苷酸之间的取代，也可以是嘌呤核苷酸之间或嘧啶核苷酸之间的取代。取代优选是保守性取代。例如，在本领域中，用性能相近或相似的核苷酸进行保守性取代时，通常不会改变多核苷酸的稳定性和功能。保守性取代例如嘌呤核苷酸之间的(A与G)的互换，嘧啶核苷酸之间的(T或U与C)的互换。因此，在本发明多核苷酸中用来自同一残基替换一个或几个位点，将不会在实质上影响其活性。此外，本发明的变体中所含有的本文所述的甲基化位点未发生突变。即本发明方法检测的是相应序列中的甲基化位点的甲基化情况，对于这些位点之外的碱基可以发生突变。

DNA或RNA的碱基之间可发生转化。本文所述“CT转化”是利用非酶促或酶促方法处理DNA，将未修饰的胞嘧啶碱基(C)转化为不与鸟嘌呤结合的碱基(例如尿嘧啶碱基(U))的过程。本文所述“AG转化”是利用非酶促或酶促方法处理DNA，将腺嘌呤(A)转化为鸟嘌呤(G)的过程。本领域周知进行转化的非酶促或酶促方法。示例性地，非酶促方法包括亚硫酸氢盐或重硫酸盐处理，例如亚硫酸氢钙、亚硫酸氢钠、亚硫酸氢钾、亚硫酸氢铵、重硫酸钠、重硫酸钾和重硫酸铵等。示例性地，酶促方法包括脱氨酶处理。经转化的DNA任选经纯化。适用于本文的DNA纯化方法本领域周知。

提及胞嘧啶时，“修饰”表示胞嘧啶碱基上的化学基团的引入或除去。在一个或多个实施方案中，修饰是指甲基化。本文所述“甲基化”或“DNA甲基化”是指在基因组DNA的CpG二核苷酸的胞嘧啶5'碳位共价结合一个甲基基团，成为5-甲基胞嘧啶(5mC)。

具体实施方案中，本发明的标志物通过如下方法筛选：(1)收集肿瘤样品和对照样品，其中滤泡状甲状腺肿瘤患者样品包括良性(FTA)和恶性甲状腺肿瘤(FTC)组织样品；(2)对样品中的基因组DNA进行甲基化检测，与147888个核酸片段(Guo et al.(2017))进行比较，找出所检测DNA样本中的所有测到的片段；(3)根据文中的mhl计算方法(包括mhl、mhl3、umhl和pdr)，以及的pdr计算方法(Landau et al.(2014))，算出每个核酸片段每种计算方法在每个样本中的甲基化评分，筛选出在良恶性甲状腺肿瘤中存在显著差异(矫正p值小于0.05)的片段及其对应算法。

在一个或多个实施方案中，与对照样品比较时，对象样品的甲基化水平升高或降低。对所测基因的甲基化水平进行数学分析，获得得分。对于检测的样品而言，当得分大于阈值，则判定结果为阳性，即为恶性结节，否则为阴性，即良性结节。本领域知晓常规数学分析的方法以及确定阈值的过程，示例性的方法是数学模型，例如支持向量机和随机森林模型。例如，对于差异NDA甲基化标志物，对两组样本构建支持向量机(SVM)和随机森林(Random forest，RF)模型，利用模型统计检测结果的准确率，敏感性和特异性以及预测值特征曲线(ROC)下面积(AUC)，统计测试集样本预测得分。

例如，本文可通过下述方法构建随机森林模型：(1)收集甲状腺肿瘤患者的肿瘤组织及其配对瘤旁样本，其中甲状腺肿瘤患者的组织样本包括良性(FTA)和恶性甲状腺肿瘤(FTC)组织样品，(2)检测样品基因组DNA中候选位点或片段的甲基化水平，使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述检测得到的甲基化水平，得到各样品的基因组DNA的甲基化信息，即经处理的甲基化水平，(3)使用R语言软件包对各样品的经处理的甲基化水平进行预处理得到横坐标为甲基化标志物，纵坐标为各样品的二维矩阵，(4)根据经处理的甲基化水平在肿瘤样品和对照样品之间是否存在显著差异筛选位点或片段和任选的其对应算法，(5)根据甲基化标志物的经处理的甲基化水平，使用R语言的caret软件包中的train函数构建滤泡状甲状腺癌良恶性鉴别模型，其为随机森林模型。在R语言环境中，载入caret模块，设置ctControl参数：ctrl＝trainControl(method＝"repeatedcv",savePredictions＝T,classProbs＝T,number＝3,repeats＝10,allowParallel＝TRUE)。

示例性的RF模型如下：

mod_rf<-train(imputed,used_pheno,method＝'rf',trControl＝ctrl)

类似地，示例性的SVM模型如下：mod_svm<-train(imputed,pheno,method＝'svmRadialSigma',trControl＝ctrl)

其中imputed为标志物二维矩阵；pheno为建模所用样本信息，FTC样本标注为p.ftc，FTA样本标注位n.fta。

上述两个模型在本文样本中的比较如下：

	准确性	灵敏性	特异性	AUC
					RF	0.9615385	0.9230769	1	0.9940828
SVM	0.6538462	1	0.3076923	0.9349112

两种模型的AUC均大于0.9，显示了较好的预测能力。其中RF的特异性更好，准确性更高。本文实施例中采用RF模型。

模型使用方法如下：使用该模型对滤泡状甲状腺肿瘤进行良恶性评估时，先对待评估的样本进行DNA抽提和甲基化测序，然后计算每个样本上述70个甲基化标志物的打分值，以形成一个横坐标为甲基化标志物，纵坐标为样本名的打分矩阵(Valmatrix)。如果矩阵中出现确失值(NA)，需先将其补齐。之后用R软件包中的randomForest模块进行预测，输入代码为：

predict(model,imputed,type＝"response")

输出结果为每个样本的良恶性预测结果，其中“n.fta”表示该样本可能为良性肿瘤，“p.ftc”表示该样本为恶性肿瘤的可能性更大。

如果想要查看每个样本为恶性肿瘤的可能性，可以用代码：

predict(model,imputed,type＝"prob")

此外，本发明还揭示了存储计算机程序的计算机可读存储介质，存储介质上所存储的计算机程序运行后执行本文所述的鉴别甲状腺良恶性的方法。结合本文中公开的实施方案描述的方法或算法的步骤可直接在硬件中、在由处理器执行的软件模块中、或在这两者的组合中体现。如果在软件中实现为计算机程序产品，则各功能可以作为一条或更多条指令或代码存储在计算机可读介质上或藉其进行传送。计算机可读介质包括计算机存储介质和通信介质两者，其包括促成计算机程序从一地向另一地转移的任何介质。存储介质可以是能被计算机访问的任何可用介质。作为示例而非限定，这样的计算机可读介质可包括RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁存储设备、或能被用来携带或存储指令或数据结构形式的合意程序代码且能被计算机访问的任何其它介质。示例性存储介质耦合到处理器以使得该处理器能从/向该存储介质读取和写入信息。在替换方案中，存储介质可以被整合到处理器。处理器和存储介质可驻留在ASIC中。ASIC可驻留在用户终端中。在替换方案中，处理器和存储介质可作为分立组件驻留在用户终端中。

实施例

下面结合附图和具体实施例对本发明作进一步详细的说明。下列实施例中，未注明具体条件的实验方法，通常按常规条件中所述的方法进行。

研究对象

良性甲状腺肿瘤患者33人，恶性甲状腺肿瘤患者33人用于筛选标志物和构建模型。26例疑似滤泡状甲状腺癌患者用于甲状腺良恶性鉴别。36例临床病理诊断为恶性潜能未定的甲状腺肿瘤(UMP)患者用于评估甲状腺肿瘤的恶性潜能。

对于这些对象，手术或穿刺活检取其肿瘤组织。其中大部分组织用作病理检测，少部分组织用Qiagen组织DNA提取试剂盒抽提全基因组DNA进行甲基化分析。

简化甲基化测序(RRBS)

组织DNA抽提

1.将组织样本剪碎，置于1.5mL的离心管中，并加入裂解缓冲液和proteinase K，充分振荡混匀后，置于56℃孵育至少一小时，直至组织接近完全裂解。

2.加入RNase A，振荡混匀后室温孵育2分钟。加入终止缓冲液，脉冲式振荡混匀15秒，随后70℃孵育10分钟。

3.加入无水乙醇，振荡混匀后短暂离心，将全部混合物转管至纯化回收柱，离心1分钟，更换新的废液收集管。

4.加入洗涤缓冲液到纯化回收柱中，离心1分钟，弃废液；再次加入洗涤缓冲液到纯化回收柱中，全速离心3分钟，弃废液；全速离心1分钟空甩后丢弃废液收集管。

5.将纯化回收柱放入新的1.5毫升离心管中，加入50-100μL洗脱液，室温孵育5分钟后离心1分钟来洗脱DNA。

6.质检：DNA不少于0.5ug

DNA甲基化文库构建

1.阳性对照：与样本DNA等量的NA12878DNA标准样

阴性对照：与样本DNA等体积的纯水

2.DNA消化：使用限制性内切酶消化基因组DNA和对照样本

3.末端修复与加dA：使用末端修复和加dA专用酶混合液和缓冲液，修复酶切产物DNA或ctDNA的末端，产生平端，并在平端的3’末端上加上一个dA，形成3’的突出，加热酶反应以灭活用于上述反应的酶。

4.连接接头序列：在灭活后加入DNA连接酶和接头序列，催化加A产物与接头序列之间的连接反应。

5.连接产物纯化：使用AMPure XP DNA纯化磁珠，纯化回收连接产物。

6.使用亚硫氰酸钠处理连接产物：在人类基因组DNA中甲基化仅发生在胞嘧啶碱基上。使用亚硫氰酸钠处理ctDNA后，未被甲基化的胞嘧啶碱基被转化为尿嘧啶(dU)，在下一步PCR扩增步骤中被进一步扩增为dT，并在测序结果中被检测出来；而甲基化的胞嘧啶不受亚硫氰酸钠处理影响，仍然为dC，因此在测序结果中仍显示为dC。因此亚硫氰酸钠处理后的连接产物将未甲基化的dC转化为dT，保留甲基化的dC不变，从而为在测序结果中区别两者奠定基础。如前文所述，接头序列中dC均已被甲基化，不受亚硫氰酸钠处理影响

7.PCR扩增和产物纯化：PCR方法扩增DNA，构建文库；使用AMPure XP DNA纯化磁珠，纯化回收扩增产物

8.质检：使用LabChip检测文库的丰度与片段分布

a)文库主要信号集中在170-400bp区间，其中峰值在250-350bp左右

b)接头二聚体产物形成的尖峰信号应显著低于文库主峰信号

c)NA12878 DNA对应的文库有较强的信号，而阴性对照中应无明显信号

9.甲基化文库测序：在使用Illumina测序平台对甲基化文库进行测序

实施例1，标志物筛选和模型构建

甲基化检测：对来自良性甲状腺肿瘤(FTC)患者33人，恶性甲状腺肿瘤(FTA)患者33人的样品利用RRBS技术进行甲基化测序和建库。

标志物筛选：将甲基化检测结果与Guo et al.(2017)中定义的147888个核酸片段进行比较，找出所检测DNA样本中的所有测到的片段。根据mhl计算方法(mhl、mhl3、umhl和pdr)和pdr计算方法(Landau et al.(2014))，算出每个核酸片段每种计算方法在每个样本中的甲基化评分，筛选出在良恶性甲状腺肿瘤中存在显著差异(矫正p值小于0.05)的片段及其对应算法，得到70组标志物和算法的组合。

数据处理：每个样本的测序结果进行质控分析后，将标志物区段的序列挑选出来，并根据每个标志物对应的甲基化值计算方式进行计算，如果有区段内某CpG位点测序深度低于10×，则该区段的甲基化值为NA。形成样本-标志物的数值矩阵，部分示例性数据见表3：

表3 10个甲状腺肿瘤样本的标志物甲基化值矩阵

构建模型

模型构建使用的是本文前面的33例FTC和33例FTA的甲基化水平矩阵(markers)，66例样本的分类信息存储为第一列为样本名，第二列为分类信息(FTC名为p.ftc，FTA名为n.fta)的矩阵(pheno)。

打开R程序包，导入66例样本的分类信息矩阵和相应的甲基化水平矩阵

pheno＝read.delim(“分类信息矩阵的存储路径”,sep＝"\t",as.is＝T,header＝T,check.names＝FALSE)

markers＝ead.delim(“甲基化水平矩阵的存储路径”,sep＝"\t",as.is＝T,header＝T,check.names＝FALSE)

将甲基化水平矩阵转置为行名为样本名称,列名为甲基化标志物名称的矩阵

imput＝t(markers)

用临近值补充法补上NA值

library(DMwR)

imputed＝knnImputation(imput)

设置建模参数

library(caret,quietly＝T)

ctrl<-trainControl(method＝"repeatedcv",savePredictions＝T,classProbs＝T,number＝3,repeats＝10,allowParallel＝TRUE)

构建随机森林模型

mod_rf<-train(imputed,pheno,method＝'rf',trControl＝ctrl)

存储模型

saveRDS(mod_rf$finalModel,file.path(“存储路径”,"RFmodelFTC.rds"))

实施例2，甲状腺良恶性鉴别

本实施例利用实施例1构建的模型对26例疑似滤泡状甲状腺癌样品进行甲状腺良恶性鉴别。过程如下：

根据实施例1的方法，构建类似表3的26个样品的样本-标志物的数值矩阵。

打开R程序包，导入26例待评估样本-标志物的数值矩阵

valdata＝read.delim(“样本-标志物矩阵的存储路径”,sep＝"\t",as.is＝T,row.names＝1,header＝T,check.names＝F)

将矩阵转置为行名为样本名称,列名为甲基化标志物名称的矩阵

imput＝t(valdata)

并用临近值补充法补上NA值

library(DMwR)

imputed＝knnImputation(imput)

再导入已建立的恶性甲状腺癌潜能的评估模型(RFmodelFTC)

model＝readRDS("RFmodelFTC存储路径")

之后开始进行模型评估，得出每一例样本是更倾向于恶性或者良性

library(randomForest)

class＝predict(model,imputed,type＝"response")

并计算出每例样本是恶性肿瘤的概率

probs＝predict(model,imputed,type＝"prob")

将最终结果与病理检测结果进行比较，发现26例样本，仅一例(14号样品)的甲基化标志物模型预测与病理不一致(见表3)，预测灵敏度92.3％，特异性100％，准确性为96.2％，ROC曲线下面积AUC为0.994(图1)。该结果展现了甲基化标志物对甲基化肿瘤预测结果与临床病理结果的高度一致性。

表4 26例甲状腺肿瘤患者的良恶性鉴别

实施例3，恶性潜能未定的甲状腺肿瘤的恶性潜能的评估

本实施例利用实施例1构建的模型对36例临床病理诊断为恶性潜能未定的甲状腺肿瘤(UMP)患者进行甲状腺肿瘤的恶性潜能的评估。

样本处理：用Qiagen组织DNA提取试剂盒抽提全基因组DNA后，一部分DNA样本进行由18个恶性肿瘤基因形成的基因panel突变检测，包括TERT,EIF1AX,HRAS,NRAS,KRAS,BRAF,TP53,PIK3CA,PTEN，GNAS,TSHR,CTNNB1,AKT1和ETV6的所有全外显子区域，部分内含子区域，以及部分基因的启动子区域的突变，以及RET,PPARG,ALK和NTRK1的融合。另外的DNA样本进行基于RRBS技术的甲基化测序。

根据实施例1的方法，构建类似表3的UMP样本的样本-标志物的数值矩阵。

打开R程序包，导入UMP样本的甲基化标志物数值矩阵

valdata＝read.delim(“样本-标志物矩阵的存储路径”,sep＝"\t",as.is＝T,row.names＝1,header＝T,check.names＝F)

将矩阵转置为行名为样本名称,列名为甲基化标志物名称的矩阵

imput＝t(valdata)

并用临近值补充法补上NA值

library(DMwR)

imputed＝knnImputation(imput)

再导入已建立的恶性甲状腺癌潜能的评估模型(RFmodelFTC)

model＝readRDS("RFmodelFTC存储路径")

进行恶性潜能打分：

probs＝predict(model,imputed,type＝"prob")

根据甲基化打分结果，可以将所有UMP样本分为三类(三个风险级别)：1——低风险，甲基化分值为0-0.4；2——中风险，甲基化分值为0.4-0.6；3——高风险，甲基化分值为0.6-1。将甲基化打分结果与基因突变检测结果进行比较(见图2)，7例低风险样本均未检测到恶性突变(0％)，14例中风险样本有5例检测到恶性突变(35.7％)，15例高风险样本有11例检测到突变(73.3％)。高风险组的突变样本比例显著高于低风险组(p＝0.004)，且显著高于中低风险组(p＝0.006)。

序列表

<110> 上海鹍远生物技术有限公司

<120> 甲状腺肿瘤良恶性鉴别模型及其应用

<130> 19A658F1

<160> 66

<170> SIPOSequenceListing 1.0

<210> 1

<211> 12

<212> DNA

<213> Homo sapiens

<400> 1

cgggcctgcg cg 12

<210> 2

<211> 77

<212> DNA

<213> Homo sapiens

<400> 2

cgcgccttgg tcgcgcggcc cccggtgcgg gcctcgtcgc tccgacagcc ctacagctgt 60

ccctacggcc cgcaccg 77

<210> 3

<211> 147

<212> DNA

<213> Homo sapiens

<400> 3

cgggcctctg gtgtccccca tggtgcaggg ggatgacaag gtgtttcgcc gcgcgccctc 60

ctggaggaag cgcttccggc cgcgggagca ccacggtcgc ggcggcatgc tcagcgcttc 120

cgcggagacc ctcccggcgg gcttccg 147

<210> 4

<211> 67

<212> DNA

<213> Homo sapiens

<400> 4

cgcgggcacc ttgctgctgg ccctcggccc cacgcacctg cccttggctg ccggcctggc 60

cctgccg 67

<210> 5

<211> 371

<212> DNA

<213> Homo sapiens

<400> 5

cggccactcc gggcgcagaa gactggaaga agggcgctga aagtccagag aagaagccgg 60

cgtgccgcaa gaagaagacg cgcacagtct tctcgcgcag ccaggtcttc cagctcgagt 120

ccaccttcga catgaagcgc tatctgagca gctcggagcg agccggcctg gccgcgtccc 180

tgcacctcac cgagacgcag gtcaagatct ggttccagaa ccgccgcaac aagtggaagc 240

ggcagctggc ggcggagctg gaggcggcca acctgagcca tgccgcggcg cagcgcatcg 300

tgcgggtgcc catcctctac cacgagaact cggcggccga gggcgcggcg gctgcagccg 360

cgggggcccc g 371

<210> 6

<211> 25

<212> DNA

<213> Homo sapiens

<400> 6

cgtcactgat ggcagcccgg ccacg 25

<210> 7

<211> 40

<212> DNA

<213> Homo sapiens

<400> 7

cgcagcgcag ctctggggac ccctaaccaa gcctgcgacg 40

<210> 8

<211> 136

<212> DNA

<213> Homo sapiens

<400> 8

cggcttagag cgcgctgagt gtccgttggg gcccctgctc ttgggggcgc ctggggctct 60

gcgcgcccgc agggcaggtt tgcgcaccga gtgtcacaga agtgatcttc ctcggtaccg 120

gcctcttagc tggccg 136

<210> 9

<211> 400

<212> DNA

<213> Homo sapiens

<400> 9

cggcaggaaa ctaaaaagga ttttctccga ttcaagcaga agttggcctc caagccggct 60

gtagatgaaa gcccagtcca cagcctccat gcccccggcc cggcccgccc cgcgcgcgtt 120

tcctgcgccg cagcccggac gtcgagaaga tacccttcgc ttaaaggccc tgcgatgtcc 180

gcggccgccc tcctgcagga ggtgctcgga ggcgcgccgc gtccctcggg cctgggtgag 240

gcggcggctc caggcaagac ccggtcgttt cgcccccggg acttttactt gcggagctcc 300

gcgttcctac ggcaccaggc cctgaaaaag cccccagtca tcgcctcggg gttcggcacg 360

gccagacccg tggtcctgct gcctccgccc gagccacccg 400

<210> 10

<211> 21

<212> DNA

<213> Homo sapiens

<400> 10

cgggcggagg gcgacggcac g 21

<210> 11

<211> 269

<212> DNA

<213> Homo sapiens

<400> 11

cgaaggctgt ttcatgagcg tcatcatcgg tcggccctgg tatgtttgtt cttccagggg 60

ctcccaggat ggatccagtc cctcccggga gccgggccac accagctgtc ccactgccag 120

gccgcatggc caagagtcgc agctgctcct acaacgtgct cacgggcccg acgccggccc 180

ccgcgccccc ccgactgccg tgctccatgg cttcccatgt ggaaagtatc acccctccca 240

cctctgagag cagccgtgtg agacccccg 269

<210> 12

<211> 297

<212> DNA

<213> Homo sapiens

<400> 12

cgcacggata agaacgagag gtgggggtga aagggcagag cgtgggtaga atagggccgc 60

gggccgcggc gggctcggct cacccacgcg tcaccctcag gttcaccatc caggacctcc 120

tggcccacgc cttcttccgc gaggagcgcg gtgtgcacgt ggaactagcg gaggaggacg 180

acggcgagaa gccgggcctc aagctctggc tgcgcatgga ggacgcgcgg cgcggggggc 240

gcccacggga caaccaggcc atcgagttcc tgttccagct gggccgggac gcggccg 297

<210> 13

<211> 673

<212> DNA

<213> Homo sapiens

<400> 13

cgatgatgat gatagtaata ataataaaac aataatttaa ctagatgcct ggggccgaaa 60

ttcctgccgc tccccttgca ataaacccca aaccatcgcg ggacggaggc caggcgagtg 120

tggaaagcga aggagccaga gacgcgatag gaaagaatgg agactccgcc ggtggtgcag 180

acagggctgg gaaggtttgt gcacccgggt agtccctggc tgctgctgcc aggctgcctc 240

agccggtcgc tgctgccgcg gcgactggcg acaagctacc agccacctac gatggcccaa 300

ggaggcgaag aagagcaagc gatcaggaca ccaaaacggt tatcggaggc aggttcccag 360

cacaccagcc ggccgagggc cgagccccgc gggcggcagc aagttttggg agctggaggt 420

aaccgaatta aaaggcgcct tagaaactcc gcttcgggac tttgctcagc agggctccgg 480

gttggagggc gccgaggcct ggcggacggg acagtgggaa gagagaaagg tgctaagggg 540

acccaagatc tgggatccag aacaagaggg ggtggggaac aactctacca agccaaacag 600

atctatttcc cttgcctcca tgttggaaaa attcaggttc catatggctc ctcgggaggg 660

agggagggag acg 673

<210> 14

<211> 13

<212> DNA

<213> Homo sapiens

<400> 14

cgttgtcgaa tcg 13

<210> 15

<211> 129

<212> DNA

<213> Homo sapiens

<400> 15

cgggcatcat cagccgcagc ccgagtggcc cccggccggg cgcctgcatc gggactggtg 60

ggaggcggcg cgcgggggtg gaggccggcg ccggcgcgag cggcgcggaa gggcgcgaag 120

gaacgcgcg 129

<210> 16

<211> 11

<212> DNA

<213> Homo sapiens

<400> 16

cggccgagtc g 11

<210> 17

<211> 185

<212> DNA

<213> Homo sapiens

<400> 17

cggcatctgc gcctccgcgg cggcaggggc cccgcagcag gcgcgctggg agcagcttct 60

ctcggggcct gccgcgttta tggcttcaat agcgccgctc cgggtccccg cggcgtcgga 120

gaggcttctg cctagcgcct ctctgctgcg cgaggtcgtc ccggccttac aaaaggagct 180

cctcg 185

<210> 18

<211> 57

<212> DNA

<213> Homo sapiens

<400> 18

cggtgaggtt gagcgcctcg cagtacagca gccgcccctc gcaccggcac agctgcg 57

<210> 19

<211> 249

<212> DNA

<213> Homo sapiens

<400> 19

cgggatgcgc agcgggacca gcgctgcccg cgcctctcac cgcactgcat ccgcctcccg 60

ccagccagga agccactgcg gcctgaggct tcccgccacc accgcgcgcg ccttcctccg 120

gggacatggg gagctggctg aaggcgtaaa ggagcgagca gaagcctcag gccagacaca 180

gcgccaccac gcgcggtagc gccgcatggc cccagccgcg ttcctcggtc tccgtctccg 240

ccgcgcccg 249

<210> 20

<211> 114

<212> DNA

<213> Homo sapiens

<400> 20

cggtcctcac cctcctctcc cgccacgcac atatccttct tgacttcgaa gtggtttgca 60

atccgaaagt gagaccttga gtcctcagat ggccggcaac gcgccgaggt cacg 114

<210> 21

<211> 24

<212> DNA

<213> Homo sapiens

<400> 21

cgggcgagcc aaatcttttt cacg 24

<210> 22

<211> 76

<212> DNA

<213> Homo sapiens

<400> 22

cgccgcctta ggcgcccgcg gcgcccggtc cgcgcattta tggcagcccc gccggaggcg 60

cccacgcgcc cacacg 76

<210> 23

<211> 18

<212> DNA

<213> Homo sapiens

<400> 23

cgcgcgcagg caggagcg 18

<210> 24

<211> 8

<212> DNA

<213> Homo sapiens

<400> 24

cgcgggcg 8

<210> 25

<211> 52

<212> DNA

<213> Homo sapiens

<400> 25

cgaccctccc cggcaggagc gcgggccgtg cagggtctcg ggggtcgaac cg 52

<210> 26

<211> 125

<212> DNA

<213> Homo sapiens

<400> 26

cgcagggccg ggcgggccca ctttcaggtg tatgatctca cctctcggaa actggcaagc 60

cgccggaggg cctcgtcaag gctttgatga ccctgtgctg aggacgcaaa tcctaggcct 120

ggccg 125

<210> 27

<211> 167

<212> DNA

<213> Homo sapiens

<400> 27

cggaggaggg gacctgctgc cctcagcctg gctggtaacc ggcctctcca tagcaacggc 60

cagcgcgcgc gtctgtgtgt gcgcgcgtgt ctgatgtgtg tgtgcccgtg gtgttcccgg 120

gactccccgc gggggctggg aggggatcgc agaaccctag ggtggcg 167

<210> 28

<211> 46

<212> DNA

<213> Homo sapiens

<400> 28

cgcgctccgg tgcctgcggc cgaggaagag gacggcctcc cgcgcg 46

<210> 29

<211> 1262

<212> DNA

<213> Homo sapiens

<400> 29

cgctgtttcc ccagcgtagc cctcctcata aattatccgc cgtgacaagc ccgattcacg 60

gctgctacag ccatcctcta cctctctgcg ccttgctcgg ctggcctgac ccgggagcgc 120

gtcccaaggc gtggggttcc agaggggttt tttgcttcct cccccttcca acgtctaaac 180

tgtcccagag aacgcccatt tcccccacta tttgtgagcg cagggtgctc gcaaagaaga 240

ggaggaagga ggaaggcagg ggagggagaa cggcaaggag agctccgcag ggctgggaga 300

aatgagacca agagagactg ggagagggcg gcagagaaga gaggggggac cgagagccgc 360

gtccccgcgg tcgcgtggat ttagaaaaag gctggcttta ccatgactta tgtgcagctt 420

gcgcatccag gggtagatct ggggttgggc gggcggcgcc gggctcggct cgctctgcgc 480

actcgcctgc tcgctgctgg caggggcgtc ctcctcggct ccggacgccg tgccaacccc 540

ctctctgctg ctgatgtggg tgctgccggc gtcggccgag gcgccgctgg agttgcttag 600

ggagtttttc ccgccgtggt ggctgtcgct gccgggcgag ggggccacgg cggagcaggg 660

cagcggatcg ggctgaggag agtgcgtgga cgtggccggc tggctgtacc tgggctcggc 720

gggcgccgcg ctggcgctgg cagcgtagct gcgggcgcgc tctccggagc caaagtggcc 780

ggagcccgag cggccgacgc tgagatccat gccattgtag ccgtagccgt acctgccgga 840

gtgcatgctc gccgagtccc tgaattgctc gctcacggaa ctatgatctc cataattatg 900

caactggtag tccgggccat ttggatagcg accgcaaaat gagtttacaa aataagagct 960

catttgtttt ttgatatgtg tgcttgattt gtggctcgcg gtcgtttgtg cgtctatagc 1020

acccttgcac aatttatgat gaattatgga aatgactggg acatgtactt ggttccctcc 1080

tacgtaggca cccaaatatg gggtacgact tcgaatcacg tgcttttgtt gtccagtcgt 1140

aaatcctgcc tgatgacctc tagaggtaaa ctcgtgcact aataggggag ttgggtggag 1200

gcgagggggg tggcgcgcgc gccccgggcg cgtgcccgcc gccagttgcc gccgttcagc 1260

cg 1262

<210> 30

<211> 692

<212> DNA

<213> Homo sapiens

<400> 30

cgtggacggg gcaggctggg gcacaaaggg cactgcttcc tctggatagc gctaactccg 60

agagctcggg ctcctatccc tagccggggc gcctggcctc ctgcagtcgg gtccccagcg 120

ggtgagccgg gcgcaggtgg acctgaccga ggaggcggcc ggagcccctg cgggtgaggc 180

gggagcggcg agcaggaggg ctgactcaga aggtttattt gcagcgctgg cggggcggac 240

aattggcggc ctcggggtgc ggcgggtgcg ggcgctcagc ggctgggcca gagcagggtg 300

aaggagccgc gagcgcagcg gaaatcgggc tctggctgct ccggcctggg ctccacggag 360

atgaggcgcc ggaccccgcg atggcttagt gggatgcagt cggagacgcg aacctccagg 420

gtgcggccct ccctgcggca ggaggagggg tcagggcggc cacgcggccg gggcttcgtg 480

cccgctaccg cccagcctgc actcacccct ggcagtgtgc ggccagcacg cccaccagct 540

ccccaacgcg gttgtccaat ggcggccggg agcggtgcag gcacaccacg aggtcgtcct 600

ggccgcgggc gtgcagcacc accagctggt ctcctccgct ggtcacgctc agccccgtca 660

cctgagcgga gcgcggggtc agagtgcagc cg 692

<210> 31

<211> 648

<212> DNA

<213> Homo sapiens

<400> 31

cgagaaggag aacgaaatga tgaagtccca cgtgatggac caagccatca acaacgccat 60

caactacctg ggggccgagt ccctgcgccc gctggtgcag acgcccccgg gcggttccga 120

ggtggtcccg gtcatcagcc cgatgtacca gctgcacaag ccgctcgcgg agggcacccc 180

gcgctccaac cactcggccc aggacagcgc cgtggagaac ctgctgctgc tctccaaggc 240

caagttggtg ccctcggagc gcgaggcgtc cccgagcaac agctgccaag actccacgga 300

caccgagagc aacaacgagg agcagcgcag cggtctcatc tacctgacca accacatcgc 360

cccgcacgcg cgcaacgggc tgtcgctcaa ggaggagcac cgcgcctacg acctgctgcg 420

cgccgcctcc gagaactcgc aggacgcgct ccgcgtggtc agcaccagcg gggagcagat 480

gaaggtgtac aagtgcgaac actgccgggt gctcttcctg gatcacgtca tgtacaccat 540

ccacatgggc tgccacggct tccgtgatcc ttttgagtgc aacatgtgcg gctaccacag 600

ccaggaccgg tacgagttct cgtcgcacat aacgcgaggg gagcaccg 648

<210> 32

<211> 54

<212> DNA

<213> Homo sapiens

<400> 32

cgccgcctag acggcgccgg gacgcgaggt gcgaatcggc ggcgtccggc tgcg 54

<210> 33

<211> 8

<212> DNA

<213> Homo sapiens

<400> 33

cgcggccg 8

<210> 34

<211> 1262

<212> DNA

<213> Homo sapiens

<400> 34

cgctgtttcc ccagcgtagc cctcctcata aattatccgc cgtgacaagc ccgattcacg 60

gctgctacag ccatcctcta cctctctgcg ccttgctcgg ctggcctgac ccgggagcgc 120

gtcccaaggc gtggggttcc agaggggttt tttgcttcct cccccttcca acgtctaaac 180

tgtcccagag aacgcccatt tcccccacta tttgtgagcg cagggtgctc gcaaagaaga 240

ggaggaagga ggaaggcagg ggagggagaa cggcaaggag agctccgcag ggctgggaga 300

aatgagacca agagagactg ggagagggcg gcagagaaga gaggggggac cgagagccgc 360

gtccccgcgg tcgcgtggat ttagaaaaag gctggcttta ccatgactta tgtgcagctt 420

gcgcatccag gggtagatct ggggttgggc gggcggcgcc gggctcggct cgctctgcgc 480

actcgcctgc tcgctgctgg caggggcgtc ctcctcggct ccggacgccg tgccaacccc 540

ctctctgctg ctgatgtggg tgctgccggc gtcggccgag gcgccgctgg agttgcttag 600

ggagtttttc ccgccgtggt ggctgtcgct gccgggcgag ggggccacgg cggagcaggg 660

cagcggatcg ggctgaggag agtgcgtgga cgtggccggc tggctgtacc tgggctcggc 720

gggcgccgcg ctggcgctgg cagcgtagct gcgggcgcgc tctccggagc caaagtggcc 780

ggagcccgag cggccgacgc tgagatccat gccattgtag ccgtagccgt acctgccgga 840

gtgcatgctc gccgagtccc tgaattgctc gctcacggaa ctatgatctc cataattatg 900

caactggtag tccgggccat ttggatagcg accgcaaaat gagtttacaa aataagagct 960

catttgtttt ttgatatgtg tgcttgattt gtggctcgcg gtcgtttgtg cgtctatagc 1020

acccttgcac aatttatgat gaattatgga aatgactggg acatgtactt ggttccctcc 1080

tacgtaggca cccaaatatg gggtacgact tcgaatcacg tgcttttgtt gtccagtcgt 1140

aaatcctgcc tgatgacctc tagaggtaaa ctcgtgcact aataggggag ttgggtggag 1200

gcgagggggg tggcgcgcgc gccccgggcg cgtgcccgcc gccagttgcc gccgttcagc 1260

cg 1262

<210> 35

<211> 156

<212> DNA

<213> Homo sapiens

<400> 35

cgaaggcctg gagcgagttg cagcgacccg gccgcagctc accactggac tagagatgcg 60

cctttgcgag gtggcagcaa gtgaccagtc ggtcgtgcgt cgccaggtcc ggagccgcgc 120

accaggttgc caggaggagg cgggagcgcg gaggcg 156

<210> 36

<211> 278

<212> DNA

<213> Homo sapiens

<400> 36

cgcagctgcc ttggcgctgc ctgtgctcct gctactgctg gtggtgctga cgccgccccc 60

gaccggcgca aggccatccc caggcccaga ttacctgcgg cgcggctgga tgcggctgct 120

agcggagggc gagggctgcg ctccctgccg gccagaagag tgcgccgcgc cgcggggctg 180

cctggcgggc agggtgcgcg acgcgtgcgg ctgctgctgg gaatgcgcca acctcgaggg 240

ccagctctgc gacctggacc ccagtgctca cttctacg 278

<210> 37

<211> 103

<212> DNA

<213> Homo sapiens

<400> 37

cgacccgcgc cgccgagagc ttgcggctac gggaacgcgg cgcgccccgg ggccacatct 60

ggcctagcgg gcgcgcgcgt gtcactaaga cgctcattca ccg 103

<210> 38

<211> 57

<212> DNA

<213> Homo sapiens

<400> 38

cgtcgccgcc ggagagggtc ccagcgcgcc caggcctccc gcaggcggga cgcggcg 57

<210> 39

<211> 486

<212> DNA

<213> Homo sapiens

<400> 39

cggccgcccg cccccgccgc cggcggtgag ggaggtgagc ggcgccgacc tgcgggacga 60

gcatcactcc gacccagccg ggggtgaggc gggtcaggat gctccggtcg caggaggaaa 120

aggaggagct ggaccaaaag cccgaagaga agaaaagggg aaggccgcgc acggagcgcg 180

gtaaaggccg gcggagctag acgccccgag gtcggagtga agcgccggga ccgagccccg 240

tctcccaggg agtccggggc gcacggcacc gaggagagcg cgggagccaa cctgggcgca 300

tcatgcgcag ggcccgggac gctgggccgg tctacaccgc cgcctgggtc acgtggcccg 360

gacgggccgg cggctgcccc ggccgggggg cgggggtcgc gccggggttg cgctggacga 420

cggagagcgg cgggcccgca gcggcctgga gcctcccaac ccgcgcgccg cgctggcccc 480

cgagcg 486

<210> 40

<211> 143

<212> DNA

<213> Homo sapiens

<400> 40

cgccgcggta gggggcggcg gtgatgcagc agcggccggg ccggggcccg caggccgaga 60

tgcagctgaa ggcgcggacc ccgcaggggg cgctgaggcg ggagtagagg caggacatgg 120

tgcggtggcg gcaggggcgg gcg 143

<210> 41

<211> 184

<212> DNA

<213> Homo sapiens

<400> 41

cgcggccgcc gcctccgagg gctgcaggga gatcagcgtc cagcaaataa gaagcaagtc 60

ctggacccgg aggaggagga gcggccgagc atctctctct gctccgccgt gtcctttaga 120

tgagcactcc cggccggagc cggaggtgga tccgcagagc tgcctctggg cgcctgaccc 180

cgcg 184

<210> 42

<211> 272

<212> DNA

<213> Homo sapiens

<400> 42

cgttcagggt gtccccaaac acctccttgg ccagcgacgt ctttccgcaa acggtgatgc 60

gcgccactcg ccggaacttg gcgtccgcct gcgcgtcccg cccgatggtg taggagccgc 120

ggtagccgat ggtgatgtag cccgagcgcc ggctgccgtc cagcgactgg gacggcgtga 180

gcagcgggcc cgccgcgccc ccggacggac tgcggctagc cagctccagc gtgggcgacg 240

gcgccccggc agaggcgccc tcctgctgtt cg 272

<210> 43

<211> 33

<212> DNA

<213> Homo sapiens

<400> 43

cgcacttcct gggcgaggcc gctgcaggca gcg 33

<210> 44

<211> 344

<212> DNA

<213> Homo sapiens

<400> 44

cgagcccccg gcggccccca ggagcatcaa ggtggaggcg gtggaggagc cggaggcggc 60

ccccatcctg ggccccggag agcctgggcc ccaggccccg tcgcggacgc cgtcgccgcg 120

cagccccgcc ccggccaggg tcaaggccga gctgtccagc cccacgccgg gctccagccc 180

ggtgcccggc gagctgggcc tggccggggc cctgttcctt ccgcagtacg tgttcgggcc 240

cgacgcggcg ccccccgcct cggagatcct ggccaagatg tccgagctgg tgcacagccg 300

gctgcagcag ggcgcgggcg cgggcgccgg cggcgcgcag accg 344

<210> 45

<211> 11

<212> DNA

<213> Homo sapiens

<400> 45

cgcaggccgc g 11

<210> 46

<211> 388

<212> DNA

<213> Homo sapiens

<400> 46

cggaggtagg ggccggggtc agcgctgggt ctctggggct ccgtgtgaga cccaggtgga 60

gcccctgagg ccgcggtggc cgcatgacga cgggaacgcc ctcggcggac aggcggaggc 120

ctccgatcca cgcccgccca agtgcggggt ccccgccagc tgcgccatcc tcggccgcgc 180

agtcgcgggc tcgggagcgc gggagggggg atgtggaaac gggccacagg cggcggcggc 240

ggggcccggc ggatcctcat gggacccagc gaagagccct cgtcccaggg tggcccatgc 300

cgctgacccc ggctccgagc ccgaaggctt ccacctccct cgcagcggct gggcgagggc 360

ccagccgcct ccacatccag gggccccg 388

<210> 47

<211> 29

<212> DNA

<213> Homo sapiens

<400> 47

cgttgggccg cgcctcctcg tgggcctcg 29

<210> 48

<211> 171

<212> DNA

<213> Homo sapiens

<400> 48

cggaggatgc tggcccagtg cagataaggg ccaggctgcc tggccgcgcc tgataaggag 60

cctggctgac cccggaggaa accagcgggc ggggacctag cagcgcgggc ctcccggccc 120

cgcccccagc cggcgctggc aggactcggt ccgccacccc tcagggcttc g 171

<210> 49

<211> 320

<212> DNA

<213> Homo sapiens

<400> 49

cgcgggcggt attagcggct ggaggaggct tcgggaggcc cggccgacgg ccgccgcctg 60

gtgctaccca cccaggggcg cgcgaccctc ccttcggtct ggctccaaag acctagcagc 120

actgacttca cccagctgtg gttccaacgg cgggtccagc ggcctcggcc cggcgccgtc 180

ctcctgctgg cccaacaggc ccgccagccc gcccctgtac gtctgtgatt ggacggcggc 240

ggccactgat gttcaagcga caggtcctgg cccgggagcc aatctgcagg tgttgaggcc 300

caggctccga gagcgggccg 320

<210> 50

<211> 76

<212> DNA

<213> Homo sapiens

<400> 50

cgacacacac cccgcgcaac cagagaccgg cgctcgggca ctgtagaccc acaaaagacg 60

caggacgcct ccaccg 76

<210> 51

<211> 64

<212> DNA

<213> Homo sapiens

<400> 51

cggtgaactg gagagaacag gctctcctca ggggacggtg gcctctgagc cagcaggcgg 60

cccg 64

<210> 52

<211> 90

<212> DNA

<213> Homo sapiens

<400> 52

cggatgtgtg ctctggggca gatgcctaat gcgaatcgct tctgtcatca ggacacctgg 60

gggcttgaca cggaccagct gccgctgccg 90

<210> 53

<211> 256

<212> DNA

<213> Homo sapiens

<400> 53

cgccctcagt cccccggcgg gccctgccac cacccctcag cgttgccgga tgtgcggctg 60

aggtccaagg tttggagagt ttccgggact ggccgggtca ccctgcggct cggccgtgtc 120

tggcagcccc acccttcctg gcagccaggc cctctgaaat tccgccagca ccgggggacg 180

aagaaggtct tgctaggatg ggggagggcg ggaagacggg gagcacgagc cggctttgcc 240

ctcccatcca gccacg 256

<210> 54

<211> 426

<212> DNA

<213> Homo sapiens

<400> 54

cgcccccgcc aaaacaaaca gacaaaaaac tgtcattctc atatgaccac ggactctgcg 60

ccccggccgc aggcggagga cagggaggag cgcacacgag aaagctccca cgcgcccgcg 120

cctcgcctcc gacgggaagg cgcctcttcc gaccgtcctg gatgcaaatt aaatacttcc 180

ctccgcagaa gacttatctc ggggtagggc cgcagccgga tttccaaatc gcgggtttga 240

ttttcgctgc gtcctctccg cgacgcgcgc agggagcgcg cccggccacc tgcagcctct 300

ggctgcggcc atctttccgc ctggggccac ctggcctgcg ccctttgtgc ccggcccgct 360

ccgtgactgt tcttcttgcc gctacttaca gaactcagcc gccttagcgt cttagagaag 420

aggacg 426

<210> 55

<211> 95

<212> DNA

<213> Homo sapiens

<400> 55

cgccacgatc tcagctcctc ggtgattggt ccatttggag aggcctacga aaaactcccg 60

gcctgagtcc gggaggccgc ggaggtttga gggcg 95

<210> 56

<211> 92

<212> DNA

<213> Homo sapiens

<400> 56

cgtcccgcac acgcgcgctc ccgccgcgtc ccgcccgggc cgctccgctg gcggcggcga 60

ctgggctgca ggtagcgagc cagccggcct cg 92

<210> 57

<211> 278

<212> DNA

<213> Homo sapiens

<400> 57

cgcagctgcc ttggcgctgc ctgtgctcct gctactgctg gtggtgctga cgccgccccc 60

gaccggcgca aggccatccc caggcccaga ttacctgcgg cgcggctgga tgcggctgct 120

agcggagggc gagggctgcg ctccctgccg gccagaagag tgcgccgcgc cgcggggctg 180

cctggcgggc agggtgcgcg acgcgtgcgg ctgctgctgg gaatgcgcca acctcgaggg 240

ccagctctgc gacctggacc ccagtgctca cttctacg 278

<210> 58

<211> 46

<212> DNA

<213> Homo sapiens

<400> 58

cggcgggcgg caggtggagg aacaaaggca gggcctggcg atgacg 46

<210> 59

<211> 132

<212> DNA

<213> Homo sapiens

<400> 59

cgtccgccat cccgccgccg ctgcctccgc ctcggccgcc tgagctgagg ggctccgagc 60

gcgacaatgg cggctgccgg ggaagcggta acggatcccg atttgaaaag gccccgagcg 120

ctcagtaacc cg 132

<210> 60

<211> 25

<212> DNA

<213> Homo sapiens

<400> 60

cgcattaggg aagccgcgtc cgccg 25

<210> 61

<211> 99

<212> DNA

<213> Homo sapiens

<400> 61

cgcccacgtg ctcttggctg cctttgaaaa ggtaagtgcc atctggccat cactttcccg 60

gggacctggg agctgggcag ggggtggcct tttgagacg 99

<210> 62

<211> 180

<212> DNA

<213> Homo sapiens

<400> 62

cgcgcggccc ccgcgccccg gcgcctgccc ccgctgtggg gcgctggcca gccgcgcgct 60

ctgccaggcc tgctccctcc tggacggcct gaaccgcggt cggccccgcc tggccatcgg 120

caagggccgc cgggggctgg acgaggaggc gacgccgggg acgcccgggg atccggcccg 180

<210> 63

<211> 85

<212> DNA

<213> Homo sapiens

<400> 63

cgggtggtaa aagtccagcc cggcggagcg cccctgcctc cttccacgcc cctggcctgc 60

cctgcgcagc tacgcccagc ccgcg 85

<210> 64

<211> 1001

<212> DNA

<213> Homo sapiens

<400> 64

cgcctgctgg gctctgtgcg cccgggaagg tgcggccacc ctcacgcgga aggcggccag 60

cggatcccgg tgcgcgcagc tcccagcgct ggggttccag cgccccgcct cttcctatag 120

caaccagcgg gacctgccgt cccccggggc accccgaggg gtctgcgccc gcttctttcc 180

gaaacgggaa ggcgctgggg gctcggcagc cagagggacg ggttcaggga gcgtccggtg 240

agcctaagac gcgcctttgc cggggttgcc gggtgtctgc ctctcactta ggtattagga 300

accgtggcac aaatctgtag gttttcctct gggggtgggc ggaggctcca aaccggacgg 360

ttttctcctg gaggactgtg ttcagacaga tactggtttc cttatccgca ggtgtgcgcg 420

gcgctcgcaa gtggtcagca taacgccggg cgaattcgga aagcccgtgc gtccgtggac 480

gacccacttg gaaggagttg ggagaagtcc ttgttcccac gcgcggacgc ttccctccgt 540

gtgtccttcg agccacaaaa agcccagacc ctaacccgct cctttctccc gccgcgtcca 600

tgcagaactc cgccgttcct gggaggggaa gcccgcgagg cgtcgggaga ggcacgtcct 660

ccgtgagcaa agagctcctc cgagcgcgcg gcggggacgc tgggccgaca ggggaccgcg 720

ggggcagggc ggagaggacc cgccctcgag tcggcccagc cctaacactc aggaccgcct 780

ccagccggag gtctgcgccc ttctgaggac cctgcctggg ggagcttatt gcggttcttt 840

tgcaaatacc cgctgcgctt ggacggagga agcgcccacg cgtcgacccc ggaaacgaag 900

gcctccctga tgggaacgca tgcgtccagg agcctttatt tactcttaat tctgcccgat 960

gcttgtacgt gtgtgaaatg cttcagatgc ttttgggagc g 1001

<210> 65

<211> 22

<212> DNA

<213> Homo sapiens

<400> 65

cgcgcacgca gcaggagcca cg 22

<210> 66

<211> 625

<212> DNA

<213> Homo sapiens

<400> 66

cgccctggta acccatggca actgcgcggc ggaaagaaaa aaaagaaccg tggtggggcg 60

gggggcggcg ggcggtgaat gggaacaacg accgcaaaga gacagcgatt taccgcgcgg 120

gcctgcggct cccggccgca cctgctgcag acggcccgcg aggccccttc ccgcgctacg 180

gtgacgattc tggctgctgg ggaaaagagc aataagcaga aagccctcgt gcagagggga 240

tccgggcggt gcacttggtg tgggaggctg ccttaaacgc cgatgaccct cccggcctcc 300

gcgtgtcccg gcttccaggg ccccgccagg ctggagagcg cgcgtggaga gggccttgcc 360

gcctgggggt tggttgaggg gcggccgcgc cggggccggg gctggactcc aggctttgtt 420

ctgcagacgc tcgcgcccgg ccggaggagg ggctacaccg tgctgccccg gccccatggg 480

gcccggcccc cgagggtccc cgcgagcgga cgcggtgggg ccgggcaaac tctacgtgcc 540

ctaaattttg tcatctgcac acgcatcgca cacattagtt aacccctttc ccttctaggc 600

ccccgagaac ctaacctgcc cggcg 625

Claims (5)

1.一种试剂盒，包含检测样品中甲基化标志物的甲基化的试剂和任选的所述甲基化标志物，所述甲基化标志物选自分离自动物的以下任意一个、多个或全部：TBX15的基因间区、LHX4的第1个内含子、MYOG的第34个外显子、MYOG的3’UTR、HMX3的启动子、KRT85的基因间区、GPC5的启动子、FOXA1的基因间区、FOXA1的的第2个外显子、SIX6的启动子、GSE1的基因间区、WNK4的启动子、HOXB3的启动子、GATA6-AS1的启动子、CBLN2的启动子、POU3F3的启动子、KIAA2012的第15个内含子、LRRTM1的启动子、FAM95A的基因间区、SIM2的第1个内含子、PITX2的第2个内含子、CRIPAK的基因间区、RBM47的启动子、TMEM174的基因间区、FRMD1的基因间区、MYL10的基因间区、MIR153-2的第12个内含子、LINC00689的启动子、HOXA5的启动子、MYO1G的3’UTR、FIGNL1的第8个外显子、SOX17的基因间区、AJM1的启动子、MIR3675的第1外显子、KAZALD1的启动子、NEURL1的启动子、LINC00173的基因间区、AQP5的启动子、KRT86的第1个内含子、FRY的启动子、KCTD12的启动子、LINC00239的基因间区、ZC3H18的的第10个外显子、RFNG的启动子、TNPO2的启动子、LYL1的启动子、EPHX3的启动子、ADAMTS10的启动子、LINC01270的基因间区、OGFRP1的基因间区、LINC00887的基因间区、LINC00884的启动子、DDAH2的启动子、SNX10的基因间区、MRPL23-AS1的基因间区、MYEF2的启动子、RORA的启动子、NUDT16L1的第7个外显子、LINC01530的基因间区、PROC的基因间区、PRMT2的基因间区、LGALS1的启动子，或包含其的基因组片段，或与其具有70％序列相同性的变体，所述片段的长度为5-1500bp，

所述甲基化标志物通过包括以下步骤的方法获得：

(1)获取肿瘤样品和对照样品的基因组DNA中候选位点或片段的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述位点或片段的甲基化水平，

(3)筛选出经处理的甲基化水平在肿瘤样品和对照样品之间存在显著差异的位点或片段和任选的其对应算法，所述位点或片段即为甲基化标志物。

2.如权利要求1所述的试剂盒，其特征在于，所述甲基化标记物是核酸分子，所述核酸分子具有(1)动物基因组的以下序列，或与其具有至少70％相同性的变体，所述变体中的甲基化位点未发生突变：chr1:119543868:119543879，chr1:180202465:180202541，chr1:203044677:203044823，chr1:203045227:203045293，chr10:124896784:124897154，chr12:52745149:52745173，chr13:92051635:92051674，chr14:38071275:38071410，chr14:38091345:38091744，chr14:60977845:60977865，chr16:85620063:85620331，chr17:40937184:40937480，chr17:46631705:46632377，chr18:19747115:19747127，chr18:70209347:70209475，chr2:105474371:105474381，chr2:203036052:203036236，chr2:80530770:80530826，chr2:95401381:95401629，chr21:38077572:38077685，chr4:111554966:111554989，chr4:1400384:1400459，chr4:40632502:40632519，chr5:72594262:72594269，chr6:168502296:168502347，chr7:101241802:101241926，chr7:157485865:157486031，chr7:158798629:158798674，chr7:27182264:27183525，chr7:45002063:45002754，chr7:50467662:50468309，chr8:55366694:55366747，chr9:139740669:139740676，chr1:16862044:16862199，chr10:102822370:102822647，chr10:105253375:105253477，chr12:116946248:116946304，chr12:50355085:50355570，chr12:52652220:52652362，chr13:32605660:32605843，chr13:77459521:77459792，chr14:102172655:102172687，chr16:88599725:88600068，chr17:80009015:80009025，chr19:12831808:12832195，chr19:13213485:13213513，chr19:13213644:13213814，chr19:15344092:15344411，chr19:8674674:8674749，chr20:48902548:48902611，chr22:42710260:42710349，chr3:193987426:193987681，chr3:194208192:194208617，chr6:31696240:31696334，chr7:26415826:26415917，chr11:2000109:2000154，chr15:48470425:48470556，chr15:60883371:60883395，chr16:4715122:4715220，chr19:52104749:52104928，chr2:128158537:128158621，chr21:48087183:48088183，chr22:38071168:38071189，和chr7:26415938:26416562(2)(1)的互补序列。

3.一种装置，其特征在于，所述装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：

(1)获取样品中权利要求1或2中所述的甲基化标志物的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，

(3)通过构建模型利用步骤(2)的经处理的甲基化水平获得评分，

(4)根据评分鉴别滤泡状甲状腺肿瘤良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

4.一种计算机可读存储介质，其特征在于，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，实现下述步骤：

(1)获取样品中权利要求1或2中所述的甲基化标志物的甲基化水平，

(2)使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，

(3)通过构建模型利用步骤(2)的经处理的各标志物的甲基化水平获得评分，

(4)根据评分鉴别滤泡状甲状腺肿瘤的良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

5.一种用于鉴别甲状腺肿瘤或评估恶性潜能未定的甲状腺肿瘤的恶性潜能的系统，其特征在于，所述系统包括：

采集装置，用于获取样品中权利要求1或2中所述的甲基化标志物的甲基化水平，

数据处理装置，使用选自mhl、mhl3、umhl和pdr中的至少一种算法处理所述各标志物的甲基化水平，并通过构建模型利用经处理的各标志物的甲基化水平获得评分，

判定装置，用于根据评分鉴别甲状腺肿瘤的良恶性或评估恶性潜能未定的甲状腺肿瘤的恶性潜能。

Images