清点一只卵巢中的卵母细胞,此前多半要把它切成连续切片再作外推。通过对整个器官成像、并以深度学习逐一分割每一枚卵母细胞,巴塞罗那基因组调控中心的团队完成了小鼠整个生育期内卵巢储备的清点——并从中发现了一种出人意料的不变性:任一时刻处于激活窗口内的储备比例纹丝不动,而储备本身却在崩塌。
来源:biorxiv.org
简明说来
女性出生时便带着一份此生不再更新的卵子储备,这份储备随年龄而消减。人们很容易把这座蓄水池想象成一只正在排空的浴缸:剩下的水越多,流量就越大。而这项在小鼠身上完成的研究观察到的却是另一回事。任一时刻,剩余储备中同样的比例——约 14 %——正处在卵母细胞脱离休眠的过渡阶段。即便储备减少到不足原先的十分之一,这一比例也不改变。换言之,受到调控的似乎是一个比例,而不是一个数目。开篇先立三条保留:这一切都是在小鼠身上测得的;研究确认了一种规律性,却没有以实验检验产生这种规律性的机制;这里既没有疗法,也没有临床承诺。
| 参数 | 数值 |
|---|---|
| 日期 | 预印本于 2025 年 11 月 8 日上线;同行评议版于 2026 年 8 月 12 日在 Nature Aging 公布 |
| 团队 | Arturo D'Angelo(第一作者)、Elvan Böke(末位作者)——CRG 巴塞罗那,并与 EMBL 巴塞罗那、MRC LMB 剑桥、DIPC/UPV-EHU 合作 |
| 模型 | 近交系 C57BL/6J 小鼠——56 只未交配雌鼠,5 至 60 周共 7 个周龄,每个周龄至少 5 只;另有一组交配队列在 31 周作对照。已在人卵巢皮质上验证可行性(2 名 25 岁供体) |
| 成像 | 对整只透明化卵巢(BABB)作光片显微成像(SPIM),DDX4 标记(卵母细胞)+ TO-PRO-3(细胞核) |
| 分割 | 基于深度学习的实例模型(BiaPy 框架)——与专家标注一致率 95,3 %,经轻度人工校正后为 96,7 % |
| 数据量 | 在约一百只完整卵巢中分割出八万五千余枚卵母细胞 |
| 关键结果 | 各周龄下处于过渡期的卵母细胞比例稳定在储备的约 14 %,而储备总量下降超过一个数量级 |
| 拟合的动力学 | 激活 k1=0,0258±0,0053 周⁻¹,即每周约有 2,5 % 的休眠储备被动员 |
| 变异度 | 同一周龄、基因型完全相同的小鼠之间,储备总量相差最多达 3 倍(5 周龄时 n = 8) |
| 状态 | bioRxiv 预印本,上线时尚未经同行评议;评议版由机构新闻稿公布 |
技术解析
卵巢为何难以清点 —— 小鼠卵巢有数毫米大小,形状并不规则。传统方法是把它切成一系列组织学切片,在其中一部分切片上清点卵泡,再乘以一个校正系数。作者指出,这类基于切片的方法在估计原始卵母细胞数目时会引入误差,并且牺牲了全部空间信息。作为方法学上的补充可以指出,其中通常潜伏着两个误差来源:被两张切片同时切到的卵母细胞可能被数两次,也可能一次都没数到;而外推则预设了分布是均匀的——恰恰是这项研究所动摇的前提。光片显微镜以一束与观察轴垂直的薄平面照亮样品,从而绕开了这一障碍,能够在不触碰组织的情况下获取光学切片。但光还得能穿透:透明化把组织的介质替换为折射率接近蛋白质的混合液,此处是苯甲醇/苯甲酸苄酯,从而消除界面处的散射,使器官变得半透明。
标记物 —— 抗 DDX4 抗体用作卵母细胞标记物;染料 TO-PRO-3 对 DNA 作复染并提供组织背景,从而能在其余细胞核之中定位卵母细胞。作者表示,其优化后的方案还能区分卵泡发生的各个阶段——依据的是大小这一判据,详见第 4 点。他们主张,在免疫标记与透明化方案共同优化之后,实现了对全部卵母细胞的原位鉴定;而这一主张所依据的正是这项优化,而非两个通道之间的交叉核对。
分割,及其一致率究竟界定了什么 —— 一只完整卵巢中含有数千枚彼此相贴的卵母细胞;把它们逐一分开属于实例分割,比单纯的检测更困难。作者手工标注了体数据,在 BiaPy 框架中训练了一个网络,随后对五只独立卵巢成像,将自动输出与专家标注作比较:校正前一致率 95,3 %,轻度校正后 96,7 %。关于这一数字的适用范围需要作一点说明:该验证是在 5 周龄卵巢上完成的,也就是密度最高的状态。没有任何依据表明它在 60 周龄仍然成立——那时局部密度已降至四分之一。因此,随周龄而变的性能漂移仍有可能存在,而它恰恰会影响对衰老本身的测量。
卵母细胞如何归类——一个以形态标定的直径阈值 —— 这是读懂主要结果最关键的一点,而且很容易说岔。作者是在共聚焦显微镜下、而不是在 SPIM 体数据中去寻找颗粒细胞的形态;主流程本身测量的是卵母细胞的直径。于是他们在这一共聚焦成像的子集上标定阈值——在那里,扁平细胞向立方细胞的转变(标志着进入生长期)可以直接观察到。取相应直径的平均值 ± 标准差,他们得到 16,30 µm 与 20,11 µm 两个边界,随后仅按大小对整个数据集分类:10 至 16,30 µm 为原始,16,30 至 20,11 µm 为过渡期,20,11 µm 以上为生长期。承载核心结果的"过渡期"这一类别,因此是一个狭窄的直径窗口,它以形态学观察为依托,却是按度量方式施用的。
核心结果:一个比例,而非一个数目 —— 从 5 周到 60 周,储备总量下降超过一个数量级。同一时期,这一池子的构成也明显移动:原始卵母细胞从约占总数的 70 % 降至 55 %,生长期卵母细胞则走了相反的路径,从约 15 % 升至 30 %。正是这一反差使第三个数字显得引人注目——处于过渡期的卵母细胞比例在各周龄下始终维持在 14 % 左右。以 R(t) 记储备、T(t) 记过渡期数目,测量给出 T(t)/R(t)≈0,14;而按绝对速率恒定的动员则要求 T(t)≈const,于是随着 R 减小,该比值将急剧上升。数据在两种机制之间作出了裁决,支持按比例动员。注意不要混淆这两个量:14 % 是某一时刻的数目占比,而不是速度。动员速率本身估计为每周约占休眠储备的 2,5 %(第 6 点)。
动力学模型 —— 该进程被描述为单向级联:已脱离休眠的卵母细胞无法再回到休眠;在每一步,它或以速率 ki 前进到下一个隔室,或以速率 di 死亡。预印本并未显式写出方程;该系统可以写成:
dtdP1=−(k1+d1)P1dtdPi=ki−1Pi−1−(ki+di)Pi (i=2,3)dtdP4=k3P3−d4P4
仅仅为了这次拟合,大小分布被重新划分为四个互不重叠的类别:原始(10–20 µm)、初级(> 20–40 µm)、次级(> 40–60 µm)、窦状(> 60 µm)。这一四隔室划分与第 4 点的三阶段方案并不相同:"过渡期"这一类别并不出现在其中。要留意术语在两套方案之间改变了含义——模型中的"原始"隔室(10–20 µm)把第 4 点的过渡窗口(16,30–20,11 µm)包含在内。因此上文所说的每周 2,5 %,是这个被扩大了的隔室的流出率,而不是严格意义上休眠池的流出率。该系统以 5 周龄测得的数目为初始条件积分,并以带约束的最小二乘法(SciPy 的 SLSQP 算法)拟合。所得速率为 k1=0,0258±0,0053、k2=0,2173±0,0400、k3=0,3375±0,0567 周⁻¹;损失一侧为 d1=0,0179±0,0074 与 d4=0,7224±0,1470 周⁻¹——后者是窦状阶段的流出率,把退化与排卵合并在一起。有一点读法值得点明,因为它改变了两个参数的性质。作者在方法部分指出,带约束的极小化使初级与次级阶段的死亡率归零,而由于文献也指向这一方向,他们随后在最终拟合中把 d2 与 d3 固定为零。因此所显示的"±0,0000"并不是测量精度:它是被冻结的参数留下的印记,而其余速率都是在这一选择既定的前提下拟合出来的。作者从中看到自身假设稳健性的佐证;人们同样可以看到一个被注入模型的假设——而人们再也无法要求模型去检验它了。
60 µm 的瓶颈 —— 一旦被激活,卵母细胞的体积会增大约 64 倍——这个数字作者取自文献,并非在此测得。倘若激活与生长都是匀速的,大小分布应当形成一条平滑的连续曲线。实际上它是双峰的:一个位于 15 µm 附近、狭窄而占主导的峰,对应休眠的储备,随后是 60 µm 附近的第二个峰。大小分布中的峰意味着放缓——细胞在推进最慢之处堆积,正如队伍会在链条最慢的环节排起。而 60 µm 恰好对应次级晚期/窦状早期的过渡,也就是卵泡开始对促性腺激素产生反应的时刻。因此所提出的解释是:随着这种激素敏感性的获得,生长进入了另一种状态。该峰在各周龄下都持续存在,说明这一现象并不取决于剩余储备的大小。
几何:卵母细胞待在哪里 —— 一条样条曲线被人工放置于器官之上,充当内部标尺:每一枚卵母细胞都以其位置及围绕该轴的角度来定位,从而使不同个体、不同周龄的卵巢之间可以比较。因此这一参照系是角度性的,而非纵向的:0° 指上表面,180° 指卵巢门,45–135° 与 225–315° 两个扇区则为中外侧区域。结果十分清晰:原始卵母细胞在卵巢门(凹面)与相对的凸极都偏少,而在中外侧区域偏多。同样的偏倚也适用于过渡期与生长期的卵母细胞,并随年龄增长而保持。样条的人工放置带入了一部分主观性,值得留心,因为全部角度定位都依赖于它。
局部密度与激活 —— 这是与上一点不同的另一项分析:对每一枚原始卵母细胞,作者清点其 100 µm 半径内的原始卵母细胞邻居。这一局部密度随年龄下降——5 周龄时平均超过 10 个邻居,60 周龄时约为 2,5 个。有意思的是相关的方向。一个由来已久的假说认为,密集的邻里会抑制激活,休眠卵泡彼此压制。数据却指向相反方向:原始密度高的区域反而容纳更多处于过渡期的卵母细胞,而反向的关系同样为正。必须按其本来面目来读这一结果——它是一种空间相关,是在未施加扰动的情况下测得的。它动摇了邻域抑制模型(诚然,该模型是针对青春期前卵巢提出的,而这些数据自 5 周龄才开始);但它并不证明高密度会导致激活。某个第三变量,例如局部血管化,可能同时支撑这两种效应。
以妊娠作检验 —— 雌鼠自 7 周龄起交配,历时 18 周,完成多次妊娠;其储备在 31 周龄时与同龄未交配个体作比较。无论在总量还是在任一子类别上,都未出现显著差异。作者由此推断,抑制排卵似乎并不减缓卵母细胞的流失,而卵泡动员应当是在组织层面受调控的,与排卵活动无关。有两点读法上的分寸:这是在单一测量时点上得到的无差异结果,并不等于证明;而且该干预并未把排卵单独分离出来,因为反复妊娠同样会改变激素轴、代谢与血管化。
Why It Worked
其贡献与其说在于某项前所未有的技术,不如说在于三个要素的组合,而每一个都弥补了另外两个的短处。体积成像消除了切片外推带来的偏倚;自动分割使得在上百只器官上作穷尽清点成为可能,而在这一规模上人工清点是行不通的;对照专家标注的验证则界定了算法的误差——至少在其所进行的那一状态下。作者自己也把这项工作置于一条脉络之中:此前已有两项近期研究把透明化、显微镜与机器学习结合起来,以自动化卵母细胞的清点。他们主张为新的是那种普适性——一套在成年卵巢上、在整个生育期内都站得住脚,并且可以移植到人体组织的流程。
最能说明问题的比较,不是方法的前后对照,而是关于卵巢排空的两种竞争模型之间的对照。按绝对速率恒定的动员预言每周被激活的卵母细胞数目固定;按比例的动员则预言比例固定。数据作出了裁决:在储备减少到不足十分之一、而池中另外两个分数各自移动约十五个百分点的同时,这一比例始终维持在 14 % 左右。
剩下的是测得之物与所断言之物之间的落差。广为流传的"卵巢在清点自己的卵子"这一说法,并非仅仅出自机构新闻稿:它是由作者本人承担的,其讨论部分提出"卵巢拥有一种机制,使其能够掌握卵母细胞的总数",而摘要则归结为一种器官尺度的调控机制。预印本以条件语气写下这些句子,而校准恰恰就在此处见分晓:所测得的是比例的不变性,以及局部密度与激活之间的相关。存在某种"读取"储备大小并据以唤醒一定百分比的机制,这是与上述观察相容的假说,但这里没有任何实验去扰动储备的大小以检验它:既没有部分切除之后再测激活速率,也没有突变体,也没有移植。这项研究是描述性和相关性的,而它的坚实之处也正在于此。
另有三项限制值得点明。初级与次级阶段的零死亡率并不是一项测量,而是作者所固定的一个参数。模型还假定生长速率与死亡速率在衰老过程中大体保持稳定——这是作者明确声明、并认为有待修正的假设,尽管所有已发表的速率都以它为条件。最后,人体方面的验证仅限于表明该流程在两名 25 岁供体的卵巢皮质上可以运作:这是技术可行性的证明,绝非对人体动态的研究。
Causal Chain
卵母细胞储备在胎儿期即已定下 → 传统上以组织学切片与外推清点,代价是原始卵母细胞计数的误差与空间信息的丧失 → 光片显微镜与折射率匹配透明化方案的发展 → DDX4 标记与细胞核复染,使每一枚卵母细胞得以在完整组织中定位 → 三维体数据的人工标注与实例分割模型的训练 → 在五只 5 周龄卵巢上对照专家取得 95,3 % 的验证 → 依据共聚焦下观察到的颗粒细胞形态标定直径阈值(16,30–20,11 µm)→ 在跨越七个周龄的 56 只小鼠中逐一清点八万五千余枚卵母细胞 → 确认在储备崩落一个数量级、另外两个分数发生移动的情况下,过渡期比例仍不变地维持在约 14 % → 拟合微分方程组,得出每周约 2,5 % 的休眠储备被动员 → 把卵巢储备重新表述为受调控的系统而非被动的蓄库 → 提出器官尺度按比例调控的假说,至今尚未经扰动检验 → 一旦获得证实,或可成为调节排空速度的靶点,但视野完全停留在临床前阶段。
Anecdote
这项研究中最令人不安的细节并不是主要结果,而是先于它的那项观察。正是抗缪勒管激素——人类临床上使用的储备标志物,其数值在不同女性之间差异极大——已知的变异度,促使作者去追问:卵母细胞的数目本身在个体之间是否也变化得这么大,哪怕是在基因型完全相同的群体内部?于是他们对八只 5 周龄的 C57BL/6J 雌鼠成像,这些都是近交系动物,因而基因型完全相同。答案是肯定的,而且差距相当可观:一只的卵母细胞最多可达另一只的三倍。这一差距在五周龄时就已可测得,提示它在性成熟之前便已确立。虽然作者并未写下这一点,但可以由此推想:初始禀赋中有一部分既不取决于遗传,也不取决于生育史。
Legacy and Current Data
作者所回顾的人类数量级勾勒出了问题所在:出生时约有 100 万至 200 万枚卵母细胞,青春期约 40 万枚,绝经时不足 1000 枚——而一生中排出的仅约 400 枚。因此,储备中绝大部分从未排卵便已流失;决定生育期长短的是这一流失的速度,而不是排卵的次数。
作者把此处提出的定量框架呈现为哺乳动物中卵泡激活速率与生长速率的首批数值估计。这些速率在女性身上究竟如何,则完全有待确立:研究中没有任何数据允许把它们直接移用。
两点透明性说明。预印本在数据量上存在内部不一致:其摘要称"八万五千余枚卵母细胞",其讨论部分则称"十万余枚"。两个数字出现在同一份文件中而未加调和;本文采用摘要中的数值八万五千,机构新闻稿沿用的也是这一数值。此外,Nature Aging 所刊版本的确切出处未能直接核实:此处所引全部数字均来自 bioRxiv 预印本,而非期刊版本。最后,这项研究是横断面的——七个不同周龄、不同个体——而非纵向的:现有技术无法对同一只卵巢作跨时间的追踪。5 周龄这一起点被选为激活与卵母细胞死亡第一波之后的最大储备,而不是生育期的真正开端。
研究者的视角 —— 未决的问题
以下实验并未出现在这项研究中;它们指出的是:要把观察到的相关变成已被证明的机制,还需要做些什么。
- 部分切除,按比例调控的决定性检验。 移除储备中已知的一部分,再测量剩余组织中处于过渡期的比例。若调控确实是按比例的,这一比例应在重新调整后回到 14 %;若是按绝对数量的,被激活的数目则应保持不变。正是这项实验能在两种模型之间作出裁决,而单侧卵巢切除提供了它的一个简便版本。
- 把密度与位置解耦。 研究中没有任何环节把密度—激活相关与中外侧富集区分开来,两项分析也从未交叉。一项控制位置的分析将能说明,两个变量中究竟是哪一个承担了效应。
- 以激素敏感性检验 60 µm 瓶颈。 若这一大小上的峰确实反映促性腺激素反应性的获得,那么操纵 FSH 信号通路应当使该峰在大小上移动,或者使其消失。一个不为所动的峰将推翻这一解释。
- 解除模型的平稳性假设。 各速率是在假定它们不随年龄漂移的前提下拟合的。按年龄区间分段拟合模型,再比较所得的 ki,就能知道这一假设是否站得住——这是最容易着手的内部对照。
- 把流程推广到人体系列。 可行性已在两名供体身上确立。一个跨越数十年的队列将能说明,比例的不变性究竟是小鼠所特有,还是哺乳动物中普遍存在:这正是决定整座建筑临床意义的问题。
来源
- D'Angelo A., Franco-Barranco D., Musy M., Sharpe J., Arganda-Carreras I., Böke E.,《3D Mapping of Intact Ovaries Reveals the Aging Dynamics of the Ovarian Reserve》,bioRxiv,2025 年 11 月 8 日上线 —— doi:10.1101/2025.11.07.686728。状态:预印本,上线时尚未经同行评议;评议版在 Nature Aging 公布。
- 基因组调控中心(CRG),《First 3D map of ovary throughout reproductive lifespan reveals mouse reproductive organ counts its own eggs》,2026 年 8 月 12 日新闻稿 —— crg.eu,亦由 EurekAlert! 发布。状态:机构新闻稿,公布在 Nature Aging 的刊出。
背景文献
以下文献并不报告所评述研究的结果;三者均见于该研究的参考文献,用以支撑上文所作的一般机制说明。
- Huisken J., Swoger J., Del Bene F., Wittbrodt J., Stelzer E. H. K.,《Optical sectioning deep inside live embryos by selective plane illumination microscopy》,Science,2004 年,第 305 卷,第 1007–1009 页 —— doi:10.1126/science.1100035。光片显微镜的原理。
- Edson M. A., Nagaraja A. K., Matzuk M. M.,《The mammalian ovary from genesis to revelation》,Endocrine Reviews,2009 年,第 30 卷第 6 期,第 624–712 页 —— doi:10.1210/er.2009-0012。卵泡发生,以及次级晚期/窦状早期促性腺激素敏感性的获得。
- Winship A. L., Sarma U. C., Alesi L. R., Hutt K. J.,《Accurate Follicle Enumeration in Adult Mouse Ovaries》,Journal of Visualized Experiments,2020 年 —— doi:10.3791/61782。连续切片法在估计原始卵泡数目时的偏倚。
