基于铁死亡相关lncRNA构建肝细胞癌预后模型

0 引言

肝细胞癌(hepatocellular carcinoma,HCC)是起源于肝脏上皮组织的恶性肿瘤,约占原发性肝癌的75%-85%,具有高致死率和高转移复发率的特性

.HCC为高度异质性肿瘤,具有高度侵袭性、表型耐药和代谢紊乱等生物学特征,多层次制约了总体临床治疗效果

.因此,迫切需要探究HCC发生发展机制、新治疗策略和潜在生存预后预测模型,以指导临床个体化治疗方案的制订.

铁死亡是近年来发现的一类新型铁离子依赖性调节性细胞死亡模式,其主要是通过铁离子的积累、自由基的产生、脂肪酸供应与脂质过氧化物增加等途径诱导细胞死亡

.目前研究表明

,铁死亡在HCC的发生发展过程中发挥着重要调控作用,铁死亡基因靶向治疗有可能成为治疗HCC的新手段.长链非编码RNA(long non-coding RNA,lncRNA)在基因表达过程中发挥着关键调控作用,与铁死亡过程存在着密切关系

.目前HCC中基于铁死亡基因作为生物标志物研究鲜见,本研究旨在通过探究铁死亡相关lncRNA与HCC患者预后之间的关系,构建生存预后预测模型,为其个体化治疗和免疫研究提供依据.

另外，质量表由按顺序排列的质量表页组成，每一页记载着一个或多个质量项。根据工序与质量项的对应关系，可以得出这样的结论：一个质量表页对应一个或多个质检类工序，而当某个质检类工序对应的质量项较多时，可能对应多个质量表页。

1 材料和方法

1.1 材料 通过癌症基因组图谱(TCGA,https://portal.gdc.cancer.gov)数据库下载HCC患者正常样本和肿瘤样本转录组数据及其临床相关信息,资料截止至2021-10-29.转录组数据共得到HCC肿瘤样本374例,正常样本50例;临床相关样本共370例,其中临床信息完整样本350例.转录组数据通过Perl软件注释并区分为mRNA和lncRNA.铁死亡基因表达数据与mRNA取交集获取,其中铁死亡基因通过铁死亡数据库(Ferrdb)获取.通过共表达分析方法获得与铁死亡基因密切相关的lncRNA(r＞0.5,

＜0.05).

1.2 方法

1.2.1 HCC预后相关铁死亡相关lncRNA筛选:利用R软件(4.1.2版)对共表达分析方法获得的铁死亡相关lncRNA进行差异分析.具体筛选标准设定为|logFC|＞1,

＜0.05.差异铁死亡相关lncRNA与临床相关信息中生存时间和生存状态合并,利用“survival”包对合并数据进行单因素Cox回归分析(

＜0.05),获得各个差异铁死亡相关lncRNA风险值(hazard ratio,HR).单因素Cox回归分析获得lncRNA进一步行多因素Cox分析,获得HCC各样本风险评分和影响HCC患者的关键铁死亡相关lncRNA.

1.2.3 预后风险评分模型验证:本研究利用三种方式对基于铁死亡相关lncRNA预后风险评分模型进行验证.利用R软件“survival”和“survminer”包对高低风险组进行Kaplan-Meier生存分析,并绘制相关生存曲线.单多因素Cox回归分析,验证风险评分能否作为评估HCC患者生存预后的独立风险因子.利用R软件“timeROC”包验证风险评分模型的特异性和灵敏度,并绘制受试者工作特征曲线(receiver operator characteristic curve,ROC).

1.2.4 免疫功能分析:整合预后风险评分模型中11个铁死亡相关lncRNA的表达量及相应样本风险评分.利用ssGSEA富集分析评估高低风险组中免疫反应情况.

采用

检验分析铁死亡相关lncRNA在正常样本和肿瘤样本中的表达差异.共表达分析采用皮尔逊相关系数(Pearson correlation coefficient)法.生存曲线通过Kaplan-Meier方法绘制,运用对数秩检验(Logrank test)进行分析.

2 结果

根据上述原理分析,结合甲烷气体的吸收曲线,本文选择中心波长为1 330 nm的DFB半导体激光器(绵阳全光通信科技有限公司的DF-14-DFB1310)作为发射光源,采用数字锁相放大器提取一次谐波和二次谐波,并利用其比值计算甲烷气体浓度,实现方案如图1所示。

A1组(浓度为1. 25%的大豆卵磷脂稀释液)低温保存绵羊精液精子活率均高于其他浓度组(P<0. 05)；A1组低温保存绵羊精液精子活率第9 d(0. 53±0. 06)、12 d(0. 46±0. 09)与C组第9 d(0. 55±0. 03)、12 d(0. 49±0. 07)差异不显著(P﹥0. 05)，且顶体完整率相当(P﹥0. 05)；A1组低温保存第9 d的精液进行人工授精,受胎率64. 3%与C组65. 6%差异不显著(P>0. 05)。表明1. 25%浓度的大豆卵磷脂稀释液能够进行绵羊精液低温长时间有效保存。

2.2 HCC预后相关关键铁死亡相关lncRNA筛选结果HCC中呈差异表达铁死亡相关lncRNA,与生存时间和生存状态合并后,纳入单因素Cox回归分析.筛选条件设置为

＜0.05,统计分析发现49个lncRNA与生存预后显著相关.多因素Cox回归将单因素Cox回归分析结果再次进行筛选,共得到11个lncRNA.结果如表1所示.

2.3 预后风险评分模型构建结果 根据多因素Cox回归分析筛选得到lncRNA构建预后风险评分模型.风险评分计算方法为11个lncRNA表达水平与风险系数的乘积.具体风险评分模型为:风险评分=LINC01224*0.569+SREBF2-AS1*0.430+NRAV*0.153+KDM4A-AS1*0.654+ZFPM2-AS1*0.043+SNHG20*(-0.287)+SNHG4*0.232+MKLN1-A S*0.982+P T O V1-A S1*(-0.502)+L I N C009-42*0.023+ZNF232-AS1*0.066.

2.4 预后风险评分模型验证结果 根据预后风险评分模型,计算各个样本的风险值.风险值中位数(0.852)作为截断值,根据截断值将370例HCC样本分为高、低风险组,每组各185例.Kaplan-Meier生存分析结果显示,高、低风险组平均生存时间分别为(1.76±1.85)年和(2.50±2.04)年,两组生存时间存在显著差异(

＜0.001),见图1.Cox回归分析评价预后风险评分模型对HCC生存状态的影响.单因素Cox回归分析结果显示,风险评分模型与生存时间和生存状态显著相关[风险比(hazard ratio,HR)=1.211;95%置信区间(confidence interval,CI)=(1.159-1.266);

＜0.001].进一步行多因素Cox回归分析结果显示,预后风险评分模型是预测生存情况的独立危险因子(HR=1.190;95%CI=1.135-1.247;

＜0.001),见表2.ROC分析结果显示,风险评分模型预测HCC患者1年、2年和3年生存率的ROC曲线下面积(area under curve,AUC)分别为0.784、0.780和0.743,如图2.ROC曲线进一步分析,临床特征与风险评分模型对HCC患者生存率的影响,结果提示风险评分模型和TNM分期对HCC患者生存具有显著影响,AUC值分别为0.784和0.672,风险评分模型优于TNM分期,如图3.

乘客使用手机模拟银行卡直接进出站并延迟扣款，不需要联机到后台进行账户认证，根据进出站交易配对结算票款。容易产生单边交易，进出站交易不匹配不能结算票款。通常采用车站级服务器或云端服务器控制用户的进出站状态，需要对所有车站进行网络改造，改造难度大、后台压力较大、票务优惠策略逻辑处理难度较大。目前广州地铁和上海磁浮线都已开通ODA应用。

3 讨论

HCC是我国常见的恶性肿瘤,目前其发生和发展的具体分子机制尚不明确.近年来研究发现,铁死亡基因和lncRNA在多种恶性肿瘤的发生发展过程中发挥着重要调控作用,如胶质母细胞瘤

、肺癌

、卵巢癌

和胰腺导管腺癌

等,并构建了相关模型预测预后.然而,铁死亡相关lncRNA在HCC分子层面中的机制尚不明确,同时预测模型构建方面研究较少.

铁死亡在HCC发生和发展过程中发挥着重要调控作用.目前研究已证实,胱氨酸-谷氨酸反向转运体(system XC)和GPX4可以调控铁死亡.Tang等研究表明,索拉菲尼是常见的铁死亡诱导剂,其可以通过抑制system XC系统和影响GSH,诱导肝癌细胞发生铁死亡

.肿瘤抑制基因p53亦可以通过调控铁死亡,影响HCC的发展.Chen等

研究提示,p53可以通过调控SLC7A11转录,影响system XC系统活性,诱发肝癌细胞铁死亡.此外,miRNA-214、BRCA1相关蛋白1等分子均可影响肝癌细胞铁死亡

.

2.1 铁死亡基因及相关lncRNA提取结果 TCGA数据库下载得到HCC肿瘤样本(374例)和正常样本(50例)的转录组数据和临床相关信息.Ferrdb数据库下载得到铁死亡基因259个,并提取铁死亡基因表达量.按照FDR＜0.05且logFC＞1标准分析,筛选得到差异表达铁死亡基因共84个,表达上调及下调基因数目分别为71个和13个.按照r＞0.5且

＜0.05筛选标准,共表达分析得到铁死亡相关lncRNA共233个,其中呈差异表达的lncRNA共147个.

本研究通过利用生物信息学技术,分析了233个铁死亡相关lncRNA.最终得到11个铁死亡相关lncRNA,并成功构建了HCC生存预后预测模型.我们得出,高风险组总体生存率显著低于低风险组(

＜0.001).预测HCC患者1年、2年及3年生存率的ROC曲线下面积均大于0.7,提示该模型具有较好的特异性和敏感性.我们进一步通过检索文献发现,11个预后lncRNA均与肿瘤的发生发展关系密切.Sun等

研究显示,LINC01224通过调控MiR-193a-5p/CDK8轴影响胃癌细胞凋亡和肿瘤进展.Chen等

通过将KDM4A-AS1过表达肝癌细胞株Hep3B移植到裸鼠皮下,发现过表达组肿瘤组织内Ki-67、N-钙粘蛋白和波形蛋白表达水平较敲减组明显升高,肿瘤负荷明显增大.进一步研究发现,KDM4A-AS1通过调控KPNA2/HIF-1α信号通路,促进肝癌的发生和发展.提示KDM4A-AS1在HCC发挥促癌作用.Ren等研究发现

,甲状腺癌组织中ZFPM2-AS1较癌旁组织显著升高;同时体外实验研究提示,ZFPM2-AS1过表达可以促进甲状腺癌细胞株SW579增殖和侵袭.Xing等

通过分析30对卵巢癌组织,发现SNHG20在卵巢癌组织中呈过表达状态;SNHG20过表达促进SKOV3和A2780细胞株增殖和侵袭.以上研究均提示,铁死亡相关lncRNA在恶性肿瘤的进展过程中发挥着重要调控作用.然而,铁死亡相关lncRNA目前在HCC中的研究罕见,其影响HCC发生及发展的分子机制和其在HCC临床实践中的指导作用有待进一步探究.

近年来研究表明,铁死亡在恶性肿瘤免疫调控过程中发挥着重要作用.铁死亡在T细胞、B细胞、粒细胞及单核巨噬细胞等免疫细胞调控中发挥着重要调控作用

.Wang等在肿瘤小鼠模型中研究发现,CD8+T细胞通过分泌干扰素γ诱导肿瘤细胞中SLC3A2 和SLC7A11的表达,从而降低肿瘤细胞对胱氨酸的摄取,致使肿瘤细胞发生铁死亡

.Muri Jonathan等实验表明,敲除B1细胞和边缘区B细胞中铁死亡关键调控基因GPX4,可以诱导细胞铁死亡,从而影响B细胞的免疫应答

.本研究分析提示,铁死亡相关lncRNA可以通过调控免疫细胞溶解活性、Ⅱ型干扰素应答及MHC Ⅱ类分子途径影响HCC免疫状态,从而影响HCC患者预后.

4 结论

综上,本研究筛选得到11个与HCC患者预后关系密切的铁死亡相关lncRNA,其构成的预后风险评分模型对HCC患者预后具有较好的预测价值.同时与其它临床指标相结合,可能为制定个体化治疗方案或随访方案提供更多的理论依据.本课题组下一步拟设计分子实验,验证铁死亡相关lncRNA在HCC发生、发展中的确切分子机制.

根据稳态分析计算结果，可以看出2种工况的管道压力均满足设计要求，但各支管处出水压力均较高，工况1管道末端压力为83 m，工况2末端出水压力水头达70.6 m。

2.5 免疫功能分析 高低风险组ssGSEA富集分析结果显示,高风险组中免疫细胞溶解活性和Ⅱ型干扰素应答呈抑制状态,而高风险组中MHC Ⅱ类分子途径呈高反应状态,均具有统计学差异(

＜0.05),如图4.

复方冬虫夏草口服液对载脂蛋白E基因敲除小鼠动脉粥样硬化的影响 …………………………………… 胡抗远等（14）：1912

肝细胞癌是我国高发恶性肿瘤之一,70%-80%首诊时已是晚期,居我国恶性肿瘤相关死亡原因第三位.尽管近年来,外科手术治疗、介入治疗、化学药物治疗等传统方法和靶向、免疫治疗等创新疗法均取得了显著进步,但是肝细胞癌患者总体生存期仍待进一步提升.因此,探究肝细胞癌发生发展机制、新治疗策略和潜在生存预后预测模型,来指导临床个体化治疗显得十分迫切.

江苏苏北有些地方水稻还在零星追肥，都在严阵以待打好3次防治总体战。第一次7月底到8月上旬，沿太湖、沿江地区以“两迁”害虫、纹枯病为主，其他地区纹枯病、钻心虫为主攻对象；第二次8月中旬，主攻纹枯病、稻纵卷叶螟、稻飞虱、螟虫等；第三次8月下旬到9月上中旬，主攻稻曲病、穗稻瘟、纹枯病及“两迁”害虫、螟虫等，3个攻坚战忙得不可开交，理“药”不理“肥”无可厚非啊！

本研究通过生物信息学方法,探究铁死亡相关lncRNA与肝细胞癌患者预后之间的关系,并构建生存预后预测模型,为其个体化治疗和免疫研究提供依据.

铁死亡是一种铁依赖性的、非细胞凋亡性的新型细胞死亡方式.研究表明,铁死亡基因及其相关lncRNA可以调控多种恶性肿瘤的发生发展过程,如胶质母细胞瘤、肺癌、卵巢癌和胰腺导管腺癌等,并初步构建了相关模型预测预后.然而,铁死亡相关lncRN肝细胞癌分子层面中的机制尚不明确,同时预测模型构建方面研究较少.

检索癌症基因组图谱数据库,获取肝细胞癌患者正常样本和肿瘤样本转录组数据及其临床相关信息.采用共表达分析方法获得与铁死亡基因密切相关lncRNA.采用单、多因素Cox回归分析方法获取与肝细胞癌患者预后密切相关铁死亡相关lncRNA,构建基于铁死亡相关lncRNA预后风险评分模型,并进行验证.

分析癌症基因组图谱数据库下载数据,筛选共获得147个铁死亡相关lncRNA.单、多因素Cox回归分析分别获得49个和11个与HCC生存预后相关的铁死亡相关lncRNA.生存分析结果显示,高风险组生存率显著低于低风险组(

＜0.001).独立预后分析结果表明,模型风险评分是预测HCC患者预后的独立风险因素.免疫功能分析结果提示,高风险组中免疫细胞溶解活性和Ⅱ型干扰素应答呈抑制状态,而高风险组中MHC Ⅱ类分子途径呈高反应状态.

铁死亡相关lncRNA模型风险评分是预测肝细胞癌患者预后的独立风险因素.机制可能时通过调控免疫细胞溶解活性、Ⅱ型干扰素应答及MHC Ⅱ类分子途径发挥调控作用.

2.2.2 噪音 Kahn［18］等认为，噪声主要来自谈话、监护报警、呼吸机等，这3类噪声分别占噪声总量的26%、20%和8%。噪音超过60分贝就会导致患者烦躁不安、刺激患者的交感神经，使心率加快、血压升高、压力感和焦虑感加重，疼痛感加剧，使患者感到抑郁、头痛、幻觉、入睡困难和昼夜睡眠节律倒转［19］。

编辑出版学学科的制度化没有达到根深蒂固的程度，加上体系重构是学科发展的必然要求，不可抗阻，因此重构的阻力是可以克服的。学科体系结构是描述一个学科体系构成的基本逻辑层级，是一种框架性描述，结构化的学科体系的优点在于可以根据学科逻辑起点知识的变迁规律与外界对知识的现实需求与时俱进的进行调整，一条线索是知识的变迁演进，一条线索是知识的规训与权利。从结构层面看，学科体系建构可秉持“三三”原则——三个方面、三个维度，即体系上可以大致划分为知识体系、学科课程体系、人才培养体系三个方面；话语生产又分为知识生产、规则生产和价值生产三个维度。

铁死亡相关lncRNA模型有望成为临床评估肝细胞癌患者预后的评分工具之一,同时该模型中相关lncRNA可能成为肝细胞癌靶向或免疫治疗的潜在靶点.

1 Rodriguez S,Skeet K,Mehmetoglu-Gurbuz T,Goldfarb M,Karri S,Rocha J,Shahinian M,Yazadi A,Poudel S,Subramani R.Phytochemicals as an Alternative or Integrative Option,in Conjunction with Conventional Treatments for Hepatocellular Carcinoma.

2021;13 [PMID:34830907 DOI:10.3390/cancers13225753]

2 Myojin Y,Kodama T,Maesaka K,Motooka D,Sato Y,Tanaka S,Abe Y,Ohkawa K,Mita E,Hayashi Y,Hikita H,Sakamori R,Tatsumi T,Taguchi A,Eguchi H,Takehara T.ST6GAL1 Is a Novel Serum Biomarker for Lenvatinib-Susceptible FGF19-Driven Hepatocellular Carcinoma.

2021;27:1150-1161 [PMID:33288659 DOI:10.1158/1078-0432.CCR-20-3382]

3 Arechederra M,Bazai SK,Abdouni A,Sequera C,Mead TJ,Richelme S,Daian F,Audebert S,Dono R,Lozano A,Gregoire D,Hibner U,Allende DS,Apte SS,Maina F.ADAMTSL5 is an epigenetically activated gene underlying tumorigenesis and drug resistance in hepatocellular carcinoma.

2021;74:893-906 [PMID:33197513 DOI:10.1016/j.jhep.2020.11.008]

4 Chen X,Yu C,Kang R,Kroemer G,Tang D.Cellular degradation systems in ferroptosis.

2021;28:1135-1148 [PMID:33462411 DOI:10.1038/s41418-020-00728-1]

5 Gao R,Kalathur RKR,Coto-Llerena M,Ercan C,Buechel D,Shuang S,Piscuoglio S,Dill MT,Camargo FD,Christofori G,Tang F.YAP/TAZ and ATF4 drive resistance to Sorafenib in hepatocellular carcinoma by preventing ferroptosis.

2021;13:e14351 [PMID:34664408 DOI:10.15252/emmm.202114351]

6 Tang R,Wu Z,Rong Z,Xu J,Wang W,Zhang B,Yu X,Shi S.Ferroptosis-related lncRNA pairs to predict the clinical outcome and molecular characteristics of pancreatic ductal adenocarcinoma.

2022;23 [PMID:34553745 DOI:10.1093/bib/bbab388]

7 Chen Q,Wang W,Wu Z,Chen S,Chen X,Zhuang S,Song G,Lv Y,Lin Y.Over-expression of lncRNA TMEM161B-AS1 promotes the malignant biological behavior of glioma cells and the resistance to temozolomide via up-regulating the expression of multiple ferroptosis-related genes by sponging hsa-miR-27a-3p.

2021;7:311 [PMID:34689169 DOI:10.1038/s41420-021-00709-4]

8 Yao J,Chen X,Liu X,Li R,Zhou X,Qu Y.Characterization of a ferroptosis and iron-metabolism related lncRNA signature in lung adenocarcinoma.

2021;21:340 [PMID:34217273 DOI:10.1186/s12935-021-02027-2]

9 Peng J,Hao Y,Rao B,Zhang Z.A ferroptosis-related lncRNA signature predicts prognosis in ovarian cancer patients.

2021;10:4802-4816 [PMID:35116333 DOI:10.21037/tcr-21-1152]

10 Tang H,Chen D,Li C,Zheng C,Wu X,Zhang Y,Song Q,Fei W.Dual GSH-exhausting sorafenib loaded manganesesilica nanodrugs for inducing the ferroptosis of hepatocellular carcinoma cells.

2019;572:118782 [PMID:31678528 DOI:10.1016/j.ijpharm.2019.118782]

11 Chen X,Kang R,Kroemer G,Tang D.Broadening horizons:the role of ferroptosis in cancer.

2021;18:280-296 [PMID:33514910 DOI:10.1038/s41571-020-00462-0]

12 Zhang Y,Shi J,Liu X,Feng L,Gong Z,Koppula P,Sirohi K,Li X,Wei Y,Lee H,Zhuang L,Chen G,Xiao ZD,Hung MC,Chen J,Huang P,Li W,Gan B.BAP1 links metabolic regulation of ferroptosis to tumour suppression.

2018;20:1181-1192 [PMID:30202049 DOI:10.1038/s41556-018-0178-0]

13 Sun H,Yan J,Tian G,Chen X,Song W.LINC01224 accelerates malignant transformation via MiR-193a-5p/CDK8 axis in gastric cancer.

2021;10:1377-1393 [PMID:33655711 DOI:10.1002/cam4.3726]

14 Chen T,Liu R,Niu Y,Mo H,Wang H,Lu Y,Wang L,Sun L,Wang Y,Tu K,Liu Q.HIF-1α-activated long non-coding RNA KDM4A-AS1 promotes hepatocellular carcinoma progression via the miR-411-5p/KPNA2/AKT pathway.

2021;12:1152 [PMID:34903711 DOI:10.1038/s41419-021-04449-2]

15 Ren R,Du Y,Niu X,Zang R.ZFPM2-AS1 transcriptionally mediated by STAT1 regulates thyroid cancer cell growth,migration and invasion via miR-515-5p/TUSC3.

2021;12:3393-3406 [PMID:33976749 DOI:10.7150/jca.51437]

16 Xing X,An M,Chen T.LncRNA SNHG20 promotes cell proliferation and invasion by suppressing miR-217 in ovarian cancer.

2021;43:1095-1104 [PMID:34302635 DOI:10.1007/s13258-021-01138-4]

17 Wang P,Lu YQ.Ferroptosis:A Critical Moderator in the Life Cycle of Immune Cells.

2022;13:877634 [PMID:35619718 DOI:10.3389/fimmu.2022.877634]

18 Wang W,Green M,Choi JE,Gijón M,Kennedy PD,Johnson JK,Liao P,Lang X,Kryczek I,Sell A,Xia H,Zhou J,Li G,Li J,Li W,Wei S,Vatan L,Zhang H,Szeliga W,Gu W,Liu R,Lawrence TS,Lamb C,Tanno Y,Cieslik M,Stone E,Georgiou G,Chan TA,Chinnaiyan A,Zou W.CD8

T cells regulate tumour ferroptosis during cancer immunotherapy.

2019;569:270-274 [PMID:31043744 DOI:10.1038/s41586-019-1170-y]

19 Muri J,Thut H,Bornkamm GW,Kopf M.B1 and Marginal Zone B Cells but Not Follicular B2 Cells Require Gpx4 to Prevent Lipid Peroxidation and Ferroptosis.

2019;29:2731-2744.e4 [PMID:31775041 DOI:10.1016/j.celrep.2019.10.070]