基于单细胞水平的孕二烯酮诱导雌性西部食蚊鱼雄性化的研究.pdf

【目的】从单细胞转录组水平研究雌性西部食蚊鱼（Gambusia affi nis）性别逆转过程中涉及到的性腺细胞分化发育和基因表达变化，探讨孕二烯酮（Gestodene，GES）作用下食蚊鱼性别反转的分子机制。【方法】将培养驯化 2 个月后的性成熟雌性食蚊鱼暴露于浓度为 500 ng/L 的孕二烯酮溶液中进行培养，于第 4、8 周取样观察其形态变化，并分别取正常雌性对照、雄性对照、暴露 4 周和 8 周样本的性腺制成单细胞悬液，进行单细胞转录组测序并分析。【结果】经孕二烯酮处理后的雌性食蚊鱼形态上出现臀鳍类生殖足现象。单细胞转录组测序结果发现，雄性食蚊鱼和雌性食蚊鱼的性腺具有完全不同的细胞构成。雄性食蚊鱼性腺主要由成纤维细胞、雄性生殖细胞、上皮细胞和原始生殖细胞构成，而雌性性腺主要由内皮细胞、成纤维细胞、卵母细胞和原始生殖细胞构成。暴露 4、8 周后雌性食蚊鱼性腺中出现少量的雄性生殖细胞，且原始生殖细胞数量明显上升。进一步分析发现，原始生殖细胞中性别决定关键基因出现差异化表达情况。细胞轨迹分析表明，部分原始生殖细胞出现向雄性分化的特征，其中有 2 477 个基因在两个不同的分化方向表现出不同的表达模式。KEGG 通路富集分析结果表明，Toll-like receptor、TGF-beta、Notch 等信号通路参与原始生殖细胞向雄性化逆转的过程，甘氨酸，丝氨酸、苏氨酸、半胱氨酸、甲硫氨酸代谢、氨基糖和核苷酸糖代谢等途径也与该过程相关。【结论】雌性食蚊鱼经孕二烯酮诱导后，在形态学上出现臀鳍雄性化特征，性腺中有部分原始生殖细胞向雄性生殖细胞转化。原始生殖细胞中性别决定关键基因 amh、fabp3、dmrt1、wt1a、wt1b、cyp19a1a 和 fgf16 出现差异化表达情况。Toll-like receptor、TGF-beta、Notch 和一些氨基酸代谢信号通路参与了原始生殖细胞向雄性化逆转的过程。关键词：孕二烯酮（GES）；西部食蚊鱼；硬骨鱼；性别逆转；单细胞转录组中图分类号：S917.4 文献标志码：A 文章编号：1004-874X（2023）06-0126-11Exploration of GES-induced masculinization of Gambusia affi nis Based on Single Cell LevelFENG Tingyin,DENG Xikai,XING Ke（College of Life Sciences,Guangzhou University,Guangzhou 510006,China）Abstract:【Objective】To investigate the molecular mechanism of sex reversal in female Western mosquitofish(Gambusia affinis)under the effect of gestodene(GES)by studying the development of gonadal cell differentiation and gene expression changes involved in the process of sex reversal in female mosquitofish at the single-cell transcriptome level.【Method】Sexually mature female mosquitofish were exposed to 500 ng/L GES after culture and domestication for 2 months,and samples were taken at the 4th and 8th week to observe morphological changes,and single-cell suspensions were made from the gonads of normal female controls,male controls,and samples exposed for 4 and 8 weeks,respectively,127【研究意义】动物的性别决定机制有两种类型，分别是遗传性别决定型（Genetic sex determination，GSD）和 环 境 依 赖 性 别 决 定 型（Environment dependent sex determination，ESD）。高等脊椎动物，如哺乳动物和鸟类的性别决定机制属于遗传性别决定型1，而鱼类、两栖类和爬行类动物的性别决定机制不完全依赖于染色体，环境温度、pH 值、激素、行为模式等因素都会改变和影响个体的性别分化2-5，而且这种影响可以贯穿个体发育的全过程。因此，在很多鱼类和爬行类动物中都可以观察到环境因素诱导的性别逆转现象。在水产养殖中，鱼类的个体大小、外形、体色等具有经济价值的性状通常与某一特定性别关联。如雄性罗非鱼（Oreochromis niloticus）和乌斑杂交鳢（Channa）都比其雌性个体生长速度快、体型大，具有更高的经济价值6-7；而大菱鲆（Scophthalmus maximus）雌鱼体型比雄鱼大50%8。西部食蚊鱼（Gambusia affi nis）为暖水性的小型硬骨鱼类，生活于水库、湖泊、坝塘、沼泽、稻田、水渠、洼地等各类静水水体中。因其嗜食蚊子幼虫及具有观赏价值而被引入我国，分布于长江以南诸多省区9。由于食蚊鱼具有容易获得、生长速度快、繁殖能力强、性别二态性，且与哺乳动物对化学物质的反应方式相似等特点，近年来食蚊鱼被当成一种很好的模式生物，广泛应用于进化和生态等研究领域10-12。因此，利用食蚊鱼研究鱼类性别逆转，有助于在水产养殖业中控制种群性别比例、培养全单性种群，具有重要的科学意义和经济价值13-15。【前人研究进展】食蚊鱼的性别决定染色体为 ZW/ZZ 型，但其性别分化同时受到环境因素的影响。研究发现，养猪场、纸浆造纸厂和城镇污水处理厂附近水域中高含量的孕激素使雌性食蚊鱼的性腺发育、生殖交配行为受到损害，并使其出现雄鱼的第二性征，进而改变该种群的雌雄比例16-17。深入研究表明由雌性食蚊鱼通过性别逆转产生的伪雄性食蚊鱼具有与野生雄性食蚊鱼相似的追逐、摇摆、交配等行为18-19，但其摆动臀鳍的频率更低。此外，这些伪雄性食蚊鱼的生殖足无法与雌性个体完全匹配，需要在交配中保持更高的交配频率。以上研究表明，雌性食蚊鱼在环境中高剂量孕激素的诱导下可发生性别逆转现象，且其产生的伪雄性个体获得了部分生殖行为能力。【本研究切入点】虽然目前针对硬骨鱼类的性别逆转已有一些研究，但针对食蚊鱼的性别逆转研究鲜有报道，从单细胞基因表达水平探讨其性别逆转机制的研究仍为空白。本研究以养殖业中广泛应用的新型孕激素孕二烯酮（Gestodene，GES）作为外源激素20，以西部食蚊鱼(以下简称食蚊鱼）为研究对象，从形态for single-cell transcriptome sequencing and analysis.【Result】The morphology of GES-treated female mosquitofish showed gonopodium,and single-cell transcriptome sequencing revealed completely different cellular composition of gonads between male and female mosquitofish.The gonads of male mosquitofish consist mainly of fibroblasts,male germ cells,epithelial cells and primordial germ cells,whereas the gonads of female mosquitofish consist mainly of endothelial cells,fibroblasts,oocytes and primordial germ cells.After 4 and 8 weeks of exposure,a small amount of male germ cells appeared in the gonads of female mosquitofish,and the number of primordial germ cells increased significantly.Further analysis revealed that several crucial genes related to the sex-differentiation were differentially expressed in primordial germ cells.Cell trajectory analysis showed that some primordial germ cells showed characteristics of differentiation towards males,in which 2 477 genes showed different expression patterns in two different different directions of differentiation.The results of KEGG pathway enrichment analysis revealed that signaling pathways,such as Toll-like receptor,TGF-beta,Notch,etc.,were involved in the process of reversal of primordial germ cells towards maleness.Glycine,serine,threonine,cysteine,methionine metabolism,amino sugar and nucleotide sugar metabolism and other pathways were also related to this process.【Conclusion】After exposure to GES solution,female mosquitofish showed morphologically masculinized anal fins,and some primordial germ cells in the gonads were transformed to male germ cells.Differential expression of amh,fabp3,dmrt1,wt1a,wt1b,cyp19a1a and fgf16,which are key genes for sex determination,was observed in the primordial germ cells.Toll-like receptor,TGF-beta,Notch,and some amino acid metabolism signaling pathways were involved in the process of sex reversal.Key words:gestodene(GES);Gambusia affi nis;teleost;sex reversal;single-cell transcriptome 128学和单细胞转录组水平研究孕二烯酮诱导下雌性食蚊鱼雄性化过程中性腺单细胞水平基因表达的变化。【拟解决的关键问题】本研究利用单细胞转录组测序技术，从单细胞水平分析孕二烯酮诱导后雌性食蚊鱼性腺基因表达谱的变化，挖掘雌性食蚊鱼雄性化的关键基因，初步解析食蚊鱼性别逆转的分子机理，丰富动物性腺单细胞数据，为后续研究提供更多有效的数据资源。1 材料与方法1.1 试验材料2022 年 10 月于广州市花地湾花鸟鱼虫市场购买 4 月龄左右的西部食蚊鱼，转至广州大学实验室水循环 养殖系统培养驯化，设置温度26(1)，光周期 14 h10h（光照黑暗），水体硬度 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,无 Ca2+、Mg2+）；30预热的 400 g/mL 1HBSS 缓冲液（含 0.04%BSA，有 Ca2+、Mg2+）；消化酶混合液：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 的 1PBS 洗涤 23 次。将组织切成小碎片，然后转移到消化酶混合液中，30 消化 3050 min，每隔 3 min 摇 1 次，进行初步显微镜检查。消化后的细胞悬液用 70 m 过滤器过滤，用 2 mL 1 PBS 清洗离心管和细胞筛。用水平离心机在4 下以 200 g 离心 5 min；取上清，加入 5 mL 1PBS 重悬细胞，再以 200 g 离心 5 min，洗涤23次；根据沉淀的量加入预冷的1PBS溶液（沉淀较大加入 1 mL，较少则加入 0.5 mL，几乎看不到沉淀则加 100 L），轻轻混匀。利用显微镜细胞计数板和 0.2%台盼蓝染色液检测细胞浓度和活性并记录。1.2.3 单细胞文库构建及测序 使用10GenomicsTM微流控芯片，将单细胞悬液和 Chromium Single Cell 3 Reagent（v3，10Genomics）试剂盒混合，在油相反应体系中形成单细胞珠状凝胶液滴22。随后，在 53 下反转录 45 min、85 下反应 5 min 后终止反应，制备出测序所需的 cDNA。将cDNA 进行 PCR 扩增，对得到的扩增 cDNA 进行片段化、末端修复和补全 polyA 尾巴，添加测序接头并进行文库扩增，最终制备出单细胞 cDNA文 库。单 细 胞 测 序 采 用 10Genomics Novaseq 6000测序平台，并用CellRanger（10GenomicsTM）比对得到原始测序数据。1.3 数据处理及分析使用 Seurat23进行主成分分析（Principle Component Analysis，PCA），选 择 适 当 的 主 成分 进 行 下 一 步 非 监 督 聚 类 和 UMAP（Uniform Manifold Approximation and Projection）分析。根据单细胞聚类结果，在 Seurat 中进行特征基因筛选，确定每个细胞亚群的特征基因，并利用 CellKB 24和 Fish Enrichr25-26数据库鉴定每个细胞亚群对应的细胞类型。根据细胞类型注释的结果，选择相关的细胞群进行细胞轨迹分析。使用 Monocle27将单细胞数据集转换成 CDS（Cell dataset object）类型，129进行预处理（使用降维函数将数据缩小到较低维度，reduction_method=“UMAP”，preprocess_method=“Aligned”），嵌入降维数据，进行拟时序计算。在过滤掉低质量数据后进行细胞聚类分析，根据不同的细胞状态分类并进行可视化分析。根据上述拟时序计算和细胞分化轨迹结果，找到关键细胞群的分化节点，通过 BEAM（Branch Expression Analysis Modeling）分析筛选对分化节点最关键的基因28-29，并使用 cluster Profiler30进行 KEGG 富集分析。2 结果与分析2.1 孕二烯酮诱导后的食蚊鱼鳍部形态变化孕二烯酮诱导后的雌性食蚊鱼臀鳍形态发生了雄性化的变化。野生雌性食蚊鱼的臀鳍为扇形或三角形，而雄性的臀鳍具有特化的生殖足结构。经孕二烯酮诱导处理 48 周后，雌性食蚊鱼臀鳍发生了不同程度的特化现象，呈现出类似于雄性生殖足的结构（图 1）。结果表明，孕二烯酮在雌性食蚊鱼的性别逆转过程中起重要的调节作用，能够促进其雄性化外部特征的出现。图 1 孕二烯酮诱导后的食蚊鱼鳍部形态变化Fig.1 Changes in fin morphology of mosquitofish exposed to GES2.2 孕二烯酮诱导后的食蚊鱼性腺细胞构成变化为进一步探究孕二烯酮诱导下雌性食蚊鱼雄性化的分子机制，本研究分别收集 MC、FC、F4和 F8 组的食蚊鱼性腺进行单细胞测序。测序产生的数据使用 CellRanger 软件进行单细胞基因表达定量分析；将表达数据进行质量控制和筛选，剔除低质量单细胞数据，最终获得每个样本的单细胞基因表达谱。MC、FC、F4 和 F8 样本的测序数据如表 1 所示，单个样本有效细胞数介于8 98213 873 之间，检测到的有效表达基因数目介于 19 89120 250 之间，有效 Barcodes 比例均在 97%以上，表明数据质量优秀。分析使用 Sctransform 对 4 个样本的数据进行归一化并整合，利用 UMAP 进行降维处理。由图2B 可见，野生雌性 FC 组和孕二烯酮处理后的雌性 F4、F8 组细胞重叠较好，而雄性 MC 组的细胞相对呈独立分布特征，说明成年雌性食蚊鱼和雄性食蚊鱼性腺细胞构成具有较大差异，且经孕二烯酮诱导后的雌性性腺细胞构成仍与雄性不同。为了确定食蚊鱼性腺的具体细胞构成类型和对应的基因表达情况，本研究进一步对细胞进行聚类，提取每个细胞类群特异高表达的 marker 基因，利用 CellKB 和 Fish Enrichr 对细胞类型进行注释分析。剔除血液细胞和免疫细胞后，结果发现食蚊鱼性腺一共包含 6 种主要的细胞类型（图2A，表 2）。不同样本间比较发现，雄性食蚊鱼性腺主要由成纤维细胞、雄性生殖细胞、上皮细胞和原始生殖细胞构成，而雌性性腺主要由内皮细胞、成纤维细胞、卵母细胞和原始生殖细胞构成。孕二烯酮诱导处理 48 周后，雌性食蚊鱼性腺中产生少量的雄性生殖细胞（F4 0.79%，F8 0.39%）。此外，诱导 8 周后性腺中原始生殖细胞比例较野生雌性样本极显著上升。2.3 原始生殖细胞的基因表达与分化2.3.1 原始生殖细胞的基因表达 提取来自 FC、F4和F8样本的原始生殖细胞进行重聚类和分析，表 1 食蚊鱼性腺单细胞测序数据概况Table 1 Overview of single-cell sequencing data on gonads of mosquitofish组别Group有效Barcodes Valid barcodes(%)测序质量高于 30 的碱基比例Q30 bases in RNA Reads(%)总细胞数Number of cells总基因数Total number of genesMC97.690.11387320250FC97.889.3898220024F498.289.71073619960F898.390.61243919891130并重点分析与性别决定相关重要基因的表达情况。结果表明，来自 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 中性别分化关键基因的表达情况，推测A：细胞类型 UMAP；B：细胞来源 UMAPA:UMAP of cell types;B:UMAP of cell sources图 2 食蚊鱼性腺单细胞 UMAP 图谱Fig.2 Single-cell UMAP diagrams of gonads in mosquitofish表 2 食蚊鱼性腺单细胞细胞类型比例Table 2 Proportion of cell types in gonads of mosquitofish(%)组别Group成纤维细胞Fibroblast内皮细胞 Endothelial cell上皮细胞 Epithelial cell雄性生殖细胞Male germ cell卵母细胞Oocyte原始生殖细胞Primordial germ cellMC55.784.7014.2622.770.002.49FC14.8352.090.050.0024.978.06F427.2137.100.020.7924.1610.72F834.4721.000.120.3918.3425.68A：细胞状态 UMAP；B：细胞来源 UMAPA:UMAP of cell state;B:UMAP of cell sources图 3 原始生殖细胞状态和细胞来源 UMAP 图谱Fig.3 UMAP diagrams of primordial germ cells state and sources131图 4 原始生殖细胞性别特征基因表达图谱Fig.4 Gene expression profiles of sex characteristics in primordial germ cellsA：不同细胞亚类型分化轨迹；B：不同细胞来源分化轨迹；C：细胞拟时序分化轨迹A:Differentiation trajectories of different seurat-clusters of cell;B:Differentiation trajectories of different sample sources;C:Differentiation trajectories of cell pseudotime 图 5 原始生殖细胞分化轨迹Fig.5 Trajectory analysis of primordial germ cells部分原始生殖细胞出现向雄性生殖细胞分化的趋势，而部分细胞维持向雌性生殖细胞分化的趋势（图 5C）。基于上述细胞分化轨迹，通过 BEAM 分析发现，有 2 477 个基因在两个不同的分化方向表现出不同的表达模式（图 6A）。对差异最大的前 200 个基因进行 KEGG 通路富集分析，结果发现 Toll-like receptor、TGF-beta、Notch 等信号通路参与了原始生殖细胞向雄性化逆转的过程（图6B）；甘氨酸，丝氨酸、苏氨酸、半胱氨酸、甲硫氨酸代谢、氨基糖和核苷酸糖代谢等途径也与该过程相关（图 6C）；不饱和脂肪酸、类固醇激素生物合成，醚酯、亚油酸代谢等途径则对维持原始生殖细胞向雌性性腺分化有重要作用（图 6B）。1323 讨论原始生殖细胞是一类生殖干细胞，具有向不同性别生殖细胞分化的潜能，是鱼类性别决定和分化过程中最重要的一类细胞。其形成方式目前存在两种学说，第一种是先成论，指受精卵带有母源生殖质，随着细胞分裂进行，生殖质被不均等地分配，获得生殖质的细胞将发育成原始生殖细胞，否则发育成体细胞，如斑马鱼（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 的抑A：原始生殖细胞 BEAM 分析结果；B、C：分别为雌性、雄性生殖细胞分化相关基因 KEGG 富集结果A:BEAM analysis of primordial germ cell;B:KEGG enrichment of female germ cell related genes;C:KEGG enrichment of male germ cell related gene图 6 原始生殖细胞 BEAM 分析及基因富集分析结果Fig.6 BEAM analysis of primordial germ cells and enrichment analysis of genes133制作用。促 卵 泡 激 素（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 孕二烯酮溶液暴露 48 周后，发现雌性食蚊鱼臀鳍出现雄性化特征；分别取野生雌性、雄性食蚊鱼和暴露4、8 周雌性食蚊鱼的性腺进行单细胞转录组测序并进行数据分析，结果发现野生雌性食蚊鱼和雄性食蚊鱼性腺细胞构成明显不同，暴露 4、8 周后的雌性食蚊鱼性腺中出现少量的雄性生殖细胞（精原细胞），且原始生殖细胞数量有明显上升；在孕二烯酮诱导下雌性食蚊鱼原始生殖细胞中性别决定关键基因 amh、fabp3、dmrt1、wt1a、wt1b、cyp19a1a 和 fgf16 出现差异化表达情况；细胞轨迹分析表明，部分原始生殖细胞出现向雄性化分化 的 特 征。Toll-like receptor、TGF-beta、Notch和一些氨基酸代谢信号通路参与了原始生殖细胞向雄性化逆转的过程。参考文献（References）：1 NAGAHAMA Y,CHAKRABORTY T,PAUL-PRASANTH B,OHTA K,NAKAMURA M.Sex determination,gonadal sex differentiation,and plasticity in vertebrate speciesJ.Physiological Reviews,2021,101(3):1237-1308.DOI:10.1152/physrev.00044.2019.2 CASTAEDA-CORTS D C,FERNANDINO J I.Stress and sex determination in fish:from brain to gonadsJ.The International Journal of Developmental Biology,2021,65(4/5/6):207-214.DOI:10.1387/ijdb.200072jf.3 YAMAMOTO Y,HATTORI R S,PATIO R,STRSSMANN C A.Environmental regulation of sex determination in fishes:Insights from AtheriniformesJ.Current Topics in Developmental Biology,2019,134:49-69.DOI:10.1016/bs.ctdb.2019.02.003.4 GARCA-CRUZ E L,YAMAMOTO Y,HATTORI R S,DE VASCONCELOS L M,YOKOTA M,STRSSMANN C A.Crowding stress during the period of sex determination causes masculinization in pejerrey Odontesthes bonariensis,a fish with temperature-dependent 134sex determinationJ.Comparative Biochemistry and Physiology.Part A,Molecular&Integrative Physiology,2020,245:110701.DOI:10.1016/j.cbpa.2020.110701.5 KOBAYASHI Y,NAGAHAMA Y,NAKAMURA M.Diversity and plasticity of sex determination and differentiation in