《Non-Coding RNA》：Integrated Transcriptomic, In Silico, and In Vitro Characterization of lncRNA ENST00000615487.1 Reveals Epithelial-Specific Expression, Differential Subcellular Distribution Between Normal and Colorectal Cancer Cells, and Potential Regulatory Functions
编辑推荐：
背景/目的：长链非编码RNA（lncRNAs）通过其与DNA、蛋白质及非编码RNA的相互作用，是肿瘤生物学的重要调控因子。尽管ENST00000615487.1（亦称CTD-2396E7.11/AC010503.4）已被发现与多种恶性肿瘤相关，但其在结直肠癌（
背景/目的：长链非编码RNA（lncRNAs）通过其与DNA、蛋白质及非编码RNA的相互作用，是肿瘤生物学的重要调控因子。尽管ENST00000615487.1（亦称CTD-2396E7.11/AC010503.4）已被发现与多种恶性肿瘤相关，但其在结直肠癌（CRC）中的生物学角色仍未被充分阐明。本研究旨在通过整合体外（in vitro）与计算机模拟（in silico）方法，探究ENST00000615487.1在CRC中的表达模式、细胞及亚细胞定位以及潜在功能角色。方法：研究人员利用CPC2和RNA Analyzer 3工具获取该转录本的分子特征。通过UCSC Xena浏览器分析ENST00000615487.1在10种肿瘤类型中的差异表达。采用PCR在正常、肿瘤及成纤维细胞结肠细胞系中实验评估转录本表达，同时通过lncATLAS、lncLocator及iLoc-LncRNA工具预测亚细胞定位，并经qRT-PCR实验验证。分析来自GSE161277数据集的单细胞RNA测序数据以确定细胞类型特异性表达模式。利用Fasim-LongTarget、miRDB及AnnoLnc2探究其与DNA、miRNA及蛋白质的潜在相互作用，随后通过STRING和Enrichr进行功能富集分析。结果：ENST00000615487.1被鉴定为一个结构稳定的非编码转录本，具有高度有序的二级结构。差异表达分析显示，与正常结肠组织相比，其在CRC中显著下调。单细胞转录组分析揭示其主要呈上皮特异性表达。计算机模拟与实验分析表明，在正常结肠细胞中主要定位于细胞核，而CRC细胞中则出现细胞质富集。功能分析鉴定出其与HIP1R、RPH3AL、特定miRNA及参与转录调控和RNA加工通路的蛋白质存在潜在相互作用，同时与参与囊泡运输的蛋白质存在功能联系。结论：ENST00000615487.1是一个结构稳定的lncRNA，在CRC中呈现环境依赖性的表达和定位模式，提示在结直肠癌变过程中，其可能从核内转录调控向细胞质转录后调控功能发生潜在转变。
**论文解读文章**
**研究背景与问题**
结直肠癌（colorectal cancer, CRC）是全球第三大常见恶性肿瘤（占9.6%），死亡率位居第二（9.3%）。2022年全球新发病例超过190万，死亡约90.4万例，每10例癌症相关死亡中即有1例由CRC导致。尽管外科手术、化疗、靶向治疗及免疫治疗取得显著进展，但晚期CRC仍具有高发病率和死亡率。早期检测可显著改善患者生存，但现有诊断和预后方法在敏感性、特异性及疾病进展预测方面仍存在重要局限。因此，当前CRC研究日益聚焦于识别新型分子标志物及参与结直肠癌变的调控机制，以期改善诊断、预后评估及治疗策略。长链非编码RNA（long non-coding RNAs, lncRNAs）作为癌症相关生物过程的重要调控因子备受关注。ENST00000615487.1（亦称CTD-2396E7.11/AC010503.4）已被发现与肾癌、卵巢癌、乳腺癌、子宫内膜癌、膀胱癌、胃癌及黑色素瘤等多种恶性肿瘤相关，并参与代谢重编程、铁死亡及自噬等关键过程，主要发挥保护性作用，是较好的预后因子。在CRC中，该转录本被定义为与疾病进展相关的四个关键“枢纽”lncRNAs之一，其在正常组织中表达最高，从腺瘤向癌变过程中逐渐降低。然而，研究人员前期基于球状体培养的下一代测序（Next-Generation Sequencing, NGS）结果却显示CRC细胞系中该转录本高表达。这些矛盾结果以及其功能研究不足，成为开展本研究的核心动机。本研究旨在通过整合体外（in vitro）与计算机模拟（in silico）分析，系统表征ENST00000615487.1在CRC及其他恶性肿瘤中的表达谱、细胞及亚细胞定位，并生成关于其潜在调控功能的假说。论文发表在《Non-Coding RNA》。
**关键技术与方法**（不超过250字）
1. **结构及编码潜能分析**：采用Coding Potential Calculator 2（CPC2）和RNA Analyzer 3评估转录本结构稳定性及编码潜能。
2. **差异表达分析**：利用UCSC Xena浏览器（TCGA Pan-Cancer队列，TOIL RSEM FPKM数据集）分析10种肿瘤类型中肿瘤与正常组织的差异表达。
3. **单细胞转录组分析**：从GSE161277数据集获取正常（11,817个细胞）、腺瘤（11,618个细胞）、癌（12,395个细胞）及血液（2,134个细胞）样本的单细胞RNA测序数据，采用Seurat包进行聚类和细胞类型注释（基于CellMarker 2.0、PanglaoDB及文献标志物）。
4. **亚细胞定位预测与验证**：运用lncATLAS、lncLocator及iLoc-LncRNA 2.0进行计算机预测；通过细胞核-细胞质分离及定量实时PCR（qRT-PCR）在HCEC-1CT、HCT116、DLD1、SW620、CCD-18Co细胞系中实验验证。
5. **分子相互作用预测**：利用Fasim-LongTarget预测DNA结合位点；miRDB预测miRNA相互作用；AnnoLnc2预测蛋白质相互作用及共表达网络；STRING和Enrichr进行功能富集分析。
6. **统计方法**：差异表达采用Welch’s t检验并经Benjamini-Krieger-Yekutieli两步法FDR校正；相关性采用Spearman检验；单细胞比较采用Wilcoxon秩和检验，p<0.05为显著。
**研究结果**
**2.1 ENST00000615487.1转录本的特征**
通过CPC2和RNA Analyzer 3分析，ENST00000615487.1被鉴定为无编码潜能的非编码RNA，由单一外显子组成，长688个碱基对（bp），具有高度有序的二级结构，结构稳定。
**2.2 ENST00000615487.1的表达**
**2.2.1 计算机模拟分析肿瘤组织中的差异表达**
利用UCSC Xena浏览器分析10种肿瘤类型中肿瘤与正常组织的差异表达。在甲状腺癌（THCA）、胆管癌（CHOL）、前列腺癌（PRAD）和肾乳头状细胞癌（KIRP）中，肿瘤组织显著上调（FDR校正q<0.001）；在结肠腺癌（COAD）、子宫内膜癌（UCEC）和肾透明细胞癌（KIRC）中显著下调（q<0.001）；在胰腺癌（PAAD）、直肠腺癌（READ）和肾嫌色细胞癌（KICH）中无显著差异。
**2.2.2 正常、肿瘤及成纤维细胞结肠细胞中ENST00000615487.1的检测**
PCR检测显示，在成纤维细胞系CCD-18Co及CRC细胞系HCT116、DLD1、SW620中检测到表达，但在正常结肠上皮细胞系HCEC-1CT中未检测到。
**2.2.3 结直肠癌中ENST00000615487.1的单细胞转录组分析**
单细胞RNA测序分析鉴定出免疫、基质及上皮细胞群。ENST00000615487.1在4.48%的正常细胞、6.14%的腺瘤细胞和8.33%的癌细胞中表达，主要局限于上皮细胞群（肠细胞样、杯状细胞样、黏液分泌型、侵袭性、应激相关及活化肿瘤上皮状态），在免疫和基质细胞中表达极低。癌样本中表达阳性细胞比例更高，但阳性正常细胞平均表达强度更高（p<0.001），提示癌变过程中转录调控发生动态重塑。
**2.3 ENST00000615487.1亚细胞定位的确定**
**2.3.1 计算机模拟预测**
lncATLAS分析显示，在多数细胞系中ENST00000615487.1的细胞质-细胞核相对浓度指数（CN RCI）为正值，提示主要定位于细胞质；仅正常人类表皮角质形成细胞（NHEK）中为负值（核定位）。lncLocator预测细胞质为最可能定位（概率0.335），iLoc-LncRNA 2.0预测核糖体相关定位（概率0.862）。
**2.3.2 正常、肿瘤及成纤维细胞结肠细胞亚细胞及胞外部分的实验验证**
qRT-PCR显示，在CRC细胞系（HCT116、DLD1、SW620）及成纤维细胞系CCD-18Co中，ENST00000615487.1主要分布于细胞质；而在正常结肠上皮细胞系HCEC-1CT中主要分布于细胞核。半定量PCR在HCEC-1CT、HCT116及CCD-18Co的外泌体RNA中均检测到该转录本，证实其存在于细胞外囊泡。
**2.4 ENST00000615487.1功能的计算机模拟预测**
**2.4.1 潜在DNA相互作用预测**
Fasim-LongTarget分析鉴定出16个与ENST00000615487.1相互作用评分>500的基因，其中HIP1R（评分1512）和RPH3AL（评分1055）最强。CRC组织中ENST00000615487.1与HIP1R（Spearman r=0.2964，p<0.0001）及RPH3AL（r=0.4748，p<0.0001）呈正相关。Enrichr功能富集分析显示这些基因参与突触囊泡胞吐、蛋白质输出及细胞内定位等过程。
**2.4.2 潜在miRNA相互作用预测**
miRDB鉴定出4个与ENST00000615487.1相互作用评分>70的miRNA：hsa-miR-4660（评分87）、hsa-miR-1295b-3p（86）、hsa-miR-4470（76）和hsa-miR-6782-5p（73）。预测结合位点位于转录本单链环状区域。进一步分析发现hsa-miR-4660靶向ZNF48，hsa-miR-1295b-3p靶向TMPRSS3，且这两个基因的表达与ENST00000615487.1正相关。
**2.4.3 潜在蛋白质关联预测**
AnnoLnc2鉴定出与ENST00000615487.1预测相互作用的蛋白质，STRING富集分析显示这些蛋白质主要参与RNA聚合酶II转录终止、复制叉逆转及SCF复合体组装等过程，涉及非同源末端连接、转录、RNA降解及核糖体生成通路。共表达分析显示与ENST00000615487.1正相关的蛋白质富集于囊泡装载、蛋白质分泌及内质网-高尔基体囊泡运输等分泌通路。
**讨论与结论**
研究表明ENST00000615487.1是一个结构稳定的非编码转录本，在CRC中表达下调，主要在上皮细胞中表达，并在正常与癌变细胞间呈现差异亚细胞分布（正常细胞核定位 vs. 肿瘤细胞质定位），提示其功能可能从核内转录调控向细胞质转录后调控发生转变。计算机预测提示其可能通过DNA相互作用（如与HIP1R、RPH3AL）参与囊泡运输和蛋白质定位，通过miRNA海绵作用（如hsa-miR-4660、hsa-miR-1295b-3p）调控靶基因，并与参与转录、复制及DNA修复的蛋白质相互作用。然而，这些发现基于计算预测和相关性分析，需通过功能实验验证。研究结论：本研究鉴定ENST00000615487.1为与CRC生物学相关的潜在重要lncRNA。编码潜能分析支持其非编码转录本注释；表达分析显示其在多种恶性肿瘤中表达失调，包括CRC中下调。单细胞转录组分析进一步证明其表达主要局限于上皮细胞群，特别是肠细胞样、杯状细胞样、黏液分泌型及肿瘤相关上皮状态，在免疫和基质细胞中表达极低。此外，癌样本中ENST00000615487.1阳性细胞比例更高，提示结直肠癌变过程中转录调控发生改变。该转录本在外泌体中的检测表明其可能被释放到细胞外环境，值得进一步探究其在细胞间通讯及肿瘤微环境调控中的潜在作用。计算机预测与亚细胞分离分析均提示ENST00000615487.1在HCEC-1CT细胞中主要呈核定位，而在CRC细胞系中主要呈细胞质定位，表明其定位在正常与恶性细胞环境中可能不同，并与环境依赖功能相关。计算分析鉴定出候选DNA、蛋白质及miRNA相关网络，重点涉及囊泡运输、蛋白质定位、转录调控、DNA复制、DNA修复及转录后调控等通路。本研究为ENST00000615487.1的表达谱、细胞及亚细胞分布以及预测分子关联提供了全面表征，为后续功能研究奠定基础。
