1. 领域背景与文献
发表期刊：未提及；影响因子：未公开；研究领域：甲状腺癌分子生物学、肿瘤微环境与分化调控
甲状腺癌是全球发病率增长最快的内分泌恶性肿瘤，BRAF V600E和RAS突变是最常见的驱动突变，分别占乳头状甲状腺癌（PTC）的40%-80%、滤泡状甲状腺癌（FTC）的30%-50%，定义了具有独特转录谱和临床病理特征的分子亚型。2022年WHO分类细化了分化型甲状腺癌（DTC）的分级体系，区分了高分化与高级别滤泡细胞来源癌，凸显了DTC进展与去分化的异质性。近年单细胞及空间转录组技术的突破，揭示了上皮细胞可塑性与肿瘤微环境（TME）重编程在甲状腺癌进展中的核心作用，但不同驱动突变对这些过程的调控差异尚未明确。当前研究热点聚焦于TME靶向治疗与精准分型策略，但未解决的核心问题包括：BRAF与RAS突变亚型的去分化轨迹差异、TME细胞互作的突变特异性机制，以及晚期去分化甲状腺癌（ATC）尤其是RAS突变亚型缺乏有效靶向治疗方案的困境。针对上述空白，本研究整合多组学技术解析两种突变亚型的去分化轨迹及TME互作模式，为甲状腺癌的精准诊疗提供理论依据与潜在靶点。
2. 文献综述解析
作者对现有研究的分类维度为是否区分驱动突变亚型，将领域内研究分为未分层的甲状腺癌研究与聚焦特定突变亚型的研究两类。现有未分层研究的优势在于揭示了上皮可塑性与TME在甲状腺癌进展中的普遍作用，局限性在于未考虑驱动突变的异质性，导致对去分化机制的理解存在偏差，无法解释临床中治疗响应的显著异质性；少量聚焦特定突变的研究仅关注单一亚型，缺乏系统的跨亚型对比。本研究的创新价值在于首次整合单细胞、空间转录组与 bulk 多组学数据，系统对比BRAF V600E与RAS突变亚型的去分化轨迹及癌相关成纤维细胞（CAF）-上皮互作差异，明确了突变特异性的肿瘤进展机制，填补了领域内对驱动突变调控TME重编程的认知空白，为亚型特异性治疗提供了潜在靶点。
3. 研究思路总结与详细解析
本研究的整体框架为：以BRAF与RAS突变驱动的甲状腺癌去分化的异质性为核心科学问题，研究目标是明确两种突变亚型的去分化轨迹及TME互作机制；技术路线遵循“样本覆盖-多组学测序-生物信息学解析-实验验证-结论提炼”的闭环逻辑，涵盖从基因组到空间层面的多维度分析。
3.1 样本收集与多组学测序平台构建
实验目的是获取覆盖不同进展阶段的BRAF与RAS突变甲状腺癌样本，构建多组学数据集以支撑后续分析。方法细节：收集16例新鲜冷冻肿瘤组织（含BRAF突变PTC（PTC-B）、RAS突变FTC（FTC-R）及对应的ATC（ATC-B、ATC-R）），进行全基因组测序（WGS）、全转录组测序（WTS）、单细胞核RNA测序（snRNA-seq）；20例福尔马林固定石蜡包埋（FFPE）样本进行10x Visium空间转录组分析；96例甲状腺癌组织构建组织芯片用于免疫组化验证；同时整合公共数据库的311例WTS、94例单细胞RNA测序（scRNA-seq）及1例空间转录组样本。结果解读：成功构建了覆盖基因组、转录组、单细胞、空间维度的多组学数据集，完整覆盖甲状腺癌从分化到去分化的进展过程，为解析突变特异性机制提供了全面的数据基础。实验所用关键产品：Qiagen DNeasy Blood & Tissue Kit、Twist NGS Library Kit、Illumina NovaSeq6000平台、10x Genomics Chromium Next GEM Single Cell 3’ Kit v3.1、Sigma-Aldrich的L1CAM抗体（货号L4543）、Abcam的ITGAV抗体（货号ab179475）等。
3.2 基因组与 bulk 转录组特征分析
实验目的是解析BRAF与RAS突变亚型的基因组及转录组差异，明确去分化过程中的分子特征变化。方法细节：通过WGS分析突变谱、拷贝数变异（CNV）及结构变异；WTS分析基因表达差异及通路富集；采用ESTIMATE算法评估TME细胞浸润水平。结果解读：基因组分析显示ATC-B/R均同时携带TERT启动子（TERTp）与TP53突变，而RAS驱动的ATC结构变异更显著；转录组分析显示BRAF突变亚型的去分化伴随TME组成重塑，RAS突变亚型则出现血管生成、上皮-间质转化（EMT）、缺氧通路的 abrupt 激活；ESTIMATE算法结果显示ATC-B的免疫与基质评分升高幅度显著高于ATC-R（n=16，P<0.01）。
3.3 单细胞转录组与细胞亚型鉴定
实验目的是解析肿瘤上皮细胞及TME细胞的亚型组成，明确去分化过程中的细胞状态变化。方法细节：对102845个细胞进行snRNA-seq聚类分析，通过Copykat进行CNV分析，采用Monocle进行伪时间轨迹分析。结果解读：鉴定出5种主要细胞类型（肿瘤上皮、T/NK、髓系、成纤维细胞、内皮细胞）；上皮细胞进一步分为7个亚型（基于甲状腺分化评分（TDS）与ERK激活评分），其中ATC-ERK-high亚型为非整倍体细胞，TDS最低且ERK激活评分最高；BRAF驱动的去分化呈现连续的轨迹变化，而RAS驱动的去分化则表现为 abrupt 的细胞状态转变，ATC-R的非整倍体细胞比例显著高于其他亚型（文献未明确提供该数据，基于图表趋势推测）。
3.4 伪时间轨迹与去分化机制解析
实验目的是明确两种突变亚型的去分化轨迹差异，解析背后的分子通路变化。方法细节：基于TDS对上皮细胞进行伪时间轨迹分析，对比BRAF与RAS突变亚型的通路富集差异。结果解读：BRAF驱动的甲状腺癌去分化为连续的渐进过程，伴随免疫通路的逐步激活；RAS驱动的去分化则为 abrupt 的转变过程，伴随缺氧、ECM重塑通路的显著激活；伪时间分析显示RAS突变亚型在去分化过程中基因组不稳定性快速升高，非整倍体细胞比例从FTC-R的低水平骤升至ATC-R的高水平。
3.5 细胞互作与空间转录组验证
实验目的是在空间层面验证细胞亚型分布及CAF-上皮细胞互作的突变特异性。方法细节：采用CellPhoneDB分析配体-受体互作，通过Visium空间转录组进行细胞亚型的空间定位，利用Stopover工具验证配体-受体对的空间共定位。结果解读：ATC中CAF与上皮细胞的互作强度显著高于DTC，且RAS突变亚型的互作强度提升幅度更大；BRAF突变ATC中以整合素相关互作（如L1CAM-ITGAV）为主，RAS突变ATC中则富集PLAU-PLAUR、TNFSF10-TNFRSF10B等互作对；空间转录组验证显示这些互作对在ATC区域特异性共定位。
3.6 CAF亚型分析与免疫组化验证
实验目的是解析CAF的亚型变化及关键互作分子的蛋白水平表达。方法细节：对成纤维细胞进行亚型聚类，分为炎性CAF（iCAF）与肌成纤维细胞样CAF（myoCAF）；通过免疫组化检测L1CAM、ITGAV、PLAU等分子的表达。结果解读：进展至ATC时，两种突变亚型的iCAF比例均显著升高；免疫组化显示L1CAM、ITGAV、PLAU等分子在ATC中特异性高表达，且呈现突变亚型差异，如PLAU在RAS突变ATC中表达更显著；生存分析显示这些互作对的高表达与较差总生存期相关（log-rank检验P<0.05，文献未明确提供HR值）。
4. Biomarker研究及发现成果解析
Biomarker定位
本研究涉及的Biomarker包括三类：一是上皮细胞亚型（如ATC-ERK-high亚型），筛选逻辑为snRNA-seq聚类分析→空间转录组验证→生存分析；二是CAF亚型（iCAF），筛选逻辑为单细胞聚类→多组学整合验证；三是配体-受体互作对（L1CAM-ITGAV、PLAU-PLAUR等），筛选逻辑为配体-受体互作分析→空间共定位验证→免疫组化及生存分析。
研究过程详述
Biomarker均来源于甲状腺癌组织的多组学及蛋白水平数据：上皮细胞亚型与CAF亚型基于snRNA-seq数据鉴定，配体-受体互作对基于单细胞互作分析筛选；验证方法包括空间转录组的空间定位、免疫组化的蛋白水平验证，以及公共数据集的生存分析；特异性与敏感性方面，ATC-ERK-high亚型在ATC中特异性富集，ROC曲线数据未明确提供，但生存分析显示该亚型比例高的患者预后显著较差；L1CAM-ITGAV互作对在BRAF突变ATC中特异性高表达，PLAU-PLAUR在RAS突变ATC中特异性高表达，Kaplan-Meier曲线显示高表达组总生存期显著短于低表达组（log-rank检验P<0.05）。
核心成果提炼
本研究的核心成果包括：首次明确BRAF与RAS突变亚型的去分化轨迹差异，BRAF驱动的去分化为连续渐进过程，RAS驱动的为 abrupt 转变过程；iCAF是两种亚型去分化的关键TME组分，且突变特异性的CAF-上皮互作通路（BRAF亚型的整合素通路、RAS亚型的PLAU-PLAUR通路）与不良预后显著相关；L1CAM、ITGAV、PLAU等分子可作为ATC的潜在预后Biomarker，其中PLAU-PLAUR互作对为RAS突变ATC提供了潜在治疗靶点。
