一种通过大样本转录组测序鉴定结直肠癌中特异分子的方法与流程

1.本发明涉及医学领域，具体涉及一种通过大样本转录组测序鉴定结直肠癌中特异分子的方法。


背景技术：

2.结直肠癌(colorectal cancer,crc)是一种始于结肠(大肠)或直肠的癌症。这两个器官都位于消化系统的下部。美国癌症协会(acs)估计，大约每23名男性和每25名女性中就有1 人会在他们的一生中罹患结肠直肠癌。国际癌症研究机构2018年全球评估显示，每年约有180 万结直肠癌的新发病例，并且90万人死于结直肠癌，使其成为第三大最常见的恶性肿瘤和第二大癌症死亡原因。某些固定的因素，如既往有结肠息肉病史、肠道疾病史、结直肠癌家族史、某些遗传综合症(如家族性腺瘤性息肉病(fap))、东欧犹太人或非洲人后裔、年龄等会增加患结肠直肠癌的风险。此外，工业化和经济增长导致西方饮食模式、久坐不动的生活方式和日益增加的肥胖，都是结直肠癌的危险因素。近年来由于结肠镜筛查的应用，结直肠癌患者的发病率和死亡率有所下降，但约25％的患者仍然存在iv期疾病，而iv期患者的 5年生存率不超过10％。
3.不同的结直肠癌分期以及患者的整体健康状况直接决定着结直肠癌的治疗方式的选择。在结直肠癌早期阶段，切除根治术是首选的治疗方式，之后使用卡培他滨(xeloda)、氟尿嘧啶、铂(eloxatin)、伊立替康(camptosar)等化疗药物或者放射对癌细胞进行后续控制。此外，靶向药物和免疫治疗也是结直肠癌治疗的有效方式，一般用于治疗转移性或晚期结肠直肠癌，根据患者靶基因的表达情况，选择合适的药物。但是目前仍有部分患者对靶向或者免疫治疗药物不敏感。
4.因此，研究结直肠癌相关的靶向分子和细胞活动，对于结直肠癌的治疗具有显著的意义。针对鉴定出来的靶向基因设计药物或者治疗方案可以为患者提供更多的治疗方案，延长患者的生存时间。


技术实现要素：

5.本发明的目的在于，基于大样本的结直肠癌患者的组织转录组数据，利用生物信息学的方法鉴定出在结直肠癌组织中显著变化的基因与信号通路，这有助于更加利用显著变化的基因对结直肠癌的个性化治疗和预后判断提供依据。为实现上述目的，本发明采用如下技术方案：
6.一种通过大样本转录组测序鉴定结直肠癌中特异分子的方法，包括如下步骤：
7.s1.数据采集
8.选取并大样本量采集诊断为结直肠癌的患者，在术中收集患者的癌组织与癌旁组织；
9.s2.转录组测序
10.使用trizol试剂提取组织样本的总rna，使用安捷伦2100rnanano6000检测试剂盒
检测总rna的完整性和浓度；总rna样本合格后，选择磁珠进行富集和纯化，目的片段经琼脂糖凝胶电泳提取，pcr扩增，然后采用pe150测序策略的illumina平台对合格的文库进行测序；
11.s3.数据分析
12.s3
‑
1.运用r语言中的“dseq2”包来筛选差异基因，筛选标准：log2foldofchange&padj＜ 0.05；然后用主成分分析方法根据差异基因表达对患者进行聚类；
13.s3
‑
2.利用gsea软件对癌组织和癌旁组织的显著变化的信号通路进行鉴定；
14.s3
‑
3.根据差异基因的表达情况，进一步筛选差异基因；
15.s3
‑
4.用验证集对筛选出来的差异基因的表达情况进行验证。
16.本发明所述方法筛选出来的特征因子有6个，分别为：ighg1、ca2、col1a1、 rna5
‑
8sn2、mir3648
‑
1与aqp8。
17.作为改进，s1数据采集中所述患者的选取标准为：排除其他重大疾病与精神疾病患者，排除样本临床信息资料不全的患者。
18.作为改进，鉴定出的特异因子应用于靶向药物的设计。
19.本发明的优点在于：
20.本发明通过大样本临床数据研究，可以鉴定出结直肠癌中显著变化的分子，对研究治疗结直肠癌的治疗方法具有深远意义。同时，针对显著变化的分子进行靶向药物设计，可以缓解患者的耐药情况，延长患者的生存时间。
附图说明
21.图1为差异基因的火山图、热图和主成分分析图；
22.图2为mtorc1、oxidative_phosphorylation、dna_repair信号通路gsea 展示图；
23.图3为ighg1、ca2、col1a1、rna5
‑
8sn2、mir3648
‑
1与aqp8表达条形图；
24.图4为验证集中ighg1、ca2、col1a1、rna5
‑
8sn2、mir3648
‑
1、aqp8的roc分析；
25.图5为aqp8、ca2和col1a1表达量与预后。
具体实施方式
26.下面通过具体实施例对本发明进行详细和具体的介绍，以使更好的理解本发明，但是下述实施例并不限定本发明的保护范围。
27.实施例1
28.本实施例公开了一种通过大样本转录组测序鉴定结直肠癌中特异分子的方法，包括如下步骤：
29.s1.数据采集
30.选取并采集大样本量诊断为结直肠癌的患者，在术中收集患者的癌组织与癌旁组织；数据采集中所述患者的选取标准为：排除其他重大疾病与精神疾病患者，排除样本临床信息资料不全的患者。
31.s2.转录组测序
32.使用trizol试剂提取组织样本的总rna，使用安捷伦2100rnanano6000检测试剂盒检测总rna的完整性和浓度；总rna样本合格后，选择磁珠进行富集和纯化，目的片段经琼脂
糖凝胶电泳提取，pcr扩增，然后采用pe150测序策略的illumina平台对合格的文库进行测序。
33.s3.数据分析
34.s3
‑
1.运用r语言中的“dseq2”包来筛选差异基因，筛选标准：log2foldofchange&padj＜ 0.05；然后用主成分分析方法根据差异基因表达对患者进行聚类；
35.s3
‑
2.利用gsea软件对癌组织和癌旁组织的显著变化的信号通路进行鉴定；
36.s3
‑
3.根据差异基因的表达情况，进一步筛选差异基因；
37.s3
‑
4.用验证集对筛选出来的差异基因的表达情况进行验证。
38.发明人按照实施例1所用方法对所在医院的结直肠癌患者进行分析鉴定：
39.s1.数据采集
40.选取2018年2月至2019年10月，在医院病理诊断为结直肠癌的144名患者，排除有其他重大疾病与精神疾病的患者，以及样本临床信息不全的患者，在术中收集患者的癌组织与癌旁组织，其中2018年2月至2019年3月的队列作为训练队列，2019年4月至2019年10 月的患者作为验证队列。
41.s2.转录组测序
42.使用trizol试剂提取组织样本的总rna；使用安捷伦2100rnanano6000检测试剂盒检测总rna的完整性和浓度；总rna样本合格后，选择磁珠进行富集和纯化。目的片段经琼脂糖凝胶电泳提取，pcr扩增。然后采用pe150测序策略的illumina平台对合格的文库进行测序。
43.s3.数据分析
44.s3
‑
1：
45.运用r语言中的“dseq2”包来筛选差异基因，根据筛选标准：log2foldofchange&padj＜0.05 共筛选出839个显著变化的基因，其中648个基因显著上调，191个基因显著下调。差异基因的热图和火山图如图1a
‑
b所示。进一步运用主成分分析方法根据差异基因表达对患者进行聚类，结果表明差异基因能显著区分癌组织与癌旁组织(图1c)。
46.s3
‑
2：
47.利用gsea软件对癌组织和癌旁组织的显著变化的信号通路进行鉴定，通过比对 hallmark数据库的50条经典信号通路，其中41信号通路在癌组织中显著激活，9条信号通路在癌旁组织中显著激活。图2展示了其中三条经典的信号通路(mtorc1、 oxidative_phosphorylation、dna_repair)的激活图。
48.s3
‑
3:
49.根据差异基因的表达情况进一步筛选出6个差异基因，ighg1、ca2、col1a1、 rna5
‑
8sn2、mir3648
‑
1与aqp8。如图3所示，在癌组织中ca2和aqp8显著降低， ighg1、col1a1、rna5
‑
8sn2与mir3648
‑
1显著增加。
50.s3
‑
4：
51.筛选出来的6个差异基因进一步用验证集对其表达情况进行验证，利用roc曲线，结果表明上述6个分子可独立筛选出癌组织与癌旁组织，并且具有良好的曲线下面积：ighg1 (0.752：0.662
‑
0.842)、ca2(0.832：0.753
‑
0.912)、col1a1(0.769：0.680
‑
0.857)、 rna5
‑
8sn2(0.902：0.849
‑
0.956)、mir3648
‑
1(0.865：0.800
‑
0.930)、aqp8(0.869： 0.799
‑
0.938)。
52.进一步利用gepia在线数据库对ca2、col1a1和aqp8表达量与预后的关系进行关联分析，结果显示高表达的aqp8患者有显著增加的总生存期(p＜0.05)，高表达的ca2患者表现出更好的总生存期(p＝0.076)，低表达的col1a1表现出更好的总生存期(p＝0.067) (图5)。
53.本次分析鉴定结果鉴定出6个在结直肠癌中显著变化的特异分子：ighg1、ca2、col1a1、 rna5
‑
8sn2、mir3648
‑
1与aqp8。对研究治疗结直肠癌的治疗方法具有深远意义。同时，针对显著变化的分子进行靶向药物设计，可以缓解患者的耐药情况，延长患者的生存时间。
54.本发明的具体实施例进行了详细描述，但其只是作为范例，本发明并不等同于以上描述的具体实施例。对于本领域技术人员而言，任何对本发明进行的等同修改和替代也都在本发明的范畴之中。因此，不脱离本发明的精神和范围下所做的均等变换和修改，都应涵盖在本发明的范围内。