基于全转录组分析的Docetaxel耐药前列腺癌细胞LNCaP与C4-2B揭示可变剪接与非编码RNA维度的Taxane耐药机制

《Cells》：Whole-Transcriptome Analysis of Docetaxel-Resistant LNCaP and C4-2B Prostate Cancer Cells Reveals Alternative Splicing and Non-Coding RNA Dimensions of Taxane Resistance

Docetaxel（DTX）是转移性去势抵抗性前列腺癌的一线化疗药物，但获得性耐药几乎普遍存在。研究人员此前构建了同基因的DTX耐药LNCaP与C4-2B细胞系衍生株，并使用Affymetrix Clariom D芯片对其进行了分析，但初步分析仅限于基因水平的差异表达，突出了ABCB1。在本研究中，研究人员采用线性模型框架重新分析了全部十二张芯片（两个细胞谱系×亲本/耐药×三个传代），以充分利用外显子/连接水平及非编码转录组数据。研究人员确认了核心耐药特征（ABCB1、VTRNA1-1、TRAIL/TNFSF10），并鉴定出一个共享的587基因耐药程序，其中丝分裂与微管运动相关基因被协同下调。剪接分析揭示了315个高置信度的交互外显子使用事件（包括NTNG1、CTNND1、MYO6、MAT2A）。按生物型解析的分析鉴定出123个耐药相关lncRNA（包括LINC00152/CYTOR），并发现共调控的反义/正义配对（如MME-AS1/MME、ATP1B3-AS1/ATP1B3），提示候选的局部（正义–反义）调控关系。增殖程序与生长动力学一致，因为细胞周期/丝分裂基因在生长较慢的C4-2BR细胞中被下调。通过充分利用Clariom D平台，研究人员将DTX耐药的分子图谱从基因表达扩展到协同的可变剪接与结构化的非编码应答。

研究背景方面，前列腺癌是全球男性中最常诊断的恶性肿瘤之一，也是男性癌症死亡的主要原因。局限性疾病通常可治愈，但晚期肿瘤会进展为转移性去势抵抗性前列腺癌（mCRPC），系统性治疗无法治愈。Taxane类药物Docetaxel（DTX）是首个被证明能延长mCRPC生存期的药物，近期研究还显示其在转移性激素敏感性疾病中联合雄激素剥夺治疗可改善结局，使其成为前列腺癌治疗全程中使用最广泛的药物之一。然而其获益有限且短暂：应答为部分性，获得性耐药几乎普遍存在，尚无经验证的生物标志物可预测Taxane获益，Taxane治疗失败后的选择有限。因此理解前列腺癌细胞如何产生DTX耐药具有直接的临床重要性。已确立的Taxane耐药机制包括ATP结合盒转运蛋白ABCB1/MDR1介导的药物外排增加、微管动力学与微管蛋白同型使用改变，以及生存与雄激素受体（AR）信号重编程。前列腺癌还提供了由信使RNA加工而非仅基因表达驱动的耐药的最清晰范例：组成性活跃的AR剪接变体AR-V7由AR转录本可变剪接产生，赋予对恩杂鲁胺和阿比特龙等AR靶向药物的耐药性，并已成为临床可操作的生物标志物。然而在AR位点之外，全转录组可变剪接对Taxane耐药的贡献在很大程度上尚未被绘制。可变剪接在癌症中普遍失调，并日益与进展和药物应答相关联，但DTX耐药模型的系统性外显子与连接水平分析基本缺失。非编码转录组同样存在盲区：长链非编码RNA（lncRNA）在前列腺癌生物学中地位突出，从诊断标志物PCA3到促进肿瘤侵袭性的SCHLAP1，lncRNA调控AR信号及对AR靶向治疗的耐药性；LINC00152/CYTOR等lncRNA进一步影响多种肿瘤类型的增殖与化疗耐药。尽管如此，按生物型解析、全基因组范围内lncRNA及其他非编码物种如何随DTX耐药变化的研究仍然稀少。

研究人员开展了以下研究：对同一Clariom D数据集进行重新分析，以表征DTX耐药的差异表达、可变剪接和按生物型解析的非编码RNA图谱，并通过lncRNA–编码基因调控配对将编码与非编码应答联系起来。研究得出以下结论：DTX耐药伴随着协调的可变剪接程序与结构化的lncRNA应答，以及局部共调控的反义配对，这些是耐药状态的关键特征；研究不仅验证了原始研究的核心耐药特征，还揭示了可能支撑前列腺癌细胞在Taxane压力下适应性生存的新转录调控层次，为实验验证和临床转化提供了优先路线图。

关键技术方法方面（不超过250字）：研究人员使用同基因DTX耐药LNCaP^R与C4-2B^R及其亲本对照，每种条件三个独立传代，共十二张Affymetrix Clariom D Human芯片（原始数据来自Lima等，GEO登录号GSE345127）。采用SST-RMA算法生成基因水平与探针组水平摘要。差异表达与差异剪接使用R语言limma包（含diffSplice），并以Affymetrix转录组分析控制台（TAC）3.1进行独立交叉验证。生物型注释、GO/KEGG/Reactome/GenMAPP富集、反义配对与共调控相关分析在R 4.3.3中完成。细胞增殖采用无透镜Cellwatcher M显微镜实时监测汇合度。候选lncRNA以逆转录定量PCR（RT-qPCR）验证，以GAPDH为内参。

研究结果方面：3.1 研究设计与数据质量——通过主成分分析（PCA）显示第一主成分（66%方差）区分两个谱系，第二主成分（13%）在两个谱系内区分耐药与亲本细胞，三个传代紧密聚类，表明DTX耐药与分析群体内一致的转录状态相关。3.2 差异表达确认并扩展已知耐药特征——采用组均值设计的limma分析，在FDR<0.05且|log2FC|>1下，共享程序包含587个基因（294上调、293下调），C4-2B^R与LNCaP^R分别有812和254个基因变化；特征由ABCB1（log2FC +3.9）、穹窿RNA VTRNA1-1（+4.2）和TNFSF10/TRAIL（+4.8）领衔，并包括SLC16A1/MCT1（+1.9）、IGFBP3（+3.2）、CLU（+2.3）和STMN1（?1.3）。3.3 DTX耐药与广泛差异外显子/连接使用相关——diffSplice在37,565个多探针组基因的600,356个检测到的PSR/JUC探针组中，鉴定出共享、C4-2B和LNCaP对比中分别为1507、1029和1976个差异剪接基因；经保守过滤排除表达驱动伪影后，获得274个不同基因中的315个高置信度交互外显子使用事件，领先候选包括NTNG1、ADGRL2、CTNND1、MYO6、MAT2A、KIAA1217和UBR5；TAC交叉验证显示75%的C4-2B和65%的LNCaP diffSplice显著基因也被TAC标记，315个高置信度事件中89%被TAC重现。3.4 通路富集——下调基因强烈富集于细胞分裂、丝分裂核分裂、染色体分离和微管运动，上调基因富集于凋亡相关条目；差异剪接基因富集于细胞周期、跨膜转运、蛋白磷酸化及mRNA加工功能。3.5 非编码转录组对耐药产生应答——芯片检测到24,867个lncRNA和1546个小RNA位点；共享对比中123个lncRNA随耐药变化（C4-2B 158个、LNCaP 62个），以LINC00152/CYTOR领衔，另有LINC01029、LINC01239和反义RNA MMP24-AS1；小RNA大多稳定，仅U2剪接体snRNA轻度增加（log2FC +1.0）。3.6 反义lncRNA与其正义基因共调控——在176个命名的反义lncRNA/正义基因配对中，共调控显著强于随机配对（中位|r| 0.33 vs 0.24；KS p=4×10^?4），且对表达匹配与DE状态匹配的零分布均稳健；60对一致、35对不一致，最强不一致配对为ATP1B3-AS1/ATP1B3（r=?0.84），一致配对包括MME-AS1/MME和SLC16A1-AS1/SLC16A1。3.7 转录增殖程序追踪活细胞生长动力学——实时汇合度监测显示C4-2B^R生长明显慢于亲本，而LNCaP^R增殖速率相当；预定义的31基因细胞周期/丝分裂特征在C4-2B^R显著降低（增殖指数?0.90，p=0.010；30/31基因下调），在LNCaP^R仅轻度非显著降低（?0.50，p=0.28）。3.8 候选转录本的实验验证——RT-qPCR显示CYTOR/LINC00152在耐药亚系中上调（LNCaP约4.2倍、C4-2B约2.3倍），TMEM132D-AS1同样上调（LNCaP约4.8倍、C4-2B约4.6倍），与其正义基因TMEM132D既无差异表达也无差异剪接，表明其诱导独立于正义基因。

讨论部分总结：研究人员指出，DTX获得性耐药是mCRPC近乎普遍的临床挑战，其分子基础远超ABCB1介导的药物外排。通过重新分析此前仅在基因水平 interrogated 的全转录组Clariom D数据集，研究揭示了同基因耐药模型中多层转录程序，确认核心耐药特征的同时揭示了三个此前未探索的维度：协调的可变剪接、结构化的lncRNA应答以及局部共调控的反义lncRNA–编码基因配对。广泛差异外显子/连接使用与候选可变剪接重塑一致，涉及细胞骨架动力学、黏附和代谢相关基因；MYO6的LI⁺同型转换和KDM5D的N端缺陷同型为谱系限制性事件，均由结构域架构预测而非证实，需同型解析实验确定因果性。非编码转录组中LINC00152/CYTOR的上调与其在增殖、化疗耐药及促进AR-V7剪接中的作用一致；核心剪接体snRNA池基本稳定，不支持snRNA水平整体改变，U2 snRNA的轻度增加可能原则上偏向外显子使用。反义配对的显著共调控提示顺式作用机制（转录干扰、染色质重塑或RNA–RNA双链形成）可能协调邻近编码与非编码转录本，但不建立因果性。谱系特异性增殖适应显示C4-2B^R进入慢周期状态并伴随丝分裂基因协同抑制，而LNCaP^R基本维持增殖，提示ABCB1介导的外排是共享核心机制而增殖表型反映谱系特异性适应。临床意义包括lncRNA特征（尤其LINC00152上调）作为DTX耐药候选生物标志物、剪接与反义lncRNA程序作为潜在治疗靶点，以及按分子亚型定制的联合治疗。局限性包括单一数据集、两个细胞系、每个谱系三个传代而非独立衍生的耐药克隆、微阵列剪接推断相对RNA-seq的固有局限，以及反义关系的相关性本质。

研究结论翻译：总之，该分析通过揭示协调的基因表达、可变剪接、结构化的lncRNA应答以及局部共调控的反义配对作为耐药状态的关键特征，扩展了对前列腺癌DTX耐药的分子理解。这些发现不仅验证了原始研究的核心耐药特征，还揭示了可能支撑前列腺癌细胞在Taxane压力下适应性生存的新转录调控层次。尽管研究存在局限，结果为实验验证和临床转化提供了优先路线图，提名了候选靶点与调控通路，这些一旦验证，可能为未来监测、预测或调控转移性前列腺癌DTX耐药的策略提供基础。