来曲唑与阿特拉津：对泥鳅性别分化基因表达的影响探究

来曲唑与阿特拉津：对泥鳅性别分化基因表达的影响探究一、引言1.1研究背景环境内分泌干扰物（EndocrineDisruptingChemicals，EDCs）是一类能够干扰内分泌系统正常功能，影响激素水平或激素受体活性，进而对生物体生长发育和繁殖产生不良影响的化学物质。其来源广泛，涵盖工业排放、农业活动以及日常生活等多个领域。在工业生产中，塑料制造、石油化工等行业排放的废水和废气含有多种EDCs；农业方面，农药、化肥和兽药的使用，也会释放出EDCs，这些物质在土壤和水体中降解，同样会生成EDCs。由于其具有持久性、生物累积性和生物放大效应，可通过大气、水体和土壤进行长距离传输，并在食物链中积累，对生态系统和人类健康构成了严重威胁。鱼类作为水生生态系统的重要组成部分，其性别分化过程极易受到环境内分泌干扰物的影响。鱼类性别分化是在性别决定机制的控制下，未分化的性腺向卵巢或精巢发育，并出现第二性征的过程，这一过程受到自身内分泌系统的精确调控。而环境内分泌干扰物能够通过多种途径干扰鱼类的内分泌机能，从而对其性别分化产生影响，具体表现为改变鱼类的性别比例、影响配子发生、干扰性腺发育以及第二性征的表达等。在英国的城市污水处理厂下游泻湖以及瑞士城市污水处理厂排水区域，均捕获到具有雌雄两性特征的斜齿鳊；在英国诺福克郡的艾尔河观测点，受工厂排水污染水域中的石斑鱼发生了严重雌化现象，60%的雄性石斑鱼出现雌性化特征，部分雄性鱼生殖器具有排卵功能，甚至转变为两性鱼或雌性鱼。由此可见，排入水中的天然或人工合成环境激素，是导致雌雄同体鱼出现的重要原因。阿特拉津（Atrazine）作为一种广泛应用的三氮苯类除草剂，主要用于林果地和玉米田中杀灭一年生禾本科杂草和阔叶杂草。自20世纪50年代在世界范围内广泛推广应用以来，其使用面积不断扩大。在我国中南部的丘陵和山区地带，玉米和水稻通常一起种植，施入玉米田的阿特拉津可能会经雨水冲刷进入稻田，对水稻造成药害。近年来，研究人员在地表水体如湖泊、池塘以及河流中均检测到阿特拉津。由于阿特拉津性质稳定，不易降解，残效期长，其可能的生态效应受到了人们的高度重视。有研究表明，阿特拉津具有内分泌干扰作用，赛林霖等用阿特拉津暴露青鳉鱼后发现，阿特拉津可导致青鳉发育畸形，精巢发育迟缓，并出现雌雄同体现象，这表明阿特拉津可能的环境激素效应，对水生生物和人类健康构成潜在威胁。来曲唑（Letrozole）是一种芳香化酶抑制剂，在临床上主要用于治疗乳腺癌等疾病。它能够抑制芳香化酶的活性，从而减少雌激素的合成。在鱼类研究中，芳香化酶在性别分化过程中起着关键作用，它能够催化雄激素转化为雌激素，对鱼类性腺的发育和分化具有重要影响。有研究通过投喂芳香化酶抑制剂来曲唑后，研究暗纹东方鲀初孵仔鱼CYP19A、DMRT1基因表达以及性腺的组织学变化，发现对照组样品CYP19A和DMRT1表达显示性二态，雌性表达CYP19A基因，雄性表达DMRT1基因，而LE处理组在性别分化上出现了明显的变化。这表明来曲唑可能通过影响芳香化酶的活性，干扰鱼类体内雌激素的合成，进而对鱼类的性别分化产生影响。泥鳅（Misgurnusanguillicaudatus）作为一种常见的淡水鱼类，具有分布广泛、适应性强等特点，是研究环境内分泌干扰物对鱼类性别分化影响的理想实验材料。目前，关于阿特拉津和来曲唑对泥鳅性别分化相关基因表达模式的影响研究较少，深入探究这两种物质对泥鳅性别分化的作用机制，对于准确评估环境内分泌干扰物的环境毒性，保护水生生物的生殖发育以及维护生态系统的平衡具有重要意义。1.2研究目的与意义本研究旨在通过实验，深入探究来曲唑和阿特拉津这两种环境内分泌干扰物对泥鳅性别分化相关基因表达模式的影响，揭示其作用机制。具体而言，本研究将以泥鳅为实验对象，设置不同浓度梯度的来曲唑和阿特拉津处理组，通过实时荧光定量PCR技术，检测泥鳅性别分化相关基因，如细胞色素P450芳香化酶基因（cyp19a）、类固醇生成因子（sf1）、抗缪勒氏管激素（amh）等在不同处理组中的表达水平，分析来曲唑和阿特拉津对这些基因表达的影响规律。同时，利用组织学切片技术，观察泥鳅性腺在不同处理条件下的形态结构变化，进一步验证基因表达变化与性腺发育之间的关联。本研究具有重要的理论意义和实际应用价值。在理论层面，有助于深化对环境内分泌干扰物影响鱼类性别分化分子机制的理解。目前，虽然已有研究表明环境内分泌干扰物能够影响鱼类性别分化，但对于具体作用机制，尤其是多种环境内分泌干扰物联合作用的机制研究还相对较少。本研究通过对来曲唑和阿特拉津这两种具有代表性的环境内分泌干扰物的研究，能够填补相关领域在泥鳅这一物种上的研究空白，为进一步揭示环境内分泌干扰物对鱼类性别分化的影响机制提供理论依据。从实际应用价值来看，本研究的结果对于评估环境内分泌干扰物的生态风险具有重要意义。阿特拉津作为一种广泛使用的除草剂，来曲唑作为临床常用药物，它们在环境中的残留可能对水生生物产生潜在威胁。通过本研究，能够准确评估这两种物质对泥鳅性别分化的影响，进而为制定合理的环境质量标准和风险评估体系提供科学依据，有助于保护水生生态系统的平衡和稳定，维护生物多样性。此外，本研究结果也为水产养殖行业提供了参考，帮助养殖者更好地了解环境因素对鱼类性别分化的影响，从而采取相应的措施，提高养殖效益和水产品质量。二、研究基础2.1动物性别决定和性别分化动物的性别决定和性别分化是一个复杂而精细的过程，受到多种因素的调控，其机制在不同物种之间存在着显著的差异。性别决定是指生物体在发育过程中确定其性别的过程，它决定了个体将来发育为雄性还是雌性；而性别分化则是在性别决定的基础上，生物体进一步发育出与性别相关的形态、生理和行为特征的过程。2.1.1遗传型性别决定遗传型性别决定是指生物体的性别由基因决定，这是一种较为常见的性别决定方式。在许多动物中，性染色体在性别决定中起着关键作用。以人类为例，正常男性的性染色体组成为XY，而女性为XX。Y染色体上的SRY基因（Sex-determiningRegionofYchromosome）是雄性性别决定的关键基因，它在胚胎发育早期启动一系列基因的表达，促使性腺发育为睾丸。若SRY基因发生突变或缺失，可能导致雄性个体发育出雌性特征。在鱼类中，虽然性染色体的类型和结构相对多样，但也存在类似的基因决定性别方式。如青鳉鱼，其Y染色体上的DMY基因（DoublesexandMab-3relatedtranscriptionfactorontheYchromosome）被证实是性别决定的关键基因，DMY基因的表达能够诱导精巢的发育，使个体发育为雄性。而在尼罗罗非鱼中，性别决定受到多个基因的调控，其中DMRT1基因（DoublesexandMab-3relatedtranscriptionfactor1）在精巢发育过程中起着重要作用，其高表达有助于精巢的形成和分化。此外，一些鱼类还存在性别决定基因的剂量效应，如某些具有ZW/ZZ性别决定系统的鱼类，W染色体上的基因剂量可能影响雌性性腺的发育。2.1.2环境型性别决定环境型性别决定是指环境因素在生物体性别决定过程中起主导作用。温度是影响性别决定的重要环境因素之一，在许多爬行动物中表现得尤为明显。例如，海龟的性别由孵化温度决定，在较低温度（25-27℃）下孵化，多发育为雄性；在较高温度（30-33℃）下孵化，则多发育为雌性。鳄鱼也是典型的温度依赖型性别决定动物，在30℃以下孵化，后代多为雌性；在34℃以上孵化，后代多为雄性。这种温度对性别决定的影响，主要是通过影响性激素合成相关基因的表达来实现的。除温度外，激素等环境因素也能对性别决定产生影响。在鱼类中，一些研究表明，暴露于雌激素或雄激素环境中的幼鱼，其性别分化可能会受到干扰。如将斑马鱼幼鱼暴露于雌激素环境中，会导致雄性个体雌性化，精巢发育受阻，甚至出现卵巢组织；相反，暴露于雄激素环境中的幼鱼，则可能出现雌性个体雄性化现象。此外，一些环境内分泌干扰物，如壬基酚、双酚A等，也具有类似激素的作用，能够干扰鱼类的内分泌系统，进而影响性别分化。这些环境内分泌干扰物可以通过与激素受体结合，激活或抑制相关信号通路，影响性别分化相关基因的表达，最终导致性别分化异常。2.2鱼类性别决定和性别分化2.2.1鱼类的性别决定类型鱼类的性别决定类型丰富多样，主要包括染色体性别决定和基因性别决定，且部分鱼类还具有独特的性别决定方式。在染色体性别决定方面，XY型是较为常见的一种类型，如大部分鲑科鱼类就属于这种类型，其雄性个体具有XY性染色体，雌性个体则为XX。这种类型与哺乳动物的性别决定方式相似，Y染色体在雄性性别决定中起着关键作用。与之相对的是ZW型，以花鲈为代表，雌性个体的性染色体组成为ZW，雄性个体为ZZ，此时雌性个体可以产生含Z或W的两种卵子，而雄性个体产生的精子只含Z染色体。此外，还有一些特殊的染色体性别决定类型，如ZO/ZZ型和XX/XO型，属于雌雄异配型，前者为雌性配子异性，后者为雄性配子异性，短颌鲚就属于这种类型。以及X1X1X2X2/X1X2Y型的复性染色体决定型，例如新月鱼，该类型主要是由于常染色体和性染色体之间以及性染色体之间发生的染色体融合所致。基因性别决定在鱼类中也较为普遍，一些基因在鱼类性别决定中发挥着核心作用。在青鳉鱼中，DMY基因位于Y染色体上，是性别决定的关键基因，其表达能够诱导精巢的发育，使个体发育为雄性；尼罗罗非鱼的性别决定则受到多个基因的调控，其中DMRT1基因在精巢发育过程中起着重要作用，其高表达有助于精巢的形成和分化。斑马鱼的性别决定与gsdf基因密切相关，该基因在雄性斑马鱼中表达，通过抑制sox102a基因的表达，促进男性生殖系统的发育。部分鱼类的性别决定还呈现出独特性和复杂性。一些鱼类没有明显的性染色体分化，其性别决定受多个基因的综合作用。银鲫的性别决定机制较为特殊，它以雌核发育为主，生殖方式独特，性别决定受到多种因素的影响，包括遗传因素和环境因素。还有一些鱼类具有性逆转现象，黄鳝在生长发育过程中会经历从雌性到雄性的性逆转，在幼体阶段，黄鳝主要为雌性，随着生长，部分个体的性腺会逐渐发生转变，最终发育为雄性，这种性逆转现象涉及到复杂的基因调控和内分泌调节机制。2.2.2鱼类性别决定与性别分化相关基因鱼类的性别决定与性别分化是一个受到众多基因精确调控的复杂过程，这些基因之间相互协作、相互制约，共同影响着性腺的发育和分化方向。细胞色素P450芳香化酶基因（cyp19a）在鱼类性别分化中占据着关键地位。它所编码的细胞色素P450芳香化酶能够催化雄激素转化为雌激素，对鱼类性腺的发育和分化具有重要影响。在许多鱼类中，cyp19a基因在卵巢中的表达水平较高，其高表达促使雌激素的合成增加，从而促进卵巢的发育和维持。在斑马鱼中，cyp19a基因的表达与卵巢的发育密切相关，敲低该基因的表达会导致卵巢发育异常，雌性个体出现雄性化特征。在尼罗罗非鱼中，cyp19a基因的表达模式也呈现出明显的性别差异，在雌性个体性腺发育过程中，该基因的表达水平逐渐升高，为卵巢的正常发育提供必要的雌激素环境。类固醇生成因子（sf1）是一种核受体转录因子，在鱼类性别决定和分化过程中发挥着重要作用。sf1基因在性腺、肾上腺等组织中表达，它参与调控类固醇激素的合成和代谢，对性腺的早期发育和分化具有关键作用。在虹鳟鱼中，sf1基因在性腺原基形成初期就开始表达，它能够激活一系列与性腺发育相关的基因，促进性腺的分化。研究还发现，sf1基因可以与其他性别决定相关基因相互作用，共同调节性腺的发育方向。在青鳉鱼中，sf1基因与DMY基因相互作用，协同调控精巢的发育。抗缪勒氏管激素（amh）也是鱼类性别决定和分化过程中的重要基因。amh基因主要在雄性性腺中表达，它能够抑制缪勒氏管的发育，促进雄性生殖器官的形成。在罗非鱼中，amh基因在精巢中的表达水平明显高于卵巢，其表达产物能够抑制雌性生殖器官相关基因的表达，推动性腺向精巢方向分化。此外，amh基因还可以通过与其他信号通路相互作用，调节鱼类的性别分化。在斑马鱼中，amh基因与转化生长因子β（TGF-β）信号通路相互关联，共同影响性腺的发育和性别分化。这些基因之间存在着复杂的相互关系和调控网络。cyp19a基因的表达可能受到sf1基因的调控，sf1可以结合到cyp19a基因的启动子区域，促进其转录表达。而amh基因与cyp19a基因之间则存在着相互抑制的关系，amh基因的表达产物能够抑制cyp19a基因的表达，从而减少雌激素的合成，有利于精巢的发育；反之，雌激素也可以通过反馈调节抑制amh基因的表达。这些基因之间的相互作用，共同构成了鱼类性别决定和分化的复杂调控机制，确保性腺能够按照正常的性别程序发育和分化。2.3环境内分泌干扰物对鱼类性别决定和性别分化的影响环境内分泌干扰物（EDCs）对鱼类性别决定和性别分化的影响机制十分复杂，涉及多个层面，主要包括诱导基因突变、影响基因表达、干扰鱼类激素合成以及与代谢途径相关酶的活性，还有不同雄性或者雌性激素对鱼类产生的雄性或雌性化作用等方面。这些影响可能导致鱼类性别比例失衡、性腺发育异常以及生殖功能受损，进而对水生生态系统的稳定性和生物多样性造成严重威胁。2.3.1诱导基因突变环境内分泌干扰物可以通过多种途径诱导鱼类基因突变，从而影响其性别决定和分化。一些具有亲电性的EDCs能够直接与DNA分子发生反应，导致DNA碱基对的替换、缺失或插入，进而引起基因突变。多环芳烃类物质中的苯并芘，它可以在体内代谢为具有强亲电性的环氧化物，这些环氧化物能够与DNA分子中的鸟嘌呤碱基结合，形成DNA加合物，阻碍DNA的正常复制和转录过程，增加基因突变的概率。如果这些突变发生在性别决定关键基因上，就可能导致鱼类性别决定和分化异常。当EDCs干扰DNA修复机制时，也会间接增加基因突变的风险。正常情况下，细胞内存在一系列DNA修复机制，能够及时识别和修复受损的DNA。然而，某些EDCs会抑制DNA修复酶的活性，或者干扰修复信号通路，使得受损DNA无法及时得到修复。双酚A可以抑制DNA聚合酶β的活性，该酶在DNA碱基切除修复过程中起着关键作用，其活性受到抑制后，DNA损伤难以修复，从而增加了基因突变的可能性。若性别决定相关基因发生突变，如青鳉鱼的DMY基因、尼罗罗非鱼的DMRT1基因等，就可能导致鱼类性别决定和分化的异常，出现性逆转、雌雄同体等现象。2.3.2影响基因表达环境内分泌干扰物主要通过干扰基因转录和翻译过程，改变鱼类性别相关基因的表达水平，进而影响其性别决定和分化。在基因转录过程中，EDCs可以与转录因子相互作用，影响转录因子与基因启动子区域的结合能力，从而调控基因的转录起始。一些具有雌激素活性的EDCs，如壬基酚、双酚A等，能够与雌激素受体结合，形成激素-受体复合物，该复合物可以结合到雌激素反应元件上，激活或抑制相关基因的转录。在斑马鱼中，壬基酚能够与雌激素受体结合，上调cyp19a基因的表达，促进雄激素向雌激素的转化，导致雄性斑马鱼出现雌性化特征。EDCs还可以影响基因转录后的加工过程，如mRNA的剪接、修饰和稳定性等。某些EDCs可能会干扰剪接体的组装或功能，导致mRNA剪接异常，产生错误的转录本。全氟化合物可以影响mRNA的甲基化修饰，改变mRNA的稳定性和翻译效率，从而间接影响基因的表达水平。在基因翻译过程中，EDCs可能干扰核糖体与mRNA的结合，或者影响翻译起始因子和延伸因子的活性，进而影响蛋白质的合成。有研究表明，有机磷农药可以抑制蛋白质合成过程中的关键酶，降低蛋白质的合成速率，影响性别相关蛋白质的表达，最终对鱼类性别分化产生影响。2.3.3干扰鱼类激素合成以及与代谢途径相关酶的活性环境内分泌干扰物能够干扰鱼类激素合成及代谢途径中酶的活性，从而对鱼类性别决定和分化产生影响。在雌激素合成过程中，细胞色素P450芳香化酶（CYP19）起着关键作用，它可以催化雄激素转化为雌激素。一些EDCs，如来曲唑、阿特拉津等，能够抑制CYP19的活性，减少雌激素的合成。来曲唑作为一种芳香化酶抑制剂，能够特异性地与CYP19的活性位点结合，阻断其催化雄激素转化为雌激素的反应，导致鱼类体内雌激素水平下降。在暗纹东方鲀的研究中发现，投喂来曲唑后，CYP19A基因的表达受到抑制，雌激素合成减少，从而影响了性腺的发育和性别分化。在雄激素合成过程中，17β-羟类固醇脱氢酶（17β-HSD）等酶参与了睾酮等雄激素的合成和代谢。某些EDCs可以抑制17β-HSD的活性，阻碍雄激素的合成。多氯联苯（PCBs）能够干扰17β-HSD的活性，使睾酮的合成减少，影响雄性鱼类的性腺发育和性别分化。EDCs还可能影响其他与激素代谢相关酶的活性，如类固醇硫酸酯酶、5α-还原酶等，这些酶参与激素的活化、灭活和代谢过程，其活性受到干扰后，会导致激素水平失衡，进而影响鱼类的性别决定和分化。2.3.4不同雄性或者雌性激素对鱼类产生的雄性或雌性化作用外源雄性或雌性激素对鱼类的性别具有显著影响，能够导致鱼类出现雄性化或雌性化现象。当鱼类暴露于雌激素环境中时，雌激素可以与雌激素受体结合，激活一系列信号通路，促进卵巢相关基因的表达，抑制精巢相关基因的表达，从而使雄性鱼类雌性化。在虹鳟鱼的研究中发现，将雄性虹鳟鱼暴露于17α-乙炔基雌二醇中，会导致其精巢发育受阻，出现卵巢样结构，同时雌性相关基因的表达上调，雄性相关基因的表达下调。雌激素还可以影响鱼类的第二性征，使雄性鱼类的体色、体型等向雌性方向转变。相反，当鱼类暴露于雄激素环境中时，雄激素会与雄激素受体结合，促进精巢相关基因的表达，抑制卵巢相关基因的表达，导致雌性鱼类雄性化。将雌性斑马鱼暴露于17α-甲基睾酮中，会使其卵巢发育受到抑制，出现精巢样结构，雄性相关基因的表达上调，雌性相关基因的表达下调。雄激素也会改变鱼类的第二性征，使雌性鱼类出现雄性特征，如体型增大、鳍条变粗等。这些外源性激素对鱼类性别的影响，不仅会影响个体的生殖能力，还可能对种群的性别比例和遗传多样性产生深远影响。2.4miRNA的研究概述2.4.1miRNA的定义miRNA（MicroRNA）是一类内源性非编码小分子RNA，长度通常在20-24个核苷酸之间。它们广泛存在于动植物、病毒等多种生物体内，虽然不编码蛋白质，但在基因表达调控中发挥着至关重要的作用。miRNA通过与靶mRNA的互补配对，对靶mRNA的翻译过程进行抑制，或者促使其降解，从而实现对基因表达的负调控。在动物中，miRNA主要通过与靶mRNA的3'-非翻译区（3'-UTR）相互作用，抑制mRNA的翻译起始或者导致mRNA的降解。在植物中，miRNA与靶mRNA的互补配对更为精确，往往能直接介导mRNA的切割降解。miRNA的这种调控作用，使得生物体能够在转录后水平对基因表达进行精细调节，从而适应不同的生理和环境需求。2.4.2miRNA的生物发生和作用机制miRNA的生物发生是一个复杂而有序的过程，涉及多个关键步骤和酶的参与。首先，在细胞核内，miRNA基因在RNA聚合酶Ⅱ的作用下转录生成初级miRNA（pri-miRNA），pri-miRNA通常具有较长的序列，并且包含多个茎环结构。随后，pri-miRNA在核酸酶Drosha及其辅助因子DGCR8的作用下，被切割成约70-100个核苷酸长度的前体miRNA（pre-miRNA），pre-miRNA呈发夹状结构。pre-miRNA通过转运蛋白Exportin-5从细胞核转运到细胞质中。在细胞质中，pre-miRNA在核酸酶Dicer的作用下，进一步被切割成约22个核苷酸的双链miRNA。双链miRNA中的一条链会被选择性地保留，形成成熟的单链miRNA，另一条链则被降解。成熟的miRNA会与AGO（Argonaute）蛋白等组成RNA诱导沉默复合体（RISC）。RISC中的miRNA通过碱基互补配对原则，识别并结合到靶mRNA的3'-UTR区域。如果miRNA与靶mRNA的互补配对程度较高，RISC中的AGO蛋白会切割靶mRNA，导致其降解；若互补配对程度较低，RISC则主要抑制靶mRNA的翻译过程，阻止蛋白质的合成。在动物细胞中，当miRNA与靶mRNA不完全互补配对时，RISC主要通过抑制翻译起始、阻碍核糖体的延伸或者促进mRNA的去腺苷酸化等方式，抑制靶mRNA的翻译。在植物细胞中，由于miRNA与靶mRNA的互补配对程度通常较高，所以RISC主要通过切割靶mRNA来实现对基因表达的调控。这种复杂而精细的生物发生和作用机制，使得miRNA能够高效、精准地调控基因表达。2.4.3miRNA的功能miRNA在生物体的细胞增殖、分化、凋亡以及个体发育等多个过程中发挥着关键的调控作用。在细胞增殖过程中，miRNA可以通过调控相关基因的表达，影响细胞周期的进程。miR-15a和miR-16-1能够通过靶向抑制细胞周期蛋白D1（CyclinD1）等基因的表达，将细胞周期阻滞在G0/G1期，从而抑制细胞增殖。在细胞分化方面，miRNA参与了多种细胞类型的分化调控。在神经干细胞分化过程中，miR-9等miRNA可以通过调控神经分化相关基因的表达，促进神经干细胞向神经元分化。在细胞凋亡过程中，miRNA也发挥着重要作用。miR-125b可以通过抑制促凋亡基因Bax的表达，减少细胞凋亡；而miR-15a和miR-16-1则可以通过靶向抑制抗凋亡基因Bcl-2的表达，促进细胞凋亡。在个体发育过程中，miRNA对生物体的正常发育至关重要。在果蝇的胚胎发育过程中，miR-34等miRNA参与了胚胎体节的形成和器官的发育调控；在小鼠的胚胎发育过程中，miR-17-92簇等miRNA对心脏、肺等器官的发育起着关键作用。在鱼类性腺发育方面，miRNA同样发挥着重要的调控作用。研究发现，一些miRNA在鱼类性腺发育的不同阶段呈现出特异性表达。在斑马鱼中，miR-218在卵巢发育过程中表达上调，通过抑制其靶基因的表达，促进卵巢的发育和成熟；而在精巢发育过程中，miR-138的表达上调，对精巢的发育和精子发生具有重要的调控作用。在尼罗罗非鱼中，miR-143等miRNA在性腺分化过程中发挥着关键作用，它们通过调控性别分化相关基因的表达，影响性腺的发育方向。这些研究表明，miRNA在鱼类性腺发育过程中通过调控相关基因的表达，维持性腺发育的正常进程，对鱼类的生殖功能具有重要意义。2.4.4miRNA的研究进展在鱼类性别分化领域，miRNA的研究取得了一系列重要成果。随着高通量测序技术和生物信息学分析方法的不断发展，越来越多与鱼类性别分化相关的miRNA被鉴定和研究。在半滑舌鳎中，通过对不同性别个体性腺的miRNA测序分析，发现了多个差异表达的miRNA，其中miR-135b等miRNA在雌性性腺中的表达显著高于雄性，进一步研究表明，miR-135b可以通过靶向抑制雄性性别决定相关基因的表达，促进性腺向雌性方向分化。在虹鳟鱼中，研究人员发现miR-202等miRNA在精巢发育过程中发挥着重要作用，它们通过调控精巢发育相关基因的表达，影响精子的发生和成熟。尽管miRNA在鱼类性别分化领域取得了一定的研究成果，但目前的研究仍存在一些问题。对于大多数已鉴定的miRNA，其具体的靶基因和调控网络尚未完全明确。在研究miRNA与靶基因的相互作用时，往往存在假阳性结果，需要进一步验证。此外，miRNA在鱼类性别分化过程中的时空表达模式和动态调控机制还需要深入研究。未来，该领域的研究方向主要包括以下几个方面。利用多种技术手段，如荧光素酶报告基因实验、RNA免疫沉淀等，深入研究miRNA与靶基因之间的相互作用，明确其调控网络。结合单细胞测序、活体成像等技术，研究miRNA在鱼类性别分化过程中的时空表达模式和动态调控机制。开展不同鱼类物种间miRNA的比较研究，探讨miRNA在鱼类性别分化过程中的进化保守性和特异性。将miRNA的研究成果应用于水产养殖实践，通过调控miRNA的表达，实现对鱼类性别比例的控制，提高养殖效益。三、实验研究3.1实验材料3.1.1实验动物本实验选用泥鳅作为研究对象，泥鳅广泛分布于各类淡水水域，对环境适应能力强，生长速度较快，性成熟周期相对较短，且在实验室条件下易于饲养和繁殖。这些特性使得泥鳅成为研究环境内分泌干扰物对鱼类性别分化影响的理想实验材料。实验所用泥鳅购自[具体供应商名称]，均为同一批次孵化的幼鱼，初始体长为[X]cm，体重为[X]g。泥鳅运回实验室后，暂养于玻璃水族箱中，水族箱规格为[长×宽×高，单位：cm]，养殖用水为曝气24h以上的自来水，水温控制在（25±1）℃，光照周期设置为12h光照：12h黑暗。每日投喂[饲料品牌及规格]饲料2次，投喂量以泥鳅在15min内吃完为宜，及时清理残饵和粪便，以保持水质清洁。暂养7d，待泥鳅适应实验室环境后，挑选健康、活力强的个体用于后续实验。3.1.2实验仪器本实验所需的主要仪器如下：仪器名称型号用途PCR仪T100（伯乐生命医学产品（上海）有限公司）用于基因扩增，通过设定特定的温度循环，实现对泥鳅性别分化相关基因的体外扩增，以便后续检测分析荧光定量PCR仪LightCycler480（德国ROCHE公司）用于对扩增后的基因进行荧光定量检测，精确测定基因的表达量，分析不同处理组间基因表达的差异切片机RM2235（徕卡显微系统（上海）有限公司）用于制作泥鳅性腺组织切片，将固定好的性腺组织切成薄片，以便进行组织学观察显微镜BX53（奥林巴斯（中国）有限公司）用于观察泥鳅性腺组织切片的形态结构，结合图像分析软件，对性腺发育情况进行评估高速冷冻离心机5424R（德国艾本德股份公司）用于离心分离细胞、蛋白质、核酸等生物样品，在RNA提取等实验步骤中，实现样品的分离和纯化移液器Researchplus（德国艾本德股份公司）准确移取不同体积的液体试剂，确保实验操作的准确性和重复性电子天平AL204（梅特勒-托利多仪器（上海）有限公司）用于精确称量实验所需的各种试剂和药品，保证实验试剂的用量准确恒温培养箱DNP-9082（上海精宏实验设备有限公司）为实验提供恒定的温度环境，满足细胞培养、细菌培养等实验需求纯水仪Milli-QAdvantageA10（美国默克密理博公司）制备高纯度的去离子水，用于实验试剂的配制和实验仪器的清洗，保证实验用水的质量3.1.3主要实验试剂本实验所用的主要试剂如下：试剂名称来源规格来曲唑上海源叶生物科技有限公司纯度≥98%阿特拉津北京索莱宝科技有限公司纯度≥99%RNA提取试剂TRIzolReagent（美国Invitrogen公司）500mL反转录试剂PrimeScriptRTreagentKitwithgDNAEraser（宝日医生物技术（北京）有限公司）50次反应荧光定量PCR试剂SYBRPremixExTaq™Ⅱ（宝日医生物技术（北京）有限公司）500μL无水乙醇国药集团化学试剂有限公司分析纯，500mL氯仿国药集团化学试剂有限公司分析纯，500mL异丙醇国药集团化学试剂有限公司分析纯，500mLDEPC水北京索莱宝科技有限公司500mL蛋白酶K北京索莱宝科技有限公司10mg/mL，1mL琼脂糖西班牙Biowest公司500gDNAMarker宝日医生物技术（北京）有限公司500μL溴化乙锭（EB）北京索莱宝科技有限公司10mg/mL，1mL3.1.4引物根据NCBI数据库中已公布的泥鳅性别分化相关基因（如cyp19a、sf1、amh等）以及内参基因β-actin的序列，利用PrimerPremier5.0软件设计引物。引物由[引物合成公司名称]合成，其序列如下：基因名称引物序列（5'-3'）cyp19aF：ATGGCAGAGCCTGTAGAGGAR：TGTGGTGGTGGTAGTGGTAGsf1F：CCGAGCTGCTGAAGAAGAGAR：TGGCTTGGCTTCTTCTCTTGamhF：AGCGACGACGACGACAAGTAR：TGCTGCTGCTGCTGCTGATCβ-actinF：CGTGACATTAAGGAGAAGCTGR：GCCAGAGGCGTACAGGGATA引物设计依据基因的保守区域，确保引物的特异性和扩增效率。引物设计完成后，通过BLAST比对分析，确认引物与其他基因无明显同源性，以保证扩增的特异性。同时，对引物进行了预实验验证，通过PCR扩增和琼脂糖凝胶电泳检测，观察扩增条带的特异性和亮度，结果表明引物能够特异性地扩增目的基因，扩增条带单一、清晰，无明显非特异性扩增产物，满足实验要求。3.2实验方法3.2.1实验设计本实验设置1个对照组和6个处理组，每组设置3个平行，每个平行放养30尾泥鳅。对照组用曝气24h以上的自来水养殖；来曲唑处理组分别设置低、中、高3个浓度梯度，即10μg/L、50μg/L、100μg/L，分别记为LT1、LT2、LT3组；阿特拉津处理组同样设置低、中、高3个浓度梯度，即50μg/L、100μg/L、200μg/L，分别记为AT1、AT2、AT3组。将挑选好的健康泥鳅随机放入各实验组的水族箱中，水族箱规格为[长×宽×高，单位：cm]，养殖用水为曝气24h以上的自来水，水温控制在（25±1）℃，光照周期设置为12h光照：12h黑暗。每日投喂[饲料品牌及规格]饲料2次，投喂量以泥鳅在15min内吃完为宜，及时清理残饵和粪便，以保持水质清洁。处理时间为60d，在处理期间，每隔7d测定一次水质指标，包括pH值、溶解氧、氨氮等，确保水质符合泥鳅生长要求。3.2.2总RNA的提取在实验结束后，从每个实验组中随机选取5尾泥鳅，用MS-222麻醉后，迅速取出脑和性腺组织，放入液氮中速冻，然后保存于-80℃冰箱备用。采用TRIzol试剂提取总RNA，具体步骤如下：将组织放入预冷的研钵中，加入液氮迅速研磨成粉末状。将研磨好的组织粉末转移至1.5mL离心管中，加入1mLTRIzol试剂，充分混匀，室温静置5min，使组织充分裂解。加入200μL氯仿，剧烈振荡15s，室温静置3min。4℃、12000r/min离心15min，此时溶液分为三层，上层为无色透明的水相，含有RNA；中层为白色的蛋白质层；下层为红色的有机相，含有DNA和蛋白质等杂质。将上层水相转移至新的1.5mL离心管中，加入等体积的异丙醇，轻轻混匀，室温静置10min。4℃、12000r/min离心10min，此时RNA会沉淀在离心管底部，形成白色或透明的沉淀。弃去上清液，加入1mL75%乙醇（用DEPC水配制），轻轻洗涤RNA沉淀，4℃、7500r/min离心5min。弃去上清液，将离心管倒扣在滤纸上，晾干RNA沉淀，注意不要让沉淀完全干燥，以免影响RNA的溶解。加入适量的DEPC水，溶解RNA沉淀，用核酸蛋白测定仪测定RNA的浓度和纯度，确保OD260/OD280在1.8-2.0之间，然后将RNA保存于-80℃冰箱备用。在提取过程中，要注意避免RNA酶的污染，所有操作均需在无RNA酶的环境中进行，使用的耗材和试剂均需经过DEPC处理。操作过程中要戴口罩、手套，避免唾液和汗液中的RNA酶污染样品。提取的RNA应尽快进行后续实验，如需长期保存，应保存在-80℃冰箱中，避免反复冻融。3.2.3普通反转录采用PrimeScriptRTreagentKitwithgDNAEraser反转录试剂盒将总RNA反转录为cDNA，其原理是利用反转录酶以RNA为模板，合成与之互补的cDNA。操作步骤如下：在冰上配制DNA去除反应体系，体系总体积为10μL，包括5×gDNAEraserBuffer2μL，gDNAEraser1μL，TotalRNA1μg，RNaseFreedH2O补齐至10μL。将上述反应体系轻轻混匀，42℃孵育2min，以去除基因组DNA。配制反转录反应体系，总体积为20μL，在上述反应体系中加入5×PrimeScriptBuffer24μL，PrimeScriptRTEnzymeMixI1μL，RTPrimerMix1μL，RNaseFreedH2O4μL。将反应体系轻轻混匀，37℃孵育15min，进行反转录反应，然后85℃孵育5s，使反转录酶失活。反转录得到的cDNA保存于-20℃冰箱备用。3.2.4q-RTPCR利用荧光定量PCR技术检测基因表达水平，其原理是在PCR反应体系中加入荧光基团，随着PCR反应的进行，荧光信号强度与PCR产物的数量成正比，通过检测荧光信号的变化，实时监测PCR反应进程，从而对基因表达水平进行定量分析。使用设计好的引物（具体引物序列见3.1.4），以β-actin作为内参基因，进行荧光定量PCR扩增。反应体系为20μL，包括SYBRPremixExTaq™Ⅱ10μL，上下游引物（10μM）各0.8μL，cDNA模板2μL，ddH2O6.4μL。反应程序设置为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。在反应过程中，利用荧光定量PCR仪实时采集荧光信号，每个样品设置3个技术重复。数据采集和分析方法如下：采用2-ΔΔCt法计算基因的相对表达量，ΔCt=Ct目的基因-Ct内参基因，ΔΔCt=ΔCt处理组-ΔCt对照组，相对表达量=2-ΔΔCt。利用GraphPadPrism8.0软件对数据进行统计分析，采用单因素方差分析（One-wayANOVA）和Dunnett's多重比较检验不同处理组与对照组之间基因表达水平的差异，以P<0.05作为差异具有统计学意义的标准。3.2.5石蜡切片在实验结束后，从每个实验组中随机选取5尾泥鳅，用MS-222麻醉后，迅速取出性腺组织，放入Bouin氏固定液中固定24h。固定后的组织依次用70%、80%、90%、95%、100%的乙醇进行脱水处理，每个浓度的乙醇处理时间为1-2h，以去除组织中的水分。脱水后的组织用二甲苯进行透明处理，二甲苯处理时间为30min-1h，使组织变得透明，便于后续包埋。将透明后的组织放入融化的石蜡中进行包埋，包埋温度为60℃左右，包埋时间为1-2h，使石蜡充分渗透到组织中。用切片机将包埋好的组织切成厚度为5μm的切片，将切片贴附在载玻片上。将载玻片放入60℃烘箱中烤片1-2h，使切片牢固地贴附在载玻片上。采用苏木精-伊红（HE）染色法对切片进行染色，染色步骤如下：将切片依次放入二甲苯Ⅰ、二甲苯Ⅱ中脱蜡10min；然后依次放入100%、95%、90%、80%、70%的乙醇中进行水化，每个浓度的乙醇处理时间为3-5min；将切片放入苏木精染液中染色5-10min，自来水冲洗10min；用1%盐酸乙醇分化数秒，自来水冲洗返蓝；将切片放入伊红染液中染色3-5min；依次用70%、80%、90%、95%、100%的乙醇进行脱水，每个浓度的乙醇处理时间为3-5min；最后用二甲苯Ⅰ、二甲苯Ⅱ透明10min，中性树胶封片。将染色后的切片置于显微镜下观察，结合图像分析软件，观察性腺的形态结构变化，如性腺的发育程度、细胞形态、组织结构等，并拍照记录。通过比较不同处理组性腺的形态结构差异，分析来曲唑和阿特拉津对泥鳅性腺发育的影响。3.3实验结果3.3.1泥鳅脑和性腺总RNA的提取使用核酸蛋白测定仪对提取的总RNA进行浓度和纯度检测，结果显示，所有样本的OD260/OD280比值均在1.8-2.0之间，表明提取的RNA纯度较高，无蛋白质和酚类等杂质污染。通过琼脂糖凝胶电泳检测RNA的完整性，结果显示，28SrRNA和18SrRNA条带清晰，亮度比值约为2:1，且无明显的弥散现象，说明提取的RNA完整性良好，无明显降解。本实验提取的泥鳅脑和性腺总RNA的纯度和完整性均满足后续实验要求，可用于反转录和荧光定量PCR等实验。3.3.2组织学分析对照组泥鳅的性腺发育正常，精巢中可见大量成熟的精子，生精小管结构清晰，管壁由精原细胞、精母细胞和精子细胞等组成；卵巢中可见不同发育阶段的卵母细胞，包括初级卵母细胞、次级卵母细胞和成熟的卵细胞，卵泡结构完整。在来曲唑处理组中，随着来曲唑浓度的升高，精巢的发育受到抑制，生精小管的管径变小，管腔内精子数量减少，部分生精小管出现萎缩现象；卵巢中卵母细胞的发育也受到影响，初级卵母细胞数量减少，成熟卵细胞的比例降低，部分卵巢组织出现退化迹象。阿特拉津处理组中，性腺发育同样出现异常。精巢中生精小管的结构紊乱，生精细胞排列不规则，精子发生受到阻碍，出现多核精子和畸形精子；卵巢中卵母细胞的形态发生改变，出现细胞核固缩、细胞质空泡化等现象，部分卵巢组织中还出现了类似精巢组织的结构，表明阿特拉津可能导致泥鳅性腺发生性逆转。3.3.3基因表达分析通过荧光定量PCR技术检测来曲唑和阿特拉津处理后泥鳅性别分化相关基因（cyp19a、sf1、amh等）在mRNA水平的表达变化，结果以图表形式呈现。从图中可以看出，与对照组相比，来曲唑处理组中cyp19a基因的表达量显著下调，且随着来曲唑浓度的升高，下调趋势更加明显，在LT3组中，cyp19a基因的表达量仅为对照组的0.3倍；sf1基因的表达量在LT1组中略有下降，但差异不显著，在LT2和LT3组中显著下降，分别为对照组的0.7倍和0.5倍；amh基因的表达量在来曲唑处理组中显著上调，在LT3组中，amh基因的表达量是对照组的2.5倍。阿特拉津处理组中，cyp19a基因的表达量同样显著下调，在AT3组中，其表达量仅为对照组的0.2倍；sf1基因的表达量在低浓度阿特拉津处理组（AT1）中无明显变化，在中高浓度处理组（AT2、AT3）中显著下降，分别为对照组的0.6倍和0.4倍；amh基因的表达量在阿特拉津处理组中显著上调，在AT3组中，amh基因的表达量是对照组的3倍。对不同浓度处理组与对照组之间基因表达的差异进行显著性分析，结果表明，来曲唑和阿特拉津处理组中cyp19a、sf1和amh基因的表达量与对照组相比，均存在显著差异（P<0.05）。四、深入分析4.1来曲唑和阿特拉津对泥鳅性腺miRNA的表达的影响及功能分析4.1.1实验方法在实验结束后，从对照组、来曲唑处理组（LT3组）和阿特拉津处理组（AT3组）中，每组随机选取5尾泥鳅，迅速解剖取出性腺组织，将性腺组织立即放入液氮中速冻，随后转移至-80℃冰箱保存，以防止RNA降解。采用TRIzol试剂提取性腺组织中的总RNA，具体操作步骤与3.2.2中总RNA提取步骤一致。在提取过程中，严格遵守无RNA酶操作规范，以确保RNA的质量。提取完成后，使用核酸蛋白测定仪测定RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.0之间，以保证RNA的高质量，满足后续实验需求。对于miRNA的反转录，选用茎环法进行。根据miRBase数据库中已有的泥鳅miRNA序列，运用PrimerPremier5.0软件设计茎环逆转录引物。茎环逆转录引物由茎环结构和miRNA的3＇末端六个反向互补的碱基组成。反转录反应使用PrimeScriptRTreagentKitwithgDNAEraser反转录试剂盒，在冰上配制反应体系。反应体系总体积为20μL，包含5×PrimeScriptBuffer24μL，PrimeScriptRTEnzymeMixI1μL，茎环逆转录引物1μL，TotalRNA1μg，RNaseFreedH2O补齐至20μL。将反应体系轻轻混匀后，置于37℃孵育15min，进行反转录反应，然后85℃孵育5s，使反转录酶失活。得到的cDNA保存于-20℃冰箱备用。使用SYBRPremixExTaq™Ⅱ试剂进行miRNA荧光定量PCR扩增。反应体系为20μL，包含SYBRPremixExTaq™Ⅱ10μL，上下游引物（10μM）各0.8μL，cDNA模板2μL，ddH2O6.4μL。反应程序设置为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3个技术重复，以提高实验的准确性和可靠性。在反应过程中，利用荧光定量PCR仪实时采集荧光信号，用于后续的数据分析。使用IlluminaHiSeq2500测序平台对性腺组织的miRNA进行测序。将提取的总RNA进行质量检测，确保RNA的完整性和纯度符合测序要求。然后，构建miRNA文库，采用TruSeqSmallRNALibraryPreparationKit试剂盒进行文库构建。具体步骤包括RNA片段化、接头连接、cDNA合成等。将构建好的文库进行质量评估，使用Agilent2100Bioanalyzer检测文库的插入片段大小和浓度，确保文库质量良好。最后，将合格的文库在IlluminaHiSeq2500测序平台上进行测序，采用PE150测序模式，获得高质量的测序数据。对测序得到的原始数据进行预处理，去除低质量reads、接头序列和长度小于18nt的reads，得到高质量的cleanreads。将cleanreads与miRBase数据库进行比对，鉴定已知的miRNA；同时，利用miRDeep2软件预测新的miRNA。对于已知miRNA的鉴定，通过比对确定其在miRBase数据库中的同源序列，并分析其表达水平。对于新miRNA的预测，miRDeep2软件根据miRNA的生物发生特征，如发夹结构、Dicer酶切割位点等，对cleanreads进行分析，预测潜在的新miRNA。利用DESeq2软件对不同样品间miRNA的表达水平进行差异分析。以|log2(FoldChange)|≥1且P<0.05作为筛选差异表达miRNA的标准。将差异表达的miRNA进行聚类分析，使用Cluster3.0软件绘制热图，直观展示不同样品间miRNA表达的差异模式，以便更清晰地了解不同处理组中miRNA表达的变化情况。随机选取5个差异表达的miRNA，利用荧光定量PCR技术对其表达量进行验证。验证过程中，反应体系和反应程序与上述miRNA荧光定量PCR扩增一致。将荧光定量PCR结果与测序数据进行相关性分析，使用GraphPadPrism8.0软件计算Pearson相关系数，以评估测序数据的准确性和可靠性。使用TargetScan和miRanda软件对差异表达的miRNA进行靶基因预测。将两个软件预测结果的交集作为最终的靶基因集合。对靶基因进行基因本体论（GO）富集分析和京都基因与基因组百科全书（KEGG）通路富集分析，使用DAVID数据库进行分析。GO富集分析从生物过程、细胞组分和分子功能三个层面分析靶基因的功能；KEGG通路富集分析确定靶基因参与的主要信号通路，从而深入了解差异表达miRNA在泥鳅性别分化过程中的潜在调控机制。4.1.2实验结果成功构建了泥鳅精巢和卵巢的miRNA文库，经Agilent2100Bioanalyzer检测，文库的插入片段大小主要分布在18-25nt之间，符合miRNA的长度特征。文库的浓度均大于10nM，且文库的质量良好，无明显降解和杂质污染，表明文库构建成功，可用于后续的测序分析。测序共获得了[X]条高质量的cleanreads，为后续的miRNA鉴定和分析提供了充足的数据基础。通过与miRBase数据库比对，在所有样品中鉴定出了[X]种已知的miRNA。不同样品中已知miRNA的表达水平存在差异，部分已知miRNA在精巢和卵巢中的表达具有明显的性别特异性。如miR-122在卵巢中的表达水平显著高于精巢，而miR-138在精巢中的表达水平明显高于卵巢，这些结果与已有的研究报道相符，进一步验证了实验方法的可靠性。利用miRDeep2软件预测得到了[X]种新的miRNA。新miRNA的发夹结构预测结果显示，它们均具有典型的miRNA前体发夹结构，茎环结构稳定，符合miRNA的生物发生特征。部分新miRNA在不同处理组中的表达也存在差异，暗示它们可能在泥鳅性别分化过程中发挥重要作用。通过DESeq2软件分析，共筛选出[X]个差异表达的miRNA。其中，与对照组相比，来曲唑处理组（LT3组）中有[X]个miRNA表达上调，[X]个miRNA表达下调；阿特拉津处理组（AT3组）中有[X]个miRNA表达上调，[X]个miRNA表达下调。聚类分析结果显示，不同处理组的miRNA表达模式存在明显差异，对照组、来曲唑处理组和阿特拉津处理组各自聚为一类，表明来曲唑和阿特拉津处理显著改变了泥鳅性腺miRNA的表达谱。随机选取的5个差异表达miRNA的荧光定量PCR验证结果与测序数据具有良好的相关性，Pearson相关系数均大于0.8，表明测序数据准确可靠，能够真实反映miRNA的表达情况。这为后续基于测序数据进行的靶基因预测和功能分析提供了有力的支持。通过TargetScan和miRanda软件预测，共得到差异表达miRNA的靶基因[X]个。GO富集分析结果显示，靶基因在生物过程中主要富集在细胞增殖、分化、信号转导等功能；在细胞组分中主要富集在细胞核、细胞膜等部位；在分子功能中主要富集在核酸结合、转录因子活性等方面。KEGG通路富集分析结果表明，靶基因主要参与了PI3K-Akt信号通路、MAPK信号通路、TGF-β信号通路等与性别分化密切相关的信号通路，这些信号通路在调控性腺发育和性别分化过程中发挥着关键作用。4.1.3讨论来曲唑和阿特拉津处理后，泥鳅性腺miRNA表达谱发生了显著变化，这与之前研究中观察到的性别分化相关基因表达变化存在一定的关联。在基因表达层面，来曲唑作为芳香化酶抑制剂，抑制了雌激素的合成，导致cyp19a基因表达下调，进而影响了性腺的发育和分化。阿特拉津也表现出对雌激素合成相关基因的干扰作用，使cyp19a基因表达下降。而在miRNA层面，来曲唑和阿特拉津处理组中多个miRNA的表达发生改变，这些miRNA可能通过调控性别分化相关基因的表达，在来曲唑和阿特拉津对泥鳅性别分化的影响中发挥重要作用。差异表达的miRNA及其靶基因在泥鳅性别分化过程中可能具有重要的潜在调控作用和