【研究意义】动物的性别决定机制有两种类型，分别是遗传性别决定型（Genetic sex determination，GSD）和环境依赖性别决定型（Environment dependent sex determination，ESD）。高等脊椎动物，如哺乳动物和鸟类的性别决定机制属于遗传性别决定型^[1]，而鱼类、两栖类和爬行类动物的性别决定机制不完全依赖于染色体，环境温度、pH值、激素、行为模式等因素都会改变和影响个体的性别分化^[2-5]，而且这种影响可以贯穿个体发育的全过程。因此，在很多鱼类和爬行类动物中都可以观察到环境因素诱导的性别逆转现象。

在水产养殖中，鱼类的个体大小、外形、体色等具有经济价值的性状通常与某一特定性别关联。如雄性罗非鱼（Oreochromis niloticus）和乌斑杂交鳢（Channa）都比其雌性个体生长速度快、体型大，具有更高的经济价值^[6-7]；而大菱鲆（Scophthalmus maximus）雌鱼体型比雄鱼大50%^[8]。西部食蚊鱼（Gambusia affinis）为暖水性的小型硬骨鱼类，生活于水库、湖泊、坝塘、沼泽、稻田、水渠、洼地等各类静水水体中。因其嗜食蚊子幼虫及具有观赏价值而被引入我国，分布于长江以南诸多省区^[9]。由于食蚊鱼具有容易获得、生长速度快、繁殖能力强、性别二态性，且与哺乳动物对化学物质的反应方式相似等特点，近年来食蚊鱼被当成一种很好的模式生物，广泛应用于进化和生态等研究领域^[10-12]。因此，利用食蚊鱼研究鱼类性别逆转，有助于在水产养殖业中控制种群性别比例、培养全单性种群，具有重要的科学意义和经济价值^[13-15]。

【前人研究进展】食蚊鱼的性别决定染色体为ZW/ZZ型，但其性别分化同时受到环境因素的影响。研究发现，养猪场、纸浆造纸厂和城镇污水处理厂附近水域中高含量的孕激素使雌性食蚊鱼的性腺发育、生殖交配行为受到损害，并使其出现雄鱼的第二性征，进而改变该种群的雌雄比例^[16-17]。深入研究表明由雌性食蚊鱼通过性别逆转产生的伪雄性食蚊鱼具有与野生雄性食蚊鱼相似的追逐、摇摆、交配等行为^[18-19]，但其摆动臀鳍的频率更低。此外，这些伪雄性食蚊鱼的生殖足无法与雌性个体完全匹配，需要在交配中保持更高的交配频率。以上研究表明，雌性食蚊鱼在环境中高剂量孕激素的诱导下可发生性别逆转现象，且其产生的伪雄性个体获得了部分生殖行为能力。

【本研究切入点】虽然目前针对硬骨鱼类的性别逆转已有一些研究，但针对食蚊鱼的性别逆转研究鲜有报道，从单细胞基因表达水平探讨其性别逆转机制的研究仍为空白。本研究以养殖业中广泛应用的新型孕激素——孕二烯酮（Gestodene，GES）作为外源激素^[20]，以西部食蚊鱼(以下简称食蚊鱼）为研究对象，从形态学和单细胞转录组水平研究孕二烯酮诱导下雌性食蚊鱼雄性化过程中性腺单细胞水平基因表达的变化。【拟解决的关键问题】本研究利用单细胞转录组测序技术，从单细胞水平分析孕二烯酮诱导后雌性食蚊鱼性腺基因表达谱的变化，挖掘雌性食蚊鱼雄性化的关键基因，初步解析食蚊鱼性别逆转的分子机理，丰富动物性腺单细胞数据，为后续研究提供更多有效的数据资源。

1　材料与方法 1.1　试验材料

2022年10月于广州市花地湾花鸟鱼虫市场购买4月龄左右的西部食蚊鱼，转至广州大学实验室水循环养殖系统培养驯化，设置温度26(±1)℃，光周期14 h: 10h（光照︰黑暗），水体硬度140(±2.5) mg/L，水体导电性20(±0.3) μS/cm，pH值7.4(±0.5)，溶解氧含量7.4(±0.2) mg/L。经过2个月驯化稳定后，食蚊鱼死亡率小于5%，已性成熟可用于试验。实验鱼的使用和实验操作符合“实验动物3R原则” ^[21]。

主要试剂包括孕二烯酮（GES, 纯度99.4%）购自上海西格玛奥德里奇贸易有限公司（Sigma-Aldrich）、二甲基亚砜（DMSO, 纯度99.9%）购自上海生工生物工程公司。500 ng/L孕二烯酮溶液：将24 μL 2 mg/mL孕二烯酮溶液溶于4 776 μL二甲基亚砜中，取1 250 μL加进25 L水中。

4 ℃预冷的400 μg/mL磷酸盐缓冲生理盐水（PBS）溶液（含0.04% BSA, 无Ca²⁺、Mg²⁺）；30℃预热的400 μg/mL 1×HBSS缓冲液（含0.04% BSA，有Ca²⁺、Mg²⁺）；消化酶混合液：1 000 μL 1× HBSS缓冲液（30℃预热）+ 1 μL DNA酶Ⅰ（2 000 U/mL）+ 10 μL胶原酶Ⅰ（10%）+ 10 μL胶原酶Ⅱ（10%）+ 10 μL胶原酶Ⅳ（10%）+ 5 μL透明质酸酶（10%）+ 10 μL中性蛋白酶（5%）+ 2 μL离子缓冲液。

1.2　试验方法

1.2.1 试验分组与处理　挑选体长为3.0（±0.5）cm、体质量为130（±50）mg的成年雌性食蚊鱼（6月龄）和体长为1.9（±0.6）cm、体质量为110（±70）mg的成年雄性食蚊鱼（6月龄）参与孕二烯酮暴露实验。设置雄性和雌性食蚊鱼的空白对照（FC，female control；MC，male control）以及雌性食蚊鱼孕二烯酮暴露4周和8周的处理组（F4，female treated for 4 weeks；F8，female treated for 8 weeks），各组设3个重复。孕二烯酮处理组为每个容积12 L的缸中放置10条雌性成鱼和10 L孕二烯酮暴露溶液。

1.2.2 单细胞悬液制备　用解剖剪刀由食蚊鱼肛门开始沿腹部剪至鳃盖下缘，再从两端分别向背部解剖，在体视镜下取出性腺组织，使用预冷至4 ℃的1×PBS洗涤2~3次。将组织切成小碎片，然后转移到消化酶混合液中，30 ℃消化30~50 min，每隔3 min摇1次，进行初步显微镜检查。消化后的细胞悬液用70 μm过滤器过滤，用2 mL 1× PBS清洗离心管和细胞筛。用水平离心机在4 ℃下以200 g离心5 min；取上清，加入5 mL 1×PBS重悬细胞，再以200 g离心5 min，洗涤2~3次；根据沉淀的量加入预冷的1×PBS溶液（沉淀较大加入1 mL，较少则加入0.5 mL，几乎看不到沉淀则加100 μL），轻轻混匀。利用显微镜细胞计数板和0.2% 台盼蓝染色液检测细胞浓度和活性并记录。

1.2.3 单细胞文库构建及测序　使用10×Genomics^TM微流控芯片，将单细胞悬液和Chromium Single Cell 3’ Reagent（v3，10×Genomics）试剂盒混合，在油相反应体系中形成单细胞珠状凝胶液滴^[22]。随后，在53 ℃下反转录45 min、85 ℃下反应5 min后终止反应，制备出测序所需的cDNA。将cDNA进行PCR扩增，对得到的扩增cDNA进行片段化、末端修复和补全polyA尾巴，添加测序接头并进行文库扩增，最终制备出单细胞cDNA文库。单细胞测序采用10×Genomics Novaseq 6000测序平台，并用CellRanger（10×Genomics^TM）比对得到原始测序数据。

1.3　数据处理及分析

使用Seurat^[23]进行主成分分析（Principle Component Analysis，PCA），选择适当的主成分进行下一步非监督聚类和UMAP（Uniform Manifold Approximation and Projection）分析。根据单细胞聚类结果，在Seurat中进行特征基因筛选，确定每个细胞亚群的特征基因，并利用CellKB^[24]和Fish Enrichr^[25-26]数据库鉴定每个细胞亚群对应的细胞类型。

根据细胞类型注释的结果，选择相关的细胞群进行细胞轨迹分析。使用Monocle^[27]将单细胞数据集转换成CDS（Cell dataset object）类型，进行预处理（使用降维函数将数据缩小到较低维度，reduction_method =“UMAP”，preprocess_method=“Aligned”），嵌入降维数据，进行拟时序计算。在过滤掉低质量数据后进行细胞聚类分析，根据不同的细胞状态分类并进行可视化分析。

根据上述拟时序计算和细胞分化轨迹结果，找到关键细胞群的分化节点，通过BEAM（Branch Expression Analysis Modeling）分析筛选对分化节点最关键的基因^[28-29]，并使用cluster Profiler^[30]进行KEGG富集分析。

2　结果与分析 2.1　孕二烯酮诱导后的食蚊鱼鳍部形态变化

孕二烯酮诱导后的雌性食蚊鱼臀鳍形态发生了雄性化的变化。野生雌性食蚊鱼的臀鳍为扇形或三角形，而雄性的臀鳍具有特化的生殖足结构。经孕二烯酮诱导处理4~8周后，雌性食蚊鱼臀鳍发生了不同程度的特化现象，呈现出类似于雄性生殖足的结构（图 1）。结果表明，孕二烯酮在雌性食蚊鱼的性别逆转过程中起重要的调节作用，能够促进其雄性化外部特征的出现。

2.2　孕二烯酮诱导后的食蚊鱼性腺细胞构成变化

为进一步探究孕二烯酮诱导下雌性食蚊鱼雄性化的分子机制，本研究分别收集MC、FC、F4和F8组的食蚊鱼性腺进行单细胞测序。测序产生的数据使用CellRanger软件进行单细胞基因表达定量分析；将表达数据进行质量控制和筛选，剔除低质量单细胞数据，最终获得每个样本的单细胞基因表达谱。MC、FC、F4和F8样本的测序数据如表 1所示，单个样本有效细胞数介于8 982~13 873之间，检测到的有效表达基因数目介于19 891~20 250之间，有效Barcodes比例均在97% 以上，表明数据质量优秀。

分析使用Sctransform对4个样本的数据进行归一化并整合，利用UMAP进行降维处理。由图 2B可见，野生雌性FC组和孕二烯酮处理后的雌性F4、F8组细胞重叠较好，而雄性MC组的细胞相对呈独立分布特征，说明成年雌性食蚊鱼和雄性食蚊鱼性腺细胞构成具有较大差异，且经孕二烯酮诱导后的雌性性腺细胞构成仍与雄性不同。

为了确定食蚊鱼性腺的具体细胞构成类型和对应的基因表达情况，本研究进一步对细胞进行聚类，提取每个细胞类群特异高表达的marker基因，利用CellKB和Fish Enrichr对细胞类型进行注释分析。剔除血液细胞和免疫细胞后，结果发现食蚊鱼性腺一共包含6种主要的细胞类型（图 2A，表 2）。不同样本间比较发现，雄性食蚊鱼性腺主要由成纤维细胞、雄性生殖细胞、上皮细胞和原始生殖细胞构成，而雌性性腺主要由内皮细胞、成纤维细胞、卵母细胞和原始生殖细胞构成。孕二烯酮诱导处理4~8周后，雌性食蚊鱼性腺中产生少量的雄性生殖细胞（F4 0.79%，F8 0.39%）。此外，诱导8周后性腺中原始生殖细胞比例较野生雌性样本极显著上升。

2.3　原始生殖细胞的基因表达与分化

2.3.1 原始生殖细胞的基因表达　提取来自FC、F4和F8样本的原始生殖细胞进行重聚类和分析，并重点分析与性别决定相关重要基因的表达情况。结果表明，来自FC、F4和F8的原始生殖细胞可以被划分为9个亚类型（图 3）。值得注意的是，不同性别分化的关键基因在不同细胞类群里出现差异表达情况。几个与雄性生殖细胞分化密切相关的基因（如amh、fabp3、dmrt1、wt1a和wt1b）在cluster1和cluster3中高表达；对雌性性腺发育和卵子形成有重要作用的cyp19a1a及fgf16基因则在cluster 0和cluster 2中高表达（图 4）。结果表明，在孕二烯酮诱导下原始生殖细胞类群中同时出现了两种具有不同性别决定潜能的细胞。

2.3.2 原始生殖细胞的分化　利用单细胞轨迹分析和拟时序分析研究原始生殖细胞的分化过程以及对应的基因表达变化序列，结果（图 5B）显示，来自FC的样本主要集中分布在原始生殖细胞分化轨迹的早期，来自F4和F8的细胞主要集中分布在分化轨迹的中晚期。孕二烯酮诱导后，原始生殖细胞出现了两种截然不同的分化方向。结合图 4中性别分化关键基因的表达情况，推测部分原始生殖细胞出现向雄性生殖细胞分化的趋势，而部分细胞维持向雌性生殖细胞分化的趋势（图 5C）。

基于上述细胞分化轨迹，通过BEAM分析发现，有2 477个基因在两个不同的分化方向表现出不同的表达模式（图 6A）。对差异最大的前200个基因进行KEGG通路富集分析，结果发现Toll-like receptor、TGF-beta、Notch等信号通路参与了原始生殖细胞向雄性化逆转的过程（图 6B）；甘氨酸，丝氨酸、苏氨酸、半胱氨酸、甲硫氨酸代谢、氨基糖和核苷酸糖代谢等途径也与该过程相关（图 6C）；不饱和脂肪酸、类固醇激素生物合成，醚酯、亚油酸代谢等途径则对维持原始生殖细胞向雌性性腺分化有重要作用（图 6B）。

3　讨论

原始生殖细胞是一类生殖干细胞，具有向不同性别生殖细胞分化的潜能，是鱼类性别决定和分化过程中最重要的一类细胞。其形成方式目前存在两种学说，第一种是先成论，指受精卵带有母源生殖质，随着细胞分裂进行，生殖质被不均等地分配，获得生殖质的细胞将发育成原始生殖细胞，否则发育成体细胞，如斑马鱼（Danio rerio）、泥鳅（Misgurnus anguillicaudatus）和大菱鲆等；另一种是后成论，指受精卵早期分裂中的体细胞与生殖细胞没有区别，但会因某些细胞间信号诱导分化成原始生殖细胞，如哺乳动物、青鳉(Oryzias latipes)、真鲷（Pagrosomus major）^[31-34]。本研究发现孕二烯酮刺激下，随着时间增加，雌性食蚊鱼中的原始生殖细胞数量不断上升，推测食蚊鱼原始生殖细胞形成方式属于后成论。通过富集分析发现Toll-like receptor、TGF-beta、Notch等信号通路参加了食蚊鱼雄性化的过程，推测以上通路也可能诱导原始生殖细胞生成。

原始生殖细胞的基因表达和细胞数量都有可能影响鱼类的性别决定。在斑马鱼和青鳉的研究中发现，原始生殖细胞数量不足的个体倾向于发育为雄鱼，而原始生殖细胞较多的个体则倾向于发育为雌鱼^[35-37]。相比之下，泥鳅的性别决定则与原始生殖细胞的数量无关^[38]。本研究发现孕二烯酮刺激下雌性食蚊鱼出现向雄性化逆转的现象，且伴随原始生殖细胞数量显著上升，这一结果与在斑马鱼、青鳉和泥鳅中观察到的情况均不相同。分析其原因可能与不同鱼类多样化的性别决定机制有关，且这一结果进一步说明原始生殖细胞数目是决定鱼类性别分化的重要因素之一。

本研究发现孕二烯酮刺激下有几个与鱼类性别分化密切相关的基因（amh、fabp3、dmrt1、wt1a、wt1b、cyp19a1a和fgf16）在原始生殖细胞内出现差异表达的情况。其中TGFβ超家族成员amh基因是硬骨鱼睾丸分化的关键基因。哺乳动物中该基因诱导缪勒氏管退化，并阻止雌性生殖管发育^[39]。而在鱼类中amh被发现可以通过调节cyp11a1、cyp17和cyp19a1a的表达进而影响性腺中的睾酮含量^[40]。研究发现，东部食蚊鱼（Gambusia holbrooki）雄性胚胎中amh基因的表达量是雌性胚胎的25倍^[41]。本研究中amh基因在孕二烯酮诱导后在部分原始生殖细胞中出现了高表达的情况，说明这部分细胞可能出现了向雄性化分化的特征。此外，本研究还发现amh和cyp19a1a基因具有截然相反的表达模式，侧面印证了性腺分化过程中amh基因对cyp19a1a的抑制作用。

促卵泡激素（Follicle-stimulating hormone，FSH）和促黄体化激素（Luteninzing hormone，LH）通过与其对应的受体FSHR和LHCGR相结合调控性腺发育。fshr和lhcgr基因在青鳉卵巢和睾丸中均具有较高的表达水平，在睾丸形成过程中，两个基因的表达均有上调现象^[42]，说明这两个基因对雄性性腺的发育有重要作用。本研究发现，在食蚊鱼原始生殖细胞向雄性化分化的过程中，fshr出现上调现象，但lhcgr则并未出现明显上调。由于在卵泡形成过程中lhcgr的表达相比fshr存在一定的滞后现象^[43]，推测这两个基因在雄性化过程中的表达同样存在表达时间差异的可能。比较野生雌性食蚊鱼和孕二烯酮处理后雌性食蚊鱼的单细胞图谱，本研究发现诱导后的原始生殖细胞开始出现向雄性化分化的特征，出现极少量的分化后雄性生殖细胞（精原细胞），这说明诱导后的雌鱼已经在一定程度上具有雄性食蚊鱼的生殖功能，但其生殖能力还远远低于正常的雄性食蚊鱼。在本试验所用孕二烯酮剂量下，完全的雄性化过程可能需要更长的时间。在前期的预实验中发现，当环境中孕二烯酮浓度达到500 ng/L时，培养2个月后的食蚊鱼就开始出现死亡情况，且随着时间增加死亡数目不断上升，3个月时70% 的食蚊鱼已经死亡。因此，为了确保获得足够的样本进行测序和分析，本研究选择8周作为最长的暴露时间来观察食蚊鱼的性别逆转现象。下一步研究将通过调节暴露浓度和延长暴露时间进一步跟踪雌性食蚊鱼完全雄性化的过程。

已有多项研究证明雌性食蚊鱼暴露在孕激素中会诱导出与雄鱼类似的生殖足，即发生臀鳍雄性化。雄鱼生殖足的形成和发育主要由类固醇激素决定，特别是雄激素^[44-45]。在一些鱼类中，雄激素尤其是11- 酮睾酮（11KT）可刺激鱼鳍的伸长，而雄激素拮抗剂或雄激素合成抑制剂会干扰睾丸功能和外生殖器发育^[46]。研究发现ar、shh、ptc1、cyp19a、er和vtgc等基因都与性别二态性鱼类臀鳍发育相关^[47-48]，本研究也发现cyp19a基因存在特异性表达，与以往研究结果符合。

单细胞测序是一项新兴的组学研究技术，凭借其高通量的准确性和特异性，广泛应用于生物的各种生理过程研究。目前利用单细胞测序分析来研究鱼类性别逆转的报道有半滑舌鳎（Cynoglossus semilaevis），其雌性个体在发育早期用高温诱导可使雌鱼的性逆转率达到73%^[49]。利用单细胞测序技术发现，在未分化的伪雄性精原细胞中，CaSR基因和MAPK信号因子上调^[50]。本研究未发现该基因和信号因子上调，可能原因是不同机制诱导的性别逆转涉及到的信号通路不同，以及不同鱼类遗传背景也会导致性别转变机制不同。鱼类的性别逆转是个复杂过程，仍需要大量的研究来填补该空白领域。

4　结论

将雌性食蚊鱼置于浓度为500 ng/L孕二烯酮溶液暴露4~8周后，发现雌性食蚊鱼臀鳍出现雄性化特征；分别取野生雌性、雄性食蚊鱼和暴露4、8周雌性食蚊鱼的性腺进行单细胞转录组测序并进行数据分析，结果发现野生雌性食蚊鱼和雄性食蚊鱼性腺细胞构成明显不同，暴露4、8周后的雌性食蚊鱼性腺中出现少量的雄性生殖细胞（精原细胞），且原始生殖细胞数量有明显上升；在孕二烯酮诱导下雌性食蚊鱼原始生殖细胞中性别决定关键基因amh、fabp3、dmrt1、wt1a、wt1b、cyp19a1a和fgf16出现差异化表达情况；细胞轨迹分析表明，部分原始生殖细胞出现向雄性化分化的特征。Toll-like receptor、TGF-beta、Notch和一些氨基酸代谢信号通路参与了原始生殖细胞向雄性化逆转的过程。