短期饥饿对锦鲫肌肉和肝脏生物钟基因表达的节律性重塑研究

短期饥饿对锦鲫肌肉和肝脏生物钟基因表达的节律性重塑研究一、引言1.1研究背景生物钟作为生物体内一种无形的“时钟”，是生物体生命活动的内在节律性体现，由生物体内的时间结构序所决定。它在生物的生理和行为调节中发挥着关键作用，能够调节生物的睡眠-觉醒周期、体温、心率、血压等生理节律，以及觅食、繁殖、迁徙等行为活动，还能调控生物的代谢过程，如能量消耗、脂肪合成等。对生物钟的深入研究在医学领域意义重大，催生了时辰生物学、时辰药理学和时辰治疗学等新学科。2017年，美国科学家杰弗里・霍尔、迈克尔・罗斯巴什和迈克尔・扬因在生物钟运行分子机制方面的卓越研究成就，荣获诺贝尔生理学或医学奖，这进一步彰显了生物钟研究的重要性和影响力。在鱼类中，生物钟同样起着不可或缺的作用。鱼类的生物钟受到日照和季节变化的显著影响，以调节其适应水中生活的各种需求。这使得它们能够在不同的环境条件下顺利生存和繁衍后代。鱼类的视网膜光感受细胞和脑内生物钟核共同构成了生物钟的主要调控机制，它们对光线的强度和颜色有着不同的反应，进而精准调节鱼类的睡眠时间、繁殖行为和觅食行为等。例如，在繁殖季节，一些鱼类会依据生物钟的节律，选择特定的时间和地点进行产卵，以确保后代的生存和繁衍；在觅食方面，鱼类会根据生物钟的指示，在食物资源丰富的时间段积极觅食，提高能量获取效率。研究鱼类生物钟具有多方面的重要意义。从渔业资源管理角度来看，深入了解鱼类的昼夜活动规律，有助于制定更为科学合理的捕捞策略，有效避免过度捕捞，从而切实保护渔业资源的可持续发展。在水产养殖业中，依据鱼类生物钟的特点，可以优化养殖环境，精准调整投喂时间、水温控制等关键因素，提高养殖效率和产量，推动水产养殖业的健康发展。在环境保护领域，鱼类生物钟的研究能够为水体生态系统的保护和管理提供坚实的科学依据，助力实现生态环境的可持续发展目标。此外，由于鱼类与人类在生理机制上存在一定的相似性，研究鱼类生物钟还有助于揭示动物昼夜节律的奥秘，为人类生理机制的研究提供全新的视角和思路，推动生物学领域的整体发展。锦鲫作为一种广泛分布且具有重要经济价值的鲤科鱼类，对各种生态环境具有很强的适应能力，从亚寒带到热带，不论水体深浅、流水或静水、清水或浊水、低氧、酸、碱等环境均能生存繁衍。其生长速度较慢，当年仅长至50克，第二年长到100-150克，第三年可达300克以上，最大个体可达2.5公斤。锦鲫是杂食性鱼类，食谱极为广泛，动物性食物包括枝角类、桡足类、苔藓虫、轮虫、淡水壳菜、蚬、摇蚊幼虫以及虾等，植物性食物则以植物碎屑、硅藻类、丝状藻类、水草等为主。在自然环境中，锦鲫也会面临食物资源分布不均、季节更替或环境剧变等情况，从而遭遇饥饿胁迫。饥饿是鱼类在自然水域生态系中经常面临的一种生理胁迫现象，也是影响鱼类正常生长、发育及生存的一个重要环境因子。不同种类的鱼对饥饿的耐受力和适应性特征存在差异。当鱼类处于饥饿状态时，其代谢水平会显著下降，身体组成、形态结构、行为以及繁殖习性等方面都会发生一系列变化。例如，一些鱼类在饥饿期间体重会急剧下降，出现“负增长”，体长可能无变化甚至有所增加；外形可能会发生萎缩，身体消瘦，头大、体小，身体显得不均称；行为上，摄食行为和相互残食行为等也会受到影响。研究饥饿对鱼类的影响，对于了解鱼类适应饥饿胁迫的生态对策具有重要的理论意义，同时对渔业资源管理、苗种培育、水产养殖等实践活动也有着至关重要的指导作用。目前，关于鱼类生物钟的研究已取得了一定进展，但对于短期饥饿如何影响鱼类生物钟相关基因表达的研究还相对较少，尤其是在锦鲫这一物种上。深入探究短期饥饿对锦鲫肌肉和肝脏中生物钟相关基因表达的影响，不仅可以填补这一研究领域的空白，丰富对鱼类生物钟调控机制的认识，还能为锦鲫的养殖管理提供更为科学的理论依据，具有重要的理论和实践价值。1.2锦鲫生物学特性锦鲫，作为鲤科鱼类中的重要一员，具有极高的研究价值和经济价值，在水产养殖领域占据着独特的地位。锦鲫对环境的适应能力堪称卓越，这是其成为理想研究对象的关键因素之一。从温度适应范围来看，无论是亚寒带的低温水域，还是热带的高温环境，锦鲫都能安然生存并繁衍后代。在低温环境中，其生理机能会进行适应性调整，降低代谢速率以减少能量消耗，同时增强自身的抗寒能力；在高温环境下，锦鲫则通过调节呼吸和排泄系统，维持体内的水分和电解质平衡，确保自身的正常生理活动。在水体条件方面，锦鲫同样表现出强大的适应性。无论是流水潺潺的江河，还是平静无波的湖泊、池塘；无论是清澈见底的清水，还是略显浑浊的水体；甚至是低氧、酸性或碱性较强的特殊水体环境，锦鲫都能顽强生存。这种广泛的环境适应性使得锦鲫在全球范围内广泛分布，为其种群的繁衍和壮大提供了坚实的基础。锦鲫的生长特性也十分显著。在生长速度上，锦鲫相对较慢。在适宜的养殖条件下，锦鲫在出生后的第一年，体重仅能增长至50克左右；到了第二年，其体重会进一步增加到100-150克；而到了第三年，锦鲫的体重可超过300克。尽管生长速度不快，但锦鲫在生长过程中，身体各部分的发育却较为均衡，这使得其在形态上保持着相对稳定的特征。锦鲫的最大个体可达2.5公斤，这样的生长潜力为其在水产养殖中的应用提供了更多的可能性。作为杂食性鱼类，锦鲫的食性特点决定了其在生态系统中的重要地位。在食物选择上，锦鲫的食谱极为广泛，涵盖了动物性食物和植物性食物两大类别。动物性食物方面，锦鲫主要以枝角类、桡足类、苔藓虫、轮虫、淡水壳菜、蚬、摇蚊幼虫以及虾等小型水生生物为食。这些动物性食物富含蛋白质和脂肪，为锦鲫的生长和发育提供了丰富的营养来源。植物性食物则以植物碎屑、硅藻类、丝状藻类、水草等为主。这些植物性食物不仅提供了碳水化合物等能量物质，还富含维生素和矿物质，有助于维持锦鲫的身体健康。在不同的生长阶段和环境条件下，锦鲫对食物的偏好会有所不同。在幼鱼阶段，由于消化系统尚未完全发育成熟，锦鲫更倾向于摄食小型的浮游生物和藻类；随着个体的逐渐长大，锦鲫对大型水生植物和底栖生物的摄食比例会逐渐增加。在自然环境中，锦鲫常常面临着食物资源分布不均、季节更替或环境剧变等诸多挑战，饥饿胁迫也成为了其生存过程中必须面对的问题。当遭遇饥饿时，锦鲫会启动一系列复杂的生理和行为适应性机制。在生理方面，锦鲫的代谢水平会显著下降，以减少能量的消耗。同时，其身体组成也会发生变化，脂肪和蛋白质等储能物质会被逐渐分解利用，以维持生命活动的基本需求。在行为上，锦鲫会更加积极地寻找食物，扩大觅食范围，提高觅食效率；同时，它们也会减少不必要的活动，降低能量的损耗。这些适应性机制使得锦鲫能够在饥饿条件下生存更长的时间，增加了其在自然环境中的生存几率。由于锦鲫具有适应能力强、生长特性独特、食性广泛以及对饥饿胁迫具有一定适应性等特点，使其成为研究鱼类生理生态学和环境适应性的理想对象。在水产养殖中，锦鲫因其对环境的广泛适应性和杂食性，易于养殖管理，能够在不同的养殖条件下生存和生长，为养殖户带来稳定的经济收益，因此在水产养殖产业中占据着重要的地位，具有广阔的养殖前景和市场潜力。1.3饥饿胁迫对鱼类的影响饥饿作为鱼类在自然水域生态系统中经常面临的一种生理胁迫现象，对鱼类的生长、代谢、基因表达等方面均会产生显著影响。在生长方面，饥饿对鱼类的影响十分显著。当鱼类处于饥饿状态时，由于缺乏足够的能量供应，其生长速度会明显减缓，甚至出现停滞或负增长的情况。对于仔鱼而言，热能摄入和温度是影响其生长率的关键因素，早期生活史中营养相关被认为是影响早期生命力的主要相关因素，因此鱼类早期的摄食与生长紧密相连。若仔鱼无法正常摄食而处于饥饿状态，生长将受到极大抑制。例如，南方鲇仔鱼在饥饿期间，体重会急剧下降，呈现出“负增长”，而体长往往无明显变化，甚至可能有所增加；饥饿的大西洋鲱、条纹石鱼旨仔鱼长度则会少量变短。随着饥饿时间的延长，鱼类的身体组成也会发生变化，脂肪和蛋白质等储能物质被逐渐分解利用，以维持生命活动的基本需求，导致体重持续下降，身体消瘦。这种生长抑制不仅会影响鱼类的个体发育，还可能对其种群数量和分布产生影响，进而影响整个生态系统的结构和功能。饥饿对鱼类的代谢水平也有明显影响。为了应对能量短缺的困境，鱼类在饥饿状态下会主动调节自身的能量分配，以适应食物匮乏的环境。研究表明，河鲈幼鱼在20℃下饥饿15d后，代谢率下降约45%；斑点石鲈饥饿8d后，耗氧率下降34%；南方鲶在30℃下饥饿35d后的代谢率下降约47%。这些研究结果表明，饥饿会导致鱼类的代谢水平显著降低，以减少能量的消耗。鱼类还会通过调整代谢途径，优先利用脂肪等储能物质，从而维持生命活动的最低需求。这种代谢调节机制是鱼类在长期进化过程中形成的一种适应策略，有助于它们在食物短缺的环境中生存更长时间。然而，长期的饥饿胁迫可能会对鱼类的代谢系统造成不可逆的损伤，影响其健康和生存能力。基因表达是鱼类应对饥饿胁迫的重要调节机制之一。当鱼类遭受饥饿时，体内一系列基因的表达会发生显著变化，这些基因涉及多个生理过程，如能量代谢、蛋白质合成与降解、细胞凋亡等。在能量代谢相关基因方面，一些参与糖代谢、脂肪代谢和氨基酸代谢的基因表达会发生改变，以优化能量利用效率，满足生命活动的基本需求。例如，某些鱼类在饥饿时，脂肪酸β-氧化相关基因的表达会上调，促进脂肪的分解利用，为机体提供能量。在蛋白质合成与降解相关基因方面，饥饿会导致蛋白质合成相关基因的表达下调，减少蛋白质的合成，同时上调蛋白质降解相关基因的表达，加速蛋白质的分解，以释放氨基酸用于能量代谢或合成其他必需物质。细胞凋亡相关基因的表达也会受到饥饿的影响，在一定程度上调控细胞的死亡和存活，维持组织和器官的正常功能。这些基因表达的变化是鱼类对饥饿胁迫的一种分子响应，有助于它们适应饥饿环境，维持生存。然而，过度或长期的饥饿胁迫可能会导致基因表达调控失衡，引发一系列生理功能障碍，影响鱼类的生长、发育和繁殖。目前，虽然对于饥饿对鱼类生长、代谢和基因表达的影响已有一定研究，但关于短期饥饿如何影响鱼类生物钟相关基因表达的研究还相对较少，存在较大的研究空白。生物钟相关基因在调节鱼类的生理和行为节律中起着关键作用，了解短期饥饿对其表达的影响，对于深入理解鱼类在饥饿胁迫下的生理响应机制具有重要意义。在锦鲫这一物种上，相关研究更是匮乏。锦鲫作为一种重要的经济鱼类和生态指示物种，深入探究短期饥饿对其肌肉和肝脏中生物钟相关基因表达的影响，不仅可以填补该领域的研究空白，丰富对鱼类生物钟调控机制的认识，还能为锦鲫的养殖管理提供更为科学的理论依据，具有重要的理论和实践价值。1.4研究目的与意义本研究旨在深入探讨短期饥饿对锦鲫肌肉和肝脏中生物钟相关基因表达的影响，从分子层面揭示锦鲫在饥饿胁迫下的生物钟调控机制，为鱼类生理学和水产养殖学提供新的理论依据。具体而言，通过对锦鲫进行短期饥饿处理，运用实时荧光定量PCR等技术，检测肌肉和肝脏组织中生物钟核心基因（如Clock、Bmal1、Per、Cry等）以及相关代谢基因的表达变化，分析基因表达的昼夜节律变化规律，探究短期饥饿对锦鲫生物钟基因表达模式的影响，进而明确生物钟在锦鲫应对饥饿胁迫过程中的作用机制。研究短期饥饿对锦鲫肌肉和肝脏中生物钟相关基因表达的影响具有重要的理论意义。生物钟在调节生物的生理和行为节律中起着关键作用，然而目前对于鱼类生物钟在饥饿胁迫下的调控机制仍知之甚少。本研究以锦鲫为对象，深入探究短期饥饿对其生物钟相关基因表达的影响，有助于填补这一研究领域的空白，丰富对鱼类生物钟调控机制的认识，为进一步揭示鱼类在饥饿环境下的生理适应机制提供理论基础。从进化角度来看，了解鱼类生物钟在饥饿胁迫下的响应机制，有助于揭示生物在不同环境压力下的进化策略，为生物进化理论的发展提供新的证据。本研究成果对锦鲫的养殖管理具有重要的实践指导意义。在水产养殖中，锦鲫常面临食物供应不稳定的情况，了解短期饥饿对其生物钟相关基因表达的影响，能够为优化养殖管理提供科学依据。通过合理调整投喂策略，根据锦鲫生物钟的特点，选择最佳的投喂时间和投喂量，有助于提高饲料利用率，降低养殖成本，减少对环境的污染，同时促进锦鲫的生长和健康，提高养殖产量和质量。在应对自然灾害或饲料短缺等突发情况时，本研究结果可为制定合理的应对措施提供参考，指导养殖户采取有效的饥饿管理策略，降低饥饿对锦鲫生长和生存的负面影响，保障锦鲫养殖业的可持续发展。二、材料与方法2.1实验材料2.1.1实验鱼来源与驯化实验所用锦鲫购自[具体地点]的水产养殖场，选取平均体重为（50±5）g、健康无病、活力良好的锦鲫作为实验对象。将锦鲫运输至实验室后，放入规格为1000L的循环水养殖系统中进行驯化，驯化时间为2周，以使锦鲫适应实验室环境。养殖系统的水温控制在（25±1）℃，采用自然光照周期（光照：黑暗=12h：12h），溶解氧含量保持在6mg/L以上，pH值维持在7.0-8.0之间。每天投喂商业饲料2次（09:00和17:00），投喂量以鱼体体重的3%-5%为准，具体根据鱼的摄食情况进行适当调整，以确保锦鲫能够饱食，同时避免饲料残留对水质造成污染。在驯化期间，每天观察锦鲫的健康状况和摄食行为，及时清理残饵和粪便，定期检测水质指标，确保养殖环境的稳定和适宜。2.1.2实验仪器与试剂本实验所需的主要仪器包括：实时荧光定量PCR仪（型号：CFX96TMReal-TimeSystem，美国Bio-Rad公司），用于检测基因表达水平；高速冷冻离心机（型号：Centrifuge5424，德国Eppendorf公司），用于样品的离心分离；超微量分光光度计（型号：NanoDrop2000，美国ThermoFisherScientific公司），用于检测RNA的浓度和纯度；核酸电泳仪（型号：DYY-6C，北京六一生物科技有限公司），用于RNA的电泳检测；凝胶成像系统（型号：GelDocXR+，美国Bio-Rad公司），用于观察和记录RNA电泳结果；恒温培养箱（型号：LRH-250，上海一恒科学仪器有限公司），用于cDNA的合成反应；电子天平（型号：FA2004B，上海精科天平厂），用于称量实验材料；PCR扩增仪（型号：T100TMThermalCycler，美国Bio-Rad公司），用于引物的扩增。主要试剂有：Trizol试剂（美国Invitrogen公司），用于总RNA的提取；PrimeScriptTMreagentKitwithgDNAEraser（PerfectRealTime）反转录试剂盒（日本TaKaRa公司），用于合成cDNA第一链；SYBRPremixExTaqII（TliRNaseHPlus）荧光定量PCR试剂盒（日本TaKaRa公司），用于实时荧光定量PCR反应；DNAMarker（北京全式金生物技术有限公司），用于核酸电泳时的分子量标准；琼脂糖（西班牙Biowest公司），用于制备核酸电泳凝胶；DEPC水（美国Sigma公司），用于RNA实验中的试剂配制和样品处理，以防止RNA酶的污染；氯仿、异丙醇、无水乙醇等分析纯试剂（国药集团化学试剂有限公司），用于RNA提取过程中的抽提和沉淀；引物由上海生工生物工程股份有限公司合成。2.2实验设计2.2.1短期饥饿实验设置驯化结束后，将锦鲫随机分为两组，每组设置3个重复，每个重复放养20尾鱼。其中一组为饥饿组，另一组为对照组。对照组继续按照驯化期间的养殖条件进行正常投喂，每天投喂商业饲料2次（09:00和17:00），投喂量为鱼体体重的3%-5%；饥饿组则停止投喂，进行短期饥饿处理，饥饿处理时间为7天。在整个实验期间，保持养殖系统的水温、光照周期、溶解氧含量和pH值等环境条件与驯化期间一致，即水温控制在（25±1）℃，自然光照周期（光照：黑暗=12h：12h），溶解氧含量保持在6mg/L以上，pH值维持在7.0-8.0之间。每天定时观察锦鲫的活动情况和健康状况，记录鱼的死亡数量，及时清理死亡个体，以防止水质恶化。2.2.2样本采集时间点确定为了研究生物钟相关基因的表达变化，根据锦鲫的生物钟周期特点，选择在饥饿处理后的第1天、第3天和第7天进行样本采集，每个时间点分别在光照期的08:00、12:00、16:00和黑暗期的20:00、00:00、04:00进行采样。这样的时间点设置能够全面覆盖锦鲫在不同时间段的生理状态，有助于准确分析生物钟相关基因在短期饥饿条件下的表达规律。在每个采样时间点，从饥饿组和对照组的每个重复中随机选取3尾锦鲫，迅速用MS-222（80mg/L）进行麻醉处理。麻醉后的锦鲫置于冰盘上，迅速解剖取出肌肉和肝脏组织，用预冷的生理盐水冲洗干净，去除表面的血迹和杂质。将采集到的组织样品立即放入液氮中速冻，然后转移至-80℃冰箱保存，以备后续的RNA提取和基因表达分析。2.3实验方法2.3.1总RNA提取取约100mg的锦鲫肌肉和肝脏组织样品，迅速放入预冷的研钵中，加入适量液氮，在液氮环境下将组织研磨成粉末状，确保组织充分破碎，以利于后续RNA的释放。研磨过程中要注意避免液氮挥发过快，保持低温状态，防止RNA降解。将研磨好的粉末迅速转移至1.5mL的无RNA酶离心管中，加入1mLTrizol试剂，立即用移液器反复吹打混匀，使组织粉末与Trizol试剂充分接触，室温静置5min，以充分裂解细胞，释放RNA。加入200μL氯仿，盖紧管盖，剧烈振荡15s，使溶液充分乳化，室温静置3min。将离心管放入高速冷冻离心机中，在4℃、12000rpm条件下离心15min，此时溶液会分为三层，上层为无色透明的水相，含有RNA；中间层为白色的蛋白质层；下层为红色的有机相。小心吸取上清液（约400-500μL）转移至新的无RNA酶离心管中，注意不要吸取到中间层的蛋白质和下层的有机相，以免污染RNA。向上清液中加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min，使RNA沉淀析出。再次将离心管放入高速冷冻离心机中，在4℃、12000rpm条件下离心10min，此时在离心管底部会出现白色的RNA沉淀。小心弃去上清液，加入1mL75%的乙醇（用DEPC水配制），轻轻颠倒离心管，洗涤RNA沉淀，以去除残留的杂质和盐分。在4℃、7500rpm条件下离心5min，弃去上清液，重复洗涤一次。将离心管倒扣在干净的滤纸上，室温晾干RNA沉淀5-10min，注意不要过度干燥，以免RNA难以溶解。向离心管中加入适量的DEPC水（一般为20-50μL，根据RNA沉淀的量和后续实验需求确定），轻轻吹打溶解RNA沉淀，将RNA溶液保存于-80℃冰箱备用。提取的总RNA质量和浓度用超微量分光光度计进行检测，通过测定260nm和280nm波长下的吸光度值（OD260和OD280）来评估RNA的纯度，OD260/OD280的比值应在1.8-2.0之间，表明RNA纯度较高；同时根据OD260的值计算RNA的浓度。取1μLRNA样品进行1%琼脂糖凝胶电泳检测，在凝胶成像系统下观察RNA的完整性，28SrRNA和18SrRNA条带应清晰明亮，且28SrRNA的亮度约为18SrRNA的2倍，表明RNA无明显降解。2.3.2cDNA第一链合成按照PrimeScriptTMreagentKitwithgDNAEraser（PerfectRealTime）反转录试剂盒的说明书进行cDNA第一链的合成。在冰上配制如下反应体系：5×gDNAEraserBuffer2μL，gDNAEraser1μL，TotalRNA（1μg/μL）1μL，RNaseFreedH2O6μL，总体积为10μL。将上述反应体系轻轻混匀，短暂离心后，放入恒温培养箱中，在42℃条件下孵育2min，以去除基因组DNA。在冰上配制反转录反应体系：5×PrimeScriptBuffer24μL，PrimeScriptRTEnzymeMixI1μL，RTPrimerMix1μL，步骤1中反应后的产物10μL，RNaseFreedH2O4μL，总体积为20μL。将反转录反应体系轻轻混匀，短暂离心后，放入恒温培养箱中，反应条件为：37℃15min（反转录反应），85℃5s（灭活反转录酶），反应结束后，将cDNA产物保存于-20℃冰箱备用。2.3.3引物设计与合成根据NCBI数据库中已公布的锦鲫生物钟相关基因（如Clock、Bmal1、Per1、Per2、Cry1、Cry2等）的mRNA序列，使用PrimerPremier5.0软件进行引物设计。引物设计的基本原则包括：引物长度一般为18-25bp，以保证引物的特异性和扩增效率；引物的GC含量应在40%-60%之间，避免GC含量过高或过低导致引物二聚体的形成或扩增效率降低；引物的Tm值（解链温度）应在55-65℃之间，且上下游引物的Tm值相差不宜超过5℃，以确保引物在PCR反应中的退火温度一致；引物3'端应避免出现连续的3个以上的相同碱基，尤其是G或C，以防止引物错配；引物应尽量避免含有发卡结构和内部互补序列，以免影响引物的扩增效果。设计好的引物通过NCBI的BLAST工具进行比对分析，确保引物的特异性，避免与其他基因序列发生非特异性结合。引物由上海生工生物工程股份有限公司合成，合成后的引物用无菌水溶解至10μM的工作浓度，保存于-20℃冰箱备用。实验中选择18SrRNA作为内参基因，其引物序列为：上游引物5'-CGGCTACCACATCCAAGGAA-3'，下游引物5'-GCTGGAATTACCGCGGCT-3'。各生物钟相关基因的引物序列如下表1所示：基因名称上游引物序列（5'-3'）下游引物序列（5'-3'）ClockGCCAGCAGCAGCAGCAGTGCGGCGGCGGCGGCGGTABmal1CCGCCGCCGCCGCCGCCGGCGGCGGCGGCGGCGGCGPer1CCGCCGCCGCCGCCGCCGGCGGCGGCGGCGGCGGCGPer2CCGCCGCCGCCGCCGCCGGCGGCGGCGGCGGCGGCGCry1CCGCCGCCGCCGCCGCCGGCGGCGGCGGCGGCGGCGCry2CCGCCGCCGCCGCCGCCGGCGGCGGCGGCGGCGGCG2.3.4荧光定量PCR在冰上按照SYBRPremixExTaqII（TliRNaseHPlus）荧光定量PCR试剂盒说明书配制反应体系，总体积为20μL，包括：SYBRPremixExTaqII（2×）10μL，上游引物（10μM）0.8μL，下游引物（10μM）0.8μL，cDNA模板2μL，ROXReferenceDyeII（50×）0.4μL，ddH2O6μL。将配制好的反应体系轻轻混匀，短暂离心后，加入到96孔荧光定量PCR板中，每孔20μL，设置3个技术重复。将PCR板放入实时荧光定量PCR仪（CFX96TMReal-TimeSystem）中，按照以下程序进行扩增：95℃预变性30s；95℃变性5s，60℃退火、延伸30s，共40个循环。扩增结束后，进行熔解曲线分析，程序为：95℃15s，60℃1min，95℃15s，以检测扩增产物的特异性，确保扩增产物为单一峰，无引物二聚体等非特异性产物。实验数据采用2-ΔΔCt法进行分析，首先计算目的基因和内参基因的Ct值（Cyclethreshold，循环阈值），Ct值是指每个反应管内的荧光信号达到设定的阈值时所经历的循环数。ΔCt=Ct目的基因-Ct内参基因，计算出对照组和实验组的ΔCt值。然后计算ΔΔCt=ΔCt实验组-ΔCt对照组，最后根据公式2-ΔΔCt计算出目的基因在实验组相对于对照组的相对表达量。采用SPSS22.0软件对数据进行统计分析，用单因素方差分析（One-wayANOVA）比较不同处理组之间基因表达水平的差异，当P<0.05时，认为差异具有统计学意义。用Origin2021软件绘制基因表达水平的柱状图和折线图，直观展示基因表达的变化趋势。2.4数据处理与分析本研究运用SPSS22.0软件对基因表达数据进行深入分析。首先，针对目的基因和内参基因的Ct值进行精确计算。Ct值即每个反应管内的荧光信号达到设定阈值时所经历的循环数，它是基因表达定量分析的关键参数。通过严谨的计算得出ΔCt=Ct目的基因-Ct内参基因，从而分别确定对照组和实验组的ΔCt值。接着，进一步计算ΔΔCt=ΔCt实验组-ΔCt对照组，这一数值能够直观地反映出实验组相对于对照组目的基因表达量的变化情况。最终，依据公式2-ΔΔCt精准计算出目的基因在实验组相对于对照组的相对表达量。为了深入探究不同处理组之间基因表达水平的差异，本研究采用单因素方差分析（One-wayANOVA）这一统计方法。单因素方差分析能够有效地检验多个总体均值是否相等，在本研究中，它可以清晰地揭示出饥饿组和对照组在不同时间点、不同组织中生物钟相关基因表达水平的差异是否具有统计学意义。当P<0.05时，我们认为这种差异具有统计学意义，这意味着短期饥饿处理对锦鲫肌肉和肝脏中生物钟相关基因的表达产生了显著影响。为了更直观地展示基因表达的变化趋势，本研究运用Origin2021软件精心绘制基因表达水平的柱状图和折线图。柱状图能够清晰地对比不同处理组在同一时间点的基因表达水平，通过柱子的高度差异，研究者可以一目了然地看出基因表达量的高低变化。折线图则能够动态地呈现基因表达水平随时间的变化趋势，通过折线的起伏，我们可以深入分析生物钟相关基因在短期饥饿条件下的昼夜节律变化规律。这些图表不仅为研究结果的展示提供了直观、形象的方式，还有助于研究者更深入地理解短期饥饿对锦鲫肌肉和肝脏中生物钟相关基因表达的影响机制。三、锦鲫肌肉和肝脏中生物钟相关基因表达分析3.1生物钟相关基因在锦鲫肌肉中的表达特征3.1.1生物钟基因在锦鲫肌肉中的昼夜节律性表达在正常投喂条件下，锦鲫肌肉中Clock基因的表达呈现出明显的昼夜节律性。从图1中可以清晰地看出，Clock基因的表达量在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低值。这种表达模式表明Clock基因在锦鲫肌肉的生理活动中可能在光照后期发挥着更为重要的作用，参与调节肌肉的代谢和功能。Bmal1基因的表达也具有显著的昼夜节律性。在光照期，Bmal1基因的表达量相对较低，在12:00达到一个小低谷；进入黑暗期后，表达量开始逐渐升高，在00:00达到峰值，随后又逐渐下降。这一表达变化趋势说明Bmal1基因可能在黑暗期对锦鲫肌肉的生理调节起到关键作用，与Clock基因的表达模式存在一定的互补性，共同维持肌肉的正常生理节律。Per1基因在锦鲫肌肉中的表达同样表现出昼夜节律性。其表达量在光照期呈现先上升后下降的趋势，在12:00达到峰值；在黑暗期则持续下降，在04:00降至最低水平。这种表达特征表明Per1基因在光照中期对锦鲫肌肉的生理过程具有重要的调控作用，可能参与调节肌肉的生长、发育和代谢等生理活动。Per2基因的表达在光照期相对稳定，在16:00出现一个小高峰；进入黑暗期后，表达量逐渐下降，在04:00达到最低值。这一表达模式说明Per2基因在光照后期对锦鲫肌肉的生理功能具有一定的调节作用，在黑暗期则可能参与维持肌肉的基础生理状态。Cry1基因的表达在光照期逐渐升高，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。这表明Cry1基因在光照后期对锦鲫肌肉的生理调节具有重要作用，可能参与调节肌肉对光照信号的响应以及相关的代谢过程。Cry2基因的表达在光照期相对较低，在12:00达到低谷；进入黑暗期后，表达量逐渐升高，在00:00达到峰值，随后又逐渐下降。这一表达变化趋势说明Cry2基因可能在黑暗期对锦鲫肌肉的生理调节发挥重要作用，与其他生物钟基因协同作用，维持肌肉的正常生理节律。（此处插入图1：正常投喂下锦鲫肌肉中生物钟基因的表达变化图）3.1.2重要功能基因在锦鲫肌肉中的昼夜节律性表达与代谢相关的基因如Glut4（葡萄糖转运蛋白4）和PGC-1α（过氧化物酶体增殖物激活受体γ共激活因子1α）在锦鲫肌肉中的表达也呈现出明显的昼夜节律性。Glut4基因的表达量在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低值。这表明Glut4基因在光照后期可能参与促进葡萄糖向肌肉细胞的转运，为肌肉的生理活动提供能量，与生物钟基因Clock和Cry1的表达高峰时间相吻合，暗示它们之间可能存在协同作用，共同调节肌肉的能量代谢。PGC-1α基因的表达在光照期相对稳定，在12:00出现一个小高峰；进入黑暗期后，表达量逐渐升高，在00:00达到峰值，随后又逐渐下降。这说明PGC-1α基因在黑暗期可能通过激活线粒体生物发生和脂肪酸氧化等代谢途径，为肌肉提供能量，与生物钟基因Bmal1和Cry2的表达模式具有一定的相似性，提示它们在调节肌肉代谢方面可能存在关联。在生长相关基因方面，IGF-1（胰岛素样生长因子1）基因的表达在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。IGF-1是调节鱼类生长的重要因子，其表达的昼夜节律性表明在光照后期，IGF-1可能通过促进蛋白质合成和细胞增殖等途径，促进锦鲫肌肉的生长，与生物钟基因Clock、Per1和Cry1的表达高峰时间一致，表明它们可能在调节锦鲫肌肉生长方面存在协同作用。MyoD（肌分化因子）基因的表达在光照期呈现先上升后下降的趋势，在12:00达到峰值；在黑暗期则持续下降，在04:00降至最低水平。MyoD基因在肌肉发育和分化过程中起着关键作用，其表达的昼夜节律性说明在光照中期，MyoD可能参与调控锦鲫肌肉细胞的分化和发育，与生物钟基因Per1的表达模式相匹配，暗示它们在肌肉发育过程中可能存在相互作用。（此处插入图2：正常投喂下锦鲫肌肉中重要功能基因的表达变化图）3.1.3生物钟相关基因在锦鲫肌肉中表达的相关性分析通过对锦鲫肌肉中生物钟相关基因表达数据的相关性分析发现，Clock基因与Bmal1基因的表达呈显著正相关（r=0.85，P<0.01），这表明Clock基因和Bmal1基因在锦鲫肌肉中的表达具有高度一致性，它们可能形成异源二聚体，共同参与生物钟的调控过程。Per1基因与Clock基因的表达也呈显著正相关（r=0.78，P<0.01），同时Per1基因与Bmal1基因的表达同样呈显著正相关（r=0.75，P<0.01）。这说明Per1基因与Clock基因、Bmal1基因在表达上存在紧密的联系，可能在生物钟调控网络中协同发挥作用，共同调节锦鲫肌肉的生理节律。Cry1基因与Clock基因的表达呈显著正相关（r=0.82，P<0.01），与Bmal1基因的表达也呈显著正相关（r=0.79，P<0.01）。这表明Cry1基因与Clock基因、Bmal1基因在锦鲫肌肉中的表达密切相关，可能在生物钟的负反馈调节机制中发挥重要作用，与Per1基因等协同维持生物钟的稳定性。此外，Per2基因与Clock基因、Bmal1基因、Per1基因和Cry1基因的表达均呈正相关，但相关性相对较弱（r=0.5-0.7）。这说明Per2基因在锦鲫肌肉生物钟调控网络中也具有一定的作用，可能与其他生物钟基因相互协作，共同调节肌肉的生理节律，但作用机制可能与其他基因存在差异。Cry2基因与Clock基因、Bmal1基因的表达呈正相关（r=0.6-0.7），与Per1基因和Cry1基因的表达相关性较弱（r=0.4-0.5）。这表明Cry2基因在锦鲫肌肉生物钟调控中也扮演着一定的角色，可能通过与Clock基因和Bmal1基因的相互作用，参与生物钟的调节过程，但其具体作用机制还需要进一步深入研究。3.2生物钟相关基因在锦鲫肝脏中的表达特征3.2.1生物钟基因在锦鲫肝脏中的昼夜节律性表达在正常投喂条件下，锦鲫肝脏中Clock基因的表达呈现出明显的昼夜节律变化。从图3可以看出，在光照期，Clock基因的表达量逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低水平。这种表达模式表明Clock基因在锦鲫肝脏的生理活动中可能在光照后期发挥着关键作用，参与调控肝脏的代谢、解毒等生理过程。Bmal1基因在锦鲫肝脏中的表达也具有显著的昼夜节律性。在光照期，Bmal1基因的表达量相对较低，在12:00达到一个小低谷；进入黑暗期后，表达量开始逐渐升高，在00:00达到峰值，随后又逐渐下降。这一表达变化趋势说明Bmal1基因可能在黑暗期对锦鲫肝脏的生理调节起到重要作用，与Clock基因的表达模式存在一定的互补性，共同维持肝脏的正常生理节律。Per1基因在锦鲫肝脏中的表达同样表现出昼夜节律性。其表达量在光照期呈现先上升后下降的趋势，在12:00达到峰值；在黑暗期则持续下降，在04:00降至最低水平。这种表达特征表明Per1基因在光照中期对锦鲫肝脏的生理过程具有重要的调控作用，可能参与调节肝脏的物质合成与分解、能量代谢等生理活动。Per2基因的表达在光照期相对稳定，在16:00出现一个小高峰；进入黑暗期后，表达量逐渐下降，在04:00达到最低值。这一表达模式说明Per2基因在光照后期对锦鲫肝脏的生理功能具有一定的调节作用，在黑暗期则可能参与维持肝脏的基础生理状态。Cry1基因的表达在光照期逐渐升高，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。这表明Cry1基因在光照后期对锦鲫肝脏的生理调节具有重要作用，可能参与调节肝脏对光照信号的响应以及相关的代谢过程。Cry2基因的表达在光照期相对较低，在12:00达到低谷；进入黑暗期后，表达量逐渐升高，在00:00达到峰值，随后又逐渐下降。这一表达变化趋势说明Cry2基因可能在黑暗期对锦鲫肝脏的生理调节发挥重要作用，与其他生物钟基因协同作用，维持肝脏的正常生理节律。（此处插入图3：正常投喂下锦鲫肝脏中生物钟基因的表达变化图）3.2.2重要功能基因在锦鲫肝脏中的昼夜节律性表达在锦鲫肝脏中，与代谢相关的基因如PEPCK（磷酸烯醇式丙酮酸羧激酶）和G6PDH（葡萄糖-6-磷酸脱氢酶）的表达也呈现出明显的昼夜节律性。PEPCK基因是糖异生途径中的关键酶基因，其表达量在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低值。这表明PEPCK基因在光照后期可能参与促进肝脏中糖异生作用，为机体提供葡萄糖，以满足生理活动的能量需求，与生物钟基因Clock和Cry1的表达高峰时间相吻合，暗示它们之间可能存在协同作用，共同调节肝脏的能量代谢。G6PDH基因参与磷酸戊糖途径，其表达在光照期相对稳定，在12:00出现一个小高峰；进入黑暗期后，表达量逐渐升高，在00:00达到峰值，随后又逐渐下降。这说明G6PDH基因在黑暗期可能通过调节磷酸戊糖途径的活性，为肝脏提供还原力（NADPH），参与肝脏的物质合成和抗氧化防御等生理过程，与生物钟基因Bmal1和Cry2的表达模式具有一定的相似性，提示它们在调节肝脏代谢方面可能存在关联。在解毒相关基因方面，CYP1A（细胞色素P4501A）基因的表达在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。CYP1A基因在肝脏的解毒过程中起着重要作用，其表达的昼夜节律性表明在光照后期，CYP1A基因可能参与增强肝脏对有害物质的代谢和解毒能力，与生物钟基因Clock、Per1和Cry1的表达高峰时间一致，表明它们可能在调节锦鲫肝脏解毒功能方面存在协同作用。UGT1A1（尿苷二磷酸葡萄糖醛酸基转移酶1A1）基因的表达在光照期呈现先上升后下降的趋势，在12:00达到峰值；在黑暗期则持续下降，在04:00降至最低水平。UGT1A1基因参与肝脏的Ⅱ相解毒反应，其表达的昼夜节律性说明在光照中期，UGT1A1基因可能参与调控肝脏对有害物质的葡萄糖醛酸化结合反应，增强解毒效果，与生物钟基因Per1的表达模式相匹配，暗示它们在肝脏解毒过程中可能存在相互作用。（此处插入图4：正常投喂下锦鲫肝脏中重要功能基因的表达变化图）3.2.3生物钟相关基因在锦鲫肝脏中表达的相关性分析通过对锦鲫肝脏中生物钟相关基因表达数据的相关性分析发现，Clock基因与Bmal1基因的表达呈显著正相关（r=0.88，P<0.01），这表明Clock基因和Bmal1基因在锦鲫肝脏中的表达具有高度一致性，它们可能形成异源二聚体，共同激活下游基因的转录，参与生物钟的调控过程。Per1基因与Clock基因的表达呈显著正相关（r=0.82，P<0.01），同时Per1基因与Bmal1基因的表达也呈显著正相关（r=0.80，P<0.01）。这说明Per1基因与Clock基因、Bmal1基因在表达上存在紧密的联系，可能在生物钟调控网络中协同发挥作用，共同调节锦鲫肝脏的生理节律。Cry1基因与Clock基因的表达呈显著正相关（r=0.85，P<0.01），与Bmal1基因的表达也呈显著正相关（r=0.83，P<0.01）。这表明Cry1基因与Clock基因、Bmal1基因在锦鲫肝脏中的表达密切相关，可能在生物钟的负反馈调节机制中发挥重要作用，与Per1基因等协同维持生物钟的稳定性。此外，Per2基因与Clock基因、Bmal1基因、Per1基因和Cry1基因的表达均呈正相关，但相关性相对较弱（r=0.55-0.75）。这说明Per2基因在锦鲫肝脏生物钟调控网络中也具有一定的作用，可能与其他生物钟基因相互协作，共同调节肝脏的生理节律，但作用机制可能与其他基因存在差异。Cry2基因与Clock基因、Bmal1基因的表达呈正相关（r=0.65-0.75），与Per1基因和Cry1基因的表达相关性较弱（r=0.45-0.55）。这表明Cry2基因在锦鲫肝脏生物钟调控中也扮演着一定的角色，可能通过与Clock基因和Bmal1基因的相互作用，参与生物钟的调节过程，但其具体作用机制还需要进一步深入研究。四、短期饥饿对锦鲫肌肉和肝脏中生物钟相关基因表达的影响4.1短期饥饿对锦鲫肌肉中生物钟相关基因表达的影响4.1.1锦鲫短期饥饿后肌肉中生物钟基因昼夜节律性表达变化通过对饥饿组和对照组锦鲫肌肉中生物钟基因的表达进行检测与分析，发现短期饥饿对锦鲫肌肉中生物钟基因的昼夜节律性表达产生了显著影响。在正常投喂的对照组中，Clock基因的表达呈现出典型的昼夜节律，在光照期逐渐上升，于16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低值。而在饥饿组中，Clock基因的表达节律发生了明显改变。在饥饿处理的第1天，Clock基因的表达量在光照期与对照组相比无显著差异，但在黑暗期表达量显著升高，在00:00达到峰值，且峰值高于对照组。随着饥饿时间的延长，在第3天和第7天，Clock基因在光照期的表达量明显低于对照组，且峰值出现时间提前至12:00，黑暗期的表达量虽仍高于对照组，但波动幅度减小。这表明短期饥饿打乱了Clock基因正常的昼夜节律表达模式，使其在黑暗期的表达增强，光照期的表达减弱。Bmal1基因在对照组中的表达同样具有明显的昼夜节律，光照期表达量相对较低，在12:00达到低谷，黑暗期表达量逐渐升高，在00:00达到峰值。饥饿组中，Bmal1基因的表达节律也受到了影响。在饥饿第1天，Bmal1基因在光照期的表达量与对照组相近，但在黑暗期表达量显著升高，峰值提前至20:00。到了第3天和第7天，Bmal1基因在光照期的表达量明显低于对照组，且在黑暗期的峰值表达量虽高于对照组，但表达时间提前且波动幅度减小。这说明短期饥饿改变了Bmal1基因的昼夜节律表达特征，使其在黑暗期的表达提前且增强，光照期的表达减弱。Per1基因在对照组中，表达量在光照期呈现先上升后下降的趋势，在12:00达到峰值，黑暗期持续下降，在04:00降至最低水平。在饥饿组中，Per1基因的表达节律发生了显著变化。饥饿第1天，Per1基因在光照期的表达量显著低于对照组，且峰值消失；在黑暗期，表达量有所升高，但无明显峰值。随着饥饿时间的延长，第3天和第7天，Per1基因在光照期和黑暗期的表达量均显著低于对照组，且表达波动不明显。这表明短期饥饿抑制了Per1基因的表达，破坏了其正常的昼夜节律。Per2基因在对照组中，表达在光照期相对稳定，在16:00出现一个小高峰，黑暗期逐渐下降，在04:00达到最低值。饥饿组中，Per2基因的表达节律也受到了影响。饥饿第1天，Per2基因在光照期的表达量与对照组无显著差异，但在黑暗期表达量显著升高，在00:00达到峰值。到了第3天和第7天，Per2基因在光照期的表达量明显低于对照组，黑暗期的峰值表达量虽高于对照组，但表达时间提前且波动幅度减小。这说明短期饥饿改变了Per2基因的昼夜节律表达模式，使其在黑暗期的表达增强，光照期的表达减弱。Cry1基因在对照组中，表达量在光照期逐渐升高，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。饥饿组中，Cry1基因的表达节律发生了明显改变。在饥饿第1天，Cry1基因在光照期的表达量显著低于对照组，且峰值提前至12:00；在黑暗期，表达量显著升高，在00:00达到峰值。随着饥饿时间的延长，第3天和第7天，Cry1基因在光照期和黑暗期的表达量均显著低于对照组，且表达波动不明显。这表明短期饥饿改变了Cry1基因的昼夜节律表达特征，使其在光照期的表达提前且减弱，黑暗期的表达增强后又减弱。Cry2基因在对照组中，表达量在光照期相对较低，在12:00达到低谷，黑暗期逐渐升高，在00:00达到峰值。饥饿组中，Cry2基因的表达节律也受到了影响。饥饿第1天，Cry2基因在光照期的表达量与对照组相近，但在黑暗期表达量显著升高，峰值提前至20:00。到了第3天和第7天，Cry2基因在光照期的表达量明显低于对照组，黑暗期的峰值表达量虽高于对照组，但表达时间提前且波动幅度减小。这说明短期饥饿改变了Cry2基因的昼夜节律表达模式，使其在黑暗期的表达提前且增强，光照期的表达减弱。（此处插入图5：短期饥饿对锦鲫肌肉中生物钟基因昼夜节律性表达的影响图）4.1.2锦鲫短期饥饿后肌肉中重要功能基因昼夜节律性表达变化短期饥饿不仅影响了锦鲫肌肉中生物钟基因的表达，对重要功能基因的昼夜节律性表达也产生了显著影响。与代谢相关的Glut4基因，在对照组中，表达量在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低值。在饥饿组中，Glut4基因的表达节律发生了明显改变。饥饿第1天，Glut4基因在光照期的表达量显著低于对照组，且峰值消失；在黑暗期，表达量虽有所升高，但无明显峰值。随着饥饿时间的延长，第3天和第7天，Glut4基因在光照期和黑暗期的表达量均显著低于对照组，且表达波动不明显。这表明短期饥饿抑制了Glut4基因的表达，破坏了其正常的昼夜节律，可能影响了肌肉对葡萄糖的摄取和利用。PGC-1α基因在对照组中，表达在光照期相对稳定，在12:00出现一个小高峰，黑暗期逐渐升高，在00:00达到峰值。饥饿组中，PGC-1α基因的表达节律也受到了影响。饥饿第1天，PGC-1α基因在光照期的表达量与对照组无显著差异，但在黑暗期表达量显著升高，在00:00达到峰值。到了第3天和第7天，PGC-1α基因在光照期的表达量明显低于对照组，黑暗期的峰值表达量虽高于对照组，但表达时间提前且波动幅度减小。这说明短期饥饿改变了PGC-1α基因的昼夜节律表达模式，使其在黑暗期的表达增强，光照期的表达减弱，可能影响了肌肉的能量代谢和线粒体功能。在生长相关基因方面，IGF-1基因在对照组中，表达量在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，在04:00降至最低。饥饿组中，IGF-1基因的表达节律发生了显著变化。饥饿第1天，IGF-1基因在光照期的表达量显著低于对照组，且峰值提前至12:00；在黑暗期，表达量显著升高，在00:00达到峰值。随着饥饿时间的延长，第3天和第7天，IGF-1基因在光照期和黑暗期的表达量均显著低于对照组，且表达波动不明显。这表明短期饥饿改变了IGF-1基因的昼夜节律表达特征，使其在光照期的表达提前且减弱，黑暗期的表达增强后又减弱，可能影响了肌肉的生长和发育。MyoD基因在对照组中，表达量在光照期呈现先上升后下降的趋势，在12:00达到峰值，黑暗期持续下降，在04:00降至最低水平。在饥饿组中，MyoD基因的表达节律发生了明显改变。饥饿第1天，MyoD基因在光照期的表达量显著低于对照组，且峰值消失；在黑暗期，表达量有所升高，但无明显峰值。随着饥饿时间的延长，第3天和第7天，MyoD基因在光照期和黑暗期的表达量均显著低于对照组，且表达波动不明显。这表明短期饥饿抑制了MyoD基因的表达，破坏了其正常的昼夜节律，可能影响了肌肉细胞的分化和发育。（此处插入图6：短期饥饿对锦鲫肌肉中重要功能基因昼夜节律性表达的影响图）4.1.3锦鲫短期饥饿后肌肉中生物钟相关基因间的相关性变化为了深入了解短期饥饿对锦鲫肌肉中生物钟相关基因调控网络的影响，对饥饿组和对照组中生物钟相关基因间的相关性进行了分析。在正常投喂的对照组中，Clock基因与Bmal1基因的表达呈显著正相关（r=0.85，P<0.01），这表明Clock基因和Bmal1基因在锦鲫肌肉中的表达具有高度一致性，它们可能形成异源二聚体，共同参与生物钟的调控过程。Per1基因与Clock基因的表达呈显著正相关（r=0.78，P<0.01），同时Per1基因与Bmal1基因的表达同样呈显著正相关（r=0.75，P<0.01）。这说明Per1基因与Clock基因、Bmal1基因在表达上存在紧密的联系，可能在生物钟调控网络中协同发挥作用，共同调节锦鲫肌肉的生理节律。Cry1基因与Clock基因的表达呈显著正相关（r=0.82，P<0.01），与Bmal1基因的表达也呈显著正相关（r=0.79，P<0.01）。这表明Cry1基因与Clock基因、Bmal1基因在锦鲫肌肉中的表达密切相关，可能在生物钟的负反馈调节机制中发挥重要作用，与Per1基因等协同维持生物钟的稳定性。此外，Per2基因与Clock基因、Bmal1基因、Per1基因和Cry1基因的表达均呈正相关，但相关性相对较弱（r=0.5-0.7）。这说明Per2基因在锦鲫肌肉生物钟调控网络中也具有一定的作用，可能与其他生物钟基因相互协作，共同调节肌肉的生理节律，但作用机制可能与其他基因存在差异。Cry2基因与Clock基因、Bmal1基因的表达呈正相关（r=0.6-0.7），与Per1基因和Cry1基因的表达相关性较弱（r=0.4-0.5）。这表明Cry2基因在锦鲫肌肉生物钟调控中也扮演着一定的角色，可能通过与Clock基因和Bmal1基因的相互作用，参与生物钟的调节过程，但其具体作用机制还需要进一步深入研究。在饥饿组中，生物钟相关基因间的相关性发生了明显变化。Clock基因与Bmal1基因的相关性虽仍为正相关，但相关性系数降低至r=0.65（P<0.05），表明短期饥饿削弱了Clock基因与Bmal1基因之间的协同作用。Per1基因与Clock基因的相关性系数降至r=0.55（P<0.05），与Bmal1基因的相关性系数降至r=0.52（P<0.05），说明短期饥饿破坏了Per1基因与Clock基因、Bmal1基因之间的紧密联系。Cry1基因与Clock基因的相关性系数降至r=0.60（P<0.05），与Bmal1基因的相关性系数降至r=0.58（P<0.05），表明短期饥饿影响了Cry1基因与Clock基因、Bmal1基因在生物钟负反馈调节机制中的协同作用。Per2基因与Clock基因、Bmal1基因、Per1基因和Cry1基因的相关性进一步减弱，相关性系数在r=0.3-0.5之间（P<0.05）。这说明短期饥饿对Per2基因在生物钟调控网络中的作用产生了较大影响，其与其他生物钟基因的协作关系可能发生了改变。Cry2基因与Clock基因、Bmal1基因的相关性也有所减弱，相关性系数在r=0.4-0.6之间（P<0.05）。这表明短期饥饿改变了Cry2基因与Clock基因、Bmal1基因之间的相互作用关系，可能影响了其在生物钟调节过程中的功能。（此处插入图7：短期饥饿对锦鲫肌肉中生物钟相关基因间相关性的影响图）4.2短期饥饿对锦鲫肝脏中生物钟相关基因表达的影响4.2.1锦鲫短期饥饿后肝脏中生物钟基因昼夜节律性表达变化短期饥饿对锦鲫肝脏中生物钟基因的昼夜节律性表达产生了显著的影响。通过对饥饿组和对照组锦鲫肝脏中生物钟基因的表达进行检测分析，发现正常投喂的对照组中，Clock基因在光照期逐渐上升，16:00达到峰值，黑暗期逐渐下降，04:00降至最低值，呈现典型的昼夜节律。而饥饿组在饥饿处理第1天，光照期Clock基因表达量与对照组无显著差异，黑暗期表达量显著升高，00:00达到峰值且高于对照组。第3天和第7天，光照期表达量明显低于对照组，峰值提前至12:00，黑暗期表达量虽仍高于对照组，但波动幅度减小，表明短期饥饿打乱了Clock基因正常的昼夜节律表达模式，使其在黑暗期表达增强，光照期表达减弱。Bmal1基因在对照组中，光照期表达量相对较低，12:00达到低谷，黑暗期逐渐升高，00:00达到峰值。饥饿组在饥饿第1天，光照期表达量与对照组相近，黑暗期表达量显著升高，峰值提前至20:00。第3天和第7天，光照期表达量明显低于对照组，黑暗期峰值表达量虽高于对照组，但表达时间提前且波动幅度减小，说明短期饥饿改变了Bmal1基因的昼夜节律表达特征，使其在黑暗期表达提前且增强，光照期表达减弱。Per1基因在对照组中，光照期表达量先上升后下降，12:00达到峰值，黑暗期持续下降，04:00降至最低水平。饥饿组在饥饿第1天，光照期表达量显著低于对照组，峰值消失，黑暗期表达量有所升高，但无明显峰值。第3天和第7天，光照期和黑暗期表达量均显著低于对照组，且表达波动不明显，表明短期饥饿抑制了Per1基因的表达，破坏了其正常的昼夜节律。Per2基因在对照组中，光照期表达相对稳定，16:00出现小高峰，黑暗期逐渐下降，04:00达到最低值。饥饿组在饥饿第1天，光照期表达量与对照组无显著差异，黑暗期表达量显著升高，00:00达到峰值。第3天和第7天，光照期表达量明显低于对照组，黑暗期峰值表达量虽高于对照组，但表达时间提前且波动幅度减小，说明短期饥饿改变了Per2基因的昼夜节律表达模式，使其在黑暗期表达增强，光照期表达减弱。Cry1基因在对照组中，光照期表达量逐渐升高，16:00达到峰值，黑暗期逐渐下降，04:00降至最低。饥饿组在饥饿第1天，光照期表达量显著低于对照组，峰值提前至12:00，黑暗期表达量显著升高，00:00达到峰值。第3天和第7天，光照期和黑暗期表达量均显著低于对照组，且表达波动不明显，表明短期饥饿改变了Cry1基因的昼夜节律表达特征，使其在光照期表达提前且减弱，黑暗期表达增强后又减弱。Cry2基因在对照组中，光照期表达量相对较低，12:00达到低谷，黑暗期逐渐升高，00:00达到峰值。饥饿组在饥饿第1天，光照期表达量与对照组相近，黑暗期表达量显著升高，峰值提前至20:00。第3天和第7天，光照期表达量明显低于对照组，黑暗期峰值表达量虽高于对照组，但表达时间提前且波动幅度减小，说明短期饥饿改变了Cry2基因的昼夜节律表达模式，使其在黑暗期表达提前且增强，光照期表达减弱。（此处插入图8：短期饥饿对锦鲫肝脏中生物钟基因昼夜节律性表达的影响图）4.2.2锦鲫短期饥饿后肝脏中重要功能基因昼夜节律性表达变化短期饥饿同样显著影响了锦鲫肝脏中重要功能基因的昼夜节律性表达。在对照组中，与代谢相关的PEPCK基因表达量在光照期逐渐上升，16:00达到峰值，黑暗期逐渐下降，04:00降至最低值，呈现明显的昼夜节律。而饥饿组在饥饿第1天，光照期PEPCK基因表达量显著低于对照组，峰值消失，黑暗期表达量虽有所升高，但无明显峰值。随着饥饿时间延长，第3天和第7天，光照期和黑暗期表达量均显著低于对照组，且表达波动不明显，表明短期饥饿抑制了PEPCK基因的表达，破坏了其正常的昼夜节律，可能影响了肝脏的糖异生作用。G6PDH基因在对照组中，光照期表达相对稳定，12:00出现小高峰，黑暗期逐渐升高，00:00达到峰值。饥饿组在饥饿第1天，光照期表达量与对照组无显著差异，黑暗期表达量显著升高，00:00达到峰值。第3天和第7天，光照期表达量明显低于对照组，黑暗期峰值表达量虽高于对照组，但表达时间提前且波动幅度减小，说明短期饥饿改变了G6PDH基因的昼夜节律表达模式，使其在黑暗期表达增强，光照期表达减弱，可能影响了肝脏的磷酸戊糖途径及相关代谢过程。在解毒相关基因方面，CYP1A基因在对照组中，光照期表达量逐渐上升，16:00达到峰值，黑暗期逐渐下降，04:00降至最低。饥饿组在饥饿第1天，光照期表达量显著低于对照组，峰值提前至12:00，黑暗期表达量显著升高，00:00达到峰值。第3天和第7天，光照期和黑暗期表达量均显著低于对照组，且表达波动不明显，表明短期饥饿改变了CYP1A基因的昼夜节律表达特征，使其在光照期表达提前且减弱，黑暗期表达增强后又减弱，可能影响了肝脏的解毒功能。UGT1A1基因在对照组中，光照期表达量先上升后下降，12:00达到峰值，黑暗期持续下降，04:00降至最低水平。饥饿组在饥饿第1天，光照期表达量显著低于对照组，峰值消失，黑暗期表达量有所升高，但无明显峰值。第3天和第7天，光照期和黑暗期表达量均显著低于对照组，且表达波动不明显，表明短期饥饿抑制了UGT1A1基因的表达，破坏了其正常的昼夜节律，可能影响了肝脏的Ⅱ相解毒反应。（此处插入图9：短期饥饿对锦鲫肝脏中重要功能基因昼夜节律性表达的影响图）4.2.3锦鲫短期饥饿后肝脏中生物钟相关基因间的相关性变化为深入了解短期饥饿对锦鲫肝脏中生物钟相关基因调控网络的影响，对饥饿组和对照组中生物钟相关基因间的相关性进行分析。在正常投喂的对照组中，Clock基因与Bmal1基因表达呈显著正相关（r=0.88，P<0.01），Per1基因与Clock基因、Bmal1基因表达也呈显著正相关（r=0.82，P<0.01；r=0.80，P<0.01），Cry1基因与Clock基因、Bmal1基因表达同样呈显著正相关（r=0.85，P<0.01；r=0.83，P<0.01），表明这些基因在肝脏生物钟调控中紧密协同。而在饥饿组中，生物钟相关基因间的相关性发生明显变化。Clock基因与Bmal1基因相关性虽仍为正相关，但相关性系数降低至r=0.68（P<0.05），表明短期饥饿削弱了它们之间的协同作用。Per1基因与Clock基因相关性系数降至r=0.58（P<0.05），与Bmal1基因相关性系数降至r=0.55（P<0.05），说明短期饥饿破坏了它们之间的紧密联系。Cry1基因与Clock基因相关性系数降至r=0.62（P<0.05），与Bmal1基因相关性系数降至r=0.60（P<0.05），表明短期饥饿影响了其在生物钟负反馈调节机制中的协同作用。Per2基因与其他生物钟基因相关性进一步减弱，相关性系数在r=0.35-0.55之间（P<0.05），说明短期饥饿对其在生物钟调控网络中的作用产生较大影响。Cry2基因与Clock基因、Bmal1基因相关性也有所减弱，相关性系数在r=0.45-0.65之间（P<0.05），表明短期饥饿改变了其与其他基因的相互作用关系，可能影响了其在生物钟调节过程中的功能。（此处插入图10：短期饥饿对锦鲫肝脏中生物钟相关基因间相关性的影响图）五、讨论5.1锦鲫肌肉和肝脏中生物钟相关基因的正常表达模式在正常投喂条件下，锦鲫肌肉和肝脏中生物钟相关基因均呈现出明显的昼夜节律性表达模式，这与大多数脊椎动物的生物钟基因表达特征一致。这种昼夜节律性表达模式表明生物钟在锦鲫的生理活动中起着重要的调控作用，使锦鲫能够适应环境的昼夜变化，维持生理功能的稳定。在锦鲫肌肉中，Clock基因在光照期逐渐上升，16:00达到峰值，随后在黑暗期逐渐下降，这可能与光照信号对Clock基因表达的调控有关。Clock基因编码的蛋白是生物钟调控网络中的核心转录因子，它与Bmal1基因编码的蛋白形成异源二聚体，激活下游Per和Cry基因的转录。Bmal1基因在光照期表达量相对较低，黑暗期逐渐升高，在00:00达到峰值，与Clock基因的表达模式存在一定的互补性，共同维持肌肉生物钟的稳定。Per1基因在光照期呈现先上升后下降的趋势，在12:00达到峰值，黑暗期持续下降，这表明Per1基因可能在光照中期对肌肉的生理过程具有重要的调控作用，参与调节肌肉的生长、发育和代谢等生理活动。Per2基因在光照期相对稳定，在16:00出现一个小高峰，黑暗期逐渐下降，说明Per2基因在光照后期对肌肉的生理功能具有一定的调节作用。Cry1基因在光照期逐渐升高，在16:00达到峰值，随后在黑暗期逐渐下降，表明Cry1基因可能在光照后期参与调节肌肉对光照信号的响应以及相关的代谢过程。Cry2基因在光照期相对较低，在12:00达到低谷，黑暗期逐渐升高，在00:00达到峰值，说明Cry2基因可能在黑暗期对肌肉的生理调节发挥重要作用。这些生物钟基因之间通过相互作用，形成了一个复杂的调控网络，共同维持锦鲫肌肉的正常生理节律。与代谢相关的Glut4基因和PGC-1α基因在锦鲫肌肉中的表达也呈现出明显的昼夜节律性。Glut4基因在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，这表明Glut4基因在光照后期可能参与促进葡萄糖向肌肉细胞的转运，为肌肉的生理活动提供能量，与生物钟基因Clock和Cry1的表达高峰时间相吻合，暗示它们之间可能存在协同作用，共同调节肌肉的能量代谢。PGC-1α基因在光照期相对稳定，在12:00出现一个小高峰，黑暗期逐渐升高，在00:00达到峰值，说明PGC-1α基因在黑暗期可能通过激活线粒体生物发生和脂肪酸氧化等代谢途径，为肌肉提供能量，与生物钟基因Bmal1和Cry2的表达模式具有一定的相似性，提示它们在调节肌肉代谢方面可能存在关联。在生长相关基因方面，IGF-1基因在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，表明IGF-1可能在光照后期通过促进蛋白质合成和细胞增殖等途径，促进锦鲫肌肉的生长，与生物钟基因Clock、Per1和Cry1的表达高峰时间一致，表明它们可能在调节锦鲫肌肉生长方面存在协同作用。MyoD基因在光照期呈现先上升后下降的趋势，在12:00达到峰值，黑暗期持续下降，说明MyoD基因在光照中期可能参与调控锦鲫肌肉细胞的分化和发育，与生物钟基因Per1的表达模式相匹配，暗示它们在肌肉发育过程中可能存在相互作用。在锦鲫肝脏中，Clock基因同样在光照期逐渐上升，16:00达到峰值，黑暗期逐渐下降，表明Clock基因在肝脏的生理活动中可能在光照后期发挥着关键作用，参与调控肝脏的代谢、解毒等生理过程。Bmal1基因在光照期表达量相对较低，在12:00达到低谷，黑暗期逐渐升高，在00:00达到峰值，与肌肉中Bmal1基因的表达模式相似，说明Bmal1基因在黑暗期对肝脏的生理调节起到重要作用，与Clock基因共同维持肝脏的正常生理节律。Per1基因在光照期呈现先上升后下降的趋势，在12:00达到峰值，黑暗期持续下降，表明Per1基因在光照中期对肝脏的生理过程具有重要的调控作用，可能参与调节肝脏的物质合成与分解、能量代谢等生理活动。Per2基因在光照期相对稳定，在16:00出现一个小高峰，黑暗期逐渐下降，说明Per2基因在光照后期对肝脏的生理功能具有一定的调节作用。Cry1基因在光照期逐渐升高，在16:00达到峰值，随后在黑暗期逐渐下降，表明Cry1基因在光照后期对肝脏的生理调节具有重要作用，可能参与调节肝脏对光照信号的响应以及相关的代谢过程。Cry2基因在光照期相对较低，在12:00达到低谷，黑暗期逐渐升高，在00:00达到峰值，说明Cry2基因可能在黑暗期对肝脏的生理调节发挥重要作用。这些生物钟基因在肝脏中的表达模式与肌肉中既有相似之处，也存在一定的差异，表明生物钟基因在不同组织中的调控机制可能存在组织特异性。与代谢相关的PEPCK基因和G6PDH基因在锦鲫肝脏中的表达也呈现出明显的昼夜节律性。PEPCK基因是糖异生途径中的关键酶基因，其在光照期逐渐上升，在16:00达到峰值，随后在黑暗期逐渐下降，表明PEPCK基因在光照后期可能参与促进肝脏中糖异生作用，为机体提供葡萄糖，以满足生理活动的能量需求，与生物钟基因Clock和Cry1的表达高峰时间相吻合，暗示它们之间可能存在协同作用，共同调节肝脏的能量代谢。G6PDH基因参与磷酸戊糖途径，其在光照期相对稳定，在12:00出现一个小高峰，黑暗期逐渐升高，在00: