《Molecular Reproduction and Development》：Sex-Specific Gene Expression Ontogeny During Gonadal Development in Post-Metamorphic Xenopus tropicalis
编辑推荐：
两栖动物性腺分化是生殖发育研究的关键环节，但变态后的成熟过程与性别特异性转录程序仍不清晰。为此，本研究整合形态学、组织学和基因表达分析，对非洲爪蟾（Xenopus tropicalis）变态后1-8周的幼蛙性腺进行追踪。研究发现，睾丸生长与体型正相关，而卵巢发育则独立于体细胞参数；基因表达谱显示cyp17、amh、amhr2与睾丸发育呈正相关，aldh1a2、ddx4则在卵巢生长和卵泡数量增加中显著上调。此外，id4在雌性下降、在雄性稳定，反映了雌雄在配子前阶段和生殖系维持的不同策略。多元分析与网络分析进一步指出dmrt1、amhr2、cyp17、esr1、ddx4是性别特异性调控网络的核心节点。该工作为脊椎动物性别分化研究提供了整合的组织学与分子终点，对理解性腺发育的分子机制具有重要意义。
在脊椎动物的生命历程中，性腺分化是一个极为关键的发育过程，它涉及从未分化性腺的性别决定，到性别特异性结构与分子特征的出现。两栖动物，尤其是非洲爪蟾（Xenopus tropicalis），因其易饲养、全年可繁殖以及相对较短的世代周期，成为研究性别决定与分化的经典模型。然而，尽管我们对幼虫期的性腺分化已有较多了解，但在蝌蚪经历变态成为幼蛙后，性腺如何继续成熟、性别特异性的转录程序如何动态变化，仍是一个未得到充分阐释的“黑箱”。这一时期——即前青春期（pre-pubertal）的发育窗口——被认为是性功能成熟前的关键过渡阶段，但其分子与组织学事件的精确时间线依然模糊。为了填补这一空白，一项发表在《Molecular Reproduction and Development》上的研究，对X. tropicalis 变态后1至8周的幼蛙进行了系统性探究，旨在揭示性腺成熟与性别特异性基因表达之间的内在联系。
研究团队综合运用了形态测量、组织学评估和基因表达分析（实时荧光定量聚合酶链式反应，qPCR）等技术手段。他们每周采集5-10只幼蛙，依据年龄划分为四个组：Met 1-2（变态后1-2周）、Met 3-4、Met 5-6和Met 7-8。在解剖时，测量了体重、体长、肝脏和脂肪体重等形态参数。左侧的性腺-肾脏复合体（Gonad-Kidney Complex, GKC）用于石蜡包埋、切片和HE（苏木精-伊红）染色，以评估性腺面积和成熟度（睾丸采用1-7分的成熟度评分系统，卵巢则计数卵泡卵母细胞数量）。右侧GKC则用于RNA提取、cDNA合成和qPCR分析，检测了涉及性别分化、类固醇生成、生殖细胞和视黄酸通路的12个候选基因的表达水平。数据分析包括广义线性模型、主成分分析、贝叶斯网络分析和斯皮尔曼相关性分析，以揭示基因表达与性腺发育之间的关系。
3.1 形态学数据
雄性幼蛙的体重和体长随年龄有增加趋势，但未达统计学显著性。肝脏和脂肪体体细胞指数在Met 3-4时升高，随后在Met 7-8时显著下降。雌性幼蛙的体重、体长和脂肪体指数在Met 5-6时达到峰值，表明其能量储存策略与雄性不同。
3.2 性腺的组织学评估
雄性睾丸的面积和成熟度评分均随年龄显著增加，并且与蛙的年龄呈显著正相关。睾丸发育与体重、体长、肢体长度等体细胞生长参数也呈正相关。雌性卵巢的面积和卵泡卵母细胞数量同样随年龄增加，但与任何体细胞测量参数均无显著相关性，提示卵巢成熟独立于整体身体生长。
3.3 基因表达
基因表达分析揭示了显著的性别和阶段依赖性模式。
- •
雄性偏向表达：cyp17、amh 和 amhr2 在雄性中表达量更高，并与睾丸发育正相关。cyp17 表达在雄性中随时间有上升趋势，反映了其对雄激素合成和睾丸生长的支持作用。
- •
雌性偏向表达：aldh1a2 和 ddx4 在雌性中表达量更高。aldh1a2（视黄酸合成酶）在雌性中从Met 5-6开始显著上调，与卵巢生长和卵泡数量正相关，提示其参与促进卵母细胞减数分裂进程。ddx4（生殖细胞标志物）在雌性中表达始终高于雄性，并与两性性腺大小和成熟度均呈强正相关，反映了生殖细胞数量的增加。
- •
性别特异性动态：id4（维持干细胞状态的因子）在雌性中有下降趋势，而在雄性中保持稳定，暗示了雌雄在配子前阶段的不同策略。esr1（雌激素受体1）在雄性后期（Met 7-8）表达更高。cyp19（芳香化酶）和sox9未表现出显著的性别二态性。cyp26b1（视黄酸降解酶）在雄性中略有上升，并与睾丸成熟度正相关，与aldh1a2的表达呈显著负相关，支持了“在雌性中促进减数分裂，在雄性中抑制减数分裂”的经典模型。
3.4 主成分分析
主成分分析清晰地将雄性和雌性个体在基因表达谱上区分开。与雌性关联较强的基因是ddx4和aldh1a2，而与雄性关联较强的基因包括cyp17、dmrt1、amh、amhr2、esr1、id4和cyp26b1。
3.5 基因表达与性腺发育的相关性
在雄性中，ddx4、dmrt1、cyp17、amhr2和cyp26b1的表达与睾丸面积和成熟度呈正相关。在雌性中，ddx4和aldh1a2与卵巢面积和卵泡数量正相关，而amhr2和3βhsd的表达与卵巢面积呈负相关。
3.6 贝叶斯网络分析
网络分析揭示了基因表达间存在依赖性的关联网络，且网络结构因发育阶段和性别而异。例如，在早期幼体（Met 1-4）中，ddx4与dmrt1正相关，与id4负相关。在后期（Met 5-8），dmrt1与amhr2、aldh1a2与cyp19呈现正相关。在雌性网络中，id4与ddx4呈强负相关，而ddx4与aldh1a2呈正相关。
3.7 基因表达水平间的相关性
跨所有样本的相关性分析显示，amh与cyp17、esr1、id4等雄性相关基因强正相关。aldh1a2与ddx4正相关，但与cyp26b1强负相关。ddx4与3βhsd负相关。esr1与id4、sox9等基因正相关。
研究结论与意义
本研究的核心结论是，在非洲爪蟾变态后的前两个月内，性腺发育在分子和组织学层面均展现出强烈的性别特异性程序。雄性睾丸的发育与体细胞生长协调一致，而雌性卵巢的成熟则相对独立。在分子层面，cyp17、amh、amhr2构成了支持睾丸发育的核心雄性基因模块；aldh1a2和ddx4的上调则是卵巢生长和卵泡发生的关键特征。id4和ddx4在雌雄中的相反动态，精准地反映了雌性从卵原细胞库向减数分裂卵母细胞过渡，以及雄性维持精原干细胞池的不同生殖策略。视黄酸通路基因（aldh1a2 vs. cyp26b1）的相反调控模式，为理解两性减数分裂启动的差异提供了分子基础。
通过多元分析和网络建模，研究进一步鉴定出dmrt1、amhr2、cyp17、esr1和ddx4等基因是潜在的性别特异性调控网络节点。这项工作不仅首次系统地描绘了X. tropicalis 前青春期性腺发育的基因表达谱系，将组织学成熟与特定分子标志物相关联，而且为利用该模式生物进行脊椎动物性别分化、环境内分泌干扰物评估以及生殖毒理学研究，提供了宝贵的整合性分子与组织学终点参考。
