簇、指状突起与单细胞：肿瘤侵袭的机械学景观

《PLOS Computational Biology》：Clusters, fingers, and singles: A mechanical landscape of tumor invasion

　　

集体侵袭是肿瘤播散和转移的关键机制，涉及异质性细胞群的协调迁移。基于球体（spheroid）的实验研究已识别出专门的领导者（leader）细胞和跟随者（follower）细胞在过程中协同作用，但支配其相互作用的生物物理规则仍不清楚。研究人员提出一种基于机制的细胞层面计算模型，使用细胞势能（Cellular Potts）框架，研究异型黏附、领导者运动性和跟随者增殖如何共同塑造侵袭。领导-跟随肿瘤在13 310组参数下被模拟，侵袭通过侵袭面积、浸润面积、指状突起、单独叛逃细胞和脱落簇来量化。从这些模拟中，研究人员识别出四种不同侵袭表型：非侵袭、整体集体、单细胞和多模态。多模态侵袭——即黏聚条带、单独细胞和小簇共存——是最普遍表型，尤其在中等黏附和高运动性下出现。增殖增加肿瘤体积而非决定侵袭模式，模式主要由黏附和领导者运动性支配。在参数空间上绘制结果揭示出侵袭模式间的急剧转变，凸显黏附与运动性在塑造侵袭复杂性中的权衡。结果表明，先前被认为罕见的混合侵袭行为可从简单机械规则中稳健涌现，并在广阔参数区域占优。该框架把侵袭的二元模型与实验观察到的异质性调和起来，为通过调节黏附和运动性改变侵袭行为提供预测性见解。

研究背景方面，转移是癌症相关死亡主要原因，局部侵袭常由单细胞独立播散或群体协调集体侵袭驱动。集体侵袭在乳腺癌、结直肠癌和鳞癌中日益被认识，其标志是迁移前沿出现领导者-跟随者功能异质性：领导者细胞具有高方向性运动性、趋化响应、基质金属蛋白酶介导的胞外基质（extracellular matrix，ECM）重塑及增殖受抑；跟随者细胞靠黏附连接维持结构、高增殖并通过力学耦合被动迁移。既有计算模型多把侵袭简化为集体或单独二元类别，少有模型显式表征领导者-跟随者异质性，也少有系统探索黏附、运动性与增殖交互如何产生真实肿瘤中条带、单细胞和脱落簇共存的混合行为。为隔离集体迁移的基本机械交互，研究人员构建理想化体外侵袭模型，而非模拟体内全复杂度，从而提供可解释的机制地图。

该研究发表于《PLOS Computational Biology》。研究人员使用软件CompuCell3D实现格点化细胞势能模型（Cellular Potts Model，CPM），在二维矩形域中初始化由领导者细胞占25%和跟随者细胞占75%组成的基底肿瘤片，施加沿y轴的静态非扩散趋化梯度，用Metropolis算法更新格点使系统能量最小化；能量项包含细胞-细胞接触能（J_LF表示领导者-跟随者接触能，低J为强黏附）、体积约束和领导者趋化项。参数扫描覆盖J_LF、迁移系数（migration coefficient，μ）和跟随者增殖概率（proliferative probability，PP）各11水平，共11³=1331组唯一参数、每组10次重复，总计13 310次模拟，运行700 MCS约36小时。表型由基于图连通分量的主瘤块识别、凸包法算浸润面积、峰检测算法数指状突起、无邻接单细胞定义singles、广度优先搜索（breadth-first search，BFS）找脱落簇，再按规则分为非侵袭、单细胞侵袭、整体侵袭、多模态侵袭四类。样本设定源自非小细胞肺癌（non-small cell lung cancer，NSCLC）球体实验系统如SaGA平台，但不含体内基质与间质微环境。

结果部分保留小标题说明。Distinct invasion metric sensitivities to adhesion, motility, and proliferation：研究人员通过全参数扫描测侵袭面积、浸润面积、指状突起数、singles数和脱落簇数，发现侵袭面积和指状突起在中等黏附与高μ下最大；浸润面积和singles随μ单调上升、弱黏附时singles剧增；脱落簇在中等至弱黏附与高μ下最频。增殖概率PP与所有形态指标几乎零相关，说明增殖只扩肿瘤块，不决定模式。Phenotypic phase-space mapping reveals prevalence of multimodal invasion：在J_LF-μ-PP三维参数空间绘表型图，非侵袭占约22%位于低μ中强黏附区，单细胞侵袭占窄带位于弱黏附低中μ区，整体集体侵袭约23%位于强黏附中μ中PP区，多模态侵袭约占54%位于中弱黏附、高μ、中高PP区，是最常见结果，证明混合侵袭是鲁棒现象。Cluster formation peaks under intermediate adhesion and strong migration：集群分析显示中间黏附（J_LF约0–3）与强μ（20–30）下集群发生率和总细胞数最高，簇多含4–8细胞、均值约7细胞，领导者富集达约60%–70%；强黏附或低μ抑制脱离，证明集群脱离由黏附-运动性平衡控制，且集群对应部分上皮-间质转化（partial epithelial–mesenchymal transition，pEMT）式力学状态。

讨论与结论翻译部分：本研究细化了集体与单独侵袭的二元区分。在1331组参数、10次重复的13 310次模拟中，涌现四种可重复表型：非侵袭、整体集体侵袭、单细胞侵袭和多模态侵袭。在均匀采样参数范围内，多模态侵袭占据最大相空间区域（约54%），而纯单细胞侵袭相对罕见；这不应被解释为体内先验概率。结果表明，领导者运动性与细胞间黏附是侵袭架构主要调节者，增殖主要作为肿瘤体积缩放因子。中间黏附与高运动性构成允许簇脱离、领导者驱动持续性和转移播散的力学许可区。该框架把实验观察到的异质前沿与机械规则联系起来，并为联合调节黏附与运动性以把肿瘤从多模态区域移开提供可检验策略。

总体而言，论文价值在于用CPM把领导者-跟随者异质性、异型黏附、方向性运动性和增殖解耦分析，证明多模态侵袭不是实验偶发现象，而是简单机械规则在广阔参数区的必然产物；这为解读肿瘤异质前沿、预测黏附或运动性干预后果提供了定量地图，对思考抗转移策略中“抑制运动性”“增强黏附以压缩播散”具有概念性意义。