CN107723352A - 一种循环肿瘤dna肝癌驱动基因高通量检测方法 - Google Patents

一种循环肿瘤dna肝癌驱动基因高通量检测方法 Download PDF

Info

Publication number
CN107723352A
CN107723352A CN201610661868.5A CN201610661868A CN107723352A CN 107723352 A CN107723352 A CN 107723352A CN 201610661868 A CN201610661868 A CN 201610661868A CN 107723352 A CN107723352 A CN 107723352A
Authority
CN
China
Prior art keywords
gene
sequence
monokaryon
liver cancer
former times
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610661868.5A
Other languages
English (en)
Inventor
张道允
巩子英
叶建伟
王伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiaxing Permitted Medical Inspection Co Ltd
Original Assignee
Jiaxing Permitted Medical Inspection Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiaxing Permitted Medical Inspection Co Ltd filed Critical Jiaxing Permitted Medical Inspection Co Ltd
Priority to CN201610661868.5A priority Critical patent/CN107723352A/zh
Publication of CN107723352A publication Critical patent/CN107723352A/zh
Pending legal-status Critical Current

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Microbiology (AREA)
  • Immunology (AREA)
  • Biotechnology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种用于循环肿瘤DNA肝癌驱动基因高通量检测方法,本方法选取目前肝癌常见的突变基因,结合第二代高通量测序技术,更全面的检测肝癌的基因突变情况。由于可以采用患者的血浆样本,因此,本方法可以极大的降低患者检测的创伤;其次、对75个肝癌基因进行全外显子检测,比基于real‑time PCR和一代Sanger测序技术平台的检测方法更全面;由于只对75个基因进行检测,能以更高的测序深度对样本进行检测,因此,检测的精度可以达到0.1%。从而为患者的更多、更优的个体化精准治疗方案提供参考。

Description

一种循环肿瘤DNA肝癌驱动基因高通量检测方法
技术领域
本发明涉及一种基因突变检测方法,尤其涉及循环肿瘤DNA肝癌驱动基因高通量检测方法。
背景技术
肝细胞癌为第五大最常见的癌症,全球癌症相关死亡率排第三。中国是肝癌大国,2012年全国肿瘤统计显示中国男性肝癌发病率为293,318人,女性肝癌发病率为101,452人,男女共有394770人发病,占到全球新发肝癌患者的52.7%,名副其实的全球第一。预计在2016年,中国男性肝癌发病率为326,698,女性肝癌发病率为111,017,两者相加437,715人,相当于中国一个县的人口,防癌与治疗形式十分严峻。如何筛选有效的肿瘤分级分子标记物并根据分子标记物进行针对性的治疗是目前迫切关注的问题。
在亚洲,包括中国、韩国等,由于地方性乙型肝炎病毒感染等原因肝癌发生率很高。目前已知的肝癌发生分子机制中,高频突变主要在维护端粒长度的信号通路上,主要是TERT启动区突变(59%)和HBV插入TERT启动区突变引起(20%);其次是Wnt/β-catenin信号通路,主要是CTNNB1(10-32%)和AXIN1(15%)突变引起;其它还有细胞周期控制途径、表观遗传修饰通路、AKT/mTOR和Ras/raf/MAP激酶信号通路、氧化胁迫通路、JAK/STAT通路等;相应发生突变的基因有TP53、CDKN2A、ATM、IRF2、MLL、FGF19、PIK3CA和JAK1等。最近,有研究表明,在59%的肝癌患者以及25%的肝硬化癌前病变患者中TERT启动子区域都发生了突变。 这些结果预示着TERT可能是一个很好的预后指标。端粒酶逆转录酶(TERT)编码催化端粒酶逆转录活性的一个亚基,目的是保持端粒的长度。TERT活性的增加是肿瘤发生的一个标志,TERT的转录调控是其癌症发生激活的主要原因。TERT启动子区域是其表达调控最重要的区域,这里结合了许多的调控因子,能够直接激活或者抑制TERT的表达。在膀胱癌的研究中,TERT启动子区域的突变使得患者预后不好,生存期缩短,复发几率远大于TERT启动子区域野生型患者。TERT启动子区域的突变与BRAF突变经常共同存在于乳头状甲状腺癌中,这部分患者预后与TERT启动子区域野生型患者也有很大的差异。TERT启动子区域野生型的肝癌患者的预后要好于突变型患者。由于TERT启动子区域突变在肝癌中频发并与肿瘤的发生密切相关。此外, CTNNB1、AXIN1、TP53、CDKN2A、ATM、IRF2、FGF19、PIK3CA和JAK1等基因的突变与肝癌的发生都有显著的相关性。这些基因突变在中国人群中的分布和预后的相关性方面都亟待研究。因此,研究这些基因突变与肝癌预后及复发的关系是一个重要的课题。
当前,癌症基因检测主要以组织样本为主。基于肿瘤组织样本的检测,存在两个问题:一、肿瘤组织大都来源于活检技术即穿刺和活检钳活检,这样会给病人带来极大的痛苦,是有创检测;二、对于某些晚期的病人,不能穿刺或者手术,病人样本就很难获取,因此就很难根据检测的结果对病人进行靶向用药。而且,采用靶向治疗的患者,都已经处于中晚期,身体状况较差。因此,尽量减少检测对身体的创伤具有非常重要意义。另一方面,随着精准治疗的不断发展,越来越多靶向治疗可供选择;同时也存在越来越多的靶向位点需要检测。因此,会产生一个问题,更多的治疗选择需要从患者身上获得更多的样品。基于上述发展过程中产生的矛盾,临床上迫切需要一种能够替代目前市场上需要较多样本的基因检测方法。近年来研究发现肿瘤患者外周血存在肿瘤循环DNA(ctDNA);而且越来越多的临床研究表明是非常好的组织样本的代替品。本专利采用血浆作为检测样本,有如下优势:1、提供肿瘤实时信息---当前肿瘤的突变状态;2、非侵入性液体活检---减少患者的成本和风险;3、没有肿瘤组织可用---仍然可以方便进行检测;4、消除组织样本偏差---评估全身肿瘤的整体状态。
因此,个体化精准治疗,能提高晚期肺癌患者的生存质量和生存时间。但是需要对靶标位点进行精准检测。而利用血浆样本,结合第二代高通量测序技术,可以更全面的检测肺癌的基因突变信息,从而获得更好的治疗方案,提高预后效果。
发明内容
本发明提供一种准确度高、临床意义明确的肝癌驱动基因高通量检测方法。
实现本发明目的的一种肝癌驱动基因高通量检测方法,包括如下步骤:
(1)根据人类基因组HG19,调取肝癌高频突变基因的外显子序列。常见肝癌基因列表如下:
调取外显子的序列包含列表基因的基因区域,以及己知各个转录本的启动子区域;
(2)对每个区域中非重复区域设计120bp的探针序列,每个序列沿着基因位置挪动设计,探针之间的挪动大小60bp;
(3)采用原位合成技术,在芯片上大量合成设计的探针,并利用多聚酶链式反应和转录的方法扩增出大量的带有生物素标记的探针,并制作肝癌基因高通量捕获试剂盒;
所述试剂盒检测方法包括如下步骤:
从患者血浆中提取100-300ng cfDNA (Cell Free DNA),然后利用肝癌常见基因探针将目的基因序列捕获出来,再用测序仪(Illumina Nextseq 500) 进行高通量测序,进而分析,找出与肝癌相关基因的所有突变信息,从而得到患者肺癌的变异情况,以达到准确基因诊断的目的。
利用测序仪(Illumina Nextseq 500)进行高通量测序,分析的过程包括单核昔酸多态性分析(SNP分析)过程、插入缺失标记分析(InDel分析)、基因融合和大片段扩增缺失分析流程;
所述单核昔酸多态性分析(SNP分析)过程包括如下步骤:
(1)测序仪(Illumina Nextseq 500)获取原始短序列并去除测序数据中的接头;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4)用GATK软件找出序列中所含有的单核昔酸多态性(SNP)的信息;
(5)统计测序结果信息,短序列数量、目标区域覆盖大小、平均测序深度等,过滤低质量值(小于25)和低覆盖度(小于10)的单核昔酸;
(6)利用ANNOVAR3、人类基因组数据库(HOMO_SAPIENS_NCBI_BUILD37.2)、dbSNP、COSMIC68信息对单核昔酸进行注释,确定突变位点发生的基因、坐标、mRNA位点、氨基酸改变、单核昔酸功能(错义突变/无义突变/可变剪切位点)、SIFT预测单核昔酸影响蛋白功能预测;
(7)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的单核昔酸作为候选的单核昔酸,在候选的单核昔酸中去除掉在dbSNP、HAPMAP、1000人类基因组、其他外显子测序项目中出现的单核昔酸。同时,过滤掉SIFT预测对蛋白功能无影响的单核昔酸作为最后疾病相关的候选单核昔酸;
所述插入缺失标记分析(InDel分析)流程包括如下步骤:
(1)把去除接头序列和低质量的测序数据用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t7 -e 40 -M 3 -f);
(2)用Pindel软件找出序列中所含有的括入/缺失(InDel)的信息(所用到参数:Pindel -k -H 50 -w 10 -E 0.6 -a 8 -M 10 -m 10 -d 100 -A 50 );
(3)利用CCDS、人类基因组数据库(HOMO_SAPIENS_NCBI_BUILD37.2)、dbSNP、COSMIC68信息对InDel进行注释,确定突变位点发生的基因、坐标、mRNA位点、编码区域序列的改变、对氨基酸的影响、InDel功能(氨基酸括入/氨基酸缺失/移码突变);
(4)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的InDel作为候选的InDels,在候选的InDels中去除掉在dbSNP、其他外显子测序项目中出现的InDel,最后筛选出疾病相关的候选InDels;
所述大片段扩增确和缺失分析流程包括如下步.:
(1)测序仪(Illumina Nextseq 500)获取原始短序列;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4) 用ONCOCNV6.4统计目标区域覆盖大小,然后根据每个目标区域位点的位置信息为横坐标,每个位置相应的覆盖度为纵坐标作图,得出扩增和缺失的分析图。
所述基因融合分析流程包括如下步.:
(1)测序仪(Illumina Nextseq 500)获取原始短序列;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4)用factera检测多个基因基因的全部或一部分的编码区首尾相连。包括染色体易位、中间缺失或染色体倒置所致的结果,得到基因融合结果。
本发明的一种循环肿瘤DNA肝癌基因的检测方法的有益效果如下:
本发明的一种循环肿瘤DNA肝癌基因的检测方法,是以DNA文库的方式进行捕获测序,因此,样品DNA的部分降解对最后的结果几乎不会产生影响,而过量的探针,保证了对目标片段的完全捕获。所以,即使可能出现的某些丰度较低的变异DNA也能被捕获并检测出来,其灵敏度相对于其他方法要高出许多。本方法以患者血浆循环肿瘤DNA为样本,采取高深度测序的方法(10000X)检测其驱动基因的突变情况。相比于临床上传统的病理学诊断而言,本方法有以下优势:
1、无创检测;本方法只需抽取患者血液(10ml)既可完成整个的诊断流程,而不需要取个体患处的组织或细胞。
2、预测性强;临床诊断通常只能在患者身体己经出现病理改变(出现肿块,蛋白改变)时做出诊断,因此,临床上肿瘤患者的发现,特别是恶性肿瘤,往往发现时病人己经处于中晚期。而基因诊断不仅仅对己患病的个体进行准确的诊断。
3、采取高深度测序,平均要达到10000X,保证低频突变能够检测出,能检测每个位点0.1%以上的及基因突变。准确度高;本方法的提供的肺癌基因捕获试剂盒覆盖了目前最新的肺癌高频突变基因的所有外显子区域。
具体实施方式
本发明的一种循环肿瘤DNA肝癌基因的检测方法,包括如下步骤:
1.样本文库制备:
(1)、吸取60ul纯化后的DNA(cfDNA)和40ul Endprep mix于PCR单管中,吹打混匀,短暂离心,进行末端修复,反应条件为30℃、30min。
(2)、将磁珠提前半个小时拿出以恢复到室温,震荡混匀,取180ul磁珠加入末端修复好的PCR产物中,抽打混匀,温室孵育5min。
(3)、将PCR单管放于磁力架上,静置5min,去上清,保持PCR管在磁力架上,加入80%乙醇(新配)200ul漂洗,静置30s后去上清,80%乙醇漂洗两遍,室温干燥10min(还有2min时用小抢头吸残液),直至无乙醇残留。
(4)、将PCR管从磁力架中取出,加入20ul超纯灭菌水洗脱,轻弹混匀,静置5min,短暂离心,置于磁力架上,待溶液澄清后吸上清17.5ul于新的PCR管中,再在PCR管中加入12.5ul dA-Tailing(解冻后放4℃,离心后再使用),吹打混匀。进行末端加尾,反应条件为37℃ 30min, 70℃ 5min,4℃ 5min。
(5)、在上步加尾后的产物中加入Ligation mix 和adapter各2.5ul,吹打混匀,短暂离心,放入PCR仪中,反应条件为30℃ 10min ,20℃ 20min。
(6)、在上步连接产物中加入15ul的超纯灭菌水,再加入50ul磁珠,震荡混匀,室温孵育5min,短暂离心,放入磁力架上,静置5min,澄清后去上清,保持PCR管在磁力架上,用200ul的80%乙醇(新配)漂洗两次,室温干燥10min,直至无乙醇残留。
(7)、将PCR管从磁力架上取出,加入25ul超纯灭菌水,轻弹混匀,静置5min,短暂离心,将PCR单管置于磁力架上,静置5min后小心吸取上清22ul于新的PCR管中。
(8)、将PCR primer mix和Amplification mix 2解冻后颠倒混匀,在上步连接产物中加入3ul的primer mix和25ul的Amplification mix 2,吹打混匀,短暂离心,进行PCR扩增,反应条件为
第一步95℃ 3min
第二步98℃ 20s
第三步60℃ 15s
第四步72℃ 30s
第五步72℃ 5min
第六步Hold at 4℃
第二----四步扩增15cycle。
(9)、用上步扩增产物进行琼脂糖凝胶电泳,观察条带大小及亮度。
(10)、将磁珠提前半个小时拿出以恢复到室温,震荡混匀,取50ul加入PCR扩增产物中,震荡混匀,室温孵育5min,短暂离心,置于磁力悬架上,静置5min,待溶液澄清后去上清,用200ul的80%乙醇(新配)漂洗两次,室温干燥10min,直至无乙醇残留。
(11)、将PCR管从磁力架中取出,加入40ul超纯灭菌水洗脱,轻弹混匀,静置5min,短暂离心,置于磁力架上,待溶液澄清后吸上清38ul于新的EP管中,用Qubit测溶液浓度并记录,-20℃保存。
2.样本富集液相杂交
(1)、用于杂交的建库后DNA总量为500ng,根据检测的溶液浓度计算杂交时加入EP管的溶液体积。
(2)、在上步的EP管中加入2.5ul的cot-1 DNA,2.5ul的SS DNA,2ul的P5(500um),2ul的P7-x(500um),短暂离心。
(3)、将上步混合好的溶液在60℃真空干燥30 min,直至无液体残留。
(4)、在上步干燥后的EP管中加入10ul超纯灭菌水,轻弹混匀,水化10 min,短暂离心,取10ul加入PCR八连排中。
(5)、将上步PCR八连排放入PCR仪中进行PCR反应,反应条件为
第一步95℃ 5min
第二步65℃ 5min
第三步65℃ 无限
在95℃结束后,将1ul RNAsin和5ul RNAbait混合(每个样品所取量),在65℃ 2.5min时将RNAsin和RNAbait混合液以及杂交液放入65℃PCR仪中,在65℃ 5min结束后先在八连排每个孔中加入6ul的混合液,加入时吹打混匀,2.5min后加入10ul(每个样品所取量)杂交液,加入后吹打混匀,65℃过夜。
3.捕获
(1)、洗磁珠,取Dnabeads磁珠震荡混匀,根据样品数取磁珠(每个样品取30ul磁珠),用200ul beads洗液洗3次。
(2)、悬浮磁珠于binding buffer中,每个样品取165ul的binding buffer。
(3)、吸取binding buffer和磁珠的结合液165ul于杂交样品中,震荡45min(每隔5min上下颠倒混匀),45min后置于磁力架上,待澄清后吸弃上清。
(4)、加165ul的洗液1于上步磁珠中,洗15min,每5min颠倒混匀一次,15min后置于磁力架上,待澄清后吸弃上清。
(5)、将洗液2提前放入70℃PCR仪预热,加165ul的洗液2于上步磁珠中,颠倒混匀,在70℃洗10min,10min后置于磁力架上,待澄清后吸弃上清,洗液2洗3次。
(6)、在上步八连排的每个孔中加入22 ul的超纯灭菌水,3ul的primer和25ul的Amplification mix2,颠倒混匀,进行PCR扩增,扩增程序为
第一步95℃ 3min
第二步98℃ 20s
第三步60℃ 15s
第四步72℃ 30s
第五步72℃ 5min
第六步Hold at 4℃
第二----四步扩增15cycle。
(7)、将上步八连排放在磁力架上,待溶液澄清后,吸上清进行琼脂糖凝胶电泳,观察条带大小及亮度。
(8)、将磁珠提前半个小时拿出以恢复到室温,震荡混匀,吸取上步八连排中的上清于新的PCR单管中,吸取50ul磁珠加入PCR单管中,震荡混匀,室温孵育5min,短暂离心,置于磁力悬架上,静置5min,待溶液澄清后去上清,用200ul的80%乙醇(新配)漂洗两次,室温干燥10min,直至无乙醇残留。
(9)、将PCR管从磁力架中取出,加入40ul超纯灭菌水洗脱,轻弹混匀,静置5min,短暂离心,置于磁力架上,待溶液澄清后吸上清38ul于新的EP管中,用Qubit测溶液浓度并记录,-20℃保存。
以上得到的DNA通过Illumina Nextseq 500测序,得到测序的数据。
4.SNP分析流程
(1)测序仪(Illumina Nextseq 500)获取原始短序列并去除测序数据中的接头;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4)用GATK软件找出序列中所含有的单核昔酸多态性(SNP)的信息;
(5)统计测序结果信息,短序列数量、目标区域覆盖大小、平均测序深度等,过滤低质量值(大于25)和低覆盖度(大于10)的单核昔酸;
(6)利用ANNOVAR3、人类基因组数据库(HOMO_SAPIENS_NCBI_BUILD37.2)、dbSNP、COSMIC68信息对单核昔酸进行注释,确定突变位点发生的基因、坐标、mRNA位点、氨基酸改变、单核昔酸功能(错义突变/无义突变/可变剪切位点)、SIFT预测单核昔酸影响蛋白功能预测;
(7)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的单核昔酸作为候选的单核昔酸,在候选的单核昔酸中去除掉在dbSNP、HAPMAP、1000人类基因组、其他外显子测序项目中出现的单核昔酸。同时,过滤掉SIFT预测对蛋白功能无影响的单核昔酸作为最后疾病相关的候选单核昔酸;
5.InDel分析流程
(1)把去除接头序列和低质量的测序数据用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t7 -e 40 -M 3 -f);
(2)用Pindel软件找出序列中所含有的括入/缺失(InDel)的信息(所用到参数:Pindel -k -H 50 -w 10 -E 0.6 -a 8 -M 10 -m 10 -d 100 -A 50 );
(3)利用CCDS、人类基因组数据库(HOMO_SAPIENS_NCBI_BUILD37.2)、dbSNP、COSMIC68信息对InDel进行注释,确定突变位点发生的基因、坐标、mRNA位点、编码区域序列的改变、对氨基酸的影响、InDel功能(氨基酸括入/氨基酸缺失/移码突变);
(4)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的InDel作为候选的InDels,在候选的InDels中去除掉在dbSNP、其他外显子测序项目中出现的InDel,最后筛选出疾病相关的候选InDels;
6.大片段扩增缺失分析流程
(1)测序仪(Illumina Nextseq 500)获取原始短序列;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4) 用ONCOCNV6.4统计目标区域覆盖大小,然后根据每个目标区域位点的位置信息为横坐标,每个位置相应的覆盖度为纵坐标作图,得出扩增和缺失的分析图。
7.基因融合分析流程
(1)测序仪(Illumina Nextseq 500)获取原始短序列;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4)用factera检测多个基因基因的全部或一部分的编码区首尾相连。包括染色体易位、中间缺失或染色体倒置所致的结果,得到基因融合结果。
上面所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神前提下,本领域普通工程技术人员对本发明技术方案做出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。

Claims (3)

1.一种用于循环肿瘤DNA肝癌驱动基因高通量检测方法,包括肝癌高频突变基因扫描试剂盒的方法和试剂盒的筛查方法。
2.肝癌高频突变基因扫描试剂盒的方法包括如下步骤:
(1)根据人类基因组HG19,调取肺癌高频突变基因的外显子序列和内含子序列;
常见肝癌基因列表如下:
TP53 CTNNB1 AXIN1 ALB ARID2 ARID1A ACVR2A NFE2L2 RPS6KA3 KEAP1 CDKN2A CDKN1A RB1 TSC2 HNF1A ZNRF3 PTEN TSC1 HGF APOB ATM SETD2 SMARCA2 JAK3 TERT DAPK1 CHD7 CREBBP CCND1 CDKN2B FGF19 FGF3 FGF4 PIK3CA KRAS NRAS HRAS EGFR MET ERBB2 BRAF KDR PDGFRB KIT FLT3 AXIN2 MYC JAK1 JAK2 STAT3 IGF1R IGF2R APC AKT1 PDGFRA VEGFA IRF2 CCNE1 MTOR FGFR1 FGFR2 FGFR3 FGFR4 ARID1B PBRM1 TCF7L2 EP300 ATR TTN FGF5 UBE3C FLT1 KMT2B KMT2C KMT2D EGFR BRAF DDR2 MEK PDGFRA ALK RET NTRK1 PIK3CA TP53 ROS1 KRAS FGFR3 PTEN MET NRAS RICTOR FGFR1 HER2 HRAS AKT1 BIM
调取的序列包含列表基因的基因区域,以及己知各个转录本的启动子区域和内含子区域;
(2)对每个区域中非重复区域设计120bp的探针序列,每个序列沿着基因位置挪动设计,探针之间的挪动大小60bp;
(3)采用原位合成技术,在芯片上大量合成设计的探针,并利用多聚酶链式反应和转录的方法扩增出大量的带有生物素标记的探针,并制作肝癌基因扫描试剂盒;
所述试剂盒筛查方法包括如下步骤:
从病人血液中提取100-300ng cfDNA (Cell Free DNA),然后利用肝癌常见基因探针将目的基因序列捕获出来,再用测序仪(Illumina Nextseq 500)进行高通量测序,进而分析,找出与肝癌相关基因的所有突变信息,从而得到患者肝癌的变异情况,以达到准确基因诊断的目的。
3.根据权利要求1所述的一种肝癌高频突变基因的筛查方法,其特征在于:采用测序仪进行高通量测序,分析的过程包括单核昔酸多态性分析过程、所述括入缺失标记分析流程、大片段扩增确实分析流程和基因融合分析流程;
所述单核昔酸多态性分析过程包括如下步骤:
(1)测序仪获取原始短序列;
(2)去除测序数据中的接头和低质量数据;
(3)把短序列定位到人类基因组数据相应的位置上;
(4)统计测序结果信息,短序列数量、目标区域覆盖大小、平均测序深度等;
(5)过滤低质量值和低覆盖度的单核昔酸;
(6)利用CCDS、人类基因组数据库、dbSNP信息对单核昔酸进行注释,确定突变位点发生的基因、坐标、mRNA位点、氨基酸改变、单核昔酸功能、SIFT预测单核昔酸影响蛋白功能预测;
(7)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的单核昔酸作为候选的单核昔酸,在候选的单核昔酸中去除掉在dbSNP、HAPMAP、1000人类基因组、其他外显子测序项目中出现的单核昔酸;同时,过滤掉SIFT预测对蛋白功能无影响的
单核昔酸作为最后疾病相关的候选单核昔酸;
所述括入缺失标记分析过程包括如下步骤
(1)把去除接头序列和低质量的测序数据比对到人类基因组上;
(2)找出序列中所含有的括入/缺失的信息;
(3)利用CCDS、人类基因组数据库、dbSNP信息对InDel进行注释,确定突变位点发生的基因、坐标、mRNA位点、编码区域序列的改变、对氨基酸的影响、InDel功能;
(4)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的InDel作为候选的InDels,在候选的InDels中去除掉在dbSNP、其他外显子测序项目中出现的InDel,最后筛选出疾病相关的候选InDels;
所述大片段扩增确实分析流程包括如下步骤:
(1)测序仪获取原始短序列;
(2)去除测序数据中的接头和低质量数据;
(3)把短序列定位到人类基因组数据相应的位置上;
(4)统计目标区域覆盖大小,然后根据每个目标区域位点的位置信息为横坐标,每个位置相应的覆盖度为纵坐标作图,得出扩增和缺失的分析图;
所述括入缺失标记分析过程包括如下步骤
(1)把去除接头序列和低质量的测序数据比对到人类基因组上;
(2)找出序列中所含有的括入/缺失的信息;
(3)利用CCDS、人类基因组数据库、dbSNP信息对InDel进行注释,确定突变位点发生的基因、坐标、mRNA位点、编码区域序列的改变、对氨基酸的影响、InDel功能;
(4)根据疾病样品和正常样品信息,选出疾病样品所共有的而在正常组中不存在的基因融合作为候选的融合基因,在候选的融合基因中去除掉在dbSNP、其他外显子测序项目中出现的融合基因,最后筛选出疾病相关的候选融合基因;
所述基因融合分析流程包括如下步.:
(1)测序仪(Illumina Nextseq 550)获取原始短序列;
(2)用Trim Galore进行QC并去除低质量数据;
(3)把短序列用BWA-Burrows-Wheeler Alignment(BWA)软件定位到人类基因组数据相应的位置上(所用到参数: bwa aln -L -l 40 -i 10 -k 2 -t 7 -e 40 -M 3 -f);
(4)用factera检测多个基因基因的全部或一部分的编码区首尾相连,包括染色体易位、中间缺失或染色体倒置所致的结果,得到基因融合结果。
CN201610661868.5A 2016-08-12 2016-08-12 一种循环肿瘤dna肝癌驱动基因高通量检测方法 Pending CN107723352A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610661868.5A CN107723352A (zh) 2016-08-12 2016-08-12 一种循环肿瘤dna肝癌驱动基因高通量检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610661868.5A CN107723352A (zh) 2016-08-12 2016-08-12 一种循环肿瘤dna肝癌驱动基因高通量检测方法

Publications (1)

Publication Number Publication Date
CN107723352A true CN107723352A (zh) 2018-02-23

Family

ID=61200890

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610661868.5A Pending CN107723352A (zh) 2016-08-12 2016-08-12 一种循环肿瘤dna肝癌驱动基因高通量检测方法

Country Status (1)

Country Link
CN (1) CN107723352A (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109295228A (zh) * 2018-10-18 2019-02-01 山东大学第二医院 一种检测血浆游离dna中肺癌相关突变基因的试剂盒及其应用
CN110669823A (zh) * 2018-07-03 2020-01-10 中国医学科学院肿瘤医院 一种同时检测多种肝癌常见突变的ctDNA文库构建和测序数据分析方法
CN110714075A (zh) * 2018-07-13 2020-01-21 立森印迹诊断技术(无锡)有限公司 一种用于检测肺肿瘤良恶性程度的分级模型及其应用
CN110964821A (zh) * 2019-12-17 2020-04-07 复旦大学附属中山医院 一种预测肝癌转移模式及风险的检测panel及其应用
CN111549131A (zh) * 2020-05-06 2020-08-18 中国科学院合肥物质科学研究院 一种基于循环肿瘤dna和肿瘤标记物的亚洲人多肿瘤早筛系统
CN111763729A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于braf基因全外显子二代测序的探针组合物、试剂及控制系统
CN111763720A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于tert基因全外显子二代测序的探针组合物、试剂及控制系统
CN111763719A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于tert基因全外显子二代测序的多重pcr引物组合物、试剂及控制系统
CN112176419A (zh) * 2019-10-16 2021-01-05 中国医学科学院肿瘤医院 一种检测ctDNA中肿瘤特异基因的变异和甲基化的方法
CN113186291A (zh) * 2021-05-26 2021-07-30 嘉兴允英医学检验有限公司 基于多重pcr的引物组和试剂盒

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103290137A (zh) * 2013-06-26 2013-09-11 北京迈基诺基因科技有限责任公司 一种肿瘤易感基因的筛查方法
CN104732116A (zh) * 2015-03-13 2015-06-24 西安交通大学 一种基于生物网络的癌症驱动基因的筛选方法
CN105063208A (zh) * 2015-08-10 2015-11-18 北京吉因加科技有限公司 一种血浆中游离的目标dna低频突变富集测序方法
WO2016095093A1 (zh) * 2014-12-15 2016-06-23 天津华大基因科技有限公司 肿瘤筛查方法、目标区域变异检测方法和装置
CN105779434A (zh) * 2014-12-15 2016-07-20 天津华大基因科技有限公司 试剂盒及其用途
CN105780129A (zh) * 2014-12-15 2016-07-20 天津华大基因科技有限公司 目标区域测序文库构建方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103290137A (zh) * 2013-06-26 2013-09-11 北京迈基诺基因科技有限责任公司 一种肿瘤易感基因的筛查方法
WO2016095093A1 (zh) * 2014-12-15 2016-06-23 天津华大基因科技有限公司 肿瘤筛查方法、目标区域变异检测方法和装置
CN105779434A (zh) * 2014-12-15 2016-07-20 天津华大基因科技有限公司 试剂盒及其用途
CN105780129A (zh) * 2014-12-15 2016-07-20 天津华大基因科技有限公司 目标区域测序文库构建方法
CN104732116A (zh) * 2015-03-13 2015-06-24 西安交通大学 一种基于生物网络的癌症驱动基因的筛选方法
CN105063208A (zh) * 2015-08-10 2015-11-18 北京吉因加科技有限公司 一种血浆中游离的目标dna低频突变富集测序方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110669823A (zh) * 2018-07-03 2020-01-10 中国医学科学院肿瘤医院 一种同时检测多种肝癌常见突变的ctDNA文库构建和测序数据分析方法
CN110669823B (zh) * 2018-07-03 2022-05-24 中国医学科学院肿瘤医院 一种同时检测多种肝癌常见突变的ctDNA文库构建和测序数据分析方法
CN110714075A (zh) * 2018-07-13 2020-01-21 立森印迹诊断技术(无锡)有限公司 一种用于检测肺肿瘤良恶性程度的分级模型及其应用
CN110714075B (zh) * 2018-07-13 2024-05-03 立森印迹诊断技术(无锡)有限公司 一种用于检测肺肿瘤良恶性程度的分级模型及其应用
CN109295228A (zh) * 2018-10-18 2019-02-01 山东大学第二医院 一种检测血浆游离dna中肺癌相关突变基因的试剂盒及其应用
CN109295228B (zh) * 2018-10-18 2021-06-22 山东大学第二医院 一种检测血浆游离dna中肺癌相关突变基因的试剂盒及其应用
CN111763720A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于tert基因全外显子二代测序的探针组合物、试剂及控制系统
CN111763719A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于tert基因全外显子二代测序的多重pcr引物组合物、试剂及控制系统
CN111763729A (zh) * 2019-04-01 2020-10-13 长沙金域医学检验实验室有限公司 用于braf基因全外显子二代测序的探针组合物、试剂及控制系统
CN112176419A (zh) * 2019-10-16 2021-01-05 中国医学科学院肿瘤医院 一种检测ctDNA中肿瘤特异基因的变异和甲基化的方法
CN110964821A (zh) * 2019-12-17 2020-04-07 复旦大学附属中山医院 一种预测肝癌转移模式及风险的检测panel及其应用
CN111549131A (zh) * 2020-05-06 2020-08-18 中国科学院合肥物质科学研究院 一种基于循环肿瘤dna和肿瘤标记物的亚洲人多肿瘤早筛系统
CN113186291A (zh) * 2021-05-26 2021-07-30 嘉兴允英医学检验有限公司 基于多重pcr的引物组和试剂盒
CN113186291B (zh) * 2021-05-26 2022-04-29 嘉兴允英医学检验有限公司 基于多重pcr的引物组和试剂盒

Similar Documents

Publication Publication Date Title
CN107723352A (zh) 一种循环肿瘤dna肝癌驱动基因高通量检测方法
CN105256051B (zh) 一种用于检测先天性巨结肠及相关综合征的致病/易感基因的探针组及试剂盒
CN107475370A (zh) 用于肺癌诊断的基因群和试剂盒及诊断方法
CN108753967A (zh) 一种用于肝癌检测的基因集及其panel检测设计方法
WO2015101175A1 (en) Method and kit for non-invasively detecting egfr gene mutaions
WO2016095093A1 (zh) 肿瘤筛查方法、目标区域变异检测方法和装置
CN111575380B (zh) 多基因检测用的探针库、杂交试剂盒和多基因检测的方法
CN107502654A (zh) 用于实体瘤靶向用药指导的多基因富集和检测方法
CN110643680B (zh) 适用于超微量dna测序的接头及其应用
CN109385666A (zh) 淋巴瘤基因捕获芯片及其应用
CN110055331B (zh) 一种用于膀胱癌辅助诊断或筛查的试剂盒及其应用
CN107723353A (zh) 一种白血病驱动基因的高通量检测方法
CN107723350A (zh) 一种遗传性乳腺癌突变基因高通量筛查方法
CN106701920A (zh) 一种预测结直肠癌肝转移的试剂盒和使用方法
CN115443341A (zh) 分析无细胞核酸的方法及其应用
CN109988838A (zh) 一种可对肿瘤靶向治疗相关靶点和免疫治疗相关tmb及msi同时进行检测的体系
CN108753954A (zh) 痴呆相关基因的捕获探针组、试剂盒、文库构建方法和用途
CN107760779A (zh) 肺动脉高压相关的突变型bmp9基因及其应用
BR112019013391A2 (pt) Adaptador de ácido nucleico, e, método para detecção de uma mutação em uma molécula de dna circulante tumoral (ctdna) de fita dupla.
CN107557461A (zh) 一种用于肝癌易感基因早筛的核酸质谱的检测方法
CN107723351A (zh) 一种循环肿瘤dna肺癌驱动基因的高通量检测方法
CN107881232A (zh) 探针组合物及基于ngs方法检测肺癌和结直肠癌基因的应用
CN110093423A (zh) 基因标志物在肺腺癌诊治中的用途
CN108753953A (zh) 一种人dmd基因外显子pcr扩增系统、检测试剂盒及其应用的方法
CN109897900B (zh) Epb42基因在肝癌sbrt疗效评估中的应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180223