Genome Biology 项目文章| 范衡宇教授团队揭示ZAR1/2调控卵母细胞母源mRNA动态变化的新机制

2025年5月9日,浙江大学生命科学研究院范衡宇教授团队在Genome Biology杂志上发表了题为“ZAR1 and ZAR2 orchestrate the dynamics of maternal mRNA polyadenylation during mouse oocyte development”的研究论文。该研究联合Smart-seq2、Total RNA-seq、PAIso-seq2和LACE-seq四种近年来开发的、适用于低起始量样本的转录组检测技术,从mRNA稳定性和聚腺苷酸化状态两个维度重新探讨了卵母细胞在减数分裂过程中的转录本动态变化,并深入分析了ZAR1在母源转录组调控中的作用机制。希望组为本研究提供了PB转录组测序服务。

研究背景

在减数分裂期间,卵母细胞基因组会长时间处于转录沉默状态,直到合子基因组激活(ZGA)才重新启动。这一阶段,母源转录组的动态变化和稳态对母源-合子转换(MZT)至关重要。然而越来越多证据表明Smart-seq2(一种低输入量建库技术)检测到的母源转录组动态可能不够准确。其原因是该技术依赖oligo-d(T)引物捕获母源mRNA的poly(A)尾可能会引入偏差。

合子停滞蛋白1(ZAR1)及其同源蛋白ZAR2是RNA结合蛋白,对母源mRNA的稳定性调控具有重要作用。早期研究推测ZAR1促进减数分裂中母源mRNA的降解,而近期研究发现ZAR1参与线粒体相关核糖核蛋白域(MARDO)的组装,推测起到稳定母源mRNA的作用。因此关于RNA结合蛋白ZAR1的功能存在的差异也说明母源mRNA的调控机制仍需进一步阐明。

结果与分析

01 Smart-seq2检测的mRNA动态可能受多聚腺苷酸化偏差影响

在生发泡破裂(GVBD)后的减数分裂过程中,母源基因组转录保持沉默,多年来研究者观察到母源转录组逐渐减少的现象(图1)。这一时期发生降解的转录本对应的基因被称为O-decay基因。通过Smart-seq2和Smart-seq3测序结果显示,与GV期相比,MII期卵母细胞中有许多基因下调,而在差异基因表达中却存在许多上调基因(图1),这些结果显示通过Smart-seq2定义的O-decay基因可能存在偏差。为探究这一问题,该研究采用PAIso-seq2分析技术检测减数分裂过程中多聚腺苷酸化状态的变化,发现GV到MII期卵母细胞中,大部分转录本poly(A)尾明显缩短,这可能导致Smart-seq2无法准确检测短poly(A)尾的mRNA。

图1 Smart-seq2和PAIso-seq2检测母源mRNA动态变化

为规避Smart-seq2的系统偏差,该研究使用Total RNA-seq数据进行分析,结果显示Total RNA-seq检测到的母源转录组在减数分裂过程中的整体下调幅度显著低于Smart-seq2的结果,且鉴定的O-decay基因数量明显更少(图2)。对Total RNA-seq鉴定的O-decay基因进行了深入分析,发现O-decay基因的转录本主要在GV期卵母细胞中发生多聚腺苷酸化和翻译,在完成生物学功能后被有序降解。因此研究推测在减数分裂过程中,母源转录组经历的多聚腺苷酸化/去腺苷酸化修饰可能比实际降解更为显著。

2 Total RNA-seq检测母源mRNA动态变化

02 Zar1/2-/-导致GV期母源mRNA稳定性下降及多聚腺苷酸化异常

早期研究推测ZAR1促进减数分裂中母源mRNA的降解,而近期研究发现ZAR1参与线粒体相关核糖核蛋白域(MARDO)的组装。为阐明这一矛盾,研究采用总RNA-seq技术对Zar1/2-/-卵母细胞进行差异表达基因分析,结果检测到GV期母源mRNA显著下调,且Zar1/2-/-GV期下调基因中61.76%与正常GV→MII期应下调或持稳的基因相关,说明Zar1/2在GV期通过稳定mRNA维持其表达。

通过对Smart-seq2数据的重新分析,研究发现Zar1/2-/- MII期卵母细胞中仍存在大量RNA转录本的异常积累现象,这与总RNA-seq的检测结果存在显著差异。为验证这一差异,研究发现Zar1/2-/-卵母细胞中的多聚腺苷酸化过程可能发生紊乱,PAIso-seq2显示大量mRNA poly(A)尾长度异常,这些异常的poly(A)尾长度导致了Smart-seq2检测偏差。
 
03 ZAR1通过结合3’UTR稳定母源转录本

ZAR1作为RNA结合蛋白(RBP),通过其C端CxxC锌指结构域与母源转录本结合。基于这一分子特征,该研究深入探究了ZAR1的调控机制。LACE-seq定位到其靶向转录本,共鉴定8,000余个潜在靶标,其中Zar1/2-/-GV期下调基因中54.8%为ZAR1靶标。此外发现3’UTR结合基因的稳定性显著高于CDS结合基因,提示ZAR1通过3’UTR结合维持mRNA稳态(图3)。

图3 ZAR1通过与3‘UTR区域结合来稳定母体转录本

04 ZAR1间接调控多聚腺苷酸化并与蛋白互作

鉴于Zar1/2-/-卵母细胞中母源转录本多聚腺苷酸化显著异常,研究进一步分析了差异多聚腺苷酸化基因(DPGs)与ZAR1靶标的关系,分析发现ZAR1本身并不直接调控多聚腺苷酸化,ZAR1更倾向于影响mRNA稳定性,而非直接影响poly(A)尾长。研究通过IP-MS技术发现ZAR1与PABP家族、IGF2BP2等RNA稳定因子互作,提示ZAR1可能与其他蛋白质相互作用以调节母源转录组的稳态(图4)。

图4  ZAR1在卵母细胞发育过程中调控mRNA多聚腺苷酸化的工作模型

05 Zar1/2-/-导致MII期染色质压缩维持失败

Zar1/2-/-卵母细胞在MII期染色质压缩维持失败,出现纺锤体错位和类原核结构等现象,伴随关键母源基因(如Lsm14b、Ccnb1)表达异常,表明ZAR1在稳定这些基因中起作用,从而有助于维持卵母细胞处于MII期。

综上所述,该研究通过多组学联合分析,揭示了聚腺苷酸化状态在减数分裂过程中主导母源转录组动态变化的关键作用。研究进一步阐明了ZAR1在卵母细胞减数分裂成熟早期维持母源转录组稳定性,以及在成熟过程中协同其他蛋白调控母源转录本聚腺苷酸化状态的功能。该研究不仅为Smart-seq2技术在卵母细胞中的的适用性提供了参考,也为进一步探索ZAR1蛋白相关的MARDO在减数分裂过程中的作用机制提供了参考资源。
浙江大学生命科学研究院博士生吴雨珂为本文第一作者。广东省第二人民医院的苏瑞宝副研究员,已毕业博士生蒋知妍和吴韵雯,浙江大学医学院附属邵逸夫医院的戎妍副研究员,浙江大学医学院的嵇姝妍研究员,中国科学院遗传与发育研究所陆发隆研究员及其课题组成员博士生刘静雯和牛卓越,以及中国科学院生物物理研究所薛愿超研究员为该研究提供了大力帮助。该研究受到国家重点研发计划、国家自然科学基金以及浙江省自然科学基金项目的资助。

文章链接:https://doi.org/10.1186/s13059-025-03593-8