长链非编码RNA LNC01309对人肝癌细胞增殖和迁移能力的影响及作用机制
DOI: 10.3969/j.issn.1001-5256.2022.03.014
Effect of long non-coding RNA LNC 01309 on proliferation and migration abilities of human hepatoma cells and its mechanism of action
-
摘要:
目的 探讨长链非编码RNA(lncRNA)LNC01309对人肝癌细胞增殖和迁移能力的影响及作用机制。 方法 收集2018年2月—2019年6月在解放军总医院第六医学中心手术治疗的12例肝细胞癌(以下简称肝癌)患者标本及对应的癌旁组织,采用荧光定量PCR方法检测LNC01309的相对表达量。利用荧光定量PCR法检测LNC01309在人肝癌细胞系(Hep G2、SNU-398和Hep 3B)和人永生化正常肝细胞系(THLE-2)中的表达水平。在肝癌细胞Hep G2中过表达LNC01309,分为质粒对照组(pEXP-control)和过表达组(pEXP-LNC01309)。采用CCK-8检测各组细胞增殖变化,采用划痕愈合试验和Transwell试验检测各组细胞迁移能力。利用RNA免疫共沉淀试验检测LNC01309与RBM38的相互作用(分组为IgG组和RBM38抗体组)。利用放线菌酮追踪分析检测LNC01309对于RBM38蛋白稳定性的影响。在肝癌细胞Hep G2中过表达RBM38进行回复试验,采用CCK-8试验、划痕愈合试验和Transwell试验,分别检测细胞增殖和迁移能力的变化。计量资料两组间比较采用t检验。 结果 LNC01309在肝癌组织中的平均表达水平高于癌旁组织(4.225±2.285 vs 1.541±0.530, t=3.618,P=0.004)。LNC01309在肝癌细胞(Hep G2、SNU-398和Hep 3B)中的相对表达水平亦均高于正常肝细胞(THLE-2)(t值分别为4.231、6.489、14.480,P值均<0.05)。过表达LNC01309的Hep G2细胞(t=9.172,P<0.001)相对于对照组,细胞的生长速度明显上升(第96小时OD450值:1.885±0.107 vs 2.527±0.234,t=4.330,P=0.012),迁移能力上调(划痕愈合试验:11.65%±2.40% vs 35.66%±4.90%,t=9.837,P<0.001;Transwell试验:100.00%±3.11% vs 161.00%±35.93%,t=4.399,P=0.005);然而过表达LNC01309诱导的肝癌细胞增殖和迁移能力的上调效应,都能够被RBM38部分抑制(第96小时OD450值:2.500±0.227 vs 1.913±0.282,t=2.812,P=0.048;168.00%±9.43% vs 117.20%±18.03%,t=6.622,P<0.001)。相对于IgG对照组,RBM38抗体能够显著富集沉淀LNC01309(t=3.846,P=0.031)。放线菌酮试验结果显示,过表达LNC01309组细胞中RBM38蛋白稳定性显著下调(t=8.038,P=0.001)。 结论 新发现的LNC01309通过与RBM38的相互作用降低了RBM38蛋白稳定性,促进肝癌细胞的增殖与迁移。 -
关键词:
- 癌, 肝细胞 /
- RNA, 长链非编码 /
- RNA结合蛋白质类
Abstract:Objective To investigate the effect of long non-coding RNA (lncRNA) LNC01309 on the proliferation and migration abilities of human hepatocellular carcinoma (HCC) cells and its mechanism of action. Methods HCC samples and corresponding adjacent tissue samples were collected from 12 patients with HCC who underwent surgical treatment in The Sixth Medical Center of PLA General Hospital from February 2018 to June 2019, and quantitative real-time PCR was used to measure the relative expression level of LNC01309. Quantitative real-time PCR was also used to measure the expression level of LNC01309 in human hepatoma cell lines (HepG2, SNU-398, and Hep3B) and the human immortalized normal liver cell line THLE-2. After LNC01309 was overexpressed in HepG2 cells, the cells were divided into plasmid control group (pEXP-control) and overexpression group (pEXP-LNC01309). CCK-8 assay was used to observe the change in cell proliferation, and wound healing assay and Transwell assay were used to observe migration ability. RNA co-immunoprecipitation was used to detect the interaction between LNC01309 with RBM38, with cells divided into IgG group and RBM38 antibody group, and cycloheximide chase assay was used to analyze the effect of LNC01309 on the stability of RBM38 protein. RBM38 was overexpressed in HepG2 cells to conduct the recovery experiment, and CCK-8 assay, wound healing assay, and Transwell assay were used to observe the changes in cell proliferation and migration abilities. The t-test was used for comparison of continuous data between two groups. Results The mean expression level of LNC01309 in HCC tissue was significantly higher than that in adjacent tissue (4.225±2.285 vs 1.541±0.530, t=3.618, P=0.004), and the relative expression level of LNC01309 in hepatoma cells (HepG2, SNU-398, and Hep3B) was also significantly higher than that in normal hepatocytes (THLE-2) (t=4.231、6.489、14.480, all P < 0.05). Compared with the control group, HepG2 cells with the overexpression of LNC01309 had significant increases in growth rate (OD450 value at 96 hours: 1.885±0.107 vs 2.527±0.234, t=4.330, P=0.012) and migration ability (11.65%±2.40% vs 35.66%±4.90%, t=9.837, P < 0.001; 100.00%±3.11% vs 161.00%±35.93%, t=4.399, P=0.005); however, the upregulated proliferation and migration abilities of hepatoma cells induced by LNC01309 overexpression were partially inhibited by RBM38 (OD450 value at 96 hours: 2.500±0.227 vs 1.913±0.282, t=2.812, P=0.048; 168.00%±9.43% vs 117.20%±18.03%, t=6.622, P < 0.001). Compared with the IgG control group, RBM38 antibody significantly enriched the precipitation of LNC01309 (t=3.846, P=0.031). The results of cycloheximide chase assay showed that the LNC01309 overexpression group had a significant reduction in the stability of RBM38 protein (t=8.038, P=0.001). Conclusion The newly identified LNC01309 reduces the stability of RBM38 protein through interaction with RBM38 and promotes the proliferation and migration of HCC cells. -
Key words:
- Carcinoma, Hepatocellular /
- RNA, Long Noncoding /
- RNA-Binding Proteins
-
肝细胞癌(以下简称肝癌)是全球发病率和病死率均较高的恶性肿瘤。近年来肝癌的诊断和治疗取得了较大进展,但晚期肝癌患者预后仍较差。目前肝癌的发病机制尚未完全阐明,深入分析肝癌的发生、发展机制,将有助于发现新的肝癌诊断标志物和治疗靶点。长链非编码RNA(long non-coding RNA,lncRNA)在肿瘤形成过程中的表达机制颇受关注,已被公认参与多种肿瘤的发生、发展及预后,与肝癌亦关系密切[1-2]。
在人类基因组中,参与编码蛋白质的基因仅占2%,超过98%基因组不编码蛋白质,其基因转录产物被称为非编码RNA(non-coding RNA,ncRNA)[3]。根据其长度,可将ncRNA分为miRNA、snoRNA、piRNA和lncRNA等[4]。ncRNA中大多数为lncRNA,lncRNA的基因长度往往>200 bp,其开放阅读框一般短于50~100 bp。根据lncRNA在基因组中与编码蛋白质基因的相对位置,可将其分成5类,即反义lncRNA、正义lncRNA、双向lncRNA、基因内lncRNA和基因间lncRNA[5]。随着二代测序技术的快速发展,成千上万的lncRNA在脊椎动物和非脊椎动物体内被快速鉴别出来。但是,二代测序技术在实际应用中也存在很多缺陷,如文库制备过程中PCR扩增很大程度会使后续结果出现误差,而且尽管已经开发出各种算法来提高测序结果的质量,但其对于组装一个横跨整个基因组的长片段仍然是一项挑战。lncRNA通常长度较长且存在多种剪接形式,单纯的RNA-Seq并不能分辨各种异构体,因此第三代测序技术的诞生就显得尤为重要[6]。
PacBio作为最新一代的测序技术,是由美国Pacific Bioscience开发的单分子实时测序(molecule real time sequencing,SMRT)技术,测序过程无需PCR扩增。其在读取长度方面远超于二代测序技术,并且无GC偏好性,能够大大提升基因组组装的质量[7]。PacBio在研究转录组学方面最大的优势在于它能够获得全长或接近全长的转录本,这对可变剪切的鉴定、lncRNA的分析[8-9]以及新基因的发现等研究有着重要意义;且能在原有研究基础上进行补充甚至发现新的机制。
本研究主要利用PacBio三代测序技术在肝癌组织中鉴定新的lncRNA,并在肝癌细胞中进行功能机制的初步探索。
1. 材料与方法
1.1 研究材料
收集2018年2月—2019年6月解放军总医院第六医学中心12例肝癌患者的肿瘤组织和对应的癌旁组织。人肝癌细胞系(Hep G2、SNU-398和Hep 3B)购自北京协和医学院基础学院基础医学细胞中心;人永生化正常肝细胞(THLE-2)购自上海泽叶生物科技有限公司;RPMI-1640及胎牛血清购自美国Thermo Fisher Scientific公司;Smarter PCR cDNA合成试剂盒购自美国Clontech公司;SYBR Green染料购自日本Takara公司;所有引物购自生工生物工程(上海)股份有限公司;CNOT4、EIF4G2、HNRNPA1、MBNL1、RBM24、RBM6、SRSF10抗体购自美国Santa Cruz Biotechnology公司;RBM38、N-cadherin、Vimentin、E-cadherin和GAPDH抗体及辣根过氧化物酶标记的二抗均购自英国Abcam公司;CCK-8试剂盒购自日本同仁化学公司;Transwell小室购自康宁公司;RNA荧光原位杂交(RNA fluorescence in situ hybridization,RNA FISH)和EdU检测试剂盒购自广州锐博生物科技有限公司;靶向lncRNA MIR205HG的siRNA由上海吉玛制药技术有限公司设计与合成;pCMV-RBM38购自北京义翘神州生物技术有限公司;pEXP-LNC01309野生型和突变型表达载体由广州锐博生物科技有限公司构建完成。
1.2 细胞培养与细胞转染
按照培养基和血清配比为9∶ 1配置完全培养基;所有细胞均采用RPMI-1640完全培养基培养;5% CO2,37 ℃恒温培养箱为细胞培养环境。细胞转染采用Lipofectamine 3000试剂。在6孔板中,接种细胞浓度为1×106/孔,培养基体积2 mL,其中siRNA用量为20 nmol/L/孔,质粒用量为8 μg/孔,共孵育24 h后,进行后续试验。
1.3 RNA提取
采用Trizol提取组织或细胞中的RNA。按100 mg组织加1 mL Trizol研磨,制备组织悬液;5×106细胞加入1 mL Trizol吹打,收获细胞悬液;且冰上裂解5~10 min后,按200 μL氯仿﹕1 mL Trizol的比例加入氯仿,混匀液体,静置15 min。然后4 ℃ 12 000 r/min离心分离RNA,吸取上清液至新的EP管中并弃去沉淀。加入等体积异丙醇,4 ℃静置10 min,12 000 r/min离心10 min,弃去上清液。加入75%酒精充分洗涤沉淀,离心去除上清液。最后将RNA样品沉淀置于超净台干燥5 min,加入DEPC水溶解,Nanodrop测量总RNA浓度,并置于-80 ℃保存。
1.4 PacBio三代测序
利用Trizol提取组织RNA。根据产品说明书,共使用4 μg混合RNA与Pacbio Sequel系统(Pacific Biosciences,CA,USA)进行全长转录序列测序。使用Clontech Smarter PCR cDNA合成试剂盒和Pacific Biosciences的luePippin Size Selection系统(PN100-092-800-03),根据异构体测序(Iso-Seq)协议制备Iso-Seq文库。后续测序与分析由北京基石生命科技有限公司完成。
1.5 荧光定量PCR
按照试剂盒说明书使用寡(dT)引物在总体积为20 μL的条件下进行cDNA合成;采用染料法,使用最终体积为25 μL的500 ng cDNA模板进行荧光定量PCR。重复3次试验,数据采用2-ΔΔCt方法处理,以GAPDH为内参。引物序列如下(5′-3′),lncRNA MIR205HG:正向引物AGCAGCAGCAGCAAGAGTAA,反向引物GAAAGATTAAGGTTCCCATC;LNC01309:正向引物AGCAGCAGCAGCAAGAGTAACT,反向引物GAAAGATTAAGGTTCCCATC;RBM38:正向引物GCTACGGCTTCGTGACCAT,反向引物GTAAGTCCGCTGGATCAAGGT;GAPDH:正向引物AGAAGGCTGGGGCTCATTTG,反向引物AGGGGCCATCCACAGTCTTC。
1.6 Western Blot
冰上裂解细胞并提取细胞中的总蛋白,BCA法检测总蛋白浓度,向50 μg蛋白样品中加入等体积的2×Laemmli上样缓冲液,总上样体积为20 μL,进行SDS-PAGE电泳。接着进行恒流(200 mA)转膜,脱脂奶粉封闭,4 ℃孵一抗过夜。次日进行二抗常温孵育,二抗孵育完成后加入化学发光液与蛋白条带反应,凝胶成像仪观察结果。
1.7 细胞增殖检测
使用CCK-8和EdU检测细胞增殖。CCK-8检测:按照1×104/孔细胞数接种至96孔细胞培养板,待细胞贴壁生长覆盖率达到50%左右,Lipofectamine 3000转染细胞并持续培养4 d。分别于不同时间点(0、24、48、72、96 h)弃去培养基,加入无血清培养基配置的10% CCK-8反应液,37 ℃恒温培养箱培养40 min,酶标仪检测波长在450 nm处的吸光度值。EdU检测:细胞接种于96孔板中,经转染处理24 h。然后将50 μmol/L EdU加入培养孔,37 ℃孵育2 h后,用4%多聚甲醛固定细胞。按照标准程序,anti-EdU工作液和DAPI染色孵育细胞。于荧光显微镜(奥林巴斯IX71)下观察,随机选取不同视野记录试验结果。
1.8 细胞划痕愈合试验
将各组Hep G2细胞按1×106/孔接种于6孔板中,置于37 ℃,5% CO2培养箱培养过夜。第2天用10 μL枪头划两条平行线,形成划痕。使用PBS轻轻漂洗2~3遍,加人无血清培养基。放入37 ℃,5% CO2培养箱继续培养,分别于0、12 h观察,拍照。用Image J分析迁移面积并计算迁移率。
1.9 Transwell迁移试验
将转染后24 h的各组细胞按照5×104细胞数分别接种于Transwell的上室,下室中加入800 μL含20%胎牛血清的培养基。继续培养24 h,取出小室,采用PBS清洗3次,甲醇固定15 min,用DAPI染色后在荧光显微镜下观察并随机选取5个以上视野拍照计数。
1.10 RNA FISH
使用预冷的Triton X-100(0.1%)使细胞膜通透,并用多聚甲醛(4%)处理载玻片。然后用LNC01309探针杂交过夜。SCC缓冲液清洗细胞。采用DAPI染细胞核后,荧光显微镜观察结果。
1.11 蛋白质的稳定性测定
用终浓度100 μg/mL的放线菌酮处理细胞,并在指定时间点(0、2、4 h)收取细胞。从细胞中提取蛋白质并进行蛋白质印迹分析。β-肌动蛋白用作参照。
1.12 统计学方法
采用Graphpad Prism 7.0进行统计学分析,Shapiro-Wilk检测试验数据均服从正态分布。计量资料用x±s表示,两组间比较采用t检验。试验数据均重复3次以上。P<0.05为差异有统计学意义。
2. 结果
2.1 lncRNA
LNC01309在肝癌组织和细胞中的表达在肝癌和对应的癌旁组织中利用PacBio测序分析发现1条高丰度且未见报道的lncRNA(图 1a),长度为309 nt,暂命名为LNC01309。利用荧光定量PCR检测分析,结果显示肝癌组织和对应的癌旁组织中LNC01309水平分别为4.225±2.285与1.541±0.530,差异有统计学意义(t=3.618,P=0.004)(图 1b)。将永生化正常肝细胞(THLE-2)和多种肝癌细胞(Hep G2、SNU-398和Hep 3B)中的LNC01309表达情况进行比较,结果显示,LNC01309在肝癌细胞中的相对含量分别为1.860±0.126、2.313±0.383和3.917±0.393,均明显高于THLE-2细胞(1.011±0.074) (t值分别为4.231、6.489、14.480,P值分别为0.004、<0.001、<0.001)(图 1c)。
2.2 LNC01309与lncRNA MIR205HG的关系分析
利用NCBI BLAST进行序列比对分析,LNC01309定位于1号染色体(209432225-209432545),与lncRNA MIR205HG存在较高的序列相似性(表 1)。lncRNA MIR205HG亦定位于1号染色体(209428820-209432549),其目前存在5种剪接体,通过序列比对分析,笔者设计了特异的PCR引物用于区分检测lncRNA MIR205HG和LNC01309(图 2a)。结果显示,lncRNA MIR205HG无论是在肝癌组织还是在肝癌细胞中均不存在异常表达现象(图 2b、c)。随后,在Hep 3B肝癌细胞中利用siRNA敲降lncRNA MIR205HG,分为NC siRNA对照组(NC si)和MIR205HG siRNA处理组(MIR205HG si)。荧光定量PCR结果显示,随着lncRNA MIR205HG的下调(1.001±0.082 vs 0.281±0.062,t=15.690,P<0.001)(图 2d),LNC01309并未明显下降(1.008±0.046 vs 1.075±0.085,t=1.534,P=0.164)(图 2e)。值得注意的是,LNC01309本身序列中并不含有miR-205序列,因此在肝癌细胞中过表达LNC01309,分为质粒对照组(pEXP-control)和过表达组(pEXP-LNC01309)(1.006± 0.099 vs 27.350±5.744,t=9.172,P<0.001)并没有导致miR-205相对含量的变化(0.971±0.083 vs 0.962±0.100,t=0.142,P=0.892)(图 2f、g)。
表 1 利用NCBI BLAST比较LNC01309与lncRNA MIR205HG序列相似性描述 种属 最高分 总分 查询覆盖值 E值 同源百分比 查询长度 登录号 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 3, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 992 NR_145434.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 1, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 895 NR_145437.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 2, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 857 NR_145433.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 4, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 940 NR_145435.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 5, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 651 NR_145436.1 图 2 LNC01309与lncRNA MIR205HG的关系分析注:a,LNC01309与lncRNA MIR205HG序列相似性,并根据序列差异性设计特异性检测引物;b,利用荧光定量PCR分析lncRNA MIR205HG在肝癌组织与癌旁组织的表达差异;c,利用荧光定量PCR分析lncRNA MIR205HG在肝癌细胞和正常肝细胞中的表达差异;d、e,Hep 3B细胞经转染siRNA(20 nmol/L)处理24 h后,利用荧光定量PCR分别检测lncRNA MIR205HG和LNC01309的相对表达含量;f、g,Hep G2细胞中转染过表达质粒24 h后,利用荧光定量PCR分别检测LNC01309和miR-205的相对表达含量。2.3 LNC01309促进肝癌细胞增殖与迁移
笔者尝试在含有较低丰度LNC01309的Hep G2细胞中过表达LNC01309(1.006±0.099 vs 27.350±5.744,t=9.172,P<0.001)。RNA FISH试验结果与在线工具Locate-R(http://locate-r.azurewebsites.net/)预测结果相一致,LNC01309主要定位于细胞质中(图 3、4)。随着LNC01309水平的升高,肝癌细胞的增殖能力(第96小时OD450值:1.885±0.107 vs 2.527±0.234,t=4.330,P=0.012)显著上调(图 5、6)。同时,划痕愈合试验(11.65%±2.40% vs 35.66%±4.90%,t=9.837,P<0.001)(图 7)和Transwell试验结果(100.00%±3.11% vs 161.00%±35.93%,t=4.399,P=0.005)(图 8)均显示肝癌细胞迁移能力增强。Western Blot结果显示,过表达LNC01309会促进上皮间质转化(EMT)的标志蛋白N- cadherin(0.466±0.067 vs 1.253±0.046,t=5.322,P=0.009)和Vimentin(0.372±0.126 vs 2.638±0.413,t=7.292,P=0.018)的表达,同时下调E-cadherin的表达(0.909±0.032 vs 0.151±0.014,t=50.450,P<0.001)(图 9)。
2.4 LNC01309与RNA结合蛋白的相互作用
利用在线工具RBPmap(http://rbpmap.technion.ac.il/index.html)和catRAPID(http://service.tartaglialab.com/page/catrapid_group)分析LNC01309的潜在结合蛋白,获得了CNOT4、EIF4G2、HNRNPA1、MBNL1、RBM24、RBM38、RBM6和SRSF10共8个排名靠前的候选蛋白(表 2)。随后,在Hep G2细胞中过表达LNC01309,利用候选蛋白的特异抗体进行RNA免疫共沉淀试验,富集结果分别为0.087±0.044、0.120±0.107、0.192±0.105、0.164± 0.093、0.153±0.108、0.502±0.213、0.124±0.085和0.168±0.068;与IgG对照组比较,仅RBM38抗体组差异有统计学意义(0.502±0.213 vs 0.060±0.029,t=3.846,P=0.031)(图 10a)。将RBM38潜在结合位点(gugugcg)缺失突变后构建突变型LNC01309。RNA免疫共沉淀实验结果显示,与IgG对照组比较,RBM38抗体组沉淀富集的野生型LNC01309相对含量显著升高(0.434±0.122 vs 0.034±0.017,t=5.123,P=0.036),而突变型LNC01309在两组中的沉淀富集丰度差异无统计学意义(0.077±0.051 vs 0.031±0.006,P=0.201)(图 10b)。
表 2 利用RBPmap和catRAPID在线预测能够与LNC01309结合的潜在靶标蛋白蛋白质 序列位置 基序 K-mer Z评分 P值 CNOT4 129 gacaga gagaga 3.554 1.90×10-4 EIF4G2 50 cgccg cgccg 4.390 5.67×10-6 HNRNPA1 66 guaguagu guagcagc 3.361 3.88×10-4 MBNL1 70 gcgcagc cagcagc 3.903 4.75×10-5 RBM24 160 wgwgugd ggagugu 3.403 3.33×10-4 RBM38 163 kkguguk gugugcg 3.333 4.30×10-4 RBM6 37 ccacc ccacc 4.228 1.18×10-5 SRSF10 70 cagcag cagcag 4.513 3.20×10-6 2.5 LNC01309影响RBM38蛋白稳定性
在相对低表达LNC01309的Hep G2细胞中过表达LNC01309时,并不会影响RBM38 mRNA的相对表达量(1.022±0.074 vs 0.953±0.076,P=0.185),但会降低RBM38蛋白的相对含量(0.960±0.057 vs 0.228±0.044,t=5.526,P=0.008)(图 11)。利用放线菌酮处理Hep G2细胞,LNC01309过表达组细胞中RBM38蛋白的稳定性显著低于对照质粒组(第4小时相对蛋白含量:0.610±0.075 vs 0.179±0.054,t=8.038,P=0.001)(图 12、13)。此外,过表达的RBM38能够有效抑制LNC01309的促细胞增殖能力(第96小时OD450值:2.500±0.227 vs 1.913±0.282,t=2.812,P=0.048)(图 14)。细胞划痕愈合试验结果(30.69%±5.80% vs 16.37%±1.66%,t=6.218,P<0.001)(图 15)和Transwell试验结果(168.00%±9.43% vs 117.20%±18.03%,t=6.622,P<0.001)(图 16)均显示过表达的RBM38能够有效抑制LNC01309诱导的肝癌细胞迁移能力的上调。
3. 讨论
大量的研究[10-12]工作已经证实lncRNA与肝癌的发生、发展及预后密切相关。例如,lncRNA ATB竞争性地结合miR-200,进而激活ZEB1和ZEB2的表达;ZEB1和ZEB2与IL-11 mRNA的相互作用增强了Stat3信号,促进肝癌转移,加速器官受累[13]。lncRNA HULC在肝癌中异常高表达[14],能够介导肿瘤抑制因子P18的下调,可促进肝癌的细胞增殖[15]。lncRNA HOTAIR降低了CREB、P300、RNApoⅢ在SETD2启动子区域的重新分配,从而抑制了SETD2的表达和磷酸化,借此促进人类肝癌干细胞的恶性生长[16-17]。目前发现lncRNA主要有以下几个方面的调控功能[18-19]。(1)表观遗传学:通过干扰转录因子与启动子结合、诱导蛋白质修饰和染色质发生重构等方式影响编码基因表达;(2)转录调控:与靶蛋白结合后影响基因转录;(3)转录后调控:能调控mRNA稳定性和翻译进程,作用蛋白结合伴侣或调节蛋白的亚细胞定位影响蛋白活性等,进而影响肿瘤进展[20];(4)编码微肽:lncRNA重新获得编码能力,产生具有功能的微肽[21]。笔者通过PacBio测序分析LNC01309显示,在肝癌组织和肝癌细胞中均呈现高表达状态。后续的试验结果发现LNC01309能够与RBM38蛋白发生相互作用,并借此下调RBM38蛋白稳定性。
RBM38是一种RNA结合蛋白,在正常细胞中通过调控RNA参与基因转录后调控,从而影响细胞增殖、细胞周期循环、肌源性分化等多种生物学过程。目前已知RBM38可与p53、mdm2、p63、p73、p21、c-Myc和MIC-1等多个目的基因的3′非翻译区结合,调节相应mRNA的稳定性,从而在肿瘤进展中发挥关键作用[22-27]。现有研究[28]显示,RBM38可能通过稳定p53-mdm2环抑制肝癌的发生。p53-mdm2通路在肝癌中经常发生改变,mdm2和p53基因突变可增加mdm2的稳定性,加速p53降解,破坏p53-mdm2平衡,从而诱导癌变,最终促进肝癌发生。Ye等[28]报道RBM38作为稳定p53-mdm2环的重要分子,在肝癌组织中RBM38失活,而mdm2表达增加,野生型p53表达随之降低,最终破坏p53-mdm2环的平衡;相反,若增加RBM38表达则可抑制肝癌细胞mdm2的表达,恢复野生型p53的表达,从而在体外抑制肝癌细胞迁移和侵袭;裸鼠体内实验也显示RBM38对肝癌的体内致瘤性具有一定抑制作用。以上体内外研究结果表明,上调RBM38可能通过抑制mdm2从而拯救野生型p53,进而改变肝癌的生物学功能和进展。本研究结果表明RBM38能够与野生型LNC01309发生结合作用,而不能够结合突变型LNC01309。LNC01309通过与RBM38的相互作用降低RBM38的蛋白稳定性,促进肝癌细胞的恶性表型特征。
LNC01309与lncRNA MIR205HG具有较高的序列相似性。然而敲降已知的lncRNA MIR205HG不同剪接体,并不会改变LNC01309的相对含量。这些结果提示,LNC01309并不是已知lncRNA MIR205HG的再加工形式或者降解片段,其可能是MIR205HG基因初始转录本的另一种剪接体形式。已知的lncRNA MIR205HG在肝癌组织和肝癌细胞中并不存在异常表达现象,这也是与LNC01309最大的区别。
LNC01309在肝癌细胞中主要定位于细胞质。由于lncRNA通常会作为内源性RNA海绵与miRNA作用并影响miRNA的表达与功能,而lncRNA定位于细胞质是这一功能发挥的前提条件,因此LNC01309理论上存在参与调控miRNA表达与功能的潜质,这也是后续应开展的研究方向。
综上,本研究团队依赖于肝癌组织的PacBio测序分析获得一条全新lncRNA LNC01309,在肝癌细胞中进行了相关功能探索,初步证明LNC01309通过与RBM38相互作用降低了RBM38的蛋白稳定性,促进肝癌细胞的增殖与迁移。
-
注:a,LNC01309与lncRNA MIR205HG序列相似性,并根据序列差异性设计特异性检测引物;b,利用荧光定量PCR分析lncRNA MIR205HG在肝癌组织与癌旁组织的表达差异;c,利用荧光定量PCR分析lncRNA MIR205HG在肝癌细胞和正常肝细胞中的表达差异;d、e,Hep 3B细胞经转染siRNA(20 nmol/L)处理24 h后,利用荧光定量PCR分别检测lncRNA MIR205HG和LNC01309的相对表达含量;f、g,Hep G2细胞中转染过表达质粒24 h后,利用荧光定量PCR分别检测LNC01309和miR-205的相对表达含量。
图 2 LNC01309与lncRNA MIR205HG的关系分析
表 1 利用NCBI BLAST比较LNC01309与lncRNA MIR205HG序列相似性
描述 种属 最高分 总分 查询覆盖值 E值 同源百分比 查询长度 登录号 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 3, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 992 NR_145434.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 1, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 895 NR_145437.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 2, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 857 NR_145433.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 4, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 940 NR_145435.1 Homo sapiens MIR205 host gene (MIR205HG), transcript variant 5, long non-coding RNA Homo sapiens 516 516 1 5×10-142 0.963 651 NR_145436.1 表 2 利用RBPmap和catRAPID在线预测能够与LNC01309结合的潜在靶标蛋白
蛋白质 序列位置 基序 K-mer Z评分 P值 CNOT4 129 gacaga gagaga 3.554 1.90×10-4 EIF4G2 50 cgccg cgccg 4.390 5.67×10-6 HNRNPA1 66 guaguagu guagcagc 3.361 3.88×10-4 MBNL1 70 gcgcagc cagcagc 3.903 4.75×10-5 RBM24 160 wgwgugd ggagugu 3.403 3.33×10-4 RBM38 163 kkguguk gugugcg 3.333 4.30×10-4 RBM6 37 ccacc ccacc 4.228 1.18×10-5 SRSF10 70 cagcag cagcag 4.513 3.20×10-6 -
[1] BHAN A, SOLEIMANI M, MANDAL SS. Long noncoding RNA and cancer: A new paradigm[J]. Cancer Res, 2017, 77(15): 3965-3981. DOI: 10.1158/0008-5472.CAN-16-2634. [2] PENG WX, KOIRALA P, MO YY. LncRNA-mediated regulation of cell signaling in cancer[J]. Oncogene, 2017, 36(41): 5661-5667. DOI: 10.1038/onc.2017.184. [3] MATTICK JS. The genetic signatures of noncoding RNAs[J]. PLoS Genet, 2009, 5(4): e1000459. DOI: 10.1371/journal.pgen.1000459. [4] KUMAR S, GONZALEZ EA, RAMESHWAR P, et al. Non-coding RNAs as mediators of epigenetic changes in malignancies[J]. Cancers (Basel), 2020, 12(12): 3657. DOI: 10.3390/cancers12123657. [5] MA L, BAJIC VB, ZHANG Z. On the classification of long non-coding RNAs[J]. RNA Biol, 2013, 10(6): 925-933. DOI: 10.4161/rna.24604. [6] MIYAMOTO M, MOTOOKA D, GOTOH K, et al. Performance comparison of second- and third-generation sequencers using a bacterial genome with two chromosomes[J]. BMC Genomics, 2014, 15: 699. DOI: 10.1186/1471-2164-15-699. [7] RHOADS A, AU KF. PacBio sequencing and its applications[J]. Genomics Proteomics Bioinformatics, 2015, 13(5): 278-289. DOI: 10.1016/j.gpb.2015.08.002. [8] XU J, ZHENG Y, PU S, et al. Third-generation sequencing found LncRNA associated with heat shock protein response to heat stress in Populus qiongdaoensis seedlings[J]. BMC Genomics, 2020, 21(1): 572. DOI: 10.1186/s12864-020-06979-z. [9] SHEN Y, LIANG W, LIN Y, et al. Single molecule real-time sequencing and RNA-seq unravel the role of long non-coding and circular RNA in the regulatory network during Nile tilapia (Oreochromis niloticus) infection with Streptococcus agalactiae[J]. Fish Shellfish Immunol, 2020, 104: 640-653. DOI: 10.1016/j.fsi.2020.06.015. [10] HUANG Z, ZHOU JK, PENG Y, et al. The role of long noncoding RNAs in hepatocellular carcinoma[J]. Mol Cancer, 2020, 19(1): 77. DOI: 10.1186/s12943-020-01188-4. [11] CHEN Y, LI Z, CHEN X, et al. Long non-coding RNAs: From disease code to drug role[J]. Acta Pharm Sin B, 2021, 11(2): 340-354. DOI: 10.1016/j.apsb.2020.10.001. [12] LI H, JI FX, XU XN, et al. Effects of long-chain non-coding RNA LINC00467 targeting miR-1247-5p on proliferation, migration and invasion of hepatocellular carcinoma cells[J]. Chin J Gerontol, 2021, 41(15): 3328-3333. DOI: 10.3969/j.issn.1005-9202.2021.15.048李豪, 姬发祥, 徐晓宁, 等. 长链非编码RNA LINC00467靶向miR-1247-5p调控肝癌细胞增殖、迁移及侵袭的影响[J]. 中国老年学杂志, 2021, 41(15): 3328-3333. DOI: 10.3969/j.issn.1005-9202.2021.15.048 [13] YUAN JH, YANG F, WANG F, et al. A long noncoding RNA activated by TGF-β promotes the invasion-metastasis cascade in hepatocellular carcinoma[J]. Cancer Cell, 2014, 25(5): 666-681. DOI: 10.1016/j.ccr.2014.03.010. [14] PANZITT K, TSCHERNATSCH MM, GUELLY C, et al. Characterization of HULC, a novel gene with striking up-regulation in hepatocellular carcinoma, as noncoding RNA[J]. Gastroenterology, 2007, 132(1): 330-342. DOI: 10.1053/j.gastro.2006.08.026. [15] DU Y, KONG G, YOU X, et al. Elevation of highly up-regulated in liver cancer (HULC) by hepatitis B virus X protein promotes hepatoma cell proliferation via down-regulating p18[J]. J Biol Chem, 2012, 287(31): 26302-26311. DOI: 10.1074/jbc.M112.342113. [16] GENG YJ, XIE SL, LI Q, et al. Large intervening non-coding RNA HOTAIR is associated with hepatocellular carcinoma progression[J]. J Int Med Res, 2011, 39(6): 2119-2128. DOI: 10.1177/147323001103900608. [17] LI G, ZHANG H, WAN X, et al. Long noncoding RNA plays a key role in metastasis and prognosis of hepatocellular carcinoma[J]. Biomed Res Int, 2014, 2014: 780521. DOI: 10.1155/2014/780521. [18] MAASS PG, LUFT FC, BÄHRING S. Long non-coding RNA in health and disease[J]. J Mol Med (Berl), 2014, 92(4): 337-346. DOI: 10.1007/s00109-014-1131-8. [19] MENG YM, JIANG J, WANG YH, et al. Mechanism of long noncoding RNA encoded peptide in tumor[J]. Clin Misdiagn Misther, 2021, 34(9): 113-116. https://www.cnki.com.cn/Article/CJFDTOTAL-LCWZ202109024.htm孟一妹, 蒋晶, 王禹涵, 等. 长链非编码RNA编码肽在肿瘤中的作用机制[J]. 临床误诊误治, 2021, 34(9): 113-116. https://www.cnki.com.cn/Article/CJFDTOTAL-LCWZ202109024.htm [20] SCHAUKOWITCH K, KIM TK. Emerging epigenetic mechanisms of long non-coding RNAs[J]. Neuroscience, 2014, 264: 25-38. DOI: 10.1016/j.neuroscience.2013.12.009. [21] HUANG JZ, CHEN M, CHEN, et al. A peptide encoded by a putative lncRNA HOXB-AS3 suppresses colon cancer growth[J]. Mol Cell, 2017, 68(1): 171-184. DOI: 10.1016/j.molcel.2017.09.015. [22] LI XX, SHI L, ZHOU XJ, et al. The role of c-Myc-RBM38 loop in the growth suppression in breast cancer[J]. J Exp Clin Cancer Res, 2017, 36(1): 49. DOI: 10.1186/s13046-017-0521-5. [23] JIANG Y, XU E, ZHANG J, et al. The Rbm38-p63 feedback loop is critical for tumor suppression and longevity[J]. Oncogene, 2018, 37(21): 2863-2872. DOI: 10.1038/s41388-018-0176-5. [24] CHO SJ, ZHANG J, CHEN X. RNPC1 modulates the RNA-binding activity of, and cooperates with, HuR to regulate p21 mRNA stability[J]. Nucleic Acids Res, 2010, 38(7): 2256-2267. DOI: 10.1093/nar/gkp1229. [25] XU E, ZHANG J, CHEN X. MDM2 expression is repressed by the RNA-binding protein RNPC1 via mRNA stability[J]. Oncogene, 2013, 32(17): 2169-2178. DOI: 10.1038/onc.2012.238. [26] YAN W, ZHANG J, ZHANG Y, et al. p73 expression is regulated by RNPC1, a target of the p53 family, via mRNA stability[J]. Mol Cell Biol, 2012, 32(13): 2336-2348. DOI: 10.1128/MCB.00215-12. [27] SUN Z, YANG GS, SIMA H, et al. Expression of RBM5 and its effect on prognosis of hepatocellular carcinoma patients after hepatectomy[J/CD]. Chin J Hepat Surg (Electronic Edition), 2021, 10(1): 93-97. DOI: 10.3877/cma.j.issn.2095-3232.2021.01.020.孙喆, 杨广顺, 司马辉, 等. RBM5在肝细胞癌中的表达及其对肝切除术后患者预后的影响[J/CD]. 中华肝脏外科手术学电子杂志, 2021, 10(1): 93-97. DOI: 10.3877/cma.j.issn.2095-3232.2021.01.020. [28] YE J, LIANG R, BAI T, et al. RBM38 plays a tumor-suppressor role via stabilizing the p53-mdm2 loop function in hepatocellular carcinoma[J]. J Exp Clin Cancer Res, 2018, 37(1): 212. DOI: 10.1186/s13046-018-0852-x. 期刊类型引用(0)
其他类型引用(1)
-