CN109097466A - 与主动脉夹层疾病相关的snp位点及其应用 - Google Patents

Abstract

本发明公开了一种基于全外显子测序技术筛选与主动脉夹层疾病相关的SNP位点，所述SNP位点包括以下主动脉夹层易感性SNP位点中的一种或两种：PTK6_c.205G>A_p.E69K；TSHZ2_c.724C>T_p.R242C。进一步，本发明提供所述SNP位点在检测主动脉夹层疾病的检测产品中的应用。本发明通过发现了两个主动脉夹层易感SNP位点并研究了其在主动脉夹层早期预测的应用前景，本发明提供的主动脉夹层相关的基因位点可能是主动脉夹层的早期生物标志物，为进一步研究主动脉夹层的遗传分子机制，探索主动脉夹层早期防治的药物靶点提供了新的方向。

Description

与主动脉夹层疾病相关的SNP位点及其应用

技术领域

本发明涉及生物医学领检测技术领域，具体涉及与主动脉夹层疾病相关的SNP位点及其应用。

背景技术

主动脉夹层(Aortic Dissection，AD)系指主动脉在一系列外在因素作用下(如高血压、外伤等因素)，基于主动脉壁本身存在或不存在病变，主动脉内膜出现破口，血液自主动脉内膜破口侵入主动脉壁的中层，主动脉中层不断撕裂、纵向分离，导致呈现出一种主动脉管腔真假两腔并存的状态。主动脉夹层起病急骤，发展迅速，预后凶险。尽管目前的各种治疗手段不断改进，但是AD的发病率和死亡率依然居高不下，并且如果不进行恰当和及时的治疗，AD死亡率极高。有报道称约有近20％的AD患者在到达医院之前死亡，而30％的患者在住院治疗期间死亡。目前临床上对AD缺乏有效的药物干预手段，主要采用外科手术进行血管置换、介入手术进行主动脉腔内支架植入以及两者结合应用，但是手术风险大，费用高，并且远期的复发率和再手术率仍然较高。伴随着影像学技术及医疗设备条件的进步，AD疾病检出率呈上升趋势，对于AD的发生及发病机制的研究也在不断深入。AD生成机制与遗传学密切相关。AD相关易感基因是目前AD遗传学研究热点，从基因角度分析与AD相关的易感基因很有意义。

伴随着人类基因组测序计划完成和高通量测序技术的发展，医学领域的科研工作者对人类基因组有了更深入的认识，对疾病相关遗传变异、相关基因的功能及调控网络的认识也逐渐加深，高通量测序得到广泛应用。全外显子组测序(Whole Exome Sequencing，WES)技术是利用外显子捕获试剂盒捕捉整个外显子区域DNA，然后进行高通量测序，是目前基因组学研究中非常成熟的一部分。近两年来，采用WES技术研究复杂疾病的文献报道开始出现，展现出该技术的巨大优势。ZiganshinBA等学者对102例胸主动脉瘤和AD患者进行WES分析，发现21.6％患者存在1个或多个其它基因的突变， 3.9％患者存在FBN1、COL5A1、MYLK、FLNA基因的有害突变。Janice L. Farlow等对45位颅内动脉瘤患者进行了WES，发现68个基因存在变异， TMEM132B基因的表达在患者和对照者中显著不同(44VS 16)。表明在常见复杂疾病中可通过高通量测序鉴别罕见变异。并且，有越来越多的证据表明稀有变异对常见疾病易患性有中到强度的效应影响，可解释很大部分的遗传缺失。罕见突变(MAF<1％或5％)的局限性，而且费用、时间成本均较高。全基因组关联分析(genome wideassociation studies，GWAS)只能识别常见变异(common variants，CV：MAF＞5％)，忽略了某些未知的与AD发病相关的功能候选基因。WES相较于全基因组测序(Whole GenomeSequencing，WGS)、GWAS具有对基因组编码区域的高深度测序，除大片段的缺失重复及SNP外，还能发现低频突变、罕见突变，识别范围更广，具有低成本、高效率的优势，是寻找疾病易感基因或位点特别是发现新易感基因或位点最常用的方法之一，可弥补GWAS不能筛选未知候选基因的缺陷。

AD血管风险高，不容易早期发现，基于已知致病基因的遗传学诊断阳性率有限，难于做早期风险评估，亟待发现新的分子标志物。

发明内容

为了解决上述问题，本发明的目的在于提供一种基于全外显子测序技术筛选与主动脉夹层疾病相关的SNP位点，所述SNP位点包括以下主动脉夹层易感性SNP位点中的一种或两种：PTK6_c.205G>A_p.E69K；TSHZ2_c.724C>T_p.R242C。

进一步地，本发明提供了所述的SNP位点在制备检测主动脉夹层疾病的试剂中的应用。

进一步地，本发明还提供了所述的SNP位点在制备预测主动脉夹层疾病早期诊断试剂盒中的应用。

优选的，所述的试剂盒包括检测主动脉夹层易感性SNP位点 PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C是否发生突变试剂。

优选的，所述试剂盒可以为采用本领域已知的任何技术检测SNP的试剂，只要其能够检出样本中易感SNP位点PTK6_c.205G>A_p.E69K和/或 TSHZ2_c.724C>T_p.R242C是否存在突变。其包括但不限于下面列举的各实施方式。

在第一实施方式中，该试剂盒包括利用测序法检测样本中易感SNP位点 PTK6_c.205G>A_p.E69K位点A等位基因存在，和/或易感SNP位点 TSHZ2_c.724C>T_p.R242C位点T等位基因存在的试剂。测序法是本领域所公知的技术，所需的引物等试剂，本领域普通技术人员均可根据需要自行选择(参见ABI、Beckman等公司测序仪相关使用说明)，在此不再赘述。利用该试剂盒，通过测序法可以直接测得样本中易感SNP位点 PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C的序列，从而判断其是否携带相应位点等位基因的变异，进而判断其主动脉夹层的易感性。

在第二实施方式中，所述试剂盒包括利用Taqman探针SNP检测法检测样本中易感SNP位点PTK6_c.205G>A_p.E69K和/或 TSHZ2_c.724C>T_p.R242C基因型的试剂。其中所用Taqman探针为针对易感SNP位点PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C设计的探针，该探针可以使用试剂公司提供的探针；也可使用软件自行设计，如PREMIERBiosoft公司的BeaconDesigner 7.5。

在第三实施方式中，所述试剂盒为利用PCR-单链构象多态性法检测样本中易感SNP位点PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C 基因型的试剂盒。该试剂盒中包括用于扩增易感SNP位点 PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C的引物，PCR试剂、对照样本和检测构象的电泳所需试剂。所述电泳优选非变性聚丙烯酰胺凝胶电泳。对照样本中包括易感SNP位点PTK6_c.205G>A_p.E69K GG纯合子的阴性对照样本和位点AA纯合子的阳性对照样本中的至少一个，另外同时还可以包括易感SNP位点TSHZ2_c.724C>T_p.R242C CC纯合子的阴性对照样本，和位点TT纯合子的阳性对照样本中的至少一个，此外还可以包括也可以不包括杂合子相应杂合子的对照样本。优选同时包括上述三类对照样本。待测样本扩增产物与对照样本的扩增产物同时电泳，比较其电泳结果可以得出待测样本中是否携带相应等位基因变异的检测结果。

更进一步地，本发明提供了所述的SNP位点在制备诊断主动脉夹层疾病装置中的应用。

优选的，所述诊断装置是测序芯片。

再进一步地，本发明提供一种筛选与主动脉夹层疾病相关的SNP位点的方法，所述方法包括以下步骤：

(1)提取主动脉夹层患者和正常对照的外周血样本的DNA；

(2)将DNA进行超声片段化,打断片段末端补平，3’端加A，连接adaptor，选择350～400bp之间的片段制备全基因组文库；

(3)应用GenCap液相捕获目标基因技术进行全外显子检测；用高通量测序仪进行双端测序,读长为100bp；

(4)外显子测序后，进行常规过滤分析，以hg19、dbSNP(v147)为参考基因组版本为筛选过滤标准；留下致病位点、正常人频率在5％以下，及有文献报道的同义突变致病位点；为了更准确的筛选疾病相关位点，过滤掉SIFT、 Polyphen2、MutationTaster、GERP++预测全为良性的SNP；过滤掉 Mutcount>5，MutRatio>30％，正常数据库中MAF>1％的Indel；

(5)进行假阳性位点的排除工作，一是通过IGV软件对去冗余duplicates 后的bam文件(rmdup.sorted.bam)进行查看验证，突变情况不符合的视为假阳性位点；同时通过samtools软件对bam文件进行查看。然后，分别以突变的基因(a)、突变的位点(b)做为maker，进行稀有突变负荷(RVB)分析，统计病例组与对照组中该基因突变的样本数并计算OR值，OR>1,P<0.01。结果发现有多个突变位点与AD有统计学关联；

(6)通过基因功能和相关通路的注释分析以及文献检索数据对测序结果进行深度分析，最终发明人从48个突变位点筛选出两个主动脉夹层易感性 SNP位点：PTK6_c.205G>A_p.E69K和TSHZ2_c.724C>T_p.R242C。

本发明有益效果：

本发明通过发现了两个主动脉夹层易感SNP位点并研究了其在主动脉夹层早期预测的应用前景，本发明提供的主动脉夹层相关的基因位点可能是主动脉夹层的早期生物标志物，为进一步研究主动脉夹层的遗传分子机制，探索主动脉夹层早期防治的药物靶点提供了新的方向。

附图说明

图1是48个位点GO通路分析(MolecμLarFunction)选取通路中包括两个及以上基因，p value小于等于0.05的通路，对富集得分作图；

图2是48个位点GO通路分析(biologicalprocess)选取通路中包括两个及以上基因，p value小于等于0.05的通路，对富集得分作图。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。若未特别指明，实施例中所用的技术手段为本领域技术人员所熟知的常规手段。

本发明中术语解释如下：

GO(Gene Ontology):描述基因功能的综合性数据库，可分为分子功能(MolecμlarFunction,MF)，生物过程(biologicalprocess,BP)，和细胞组成(cellμl arcomponent,CC)三个部分。GO富集以p小于0.05为显著富集。

实施例1样品收集

自2017年1月至2018年4月在深圳市孙逸仙心血管医院通过主动脉 CTA确诊StanfordA型主动脉夹层的99例患者，收集全血2mL和590例正常对照样本来自迈基诺测序公司检索的数据库。已取得患者的知情同意，并通过伦理委员会审核。

样本处理：将EDTA抗凝全血按1:1比例与Trizol混合，充分混匀后置于1.8mL细胞冻存管中，在液氮中迅速冷却30s后置于-80℃的冰箱保存。

实施例2提取血液样品DNA

(1)将1mL通过EDTA(0.01M，China,Huamei Bioengineer)抗凝处理的血液中加入1mLCL细胞裂解液，轻柔地颠倒混匀6次，用3600rpm的转速离心5min，弃除上清液；

(2)像离心管中再次倒入1mLCL细胞裂解液，轻柔地颠倒混匀6次，用3600rpm的转速离心5min，弃除上清液；确保沉淀保留在管内的前提下，将离心管倒置在清洁的吸水纸上静置2min；

(3)配置蛋白酶K和缓冲液FG的混合液；

(4)加入500μL蛋白酶K和缓冲液FG的混合液，随即混匀至溶液无团块；

(5)65℃水浴30min，期间颠倒混匀数次；

(6)加入1mL异丙醇，随即颠倒混匀至出现簇状或丝状基因组DNA；

(7)用3600rpm的转速离心8min，弃除上清液；确保沉淀保留在管内的前提下，将离心管倒置在清洁的吸水纸上静置2min；

(8)加入1mL70％乙醇，振荡5sec，用3600rpm的转速离心3min，弃除上清液；

(9)重复步骤(8)；

(10)确保沉淀保留在管内的前提下，将离心管倒置在清洁的吸水纸上静置5min；

(11)常温状态下，空气干燥基因组DNA沉淀至所有液体完全挥发(至少5min)振荡5sec,用3600rpm的转速离心3min，弃除上清液；

(12)加入200μLTB缓冲液，低速涡旋振荡5sec，65℃水浴加热1h，期间轻弹助溶数次；

(13)使用NanoDrop ND8000(THERMO，USA)鉴定所提取的基因组 DNA浓度和纯度。浓度＞30ng/μL，1.8＜OD260/OD280＜2.0，且总量＞3μg 的基因组DNA样本视为合格，保存于-80℃的冰箱备用；

(14)在DNA初步定量基础上，进一步采用琼脂糖凝胶电泳检测(胶浓度：0.8％，电压：120V，时间：20min)精确定量DNA。电泳后出现弱带或主带之外存在杂带，说明DNA可能有降解现象或存在杂质，如果电泳后出现一条强带，则说明提取的DNA质量较高。全外显子测序要求DNA浓度≥50ng/μL，总量不低于5μg，OD260/280在1.8～2.0之间，OD260/230在2.0左右，样本无RNA污染，无降解或者轻微降解。

实施例3全外显子测序

1、DNA文库构建

将3μg DNA进行超声片段化,打断片段末端补平，3’端加A，连接 adaptor，选择350～400bp之间的片段制备全基因组文库。文库样本使用 Agilent2100生物分析仪(安捷伦科技公司，美国)进行质控。

2、目标区域捕获测序

应用GenCap液相捕获目标基因技术(北京迈基诺公司)进行全外显子检测。将1μgDNA文库与BL缓冲液和探针混合，95℃加热7min，65℃加热2min加入23μL预热至65℃的HY缓冲液，65℃杂交22h。用500μL 1X 结合缓冲液清洗50μL MyOne磁珠(美国Life Technology公司)3次，重悬于 80μL 1X结合缓冲液。加入64μL 2X结合缓冲液至杂交混合物中，转移至含有80μL MyOne磁珠的试管中。旋转混匀1h。用WB1缓冲液室温清洗磁珠 15min，用WB3缓冲液65℃清洗3次，每次15min。然后用洗脱缓冲液将结合的DNA洗脱下来。洗脱的DNA进行PCR反应，反应条件如下：98℃预变性30s；98℃变性25s，65℃退火30s，72℃延伸30s，共进行15个循环；最后72℃延伸5min。根据产品说明书用SPRIbeads(美国Beckman CoμLter) 纯化PCR产物。富集的文库用illumina HiSeq 2000测序仪进行双端测序,读长为100bp。

3、数据分析

测序得到的原始测序序列，里面含有带接头的、低质量的reads，会对后续信息分析造成干扰。为了保证信息分析质量，对原始数据用Cutadapt软件进行过滤，将去除接头序列和低质量碱基的测序数据用Burrows-Wheeler Aligner(BWA)比对到人类基因组(版本GRCh37/hg19)上，并用Picard工具去除冗余。目标区域的深度用coverageBed工具计算。通过GATK软件对 indel位点附近进行重比对从而提高序列的质量。GATK HaplotypeCaller检测SNP与indel突变信息，并用ANNOVAR对突变注释，注释数据库包括 dbSNP147,1000Genomes project,Exome sequencing project(ESP6500), Inhouse database(MyGenostics),gnomAD_genome_EAS.Online Mendelian Inheritance in Man(OMIM)andHuman Gene Mutation Database(HGMD Professional 2016.10)。致病性预测数据库包括SIFT,PolyPhen-2and MutationTaster，以及通过ACMG标准指南对致病性进行划分。

4、变异过滤

为了找出潜在的致病性突变，需要对数据进行过滤筛选。主要针对外显子区和可变剪切区域的突变，第一步：留下致病性分析(pathogenic_analysis) 中的致病性突变(pathogenic)位点。过滤筛选标准如下：

第二步：筛选突变碱基测序次数大于5，突变频率大于等于30％以及1000g2015apr、ESP6500si、Inhouse、ExAC_ALL、ExAC_EAS五个正常人突变数据库中保留没有出现或小于5％的位点；去除数据集中的同义突变位点，另外有文献报道的突变位点留下。

5、统计分析

对上述过滤后的突变结果进行稀有突变负荷分析，选择迈基诺590例数据作为对照样本，利用R语言中的Metabin分析OR值与95％置信区间，并用Mantel-Haenszel’s方法进行统计检验，p值小于0.01的认为显著。

6、结果

外显子测序后，进行常规过滤分析，筛选过滤标准：以hg19、dbSNP(v147) 为参考基因组版本。留下致病(pathogenic)位点、正常人频率在5％以下，及有文献报道的同义突变致病位点，得到59351个点突变。为了更准确的筛选疾病相关位点，进一步过滤掉SIFT、Polyphen2、MutationTaster、GERP++ 预测全为良性的SNP；过滤掉Mutcount>5，MutRatio>30％，正常数据库中 MAF>1％的Indel；经上述筛选得到29979个点突变。

之后，发明人进行了假阳性位点的排除工作，一是通过IGV软件对去冗余duplicates后的bam文件(rmdup.sorted.bam)进行查看验证，突变情况不符合的视为假阳性位点；同时通过samtools软件对bam文件进行查看。然后，分别以突变的基因(a)、突变的位点(b)做为maker，进行稀有突变负荷 (RVB)分析，统计病例组与对照组中该基因突变的样本数并计算OR值， OR>1,P<0.01。结果发现有48个突变位点与AD有统计学关联。

进一步，发明人对筛选出的48个突变位点进行进行基因功能和相关通路的注释。

48个突变位点所在基因的通路注释中选取count大于等于2的作图，如图1、2所示，在BP分析中，这些位点所在基因显著富集于氧化还原生物过程。

在48个突变位点的筛选结果中发现：

PTK6_c.205G>A_p.E69K(OR＝55.65，95％CI＝2.97～1041.89)与细胞迁移分化相关的位点突变。血管内皮细胞的增殖、迁移、分化可导致血管内皮形态变化，可能在AD的血管内膜形态变化相关。

TSHZ2_c.724C>T_p.R242C(OR＝24.8，95％CI＝2.74～224.27)与金属离子结合相关。金属离子与血管张力有密切联系，血管张力通过影响血压，可能与AD存在联系。

发明人结合基因功能和相关通路的注释分析以及文献检索数据对测序结果进行深度分析，最终发明人筛选出在病例组中发生突变的两个基因PTK6、 TSHZ2，所述两个基因携带的突变SNP位点如表1所示。 PTK6_c.205G>A_p.E69K和TSHZ2_c.724C>T_p.R242C突变位点可能是导致 AD的关键性分子，其可作为AD新型生物标志物，为指导其临床早期干预及靶向治疗提供了重要理论依据。

表1主动脉夹层SNP位点信息

注：Gene，基因名称；Chr，染色体；Pos，SNP位置；AA_change，氨基酸改变；ID，dbSNP注释ID。

实施例4用于预测主动脉夹层疾病早期诊断试剂盒的制作

组装本发明所述的用于主动脉夹层的试剂盒，具体有如下三种方案：

剂盒1：所述试剂盒包括检测PTK6_c.205G>A_p.E69K位点是否发生突变的试剂。

试剂盒2：所述试剂盒包括检测TSHZ2_c.724C>T_p.R242C位点是否发生突变的试剂。

试剂盒3：所述试剂盒包括检测PTK6_c.205G>A_p.E69K位点和 TSHZ2_c.724C>T_p.R242C位点是否发生突变的试剂。

所述试剂包括用于扩增PTK6_c.205G>A_p.E69K位点和/或 TSHZ2_c.724C>T_p.R242C位点的特异性引物，该引物可以采用软件来设计，如使用Primer5、Oligo6等；所述试剂还可以有相应PCR技术所需的常用试剂，如：dNTPs，MgCl₂，双蒸水，Taq酶等，这些常用试剂都是本领域技术人员熟知的，另外还可以有标准品和对照(如确定基因型的标准品和空白对照等)。此试剂盒的价值在于只需要外周血而不需要其他组织样品，通过最精简和特异的引物对检测SNP，再通过SNP谱辅助判断主动脉夹层，不仅稳定，检测方便，且精确，大大提高疾病诊断的敏感性和特异性，因此将此试剂盒投入实践，可以帮助指导诊断和更有效的个体化治疗。

虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (7)

1.一种基于全外显子测序技术筛选与主动脉夹层疾病相关的SNP位点，其特征在于，所述SNP位点包括以下主动脉夹层易感性SNP位点中的一种或两种：PTK6_c.205G>A_p.E69K；TSHZ2_c.724C>T_p.R242C。

2.权利要求1所述的SNP位点在制备检测主动脉夹层疾病的试剂中的应用。

3.权利要求1所述的SNP位点在制备预测主动脉夹层疾病早期诊断试剂盒中的应用。

4.如权利要求3所述的应用，其特征在于，所述的试剂盒包括检测主动脉夹层易感性SNP位点PTK6_c.205G>A_p.E69K和/或TSHZ2_c.724C>T_p.R242C是否发生突变试剂。

5.权利要求1所述的SNP位点在制备检测主动脉夹层疾病装置中的应用。

6.如权利要求5所述的应用，其特征在于，所述检测装置是测序芯片。

7.一种筛选与主动脉夹层疾病相关的SNP位点的筛选方法，其特征在于，所述方法包括以下步骤：

(1)提取主动脉夹层患者和正常对照的外周血样本的DNA；

(2)将DNA进行超声片段化,打断片段末端补平，3’端加A，连接adaptor，选择350～400bp之间的片段制备全基因组文库；

(3)应用GenCap液相捕获目标基因技术进行全外显子检测；用高通量测序仪进行双端测序,读长为100bp；

(4)外显子测序后，进行常规过滤分析，以hg19、dbSNP(v147)为参考基因组版本为筛选过滤标准；留下致病位点、正常人频率在5％以下，及有文献报道的同义突变致病位点；为了更准确的筛选疾病相关位点，过滤掉SIFT、Polyphen2、MutationTaster、GERP++预测全为良性的SNP；过滤掉Mutcount>5，MutRatio>30％，正常数据库中MAF>1％的Indel；

(5)进行假阳性位点的排除工作，一是通过IGV软件对去冗余duplicates后的bam文件(rmdup.sorted.bam)进行查看验证，突变情况不符合的视为假阳性位点；同时通过samtools软件对bam文件进行查看。然后，分别以突变的基因(a)、突变的位点(b)做为maker，进行稀有突变负荷(RVB)分析，统计病例组与对照组中该基因突变的样本数并计算OR值，OR>1,P<0.01。结果发现有多个突变位点与AD有统计学关联；

(6)通过基因功能和相关通路的注释分析以及文献检索数据对测序结果进行深度分析，最终发明人从48个突变位点筛选出两个主动脉夹层易感性SNP位点：PTK6_c.205G>A_p.E69K和TSHZ2_c.724C>T_p.R242C。