阅读说明：本技术 一种基于检测LncRNA的癌症预后评估试剂盒 (Cancer prognosis evaluation kit based on LncRNA detection ) 是由 叶定伟 于 2018-07-12 设计创作，主要内容包括：本发明公开了一种基于检测LncRNA的癌症预后评估试剂盒,试剂盒包括可检测RP11-783K16.13、RP11-727F15.11、PRKAG2-AS1、AC013460.1、CRNDE表达水平的检测试剂；本发明提供了可高效预测癌症的预后的试剂盒。(The invention discloses a cancer prognosis evaluation kit based on LncRNA detection, which comprises a detection reagent capable of detecting expression levels of RP11-783K16.13, RP11-727F15.11, PRKAG2-AS1, AC013460.1 and CRNDE; the present invention provides a kit that can efficiently predict the prognosis of cancer.)

一种基于检测LncRNA的癌症预后评估试剂盒

技术领域

本发明涉及医药卫生领域，尤其涉及一种基于检测LncRNA的癌症预后评估试剂盒。

背景技术

***癌是指发生在***的上皮性恶性肿瘤。2012年我国肿瘤登记地区***癌发病率为9.92/10万，列男性恶性肿瘤发病率的第6位。发病年龄在55岁前处于较低水平，55岁后逐渐升高，发病率随着年龄的增长而增长，高峰年龄是70～80岁。家族遗传型***癌患者发病年龄稍早，年龄≤55岁的患者占43％。

***癌根治术是治疗局限性***癌的重要手段之一。但是，***癌根治术后，仍有相当多的患者会出现生化复发。当患者在发生生化复发后，需要明确是否已经发生临床复发，以及是局部复发，区域***转移还是远处转移。在临床工作中，***癌根治术后生化复发的评估常常是令人头疼的，从而导致其治疗也是缺乏精准的。

长链非编码RNA(Long non-coding RNA,LncRNA)是长度大于200个核苷酸的非编码RNA。研究表明,LncRNA在剂量补偿效应、表观遗传调控、细胞周期调控和细胞分化调控等众多生命活动中发挥重要作用，成为遗传学研究热点。研究发现,长链非编码RNA的表达或功能异常与人类疾病的发生密切相关,其中就包括癌症、退行性神经疾病在内的多种严重危害人类健康的重大疾病,具体表现为长链非编码RNA在序列和空间结构上的异常、表达水平的异常、与结合蛋白相互作用的异常等。

现有研究已经发现了大量的与***癌相关基因，2018年发表于癌症研究和临床肿瘤学杂志(Cancer Research and Clinical Oncology)的一篇论文，Validation of a10-gene molecular signature for predicting biochemical recurrence andclinical metastasis in localized prostate cancer，公开了FRZB,LEF1,SDCBP,WNT2,ING3,ANK3,MEIS2,ANXA4,PLA2G7和CHD5共10个基因的表达用于预测***癌生化复发的模型。该论文还公开了采用该模型预测***癌的结果：AUC值为0.65，HR值为0.24，95％CI：0.09-0.59。

虽然现有技术公开了大量的与***癌相关基因，但是***癌相关基因在实际应用中用于评估***癌却受限。现有模型的预测效果仍有待增强，需要从大量的***癌基因中建立能够更加精准预测***癌复发风险的模型。

发明内容

针对上述技术问题，本发明提供了一种基于检测LncRNA的癌症预后评估试剂盒。

本发明的一个方面，公开了一种基于检测LncRNA的癌症预后评估试剂盒，所述试剂盒包括可检测RP11-783K16.13、RP11-727F15.11、PRKAG2-AS1、AC013460.1、CRNDE表达水平的检测试剂。

优选地，所述检测试剂为多核苷酸引物或探针。

优选地，多核苷酸引物的序列如SEQ ID NO:1-10所示。

优选地，其中所述癌症选自***癌，优选的，所述试剂盒用于***癌根治术后预后的评估。

优选地，本发明的癌症预后方法包括以下步骤：

a)检测样本中RP11-783K16.13、RP11-727F15.11、PRKAG2-AS1、AC013460.1、CRNDE的mRNA表达水平；

b)根据步骤a)获得的表达水平数据评估患者的癌症复发风险。

优选地，步骤b)中通过以下方法评估患者的癌症复发风险：

风险值＝(0.0766×RP11-783K16.13表达水平)+(0.2443×RP11-727F15.11表达水平)+(0.0042×PRKAG2-AS1表达水平)+(2.8117×AC013460.1表达水平)+(0.0162×CRNDE表达水平)，其中所述的表达水平为步骤a)中检测获得的mRNA表达值。

优选的，所述检测mRNA表达水平的方法包括：Affymetrix/Illumina的芯片检测、全转录组鸟枪法测序、RT-PCR。

本发明的另一个方面，检测RP11-783K16.13、RP11-727F15.11、PRKAG2-AS1、AC013460.1、CRNDE表达水平的检测试剂在制备用于癌症预后的产品中的用途。

优选地，本发明可用于***癌生化复发的早期预测，优选的，本发明可用于***癌根治术后预后的评估。

优选地，所述癌症预后方法包括以下步骤：

a)检测样本中RP11-783K16.13、RP11-727F15.11、PRKAG2-AS1、AC013460.1、CRNDE的mRNA表达水平；

b)根据步骤a)获得的表达水平数据评估患者的癌症复发风险。

优选地，步骤b)中通过以下方法评估患者的癌症复发风险：

风险值＝(0.0766×RP11-783K16.13表达水平)+(0.2443×RP11-727F15.11表达水平)+(0.0042×PRKAG2-AS1表达水平)+(2.8117×AC013460.1表达水平)+(0.0162×CRNDE表达水平)，其中所述的表达水平为步骤a)中检测获得的mRNA表达值。

优选的，所述检测mRNA表达量的方法包括：Affymetrix/Illumina的芯片检测、全转录组鸟枪法测序、RT-PCR。

与现有技术相比，本发明的技术方案具有以下优点：本发明的试剂盒和方法能够准确预测癌症预后后果，尤其是对于***癌生化复发的诊断具有较高的预测价值和重要的基础及临床价值和广阔的应用前景。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例的5个LncRNA为基础的基因决策模型；

图2是本发明实施例的343例患者ROC曲线图；

图3是本发明实施例的343例患者生化复发曲线图；

图4是本发明实施例的114例患者ROC曲线图；

图5是本发明实施例的114例患者生化复发曲线图。

具体实施方式

以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。

术语“癌症”和“癌”是指或者描述哺乳动物中通常以异常或失控的细胞生长为特征的生理状态。癌症和癌症病理可以伴随着例如转移、干扰邻近细胞的正常机能、以异常水平释放细胞因子或其他分泌产物、抑制或加重炎性或免疫应答、瘤形成、癌前病变、恶性肿瘤、周围或远距离组织或器官如***的浸润等。尤其包括的是***癌。

术语“预后”指对医学结果(medical outcome)，例如不良或良好结果的预测(例如长期存活的可能性)；消极的预后或不良结果包括复发、疾病进展(例如肿瘤生长或转移或耐药性)或死亡的预测。积极的预后或良好结果包括疾病好转(例如无疾状态)、改善(例如肿瘤消退)或稳定的预测。

本发明所述的RP11-783K16.13同The Cancer Genome Atlas(TCGA)数据库中的RNA-seq数据中的RP11-783K16.13。

本发明所述的RP11-727F15.11同The Cancer Genome Atlas(TCGA)数据库中的RNA-seq数据中的RP11-727F15.11。

本发明所述的PRKAG2-AS1同The Cancer Genome Atlas(TCGA)数据库中的RNA-seq数据中的PRKAG2-AS1。

本发明所述的AC013460.1同The Cancer Genome Atlas(TCGA)数据库中的RNA-seq数据中的AC013460.1。

实施例1***癌复发和/或存活后果或预后的评估模型的建立

通过从The Cancer Genome Atlas(TCGA)数据库检索并通过倾向性匹配分析，筛选了343例***癌患者作为发现集，从其中提取了相应的RNA-seq数据，应用Lasso Cox回归模型分析预测这些患者生化复发密切相关的LncRNA。最终确定和构建了以5个LncRNA为基础的基因决策模型(见图1)。并根据决策模型的风险，采用以下评分公式，得出风险值，将患者分为***癌术后生化复发高危组和低危组。

风险值＝(0.0766×RP11-783K16.13表达水平)+(0.2443×RP11-727F15.11表达水平)+(0.0042×PRKAG2-AS1表达水平)+(2.8117×AC013460.1表达水平)+(0.0162×CRNDE表达水平)。

决策LncRNA表达值代入模型中计算出风险值，风险值大于0.55判定为***癌生化复发高危组，风险值小于0.55判定为***癌生化复发低危组。

将343例***癌患者的风险值分别计算出，结合患者的临床资料，绘制ROC曲线，见图2，可见两年内，***癌患者的早期生化复发诊断AUC值为0.722,推衍至五年生化复发的AUC值为0.704。并根据risk score分为的***癌术后生化复发高危组和低危组，绘制生化复发曲线，见图3，可见两组术后生化复发存在着显著差异，HR＝0.44,95％CI：0.27-0.72，C-index值为0.63。

实施例2用于癌症预后评估的试剂盒

本实施例中用于癌症预后的试剂盒包括：

总RNA提取试剂Trizol；

氯仿(三氯甲烷)；

异戊醇；

无水乙醇；

DEPC水(DD1005)；

焦磷酸乙二酯(DEPC)；

抗RNA酶液(RNaseZap)；

反转录试剂盒；

iQ SYBR Green Supermix；

序列为SEQ ID NO:1-10的多核苷酸引物。

利用试剂盒对各LncRNA进行实时定量PCR(qRT-PCR)检测后，初始数据结果以Ct值(cycle threshold)表示，即每个反应体系内的荧光信号达到设定的阈值所需要的循环数。每一样品的Ct值与该样品起始拷贝数的对数存在线性关系，起始拷贝数越多，则Ct值越小。采用ΔΔCT法计算各LncRNA表达水平，并进行标准化处理。

根据上述步骤获得的各LncRNA表达水平，代入以下***癌生化复发预测模型计算***癌生化复发的风险值，风险值大于0.55判定为***癌生化复发高危组，风险值小于0.55判定为***癌生化复发低危组。

风险值＝(0.0766×RP11-783K16.13表达水平)+(0.2443×RP11-727F15.11表达水平)+(0.0042×PRKAG2-AS1表达水平)+(2.8117×AC013460.1表达水平)+(0.0162×CRNDE表达水平)。

实施例3样本RNA的制备和预处理

(1)组织RNA提取实验用枪头、镊子在DEPC水中浸泡过夜后高温高压灭菌，预冷4℃离心机，实验台、移液枪、手套等用抗RNA酶液(RNase Zap)擦拭。

(2)组织匀浆取50～100mg组织样品，用无菌手术刀片稍切碎后置于1.5mL EP管中，加入500μL Trizol试剂，用电动组织研磨器充分匀浆，然后补充加入500μL Trizol试剂。

(3)每1mL Trizol试剂匀浆样品中加入0.2mL氯仿，盖紧EP管盖。剧烈振荡15秒，室温放置3min。4℃12,000rpm离心15min(提前预冷离心机)。离心后混合体系将分为上层的无色水相，中层的蛋白质和下层的红色酚氯仿相。DNA溶于氯仿分布于下层，RNA溶于水相分布于上层。水相的体积为匀浆时加入Trizol试剂体积的60％左右。

(4)RNA沉淀，将上层水相转移至新的EP管中。每1mL Trizol试剂匀浆样品中加入0.5mL异丙醇。混匀后室温放置10min，4℃下12,000rpm离心10min。离心后将在管壁底部和侧壁上可见白色沉淀物。

(5)RNA清洗，轻轻弃去上清，向每1mL Trizol试剂匀浆样品中加入1mL 75％的乙醇，清洗RNA沉淀。振荡，4℃下10,000rpm离心5min。

(6)RNA溶解，轻弃乙醇溶液，空气中干燥RNA沉淀约5～10min。注意切勿完全干燥RNA沉淀，否则将大大降低RNA的可溶性。溶解RNA时，加入适量无RNase的DEPC水，用移液枪反复吹打，确保RNA充分溶解后，将RNA溶液保存于-80℃冰箱。

(7)RNA浓度测定和质控，使用

ND-2000紫外分光光度计测定RNA溶液的浓度和纯度。

1)测量前先用溶解RNA用的DEPC水调零；

2)用移液枪吸取2μL RNA样品滴加至测量基座的表面；

3)轻轻合上基座后，液滴会自动在上下基座之间形成液柱，完成测定后电脑中即显示RNA溶液的各种参数包括RNA浓度和纯度等。A260/A280的比值是一种常用的评估RNA纯度的参数，一般认为其比值范围1.8～2.1表示RNA纯度较好。

4)一次检测完成后，用擦镜纸轻轻擦去上下基座表面液体，即可进行下一个样品的检测

ii)琼脂糖凝胶电泳

1)制胶：称取1g琼脂糖加入100mL 1×TAE buffer中，置于微波炉加热至沸腾，灌制凝胶板，待胶凝后取下梳子，将凝胶板放入电泳槽内，加入适量的1×TAE buffer至液面完全覆盖胶面。

2)准备RNA样品：取3μg RNA，加3倍体积的甲醛上样染液，再加EB于甲醛上样染液中至EB终浓度为10ug/mL，将体系加热至70℃孵育5min使样品变性。

3)电泳：上样后，5～6V/cm电压下电泳，至溴酚蓝指示剂进入胶内至少2～3cm。

4)紫外透射光下观察结果：经过变性RNA电泳后，在凝胶成像系统上可见28SrRNA、18S rRNA和5S rRNA三条带。观察到28S rRNA条带的强度约为18S rRNA的2倍，而且5SrRNA条带较弱，说明总RNA未发生明显降解。

实施例4决策LncRNA检测

(1)采用全式金公司的1st Strand cDNA Synthesis SuperMix进行第一链cDNA的合成，反应体系的配制在冰上进行，反应体系如下：

(2)采用Bio-Rad公司的iQ SYBR Green Supermix进行实时定量PCR反应。一般引物终浓度为0.2μM可以得到较好的结果，反应性能不佳时，可在0.1～1.0μM范围内调整引物浓度。

(3)反应体系

PCR引物分别为：RP11-783K16.13的扩增引物为：SEQ ID NO:1、SEQ ID NO:2；

RP11-727F15.11的扩增引物为：SEQ ID NO:3、SEQ ID NO:4；

PRKAG2-AS1的扩增引物为：SEQ ID NO:5、SEQ ID NO:6；

AC013460.1的扩增引物为：SEQ ID NO:7、SEQ ID NO:8；

CRNDE的扩增引物为：SEQ ID NO:9、SEQ ID NO:10。

(4)反应程序

(5)对各LncRNA的表达进行定量

实时定量PCR反应结束后对Real time PCR扩增曲线及熔解曲线进行数据分析处理，初始数据结果以Ct值(cycle threshold)表示，即每个反应体系内的荧光信号达到设定的阈值所需要的循环数。每一样品的Ct值与该样品起始拷贝数的对数存在线性关系，起始拷贝数越多，则Ct值越小。采用ΔΔCT法计算各LncRNA表达水平。

实施例5***癌生化复发模型与***癌生化复发的相关性分析

采用本发明的试剂盒进行预测，对114例验证集的***癌患进行了验证，通过上述实施例的方法，将验证人群风险值计算出并分为高危组和低危组，绘制ROC曲线及生化复发曲线，见图4和图5，结果显示在验证人群中本发明的方法对***癌术后患者生化复发有很强的预测价值，两年及五年的生化复发ACU值分别为0.680和0.702，HR＝0.22,95％CI：0.09-0.56，C-index值为0.65。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

序列表

<110> 复旦大学附属肿瘤医院

<120> 一种基于检测LncRNA的癌症预后评估试剂盒

<141> 2018-06-20

<160> 10

<170> SIPOSequenceListing 1.0

<210> 1

<211> 21

<212> DNA

<213> Artificial Sequence

<400> 1

cctctgccca gcagacagaa c 21

<210> 2

<211> 19

<212> DNA

<213> Artificial Sequence

<400> 2

accccgcttt gggcagttt 19

<210> 3

<211> 21

<212> DNA

<213> Artificial Sequence

<400> 3

agagaagacc agcaccgagt a 21

<210> 4

<211> 20

<212> DNA

<213> Artificial Sequence

<400> 4

cagtgaaaac ccctcacccc 20

<210> 5

<211> 21

<212> DNA

<213> Artificial Sequence

<400> 5

agaaaagggg agggggtagt a 21

<210> 6

<211> 21

<212> DNA

<213> Artificial Sequence

<400> 6

tgttttatgc accctggcaa a 21

<210> 7

<211> 20

<212> DNA

<213> Artificial Sequence

<400> 7

gggcatctgt gatctggcat 20

<210> 8

<211> 20

<212> DNA

<213> Artificial Sequence

<400> 8

gtcaaccaca caggccatct 20

<210> 9

<211> 18

<212> DNA

<213> Artificial Sequence

<400> 9

cacgcctgtt cccttcgt 18

<210> 10

<211> 20

<212> DNA

<213> Artificial Sequence

<400> 10

ttcaaagacc aacggctgaa 20