《Antioxidants》：The Yeast Metabolic Cycle as a Tractable Cellular Framework for Redox Timing, Redox Buffering, and Transcriptome Fidelity
编辑推荐：
酿酒酵母中的酵母代谢周期（YMC）提供了一个易处理模型，用于研究线粒体呼吸、氧化还原计时和代谢阶段如何塑造转录组丰度和保真度。研究人员在低溶解氧（Low-DO）/高呼吸和高溶解氧（High-DO）/低呼吸阶段进行了核糖体RNA去除的全转录组RNA测序（WRS）
酿酒酵母中的酵母代谢周期（YMC）提供了一个易处理模型，用于研究线粒体呼吸、氧化还原计时和代谢阶段如何塑造转录组丰度和保真度。研究人员在低溶解氧（Low-DO）/高呼吸和高溶解氧（High-DO）/低呼吸阶段进行了核糖体RNA去除的全转录组RNA测序（WRS）和RNA-seq衍生的错配分析。在1505个阶段差异表达基因中，Low-DO富集了核糖体生物合成、rRNA加工、翻译、硫代谢和蛋白质合成，而High-DO富集了氧化剂解毒、氧化还原酶活性和氧化还原缓冲相关通路。广义线性混合模型识别出一个替换类依赖的Low-DO相关RNA-seq错配响应。最强的错配水平增加发生在合并的C>T/G>A兼容类中，而C>A/G>T并未增加。这种模式与简单的单一损伤模型不一致，而是支持一个混合的Low-DO相关RNA-seq序列不一致景观。变异率建模还检测到T>C/A>G和T>A/A>T的增加，表明多个生物学和技术过程可能促成了观察到的谱。复发分析显示，大多数被调用的变异是样本特异性的，支持一个瞬时的RNA-seq错配景观，而非稳定的DNA突变。这些发现将YMC建立为一个简化的真核框架，用于研究代谢阶段和氧化还原状态如何塑造转录组保真度。
**研究背景与问题**
生物信息流（DNA→RNA→蛋白质）并非线性生化序列，而是受时空组织调控的过程。在真核细胞中，转录、RNA加工、核输出、翻译、降解和质量控制发生在动态的核质结构中。同时，细胞代谢并非恒定不变：能量生产、氧化还原状态、氧化还原缓冲需求和生物合成活性随生物周期波动，导致大分子在不同时间窗口受到差异化保护或易受损伤。酵母代谢周期（YMC）在葡萄糖限制恒化器培养中通过同步的溶解氧（DO）、呼吸、ATP产生、氧化还原状态和基因表达振荡，为研究生物计时、线粒体活性、氧化还原缓冲和转录组调控的协调提供了实验模型。关键特征是生物合成与基因组保护功能的时序分离：低DO阶段（高呼吸/高耗氧）对应核糖体生物合成和较氧化环境，高DO阶段（低呼吸）富集解毒和氧化还原缓冲程序。然而，代谢阶段是否以及如何影响转录组保真度（即RNA序列一致性）尚不清楚。本研究旨在利用YMC的简化系统，检验代谢阶段（尤其是低DO高呼吸期）是否与RNA-seq衍生的序列不一致（错配）景观相关，并探索其替换类依赖性。
**研究内容与结论**
研究人员使用非工程酵母菌株DBY12007（来源于Hickman和Winston，后经Botstein实验室表征，由Gibney实验室提供），在葡萄糖限制恒化器中培养，并基于DO振荡在低DO（高呼吸相）和高DO（低呼吸相）阶段采集三个生物学重复样本。通过rRNA去除的全转录组RNA测序（WRS）和RNA-seq衍生的错配分析，研究人员发现低DO阶段与高DO阶段相比，存在1505个差异表达基因，且低DO富集核糖体生物合成、rRNA加工、翻译、硫代谢等，高DO富集氧化剂解毒、氧化还原酶活性和氧化还原缓冲通路。错配率分析显示，低DO相关错配响应具有替换类依赖性：合并的C>T/G>A类（包括互补的G>A事件）表现出最强的错配率增加（约1.074倍），同时T>C/A>G和T>A/A>T也显著增加，而C>A和T>G轻微下降。变异率分析进一步发现T>C/A>G和T>A/A>T的相对增加幅度最大（分别约13.8%和12.0%），但C>T/G>A类仍是最丰富的。复发分析显示大多数变异是样本特异性的，而非跨三个生物学重复稳定存在，支持瞬时的RNA-seq错配景观而非稳定的DNA突变。这些发现将YMC确立为研究代谢阶段、氧化还原状态和转录组保真度相互作用的简化真核框架，论文发表在《Antioxidants》。
**主要关键技术方法**
研究人员采用连续葡萄糖限制恒化器培养（菌株DBY12007，S288C背景，野生型HAP1
+），在低DO和高DO阶段采集三个生物学重复样本。主要技术包括：rRNA去除的全转录组RNA测序（WRS，使用RiboCop rRNA去除试剂盒和CORALL总RNA-seq V2试剂盒，Element AVITI24平台双端测序2×150 bp），STAR v2.7.11b比对至sacCer3参考基因组，featureCounts基因定量，DESeq2差异表达分析，clusterProfiler进行基因集富集分析（GSEA，基于GO生物学过程和分子功能）。错配分析使用GATK CallableLoci定义可调用位点，SAMtools mpileup提取碱基计数，beta-binomial广义线性混合模型（GLMM）评估错配率，负二项GLMM评估变异率（变异由Mutect2调用），复发分析通过交叉表卡方检验。所有分析在64位Linux工作站完成，使用R包glmmTMB、emmeans、ggplot2等。
**研究结果**
- **3.1 实验设计与数据质量**：通过连续葡萄糖限制恒化器培养，酵母菌株DBY12007在低DO和高DO阶段采样，DO振荡确认同步呼吸周期。WRS样本每个得到70-100百万唯一比对reads，主成分分析（PCA）显示低DO和高DO样本沿第一主成分分离，生物学重复紧密聚类，表明高重复性和干净阶段分配。
- **3.2 阶段特异性转录组**：DESeq2分析发现1505个基因在低DO和高DO间差异表达（调整p<0.05），|log2FC|>1阈值下487个低DO富集基因和574个高DO富集基因。GSEA显示低DO富集RNA加工与成熟、rRNA甲基化、tRNA转录与修饰、细胞质翻译、核糖体大亚基输出、硫代谢等；高DO富集醛和甲基乙二醛代谢、活性氧代谢、细胞氧化剂解毒、细胞对氧化应激和热响应、糖原代谢、线粒体呼吸链复合体组装。分子功能层面，高DO富集硫氧还蛋白依赖的过氧化物酶活性和谷胱甘肽过氧化物酶活性，低DO富集mRNA结合、RNA解旋酶活性、U3 snoRNA结合、rRNA结合、tRNA结合、甲基转移酶活性等。与独立YMC数据集（Murray等）高度一致（全局相关系数约0.78），支持阶段分配和表达架构的稳健性。
- **3.3 全基因组碱基替换分析**
- **3.3.1 错配水平分析**：基于可调用位点，每位点错配等位基因频率（MAF）的平均变化显示低DO相关响应具有替换类依赖性。C>T类（合并G>A）正向变化最大，而C>A和T>G下降。样本水平错配率（MR）通过beta-binomial GLMM估计，C>T/G>A兼容类在低DO中MR为1.519×10
-3，高DO为1.414×10
-3，倍数为1.074（显著增加）。C>G、T>A、T>C也显著增加，T>G和C>A轻微下降。
- **3.3.2 变异水平分析**：变异率（VR）通过负二项GLMM估计，C>T/G>A类在低DO中VR为1.826×10
-3，高DO为1.680×10
-3，增加约8.5%。T>C和T>A增加最大（约13.8%和12.0%），C>G和T>G无显著变化，C>A仅微弱增加且未达显著性。复发分析显示大多数变异为样本特异性，C>T类在低DO中单样本变异比例从78.2%略升至79.1%，三样本变异从6.2%降至5.4%（校正后p=0.0259），其他类无显著差异。
- **3.3.3 对比配对变异分析**：去除高DO和低DO共享的候选位点后，每样本变异率仍保持与主模型一致的模式，C>T/G>A类仍为优势信号，T>A和T>C保留低DO相关增加。拆解C>T类为C>T和G>A方向，各自占比约一半，条件间变化可忽略（Cramer's V=0.0009），说明碰撞类平衡稳定。
**讨论与结论**
讨论部分总结：YMC在时间上分离了生物合成与氧化还原缓冲程序，低DO代表高呼吸/高耗氧状态，但高DO富集的解毒通路并非意味着高DO更具氧化损伤性，而是反映恢复期。错配和变异率分析支持低DO相关替换类依赖的RNA-seq错配响应，其中C>T/G>A类最显著，但T>C/A>G和T>A/A>T增加表明机制复杂，不能归因于单一损伤（如8-氧鸟嘌呤或胞嘧啶脱氨）。复发分析支持瞬时、异质性景观而非稳定DNA突变。研究将YMC作为简化模型，提出代谢阶段可结构化转录组保真度，与氧化还原周期共振模型（Stolc等）和RNA-DNA差异（RDD）框架一致。未来需直接测量ROS、匹配DNA/RNA测序、直接RNA测序和化学鉴定来明确分子起源。结论部分翻译如下：
本研究证明酵母代谢周期在定义代谢阶段中同时组织转录本丰度和RNA-seq序列不一致。低DO/高呼吸期富集核糖体生物合成、RNA加工与修饰、翻译、硫代谢和蛋白质合成；高DO/低呼吸期富集氧化剂解毒、氧化还原酶活性和氧化还原缓冲相关通路，符合生物合成需求与细胞恢复的时序分离。错配率和候选变异率分析识别出低DO相关的替换类依赖性响应，其中合并的C>T/G>A兼容类错配率增加最大，而T>C/A>G和T>A/A>T在变异率分析中也增加。样本特异性候选位点占主导和复发差异小支持瞬态、异质的RNA-seq错配景观，而非广泛的稳定DNA突变。综上，结果支持氧化还原周期共振模型的核心预测，即序列变异可源于转录本合成和加工时的代谢状态。阶段依赖性和替换类依赖性的RNA-seq错配景观表明，代谢胁迫不仅是背景条件，更可结构化产生RDD兼容序列变异的生化环境。因此，YMC可作为简化真核模型，用于研究代谢计时、呼吸状态、氧化还原缓冲、RNA加工和转录组保真度如何相互作用。进一步研究需要直接ROS测量、匹配DNA/RNA测序、化学定量修饰或氧化核苷、直接RNA测序以及遗传或药理学扰动氧化还原通路，以确定观察到的错配模式的分子起源和功能后果。
