本发明属于分子生物学的应用领域。具体地说,本发明涉及mrnaadipor1、与其对应的分离的多核苷酸、与其特异性结合的寡核苷酸引物对、mrna芯片、试剂盒,还涉及用所述试剂盒定量检测血液中对应mrna的方法。
背景技术:
原发性肝癌是我国最常见的恶性肿瘤之一,其发病率、死亡率分别排名全国恶性肿瘤第4位和第2位。乙型、丙型肝炎病毒感染、黄曲霉素、饮酒、非酒精性脂肪肝肝、肥胖等因素是肝癌的危险因素。目前,原发性肝癌的治疗方法多样,其中手术治疗包括根治性外科切除、肝脏移植;非手术治疗,包括局部消融治疗、动脉化疗栓塞、基因分子靶向治疗、系统性化疗、放疗等。手术切除已被证实是使肝癌患者长期存活的最优治疗方案。然而,超过70%的原发性肝癌患者由于肿瘤的部位、大小、数量以及肝功能受损而无法接受肝脏切除手术。因此非手术疗法在肝癌治疗中的地位不言而喻。对于不可手术切除或无法接受外科手术的肝癌患者,nccn指南推荐外照射放疗作为治疗手段之一。随着计算机、放疗和影像技术的发展,放疗技术的飞速发展使精确放疗成为可能,体现在靶区足够精准、射线剂量集中释放于肿瘤靶区消灭肿瘤,其清除肿瘤的效果近似于“手术刀”般精准。从三维适形放射治疗开始,放射治疗越来越多地被用于原发性肝癌的治疗。目前,原发性肝癌的放射治疗包括一系列先进的技术,如调强放射治疗,体部立体定向放射治疗,粒子治疗等。目前的精准外照射技术可以确保肿瘤局部给予高剂量照射,同时保护剩余正常肝脏组织免受或仅受到低剂量照射。此外,外照射放疗适用于肝脏几乎所有位置的肿瘤。体部立体定向放射治疗(sbrt),也可作为消融/tace等治疗手段的替代方案,消融/tace等治疗失败后的选择方案,或者消融/tace禁忌证患者的治疗选择。
目前临床应用的超声检查、ct检查、核磁共振检查等手段能从不同角度对肝脏肿瘤的放疗疗效进行判断,然而,到目前为止,放疗疗效的评估主要依赖于局部控制率和转归,如完全治愈cr、部分治愈pr、稳定sd、进展pd等,总生存期(os),无病生存期(dfs),无进展生存期(pfs)等指标,目前尚没有基于个体的、实时评估放疗疗效的生物学指标。因而,寻找出若干具有临床价值的生物标志物,判定其用于评价的效能,是研究人员关注已久的问题。
肿瘤标志物(tumormarker,tm)是指在恶性肿瘤的发生和增殖过程中,因肿瘤细胞的相关基因表达或机体对肿瘤产生反应而异常变化的一类物质。其随着肿瘤的发生和发展进行变化,主要表现为体液中某种正常活性物质的异常增多或者异常物质的出现。因此,肿瘤标志物水平能在一定程度上反映出肿瘤细胞的增殖程度,并在肿瘤的诊断、治疗、预后判断中起到重要的作用。肿瘤标志物检查简便易行,对身体伤害小,仅需要少量血液或者其他体液即可做相关检测。
肿瘤是由环境因素和遗传因素相互作用所导致的一类疾病,从正常组织、癌前病变到恶性肿瘤,多种mrna参与了肿瘤的发生和发展。mrna的作用在多种肿瘤中已被证实,但是目前关于mrna的研究主要集中于肿瘤与正常组织的比较,而已有的文献报道显示组织中的mrna表达水平与血液中的mrna表达水平并不相同,两者并非都存在相关关系,这是因为组织中的mrna表达水平主要反映的是肿瘤局部组织的表达量,而血液中的mrna则是机体整体的表达情况。对于实施治疗后静脉血中mrna水平变化是否对癌症预后以及治疗方案的确定有指示作用仍有待研究,尤其是在sbrt治疗后肝癌患者血液中的mrna的变化尚未见报道。
技术实现要素:
本发明的目的:为了弥补现有技术中的不足,本发明的目的在于提供一种可用于立体定向放疗后迅速判断肝癌预后的分子标志物,并辅助确定后续治疗方案。
技术实现要素:
本发明实验研究发现肝癌患者sbrt治疗后adipor1基因表达升高,且升高比例不高于0.5838时,患者预后较好。据此可认为adipor1基因可作为肝癌sbrt治疗的预后标志物,可以开发判断肝癌sbrt预后的工具。
通过使用上述定量检测血液样本中mrnaadipor1的方法,可分别采集肝癌患者sbrt治疗前和治疗后的血样,检测adipor1的表达水平。筛选出放疗后较放疗前adipor1表达上调,且变化比例高于0.5838的样本,可视为sbrt治疗后肝癌预后差的对象,建议改变sbrt治疗策略,及时行进一步放疗疗程,或建议进行sbrt之外的其它治疗方案。若放疗后较放疗前adipor1升高比例低于0.5838,可认为sbrt疗效较好。
此方法方便简捷,通过对定量检测出的患病个体在sbrt治疗前后血液样本中上述mrna的表达量进行比较,利于在sbrt治疗后早期判断预后并予以积极预防和及时治疗,对提高肝癌的生存率,降低死亡率有重要意义。
本发明所使用的肝癌血液样本有以下纳入标准:
1.小肝癌(≤5cm)或大肝癌(5-10cm),肝功a级或b级;
2.未经任何治疗的首诊sbrt患者,80gy<bed<100gy;
3.无化疗等严重影响血液指标的其他综合治疗。
所有研究对象对本研究均知情并签署了知情同意书。
本发明包括下述技术方案:
第一方面的技术方案提供了一种mrna,所述的mrna具有如seqidno:1所示的序列。
另外,还提供了第一方面技术方案所述的mrna的用途,用于制备用于定量检测血液样本中具有如seqidno:1所示的序列的mrna的寡核苷酸引物对、mrna芯片或试剂盒。
第二方面的技术方案提供了一种分离的多核苷酸,所述多核苷酸能转录成第一方面技术方案所述的具有如seqidno:1所示的序列的mrna。
第三方面的技术方案提供了寡核苷酸引物对,所述的寡核苷酸引物对特异性地结合于第一方面技术方案所述的具有seqidno:1所示的序列的mrna。
优选地,所述的寡核苷酸引物对具有如seqidno:3-4所示的序列。所述seqidno:3-4分别对应为seqidno:1所示的mrna的上游和下游引物的序列。
第四方面的技术方案提供了一种mrna芯片,所述的mrna芯片包括:固相载体;以及固定在所述固相载体上的第三方面技术方案所述的寡核苷酸引物对。优选地,所述的寡核苷酸引物对具有如seqidno:3-4所示的序列。
另外,还提供了所述mrna芯片的用途,用于制备用于定量检测血液样本中具有如seqidno:1所示的序列的mrna的试剂盒。
第五方面的技术方案提供了一种用于定量检测血液样本中具有如seqidno:1所示的序列的mrna的试剂盒,所述的试剂盒中含有第三方面技术方案所述的寡核苷酸引物对或第四方面技术方案所述的mrna芯片。优选地,所述的寡核苷酸引物对具有如seqidno:3-4所示的序列。
第六方面的技术方案提供了一种定量检测血液样本中具有如seqidno:1所示的序列的mrna的方法,使用第三方面技术方案所述的寡核苷酸引物对、第四方面技术方案所述的mrna芯片、或第五方面技术方案所述的试剂盒。优选地,所述的寡核苷酸引物对具有如seqidno:3-4所示的序列。
在本发明中,筛选mrna的研究步骤简述为在得到肝癌样本原始测序数据并预处理后,筛选出肝癌sbrt治疗前后基因表达变化一致的mrna,结合治疗后3个月的随访数据,从中筛选sbrt治疗后差异表达且有治疗价值的mrna,再通过绘制roc曲线评价mrna作为肝癌患者sbrt预后标志物的能力。
另外,本发明针对上述筛选出的mrnaadipor1,设计了可与其特异性结合的具有seqidno:3-4的序列的寡核苷酸引物对,进而获得包括固相载体以及固定在所述固相载体上的寡核苷酸引物对的mrna芯片、含有上述寡核苷酸引物对或上述mrna芯片的试剂盒。
通过使用上述的寡核苷酸引物对、mrna芯片、或试剂盒,本发明通过对血液样本提取的rna进行逆转录pcr和pcr,提供了一种定量检测血液样本中adipor1基因表达水平改变的便捷方法。此方法可用于sbrt治疗后迅速检测肝癌预后以及后续治疗方案的辅助确定,具体参见下文介绍。
具体实施方式
以下是本发明的具体实施例,对本发明的技术方案做进一步的描述,但是本发明的保护范围并不限于这些实施例。凡是不背离本发明构思的改变或等同替代均包括在本发明的保护范围之内。本发明属于分子生物学领域,在实施例中,包括了很多常用试剂、常规方法,故实施例中出现的试剂不仅限于列出的试剂供应商品牌。下述实施例中所使用的实验方法如无特殊说明,均为本领域的常规方法,或按照制造厂商所建议的条件与实验步骤。
实施例1.筛选sbrt治疗前后具有共同变化特征的mrna
通常地,采集临床血液样本用于rna检测使用的是edta抗凝管。但由于rna的不稳定性,在采血后的几分钟内,血液中rna转录水平急剧降低,处理过程中的变异因素会导致试验结果不稳定。本研究中,我们使用了paxgenebloodrnatube(qiagen,no./id:762165)(简称brt管)采集临床血样,该采血管中添加了特殊试剂,可以快速地保护细胞内的rna不被降解,在18~25℃可以保存3天,2~8℃可以保存5天,-20/-80℃条件下保存8年。因此,可以为样本分析前处理提供准确的测试,并为实验使用提供有效性。此外,为保证实验结果的准确性和可靠性,我们使用了与之配套的试剂盒(paxgenebloodrnakit;catno./id:762174)进行后续实验研究,而非常规使用的trizol法提取血液白细胞中的总rna。
我们采集肝癌患者sbrt治疗前、治疗后1(出院前)和治疗后2(1.5-2个月)血样进行高通量测序,筛选出sbrt治疗后具有共同变化特征的mrna,筛选标准为qvalue(即校正后的pvalue)<0.05且|log2fc|>1。
一.高通量测序实验步骤:
1.样品采集和保存
1.1采血纳入标准:
(1)小肝癌(≤5cm)或大肝癌(5-10cm),肝功a级或b级;
(2)未经任何治疗的首诊sbrt患者,80gy<bed<100gy;
(3)无化疗等严重影响血液指标的其他综合治疗。
1.2采血流程
(1)使用前,将bd公司的paxgenebloodrnatube(货号:762165)正立放置于室温;
(2)采集血样,2管/人,2.5ml/管;
(3)血样采集后,立即将采血管上下颠倒混匀(十次),将采血管正立置于塑料试管架上,室温放置2h;
(4)室温放置2h后先于-20℃放置24h,然后转移至-80℃长期保存。
2.rna分离和质控(由novogene实验部完成)
2.1rna分离
实验步骤:(试剂盒参考:paxgenebloodrnakit;catno./id:762174)
(1)离心收集,将全血转入brt管中,3000-5000g,室温,离心10min(注意:确保全血在brt管中室温下孵育至少2h,使其充分裂解)
(2)去上清,加入4mlrnase-freewater到中心沉淀中,用新的bdhemogardclosure盖上盖子(注意:弃上清时不要破坏沉淀,并用干净纸干燥管壁)
(3)涡旋振荡直到沉淀完全溶解,3000-5000g,室温,离心10min,去上清(注意:上清清除不完全会抑制裂解并且会稀释裂解所得产物,因此会影响rna与膜结合。)
(4)加350μlresuspensionbuffer(br1),涡旋振荡直至沉淀消失
(5)将样品转移到1.5ml离心管(mct),加300ulbindingbuffer(br2)加40μlprotrninasek(pk)涡旋混匀5秒,在恒温振荡器中55度孵育10min,转数400-1400rpm,孵育完成后提前准备65度金属浴为第20步做准备(注意:br2和pk要分开加)
(6)将裂解产物转入paxgeneshredderspincolumn(psc)中,置于2mlprocessingtube(pt),12000g离心3min(注意:产物因全部转移,并所用转速不要高于20000g)
(7)取上清到新1.5ml离心管中,不要吸到沉淀
(8)加350μl无水乙醇,涡旋混匀,微型离心机离心1-2秒使盖子管壁上液体落到底部。(注意:离心不要超过1-2秒,不然会影响总产量)
(9)取700μl混合物到paxgenernaspincolumn(prc,红色)装到2mlpt中,12000g离心1min,将pcr放到新pt中,弃掉旧管。
(10)重复9,将剩下的混合物都通过prc
(11)取350μlwashbuffer1(br3)到prc中,12000g离心1min,将pcr放到新pt中,弃掉旧管。
(12)加10μl预先配置dnasei(rnfd)和70uldnadigestionbuffer(rdd)到1.5离心管中,轻轻混匀(注意:dnasei对物理损伤敏感,应轻轻混匀,避免涡旋振荡。此步骤可根据样品数做相应倍数混合母液)
(13)取上述混合液80μl到prc中,室温孵育15min(注意:上述混合液应该完全加入到prc管中心的膜上,避免加到管壁否则会影响dna的消化)
(14)加350μlwashbuffer1(br3)到prc中,12000g离心1min,将pcr放到新pt中,弃掉旧管。
(15)加500μlwashbuffer2(br4)到prc中,12000g离心1min,将pcr放到新pt中,弃掉旧管。(注意:br4使用前确保乙醇已经加入)
(16)加500μlwashbuffer2(br4)到prc中,12000g离心3min
(17)弃掉有液体的pt旧管,将pcr放到新pt中12000g离心1min
(18)弃掉有液体的pt旧管,将pcr放到新1.5ml离心管mct中,加入40μlelutionbuffer(br5)到prc管中心膜上,12000g离心1min洗脱rna(注意:此步骤很关键,影响整个rna产量,因此要确保buffer完全加到膜上)
(19)重复步骤18中的洗脱步骤,即再加入40ul的br5离心到同一离心管中
(20)将洗脱下的产物在65℃下孵育5min后立即放冰上。(注意:此步骤将rna变性为下游实验做准备,时间温度不要随意改变)
(21)如果rna产物不立即使用须放置于-80℃保存。
2.2rna定量和质控
(1)在1%琼脂糖凝胶上监测rna降解和污染。
(2)使用
(3)使用qubit
(4)使用bioanalyzer2100系统的rnanano6000assaykit(agilenttechnologies,ca,usa)评估rna完整性。
3.文库制备和测序
3.1测序的文库制备
每个样品的总量为3μgrna用作rna样品制备的输入材料。首先,通过epicenterribo-zerotmrrna去除试剂盒(epicentre,usa)除去核糖体rna,并通过乙醇沉淀清除无rrna残基。随后,按照制造商的推荐,使用
3.2聚类和测序
根据制造商的说明,使用truseqpeclusterkitv3-cbot-hs(illumia)在cbotclustergenerationsystem上进行索引编码样本的聚类。在簇生成后,将文库在illuminahiseq4000平台上测序,并产生150bp的配对末端读数。
二.高通量测序数据分析(由novogene基因调控部完成)
1.质量控制
fastq格式的原始数据(原始读取)首先通过内部perl脚本处理。在此步骤中,通过删除包含适配器的读取,读取包含ploy-n和从原始数据读取的低质量读取来获得干净数据(干净读取)。同时,计算了清洁数据的q20,q30和gc含量。所有下游分析均基于高质量的清洁数据。
2.比对到参考基因组
参考基因组和基因模型注释文件直接从基因组网站下载。使用bowtie2v2.2.8构建参考基因组的索引,并使用hisat2(langmead,b,etal.)v2.0.4将配对末端清洁读数与参考基因组比对。hisat2使用′--rna-strandnessrf′运行,其他参数设置为默认值。
3.转录组装配
每个样品的映射读数由stringtie(v1.3.1)(mihaelapertea,etal.2016)在基于参考的方法中组装。stringtie使用新的网络流算法以及可选的从头组装步骤来组装和定量代表每个基因座的多个剪接变体的全长转录物。
4.编码潜力分析
4.1cnci(编码-非编码-索引)(v2)描述邻接的核苷酸三联体,以有效区分蛋白质编码和非编码序列,而不依赖于已知的注释(sun,etal.2013)。我们使用cnci和默认参数。
4.2cpc(编码潜力计算器)(0.9-r2)主要通过评估转录物中orf的范围和质量,并用已知的蛋白质序列数据库搜索序列以阐明编码和非编码转录物(kong,etal.2007)。我们使用了ncbi真核生物的蛋白质数据库,并在我们的分析中设置了电子值′1e-10′。
4.3pfam-sca。在所有三个可能的框架中翻译了每个转录本并使用pfamscan(v1.3)来识别pfam数据库中记录的任何已知蛋白质家族域的出现(第27版;使用pfama和pfamb)(punta,etal.2012)。任何具有pfam命中的抄本将在以下步骤中排除。pfam搜索使用-e0.001--dome0.001的默认参数(bateman,etal.2002)。
4.4phylocsf(系统发育密码子替换频率)(v20121028)检查了保守编码区比对特征的进化特征,例如同义密码子取代和保守氨基酸取代的高频率,以及其他错义和无意义取代的低频率来区分蛋白质编码和非编码转录本(lin,etal.2011)。构建多物种基因组序列比对,并使用默认参数运行phylocsf。
4.5phast(v1.3)是一个包含大量统计程序的软件包,大多数用于系统发育分析(siepel,etal.2005),而phastcons是保守元素的保护评分和识别程序。我们使用phylofit计算物种间保守和非保守区域的系统发育模型,然后将模型和hmm转换参数提供给phylop,以计算编码基因的一组保守分数。
4.6基因表达水平的定量。使用cuffdiff(v2.1.1)计算每个样品中编码基因的fpkm(trapnell,c.etal.2010)。通过对每个基因组中转录物的fpkm求和来计算基因fpkm。fpkm表示基于片段长度计算的每百万个外显子每百万个片段的片段,并且读取计数比对到该片段。
4.7差异表达分析。ballgown套件包括用于交互式探索转录组装配的功能,转录物结构的可视化和每个基因座的特征特异性丰度,以及对注释特征的组装特征的事后注释(alyssac.frazee,etal.2014)。具有p-adjust<0.05的转录物被指定为差异表达的。cuffdiff提供了统计程序,用于使用基于负二项分布的模型确定数字转录物或基因表达数据中的差异表达(trapnell,c.etal.2010)。具有p-adjust<0.05的转录物被指定为差异表达的。
5.高通量测序差异基因筛选结果
筛选qvalue(即校正后的pvalue)<0.05且sbrt治疗后|log2fc|>1的差异表达的基因。共得到28个放疗后1(出院前)升高/降低,并在放疗后2(1.5-2个月)仍保持升高/降低的mrna。其中,保持升高的mrna共16个,保持降低的mrna共12个。见表1。
表1.肝癌sbrt治疗后1(出院前)和治疗后2(1.5-2个月)较治疗前保持升高/降低的mrna
实施例2.对sbrt治疗后血液样本的实时pcr定量检测
对接受sbrt治疗的肝癌患者(50%为大肝癌,50%为小肝癌),取其sbrt治疗前、治疗后1(出院前)和治疗后2(1.5-2个月)的血样,采血纳入标准和采血流程同实施例一,对血样进行实时pcr定量检测。
实验步骤具体为:
一.全血rna的提取
因所有血样的采集均使用paxgenebloodrnatube(货号:762165),为保证实验结果的准确性和可靠性,我们使用了与之配套的试剂盒(paxgenebloodrnakit;catno./id:762174)进行全血rna的提取,而非常规使用的trizol法提取血液白细胞中的总rna。
方法同实施例1。
二.逆转录反应:
1.按下列组份配制rt反应液(反应液配制请在冰上进行)。
2.轻柔混匀后进行反转录反应,条件如下:
37℃15min(反转录反应)
85℃5sec(反转录酶的失活反应)
4℃(保存)
三.高通量荧光定量qpcr检测
1.上机使用的模板稀释倍数:10倍
2.检测仪器:steponeplustm实时荧光定量pcr
3.qpcr检测试剂盒参考:tbgreentmpremixextaqtmii(tlirnasehplus)(takara,codeno.rr820a)
4.引物设计:根据genbank中adipor1基因和gapdh基因的编码序列设计qpcr扩增引物,由上海启因生物科技有限公司合成。具体引物序列如下:
adipor1基因:
正向引物为5’-acttgttgaaccagaatggtctc-3’(seqidno.3);
反向引物为5’-tccacttctctacctgcttgtc-3’(seqidno.4),
gapdh基因:
正向引物为5’-caatgaccccttcattgacc-3’(seqidno.5);
反向引物为5’-gacaagcttcccgttctcag-3’(seqidno.6)。
5.pcr反应体系:
pcr反应条件:
95℃30s1cycle
95℃5s40cycles
60℃30s40cycles
ct检测限:40个循环
6.结果计算:采用比较阈值法计算结果,即目的基因的量=2-δδct,在该公式中,ct指的是热循环仪检测到反应体系中荧光信号的强度值,首先将所有基因ct值整理好,用每一组样本自身的目的基因ct值减去自身内参基因ct值,得到的数就是δct,然后将每一组样本每一个目的基因的δct减去对照组样本的δct,并同时对所有结果取相反数(就是加一个负运算)该步运算得到的结果就是-δδct,即目的基因的表达量。
本研究qpcr技术创新点在于我们未采用常规的“三步法”,而使用了“两步法”进行pcr反应:95℃预变性,酶激活30s;95℃变性,5s;60℃退火延伸,30s,变性退火延伸共40个循环。这样既能节省pcr反应时间,又能降低引物错配率,提高引物特异性。
实施例3.结合临床病历,评价预后价值
1.收集临床采血的sbrt治疗的肝癌患者资料,随访3个月,复查经sbrt治疗后3个月的肝脏病灶大小,判断肝癌三个月转归情况(cr,pr,sd或pd),将cr与pr归为sbrt有效组,sd与pd归为sbrt无效组。
2.结合qpcr定量结果,将疗效分组与adipor1变化差值倍数做关联性分析,观察sbrt治疗后1(出院前)vs治疗前和治疗后2(1.5-2个月)vs出院前,adipor1的比例变化在有效组和无效组中是否有差异。计算公式为:治疗后vs治疗前的比例变化=(sbrt治疗后1-治疗前)/治疗前或(sbrt治疗后2-治疗前)/治疗前。(结果见表2)。
表2.adipor1在sbrt治疗后两个不同阶段的有效组和无效组中的差值变化情况
结果显示adipor1放疗后1(出院前)vs放疗前的差值变化比在有效组和无效组中存在显著差异(p<0.05),而adipor1在放疗后2(1.5-2个月)vs放疗前的变化比例在有效组和无效组中差异无统计学意义(p>0.05)。表明检测sbrt治疗后1(出院前)的血样和sbrt治疗前的血样中adipor1的变化比例,对判断肝癌sbrt治疗后3个月的转归情况,有一定的指导意义。
依据疗效分组和肝癌患者sbrt治疗前、治疗后1(出院前)的qpcr检测结果,使用medcalc软件绘制roc曲线,根据灵敏度、特异度、约登指数和roc曲线下面积auc值评价adipor1作为肝癌sbrt预后生物标志物的能力。roc分析结果显示:sbrt治疗后1(出院前)vs治疗前的变化比例在有效组无效组的区分中,adipor1的曲线下面积auc为0.896,灵敏度100%,特异度为83.33%,约登指数0.83。(结果见表3)。
上述结果提示肝癌患者sbrt治疗后1(出院前)血样中的mrnaadipor1作为肝癌患者sbrt三个月转归情况预后标志物的能力较强,提示,adipor1可用作sbrt治疗后出院前血液检测用的mrna以预测肝癌sbrt治疗后转归情况。
表3adipor1在肝癌sbrt治疗后1vs治疗前的roc分析
结合表2和表3可知,adipor1在肝癌患者sbrt治疗完成后出院前较治疗前表达上调,且在无效组中上调比例高于有效组,使用adipor1,区分sbrt治疗后有效组无效组的最佳阈值是0.5838,其意义在于,患者sbrt治疗后出院前血样检查若提示adipor1升高,且升高比例高于0.5838,则sbrt治疗后三个月转归情况较差;若升高比例低于0.5838,则疗效较好。
另外,患者出院前医院会例行抽血做出院前检查,监测sbrt治疗后出院前这个时间点上adipor1的变化情况可及时有利地帮助确定是否患者预后较好,并且辅助判断患者是否需要进一步放疗疗程,为临床总体治疗方案的确定提供可靠的辅助参考。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。