《Nature Communications》：Single-cell mapping of chromosome breaks identifies multiple fragile site classes with distinct DNA replication timing landscapes
编辑推荐：
常见脆性位点（common fragile sites, CFSs）是在复制压力条件下易发生染色体断裂的基因组区域。由于以往鉴定主要依赖细胞群体水平分析，低频非常见脆性位点（non-CFSs）可能被忽视。研究人员利用单细胞全基因组测序（single-cell
常见脆性位点（common fragile sites, CFSs）是在复制压力条件下易发生染色体断裂的基因组区域。由于以往鉴定主要依赖细胞群体水平分析，低频非常见脆性位点（non-CFSs）可能被忽视。研究人员利用单细胞全基因组测序（single-cell whole-genome sequencing, scWGS）对人骨肉瘤U2OS细胞进行aphidicolin诱导的染色体断裂全景分类。结果显示，42%的断裂具有CFS特征，其余58%为在整个S期均可出现的罕见non-CFSs。在non-CFSs中，研究人员发现与转录-复制冲突相关的早S期复制断裂，以及两类新断裂：一类依赖转录、定位于早到晚S期复制时序（replication timing, RT）转换区；另一类呈晚S期RT，且在抑制转录后断裂频率升高。研究人员进一步发现，未处理细胞中不同的RT景观与aphidicolin处理后不同的RT行为和断裂类别相关联。该工作揭示了不同断裂类别的频率、RT景观及机制差异，为理解轻度复制压力下基因组不稳定区域提供了全景视图。
研究背景方面，基因组不稳定性是癌症发生与衰老的重要驱动因素。哺乳动物基因组中存在天然易于不稳定的区域，其中常见脆性位点（common fragile sites, CFSs）研究最充分。CFSs最早由低剂量aphidicolin（真核DNA聚合酶抑制剂）处理后在中期染色体上观察到的裂隙和断裂来定义，并被认为多在mid-late S期复制、与大基因和延迟复制相关。群体水平的MiDAS-seq显示G2/M期DNA合成位点与CFSs重合；另有一类早复制脆性位点（early replicating fragile sites, ERFSs）在小鼠B细胞中被鉴定。但CFS与ERFS实验条件不同，且其相对频率与低频non-CFSs长期未被系统研究。由于群体分析会漏掉罕见位点，且已有研究提示小拷贝数变异（copy number variations, CNVs）可能来自尚未定义的脆性位点类别，因此研究人员开展单细胞水平全基因组断裂制图。
研究人员以人骨肉瘤U2OS细胞为模型，进行低剂量aphidicolin单周期处理，分选单个G1期细胞做scWGS，用AneuFinder在500-kb分辨率下做拷贝数核型，定义染色体断裂区域并剔除固有断裂，得到aphidicolin诱导断裂集；同时用10细胞Repli-seq构建假批量（pseudobulk）复制时序（replication timing, RT）曲线，并结合aphidicolin诱导的MiDAS-seq信号，按RT与CFS/non-CFS两个轴将断裂分为五类。结论上，研究人员证明除经典CFS外，还存在多类non-CFS断裂，其RT景观、转录依赖性和复制压力响应均不同；未处理细胞的RT地形可预测aphidicolin后的RT改变与断裂类别。该研究发表于《Nature Communications》，意义在于把脆性位点研究从“群体共性”推进到“单细胞异质性与类别全景”，为癌症基因组不稳定性机制提供新框架。
关键技术方法上，研究人员使用人骨肉瘤U2OS细胞系作为样本来源；采用scRepli-seq框架下的单细胞全基因组测序（scWGS）做G1期拷贝数断裂识别；用多分数10细胞Repli-seq获得pseudobulk RT；用MiDAS-seq定义CFS断裂；用BrdU-IP Repli-seq与very-late-S/G2/M BrdU-IP校正RT判断偏差；用GRO-seq评估转录活性；用Hi-C做TAD边界与A/B区室分析；用immuno-FISH和γH2AX、FANCD2免疫荧光分析转录抑制下损伤位点；用DRIP-seq、G4-seq、NS-seq等公共数据注释R-loop、G-quadruplex和复制起始。
研究结果如下。
Copy number-based karyotyping in single mammalian cells：研究人员通过单G1细胞拷贝数核型发现，aphidicolin使U2OS细胞间核型异质性增加。
Identification of chromosome breaks induced by low-dose aphidicolin：研究人员定义并筛选出274个aphidicolin诱导断裂区域，其中相当比例与已知CFS数据库重合，证明实验体系可捕获真实CFS。
Identification of five classes of chromosome breaks based on CFS and RT：研究人员用pseudobulk RT与MiDAS信号建立RT-MiDAS轴，将断裂分为CFS-L、CFS-M、non-CFS-L、non-CFS-M、non-CFS-E五类，且未观察到早复制CFS型。
CFS-L and CFS-M breaks show classic CFS features：研究人员通过RT延迟、G2/M合成增强、靠近长基因等，确认CFS-L与CFS-M具经典CFS特征，且不偏好TAD边界。
Non-CFS-L RT-delayed breaks represent a category distinct from CFS-L：研究人员通过GRO-seq、immuno-FISH与DRB转录抑制实验，发现non-CFS-L RT-delayed断裂转录活性低，且抑制转录后γH2AX升高，不同于CFS-L，代表一类晚复制但非CFS的脆性区域。
Steep late-S RT valleys are associated with CFS-L breaks：研究人员对终止区（termination zones, TZs）聚类，发现陡的深late-S RT谷对应慢叉移动和G2/M期低复制，是CFS-L的特征性RT地形。
Non-CFS-M RT-delayed breaks are distinct and coincide with early/late-S TTRs：研究人员发现non-CFS-M RT-delayed断裂位于早/晚S期复制转换区（transition regions, TTRs），抑制转录后损伤下降，说明其依赖转录且位于早期到晚期RT边界。
Early/late-S TTRs exhibit aphidicolin-induced under-replication in G2/M：研究人员对TTRs聚类显示，早/晚S边界TTRs在aphidicolin后于较晚RT侧发生G2/M低复制，伴随A/B区室边界特征。
Non-CFS-E RT-advanced breaks show features of transcription-replication conflicts：研究人员用FANCD2、DRIP-seq、G4-seq、NS-seq等显示，non-CFS-E RT-advanced断裂富集于高转录、R-loop、G-quadruplex和复制起始区，符合转录-复制冲突（transcription-replication conflicts, TRCs）机制。
DNA damage at non-CFS-E RT-advanced breaks decreased upon transcription inhibition：研究人员用γH2AX-FANCD2共定位分析发现，仅早S期在转录抑制后共定位下降，支持non-CFS-E RT-advanced断裂由TRCs驱动。
RT landscapes in untreated cells can distinguish different break classes：研究人员总结未处理细胞中的steep RT valley、TTR、early-S转录活跃区等RT景观可区分不同断裂类。
讨论与结论部分中，研究人员指出群体水平CFS制图会低估脆性位点多样性；pseudobulk 10-cell Repli-seq比传统early-versus-late BrdU-IP更不易被低复制偏差误导。结论为：轻度复制压力下染色体断裂不是单一CFS现象，而是由late-S陡谷、early/late-S TTR、early-S转录活跃区等不同RT景观所塑造的多类脆性位点；其中CFS-L/CFS-M源于under-replication与G2/M MiDAS，non-CFS-L晚复制断裂与转录抑制敏感相关，non-CFS-M RT-delayed依赖转录并位于TTR，non-CFS-E RT-advanced由TRCs、R-loop和G-quadruplex促成。该分类框架揭示了脆性位点的频率、RT景观与机制异质性，为癌症基因组不稳定性及复制压力相关断裂研究提供了单细胞尺度的全景基础。
