整合单细胞转录组与染色质可及性图谱解析肺腺癌（LUAD）细胞放射反应异质性及放疗增敏候选靶点

《MedComm》：Dissecting the Radioresponse Heterogeneity of Lung Adenocarcinoma Cells From the Perspective of Transcriptional and Regulatory Landscape

　　

放射抵抗严重削弱肺腺癌（LUAD）的放疗疗效，但同基因肿瘤细胞间细胞放射反应的内在异质性仍未被充分表征。在此，研究人员整合单细胞RNA测序（scRNA-seq）与单细胞转座酶可及染色质高通量测序（scATAC-seq），解析单细胞来源LUAD细胞接受6 Gy γ射线照射后异质性放射反应背后的转录与调控景观。研究人员发现，同基因LUAD细胞具有决定放射敏感性差异的内在分子因素，并在照射后表现出异质性放射反应性。照射后6 h出现一种与放射相关的转录细胞状态，并伴随DNA损伤修复、凋亡及细胞周期调控。抑制赖氨酸去甲基化酶4B（KDM4B）可通过加剧DNA损伤增强LUAD放射敏感性，而细胞周期蛋白B1（CCNB1）和DNA拓扑异构酶IIα（TOP2A）表达降低可能促进放射触发的G₂/M期阻滞。值得注意的是，放射相关调控细胞状态出现早于放射诱导的转录细胞状态。RNA聚合酶II亚基A（POLR2A）在高反应性LUAD细胞中发挥核心调控作用。基于这些发现，研究人员提出LUAD细胞放疗增敏的假说。总之，本研究为LUAD放射反应异质性提供多层机制见解，并为开发改进的放疗增敏策略、推进LUAD精准放疗提供可操作的分子候选物。

放射抵抗是影响肺癌放疗效果的主要障碍，尤其对最常见的组织学亚型肺腺癌（LUAD）而言，其可导致治疗失败、复发和转移。现有研究提示，放射敏感性在患者、肿瘤乃至同基因癌细胞间存在显著差异，但具体分子机制尚未充分阐明；传统免疫荧光和流式细胞术虽可评估DNA损伤、细胞周期阻滞和凋亡，却难以同时刻画癌细胞放射敏感性与多基因调控效应。单细胞测序已在肿瘤异质性等研究中广泛应用，但尚未充分用于解析癌细胞放射反应异质性。基于此，研究人员对单细胞来源LUAD细胞行6 Gy γ射线照射，整合scRNA-seq与scATAC-seq，从转录和调控景观解析其放射反应异质性，提出放疗增敏假说并筛选候选靶点。该研究发表于《MedComm》，为LUAD放射反应异质性提供了多层机制见解和可操作分子候选，对推进精准放疗具有意义。

主要技术方法（约250字）：研究采用单细胞源性A549肺腺癌细胞（TP53野生型，购自GeneChem）单克隆培养，以6 Gy ⁶⁰Co γ射线（61.42 cGy/min）照射；分别收集未照射、照射后2 h和6 h样本行scRNA-seq，收集未照射、照射后即刻和2 h样本行scATAC-seq。生物信息学使用scCANCER、MELD、Signac、cisTopic、SCENIC、CellChat、PAGA-tree、dynamo RNA velocity等。实验采用CCK-8、γH2AX免疫荧光、7-AAD/Annexin V凋亡检测和PI细胞周期流式分析，并以NCGC00244536抑制KDM4B。临床分析使用TCGA中107例接受放疗的LUAD患者数据；TP53缺失/突变分析使用GEO数据集。

研究结果如下。2.1 同基因LUAD细胞具有异质性放射敏感性的分子基础：对单细胞源性A549细胞行scRNA-seq（13 619个细胞）和scATAC-seq（8 927个细胞），scRNA-seq识别10个簇，功能富集主要涉及细胞周期和DNA构象变化；簇10 ATP合成活跃，提示放射抵抗特征；scATAC-seq亦识别10个染色质可及性不同的簇，包括RAD51B等异质性。结论是即便同基因LUAD细胞也存在与细胞周期和基因组构象相关的转录与调控异质性。

2.2.1 利用scRNA-seq数据量化受照LUAD细胞的放射反应：照射后2 h和6 h行scRNA-seq，以MELD计算电离辐射相关似然并分为高、中、低反应；2 h以中反应为主，6 h高反应细胞最多；高反应细胞2 h富集DNA损伤过程，6 h富集凋亡和细胞周期调控；实验显示6 h增殖下降、γH2AX信号2 h更明显且分布不均、6 h凋亡和G₂/M阻滞增加。结论是MELD分类与实验结果一致，放射反应具时间性和异质性。

2.2.2 照射后6 h出现放射相关转录细胞状态：整合照射与对照scRNA-seq，2 h未见新转录状态，6 h出现放射相关转录细胞状态且受放射影响最强；该状态高表达基因富集TP53介导DNA修复、凋亡和细胞周期调控，主要为G₁/S期且DNA损伤不可修复；低反应簇高表达基因富集ATR通路DNA损伤修复。结论是6 h细胞命运取决于DNA损伤与修复平衡。

2.2.3 放射相关转录细胞状态的关键基因：重建TP53信号通路；TCGA中107例放疗LUAD按放射敏感性指数（RSI，中位0.225）分组，KDM4B和CDKN1A在放射抵抗病例中上调；NCGC00244536抑制KDM4B后6 h放射敏感性和DNA损伤增加；GEO数据示TP53敲除/突变下PHLDA3、BAX、PURPL下调，而CCND1、ELF3、PFKFB3不依赖TP53。结论是KDM4B为候选增敏靶点，TP53状态显著调节放疗反应，筛选靶点需考虑TP53背景。

2.3.1 G₂/M期阻滞LUAD细胞中细胞周期转换关键基因降低：通过PAGA-tree轨迹、dynamo RNA velocity和细胞周期注释分析，6 h出现明显G₂/M阻滞，CCNB1和TOP2A在阻滞细胞中下调，RNA velocity示6 h G₂/M至G₁转录转变改变。结论是放射诱导转录变化与G₂/M阻滞密切相关，CCNB1/TOP2A降低与阻滞相关，但原文指出其确切关系仍需进一步阐明。

2.3.2 细胞间通讯调控受照LUAD细胞G₂/M阻滞：CellChat分析显示6 h G₂/M细胞中TGFβ、GDF和SPP1介导入信号；TGFβ2和GDF15相关配受体通讯概率最高，TGFβ2信号主要来自G₂/M细胞，SPP1主要来自G₁细胞并由ITGB1-ITGA4/5调控；原文提示TGFβ2与SPP1信号可能对G₂/M阻滞起相反作用。结论是细胞间通讯参与放射诱导G₂/M阻滞调控。

2.4.1 受照LUAD细胞的特征性调控变化：照射后即刻和2 h行scATAC-seq，Signac/chromVAR示TEAD、FOSL等转录因子染色质可及性变化；照射后失去可及性的motif富集5′-C*C或5′-CC*二核苷酸序列，提示与放射诱导DNA损伤位点相关。结论是放射早期可检测调控层染色质变化，提示相关转录因子在基因组脆弱性中的潜在作用。

2.4.2 放射相关调控细胞状态富集放射相关转录因子：整合3个scATAC-seq数据集，cisTopic识别12个topic，topic 1特异贡献于照射细胞；topic 1区域主要位于远端基因间和内含子等调控区，富集POLR2A、E2F1、YY1、POU4F1、DEAF1、GATA1、NKX2-8等motif。结论是存在放射相关调控细胞状态并富集放射相关转录因子。

2.4.3 POLR2A是高反应LUAD细胞中调控放射反应的关键转录因子：SCENIC从2 h scRNA-seq识别331个regulon，POLR2A regulon在高反应细胞中活跃且在放射相关调控状态富集；POLR2A调控多个基因，CYB561D1推断调控重要性最高，参与铁死亡调控。结论是POLR2A为高反应LUAD细胞的关键调控因子和候选核心靶点。

2.5 LUAD细胞放疗增敏假说：研究人员提出三阶段假说，即初始阶段由细胞周期状态和染色质结构等内在特性决定放射敏感性；反应阶段以DNA损伤与修复平衡、细胞周期检查点、CCNB1/TOP2A降低及TGFβ2通讯参与G₂/M阻滞为核心；命运决定阶段中6 h放射相关转录状态内少量G₁/S期细胞DNA损伤严重并发生凋亡，且调控状态早于转录状态，POLR2A为核心转录因子。

讨论部分总结认为，LUAD细胞存在与细胞周期动态和3D基因组组织密切相关的内在转录与调控异质性，可能促成同基因癌细胞间放射反应差异；放射相关调控状态早于转录状态出现；TP53信号激活是放射相关转录状态特征，但并非所有关键基因均依赖TP53；KDM4B为候选增敏靶点，CCNB1/TOP2A降低及TGFβ2通讯与G₂/M阻滞相关，POLR2A为早期高反应细胞核心调控因子。三阶段模型强调初始阶段干预机会。局限包括仅使用TP53野生型A549细胞、二维单层培养缺乏肿瘤微环境，需遗传工程细胞、患者来源类器官、单细胞/空间转录组及单细胞DSB定位技术进一步验证。研究结论翻译：总之，本研究通过整合转录与调控景观，为同基因LUAD细胞间放射反应异质性提供新见解。基于这些发现，研究人员提出一种放疗增敏策略并鉴定若干候选分子靶点，为推进肺腺癌精准放疗的机制与转化研究奠定基础。