本发明属于生物分子检测领域，更具体地说，涉及一种基于高通量测序技术的肺癌靶向治疗分子诊断引物、试剂盒及诊断方法。
背景技术：
根据国家癌症中心公布的最新数字显示，我国恶性肿瘤发病及死亡居首位的是肺癌，每年新发病例约73.3万，死亡约59.1万，占所有因癌死亡人数的35％以上，超越了交通事故、自然灾害等，成为了影响我国人口死亡的最主要因素。对比过去三十年统计数据，我国肺癌的死亡率上升了465％。而根据美国癌症协会公布的《2017年度美国癌症数据调查报告》显示，美国肺癌发病率和死亡率持续显著下降。该报告分析认为，美国肺癌发病率和死亡率呈现双下降主要得益于三大因素：积极控烟、推广筛查和新型疗法的应用。大量学者和医务人员也认为推广筛查和新型疗法的应用将大大推进我国的肺癌筛查、诊断和治疗水平。
目前，我国临床上对肺癌仍以组织病理学诊断为确诊和治疗的依据，使用的治疗方法是根据其分期采取手术或手术结合放化疗进行治疗。但由于肺癌缺乏有效的早期诊断手段以及受国民健康意识的影响，70％以上的患者在确诊时已是中晚期，错过了最佳手术时期。而肺癌各种化疗方案的总体有效性仅为30％左右，同时有些患者无法耐受化疗和放疗以及化疗后很快耐药，从而导致肺癌患者确诊后5年生存率很低，晚期患者五年生存期小于15％。
近年来，分子靶向治疗逐渐受到人们的关注且在欧美国家已广泛应用于临床。分子靶向治疗是使用药物针对已经明确的致癌基因，特异地杀死肿瘤细胞，因此分子靶向治疗的副作用更小和特异性更高，也为晚期肺癌患者提供了新的有效治疗手段。这其中最典型的例子是肺腺癌的egfr突变靶向治疗。在我国，约30％非小细胞肺癌患者存在egfr突变。对这类患者采用egfr激酶抑制剂-吉非替尼、厄洛替尼、阿法替尼，有效率可以达到70％～80％。过去十年针对肺癌的特定基因突变开展的靶向治疗就已让晚期肺癌患者中位生存期显著增加2.4倍，从此前的14.1个月延长至33.5个月，相较于传统的治疗方法优势明显。目前，临床上关于肺癌的靶向药物，然而除激酶抑制剂-吉非替尼、厄洛替尼、阿法替尼以外，还有多种基因特异位点的抑制剂对肺癌显示出较好的抑制效果。然而，并非所有的肺癌患者对分子靶向治疗均表现出较好的疗效，比如现有资料表明egfr基因突变与egfr激酶抑制剂敏感性相关，egfr基因突变患者对egfr激酶抑制剂吉非替尼的有效率达60％以上，而egfr基因野生型(即未突变型)患者的有效率仅为10％～15％。因此，明确肺癌组织的突变状况对靶向药物的临床使用和疗效预测有重要的指导意义。而分子靶向药物的使用需要明确肿瘤患者存在可靶向治疗的基因突变，所以基因突变的检测是确定治疗方案的关键。
临床中肺癌的诊断样品70％是小的活检样品，其首先被用于临床病理学诊断，剩下的有限组织若要用作分子诊断，则对其核酸含量及提取方法都有较高的技术操作要求。且病理形态相同的肺癌，由于分子遗传学改变，在分子水平上呈现高度异质，这种组织的异质性使得一些活检样本无法反应肿瘤基因突变的全貌，从而导致对靶向治疗的反应差别很大。同时，肿瘤异质性也是不断变化的，无法通过一两次的活检或组织样本来检测患者的病情。
高通量测序技术通过多基因、多位点的集成式检测，可一次性检测与肿瘤靶向治疗密切相关的基因位点。它所需的核酸要求不高，是应用于肿瘤基因检测的一个很好的选择。此外，肿瘤是基因组疾病，单一基因的检测难以满足指导个体化治疗的需要。而以高通量测序和大数据分析为基础，对肺癌进行多基因突变检测，能够精确的地给出患者癌组织的基因突变信息，从而指导医生为病人制定针对性的个性化靶向治疗方案。同时，也可用于肺癌病人的早期诊断和筛查。目前，在美国高通量测序技术已广泛用于癌症的筛查、诊断和治疗过程中，对美国癌症的死亡率下降做出了突出的贡献。
然而，在国内市场上对于肺癌的检测主要停留在单个或几个位点的检测水平，检测方法一般使用荧光pcr、芯片分析等方法，例如，中国专利申请号为201310200169.7，申请公布日为2013年9月11日的专利文件公开了一种用于检测alk基因表达的探针、引物及试剂盒，包括以下序列：seqidno1-seqidno12，该发明的引物和探针可以特异检测alk基因表达差异，以确定是否存在alk基因融合，该发明建立的用于检测alk基因表达差异的实时荧光pcr体系，通过检测5’端和3’端基因表达量的差异，提供alk基因表达差异外显子的定性评估，适用于肺癌患者在进入个性化靶向治疗之前。中国专利申请号为201310068052.8，申请公布日为2013年6月5日的专利文件公开了一种早期肺癌多位点关联基因的引物、探针、试剂盒及其制备方法和检测方法。该发明针对肺癌的早期检测与之相关联的egfr，chrna3，chrna5及ercc1基因内部的rs2239680，rs16969968，rs12910984及rs2298881位点的基因型进行分型，利用以上4个位点联合作为肺癌早期诊断的标志物，提高了检测的真实性和准确性。然而，由于癌症是基因组疾病，单一或几个基因突变位点的检测难以满足指导个体化治疗的需要。已经报道的几种检测肺癌的方法只包含了极少的几个基因突变位点，不能有效的反应患者癌症突变的全貌，难以用于临床指导，且临床发现单一基因的靶向用药多数存在治疗无效或很快产生药物抗性等不利反应，因此，我们需要一种更加全面有效的检测手段。同时，对癌症进行全基因组或全外显子测序可以获得海量的基因突变数据。然而，这些测序费用极其昂贵，后期数据分析要求极高、耗时长达数月，且由于我们对绝大多数基因突变在癌症中的功能认识有限，目前仅针对极少的一部分基因研发出了靶向药物。因此，全基因组或全外显子测序检测肺癌的方法无法应用于临床需求。根据前期研究表明，肺癌存在多种基因突变，其中最常见的突变包括致癌基因akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten，更重要的是这些突变临床上均已开发或正在测试靶向药物，对于这些基因突变的检测将更好的指导医生为病人制定针对性的个性化靶向治疗方案。
综上所述，现有的肺癌的检测方法都存在一定的缺陷，例如结果不易判读、重复性差、假阳性和假阴性多、检测基因少、通量低等问题。而现在市场上存在的一些基于高通量测序技术检测肺癌的方法，存在对样品的要求较高、目的性不强、操作繁琐、引物或探针设计复杂、检测费用高等缺点，特别是不能指导靶向用药，因此急需一种新型的检测方法以满足市场和医疗的需求。
技术实现要素：
1.要解决的问题
针对现有的肺癌基因突变位点难以准确确定、重复性差、假阳性和假阴性多、通量低等问题，本发明提供一种基于高通量测序技术确定肺癌基因突变位点的引物、试剂盒及方法，用于检测确定16个致癌基因的316个突变的方法，所述的致癌基因为akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten，所述突变可以是一个或多个碱基的置换、插入和/或缺失。
2.技术方案
为了解决上述问题，本发明所采用的技术方案如下：
本发明针对肺癌靶向治疗分子诊断的16个致癌基因的316个突变位点的特定序列设计了特异性引物(seqidno.1～seqidno.82)，检测的基因为akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten；可以扩增致癌基因akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten的突变位点区域，扩增产物经高通量测序获得其序列信息，从而检测确定出样品中这些基因的突变状况。本发明所用引物(seqidno.1～seqidno.82)，可使用多重pcr扩增技术在一管中同时有效的扩增目标序列，所得pcr产物纯化后经二轮pcr引入测序接头和barcode后，用于制备测序文库和上机进行高通量测序，所构建的测序文库经nanodrop2000、agilentbioanalyzer2100、labchipgxtouch、qubit3.0和qpcr检测合格后适用于多种高通量测序平台，包括illuminanextseq500/550测序仪、illuminamiseq测序仪、illuminaminiseq测序仪、illuminahiseq2500/3000/4000测序仪和illuminanovaseq等测序仪。
一种确定肺癌基因突变位点的特异性引物，引物序列为seqidno.1～seqidno.82中的一种或多种组合。引物序列的设计要结合特定基因片段，检测目的突变位点等因素，采用合适的算法设计得到。
一种肺癌靶向治疗分子诊断试剂盒，包括引物序列为seqidno.1～seqidno.82中的一种或多种组合。
表1肺癌16个靶向治疗分子诊断基因突变位点检测seqno.1～seqno.82引物序列
备注：序列号代表引物序列的编号，forwardprimer表示正向引物序列，reverseprimer代表反向引物序列。
一种基于高通量测序技术确定肺癌基因突变位点的方法，包括以下步骤：
a.从待测样品中提取基因组dna；
b.使用上述的基因特异性引物seqidno.1～seqidno.82，对步骤a中提取的基因组dna上的目的片段进行多重pcr扩增；
c.对步骤b中得到的多重pcr扩增产物进行纯化；
d.对步骤c中的纯化产物进行二轮pcr扩增，引入测序接头和barcode；
e.对步骤d中的二轮pcr产物进行纯化和质控，完成测序文库的构建；
f.对步骤e中构建的测序文库进行高通量测序，获得目的基因片段的序列信息；
g.将步骤f中的获得的目的基因片段的序列信息与正常的基因序列进行比对分析，得到目的基因片段的突变信息。
更进一步地，步骤a中的待测样品包括但不限于新鲜肿瘤组织、新鲜活检组织和石蜡包埋肿瘤组织；提取基因组dna时，使用但不限于组织dna提取试剂盒和石蜡包埋组织dna提取试剂盒。
步骤a中提取基因组dna后对dna进行质量检测，用于质量检测的仪器优选为nanodrop2000和qubit3.0，对dna的质量要求为：dna浓度大于5ng/μl，a260/a280>1.8，a260/a230>2.0；步骤b中的特异性引物seqidno.1～seqidno.82进行多重pcr扩增时同时加入
更进一步地，步骤b中所述的基因组dna上的目的片段是指分别包含肺癌16个靶向治疗分子诊断基因的316个突变位点的基因片段，16个靶向治疗分子诊断基因分别为akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten。16个靶向治疗分子诊断基因的316个突变位点相关信息如表2所示。
表216个肺癌靶向治疗分子诊断基因的316个突变位点和相应靶向药物信息
备注：基因表示基因名称，如akt1代表akt1基因；cosmicid是指基因突变在癌症基因组突变数据库中的编号，如肺癌中akt1基因一个最常见的突变在癌症基因组突变数据库中的编号为cosm33765；核苷酸变异代表编码序列的突变，如c.49g>a表示编码区第49位鸟苷酸突变成为腺苷酸；氨基酸变异表示氨基酸序列的突变，如p.e17k表示蛋白序列中第17位的谷氨酸突变为赖氨酸；位置则指该突变在染色体上的位置，如14:105246551-105246551表示突变发生在14号染色体105246551位点；此外，ins代表插入，del代表缺失，*表示终止子，？代表对编码蛋白影响未知。
更进一步地，步骤b中多重pcr扩增的反应体系为：2xmastermix(包含dna聚合酶、mgcl2和dntps)25μl、包含所述seqno.1～seqno.82引物的引物panel(10nmeach)10μl、待测dna样品10～20ng，用ddh2o加至50μl；步骤b中多重pcr扩增的反应条件为：95℃预变性3min；然后进行聚合酶链式扩增：95℃变性20s，60℃退火和延伸4min，共进行20～23个循环；最终72℃延伸4min。
步骤b中使用肺癌突变基因特异性引物seqidno.1～seqidno.82，对步骤a中提取的基因组dna上的目的片段进行多重pcr扩增，扩增产物经处理后进行高通量测序，用以检测表2中的基因突变状况。检测每个突变位点的引物对分别为：seqno.1和seqno.2检测akt1基因cosm33765突变，seqno.3和seqno.4检测alk基因cosm28056突变，seqno.5和seqno.6检测alk基因cosm28055、cosm28061和cosm28057突变，seqno.7和seqno.8检测araf基因cosm1742787和cosm5044705突变，seqno.9和seqno.10检测braf基因cosm308550、cosm475、cosm477、cosm478、cosm1127、cosm1583011、cosm473、cosm474、cosm476、cosm6137、cosm1130、cosm1126、cosm1125、cosm471、cosm470、cosm1121、cosm253330、cosm211600和cosm27639突变，seqno.11和seqno.12检测braf基因cosm460和cosm451突变，seqno.13和seqno.14检测cdkn2a基因cosm3092239、cosm13670、cosm13221、cosm28562、cosm1247967、cosm3906638、cosm13294、cosm13698、cosm1247961、cosm12479、cosm753746、cosm4774787、cosm13830、cosm3952632、cosm13811、cosm3433073、cosm3092261、cosm13701、cosm12547、cosm48297、cosm12481、cosm753741、cosm12484、cosm13489、cosm3092274、cosm13618、cosm12751和cosm4604974突变，seqno.15和seqno.16检测egfr基因cosm1451536和cosm21683突变，seqno.17和seqno.18检测egfr基因cosm21686、cosm21685、和cosm21687突变，seqno.19和seqno.20检测egfr基因cosm18441、cosm6252、cosm6253、cosm18425、cosm6239和cosm747426突变，seqno.21和seqno.22检测egfr基因cosm3734668突变，seqno.23和seqno.24检测egfr基因cosm4745530、cosm291998、cosm6241、cosm48921、cosm48923和cosm6240突变，seqno.25和seqno.26检测egfr基因cosm13553、cosm12429、cosm133630、cosm6224、cosm12374和cosm6213突变，seqno.27和seqno.28检测erbb2基因cosm1205571、cosm14060、cosm436499、cosm14062、cosm436500和cosm5802314突变，seqno.29和seqno.30检测fgfr3基因cosm24842突变，seqno.31和seqno.32检测fgfr3基因cosm719、cosm726、cosm720、cosm731、cosm1428730和cosm3993568突变，seqno.33和seqno.34检测kit基因cosm6142、cosm1177、cosm1180、cosm1328、cosm1190、cosm133763、cosm1192、cosm1329、cosm1198、cosm23418、cosm1210、cosm1211、cosm1217、cosm1216、cosm1219、cosm1221、cosm1235、cosm1239、cosm18896、cosm1246、cosm1247、cosm1248、cosm19216、cosm1252、cosm1253、cosm1255、cosm1256、cosm1258、cosm1257和cosm1260突变，seqno.35和seqno.36检测kit基因cosm1304、cosm19326和cosm12706突变，seqno.37和seqno.38检测kit基因cosm36053突变，seqno.39和seqno.40检测kit基因cosm1312、cosm1310、cosm1314、cosm12710、cosm1316、cosm133670、cosm12709、cosm19280、cosm1318、cosm19109、cosm3604523、cosm1320、cosm1321和cosm1322突变，seqno.41和seqno.42检测kras基因cosm19900、cosm19404和cosm19905突变，seqno.43和seqno.44检测kras基因cosm1212768、cosm1168052、cosm554、cosm555、cosm552、cosm553、cosm87298、cosm549、cosm550、cosm547和cosm546突变，seqno.45和seqno.46检测kras基因cosm1511785、cosm4384683、cosm543、cosm12703、cosm20818、cosm12721、cosm12722、cosm526、cosm530、cosm531、cosm4745557、cosm532、cosm534、cosm87281、cosm14209、cosm249888、cosm4387522、cosm513、cosm515、cosm527、cosm529、cosm5413585、cosm25801、cosm512、cosm520、cosm521、cosm522、cosm516、cosm517和cosm516突变，seqno.47和seqno.48检测map2k1基因cosm1235485、cosm2257207、cosm3690509、cosm3503329、cosm555604、cosm1562837、cosm5588244、cosm1725008、cosm2257208、cosm1235481、cosm5369532、cosm4756761、cosm1235478、cosm5520914、cosm1235479、cosm1678546和cosm5732210突变，seqno.49和seqno.50检测map2k1基因cosm4166154、cosm4166152、cosm4166153、cosm555601、cosm1315829、cosm235614、cosm1167912、cosm1315861、cosm1374186和cosm1235480突变，seqno.51和seqno.52检测mtor基因cosm1686998突变，seqno.53和seqno.54检测nras基因cosm30646、cosm33693、cosm53223、cosm12725、cosm579、cosm585、cosm586、cosm12730、cosm583、cosm584和cosm580突变，seqno.55和seqno.56检测nras基因cosm572、cosm574、cosm569、cosm564、cosm565、cosm566、cosm561、cosm562和cosm563突变，seqno.57和seqno.58检测pdgfra基因cosm133645和cosm739突变，seqno.59和seqno.60检测pdgfra基因cosm22414和cosm22415突变，seqno.61和seqno.62检测pik3ca基因cosm759、cosm17442、cosm760、cosm1041495、cosm761、cosm27133、cosm763、cosm12458、cosm764、cosm27374、cosm765、cosm766、cosm12459、cosm767、cosm24712和cosm5881065突变，seqno.63和seqno.64检测pik3ca基因cosm17444、cosm12461、cosm12590、cosm771、cosm36286、cosm12591、cosm29313、cosm773、cosm94984、cosm36288、cosm12592、cosm27504、cosm774、cosm775、cosm776、cosm1041525、cosm24714和cosm12597突变，seqno.65和seqno.66检测pten基因cosm428078、cosm3967199、cosm5020、cosm5270、cosm5275、cosm5878、cosm28916、cosm921056和cosm5232429突变，seqno.67和seqno.68检测pten基因cosm5311和cosm5328突变，seqno.69和seqno.70检测pten基因cosm5313和cosm685105突变，seqno.71和seqno.72检测pten基因cosm27933、cosm1474840、cosm1474841、cosm23629、cosm4948、cosm5322和cosm5960388突变，seqno.73和seqno.74检测pten基因cosm5149、cosm13116和cosm35665突变，seqno.75和seqno.76检测pten基因cosm5154、cosm1297486、cosm88109、cosm5810、cosm5159、cosm1173612、cosm428088、cosm921131、cosm5231546和cosm428089突变，seqno.77和seqno.78检测pten基因cosm86061、cosm4613228、cosm5312和cosm28899突变，seqno.79和seqno.80检测pten基因cosm921150、cosm921151、cosm5008和cosm5892突变，seqno.81和seqno.82检测pten基因cosm5310、cosm685100和cosm3807907突变。
更进一步地，步骤c中对多重pcr扩增产物纯化的方法为：使用agencourtampurexp磁珠进行纯化，50μlpcr反应体系中首先加入25μlagencourtampurexp磁珠除去基因组dna，然后在上清中加入30μlagencourtampurexp磁珠，经200μl70％乙醇洗涤两次后，彻底挥发乙醇。
更进一步地，步骤d中二轮pcr扩增引入测序接头和barcode的反应体系为：干燥的磁珠中加入2xmastermix25μl和包含不同barcode的接头引物3μl，加ddh2o至50μl；2xmastermix中包含dna聚合酶、mgcl2和dntps，接头引物的浓度为10nm，步骤d中二轮pcr扩增的反应条件为：95℃预变性3min；然后进行聚合酶链式扩增：95℃变性15s，58℃退火15s，72℃延伸30s，共进行6～8个循环；最终72℃延伸5min。设计不同的barcode，目的用来区分不同的样本，可以根据设计者的意愿随意加不同长度和序列来进行区分。
更进一步地，包含barcode的接头引物序列为seqidno.83～seqidno.95(序列如表3所示)，适用于illuminate测序仪。
表3用于illuminate测序仪测序文库构建的接头引物序列(seqidno.83～seqidno.95)
备注：adaptername代表接头的名称，sequence表示接头的核苷酸序列信息，所有接头的5’端都经过磷酸化修饰，*代表修饰引物中的磷酸酯键，斜体黑色字母标示barcode序列。
更进一步地，步骤e中二轮pcr扩增产物使用agencourtampurexp磁珠进行纯化，50μlpcr反应体系中加入45μlagencourtampurexp磁珠，经200μl70％乙醇洗涤两次后，彻底挥发乙醇后，加入20μlddh2o进行洗脱，获得测序文库；质控是指使用qubit3.0、nanodrop2000、agilentbioanalyzer2100、labchipgxtouch或qpcr对测序文库的dna进行质量检测，测序文库的dna需要满足以下条件：浓度大于20ng/μl，a260/a280>1.8，a260/a230>2.0，片段大小位于250bp～500bp之间。
本发明中构建的测序文库适用于多种高通量测序平台，包括illuminanextseq500/550测序仪、illuminamiseq测序仪、illuminaminiseq测序仪、illuminahiseq2500/3000/4000测序仪和illuminanovaseq等测序仪；步骤e中所得的测序文库按照不同测序方法的要求进行上机前处理，使得上样量相当的、不同barcode标记的测序文库处于测序方法需要的缓冲液体系中，能够直接用于上机测序。
3.有益效果
相比于现有技术，本发明的有益效果为：
(1)本发明提供的特异性引物，可以一次性检测确定肺癌靶向治疗分子诊断16个基因的316个突变位点的突变情况，具有特异性强、非特异性扩增程度低、各引物间相互干扰小、扩增效果好等优点，在进行多重pcr扩增时，这些引物混合溶解于同一管中同时加入，效率非常高；
(2)本发明提供了一种确定肺癌基因突变位点的方法，样品要求低，起始样品dna的量可以低至1ng，可以使用各种来源的样品进行dna的提取，包括且不限于新鲜肿瘤组织、新鲜活检组织和石蜡包埋肿瘤组织等，而且可以使用已发生降解的dna样品；
(3)本发明的检测通量高，一次可以检测十几个到几百个样品的16个基因的316个突变位点状况，高效便捷；
(4)本发明的确定方法，灵敏度高，检测灵敏度高达1％，检测结果明确客观，可以直接确定相关基因的具体突变位点，对于指导临床肺癌靶向治疗分子诊断、靶向用药以及用于癌症早期诊断或辅助诊断及筛查和癌症愈后监测有重要意义；
(5)本发明的检测确定方法安全快速，不含有毒有害物质，对实验人员和环境都无危害，检测速度快，仅需两到三天就可以获得检测结果；
(6)本发明的检测确定方法基于多重pcr技术和高通量测序技术，可以同时检测16个临床已有或正在测试靶向药物的基因的316个突变位点，使用的特异性引物和接头及优化的实验体系，使得检测方法简单易行、特异性高、重复性好、准确度高，对肺癌靶向治疗分子诊断基因的突变位点具有检测通量高、检测速度快、准确性和重复性好、安全性高等众多优点，为指导肺癌患者的个性化治疗，辅助肺癌的早期诊断及癌症愈后检测提供一种高效的方法。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施例对本发明进一步详细说明。应该理解，这些实施例仅用于说明本发明，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知技术的描述，以避免不必要的混淆本发明的概念。在阅读了本发明的内容后，本领域技术人员可以对本发明做各种改动，这些等价形式同样落入本发明所限定的范围。
实施例1
对某t2bn1m1期肺鳞癌活检组织进行检测确定基因突变的相关信息，步骤如下：
所用样本为t2bn1m1期肺癌新鲜活检组织，采用dna提取试剂盒提取肺癌新鲜活检组织中的基因组dna，应用seqno.1～seqno.82对提取得到的基因组dna进行多重pcr扩增后，采用illuminanextseq500测序仪进行测序分析，具体实施方法如下：
a.基因组dna提取：使用qiagendneasyblood&tissuekit，参照试剂盒说明书进行基因组dna的提取。获得的基因组dna进行凝胶电泳和使用nanodrop2000进行质量检测，要求基因组dna没有明显的降解，浓度大于10ng/μl，a260/a280>1.8，a260/a230>2.0，并使用qubit3.0进行精确定量。
b.多重pcr扩增：使用引物seqno.1～seqno.82对10ng基因组dna进行多重pcr扩增，获得akt1、alk、araf、braf、cdkn2a、egfr、erbb2、fgfr3、kit、kras、map2k1、mtor、nras、pdgfra、pik3ca和pten这16个目的基因的316个突变位点区域的目标扩增片段。
所用的多重pcr反应体系为：genereaddnaseqpanel5xpcrbuffer10μl、genereadhotstartaqdnapolymerase(6u/μl)2.5μl、含所述seqno.1～seqno.82引物的引物panel(10nmeach)10μl、待测dna样品10～20ng，用ddh2o加至50μl；多重pcr扩增的反应条件为：95℃预变性3min；然后进行聚合酶链式扩增：95℃变性20s，60℃退火和延伸4min，共进行22个循环；最终72℃延伸4min。
c.pcr产物纯化：步骤b中多重pcr扩增产物使用agencourtampurexp磁珠进行纯化，50μlpcr反应体系中首先加入25μlagencourtampurexp磁珠除去基因组dna，然后在上清中加入30μlagencourtampurexp磁珠，经200μl70％乙醇洗涤两次后，彻底挥发乙醇。
d.二轮pcr引入测序接头和barcode：所用的二轮pcr反应体系为：干燥的磁珠中加入genereaddnaseqpanel5xpcrbuffer10μl、genereadhotstartaqdnapolymerase
(6u/μl)2.5μl、含不同barcode的接头引物(10nm)3μl，加ddh2o至50μl；二次pcr扩增的反应条件为：95℃预变性3min；然后进行聚合酶链式扩增：95℃变性15s，58℃退火15s，72℃延伸30s，共进行7个循环；最终72℃延伸5min。
e.二轮pcr产物纯化和质控：使用agencourtampurexp磁珠进行纯化，50μlpcr反应体系中加入45μlagencourtampurexp磁珠，经200μl70％乙醇洗涤两次后，彻底挥发乙醇后，加入20μlddh2o进行洗脱，获得测序文库。测序文库的dna浓度大于20ng/μl，a260/a280>1.8，a260/a230>2.0，片段大小位于250bp～500bp之间(使用qubit3.0、nanodrop2000和labchipgxtouch对测序文库的dna进行质量检测)。
f.高通量测序：等量的dna测序文库使用illuminatelibrarynormalizationsolutions进行平衡后按照illuminatenextseq500测序仪操作进行上样测序，测序结束后，对获得的序列进行拼接后使用软件进行比对分析，筛选出发生突变的基因位点。
检测信息如下表4所示
检测结果如下表5所示
根据检测结果，确定了该t2bn1m1期肺鳞癌是由于egfr基因的编码区第2573位胸腺嘧啶核苷酸突变为鸟嘌呤核苷酸，导致蛋白序列中第858位的亮氨酸突变为精氨酸。本实施例中确定了基因以及氨基酸的突变位点，有利于t2bn1m1期肺鳞癌的靶向治疗方案确定。
实施例2
某转移性肺腺癌，经fish检测结果为egfr拷贝数扩增，erbb2拷贝数扩增。使用本发明的方法对其进行分子检测，确定具体的突变位点。
所用样本为石蜡包埋的肺癌组织，采用dna提取试剂盒提取石蜡包埋组织中的基因组dna，应用seqno.1～seqno.82对提取得到的基因组dna进行多重pcr扩增后，采用illuminanextseq500测序仪进行测序分析，具体实施方法如下：
a.基因组dna提取：使用qiagengenereaddnaffpekit，参照试剂盒说明书进行基因组dna的提取。获得的基因组dna使用nanodrop2000进行质量检测，要求dna浓度大于5ng/μl，a260/a280>1.8，a260/a230>2.0，并使用qubit3.0进行精确定量。
步骤b～f的操作与实施例1相同。
检测信息如下表6所示
检测结果如下表7所示
根据检测结果，确定了该转移性肺腺癌是由于egfr基因的编码区第2310-2311位之间插入了gggtta核苷酸序列，导致蛋白序列中第770-771位插入了甘氨酸和亮氨酸两个氨基酸；erbb2基因的编码区第2236位发生g核苷酸缺失并插入了ttat序列，导致蛋白序列中第746位甘氨酸的缺失并发生亮氨酸和半胱氨酸的插入。本实施例中确定了基因以及氨基酸的突变位点，有利于转移性肺腺癌的靶向治疗方案确定。
本发明对肺癌靶向治疗分子诊断16个基因的316个突变位点具有检测通量高、特异性强、样本来源广泛、不易发生污染且安全性高等众多优点，检测结果具有很好的重复性和准确性，对肺癌具有较好的辅助诊断和指导治疗作用。