m6A修饰调控的LDB2抑制肺腺癌细胞增殖

翟东凤，王 鸽，李 蕾，贾小婷，郑国沛，尹 江

广州医科大学附属肿瘤医院肿瘤研究所//广州恶性肿瘤治疗转化医学重点实验室，广东 广州510095

肺癌是全球发病率和死亡率均居首位的恶性肿瘤。根据组织学分类，肺癌分为小细胞肺癌和非小细胞肺癌［1-2］。非小细胞肺癌中最常见的亚型为肺腺癌，约占肺癌总病例的40%［3］。与其他肺癌亚型相比，肺腺癌的发病率近年来明显增加［4］。肺腺癌由小的气道上皮，Ⅱ型肺泡细胞发展而来［2］。目前已鉴定了若干高频遗传变化（诸如EGFR、KRAS突变和ALK重排）与肺腺癌的发生发展相关，且针对这些靶点开发的药物能够有效抑制肺腺癌的进程［5-7］。近年来，在肺腺癌中针对程序性死亡配体1（PD-L1）和细胞毒性T 淋巴细胞相关抗原4（CTLA-4）的免疫治疗备受瞩目［8-9］。但是，还有更多的患者既无上述遗传突变，且对免疫检查点治疗不敏感。因此，探究肺腺癌发生和进展过程中新的关键分子，有望为临床肺腺癌提供潜在的诊断标志物或有效的治疗靶标。

LIM域结合蛋白2（LDB2）属于LIM域家族蛋白。LIM结构域是一种特殊的双锌指基序，存在于多种蛋白质中。LIM域介导功能复合物成员之间的特定接触并调节某些组成蛋白的活性［10］。目前公认的几种重要的LIM蛋白均与细胞骨架有关，在粘附斑块和肌动蛋白微丝组织中起作用［10］。LDB家族蛋白包括LDB1、LDB2和LDB3。LDB1 定位于细胞核，在细胞中组成型表达。LDB1在连接增强子和基因方面起着关键作用，它通过跨神经发育、心脏发生、视网膜发生和造血等多种发展途径形成具有细胞类型特异性的复合物［11］。LDB3特异性剪接会引起特发性扩张型心肌病［12］。许多研究表明，LDB2基因在视网膜发育和细胞周期中起调节作用［13］。LDB2通过与DLL4启动子区域结合调节该基因启动子活性，并且LDB2还介导VEGF诱导的内皮细胞中的DLL4表达，从而参与血管生成发芽和血管重塑等生理过程［14］。此外，LDB2还可驱动白细胞跨上皮迁移和动脉粥样硬化［13］。研究认为在乳腺癌中LDB2对ER受体活性有负调节作用，并破坏ERE依赖的转录［15］。在肝细胞癌中LDB2通过下调HEY1的表达抑制肿瘤细胞的增殖和迁移［16］。但是LDB2在肺腺癌中的作用尚未被阐明。

在本研究中，我们从多个数据库中预测LDB2在肺腺癌中低表达，且与患者不良预后负相关。在此基础上，本课题拟验证LDB2在肺腺癌中的表达模式，明确其在肺腺癌中的表达意义，探索其在肺腺癌中发挥的作用及潜在的分子机制。

1 资料和方法

1.1 细胞培养

正常肺上皮细胞16HBE 和肺腺癌细胞（H1299、H1975、A549、H1650和PC9）均由本实验室自主保存提供，上述细胞均在含有10%胎牛血清（Gibco）的1640 完全培养基中培养，将细胞置于37 ℃，5%培养箱中培养。

1.2 标本收集

收集2013年3月～2018年12月来本院就诊的肺腺癌患者的癌组织80例和癌旁组织17例，置于RNA保存液中，在液氮中保存备用。

1.3 转染

Lipofectamine 3000转染试剂购自Fermentas。将H1299细胞接种6孔板，次日进行脂质体转染。6孔板1孔的转染量如下：（1）质粒4 μg 和5 μL p3000 加入250 μL无血清培养基，（2）5 μL Lipofectamine 3000 加入250 μL无血清培养基，分别混匀，将（1）加入（2）混合液中，轻柔混匀后室温静置15 min，将上述混合液加入6孔板中，48 h后用于后续研究。

1.4 qRT-PCR

使用Omega RNA提取试剂盒来获得目标细胞及组织总RNAs，并将其逆转录成cDNAs。SYBR法用于检测各基因mRNA水平的表达量。各基因引物序列见表1。

1.5 Western blot检测

使用含有蛋白酶抑制剂的细胞裂解液处理目标细胞，BCA法测定蛋白浓度。蛋白样品上样量是30 μg，蛋白样品经SDS-PAGE电泳后，转至PVDF膜上。5%BSA室温封闭2 h，一抗孵育过夜，HRP偶联二抗室温孵育2 h，ECL底物检测蛋白表达量。

1.6 细胞计数

将5×103细胞种于12孔板中，混匀后置于37 ℃，5%培养箱中培养，此后连续7 d用胰酶消化不同孔的细胞，混匀后吸出同等体积的细胞混合液，使用细胞计数仪计算细胞数目。

1.7 软琼脂克隆

在6孔板中均匀铺入含有0.75%琼脂糖的完全培养基的下层胶，室温下冷却后，将2×103细胞加入含有0.375%琼脂糖的完全培养基中，混匀后铺入凝固的下层胶上，室温下冷却后置于37 ℃，5%培养箱中培养2周，显微镜下观察克隆大小并拍照。

1.8 流式细胞检测

将5×105细胞均匀铺于六孔板中，使用含10%胎牛血清的1640完全培养基培养，并置于37 ℃，5%培养箱中。24 h后将完全培养基置换成无血清的1640进行饥饿处理24 h，随后换成完全培养基培养，按照说明书对处理后的细胞进行固定及染色，最后使用流式细胞仪分析细胞周期。

1.9 MeRIP-qPCR实验

购买Millipore 的Magna RIPTMRNA-Binding Protein Immunoprecipitation Kit，按照说明书开展实验。RIP裂解液处理细胞后，使用anti-m6A抗体（RIPqPCR实验中用到的是YTHDC2抗体）孵育细胞裂解液，磁珠与上述混合物孵育，蛋白酶K处理相关蛋白，纯化RNA并逆转录成cDNA，qRT-PCR检测目标RNA的含量。

1.10 双荧光素酶报告基因实验

购买TaKaRa的定点突变试剂盒MutanBest Kit，将LDB2转录本上“AGA（274 bp）CT”序列缺失，将野生型的序列（转录本1～400 bp）和缺失型的序列构建到报告基因载体pGL4 Basic，分别命名为“LDB2-WT”和“LDB2-DEL”。分别将野生型报告基因/缺失型报告基因与pRLTK（海肾荧光素酶报告基因载体，作为内参）、YTHDC2过表达质粒共转染H1299细胞。48 h后，1×passive lysis buffer裂解各组细胞30 min，将裂解液转移至白板，依次加入萤火虫荧光素酶底物和海肾荧光素酶底物（含有荧光猝灭剂），多功能酶标仪检测各孔化学发光值，计算各孔荧光素酶活性。

1.11 生物信息学预测网址

Prognoscan(http://dna00.bio.kyutech.ac.jp/PrognoScan/)

GEPIA(http://gepia.cancer-pku.cn/)

cBioPortal for Cancer Genomics (https://www.cbioportal.org/)

TCGA(http://www.tcga.org/)

SRAMP(http://www.cuilab.cn/sramp)

1.12 统计处理

所有实验重复3次以上，数据均用均数±标准差表示。两对基因表达的相关性通过SPSS 19.0软件包进行Pearson相关性分析。两组间比较采用独立样本t检验和卡方检验，P＜0.05认为差异有统计学意义。

2 结果

2.1 LDB2在肺腺癌中的表达模式

下载GEO、TCGA及CPTAC 3个数据库的肺腺癌中下调的基因表达信息，取交集分析后发现肺腺癌中有25 个表达下调的基因（图1A）。GEPIA 数据库显示LDB2在483例肺腺癌组织中的表达量远低于其在347例肺正常组织中的表达量（图1B）。GSE10072 和GSE43767数据集均展示LDB2在肺腺癌组织中低表达；而在GSE32863数据集中发现LDB2在58例肺腺癌中的表达量显著低于其在配对癌旁组织中的表达量（图1C）。CPTAC 数据库证实LDB2 蛋白在111 例肺腺癌组织中的表达量低于其在111例正常组织中的表达量（图1D）。我们通过免疫组化方法检测LDB2在80例肺腺癌组织和17例正常组织的表达量，亦获得相似的结果（图1E）。最后，qRT-PCR方法显示LDB2在肺腺癌细胞株A549、H1299、H1975、H1650和PC9中的表达量均显著低于其在正常肺上皮细胞16HBE中的表达量（图1F）。

2.2 LDB2在肺腺癌中低表达的临床意义

在GEPIA数据库中统计了273例LDB2高表达肺腺癌患者和144例LDB2低表达的肺腺癌患者的总生存率，结果发现LDB2高表达与肺腺癌患者较好的预后正相关（图2A，P=0.049）。在线数据库Prognoscan 预测到LDB2在多个数据集中与肺腺癌患者不良预后负相关。最后，将收集到的80 例肺腺癌组织分为51 例LDB2低表达的肺腺癌组织和29例高表达LDB2的肺腺癌组织，经统计分析：随着LDB2表达量越低，肺腺癌患者总的生存率越差，LDB2表达量与肺腺癌患者总的生存率表达负相关（P=0.0087，图2C）。

2.3 LDB2抑制肺腺癌细胞增殖

H1299中过表达LDB2，通过qRT-PCR及Western blot验证，H1299细胞株中过表达LDB2构建成功（图3A）。体外细胞计数实验结果显示：与对照组对比，过表达LDB2 可明显抑制H1299 细胞的增殖能力（图3B）。软琼脂克隆实验发现：与对照组相比，过表达LDB2组可显著抑制H1299细胞的克隆形成能力（图3C）。流式细胞术证实：与对照组相比，过表达LDB2后将诱导肺腺癌细胞发生S期阻滞（图3D）。

2.4 m6A修饰调控LDB2表达

我们首先通过在线数据库SRAMP预测到LDB2转录本上存在多个m6A结合位点，其中274位点的m6A具有高度可靠性，后期重点验证该位点（图4A）。在H1299 细胞中开展MeRIP-qPCR 实验，发现与IgG 相比，其m6A组的LDB2含量显著升高（图4B）。GEPIA数据库显示YTHDC2在肺腺癌中低表达，且与LDB2表达正相关（图4C，r=0.22，P=1.2e-06）。在H1299细胞中过表达YTHDC2，qRT-PCR和WB检测结果均证实过表达细胞系构建成功（图4D）。通过qRT-PCR 和Western blot检测发现H1299细胞中过表达YTHDC2可增加LDB2的表达（图4E）。H1299过表达YTHDC2的细胞中开展RIP 实验，定量结果显示LDB2 在YTHDC2 组的含量是其在IgG 组含量的19.35 倍（图4F）。双荧光素酶报告基因试验发现：与对照组相比，YTHDC2可上调野生型报告基因活性13.04倍；而当该位点m6A发生缺失后，YTHDC2对其活性调控也大幅度降低（图4G）。

图1 LDB2 在肺腺癌中的表达模式Fig.1 Expression pattern of LDB2 in LUAD.A:TCGA,GEO and CPTAC databases were searched for genes with low expressions in LUAD. B: Expression levels of LDB2 in LUAD and normal tissues were evaluated by GEPIA databank, **P＜0.01. C: GEO database analysis of LDB2 expression in LUAD and normal tissues, **P＜0.01. D: CPTAC databank analysis of the protein expression levels of LDB2 in LUAD and normal tissues, **P＜0.01. E: Immunohistochemical method was used to detect LDB2 expression in 80 LUAD tissues and 17 normal lung tissues,**P＜0.01.F:qRT-PCR for to estimating LDB2 expression in 16HBE and several LUAD cells,**P＜0.01.

图3 LDB2 抑制肺腺癌细胞的增殖能力Fig.3 LDB2 inhibits proliferation ability of LUAD cells.A:qRT-PCR and Western blot were performed to verify the expression level of LDB2 in H1299 cells. **P＜0.01 vs H1299-con. (B) Cell counting and (C) soft agar clone formation assays were used to analyze the cell proliferation ability.**P＜0.01 vs H1299-con.D:Cell cycle was detected by flow cytometry assay in H1299 cells.

3 讨论

尽管在研究肺腺癌的相关发病机制及临床治疗上取得了新的进展，但是肺腺癌仍然是最具侵袭性及快速致死性的肿瘤类型之一，总生存率不到5年［3］。因此探讨肺腺癌发生发展的分子机制是目前亟待解决的重大课题。在本研究中，我们前期通过数据库预测、组织学及细胞学水平验证，发现LDB2在肺腺癌中显著低表达。并且LDB2的表达与肺腺癌患者良好的临床预后呈正相关。细胞功能实验发现，过表达LDB2后可抑制肺腺癌细胞的增殖能力。深入探究后发现m6A修饰调控LDB2在肺腺癌中低表达。

m6A修饰是真核mRNA最常见的内部修饰，广泛参与许多细胞过程，例如RNA转录、剪接、核转运、降解和翻译［17］。大量研究发现m6A 扮演了肿瘤启动子或抑制子的角色，通过调节肿瘤相关基因的表达，而参与肿瘤细胞的启动、增殖、分化、转移和代谢重编程等生命过程［17］。m6A修饰主要发生在DRACH基序上（D对应于A，G或U；R对应于G或A；H对应于A，C或U），通常在转录本的5'-UTR、3'-UTR 和终止密码子附近［18］。m6A的效应器包括编码器、读码器和消码器。编码器（如METTL3、METTL14和WTAP等）通过S-腺苷甲硫氨酸（SAM）转移酶将甲基添加到目标RNA 上，形成m6A 甲基化位点［19］。FTO 和ALKBH5 是两种主要的m6A销码器，它们以Fe(II)/α-酮戊二酸依赖性方式催化去除m6A修饰［20］。读码器（如YTHDF1-3、YTHDC1-2和HNRNPA2B1等）能够专门识别并结合带有m6A修饰的RNA，以控制其命运并调节下游功能［21］。

YTHDC2的YTH结构域优先通过保守的疏水口袋与含m6A的RNA结合，而锚蛋白重复序列介导与5'-3'外核糖核酸酶XRN1的RNA相互作用。YTHDC2与小核糖体亚基相互作用，接近mRNA进入/退出位点，从而调控特定基因的表达［22］。目前关于YTHDC2对转录本的调控尚未有定论。有报道指出YTHDC2能够结合到SLC7A11转录本上促进其降解［23］。也有文章发现YTHDC2在耐辐射的鼻咽癌细胞中高表达。生物信息学和机制研究发现，YTHDC2能结合到胰岛素样生长因子1受体（IGF1R）mRNA上，促进该转录本的翻译起始，进而激活IGF1R-AKT/S6信号通路，从而实现鼻咽癌放疗抵抗［24］。YTHDC2在肺腺癌中低表达，且与肺腺癌的低分化、淋巴结转移、肿瘤大小和临床分期等特征相关。YTHDC2能够抑制肺腺癌细胞的增殖和迁移能力［25］。与此不同的是，我们在研究中发现YTHDC2可以调控LDB2的mRNA丰度，也调控其蛋白丰度，我们优先考虑YTHDC2 对LDB2 转录活性有影响。由于LDB2的转录起始位点位于该转录本161位，274位点与转录起始位点比较接近，仍可作为启动子区域，因此构建启动子报告基因载体，双荧光素酶报告基因证实YTHDC2可识别LDB2上的m6A位点，从而转录活化该转录本。与这些数据相吻合的是，GEPIA数据库数据亦显示YTHDC2与LDB2在肺腺癌中表达正相关。在后续研究中，我们将进一步寻找LDB2转录本上m6A修饰的编码器和销码器。

图4 YTHDC2通过m6A修饰上调LDB2的表达Fig.4 M6A modifies the expression of LDB2.A: The distribution of m6A sites on LDB2 RNA was predicted via online software SRAMP;B:MeRIP-qPCR assay was carried out to examine m6A sites on LDB2,**P＜0.01;C:Expression pattern of YTHDC2 in LUAD.Pearson correlation analysis was performed to discover the expression relationship between LDB2 and YTHDC2 in LUAD; YTHDC2 overexpression plasmid was transfected in H1299 cells, qRT-PCR and Western blot assays were implemented to detect the expression of YTHDC2(D)and LDB2(E),**P＜0.01;F:We performed RIP-qPCR assays to evaluate whether YTHDC2 could bind to the transcript of LDB2,**P＜0.01;G:Dual-luciferase report assay in 293T cells with wildtype or deletion reporter plasmid,**P＜0.01.

因此，本研究阐述了LDB2在肺腺癌中表达量显著下降，其表达受m6A修饰调控抑制，且LDB2可诱导肺腺癌细胞发生S期阻滞，抑制其增殖能力。LDB2调控肺腺癌细胞周期的分子机制有待于进一步研究。