TUBB2A和TUBB3在骨转移性乳腺癌中调控微管动力学、内体运输与转移表型的功能研究

《Cancers》：A Novel Function of Tubulin β2a and Tubulin β3 in Regulating Metastatic Properties in Bone-Metastatic Breast Cancer Cells

　　

背景/目的：骨转移（bone metastasis, BM）发生于约70%的晚期乳腺癌（breast cancer, BCa）患者，治疗选择有限，中位生存期仅19–25个月。尽管紫杉醇等微管（microtubule, MT）靶向药物仍是标准治疗，但其疗效常受化疗耐药限制，部分与β-tubulin同种型（β-tubulin isotype）表达改变有关。其中，tubulin β2a和tubulin β3（分别由TUBB2A和TUBB3编码）与不良临床结局相关；然而，它们在骨转移性乳腺癌（bone-metastatic breast cancer, BMBC）中的作用仍不清楚。本研究中，研究人员在骨来源BCa细胞中研究了TUBB2A和TUBB3调控MT动力学和转移行为的功能意义。方法：研究人员评估了原发乳腺肿瘤和骨转移灶中TUBB2A和TUBB3的表达水平。在骨来源MDA-MB-231细胞中实施CRISPR介导的TUBB2A敲除和TUBB3敲低，以评估微管动力学、运输、转移表型和基因表达。结果：TUBB2A和TUBB3表达升高均与患者生存率降低相关。免疫组织化学（immunohistochemistry, IHC）分析显示，与配对原发乳腺肿瘤相比，BM病灶中tubulin β2a丰度增加，同时tubulin β3中等至强阳性细胞比例在BM中也增加。利用CRISPR/Cas9介导的敲除/敲低骨来源MDA-MB-231细胞模型，研究人员发现TUBB2A缺失或TUBB3部分敲低会损害MT再生长、改变细胞周期并破坏内体运输，其特征为核周Rab5阳性囊泡增大。功能上，任一isoform（同种型）的耗竭均降低细胞迁移和侵袭，同时增强黏附，并伴随整合素（integrin）β4、β5、α4、α5和talin-1表达增加。探索性RNA测序（RNA sequencing, RNA-seq）分析提示细胞因子—细胞因子受体信号通路富集，与MT重塑和炎症/分泌程序之间的潜在联系一致。结论：综上，这些发现表明TUBB2A和TUBB3是转移表型的调控因子，也是BMBC的潜在生物标志物和治疗靶点。

研究背景方面，骨转移（bone metastasis, BM）在晚期乳腺癌（breast cancer, BCa）中发生率约65–75%，是导致癌症相关发病和死亡的重要因素，且骨转移性乳腺癌（bone-metastatic breast cancer, BMBC）仍大多不可治愈。微管（microtubule, MT）由α-和β-tubulin异二聚体构成，参与有丝分裂、胞内运输、迁移和黏附等过程；MT动力学失调与侵袭、转移能力和化疗耐药相关。“tubulin code”提示不同β-tubulin同种型（β-tubulin isotype）可通过与微管相关蛋白和马达蛋白差异互作赋予MT网络功能异质性。TUBB2A和TUBB3分别编码tubulin β2a和tubulin β3，TUBB3与紫杉类/长春碱类耐药及不良预后相关，TUBB2A与TGF-β信号和上皮—间质转化相关；但二者在BMBC中的机制作用尚不明确。骨转移微环境中MT稳定性对基质修饰酶、细胞因子和囊泡货物运输尤为重要，α-tubulin乙酰化可增强MT稳定性并促进RANKL等促骨溶解因子分泌；β-tubulin同种型组成改变可能帮助肿瘤细胞适应骨微环境，但特定同种型的作用尚未系统研究。基于此，研究人员分析TUBB2A和TUBB3的临床相关性及其在骨来源BCa细胞中的功能，发现二者高表达与患者生存降低相关，并在骨转移灶中富集；功能缺失会改变MT动力学、内体组织、细胞因子信号、黏附、迁移和侵袭，提示TUBB2A和TUBB3是BMBC转移表型的调控因子及潜在生物标志物和治疗靶点。
关键技术方法方面，研究人员利用Kaplan–Meier Plotter数据库分析乳腺癌患者复发无生存（n=4929），并使用TIMER 2.0数据集比较乳腺浸润癌与邻近正常组织；免疫组织化学分析Rush Pathology来源的配对原发乳腺肿瘤与骨转移活检组织（n=16例患者）；在骨来源MDA-MB-231细胞中采用CRISPR/Cas9进行TUBB2A敲除（KO）和TUBB3敲低（KD），以二代测序（next-generation sequencing, NGS）和Western blot验证；开展微管再生长实验、流式细胞周期分析、Transwell迁移与Matrigel侵袭实验、细胞黏附实验、ELISA检测IL-6分泌，以及每组n=1的探索性RNA-seq（RNA sequencing）转录组分析。
研究结果如下。3.1 乳腺癌中TUBB2A和TUBB3表达升高与不良临床结局相关：通过Kaplan–Meier Plotter分析显示，TUBB2A（HR=1.39）和TUBB3（HR=1.57）高表达与复发无生存降低显著相关（logrank p<0.0001），TUBB1高表达与预后改善相关（HR=0.66），TUBB2B关联较弱（HR=0.90）；IHC显示正常乳腺组织中tubulin β2染色强于β3，而TIMER 2.0数据显示原发乳腺肿瘤中TUBB2A表达降低、TUBB3表达显著升高，提示二者为临床相关且动态调控的β-tubulin同种型。3.2 乳腺癌骨转移肿瘤中tubulin β2和tubulin β3丰度增加：配对患者样本IHC分析（n=16）显示，骨转移灶中tubulin β2强阳性细胞较原发肿瘤增加约1.7倍（9.4% vs. 16.2%）；tubulin β3中等阳性细胞增加1.60倍（5.0% vs. 7.9%），强阳性细胞增加2.22倍（2.6% vs. 5.9%），弱阳性比例和总阳性细胞负荷无实质变化，提示骨转移灶中β2和β3蛋白负荷更高，可能参与肿瘤细胞对骨微环境的适应。3.3 TUBB2A敲除和TUBB3敲低改变有丝分裂进程并延迟微管再生长：NGS证实TUBB2A KO中96%读段携带移码插入缺失，TUBB3 KD中约33%读段携带插入缺失；Western blot显示tubulin β2蛋白降低约75%，tubulin β3降低约25%；qPCR显示TUBB2A KO中TUBB4A表达升高1.8倍，其他同种型无显著变化，γ-tubulin在TUBB3 KD中轻度升高1.36倍。流式细胞术显示TUBB3 KD使G₂/M比例由14.7%降至2.3%，G₁由49.4%升至57.7%，sub-G₁由3.0%升至12.5%，S期无变化；TUBB2A KO无G₂/M降低，仅有较轻的细胞周期重分布。微管再生长实验显示，5 min时各组aster面积相当，30 min时TUBB2A KO和TUBB3 KD再生长严重受损（分别约为scramble的5%和4%）；TUBB3 KD细胞长宽比增加（2.24±0.11 vs. 1.74±0.06），提示细胞更细长。结论为二者缺失改变有丝分裂进程并破坏MT再生长。3.4 TUBB2A KO和TUBB3 KD破坏内体亚细胞分布：通过免疫荧光检测Rab5和caveolin-1（Cav1），scramble细胞Rab5呈核周内体云状点状分布；TUBB2A KO出现核周聚集增强，TUBB3 KD出现向细胞周边延伸的点状梯度。定量显示Rab5强度在TUBB2A KO和TUBB3 KD中分别升高1.6倍和2.0倍，大直径Rab5阳性puncta分别增加6倍和4倍，小puncta数量相对不变；Cav1荧光强度在两种模型中轻度升高，提示MT依赖的内体分布受损、Rab5阳性囊泡增大和空间重分布。3.5 TUBB2A敲除或TUBB3敲低诱导转录重编程：探索性RNA-seq鉴定TUBB2A KO中121个差异表达基因（50个上调、71个下调），TUBB3 KD中122个（45个上调、77个下调），阈值|log2FC|≥1且校正p≤0.05；qRT-PCR验证IL1B、ATP9A、HAS2、EPHB3、GABRA3变化方向与幅度。Reactome和基因本体（Gene Ontology, GO）分析提示细胞因子—细胞因子受体互作通路及IL-10相关炎症信号富集；TUBB2A KO呈现NF-κB/IRF驱动炎症模块，包括IL-1β、IL-6及CXCR2配体趋化因子升高；EMC10为两种模型中改变最显著转录本，MMP1在top100差异基因中表达最高且二者均升高，CSF3、MUC2、CCL2也共同上调。TUBB2A KO下调骨定植/微环境整合相关程序，包括CDH11、Notch/Wnt组分NRARP、NKD1、ANKRD6、TCF7L1，以及OSCAR、CLEC11A、MSR1、ENG；MAGE家族呈同种型特异性，MAGEA1和MAGEA12在TUBB2A KO中下调，MAGEA6在TUBB3 KD中上调。ELISA显示IL-6分泌在TUBB2A KO中升高1.7倍，在TUBB3 KD中降低至0.4倍，提示二者对细胞因子分泌具有不同且相反的作用。3.6 TUBB2A敲除和TUBB3敲低降低细胞迁移和侵袭能力：Transwell迁移实验显示TUBB2A KO迁移降低约20%，TUBB3 KD降低约30%；Matrigel包被Transwell侵袭实验显示TUBB2A KO侵袭降低约20%，TUBB3 KD降低约60%（n=2），提示二者是迁移和侵袭行为的重要调控因子。3.7 TUBB2A缺失和TUBB3部分耗竭改变细胞黏附及黏着斑信号：定量黏附实验显示，在塑料表面TUBB2A KO黏附增加约15%，TUBB3 KD降低约13%；在fibronectin上两种模型黏附均显著增加，在collagen I上呈增加趋势（TUBB2A KO 1.24倍，TUBB3 KD 1.12倍）但未达统计学显著。Western blot显示黏着斑相关蛋白integrin β5、β4、α4、α5和talin-1水平升高，提示TUBB2A和TUBB3通过调节黏着斑组成和细胞骨架动力学影响细胞—基质黏附。
讨论部分总结：研究人员认为β-tubulin同种型并非被动结构成分，TUBB2A和TUBB3高表达与复发无生存降低及骨转移灶富集相关，提示其可能参与骨转移而非仅反映MT结构。不同β-tubulin同种型功能不等价：TUBB2A、TUBB3、TUBB4B高表达与不良结局相关，TUBB2B关联较弱，TUBB1与预后改善相关。既往研究提示TUBB3促进迁移和转移潜能，尤其在脑转移乳腺癌中；本研究扩展至骨转移，显示tubulin β2阳性细胞及中等至强tubulin β3阳性细胞在骨转移中富集，可能促进肿瘤细胞适应骨微环境或与微环境互作。研究队列富集HR⁺/HER2^?肿瘤，许多患者接受过紫杉类和骨靶向治疗，可能影响β-tubulin表达并混杂解释；需更大、更临床多样队列验证。骨微环境具有坚硬矿化基质、高机械应力和与基质/免疫细胞广泛互作；骨灶中β2和β3相对原发肿瘤富集，正常乳腺中β2高于β3而原发肿瘤和骨转移中关系逆转，提示tubulin组成重编程可能改变MT动力学以利于转移行为。TUBB2A或TUBB3耗竭损害MT再生长并改变有丝分裂各阶段比例；TUBB3 KD显著降低G₂/M细胞，支持TUBB3参与细胞周期控制。除细胞周期外，研究揭示β-tubulin同种型在胞内运输和细胞因子信号中的作用：TUBB2A缺失和TUBB3部分耗竭破坏内体组织，导致Rab5阳性囊泡增大和Cav1阳性囊泡增加；二者对IL-6分泌作用相反，TUBB2A耗竭增加IL-6输出，TUBB3敲低降低IL-6，鉴于IL-6在破骨细胞活化和骨转移“恶性循环”中的作用，提示β-tubulin同种型通过MT依赖运输微调炎症信号。研究还显示二者影响细胞—基质互作和侵袭行为：任一isoform耗竭均增加fibronectin黏附、降低迁移能力并呈侵袭降低趋势，伴随黏着斑相关蛋白改变，提示黏附动态失衡稳定细胞—基质互作并限制运动和侵袭潜能。探索性RNA-seq提示β-tubulin同种型为转移转录程序候选调控因子：TUBB2A缺失伴NF-κB/IRF/AP-1炎症网络及促破骨细胞因子/趋化因子升高，同时下调CDH11、Notch/Wnt组分及OSCAR、CLEC11A、MSR1、ENG等骨定植/微环境适应程序，提示TUBB2A缺失可能倾向于炎症性微环境失稳而非稳定骨定植；TUBB3敲低伴MAGEA6等MAGE家族升高，提示TUBB3可能 normally restrain 支持微环境生存和转移可塑性的程序。局限包括缺乏复杂骨微环境、临床队列较小、TUBB3 CRISPR/Cas9敲低效率随传代变异，以及患者特异因素可能影响tubulin同种型表达；需患者来源模型和临床前小鼠模型验证。结论部分翻译：研究结果表明，TUBB2A和TUBB3是骨转移性乳腺癌转移行为的关键调控因子。这些β-tubulin同种型表达升高与患者不良结局相关，并在骨转移灶中丰度增加。功能上降低TUBB2A或TUBB3会破坏微管动力学、细胞周期进程、内体运输、迁移和侵袭，同时增强细胞黏附。此外，转录组分析将这些变化与细胞因子和炎症信号通路联系起来，突出这些同种型在肿瘤—微环境互作中的作用。综上，这些结果确定TUBB2A和TUBB3为转移表型的重要贡献者，并支持其作为乳腺癌骨转移生物标志物和治疗靶点的潜在价值。