视黄酸促进黄鳍鲷精巢卵巢中减数分裂与生殖细胞分化程序

《Aquaculture Reports》：Retinoic acid promotes meiotic and germ cell differentiation programs in protandrous yellowfin seabream ( Acanthopagrus latus)

　　

视黄酸（retinoic acid, RA）是生殖细胞分化的保守调节因子，但其在具有解剖学上分离的卵巢区和精巢区的雄性先熟雌雄同体鱼类中的作用仍不清楚。本研究对黄鳍鲷（Acanthopagrus latus）幼鱼进行RA处理，并在处理后第30、60和90天检测性腺发育情况。组织学分析显示，在第90天时，对照组和RA处理组鱼的精巢区均出现初级精母细胞，而仅在RA处理组的卵巢区检测到少量初级卵母细胞，表明RA处理优先促进了卵巢生殖细胞的进程。转录组学分析在第30、60和90天的对照组与RA处理组性腺之间分别鉴定出1276、8172和989个差异表达基因（differentially expressed genes, DEGs），表明存在强烈的阶段特异性反应，其中第60天时重编程最为广泛。RA处理诱导了类视黄醇、钙离子和MAPK相关通路的早期激活，随后启动了神经内分泌、类固醇生成和减数分裂程序。在第60天时，RA处理组鱼的关键减数分裂基因（sycp3、dmc1、dazl和rec8）和性别分化基因（foxl2、dmrt1、sox9a和cyp19a1a）上调，支持生殖发育程序的更早分子激活。同时，RA降解酶（cyp26a1和cyp26b1）在第30天被强烈诱导，而aldh1a2在对照组中表达更高，提示RA可用性存在反馈调节。热图和蛋白互作网络（protein-protein interaction, PPI）分析进一步将RA代谢与类固醇生成、减数分裂进程和性别分化联系起来。这些发现表明，RA处理促进了黄鳍鲷精巢卵巢中的生殖细胞分化并提前了生殖发育程序，其中卵巢谱系表现出最明显的组织学响应，并提示RA生物活性在精巢卵巢发育过程中受到代谢反馈检查点的限制。

研究背景与科学问题

视黄酸（retinoic acid, RA）是维生素A的活性衍生物，作为保守的形态发生素调控脊椎动物的细胞分化、模式形成和器官发生。在生殖组织中，RA作为性腺环境信号连接生殖细胞命运与减数分裂起始时间。哺乳动物研究中已明确RA在促进精原细胞分化及胎儿卵巢减数分裂启动中的关键作用，其合成酶ALDH1A与降解酶CYP26B1的平衡决定了性腺局部RA浓度，进而影响减数分裂的性别二态性时序。然而，在硬骨鱼类中，RA的作用机制呈现物种特异性差异：斑马鱼缺乏STRA8但仍保留MEIOSIN直系同源物介导的不依赖STRA8的减数分裂活性；青鳉中RA对首次减数分裂启动作用有限；而在南方鲶中RA则激活经典的STRA8依赖性减数分裂通路。黄鳍鲷（Acanthopagrus latus）作为一种经济价值重要的鲷科鱼类和雄性先熟雌雄同体模型，其性腺同时发育出卵巢区和精巢区，为评估同一内分泌信号如何在同一生理环境中作用于两条生殖细胞谱系提供了独特机会。然而，RA处理是否能同步加速两条谱系的减数分裂，以及RA代谢如何与精巢卵巢中的性别分化程序整合，尚不清楚。

研究设计与主要结论

研究人员对黄鳍鲷幼鱼（体重0.99±0.25 g，来源：珠海长丰水产科技有限公司）进行为期90天的RA饲喂处理（300 mg/kg饲料），分别在处理后第30、60和90天采样，结合组织学分析、RNA测序、短时间序列表达挖掘器（Short Time series Expression Miner, STEM）聚类、热图分析、蛋白质-蛋白质相互作用（protein-protein interaction, PPI）网络分析和实时荧光定量PCR（RT-qPCR）验证，系统研究了RA对精巢卵巢发育的影响。主要结论包括：（1）RA处理促进了生殖细胞分化，尤其优先推动卵巢谱系向初级卵母细胞阶段进展；（2）转录组分析揭示了强烈的阶段特异性响应，第60天时重编程最为广泛，涉及1276、8172和989个差异表达基因（differentially expressed genes, DEGs）；（3）RA处理诱导了从早期信号通路（类视黄醇、钙离子和MAPK相关通路）到后期神经内分泌、类固醇生成和减数分裂程序的顺序性转录激活；（4）RA生物活性受到由CYP26酶和ALDH1A2组成的代谢反馈检查点的严格约束。这些发现为理解雌雄同体硬骨鱼类中RA调控生殖细胞分化的分子机制提供了新见解，并对水产养殖中生殖调控技术开发具有潜在应用价值。

关键技术方法

1. 1.
   实验设计：200尾黄鳍鲷幼鱼随机分为对照组和RA处理组（每组三个重复缸），RA处理组饲料喷洒全反式RA（MCE, HY-14649，300 mg/kg），对照组喷洒等体积无水乙醇。
2. 2.
   组织学分析：性腺样品经Bouin氏液固定、石蜡包埋、切片（5 μm）、苏木精-伊红（H&E）染色后光镜观察。
3. 3.
   RNA测序与数据分析：TRIzol法提取总RNA，构建Illumina文库后在NovaSeq X Plus平台进行150 bp双端测序；fastp过滤原始数据，HISAT2比对至黄鳍鲷参考基因组，RSEM计算FPKM值；DESeq2鉴定DEGs（|log₂FC|≥1且校正P<0.05）。
4. 4.
   生物信息学分析：STEM分析时间序列表达模式；KEGG富集分析（超几何检验，FDR校正q≤0.05）；STRING数据库（v11.5）构建PPI网络（置信度评分>700），Cytoscape可视化。
5. 5.
   RT-qPCR验证：HiScript III RT SuperMix反转录，LightCycler 480 SYBR Green I Master系统检测，β-肌动蛋白（β-actin）为内参，2^-ΔΔCt法计算相对表达量。

研究结果

3.1 RA处理增强了体细胞生长：在第90天时，RA处理组鱼的体重（23.40±6.08 g vs. 17.43±4.43 g, P<0.001）和体长（8.90±0.87 cm vs. 8.07±0.95 cm, P<0.01）显著高于对照组，但性腺指数无显著差异。

3.2 RA处理促进了生殖细胞分化：H&E染色显示，第90天时对照组精巢区出现初级精母细胞但卵巢区罕见初级卵母细胞，而RA处理组精巢区出现初级精母细胞的同时卵巢区可见少量初级卵母细胞（数量：8.17±3.43 vs. 0.50±0.84），表明RA优先推动了卵巢生殖细胞向初级卵母细胞阶段的进展。

3.3 转录组分析揭示了对RA处理的阶段特异性响应：PCA显示不同发育阶段及处理组间转录状态明显分离。RA处理组与对照组比较，第30天鉴定出1276个DEGs（210上调，1066下调），第60天为8172个DEGs（7230上调，942下调），第90天为989个DEGs（725上调，264下调）。维恩图显示仅25个DEGs为三时间点共有。KEGG富集分析表明，第30天主要涉及Rap1信号、Ras信号、MAPK信号、钙信号和视黄醇代谢等早期信号通路；第60天富集于谷氨酸能突触、GABA能突触、GnRH分泌、GnRH信号和雌激素信号等神经内分泌通路；第90天则富集于细胞黏附分子互作、ECM-受体互作、卵母细胞减数分裂和卵巢类固醇生成通路。

3.4 RA处理重塑了时间转录轨迹：时间序列比较显示RA处理组在第30-60天间的转录变化远大于对照组（7784上调 vs. 1894上调），且RA处理组在第30-60天即已显著富集GnRH分泌、卵巢类固醇生成、催产素信号和GnRH信号通路，而对照组这些通路至第30-90天才出现，证明外源性RA将减数分裂启动和类固醇生成的转录程序提前到了更早的发育窗口。

3.5 STEM分析揭示了RA相关的发育表达谱重编程：对照组主要富集Profile 4（神经活性配体-受体互作、钙信号）、Profile 7（DNA复制、细胞周期）和Profile 6（GABA能突触、谷氨酸能突触）。RA处理组上调的Profile 6显著富集于谷氨酸能突触、GABA能突触、GnRH分泌和催产素信号；上调的Profile 7富集于核质转运、剪接体和DNA复制；而下调的Profile 1富集于核糖体、糖酵解和视黄醇代谢。

3.6 热图和PPI分析将RA代谢与减数分裂和性别分化联系起来：层级聚类显示，减数分裂基因（sycp1、sycp2、sycp3、dmc1、rec8、dazl、piwil1/2）和性别分化/类固醇生成基因（foxl2、dmrt1、sox9a、cyp19a1a、cyp17a1、cyp11a1、esr1/2、pgr）在RA处理第60天和/或第90天表达增强。值得注意的是，RA降解酶cyp26a1和cyp26b1在第30天急剧诱导，而RA合成酶aldh1a2在对照组中高表达。PPI网络揭示了以Cyp19a1为中心的类固醇生成模块、包含Sycp1/2/3、Rec8、Dazl、Piwil1/2和Dmc1的减数分裂模块以及包含Cyp26a1、Cyp26b1、Aldh1a1、Aldh1a2、Rarab和Rxraa的RA代谢/受体模块之间的功能关联。

3.7 RT-qPCR验证了关键的RNA-seq表达模式：RT-qPCR证实了RA代谢基因（aldh1a2早期降低，cyp26a1和cyp26b1第30天强烈诱导）、核受体（rarab和rxraa呈阶段依赖性变化）、性别分化基因（sox9a、cyp19a1a、foxl2在第60天和/或第90天上调，dmrt1在第90天强增加）和减数分裂标志基因（dmc1、sycp3、dazl在第60天和/或第90天上调）的表达趋势与RNA-seq一致。

讨论与结论

讨论部分：RA处理对黄鳍鲷生殖细胞分化的促进作用需置于精巢卵巢正常发育序列的背景下理解。组织学上，最明显的效应并非两个谱系同步提前30天，而是优先推动卵巢生殖细胞向初级卵母细胞形成。这种模式可能反映了鲷科精巢卵巢的结构特点——卵巢区和精巢区共存但未必同步进展，RA可作用于两条谱系但在尚未达到初级卵母细胞阶段的谱系中产生更强的可见效应。类似的发育背景依赖性在其他雌雄同体硬骨鱼类中也有报道，如石斑鱼和稻田鳗中RA合成与降解基因在减数分裂卵母细胞出现前呈现阶段性动态变化。RA处理后的早期显著特征是cyp26a1和cyp26b1强烈诱导伴随aldh1a2表达降低，这与负反馈代谢响应一致——通过CYP26酶限制过量RA信号，同时抑制内源性RA合成。核受体rarab和rxraa的阶段依赖性变化提示受体水平也可能参与细胞对RA敏感性的微调。RA诱导的性腺响应呈现发育有序性而非同步激活：第60天foxl2和cyp19a1a的显著增加更符合卵巢体细胞/类固醇生成程序的激活，而非明确的早熟减数分裂启动。值得注意的是，RA处理同时影响了卵巢（foxl2、cyp19a1a）和精巢（dmrt1、sox9a）分化相关基因，这在精巢卵巢中具有生物学合理性——两条谱系共存，性别命运依赖于竞争性性腺程序之间的动态平衡。PPI网络进一步支持了这些过程之间的功能汇聚，但需注意该网络基于预测的功能关联而非实验验证的蛋白质互作，且全精巢卵巢RNA-seq无法区分细胞类型特异性起源。

研究结论：本研究表明RA处理促进了黄鳍鲷幼鱼的生殖细胞分化和减数分裂程序的更早分子激活。组织学证据显示最明显的效应出现在卵巢区域——RA处理组鱼在第90天出现少量初级卵母细胞而对照组未出现；两组鱼在第90天均存在精巢初级精母细胞。转录组学、热图、PPI和RT-qPCR分析表明，RA处理激活了减数分裂标志物和性别分化基因，同时诱导了涉及cyp26a1、cyp26b1和aldh1a2的反馈响应。这些发现提示RA生物活性受内部代谢检查点调控，并与雄性先熟精巢卵巢的区室化发育序列相互作用。研究结果为理解雌雄同体硬骨鱼类生殖细胞分化提供了分子基础。