《Biomolecules》：Detection of N4-Acetylcytosine in Escherichia coli Using APOBEC3A-Mediated Deamination Sequencing
编辑推荐：
N?-乙酰胞嘧啶（4acC）是一种天然存在的DNA修饰，然而其在原核生物中的存在及分布特征此前尚不明确。本研究基于APOBEC3A（A3A）对胞嘧啶（C）与4acC差异性的脱氨活性，开发了一种单碱基分辨率的测序方法，命名为4acC-A3A-seq（4acC-A
N?-乙酰胞嘧啶（4acC）是一种天然存在的DNA修饰，然而其在原核生物中的存在及分布特征此前尚不明确。本研究基于APOBEC3A（A3A）对胞嘧啶（C）与4acC差异性的脱氨活性，开发了一种单碱基分辨率的测序方法，命名为4acC-A3A-seq（4acC-APOBEC3A测序）。研究人员在利用模板链进行验证并优化反应条件后，将该方法应用于大肠杆菌（E. coli）基因组，成功鉴定出2271个高置信度候选4acC位点，构建了原核生物中首个单碱基分辨率的全基因组4acC图谱。该研究为在单碱基分辨率水平上进一步探索4acC的生物学功能及进化保守性提供了关键技术途径和基础数据。
论文解读文章
一、研究背景与问题提出
DNA化学修饰作为表观遗传调控的核心机制，在基因表达调控、基因组稳定性维持及细胞命运决定等关键生物学过程中发挥重要作用。随着高通量测序技术和质谱分析方法的突破性进展，研究人员已鉴定出多种DNA修饰，显著扩展了遗传密码的复杂性。其中，5-甲基胞嘧啶（5mC）及其衍生物如5-羟甲基胞嘧啶（5hmC）和5-甲酰胞嘧啶（5fC）已被最系统地研究，并在胚胎发育、肿瘤发生及神经退行性疾病等生理和病理过程中展现出关键调控功能。
值得注意的是，近期发现的N位修饰的胞嘧啶，如N?-甲基胞嘧啶（4mC）和N?-乙酰胞嘧啶（4acC），正在重塑研究人员对修饰多样性的认知。与已被充分研究的C5位修饰不同，这些N4位衍生物代表了一个独特且相对未探索的表观遗传学前沿领域。尽管4acC具有重要的科学潜力，但当前对其认知仅局限于单一标志性出版物，其分布和功能的全貌尚未被揭示。
目前，4acC研究面临若干关键挑战。首先，由于其在基因组中丰度极低且固有的化学不稳定性（易发生脱乙酰化），常规检测方法如质谱法和抗体富集测序难以实现高灵敏度和高特异性。其次，缺乏能够在全基因组范围内精确绘制该修饰的单碱基分辨率技术，严重限制了对其他生物学功能的深入分析。此外，4acC在各种生理和病理条件下的动态变化及其调控机制仍不清楚。这些技术瓶颈直接阻碍了该领域的研究进展。
近年来，基于APOBEC（载脂蛋白B mRNA编辑酶催化多肽样）家族脱氨酶的测序技术为DNA修饰研究带来了新机遇。特别是APOBEC3A（A3A）蛋白因其独特的底物特异性而受到广泛关注。研究表明，A3A能够特异性识别并脱氨未修饰的C和5mC，而对4mC保持惰性。这一特性为开发新型DNA修饰检测方法提供了关键思路。从结构上看，4acC和4mC均为N4修饰的C衍生物，4acC带有乙酰基而4mC带有甲基基团。鉴于A3A对4mC表现出脱氨抗性，研究人员推测4acC可能同样对A3A介导的脱氨具有惰性。同时，碱性处理已被确立为特异性去除乙酰修饰的化学方法。基于这些发现，研究人员开发了4acC-A3A-seq方法，以实现大肠杆菌基因组DNA（gDNA）中4acC的单碱基分辨率定位。
二、研究内容与主要结论
研究人员通过系统的实验设计和验证，成功开发并应用了4acC-A3A-seq方法。该方法的核心原理是：经A3A处理后，DNA链中未修饰的C被脱氨为尿嘧啶（U），而4acC由于其与4mC的结构相似性，预期保持不变。经PCR扩增后，原始“C”位点被读取为“T”，而4acC位点被读取为“C”。为消除实验背景并确保高特异性，研究人员引入了碱性预处理作为阴性对照。碱性处理可将4acC脱乙酰化为未修饰的C，使其随后易受A3A介导的脱氨作用。通过比较碱性处理前后的测序结果，可生成全基因组4acC位点图谱。
研究首先通过高效液相色谱-串联质谱法（HPLC-MS/MS）验证了大肠杆菌gDNA中4acC的存在。在约7.5–8.0分钟的保留时间处鉴定出与4acC标准品对应的峰，特征质量跃迁m/z 270.1 → 154.0和270.1 → 112.0明确证实了大肠杆菌基因组中4acC的存在。定量分析显示，大肠杆菌gDNA样本中4acC/dC比率约为0.01%，估计每个大肠杆菌gDNA中约有120个4acC分子。
在方法验证方面，研究人员比较了三种DNA聚合酶的延伸性能，发现Deep Vent（exo-）DNA聚合酶表现最佳，能够成功合成含4acC的链且未检测到脱乙酰化。通过Sanger测序验证，A3A处理后所有未修饰胞嘧啶均完全脱氨并被读取为“T”，而4acC位点保持不变并被读取为“C”。在TC、CC、GC和AC四种序列背景下，A3A均未对4acC表现出明显的脱氨活性。在碱性试剂筛选中，NaHCO?表现出最高的脱乙酰化效率。经NaHCO?预处理后再经A3A处理，原始4acC位点出现清晰的“C→T”转换，验证了NaHCO?介导的化学脱乙酰化的可行性。
在应用于大肠杆菌gDNA时，研究人员发现A3A对gDNA中未修饰C的脱氨效率显著低于合成DNA片段。为解决这一问题，研究人员将A3A孵育时间延长至20小时，并引入已知序列的dsDNA-4acC和dsDNA-C作为spike-in对照，以客观评估每次实验中A3A对未修饰C和4acC的实际脱氨效率。结果显示，dsDNA-C中未修饰C的C-to-T转换率在单独A3A处理后为99.8%，在NaHCO?和A3A顺序处理后为97%。dsDNA-4acC中4acC的转换率在单独A3A处理后约为23%，在NaHCO?和A3A顺序处理后增加至约94%。
通过分析两个技术重复组，研究人员在大肠杆菌gDNA中鉴定出2271个高置信度候选4acC位点。分布分析显示，4acC位点主要位于基因区域内，在基因起始和终止位点密度较高，特别是在基因下游100 bp区域，呈现马鞍形分布模式。乙酰化化学计量学分析显示，鉴定位点的乙酰化率约为20%，考虑约23%的C-to-T转换率后，基因组中4acC修饰的实际乙酰化率估计约为25%。序列偏好性分析显示，大肠杆菌中4acC修饰主要富集于CHH序列背景中，这一分布模式与主要存在于CG背景中的DNA 5mC甲基化显著不同。基序富集分析鉴定出两个高置信度序列基序：WGACCMA和NMGCGNN。其中，WGACCMA基序位点的乙酰化水平显著高于非WGACCMA序列，表明4acC修饰对WGACCMA基序存在偏好性。
三、关键技术方法
研究人员为开展本研究主要运用了以下关键技术方法：高效液相色谱-串联质谱法（HPLC-MS/MS）用于验证和定量检测大肠杆菌gDNA中4acC的存在；基于APOBEC3A（A3A）脱氨酶的酶促脱氨反应，利用其对未修饰C和4acC的差异性脱氨活性实现修饰位点的区分；NaHCO?碱性脱乙酰化处理作为阴性对照，将4acC转化为未修饰C以验证检测特异性；spike-in对照策略，将已知序列的含4acC或C的DNA整合至片段化gDNA中，用于客观评估脱氨效率；高通量测序文库构建及生物信息学分析，包括fastp、bowtie2、bismark等工具进行数据处理和甲基化提取分析。
四、结论与讨论
本研究不仅证实了原核生物基因组中4acC的存在及其低丰度特征，更重要的是建立并验证了一种能够在全基因组范围内以单碱基分辨率检测该修饰的方法学框架。这些发现为后续深入研究原核生物基因组中4acC的分布模式、进化意义和生物学功能奠定了基础。研究中鉴定的2271个位点因此被报告为高置信度候选4acC位点而非最终确认位点。
尽管该方法已实现大肠杆菌全基因组4acC的单碱基测序，但仍存在一定局限性：（1）在A3A处理过程中，部分4acC会发生C到T的转换（通过spike-in数据分析，该部分转换率约为25%），导致测得的4acC乙酰化修饰水平偏低。对于丰度相对较低的4acC修饰，这可能导致平行组间信号重叠减少，因此检测到的位点数量可能少于实际位点。（2）该方法需要较高的位点测序深度以确保不同处理方式下4acC位点C到T信号差异的产生，这使得该方法难以应用于基因组更大的高等生物（如人类、小鼠等）。因此，未来的改进方法可侧重于开发更好的脱氨方法或通过富集信号来降低测序深度要求。
