《Oncogene》：Phylogenetic and phylodynamic approaches to understanding the evolution of non-small cell lung cancer using single-cell sequencing
编辑推荐：
刻画肿瘤细胞的谱系关系及其种群动态对于理解肿瘤演化至关重要。由于直接监测肿瘤演化过程不可行，研究人员利用系统发育学（phylogenetics）和系统动力学（phylodynamics）的进展，基于来自11例非小细胞肺癌（non-small cell lung
刻画肿瘤细胞的谱系关系及其种群动态对于理解肿瘤演化至关重要。由于直接监测肿瘤演化过程不可行，研究人员利用系统发育学（phylogenetics）和系统动力学（phylodynamics）的进展，基于来自11例非小细胞肺癌（non-small cell lung cancer, NSCLC）患者共351个肿瘤细胞的单细胞全外显子组测序数据，探索NSCLC的演化动态，重点关注原发肿瘤及NSCLC早期阶段。研究结果揭示了NSCLC不同演化阶段中系统发育树形状和演化轨迹的显著差异。然而，当前肿瘤模型尚无法完全阐明这一复杂且多阶段的演化过程。为此，研究人员提出假设：肿瘤依次经历以下阶段：（1）通过类似大爆炸（Big Bang）式生长实现初始快速扩增；（2）长期的不进展平台期，伴随可能的中性替换积累；（3）随后适应性替换逐渐积累的阶段；（4）瓶颈事件，即具有选择性优势的亚克隆获得侵袭能力（以间断演化模式为特征），并可能发生导致转移的快速扩增。本工作凸显了单细胞系统发育学和系统动力学在解读NSCLC演化中的应用价值，并有望用于预测肿瘤进展、耐药机制及治疗应答。
**研究背景与意义**
癌症目前位居全球死亡原因首位，非小细胞肺癌（non-small cell lung cancer, NSCLC）是最常见的肺癌组织学类型，约占所有肺癌的85%，总体5年生存率仅28%，约70%的新诊断患者就诊时已处于晚期，预后不良。肿瘤发展最早由Nowell描述为通常跨越数十年的演化过程：肿瘤从单个正常细胞起始，在诸如免疫识别等选择压力下演化，突变在细胞增殖过程中不断积累，导致肿瘤细胞间存在遗传异质性。近年来的TRACERx（TRAcking non-small cell lung Cancer Evolution through therapy (Rx)）等大规模多区域及纵向研究已确立NSCLC演化的多项关键原则，包括早期克隆多样化、广泛瘤内异质性（intratumor heterogeneity, ITH）、亚克隆选择、免疫介导选择和分支转移播散。然而，这些研究大多依赖批量多区域测序，只能从混合细胞群体中推断克隆结构，无法直接重建单个肿瘤细胞的谱系关系。单细胞测序技术为解析ITH和理解肿瘤发生过程中的细胞谱系演化提供了新途径，但该技术存在全基因组扩增（whole-genome amplification, WGA）引入的高扩增错误率、等位基因脱扣和覆盖深度不均等挑战。为此，研究人员利用单细胞全外显子组测序（whole exome sequencing, WES）数据，结合系统发育学（phylogenetics，研究演化单元间谱系关系）和系统动力学（phylodynamics，利用系统发育树推断种群大小变化及影响因素）方法，探索NSCLC的复杂演化动态及转移播散背后的演化机制。
**研究内容与结论概述**
研究人员对来自11例NSCLC患者原发肿瘤的351个单细胞及匹配正常组织bulk样本进行WES，运用双单细胞变异检出器交叉验证获得高置信度单核苷酸变异（single-nucleotide variant, SNV），基于CellPhy构建单细胞系统发育树，采用贝叶斯天际线图（Bayesian skyline plot）推断种群动态历史，结合树形统计与主成分分析（principal component analysis, PCA）聚类，提出NSCLC演化的多阶段假设模型，并讨论了该模型对肿瘤进展预测、耐药机制及治疗应答的潜在意义。
**关键技术方法**
研究样本来源于美国杜克大学健康系统11例接受手术治疗的NSCLC患者（IA–IV期）。单细胞经荧光激活细胞分选（fluorescence-activated cell sorting, FACS）分离后进行全基因组扩增，与匹配正常组织bulk样本一同进行全外显子组测序。SNV检测采用Monovar与SCIΦ双工具交叉验证。系统发育树构建采用CellPhy（最大似然框架，含双倍体基因型错误模型）。种群动态推断采用BEAST2中的phylonco（双倍体基因型GT16错误模型加松弛时钟）。树形评估采用外部/内部支长比（EI-ratio）及三种树平衡指数（Colless、Cophenetic、scaled Sackin），并结合PCA进行分组。
**研究结果**
*临床特征与样本概况*：11例患者中，IA期4例、IB期3例、IIA期1例、IIB期2例（均伴淋巴结转移），另1例为IV期梭形细胞癌。共获得351个单细胞（中位23个/患者，范围7–71个），正常bulk组织平均覆盖深度198.1X，肿瘤单细胞平均覆盖深度101.5X。
*SNV鉴定*：经Monovar和SCIΦ双工具交叉验证，每例患者获得中位4819个SNV（范围1571–7027），其中57.1%–68.6%存在于COSMIC数据库。各患者内部均观察到广泛ITH。肿瘤分期与非同义突变比例呈正相关但未达统计学显著（
