m6A相关lncRNA在预测小细胞肺癌预后和化疗反应的应用的制作方法

m6a相关lncrna在预测小细胞肺癌预后和化疗反应的应用
技术领域
1.本发明涉及生物医学领域，具体涉及m6a相关lncrna在预测小细胞肺癌预后和化疗反应的应用。


背景技术：

2.小细胞肺癌(small cell lung cancer,sclc)是高致死性的高级别神经内分泌肿瘤，其特点在于肿瘤倍增时间短、生长迅速和早期转移扩散。sclc约占肺癌的15％，患者五年生存率不足7％。尽管分子靶向药物、免疫检查点抑制剂等新的治疗措施不断发展，小细胞肺癌患者的治疗策略近几十年仍未有明显的突破，化疗仍然是小细胞肺癌患者不可替代的一线治疗方案。然而，大多数小细胞肺癌患者即使对化疗敏感，但治疗后快速产生耐药。且患者病情进展迅速，极易发生转移，治疗手段十分有限。因此，临床亟需精准筛选预测化疗获益及预后的标志物，以便提高患者治疗效果，提高患者的预后。
3.n6-甲基腺苷(n6-methyladenosine，m6a)是真核生物rna中最丰富和最普遍的rna修饰，是癌症生物学的重要组成部分。m6a相关的生物学过程是动态、多层面、可逆的过程，主要由甲基化酶、甲基转移酶和结合蛋白介导功能的发挥。该修饰方式可以调控多种rna相关的生物学进程，包括rna降解、稳定、翻译、剪切和运输，最终调节靶基因的表达。
4.长链非编码rna(long-coding rnas，lncrnas)是核苷酸长度长于200的rna亚群。lncrnas参与m6a调节的肿瘤的发生和发展。m6a和lncrna都参与多项生物学进程，具有预测患者预后和治疗靶点的潜力，而二者之间相关性的研究少之又少。综上所述，鉴于小细胞肺癌的高度恶性、有限的治疗措施和较差的预后，识别并建立小细胞肺癌的预后预测和化疗获益预测标志物意义重大。


技术实现要素：

5.本发明所要解决的技术问题是如何预测小细胞肺癌预后和/或化疗反应和/或确定小细胞肺癌的预后预测和/或化疗反应预测的标志物。
6.为了解决上述技术问题，本发明首先提供了检测wars2-it1、ac005162、ac130352、ac015971、foxp1-it1、ac040963和linc00355这七种lncrna的表达量的物质在制备预测小细胞肺癌患者预后的系统和/或预测小细胞肺癌患者化疗获益的系统中的应用。
7.上文所述应用中，所述wars2-it1可为wars2-it1(gene id:104472716,updated on 24-jul-2021)。所述ac005162可为ac005162(genbank:ac005162.1，updated on 31-jan-2004)。所述ac130352可为ac130352(genbank:ac130352.6，updated on 13-dec-2002)。所述ac015971(genbank:ac015971.4，updated on 16-apr-2005)。所述foxp1-it1可为foxp1-it1(gene id:100506815,updated on 24-jul-2021)。所述ac040963可为ac040963(genbank:ac040963.6，updated on jun-26-2002)。所述linc00355可为linc00355(gene id:144766,updated on 16-aug-2021)。
8.上文所述应用中，所述物质可包括检测所述七种lncrna基因的表达量所需的试剂
和/或仪器。
9.上文所述应用中，所述预测小细胞肺癌患者预后的系统可包括所述物质和lncrna基因表达量数据处理模块。所述lncrna基因表达量数据处理模块可用于将来自待预测小细胞肺癌患者的分离的小细胞肺癌组织中所述七种lncrna基因的表达量转换为风险值，根据所述风险值预测所述待预测小细胞肺癌患者的预后。
10.上文所述应用中，所述预测小细胞肺癌患者化疗获益的系统也可包括所述物质和lncrna基因表达量数据处理模块。所述lncrna基因表达量数据处理模块用于将来自待预测小细胞肺癌患者的分离的小细胞肺癌组织中所述七种lncrna基因的表达量转换为风险值，根据所述风险值预测所述待预测小细胞肺癌患者是否从化疗中获益。
11.上文所述应用中，所述lncrna基因表达量数据处理模块可根据公式i得到所述风险值。
12.所述公式i可为：风险值＝(0.5622
×
wars2-it1基因表达量) (1.0842
×
ac005162基因表达量) (1.1170
×
ac130352基因表达量) (1.5938
×
ac015971基因表达量)-(0.6460
×
foxp1-it1基因表达量)-(0.0665
×
ac040963基因表达量)-(0.5835
×
linc00355基因表达量)。
13.上文所述应用中，所述lncrna基因表达量数据处理模块可利用r语言软件中的survminer软件包中的surv_cutpoint确定所述风险值的阈值。
14.为了解决上述技术问题，本发明还提供了七种lncrna作为标志物在预测小细胞肺癌患者的预后中的应用。所述七种lncrna可为wars2-it1、ac005162、ac130352、ac015971、foxp1-it1、ac040963和linc00355。
15.为了解决上述技术问题，本发明还提供了七种lncrna作为标志物在预测小细胞肺癌患者的化疗获益中的应用。所述七种lncrna可为wars2-it1、ac005162、ac130352、ac015971、foxp1-it1、ac040963和linc00355。
16.上文所述应用中，所述wars2-it1可为wars2-it1(gene id:104472716,updated on 24-jul-2021)。所述ac005162可为ac005162(genbank:ac005162.1，updated on 31-jan-2004)。所述ac130352可为ac130352(genbank:ac130352.6，updated on 13-dec-2002)。所述ac015971(genbank:ac015971.4，updated on 16-apr-2005)。所述foxp1-it1可为foxp1-it1(gene id:100506815,updated on 24-jul-2021)。所述ac040963可为ac040963(genbank:ac040963.6，updated on jun-26-2002)。所述linc00355可为linc00355(gene id:144766,updated on 16-aug-2021)。
17.为了解决上述技术问题，本发明还提供了预测小细胞肺癌患者预后和/或预测小细胞肺癌患者化疗获益的装置。所述装置可包括用于获取wars2-it1基因、ac005162基因、ac130352基因、ac015971基因、foxp1-it1基因、ac040963基因和linc00355基因这七种lncrna基因的表达量的模块。
18.上文所述装置中，所述wars2-it1可为wars2-it1(gene id:104472716,updated on 24-jul-2021)。所述ac005162可为ac005162(genbank:ac005162.1，updated on 31-jan-2004)。所述ac130352可为ac130352(genbank:ac130352.6，updated on 13-dec-2002)。所述ac015971(genbank:ac015971.4，updated on 16-apr-2005)。所述foxp1-it1可为foxp1-it1(gene id:100506815,updated on 24-jul-2021)。所述ac040963可为
ac040963(genbank:ac040963.6，updated on jun-26-2002)。所述linc00355可为linc00355(gene id:144766,updated on 16-aug-2021)。
19.上文所述装置还可包括所述lncrna基因表达量数据处理模块。所述lncrna基因表达量数据处理模块可用于将来自待预测小细胞肺癌患者的分离的小细胞肺癌组织中所述八种lncrna基因的表达量转换为风险值，根据所述风险值预测所述待预测小细胞肺癌患者的预后和/或是否从化疗中获益。
20.上文所述的预测小细胞肺癌患者预后的系统和/或预测小细胞肺癌患者化疗获益的系统也属于本发明的保护范围。
21.上文所述的小细胞肺癌组织可来自所述待预测小细胞肺癌患者的分离的小细胞肺癌组织经过福尔马林固定石蜡包埋制备的样本。
22.上文所述化疗可为辅助化疗。
23.上文所述七种lncrna可均来自于人(homo sapiens)。
24.本发明从geo数据库(https://www.ncbi.nlm.nih.gov/geo/)中下载gse65002队列作为训练队列建立了基于sclc中7个免疫相关lncrna(包括wars2-it1、ac005162、ac130352、ac015971、foxp1-it1、ac040963和linc00355)的分子模型(风险值预测模型)来预测sclc患者的化疗反应和预后风险。并在国家癌症中心收集158例小细胞肺癌患者福尔马林固定石蜡包埋的ffpe组织样本作为验证队列，用于模型验证。结果表明本发明建立的n6-甲基腺苷(m6a)相关lncrna的风险值预测模型在对不同生存时间、不同临床亚组的小细胞肺癌患者的总生存期或无复发生存期的预后预测和化疗获益预测中均得到了很好的验证，预测能力都是最好的；7个免疫相关lncrna可作为小细胞肺癌的预后预测和辅助化疗获益预测的标志物。该发明可为未来化疗在sclc患者中的临床应用提供参考依据，可用于指导可进行化疗的患者，使不受益者免于化疗，避免过度治疗，将有助于提高临床小细胞肺癌的治疗疗效。
附图说明
25.图1为sclc中最具预后意义的m6a相关lncrna的筛选。(a)为基于gse40275数据库中30个m6a调节元件的表达谱的sclc和正常肺组织的主成分分析。(b)为来自gse40275的30个m6a调节元件表达的热图。(c)为训练队列中m6a相关lncrna与sclc预后之间关联的森林图。
26.图2为训练队列中m6a相关lnrna预测模型的构建。(a)为模型中m6a相关lncrna的lasso系数分布。横坐标为系数(b)为m6a相关lnrna与风险值(risk score)之间的相关性矩阵图。横坐标的颜色代表相关性的高低。(c)为训练队列中患者生存状态的风险值分布，上图为患者风险值，纵坐标为风险值，横坐标为不同的患者；中图为训练队列患者预后情况，纵坐标为生存时间，横坐标为不同的患者，深色点表示患者死亡，浅色点表示存活；下图为训练队列中七个m6a相关lncrna的表达分布。(d)为训练队列患者的os的kaplan-meier曲线。纵坐标为总生存率，横坐标为患者术后生存时间(月)。(e)为m6a相关lncrna特征的roc分析。
27.图3为预测模型在独立队列的验证和临床应用。(a)为独立队列中158名患者os的kaplan-meier曲线。纵坐标为总生存率，横坐标为患者术后生存时间(月)(b)为独立队列中
os的风险值和多个临床特征的roc分析。纵坐标为灵敏性，横坐标为1-特异性(c)独立队列中患者os的风险值和不同临床参数的c指数分析。纵坐标为不同临床参数，横坐标为c指数值。(d)为独立队列中158名患者rfs的kaplan-meier曲线。纵坐标为无复发生存率，横坐标为患者术后生存时间(月)(e)为独立队列中rfs的风险值和多个临床特征的roc分析。纵坐标为灵敏性，横坐标为1-特异性(f)为独立队列患者rfs的风险值和不同临床参数的c指数分析。纵坐标为不同临床参数，横坐标为c指数值。(g)为独立队列接受辅助化疗的亚组中os的kaplan-meier曲线。纵坐标为无复发生存率，横坐标为患者术后生存时间(月)(h)为独立队列接受辅助化疗的亚组中os的roc分析。纵坐标为灵敏性，横坐标为1-特异性(i)为独立队列接受辅助化疗的亚组的os的风险值和不同临床参数的c指数分析。纵坐标为不同临床参数，横坐标为c指数值。(j)为独立队列接受辅助化疗的亚组中rfs的kaplan-meier曲线。(k)为独立队列接受辅助化疗的亚组中rfs的roc分析。纵坐标为灵敏性，横坐标为1-特异性(l)为独立队列接受辅助化疗的亚组的rfs的风险值和不同临床参数的c指数分析。纵坐标为不同临床参数，横坐标为c指数值。
28.图4为各临床亚组中预测模型的os预测性能。(a)验证队列中男性患者os的kaplan-meier曲线。(b)验证队列中女性患者os的kaplan-meier曲线。(c)验证队列中年老患者os的kaplan-meier曲线。(d)验证队列中青年患者os的kaplan-meier曲线。(e)验证队列中吸烟患者os的kaplan-meier曲线。(f)验证队列中非吸烟患者os的kaplan-meier曲线。纵坐标为总生存率，横坐标为术后时间。
29.图5为各临床亚组中预测模型的rfs预测性能。(a)验证队列中男性患者rfs的kaplan-meier曲线。(b)验证队列中女性患者rfs的kaplan-meier曲线。(c)验证队列中年老患者rfs的kaplan-meier曲线。(d)验证队列中青年患者rfs的kaplan-meier曲线。(e)验证队列中吸烟患者rfs的kaplan-meier曲线。(f)验证队列中非吸烟患者rfs的kaplan-meier曲线。纵坐标为无复发生存率，横坐标为术后时间。
30.图6为预测模型对生存期和化疗获益的预测能力分析。(a)训练队列和独立队列中风险值和多个临床特征的单变量cox回归分析。(b)训练队列和独立队列中风险评分和多个临床特征的多变量cox回归分析(c)独立队列接受辅助化疗的亚组中风险评分和多个临床特征的单变量cox回归分析。(d)独立队列接受辅助化疗的亚组中风险评分和多个临床特征的多变量cox回归分析。纵坐标为hr的对数值，横坐标为不同临床特征。
具体实施方式
31.下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。以下提供的实施例可作为本技术领域普通技术人员进行进一步改进的指南，并不以任何方式构成对本发明的限制。
32.下述实施例中的实验方法，如无特殊说明，均为常规方法，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
33.实施例1、小细胞肺癌患者m6a相关lncrna预测模型的建立和应用
34.从geo数据库(https://www.ncbi.nlm.nih.gov/geo/)中下载gse65002队列作为训练队列，用于构建小细胞肺癌预测模型；并在国家癌症中心(ncc)收集小细胞肺癌患者福
尔马林固定石蜡包埋的ffpe组织作为验证队列，用于模型验证。所有患者的临床特征如表1所示。
35.表1小细胞肺癌患者的临床特征
[0036][0037]
注：sclc：小细胞肺癌；os：总生存期。
[0038]
验证队列中，小细胞肺癌患者的入选标准均如下：1)患者经过国家癌症中心诊断为小细胞肺癌；2)患者诊疗信息完善；3)患者术前未接受辅助治疗；4)患者病理组织经过2位临床病理医师的独立评估，均确认为小细胞肺癌。
[0039]
辅助化疗以铂类药物为基础。
[0040]
一、小细胞肺癌患者m6a相关lncrna预测模型的建立
[0041]
1.1m6a相关lncrna预测模型的建立
[0042]
(1)在训练队列中，为了建立小细胞肺癌患者m6a分子的治疗疗效及预后预测模型，采用单因素cox比例回归模型，研究m6a相关基因对总生存期(os)预后指标的影响。
[0043]
为了使预测模型更加优化和实用，采用逐步cox比例风险回归模型，最终构建出一个包括如下筛选获得的7个m6a相关的lncrna基因的预后模型(m6a相关的lncrna筛选过程见实施例2)：wars2-it1(gene id:104472716,updated on 24-jul-2021)、ac005162(genbank:ac005162.1，updated on 31-jan-2004)、ac130352(genbank:ac130352.6，updated on 13-dec-2002)、ac015971(genbank:ac015971.4，updated on 16-apr-2005)、foxp1-it1(gene id:100506815,updated on 24-jul-2021)、ac040963(genbank:ac040963.6，updated on jun-26-2002)和linc00355(gene id:144766,updated on 16-aug-2021)。
[0044]
(2)根据每个患者目的基因的相对表达量，通过lasso分析，得出如下公式i用于计算每个患者的风险值(图2中a)：
[0045]
风险值＝(0.5622
×
wars2-it1基因表达量) (1.0842
×
ac005162基因表达量)
(1.1170
×
ac130352基因表达量) (1.5938
×
ac015971基因表达量)-(0.6460
×
foxp1-it1基因表达量)-(0.0665
×
ac040963基因表达量)-(0.5835
×
linc00355基因表达量)
[0046]
下文中，将公式i简称为风险值预测模型。
[0047]
(3)风险值预测模型中，各m6a相关lncrna和风险值的关系如图2中b所示。
[0048]
1.3m6a相关lncrna预测模型的验证
[0049]
(1)通过r语言软件的“survminer”软件包的“surv_cutpoint”确定阈值，具体方法如下：将待预测小细胞肺癌患者的风险值与匹配的预后信息，输入至r语言软件中，在“survminer”软件包的“surv_cutpoint”的算法下，软件会自动计算出p值最小的分割点，该分割点即为高风险组和低风险组的阈值(最优cutoff点)。
[0050]
结果显示，训练队列中，计算出的高风险组和低风险组的阈值为3.4995，患者风险值大于等于3.4995的为高风险组，患者风险值小于3.4995的为低风险组。(图2中c)
[0051]
(2)采用kaplan-meier生存分析法分析训练队列的高风险组和低风险组患者的总生存率os差异
[0052]
利用kaplan-meier分析上述(1)得到的高风险组和低风险组的风险值和生存数据。kaplan-meier生存分析结果显示，训练队列中高风险组患者的os比低分患者短(p《0.001)(图2中d)。
[0053]
(3)风险值预测模型的roc分析
[0054]
将训练队列患者的1年、3年和5年预后情况与风险值预测模型进行roc检验，结果显示，该风险值预测模型对患者多时间段的预后均具有很好的预测价值，具体的，1年：auc＝0.775；2年：auc＝0.826；5年：auc＝0.86(图2中e)。
[0055]
二、小细胞肺癌患者m6a相关lncrna预测模型对预后预测能力的验证
[0056]
为了验证步骤一中建立的风险值预测模型的可重复性和效果，用验证队列对风险值预测模型的有效性进行验证。
[0057]
2.1收集158例sclc患者(表1)的ffpe组织作为验证队列，提取ffpe组织样本的rna，通过pcr检测基因的表达量。
[0058]
具体检测方法如下：
[0059]
1)组织样本处理
[0060]
a.取小细胞肺癌和癌旁组织，各100mg，放入2ml无酶ep管中，加入1ml rnaiso plus(中国大连宝生物(takara)公司)；
[0061]
b.加4粒灭菌后的钢珠，放入高速低温组织研磨仪，参数设置为50hz，匀浆处理5min，去除组织匀浆液；
[0062]
c.用4℃离心机，用12000rpm/min离心10min，吸取上清液加入新的ep管中，置于冰上。
[0063]
2)rna浓度测定
[0064]
a.打开nanodrop，选择rna检测模式，用1μl无酶水清洗探头3次，吸水纸擦干；
[0065]
b.用1μl无酶水，校准仪器，调零；
[0066]
c.加1μl待测样品至探头，检测rna浓度，检测完毕后，用吸水纸擦干；
[0067]
d.重复上一步骤，至所有样品检测完成；
[0068]
e.用无酶水清洗探头3次。
[0069]
3)rna质控
[0070]
通过nanodrop检测的rna浓度，同时观察以下两个值：
[0071]
a.a260/a280比值：为rna浓度与蛋白浓度的比值，1.8-2.0间表示质控合格；
[0072]
b.a260/a230比值：为rna浓度与共提取污染的比值，1.8-2.2间表示质控合格。
[0073]
4)反转录
[0074]
用cdna反转录试剂盒(中国大连宝生物(takara)公司)，具体如下：
[0075]
a.去除基因组dna：体系详见表2.
[0076]
表2去除基因组dna体系
[0077][0078]
在冰上进行配置反应体系，混匀并短暂离心，42℃反应2min。
[0079]
b.反转录反应：体系详见表3。
[0080]
表3反转录体系
[0081][0082]
将反应体系短暂离心并放置于pcr仪中，程序为37℃，15min，85℃反应5s，获得ffpe组织样本的cdna。
[0083]
4)pcr扩增
[0084]
采用sybr green试剂(中国大连宝生物(takara)公司)，通过循环阈值(cycle threshold valve,ct)和标准曲线对起始模板进行定量分析。所需要的基因特异性引物均由捷瑞科技有限公司合成，检测各个目的基因和gapdh基因的引物序列如表4(序列1-序列16)所示。
[0085]
以3)获得的cdna为模板，每个基因表达的检测需要设置三个复孔，其反应体系如下(表5)：
[0086]
表4目的基因引物序列
[0087][0088]
表5 pcr扩增体系
[0089][0090]
按照上述体系进行配置，并在冰上避光加入到pcr仪器专用的八联排离心管中，涡旋混匀。注意轻弹无气泡后离心。
[0091]
在abi 7900ht荧光定量pcr仪中，设置反应条件如下：95℃10min；95℃10s，60℃10s，72℃10s，共40个循环，溶解曲线。以gapdh作为内参，按公式δct＝ct gene-ct gapdh进行数据分析，folds＝2-δδct
公式计算目的基因的相对表达量。
[0092]
2.2按照步骤1.2中(2)的方法分别对验证队列的患者ffpe样本的7个基因的相对表达量进行检测，计算风险值。
[0093]
2.3风险值预测模型对验证队列总生存率(overall survival，os)的预后预测
[0094]
(1)按照步骤1.3中(1)的方法，将验证队列的患者样本的风险值与os，输入至r语言软件中确定os阈值。结果确定阈值为0.1220，根据阈值，将患者分为高风险组和低风险组，患者风险值大于等于0.1220的为高风险组，患者风险值小于0.1220的为低风险组。
[0095]
(2)采用kaplan-meier生存分析法分析高风险组和低风险组患者的总生存率os差异
[0096]
利用kaplan-meier分析上述(1)得到的高风险组和低风险组的风险值和生存数据。kaplan-meier生存分析os,结果显示，验证队列中高风险组(图3中a的high risk)患者的os比低分(图3中a的low risk)患者短(p《0.001)(图3中a)。
[0097]
(3)不同临床亚组中风险值预测模型的效用验证
[0098]
在验证队列不同的临床亚组(性别、是否吸烟)中，风险值预测模型效用也得到了很好的验证。
[0099]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行对5年os的roc分析，结果证实，与患者年龄、性别、是否吸烟和肿瘤分期等比较，预测模型预测能力最好(图3中b)。具体的，风险值预测模型(图3中b的risk score)：auc＝0.645、
性别(图3中b的sex)：auc＝0.548、年龄(图3中b的age)：auc＝0.582、是否吸烟(图3中b的smoking)：auc＝0.579和分期(图3中b的staging)：auc＝0.593。对临床常用的病理参数，年龄(图3中c的age)、性别(图3中c的sex)、是否吸烟(图3中c的smoking)和肿瘤分期(图3中c的sclc_staging)与预测模型(图3中c的risk score)同时进行os的多因素cox分析，结果显示风险值预测模型对os的预测能力最好(图3中c)。
[0100]
2.4预测模型对验证队列无复发生存期(relapse-free survival,rfs)的预后预测
[0101]
(1)按照步骤1.3中(1)方法，将验证队列的患者样本的风险值与rfs预后信息，输入至r语言软件中确定rfs阈值。结果确定阈值为0.1220，根据阈值，将患者分为高风险组和低风险组，患者风险值大于等于0.1220的为高风险组，患者风险值小于0.1220的为低风险组。
[0102]
(2)采用kaplan-meier生存分析法分析高风险组和低风险组患者的总生存率os差异
[0103]
利用kaplan-meier分析上述(1)得到的高风险组和低风险组的风险值和生存数据。kaplan-meier生存分析os,结果显示，ncc验证队列中高风险组(图3中d的high risk)患者的os比低分(图3中d的low risk)患者短(p《0.001)(图3中d)。
[0104]
(3)不同临床亚组中风险值预测模型的效用验证
[0105]
在验证队列不同的临床亚组(性别、是否吸烟)中，风险值预测模型效用也得到了很好的验证。
[0106]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行对5年os的roc分析，结果证实，与患者年龄(图3中e的age)、性别(图3中e的sex)、是否吸烟(图3中e的smoking)和肿瘤分期(图3中e的staging)等比较，风险值预测模型(图3中e的risk score)预测能力最好(图3中e)。具体的，风险值预测模型：auc＝0.666、性别：auc＝0.581、年龄：auc＝0.54、是否吸烟：auc＝0.591和肿瘤分期：auc＝0.608。对临床常用的病理参数，年龄(图3中f的age)、性别(图3中f的sex)、是否吸烟(图3中f的smoking)和肿瘤分期(图3中f的sclc_staging)与风险值预测模型(图3中f的risk score)同时进行os的多因素cox分析，结果显示风险值预测模型对os的预测能力最好(图3中f)。
[0107]
三、小细胞肺癌患者m6a相关lncrna预测模型对辅助化疗获益程度预测能力的验证
[0108]
为了验证预测模型对辅助化疗获益程度预测的能力，选取表1的验证队列中接受辅助化疗的138例患者对预测模型的有效性进行验证。
[0109]
3.1风险值预测模型对辅助化疗后os获益的验证
[0110]
(1)按照步骤1.3中(1)方法，将验证队列138例接受辅助化疗的的患者样本的风险值与辅助化疗后os，输入至r语言软件中确定辅助化疗后os阈值。结果确定阈值为0.0143，根据阈值，将患者分为高风险组和低风险组，患者风险值大于等于0.0143的为高风险组，患者风险值小于0.0143的为低风险组。
[0111]
(2)采用kaplan-meier生存分析法分析高风险组和低风险组患者的总生存率os差异
[0112]
利用kaplan-meier分析上述(1)得到的高风险组和低风险组的风险值和生存数
据。kaplan-meier生存分析os,结果显示，验证队列中高风险组(图3中g的high risk)患者的os比低分(图3中g的low risk)患者短(p《0.001)(图3中g)。
[0113]
(3)不同临床亚组中预测模型的效用验证
[0114]
在验证队列不同的临床亚组(性别、是否吸烟)中，风险值预测模型效用也得到了很好的验证。
[0115]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行对5年os的roc分析，结果证实，与患者年龄(图3中h的age)、性别(图3中h的sex)、是否吸烟(图3中h的smoking)和肿瘤分期(图3中h的staging)等比较，风险值预测模型(图3中h的risk score)预测能力最好(图3中h)。具体的，风险值预测模型：auc＝0.629、性别auc＝0.525、年龄auc＝0.589、是否吸烟auc＝0.569和分期auc＝0.596。对临床常用的病理参数，年龄(图3中i的age)、性别(图3中i的sex)、是否吸烟(图3中i的smoking)和肿瘤分期(图3中i的sclc_staging)与风险值预测模型(图3中i的risk score)同时进行os的多因素cox分析，结果显示风险值预测模型对os的预测能力最好(图3中i)。
[0116]
3.2风险值预测模型对辅助化疗后rfs获益的验证
[0117]
(1)按照步骤1.3中(1)方法，将验证队列138例接受辅助化疗的患者样本的风险值与辅助化疗后rfs，输入至r语言软件中确定rfs阈值。结果确定阈值为0.0143，根据阈值，将患者分为高风险组和低风险组，患者风险值大于等于0.0143的为高风险组，患者风险值小于0.0143的为低风险组。
[0118]
(2)采用kaplan-meier生存分析法分析高风险组和低风险组患者的总生存率os差异
[0119]
利用kaplan-meier分析上述(1)得到的高风险组和低风险组的风险值和生存数据。kaplan-meier生存分析os,结果显示，验证队列中高风险组(图3中j的high risk)患者的os比低分(图3中j的low risk)患者短(p《0.001)(图3中j)。
[0120]
(3)不同临床亚组中预测模型的效用验证
[0121]
在验证队列不同的临床亚组(性别、是否吸烟)中，风险值预测模型效用也得到了很好的验证。
[0122]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行对5年os的roc分析，结果证实，与患者年龄(图3中k的age)、性别(图3中k的sex)、是否吸烟(图3中k的smoking)和肿瘤分期(图3中k的staging)等比较，风险值预测模型(图3中k的risk score)预测能力最好(图3中k)。具体的，风险值预测模型：auc＝0.635、性别：auc＝0.55、年龄：auc＝0.565、是否吸烟：auc＝0.59和分期：auc＝0.631。对临床常用的病理参数，年龄(图3中l的age)、性别(图3中l的sex)、是否吸烟(图3中l的smoking)和肿瘤分期(图3中l的sclc_staging)与风险值预测模型(图3中l的risk score)同时进行os的多因素cox分析，结果显示预测模型对os的预测能力最好(图3中l)。
[0123]
3.3预测模型对不同临床亚组的kaplan-meier分析
[0124]
(1)通过对sclc患者的不同亚组的os进行kaplan-meier分析，结果显示，在年龄、性别和是否吸烟等临床亚组中，高风险组患者的os比低分患者短。具体的，在男性患者中，高风险组(图4中a的high risk)患者的os比低分(图4中a的low risk)患者短(p《0.001)；在女性患者中，高风险组(图4中b的high risk)患者的os比低分(图4中b的low risk)患者短
(p＝0.013)；在老年患者中，高风险组(图4中c的high risk)患者的os比低分(图4中c的low risk)患者短(p《0.001)；在青年患者中，高风险组(图4中d的high risk)患者的os比低分(图4中d的low risk)患者短(p＝0.001)；在吸烟患者中，高风险组(图4中e的high risk)患者的os比低分(图4中e的low risk)患者短(p＝0.001)；在非吸烟患者中，高风险组(图4中f的high risk)患者的os比低分(图4中f的low risk)患者短(p＝0.002)(图4)。
[0125]
(2)通过对sclc患者的不同亚组的rfs进行kaplan-meier分析，结果显示，在年龄、性别和是否吸烟等临床亚组中，高风险组患者的os比低分患者短。具体的，在男性患者中，高风险组患者(图5中a的high risk)的os比低分(图5中a的low risk)患者短(p《0.001)；在女性患者中，高风险组(图5中b的high risk)患者的os比低分(图5中b的low risk)患者短(p＝0.001)；在老年患者中，高风险组(图5中c的high risk)患者的os比低分(图5中c的low risk)患者短(p＝0.001)；在青年患者中，高风险组(图5中d的high risk)患者的os比低分(图5中d的low risk)患者短(p＝0.001)；在吸烟患者中，高风险组(图5中e的high risk)患者的os比低分(图5中e的low risk)患者短(p《0.001)；在非吸烟患者中，高风险组(图5中f的high risk)患者的os比低分(图5中f的low risk)患者短(p＝0.002)(图5)。
[0126]
四、小细胞肺癌患者m6a相关lncrna风险值预测模型预测能力分析
[0127]
4.1、风险值预测模型对sclc患者预后预测能力验证
[0128]
(1)风险值预测模型对训练队列中os的预测能力
[0129]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行os的单因素cox分析，结果显示，相比于性别(图6中a的traning cohort组sex行：p＝0.947，95％ci：0.287-3.341)、年龄(图6中a的traning cohort组age行：p＝0.890，95％ci：0.407-2.181)、是否吸烟(图6中a的traning cohort组smoking行：p＝0.389,95％ci：0.593-3.831)和肿瘤分期(图6中a的traning cohort组sclc staging行：p＝0.002,95％ci：1.706-10.374)，风险值预测模型对os的预测能力最好(图6中a的traning cohort组risk score行：p《0.001，95％ci：3.166-58.485)。
[0130]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行os的多因素cox分析，结果显示，相比于性别(图6中b的traning cohort组sex行：p＝0.426，95％ci：0.337-13.145)、年龄(图6中b的traning cohort组age行：p＝0.228，95％ci：0.709-4.232)、是否吸烟(图6中b的traning cohort组smoking行：p＝0.501,95％ci：0.162-2.432)和肿瘤分期(图6中b的traning cohort组sclc staging行：p＝0.095,95％ci：0.857-6.887)，风险值预测模型对os的预测能力最好(图6中b的traning cohort组risk score行：p＝0.004，95％ci：2.093-51.652)。
[0131]
(2)预测模型对验证队列中os的预测能力
[0132]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行os的单因素cox分析，结果显示，相比于性别(图6中a的validation cohort(os)组sex行：p＝0.615，95％ci：0.684-1.898)、年龄(图6中a的validation cohort(os)组age行：p＝0.105，95％ci：0.929-2.156)、是否吸烟(图6中a的validation cohort(os)组smoking行：p＝0.334,95％ci：0.799-1.933)和肿瘤分期(图6中a的validation cohort(os)组sclc staging行：p＝0.001,95％ci：1.183-2.007)，预测模型对os的预测能力最好(图6中a的validation cohort(os)组risk score行：p《0.001，95％ci：1.871-4.737)。
[0133]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与预测模型同时进行os的多因素cox分析，结果显示，相比于性别(图6中b的validation cohort(os)组sex行：p＝0.891，95％ci：0.489-1.862)、年龄(图6中b的validation cohort(os)组age行：p＝0.094，95％ci：0.939-2.214)、是否吸烟(图6中b的validation cohort(os)组smoking行：p＝0.210，95％ci：0.810-2.611)和肿瘤分期(图6中b的validation cohort(os)组sclc staging行：p＝0.008,95％ci：1.101-1.882)，风险值预测模型对os的预测能力最好(图6中b的validation cohort(os)组risk score行：p＝0.004，95％ci：2.093-51.652)。
[0134]
(3)预测模型对验证队列中rfs的预测能力
[0135]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行rfs的单因素cox分析，结果显示，相比于性别(图6中a的validation cohort(rfs)组sex行：p＝0.153，95％ci：0.873-2.379)、年龄(图6中a的validation cohort(rfs)组age行：p＝0.453，95％ci：0.785-1.721)、是否吸烟(图6中a的validation cohort(rfs)组smoking行：p＝0.140,95％ci：0.902-2.076)和肿瘤分期(图6中a的validation cohort(rfs)组sclc staging行：p＝0.003,95％ci：1.132-1.865)，风险值预测模型对os的预测能力最好(图6中a的validation cohort(rfs)组risk score行：p《0.001，95％ci：1.755-4.070)。
[0136]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行rfs的多因素cox分析，结果显示，相比于性别(图6中b的validation cohort(rfs)组sex行：p＝0.407，95％ci：0.694-2.462)、年龄(图6中b的validation cohort(rfs)组age行：p＝0.463，95％ci：0.778-1.737)、是否吸烟(图6中b的validation cohort(rfs)组smoking行：p＝0.266，95％ci：0.793-2.312)和肿瘤分期(图6中b的validation cohort(rfs)组sclc staging行：p＝0.035,95％ci：1.020-1.689)，风险值预测模型对os的预测能力最好(图6中b的validation cohort(rfs)组risk score行：p《0.001，95％ci：1.763-4.165)。
[0137]
4.2、预测模型对接受辅助化疗sclc患者获益程度的预测能力验证
[0138]
(1)预测模型对验证队列中辅助化疗后os的预测能力
[0139]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与预测模型同时进行os的单因素cox分析，结果显示，相比于性别(图6中c的validation cohort(os)组sex行：p＝0.982，95％ci：0.594-1.705)、年龄(图6中c的validation cohort(os)组age行：p＝0.115，95％ci：0.916-2.247)、是否吸烟(图6中c的validation cohort(os)组smoking行：p＝0.439,95％ci：0.757-1.900)和肿瘤分期(图6中c的validation cohort(os)组sclc staging行：p＝0.002,95％ci：1.177-2.087)，风险值预测模型对os的预测能力最好(图6中c的validation cohort(os)组risk score行：p《0.001，95％ci：1.651-4.288)。
[0140]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与预测模型同时进行os的多因素cox分析，结果显示，相比于性别(图6中d的validation cohort(os)组sex行：p＝0.596，95％ci：0.414-1.660)、年龄(图6中d的validation cohort(os)组age行：p＝0.180，95％ci：0.866-2.152)、是否吸烟(图6中d的validation cohort(os)组smoking行：p＝0.258，95％ci：0.773-2.609)和肿瘤分期(图6中d的validation cohort(os)组sclc staging行：p＝0.014,95％ci：1.077-1.934)，风险值预测模型对os的预测能力最好(图6中
d的validation cohort(os)组risk score行：p＝0.004，95％ci：1.499-3.996)。
[0141]
(2)预测模型对验证队列中辅助化疗后rfs的预测能力
[0142]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与预测模型同时进行rfs的单因素cox分析，结果显示，相比于性别(图6中c的validation cohort(rfs)组sex行：p＝0.320，95％ci：0.775-2.183)、年龄(图6中c的validation cohort(rfs)组age行：p＝0.396，95％ci：0.788-1.824)、是否吸烟(图6中a的validation cohort(rfs)组smoking行：p＝0.130,95％ci：0.905-2.169)和肿瘤分期(图6中c的validation cohort(rfs)组sclc staging行：p＝0.005,95％ci：1.127-1.937)，风险值预测模型对os的预测能力最好(图6中c的validation cohort(rfs)组risk score行：p《0.001，95％ci：1.605-4.063)。
[0143]
对临床常用的病理参数，年龄、性别、是否吸烟和肿瘤分期与风险值预测模型同时进行rfs的多因素cox分析，结果显示，相比于性别(图6中d的validation cohort(rfs)组sex行：p＝0.822，95％ci：0.557-2.091)、年龄(图6中d的validation cohort(rfs)组age行：p＝0.558，95％ci：0.741-1.741)、是否吸烟(图6中d的validation cohort(rfs)组smoking行：p＝0.302，95％ci：0.766-2.360)和肿瘤分期(图6中d的validation cohort(rfs)组sclc staging行：p＝0.064,95％ci：0.985-1.725)，预测模型对os的预测能力最好(图6中d的validation cohort(rfs)组risk score行：p《0.001，95％ci：1.463-3.792)。
[0144]
综上所述，本发明通过确定sclc中m6a相关lncrna的表达谱，并建立了基于7个m6a相关lncrna(包括wars2-it1、ac005162、ac130352、ac015971、foxp1-it1、ac040963和linc00355)的分子模型来预测sclc患者的化疗反应和预后风险，结果表明本发明所建立的m6a相关lncrna的风险值预测模型可以有效对小细胞肺癌的预后和化疗获益进行预测；7个m6a相关lncrna可作为小细胞肺癌的预后预测和化疗获益预测的标志物。该发明可能会为未来化疗在sclc患者中的临床应用提供参考依据，将有助于提高临床小细胞肺癌的治疗疗效。
[0145]
实施例2、m6a相关lncrna的筛选过程
[0146]
1.1m6a相关lncrna的筛选
[0147]
(1)m6a各调节元件在肿瘤组织和正常组织中的表达
[0148]
本发明从已发表文献中((1).li y,xiao j,bai j,tian y,qu y,chen x,wang q,li x,zhang y,and xu j,molecular characterization and clinical relevance of m(6)a regulators across 33cancer types.mol cancer,2019.18(1):137.(2).liu j,harada bt,and he c,regulation of gene expression by n(6)-methyladenosine in cancer.trends cell biol,2019.29(6):487-499.(3).huang h,weng h,and chen j,m(6)a modification in coding and non-coding rnas:roles and therapeutic implications in cancer.cancer cell,2020.37(3):270-288.(4).nombela p,miguel-l
ó
pez b,and blanco s,the role of m(6)a,m(5)candψrna modifications in cancer:novel therapeutic opportunities.mol cancer,2021.20(1):18.)，总结了30个m6a的调节元件，包括11个writer调节元件(mettl3,mettl14,mettl16,mettl5,wtap,virma,rbm15,rbm15b,zc3h13,cbll1,和zcchc4),2个eraser调节元件(fto和alkbh5),以及17各reader调节元件(ythdf1,ythdf2,ythdf3,ythdc1,ythdc2,hnrnpa2b1,hnrnpc,fmr1,eif3a,igf2bp1,igf2bp2,igf2bp3,elavl1,g3bp1,g3bp2,prrc2a和rbmx)。
[0149]
通过从geo数据库(https://www.ncbi.nlm.nih.gov/geo/)中下载gse40275数据，对比m6a调节元件在肿瘤组织和正常组织中的表达。结果显示，m6a调节元件在肿瘤组织和正常组织中表达差异明显(图1中a)。通过热图分析发现，m6a调节元件在肿瘤组织中呈现高表达(图1中b)。
[0150]
(2)m6a相关lncrna的筛选
[0151]
从geo数据库(https://www.ncbi.nlm.nih.gov/geo/)中下载gse60052数据作为训练队列进行分析。首先排除了训练队列中表达值为0的lncrna，之后将剩余的lncrna与m6a调节元件分别进行pearson分析。同时满足|pearson r|》0.5和p《0.0001条件的lncrna被筛选为m6a相关lncrna。最终，有289个lncrna被选中。
[0152]
(3)m6a相关lncrna与预后的相关性
[0153]
选取训练队列中有预后信息的患者，分别分析上步筛选出的m6a相关lncrna对预后的影响。结果显示，有19个m6a相关lncrna对患者总生存期存在相关性(p《0.2)(图1中c)。
[0154]
为了使构建的小细胞肺癌预测模型更加优化和实用，最终筛选出如下7个lncrna基因用于构建模型：wars2-it1(gene id:104472716,updated on 24-jul-2021)、ac005162(genbank:ac005162.1，updated on 31-jan-2004)、ac130352(genbank:ac130352.6，updated on 13-dec-2002)、ac015971(genbank:ac015971.4，updated on 16-apr-2005)、foxp1-it1(gene id:100506815,updated on 24-jul-2021)、ac040963(genbank:ac040963.6，updated on jun-26-2002)和linc00355(gene id:144766,updated on 16-aug-2021)。
[0155]
以上对本发明进行了详述。对于本领域技术人员来说，在不脱离本发明的宗旨和范围，以及无需进行不必要的实验情况下，可在等同参数、浓度和条件下，在较宽范围内实施本发明。虽然本发明给出了特殊的实施例，应该理解为，可以对本发明作进一步的改进。总之，按本发明的原理，本技术欲包括任何变更、用途或对本发明的改进，包括脱离了本技术中已公开范围，而用本领域已知的常规技术进行的改变。按以下附带的权利要求的范围，可以进行一些基本特征的应用。