探秘鳗弧菌感染下鮸鱼基因表达与miRNA调控网络

探秘鳗弧菌感染下鮸鱼基因表达与miRNA调控网络一、引言1.1研究背景在全球范围内，水产养殖业是食品生产领域中增长最为迅速的部分之一，对保障食物安全和促进经济发展具有重要意义。然而，随着养殖规模的不断扩大和集约化程度的提高，鱼类面临着日益严峻的病害威胁。鳗弧菌（Vibrioanguillarum）作为一种革兰氏阴性菌，是海水养殖鱼类的主要病原菌之一，能感染多种鱼类，如鲑鱼、虹鳟、鳗鲡、香鱼、鲈鱼、鳕鱼、大菱鲆、牙鲆、黄鱼等，给水产养殖业带来了巨大的经济损失。鳗弧菌感染会引发鱼类一系列严重的症状。感染初期，鱼体体表局部褪色，鳍条、鳍基部及鳃骨下部充血发红，肛门红肿；随着病情发展，肌肉组织出现弥撒性或点状出血，体表发黑，鳍部溃烂。解剖检验可见明显的黄色粘稠腹水，肠粘膜组织腐烂脱落，部分鱼肝脏坏死。不同鱼类感染鳗弧菌后的症状虽有所差异，但主要症状均表现为体表出血。例如，感染鳗弧菌的养殖牙鲆最明显的症状是鳍部严重出血以及体表溃烂；大菱鲆感染后则出现全身性的出血性败血症，鳍基部出血，眼球突出，角膜混浊，肝脏苍白，有时伴有腹水症状。鱼类在面对鳗弧菌感染时，会启动自身的免疫系统来抵御病原菌的入侵。免疫系统中的细胞免疫和体液免疫协同作用，免疫器官、免疫细胞和免疫分子共同参与免疫应答过程。在细胞免疫方面，吞噬细胞如中性粒细胞、巨噬细胞等负责吞噬和消化病原体；淋巴细胞包括T细胞和B细胞，参与特异性免疫应答。体液免疫中，抗体、溶菌酶、吞噬细胞等发挥重要作用。同时，鱼类免疫系统中还存在一些独特的免疫分子和通路，如Toll/IL-1受体超家族（Toll/IL-1receptorsuperfamily）相关的信号通路，在对抗外来入侵病原体的天然免疫中发挥着关键的作用。MicroRNA（miRNA）作为一类约含22核苷酸的内生非编码小分子RNA，在基因调控中扮演着重要角色。它通过与mRNA的3'非翻译区（3'UTR）以不完全互补配对的方式结合，从而引起该基因翻译受阻或者直接降解mRNA，实现对基因转录后水平的调控。在鱼类免疫领域，miRNA参与了鱼类对病原体感染的免疫应答过程，对鱼类的免疫调节起着重要作用。鮸鱼（Miichthysmiiuy）是一种重要的海水经济鱼类，在我国沿海地区广泛养殖。鳗弧菌感染对鮸鱼的养殖造成了严重的影响，制约了鮸鱼养殖业的发展。研究鮸鱼在鳗弧菌感染后的基因表达变化以及miRNA对这些基因的调控机制，对于深入了解鮸鱼的免疫防御机制、开发有效的病害防治策略具有重要的理论和实践意义。一方面，通过揭示miRNA与差异表达基因之间的相互作用关系，有助于从分子层面解析鮸鱼抵抗鳗弧菌感染的免疫调控网络；另一方面，为开发基于miRNA的新型鱼类病害防治技术提供理论依据，有望通过调节miRNA的表达来增强鮸鱼的免疫力，减少鳗弧菌感染带来的损失，促进鮸鱼养殖业的健康可持续发展。1.2研究目的与问题提出本研究旨在深入揭示miRNA对鳗弧菌感染鮸鱼基因的调控机制，为理解鮸鱼免疫防御机制以及开发有效的病害防治策略提供关键的理论基础。通过对鳗弧菌感染鮸鱼后的基因表达谱进行全面分析，鉴定出差异表达基因，并运用生物信息学和实验技术，系统地探究miRNA在其中所扮演的调控角色。基于此，本研究拟解决以下关键问题：首先，在鳗弧菌感染鮸鱼后，哪些miRNA的表达发生了显著变化，这些变化与鮸鱼的免疫应答过程存在怎样的关联？其次，这些差异表达的miRNA具体调控了哪些基因，它们之间的相互作用关系如何？进一步地，这些miRNA-基因调控网络在鮸鱼抵御鳗弧菌感染的免疫防御机制中发挥着怎样的功能？此外，能否通过人为调控这些关键的miRNA及其靶基因，来增强鮸鱼对鳗弧菌感染的抵抗力，从而为水产养殖中的病害防治提供新的策略和方法？二、研究综述2.1鱼类免疫学研究进展2.1.1鱼类免疫系统构成鱼类的免疫系统是其抵御病原体入侵、维持自身健康的重要防御体系，由免疫器官、免疫细胞和免疫分子共同构成。免疫器官是免疫细胞产生、分化、成熟、定居和增殖以及产生免疫应答的场所。胸腺是鱼类的中枢免疫器官，由淋巴细胞、淋巴母细胞、浆母细胞、分泌样细胞以及其他游离间充质细胞（如巨噬细胞、肌样细胞、肥大细胞等）组成，分布于由网状上皮细胞形成的基质网孔内，是T细胞的发源地，主要承担细胞免疫功能，硬骨鱼类胸腺中存在形态学上与高等脊椎动物相似的“血胸屏障”。肾脏在鱼类免疫中也起着关键作用，硬骨鱼类的肾脏分头肾、中肾和后肾三部分，头肾是鱼类继胸腺之后第二个发育的免疫器官，同时具有造血功能，在不依赖抗原刺激时，头肾可以产生红细胞和B淋巴细胞等细胞，相当于哺乳动物的骨髓，受抗原刺激后，头肾和后肾造血实质细胞出现增生，且存在抗体产生细胞。脾脏则是各种中性粒细胞、红细胞等产生、储存和成熟的主要场所，不仅具有产生各种免疫细胞的功能，还具备造血功能，对于软骨鱼和硬骨鱼而言，虽然脾脏结构存在一定差异，如软骨鱼的脾脏分为红髓和白髓，硬骨鱼无明显的红髓、白髓之分，但都具有免疫功能。此外，粘膜淋巴组织（Mucosa-associatedlymphoidtissue,MALT）广泛分布于鱼类的消化道、呼吸道和泌尿生殖道等粘膜表面，是鱼类抵御病原体入侵的第一道防线，在鱼类的免疫防御中发挥着重要作用。免疫细胞是免疫系统的核心组成部分，包括巨噬细胞、淋巴细胞、中性粒细胞、嗜酸性粒细胞等。巨噬细胞具有强大的吞噬和消化病原体的能力，同时能够分泌多种细胞因子，参与免疫调节和炎症反应；淋巴细胞主要分为T淋巴细胞和B淋巴细胞，T淋巴细胞主要介导细胞免疫，通过分泌细胞因子参与机体的非特异性免疫反应，B淋巴细胞主要介导体液免疫，以分泌抗体和表面抗体来参与机体的特异性免疫反应。中性粒细胞能够快速迁移到感染部位，通过吞噬和释放杀菌物质来抵御病原体；嗜酸性粒细胞则在抗寄生虫感染和过敏反应中发挥重要作用。免疫分子是免疫应答过程中的效应物质，包括细胞因子、补体、免疫球蛋白、溶菌酶等。细胞因子是一类由免疫细胞和某些非免疫细胞分泌的小分子蛋白质，具有调节免疫细胞功能、介导炎症反应等多种生物学活性，如白细胞介素、干扰素、肿瘤坏死因子等；补体系统是鱼类免疫系统中重要的组成部分，能同时参与特异性免疫和非特异性免疫反应，由大约30余种酶、酶抑制因子和受体构成，具有较强的调理作用，能增强体液和细胞介导的特异性免疫；免疫球蛋白是鱼类特异性免疫的关键分子，虽然鱼类的免疫球蛋白种类相对哺乳动物较少，目前从鱼类中分离得到的免疫球蛋白分子主要为IgM，但它能参与鱼类的特异性免疫反应；溶菌酶存在于许多鱼类的体表、肠道黏液、血清和巨噬细胞中，能作用于病原细胞壁，进行有效的攻击，从而发挥免疫防御作用。2.1.2鱼类免疫机制概述鱼类的免疫机制分为先天性免疫和获得性免疫，二者相互协作，共同抵御病原体的入侵。先天性免疫是鱼类与生俱来的防御机制，是鱼类抵御病原体的第一道防线，具有快速响应、非特异性等特点，在病原体入侵的早期发挥关键作用。模式识别受体（PatternRecognitionReceptors，PRRs）是先天性免疫的重要组成部分，能够识别病原体表面的保守分子结构，即病原体相关分子模式（Pathogen-AssociatedMolecularPatterns，PAMPs），如细菌的脂多糖、肽聚糖，病毒的双链RNA等。常见的模式识别受体包括Toll样受体（Toll-likeReceptors，TLRs）、NOD样受体（NOD-likeReceptors，NLRs）、RIG-I样受体（RIG-I-likeReceptors，RLRs）等。以Toll样受体为例，它是一类Ⅰ型跨膜蛋白，可分为胞膜外区、胞浆区和跨膜区三部分，胞膜外区主要行使识别受体及与其他辅助受体结合形成受体复合物的功能，其胞浆区与IL-1R家族成员胞浆区高度同源，称为Toll-IL-1受体结构域（Toll-IL-1receptordomain，TIR结构域），TIR具有嗜同性相互作用，藉此来募集下游含有TIR的信号分子，组成信号复合体。当Toll样受体识别到病原体相关分子模式后，会激活一系列下游信号通路，如MyD88依赖的信号通路和MyD88非依赖的信号通路，最终导致核因子κB（NF-κB）等转录因子的激活，促使细胞因子、趋化因子等免疫分子的表达，引发炎症反应，招募免疫细胞到感染部位，以清除病原体。获得性免疫是鱼类在接触病原体后逐渐产生的特异性免疫应答，具有特异性、记忆性等特点，在病原体入侵的后期发挥重要作用。当病原体突破先天性免疫防线后，抗原呈递细胞（如巨噬细胞、树突状细胞等）会摄取、加工病原体抗原，并将其呈递给T淋巴细胞和B淋巴细胞。T淋巴细胞被激活后，会分化为效应T细胞和记忆T细胞，效应T细胞能够直接杀伤被病原体感染的细胞，或者分泌细胞因子来调节免疫应答；B淋巴细胞被激活后，会分化为浆细胞，浆细胞分泌特异性抗体，抗体与病原体结合，通过中和、凝集、调理等作用，清除病原体。此外，记忆T细胞和记忆B细胞能够记住病原体的抗原信息，当再次接触相同病原体时，能够迅速启动免疫应答，产生更强的免疫反应，这就是获得性免疫的记忆性。2.1.3鮸鱼的生物学特性与免疫特点鮸鱼（Miichthysmiiuy）属于鲈形目石首鱼科，是一种重要的海水经济鱼类，在我国沿海地区广泛养殖。鮸鱼身体较长，侧扁，背部和腹部呈浅弧形，头中等大，吻短而钝尖，无须。其体色呈较暗的灰褐色带部分紫绿色，腹部为灰白色，体表无明显条纹。鮸鱼属于近海温水性的近底层鱼类，喜欢栖息在15-75米深且底质为泥沙或者泥的海底，是广温广盐性鱼类，以鮸鱼幼鱼为例，可适应温度范围在4-32摄氏度，盐度的范围在12.1-52.6。鮸鱼具有垂直活动现象，白天下沉，夜间上浮到水面，同时它还是一种小区域性的洄游鱼类，每到产卵季节临近，鱼群才会较为集中。在觅食方面，鮸鱼属于底栖肉食性鱼类，其摄取的食物等级随着个体发育而增加，体重小于750克的仔鱼以摄食轮虫、小球藻、小型桡足类等为主，750克以上的可摄食小带鱼、沙丁鱼等，成鱼的食物主要以黄鲫、青鳞鱼等以及对虾、鼓虾、毛虾等。由于鮸鱼长期生活在水下光线较弱的地方，视觉系统退化，但嗅觉灵敏，所以对于腥味重的食物摄食率较高。在免疫方面，鮸鱼具有与其他鱼类相似的免疫器官和免疫细胞，其免疫机制也包括先天性免疫和获得性免疫。然而，鮸鱼对鳗弧菌感染具有较高的易感性。鳗弧菌感染会导致鮸鱼出现一系列严重的症状，如体表出血、鳍条溃烂、内脏器官病变等，严重影响鮸鱼的生长和生存，给鮸鱼养殖业带来巨大的经济损失。研究表明，鳗弧菌感染鮸鱼后，会引发鮸鱼免疫系统的一系列应答反应，包括免疫细胞的活化、免疫分子的表达变化等，但鮸鱼在抵御鳗弧菌感染过程中，其免疫调节机制尚不完全清楚，尤其是miRNA在其中的调控作用，有待进一步深入研究。2.2Toll/IL-1R家族及信号通路2.2.1Toll/IL-1R家族成员与结构Toll/IL-1R家族在免疫调控中发挥着关键作用，是先天性免疫研究的重要领域。该家族成员在进化上高度保守，从果蝇到哺乳动物都存在其同源分子。在哺乳动物中，Toll样受体（Toll-likeReceptors，TLRs）是Toll/IL-1R家族的重要组成部分，目前已发现人类有10个TLR家族成员，如TLR1位于4p14、TLR2位于4q32、TLR3位于4q35、TLR4位于9q32-33等。Toll/IL-1R家族成员均为Ⅰ型跨膜蛋白，其结构可分为胞膜外区、胞浆区和跨膜区三部分。胞膜外区是受体与配体结合的关键部位，主要行使识别受体及与其他辅助受体结合形成受体复合物的功能。以TLRs为例，其胞膜外区含有17-31个亮氨酸富集的重复序列（Leucinerichrepeats，LRRs），这些LRRs结构域赋予了TLRs识别多种病原体相关分子模式（Pathogen-AssociatedMolecularPatterns，PAMPs）的能力，如细菌的脂多糖（LPS）、肽聚糖，病毒的双链RNA等。不同TLR成员的胞膜外区结构存在一定差异，这决定了它们对不同PAMPs的特异性识别，例如TLR4主要识别革兰氏阴性菌的脂多糖，而TLR2则可识别革兰氏阳性菌的肽聚糖等多种PAMPs。胞浆区是信号传导的关键区域，Toll/IL-1R家族成员的胞浆区与IL-1R家族成员胞浆区高度同源，该区称为Toll-IL-1受体结构域（Toll-IL-1receptordomain，TIR结构域）。TIR结构域具有嗜同性相互作用，能够募集下游含有TIR的信号分子，组成信号复合体，从而启动下游信号转导通路。在哺乳动物中，一些胞浆蛋白也存在TIR域，如2种信号接受蛋白，MyD88和TIR域的接受子蛋白（TIRdomain-containingadaptorprotein，TIRAP），它们在TIR信号传导中发挥着重要作用。MyD88是TLRs信号通路中重要的接头蛋白，几乎参与所有TLRs介导的信号转导，通过其TIR结构域与TLRs的TIR结构域相互作用，招募下游的白细胞介素-1受体相关激酶（IL-1receptor-associatedkinases，IRAKs）等信号分子，激活下游信号通路。2.2.2Toll/IL-1R信号通路解析Toll/IL-1R信号通路在鱼类抵御病原体入侵的先天性免疫应答中起着核心作用，其激活过程涉及一系列复杂的分子事件和信号转导途径。当病原体入侵鱼类机体时，其表面的病原体相关分子模式（PAMPs）首先被Toll/IL-1R家族成员识别。以TLRs为例，不同的TLR成员识别不同的PAMPs，如TLR4识别革兰氏阴性菌的脂多糖（LPS），TLR2识别革兰氏阳性菌的肽聚糖等。这种识别过程具有高度的特异性，是Toll/IL-1R信号通路激活的起始步骤。当TLR识别到PAMP后，其胞内的Toll-IL-1受体结构域（TIR结构域）发生构象变化，从而招募含有TIR结构域的接头蛋白，启动下游信号转导。在Toll/IL-1R信号通路中，主要存在MyD88依赖和MyD88非依赖两条信号转导途径。MyD88依赖的信号通路是Toll/IL-1R信号通路的经典途径，大多数TLRs激活后都会通过该途径传递信号。在这条途径中，MyD88作为关键的接头蛋白，通过其TIR结构域与TLRs的TIR结构域相互作用，招募白细胞介素-1受体相关激酶（IL-1receptor-associatedkinases，IRAKs）家族成员。IRAKs被招募后，会发生磷酸化激活，进而激活肿瘤坏死因子受体相关因子6（TNFreceptor-associatedfactor6，TRAF6）。TRAF6激活后，会通过一系列激酶级联反应，激活转化生长因子β激活激酶1（TGF-β-activatedkinase1，TAK1）。TAK1进一步激活核因子κB（NF-κB）抑制蛋白激酶（IκBkinases，IKKs），使IκB磷酸化并降解，从而释放NF-κB。NF-κB进入细胞核，与相关基因的启动子区域结合，启动炎症细胞因子（如白细胞介素-1β、白细胞介素-6、肿瘤坏死因子α等）、趋化因子等免疫相关基因的转录，引发炎症反应，招募免疫细胞到感染部位，以清除病原体。MyD88非依赖的信号通路主要由TLR3和TLR4激活。在这条途径中，接头蛋白TRIF（TIR-domain-containingadapter-inducinginterferon-β）发挥关键作用。当TLR3或TLR4识别PAMP后，招募TRIF，TRIF通过与TRAF3相互作用，激活干扰素调节因子3（Interferonregulatoryfactor3，IRF3）。IRF3磷酸化后形成二聚体，进入细胞核，与干扰素刺激反应元件（Interferon-stimulatedresponseelement，ISRE）结合，启动干扰素（如IFN-β）等抗病毒基因的转录，发挥抗病毒免疫作用。同时，TRIF还可以通过激活RIP1（Receptor-interactingprotein1），进一步激活NF-κB，引发炎症反应。Toll/IL-1R信号通路的激活对于鱼类免疫应答具有重要意义。它能够迅速启动先天性免疫反应，在病原体入侵的早期阶段发挥关键的防御作用。通过激活NF-κB等转录因子，上调炎症细胞因子和趋化因子的表达，招募巨噬细胞、中性粒细胞等免疫细胞到感染部位，增强吞噬作用和杀菌能力，从而有效清除病原体。同时，Toll/IL-1R信号通路的激活还可以调节获得性免疫反应，为后续的特异性免疫应答奠定基础。然而，Toll/IL-1R信号通路的过度激活也可能导致炎症反应失控，引发组织损伤和免疫病理反应，因此其激活过程受到严格的调控。2.3miRNA研究进展2.3.1miRNA的生物合成miRNA的生物合成是一个复杂且精细的过程，涉及多个关键步骤和酶的参与。首先，miRNA基因通过RNA聚合酶Ⅱ转录生成具有发夹结构的初级转录本（primarymiRNA，pri-miRNA）。pri-miRNA通常长度可达数百至数千个核苷酸，其结构包含一个或多个茎环结构以及侧翼序列。在细胞核内，pri-miRNA被一种由Drosha酶和DGCR8蛋白组成的Microprocessor复合物识别并切割。Drosha是一种RNaseⅢ酶，它能够在pri-miRNA的茎环结构处进行特异性切割，将pri-miRNA加工成约70-100个核苷酸的前体miRNA（precursormiRNA，pre-miRNA）。pre-miRNA依然保持着发夹结构，其5'端带有磷酸基团，3'端则有2个核苷酸的突出。随后，pre-miRNA在转运蛋白exportin-5的作用下，从细胞核被转运至细胞质中。exportin-5与pre-miRNA结合形成复合物，通过核孔复合体进入细胞质。在细胞质中，pre-miRNA被另一种RNaseⅢ酶Dicer识别并进一步切割。Dicer能够识别pre-miRNA的发夹结构，将其切割成双链小RNA，长度约为21-25个核苷酸。这一双链小RNA的两条链分别被称为引导链（guidestrand）和过客链（passengerstrand）。接下来，双链小RNA与Argonaute（Ago）蛋白结合形成RNA诱导沉默复合体（RNA-inducedsilencingcomplex，RISC）。在RISC组装过程中，过客链逐渐被降解，而引导链则保留在Ago蛋白中，形成成熟的miRNA复合物。成熟的miRNA复合物通过其5'端的种子序列（seedsequence）与靶mRNA的3'非翻译区（3'UTR）进行互补配对，从而实现对靶基因的调控作用。除了上述经典的miRNA生物合成途径外，还存在一些非经典途径。例如，某些miRNA可以不依赖Drosha或Dicer进行加工生成。在一些情况下，pri-miRNA可以直接被Dicer切割生成成熟的miRNA，这种途径被称为Drosha-independent途径。此外，还有一些miRNA是由mRNA的内含子加工而来，这些内含子在mRNA剪接过程中形成特殊的结构，直接被Dicer识别并加工成miRNA，而不需要经过Drosha的切割步骤。2.3.2miRNA的作用机制miRNA主要通过与靶mRNA的3'非翻译区（3'UTR）互补配对，来实现对基因表达的调控，其作用方式主要包括翻译抑制和mRNA降解两种。当miRNA与靶mRNA的3'UTR以不完全互补配对的方式结合时，主要抑制靶mRNA的翻译过程。在翻译起始阶段，miRNA与靶mRNA结合后，会阻碍核糖体与mRNA的结合，从而抑制翻译起始复合物的形成。此外，miRNA还可能影响翻译延伸和终止过程。例如，miRNA与靶mRNA结合后，会导致核糖体在mRNA上的移动速度减慢，从而抑制翻译延伸；或者miRNA会使翻译提前终止，导致产生不完整的蛋白质。在哺乳动物细胞中，let-7家族的miRNA可以通过与靶mRNA的3'UTR结合，抑制其翻译过程，从而调控细胞的增殖、分化和凋亡等过程。当miRNA与靶mRNA的3'UTR存在较高程度的互补配对时，尤其是在种子序列区域完全互补时，miRNA会介导靶mRNA的降解。在这种情况下，miRNA与靶mRNA结合后，会招募核酸内切酶等相关蛋白，对靶mRNA进行切割，使其降解为小片段。这些小片段随后会被细胞内的核酸外切酶进一步降解，从而实现对靶mRNA的完全清除。在植物中，miRNA介导的mRNA降解是一种常见的调控方式。例如，拟南芥中的miR164可以与NAC1基因的mRNA互补配对，介导其降解，从而调控植物的侧根发育。miRNA对靶基因的调控具有多效性和复杂性。一个miRNA可以同时调控多个靶基因的表达，而一个靶基因也可能受到多个miRNA的调控。这种复杂的调控网络使得miRNA能够在细胞内发挥广泛而精细的调控作用，参与细胞的各种生理过程，如细胞增殖、分化、凋亡、代谢等。2.3.3miRNA在免疫系统中的功能miRNA在免疫系统中发挥着至关重要的调控作用，参与免疫细胞的发育、分化和免疫应答等多个环节。在免疫细胞发育方面，miRNA对免疫细胞的分化和成熟起着关键的调节作用。以T淋巴细胞为例，miR-181a在T淋巴细胞发育过程中表达上调，它可以通过抑制多个负调控因子的表达，促进T淋巴细胞的分化和成熟。研究表明，miR-181a过表达的小鼠中，T淋巴细胞的数量和功能都得到显著增强。在B淋巴细胞发育过程中，miR-150的表达变化对B淋巴细胞的分化和成熟也具有重要影响。miR-150在未成熟B淋巴细胞中表达较低，而在成熟B淋巴细胞中表达升高。miR-150通过靶向抑制转录因子c-Myb的表达，促进B淋巴细胞从幼稚阶段向成熟阶段的分化。在免疫应答过程中，miRNA参与调节免疫细胞的活化和细胞因子的产生。当免疫细胞受到病原体刺激时，miRNA的表达会发生显著变化。在巨噬细胞中，miR-155在受到脂多糖（LPS）刺激后表达上调，它可以通过调节多个靶基因的表达，促进巨噬细胞的活化和炎症细胞因子（如肿瘤坏死因子α、白细胞介素-6等）的产生。miR-155通过靶向抑制SHIP1基因的表达，激活PI3K-Akt信号通路，从而增强巨噬细胞的免疫应答能力。相反，一些miRNA则起到负调控免疫应答的作用。miR-146a在免疫细胞受到刺激后表达上调，它可以通过靶向抑制肿瘤坏死因子受体相关因子6（TRAF6）和白细胞介素-1受体相关激酶1（IRAK1）等信号分子，抑制NF-κB信号通路的激活，从而负调控炎症细胞因子的产生，防止免疫应答过度激活。miRNA还参与调节免疫细胞之间的相互作用和免疫记忆的形成。在T淋巴细胞和B淋巴细胞的相互作用中，miRNA可以通过调节细胞表面分子的表达，影响两者之间的信号传递和协同作用。在免疫记忆形成过程中，miRNA也发挥着重要作用。研究发现，miR-125b在记忆性T淋巴细胞中表达上调，它可以通过调节相关基因的表达，维持记忆性T淋巴细胞的存活和功能，从而参与免疫记忆的形成和维持。三、研究设计与方法3.1实验材料实验所用鮸鱼购自[具体养殖场名称]，为健康的当年生幼鱼，体长约[X]cm，体重约[X]g。将鮸鱼暂养于实验室循环水养殖系统中，水温控制在[20±1]℃，盐度为[28±2]‰，pH值维持在[7.8-8.2]，每天投喂适量的商业饲料，暂养一周使其适应实验室环境后进行后续实验。鳗弧菌菌种（Vibrioanguillarum）购自[菌种保藏中心名称]，编号为[具体编号]。菌种保存于-80℃冰箱中，使用前将其从冰箱取出，在冰上缓慢融化。然后在无菌条件下，将菌种接种于含有3%NaCl的LB固体培养基上，置于28℃恒温培养箱中培养24-48h，待长出单菌落。挑选单菌落接种于含有3%NaCl的LB液体培养基中，在28℃、200rpm的摇床中振荡培养至对数生长期，备用。实验中使用的细胞系为鮸鱼头肾巨噬细胞系（Miiuycroakerhead-kidneymacrophagecellline，MHK），由本实验室前期建立并保存。细胞培养于含有10%胎牛血清（Fetalbovineserum，FBS）、100U/mL青霉素和100μg/mL链霉素的RPMI-1640培养基中，置于37℃、5%CO₂的细胞培养箱中培养，定期换液传代。转染试剂选用Lipofectamine3000（Invitrogen公司），按照其说明书进行操作，用于将miRNA模拟体、抑制体和载体转染至细胞中。载体选用pGL3-basic载体（Promega公司），用于构建荧光素酶报告基因载体。引物由[引物合成公司名称]合成，包括用于荧光定量PCR的特异性引物、用于构建载体的引物等，引物序列通过NCBI数据库和相关文献进行设计，并经过BLAST比对验证其特异性。miRNA模拟体（miRNAmimics）和miRNA抑制体（miRNAinhibitors）均购自[生物公司名称]，分别用于模拟miRNA的过表达和抑制miRNA的表达。miRNA模拟体是与成熟miRNA序列相同的双链RNA，能够增强细胞内miRNA的表达水平；miRNA抑制体是经过化学修饰的单链RNA，能够特异性地与细胞内的miRNA结合，从而抑制其功能。3.2实验方法3.2.1鮸鱼巨噬细胞的分离与培养在无菌条件下，从暂养适应后的鮸鱼中选取个体，使用10mg/L的MS-222（三卡因甲烷磺酸盐）溶液进行麻醉处理。将麻醉后的鮸鱼置于超净工作台中，用75%酒精棉球擦拭鱼体表面进行消毒。随后，迅速剪开鱼体腹腔，取出脾脏，将其放入含有预冷的无菌PBS（磷酸盐缓冲液，pH7.4）的培养皿中，轻轻漂洗，以去除脾脏表面的血液和杂质。将漂洗后的脾脏转移至另一无菌培养皿中，加入适量含有10%胎牛血清（Fetalbovineserum，FBS）、100U/mL青霉素和100μg/mL链霉素的RPMI-1640培养基，用眼科剪将脾脏剪成1mm³左右的小块。接着，将剪碎的脾脏组织转移至15mL离心管中，加入适量的0.25%胰蛋白酶-EDTA溶液，置于37℃水浴锅中消化15-20min，期间每隔5min轻轻振荡一次，使组织块与胰蛋白酶充分接触。消化结束后，向离心管中加入等体积的含有10%FBS的RPMI-1640培养基，以终止胰蛋白酶的消化作用。然后，将离心管以1000rpm的转速离心5min，弃去上清液。用含有10%FBS的RPMI-1640培养基重悬细胞沉淀，并用200目细胞筛过滤，去除未消化的组织块，得到单细胞悬液。将单细胞悬液转移至25cm²细胞培养瓶中，置于37℃、5%CO₂的细胞培养箱中培养。培养2h后，轻轻晃动培养瓶，去除未贴壁的细胞，然后加入新鲜的含有10%FBS的RPMI-1640培养基，继续培养。此后，每隔2-3天更换一次培养基，待细胞融合度达到80%-90%时，进行传代培养。传代时，先用PBS轻轻漂洗细胞两次，然后加入适量的0.25%胰蛋白酶-EDTA溶液，消化2-3min，待细胞变圆并脱离瓶壁后，加入含有10%FBS的RPMI-1640培养基终止消化，吹打细胞使其均匀分散，按1:2或1:3的比例接种到新的细胞培养瓶中继续培养。3.2.2病原菌感染实验设计将培养至对数生长期的鳗弧菌接种于含有3%NaCl的LB液体培养基中，在28℃、200rpm的摇床中振荡培养至OD₆₀₀值为0.6-0.8，此时的鳗弧菌处于对数生长期，活力较强。然后，将菌液以5000rpm的转速离心5min，弃去上清液，用无菌PBS重悬菌体，并调整菌液浓度至1×10⁷CFU/mL（菌落形成单位/毫升）备用。选取生长状态良好、大小均匀的鮸鱼头肾巨噬细胞，以1×10⁶个/孔的密度接种于6孔细胞培养板中，每孔加入2mL含有10%FBS的RPMI-1640培养基，置于37℃、5%CO₂的细胞培养箱中培养，待细胞贴壁生长至融合度达到80%左右时，进行感染实验。实验共设置3组，分别为感染组、对照组和空白组，每组设置3个复孔。感染组每孔加入100μL浓度为1×10⁷CFU/mL的鳗弧菌菌液，使感染复数（MOI）为10:1（细菌与细胞数量之比）；对照组每孔加入100μL无菌PBS，以排除PBS对实验结果的影响；空白组不做任何处理，仅加入2mL含有10%FBS的RPMI-1640培养基，用于观察细胞的正常生长状态。将接种后的6孔细胞培养板置于37℃、5%CO₂的细胞培养箱中继续培养，分别在感染后0h、3h、6h、12h、24h等时间点收集细胞，用于后续的RNA提取和相关分析。在每个时间点，小心吸去培养孔中的培养基，用预冷的PBS轻轻漂洗细胞3次，以去除未感染的细菌和杂质，然后按照RNA提取试剂盒的说明书进行RNA提取。3.2.3RNA提取与质量检测采用Trizol试剂法从感染和未感染鳗弧菌的鮸鱼头肾巨噬细胞以及鮸鱼组织（如脾脏、肝脏等）中提取总RNA。在超净工作台中，将收集的细胞或组织样品加入适量的Trizol试剂，用移液器反复吹打，使样品充分裂解。将裂解液转移至1.5mL无RNA酶的EP管中，室温静置5min，使核酸蛋白复合物完全解离。向EP管中加入1/5体积的氯仿，盖紧管盖，手动剧烈振荡15s，使溶液充分混匀，然后室温静置5min。将EP管置于4℃、12000g的离心机中离心15min，此时溶液会分为三层，上层为无色透明的水相，含有RNA；中间层为白色的蛋白质层；下层为红色的酚氯仿相，含有DNA和蛋白质等杂质。小心吸取上层水相转移至新的1.5mL无RNA酶的EP管中，加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min，使RNA沉淀。将EP管置于4℃、12000g的离心机中离心10min，此时RNA会沉淀在管底，形成白色或透明的胶状沉淀。小心弃去上清液，加入1mL预冷的75%乙醇（用DEPC水配制），轻轻洗涤RNA沉淀，以去除杂质和残留的异丙醇。将EP管置于4℃、7500g的离心机中离心5min，弃去上清液，尽量吸尽乙醇。将EP管敞口置于超净工作台中干燥5-10min，使RNA沉淀中的乙醇完全挥发。向干燥后的RNA沉淀中加入适量的无RNA酶的DEPC水，用移液器反复吹打，使RNA充分溶解。将提取的RNA样品保存于-80℃冰箱中备用。使用超微量紫外分光光度计测定RNA样品在260nm和280nm处的吸光值，计算OD₂₆₀/OD₂₈₀比值，以评估RNA的纯度。一般来说，纯RNA样品的OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，如果比值低于1.8，表明RNA样品可能存在蛋白质污染；如果比值高于2.0，可能存在RNA降解或DNA污染。同时，根据OD₂₆₀值计算RNA样品的浓度，公式为：RNA浓度（μg/μL）=OD₂₆₀×稀释倍数×40/1000。采用1%琼脂糖凝胶电泳检测RNA的完整性。配制1×TAE电泳缓冲液，称取适量的琼脂糖，加入到1×TAE电泳缓冲液中，加热溶解，待溶液冷却至60℃左右时，加入适量的GoldView核酸染料，充分混匀。将温热的琼脂糖凝胶溶液倒入凝胶模具中，插入梳子，待凝胶凝固后，小心拔出梳子。将凝胶放入电泳槽中，加入1×TAE电泳缓冲液，使液面没过凝胶。取适量的RNA样品与上样缓冲液混合，然后加入到凝胶的加样孔中，同时加入RNAMarker作为分子量标准。在100V的电压下电泳30-40min，电泳结束后，在紫外凝胶成像系统下观察并拍照。如果RNA样品的28S和18SrRNA条带清晰、明亮，且28S条带的亮度约为18S条带的2倍，表明RNA的完整性较好；如果条带模糊或出现拖尾现象，说明RNA可能存在降解。3.2.4cDNA合成与文库构建以提取的总RNA为模板，使用反转录试剂盒进行cDNA合成。在冰上配制反转录反应体系，体系中包含5×反转录缓冲液、dNTP混合物、随机引物、逆转录酶、RNA酶抑制剂和总RNA模板等。将反应体系轻轻混匀后，短暂离心，使液体集中在管底。将反应管置于PCR仪中，按照以下程序进行反转录反应：37℃孵育15min，85℃加热5s，然后4℃保存。反应结束后，得到的cDNA产物可用于后续的实验分析。构建cDNA文库时，首先使用FragmentationBuffer将cDNA进行片段化处理，使其长度分布在200-500bp之间。然后，对片段化的cDNA进行末端修复，在DNA聚合酶、dNTP等作用下，将cDNA片段的末端修复为平端，并在3'端加上一个“A”碱基。接着，将带有“A”碱基的cDNA片段与接头进行连接，接头包含了用于PCR扩增和测序的引物结合位点。连接反应完成后，通过PCR扩增富集带有接头的cDNA片段，得到cDNA文库。在PCR扩增过程中，需要优化引物浓度、退火温度等条件，以保证扩增效率和文库质量。扩增结束后，使用磁珠纯化试剂盒对PCR产物进行纯化，去除未反应的引物、dNTP等杂质，得到高质量的cDNA文库。最后，使用Qubit荧光定量仪对文库浓度进行精确测定，并通过Agilent2100生物分析仪检测文库的插入片段大小和质量，确保文库符合高通量测序的要求。3.2.5高通量测序与数据分析选择IlluminaHiSeq测序平台对构建好的cDNA文库进行高通量测序。将cDNA文库与测序芯片进行杂交，在测序过程中，DNA聚合酶会根据碱基互补配对原则，依次添加dNTP，同时释放出荧光信号，测序仪通过检测荧光信号来确定每个碱基的序列。测序完成后，得到的原始数据以FASTQ格式存储，包含了测序序列（reads）以及对应的质量分数信息。对原始测序数据进行预处理，使用FastQC软件对原始数据进行质量评估，检查数据的质量分布、碱基组成、GC含量等指标。使用Trimmomatic软件去除低质量的reads、接头序列以及含有过多N（未知碱基）的reads。通过这些预处理步骤，得到高质量的cleanreads，为后续的数据分析提供可靠的数据基础。利用TopHat软件将cleanreads比对到鮸鱼的参考基因组上，确定每个read在基因组上的位置。使用Cufflinks软件根据比对结果计算基因的表达量，通常以每百万映射reads中来自某基因每千碱基长度的reads数（FPKM，FragmentsPerKilobaseofexonperMillionreadsmapped）来表示基因的表达水平。通过比较感染组和对照组中基因的表达量，使用DESeq2软件进行差异表达分析，筛选出在鳗弧菌感染后表达发生显著变化的基因。设定差异表达基因的筛选标准为：|log₂(FoldChange)|≥1且FDR（FalseDiscoveryRate）＜0.05。其中，FoldChange表示感染组与对照组基因表达量的比值，FDR用于校正多重检验中的假阳性率。对筛选出的差异表达基因进行功能注释和富集分析，使用DAVID数据库对差异表达基因进行GO（GeneOntology）功能富集分析和KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析，以了解差异表达基因参与的生物学过程、细胞组成和分子功能，以及相关的信号通路，从而深入探究鮸鱼在鳗弧菌感染后的免疫应答机制。3.2.6miRNA筛选与靶基因预测运用生物信息学工具，如miRDeep2、miRanda、TargetScan等，对高通量测序数据进行分析，筛选出与差异表达基因相关的miRNA。首先，使用miRDeep2软件对测序数据进行处理，识别出已知的miRNA和预测新的miRNA。根据miRNA的表达量变化，筛选出在鳗弧菌感染后表达显著上调或下调的miRNA。利用miRanda和TargetScan等软件预测miRNA的靶基因。这些软件基于miRNA与靶基因3'UTR互补配对的原理，通过计算miRNA与靶基因之间的结合能、种子序列匹配等参数，预测可能的靶基因。在预测过程中，设定一定的筛选标准，如结合能阈值、种子序列匹配长度等，以提高预测的准确性。对预测得到的靶基因进行功能注释和富集分析，使用DAVID数据库进行GO功能富集分析和KEGG通路富集分析，了解靶基因参与的生物学过程和信号通路，从而初步推断miRNA在鮸鱼免疫应答中的调控作用。通过构建miRNA-靶基因调控网络，直观地展示miRNA与靶基因之间的相互作用关系，进一步分析miRNA在免疫调控中的核心作用和关键节点。3.2.7实时荧光定量PCR验证根据高通量测序结果，选取部分差异表达基因和与之相关的miRNA进行实时荧光定量PCR验证。使用PrimerPremier5.0软件设计特异性引物，引物设计时遵循以下原则：引物长度为18-25bp，GC含量在40%-60%之间，引物3'端避免出现连续的3个以上相同碱基，引物之间避免形成二聚体和发夹结构。将设计好的引物序列在NCBI数据库中进行BLAST比对，确保引物的特异性。以提取的总RNA为模板，按照反转录试剂盒的说明书合成cDNA。使用SYBRGreen荧光染料法进行实时荧光定量PCR反应，反应体系包含2×SYBRGreenPCRMasterMix、上下游引物、cDNA模板和无菌水。将反应体系轻轻混匀后，加入到96孔PCR板中，每个样品设置3个技术重复。在PCR仪上按照以下程序进行反应：95℃预变性30s，然后进行40个循环，每个循环包括95℃变性5s、60℃退火30s。在每个循环的退火阶段，检测荧光信号的强度，根据荧光信号的变化绘制扩增曲线。采用2⁻ΔΔCt法计算基因的相对表达量。其中，ΔCt=Ct（目的基因）-Ct（内参基因），ΔΔCt=ΔCt（感染组）-ΔCt（对照组）。内参基因选择在不同组织和处理条件下表达相对稳定的基因，如β-actin、GAPDH等。通过比较实时荧光定量PCR结果与高通量测序数据，验证差异表达基因和miRNA的表达变化趋势是否一致，以确保高通量测序结果的可靠性。3.2.8双荧光素酶报告基因实验构建荧光素酶报告载体，验证miRNA与靶基因3'UTR的相互作用。从鮸鱼基因组DNA中扩增出靶基因3'UTR序列，使用限制性内切酶将其克隆到pGL3-basic载体的多克隆位点中，构建成pGL3-target载体。同时，构建突变型的pGL3-mut-target载体，即将靶基因3'UTR中与miRNA结合的种子序列区域进行突变。将构建好的载体进行测序验证，确保序列的准确性。将miRNA模拟体（mimics）或阴性对照（NC）与pGL3-target载体或pGL3-mut-target载体共转染至HEK293T细胞中。使用Lipofectamine3000转染试剂，按照其说明书进行操作。转染后48h，收集细胞，使用双荧光素酶报告基因检测试剂盒检测荧光素酶活性。该试剂盒利用萤火虫荧光素酶和海肾荧光素酶的不同底物，分别检测萤火虫荧光素酶活性和海肾荧光素酶活性。以海肾荧光素酶活性作为内参，校正转染效率，计算萤火虫荧光素酶活性与海肾荧光素酶活性的比值。如果miRNA能够与靶基因3'UTR相互作用，抑制荧光素酶的表达，那么转染miRNA模拟体和pGL3-target载体的细胞中荧光素酶活性将显著低于转染阴性对照和pGL3-target载体的细胞；而转染miRNA模拟体和pGL3-mut-target载体的细胞中荧光素酶活性与转染阴性对照和pGL3-mut-target载体的细胞相比，无明显差异。通过双荧光素酶报告基因实验，验证miRNA与靶基因3'UTR的靶向调控关系。3.2.9数据分析方法运用统计学方法对实验数据进行分析，以揭示实验结果的显著性差异和相关性。对于差异表达基因和miRNA的表达量数据，使用DESeq2软件进行差异显著性检验，计算出每个基因或miRNA在感染组和对照组之间的差异倍数（FoldChange）和校正后的P值（FDR）。设定|log₂(FoldChange)|≥1且FDR＜0.05为差异表达的筛选标准，筛选出在鳗弧菌感染后表达发生显著变化的基因和miRNA。对实时荧光定量PCR验证的数据，采用SPSS软件进行统计学分析。使用单因素方差分析（One-wayANOVA）比较不同处理组之间基因相对表达量的差异，当P＜0.05时，认为差异具有统计学意义。如果存在多个处理组，进一步使用LSD（Least-SignificantDifference）法进行多重比较，以确定具体哪些组之间存在显著差异。在研究miRNA四、实验结果4.1转录组测序与基因表达分析4.1.1测序数据质量评估利用IlluminaHiSeq测序平台对构建好的cDNA文库进行高通量测序，获得了大量的原始数据。原始数据以FASTQ格式存储，包含了测序序列（reads）以及对应的质量分数信息。通过FastQC软件对原始数据进行质量评估，结果显示，原始数据的碱基质量分布较为均匀，大部分碱基的质量分数在30以上，表明测序质量较高。在碱基组成方面，A、T、C、G四种碱基的含量相对均衡，符合正常的测序数据特征。同时，GC含量约为[X]%，处于合理范围内，进一步验证了测序数据的可靠性。使用Trimmomatic软件对原始数据进行预处理，去除低质量的reads、接头序列以及含有过多N（未知碱基）的reads。经过预处理后，得到了高质量的cleanreads，其数据量为[X]条，占原始数据的[X]%。Cleanreads的质量评估结果显示，其碱基质量分数进一步提高，平均质量分数达到了[X]以上，为后续的数据分析提供了可靠的数据基础。4.1.2序列组装与基因注释利用TopHat软件将cleanreads比对到鮸鱼的参考基因组上，确定每个read在基因组上的位置。比对结果显示，cleanreads与参考基因组的比对率为[X]%，表明大部分reads能够准确地比对到参考基因组上。使用Cufflinks软件根据比对结果进行序列组装，共得到了[X]条unigenes。这些unigenes的长度分布范围较广，从[X]bp到[X]bp不等，平均长度为[X]bp。其中，长度大于1000bp的unigenes有[X]条，占总unigenes的[X]%，这些较长的unigenes可能包含了完整的基因编码序列，对于后续的基因功能研究具有重要意义。将组装得到的unigenes在多个数据库中进行比对注释，包括NT、NR、Swiss-Prot和KEGG等数据库。通过BLAST相似性比对方法，共有[X]条unigenes得到了注释，注释率为[X]%。在GO功能注释中，有[X]个unigenes获得了GO注释，并被分类到生物过程、细胞组成和分子功能等51个功能类别中。例如，在生物过程类别中，参与细胞代谢过程的基因有[X]个，参与免疫应答过程的基因有[X]个；在细胞组成类别中，与细胞膜相关的基因有[X]个，与细胞核相关的基因有[X]个；在分子功能类别中，具有酶活性的基因有[X]个，具有转录因子活性的基因有[X]个。在COG功能分类中，有[X]个unigenes被归纳到25个功能分类中，进一步揭示了unigenes的功能分布情况。在KEGG注释中，共有[X]条unigenes获得了KEGG注释，共参与123条代谢通路中，这些代谢通路涉及到细胞信号传导、能量代谢、物质合成与代谢等多个方面，为深入了解鮸鱼的生物学过程提供了重要线索。4.1.3差异表达基因筛选与鉴定使用DESeq2软件对感染组和对照组的基因表达量数据进行差异表达分析，筛选出在鳗弧菌感染后表达发生显著变化的基因。设定差异表达基因的筛选标准为：|log₂(FoldChange)|≥1且FDR（FalseDiscoveryRate）＜0.05。经过筛选，共得到差异表达基因[X]个，其中上调基因有[X]个，下调基因有[X]个。上调基因的表达量在感染组中显著高于对照组，而下调基因的表达量在感染组中显著低于对照组。对差异表达基因进行进一步分析，发现其中一些基因与免疫相关。例如，Toll样受体1（TLR1）基因在鳗弧菌感染后表达显著上调，其log₂(FoldChange)值为[X]，FDR值为[X]。TLR1是Toll/IL-1R家族的重要成员，在鱼类的先天性免疫中发挥着关键作用，能够识别病原体相关分子模式，激活下游信号通路，引发免疫应答。白细胞介素-1受体1（IL-1R1）基因的表达也显著上调，log₂(FoldChange)值为[X]，FDR值为[X]。IL-1R1参与白细胞介素-1介导的信号传导通路，在炎症反应和免疫调节中具有重要作用。此外，还筛选到一些与免疫细胞功能相关的基因，如巨噬细胞迁移抑制因子（MIF）基因，其在感染后表达上调，可能参与了巨噬细胞的活化和迁移过程，增强了免疫细胞对病原体的清除能力。这些与免疫相关的关键差异表达基因的筛选，为进一步研究鮸鱼在鳗弧菌感染后的免疫应答机制提供了重要的基因靶点。4.2差异表达基因的功能富集分析4.2.1GO功能富集分析结果利用DAVID数据库对筛选出的差异表达基因进行GO（GeneOntology）功能富集分析，以了解这些基因在生物学过程、细胞组分和分子功能等方面的富集情况。分析结果显示，在生物学过程类别中，差异表达基因显著富集在免疫应答（immuneresponse）、炎症反应（inflammatoryresponse）、细胞对细菌的防御反应（cellulardefenseresponsetobacterium）等条目。其中，参与免疫应答过程的基因有[X]个，占差异表达基因总数的[X]%；参与炎症反应的基因有[X]个，占比为[X]%；参与细胞对细菌防御反应的基因有[X]个，占比[X]%。这表明在鳗弧菌感染后，鮸鱼的免疫应答和炎症反应相关的生物学过程被显著激活，体现了鮸鱼免疫系统对病原菌入侵的积极响应。在细胞组分类别中，差异表达基因主要富集在细胞膜（cellmembrane）、细胞外基质（extracellularmatrix）、溶酶体（lysosome）等条目。与细胞膜相关的基因有[X]个，占比[X]%，细胞膜是细胞与外界环境进行物质交换和信号传递的重要界面，其相关基因的差异表达可能影响细胞对病原菌的识别和摄取；与细胞外基质相关的基因有[X]个，占比[X]%，细胞外基质在维持细胞形态、调节细胞行为等方面发挥重要作用，其相关基因的变化可能与免疫细胞的迁移和组织修复有关；与溶酶体相关的基因有[X]个，占比[X]%，溶酶体含有多种水解酶，参与细胞内的物质降解和病原体清除，其相关基因的差异表达可能影响免疫细胞对病原菌的消化和清除能力。在分子功能类别中，差异表达基因显著富集在蛋白结合（proteinbinding）、酶活性（enzymeactivity）、细胞因子活性（cytokineactivity）等条目。具有蛋白结合功能的基因有[X]个，占比[X]%，蛋白结合是许多生物学过程的基础，如信号传导、蛋白质相互作用等，这些基因的差异表达可能影响免疫相关信号通路的传递；具有酶活性的基因有[X]个，占比[X]%，酶在生物化学反应中起催化作用，参与免疫应答和炎症反应中的多种代谢过程；具有细胞因子活性的基因有[X]个，占比[X]%，细胞因子在免疫调节和炎症反应中发挥关键作用，能够调节免疫细胞的活性和功能，促进免疫细胞的增殖、分化和迁移。4.2.2KEGG通路富集分析结果对差异表达基因进行KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析，以揭示这些基因参与的主要信号通路。分析结果显示，差异表达基因显著富集在多条与免疫相关的信号通路中，如Toll/IL-1R信号通路（Toll/IL-1Rsignalingpathway）、NF-κB信号通路（NF-κBsignalingpathway）、MAPK信号通路（MAPKsignalingpathway）等。在Toll/IL-1R信号通路中，有[X]个差异表达基因参与，如Toll样受体1（TLR1）、白细胞介素-1受体1（IL-1R1）等基因在该通路中表达显著上调。Toll/IL-1R信号通路是鱼类先天性免疫的关键信号通路，能够识别病原体相关分子模式，激活下游信号转导，引发免疫应答。当鳗弧菌感染鮸鱼后，TLR1和IL-1R1等受体识别鳗弧菌表面的病原体相关分子模式，通过MyD88依赖和MyD88非依赖两条信号转导途径，激活NF-κB等转录因子，上调炎症细胞因子和趋化因子的表达，招募免疫细胞到感染部位，增强吞噬作用和杀菌能力，从而有效清除病原体。在NF-κB信号通路中，有[X]个差异表达基因参与。NF-κB是一种重要的转录因子，在免疫应答和炎症反应中发挥核心调控作用。在正常情况下，NF-κB与抑制蛋白IκB结合，以无活性的形式存在于细胞质中。当细胞受到病原体刺激时，IκB被磷酸化并降解，释放NF-κB，使其进入细胞核，与相关基因的启动子区域结合，启动炎症细胞因子、趋化因子等免疫相关基因的转录，促进炎症反应和免疫应答。在鳗弧菌感染鮸鱼后，NF-κB信号通路被激活，相关差异表达基因的变化可能影响NF-κB的激活和调控，进而影响免疫应答的强度和进程。在MAPK信号通路中，也有[X]个差异表达基因参与。MAPK信号通路在细胞增殖、分化、凋亡和应激反应等过程中发挥重要作用，参与免疫细胞的活化和细胞因子的产生。当免疫细胞受到病原体刺激时，MAPK信号通路被激活，通过一系列激酶级联反应，激活下游的转录因子，调节免疫相关基因的表达。在鳗弧菌感染鮸鱼后，MAPK信号通路的激活可能导致免疫细胞的活化和功能增强，促进免疫应答的发生。4.3免疫相关差异表达基因的miRNA调控分析4.3.1鮸鱼IL-1RI的miRNA调控通过生物信息学分析，筛选出了可能与IL-1RI3'UTR区域结合的miRNA，其中miR-192被重点关注。对miR-192的特征进行深入分析，发现其成熟序列长度为22个核苷酸，具有典型的miRNA二级结构，呈现出茎环状。通过对miR-192的启动子区域分析，预测出其可以被转录因子AP-1所调控，这表明miR-192的表达可能受到AP-1的精确调控，进而影响其对靶基因的作用。为了确定miR-192与IL-1RI3'UTR的作用位点，构建了含有IL-1RI3'UTR野生型序列的荧光素酶报告载体（pGL3-IL-1RI-3'UTR-WT）以及将miR-192与IL-1RI3'UTR结合的种子序列区域进行突变后的荧光素酶报告载体（pGL3-IL-1RI-3'UTR-MUT）。将miR-192模拟体（mimics）或阴性对照（NC）分别与上述两种报告载体共转染至HEK293T细胞中。双荧光素酶报告基因实验结果显示，转染miR-192模拟体和pGL3-IL-1RI-3'UTR-WT载体的细胞中，荧光素酶活性显著低于转染阴性对照和pGL3-IL-1RI-3'UTR-WT载体的细胞，这表明miR-192能够与IL-1RI3'UTR野生型序列相互作用，抑制荧光素酶的表达。而转染miR-192模拟体和pGL3-IL-1RI-3'UTR-MUT载体的细胞中，荧光素酶活性与转染阴性对照和pGL3-IL-1RI-3'UTR-MUT载体的细胞相比，无明显差异，这证实miR-192对IL-1R的抑制效应是严格依赖于“Seedregion”与IL-1RImRNA之间的相互作用。进一步对miR-192与IL-1RI在鳗弧菌感染过程中的表达变化进行分析。荧光定量分析结果显示，在鳗弧菌感染的鮸鱼脾脏组织和LPS刺激的鮸鱼巨噬细胞中，IL-1RI的表达量显著升高，而miR-192的表达量变化趋势则与IL-1RI相反，呈现出显著下降的趋势。这一结果进一步证实了miR-192对IL-1RI表达的调节作用，即miR-192通过与IL-1RI3'UTR的相互作用，在转录后水平抑制IL-1RI的表达，从而在鮸鱼免疫及炎症反应中发挥重要的调控作用。4.3.2鮸鱼TLR1的miRNA调控通过生物学信息学方法及双荧光素酶报告基因实验，筛选出miR-217可靶作用于TLR1的3'UTR。miR-217的成熟序列长度为21个核苷酸，其种子序列与TLR13'UTR的特定区域具有高度的互补性。构建含有TLR13'UTR野生型序列的荧光素酶报告载体（pGL3-TLR1-3'UTR-WT）以及将miR-217与TLR13'UTR结合的种子序列区域进行突变后的荧光素酶报告载体（pGL3-TLR1-3'UTR-MUT）。将miR-217模拟体（mimics）或阴性对照（NC）分别与这两种报告载体共转染至HEK293T细胞中。双荧光素酶报告基因实验结果表明，miR-217可以抑制TLR13'UTR野生型报告质粒的荧光活性，而对突变型报告质粒的荧光活性无明显影响，这验证了miR-217可以直接靶作用于TLR13'UTR的“Seedregion”，通过与TLR13'UTR的互补配对，抑制荧光素酶的表达，进而影响TLR1的表达水平。为了进一步探究miR-217前体对TLR1表达的影响，构建了miR-217前体质粒。将miR-217前体质粒或阴性对照质粒与pGL3-TLR1-3'UTR-WT载体共转染至HEK293T细胞中。双荧光素酶报告基因实验发现，miR-217前体具有类似于miR-217模拟体的活性，也可以靶作用于TLR13'UTR而抑制荧光活性。这表明miR-217前体在细胞内能够被加工成成熟的miR-217，进而发挥对TLR1表达的调控作用。对miR-217在LPS刺激的脾脏组织和LPS刺激的鮸鱼巨噬细胞中的表达量进行荧光定量分析。结果显示，miR-217的表达量在LPS刺激后显著上调，表明miR-217可能参与调节机体对LPS刺激的反应。结合之前的实验结果，miR-217通过与TLR13'UTR的相互作用，在转录后水平对TLR1的表达进行调控，在鮸鱼应对革兰氏阴性菌感染的免疫应答过程中发挥重要的作用。五、分析与讨论5.1鳗弧菌感染对鮸鱼基因表达的影响通过对鳗弧菌感染鮸鱼后的转录组测序分析，我们发现了大量差异表达基因，这些基因在免疫应答、代谢、细胞凋亡等多个方面发生了显著变化，对鮸鱼抵抗鳗弧菌感染具有重要意义。在免疫应答方面，许多与免疫相关的基因表达发生改变。Toll样受体1（TLR1）基因表达显著上调，TLR1作为Toll/IL-1R家族的重要成员，在鱼类先天性免疫中发挥关键作用。它能够识别病原体相关分子模式，激活下游信号通路，如MyD88依赖和MyD88非依赖信号通路，最终导致核因子κB（NF-κB）等转录因子的激活，启动炎症细胞因子（如白细胞介素-1β、白细胞介素-6、肿瘤坏死因子α等）、趋化因子等免疫相关基因的转录。TLR1基因的上调表达，表明鮸鱼在鳗弧菌感染后，能够通过激活Toll/IL-1R信号通路，增强免疫应答，以抵御病原菌的入侵。白细胞介素-1受体1（IL-1R1）基因表达也显著上调，IL-1R1参与白细胞介素-1介导的信号传导通路，在炎症反应和免疫调节中具有重要作用。其表达上调可能促进了炎症反应的发生，吸引更多免疫细胞到感染部位，增强免疫防御能力。在代谢方面，一些参与能量代谢和物质合成代谢的基因表达发生变化。与糖代谢相关的基因，如己糖激酶基因表达上调，这可能是为了满足鮸鱼在免疫应答过程中对能量的需求。在免疫应答过程中，免疫细胞的活化、增殖以及炎症细胞因子的合成和分泌等都需要消耗大量能量，上调糖代谢相关基因的表达，有助于提高糖的分解代谢，产生更多的能量，为免疫反应提供支持。此外，一些参与脂质代谢和蛋白质合成代谢的基因表达也发生改变，这可能与免疫细胞的功能和结构变化有关。免疫细胞在激活后，需要合成更多的蛋白质和脂质，以构建细胞膜、细胞器等结构，同时满足细胞内信号传导和免疫分子合成的需求。细胞凋亡相关基因的表达变化在鮸鱼抵抗鳗弧菌感染中也具有重要意义。一些促凋亡基因，如Bax基因表达上调，而抗凋亡基因，如Bcl-2基因表达下调。细胞凋亡是一种程序性细胞死亡过程，在免疫应答中，适当的细胞凋亡可以清除被病原体感染的细胞，防止病原体在细胞内大量繁殖，同时也有助于调节免疫细胞的数量和功能，维持免疫平衡。Bax基因表达上调和Bcl-2基因表达下调，可能导致细胞凋亡的增加，从而促进被鳗弧菌感染细胞的清除，减少病原体的扩散，保护机体免受进一步的损害。这些差异表达基因在免疫应答、代谢、细胞凋亡等方面的变化，是鮸鱼在鳗弧菌感染后启动的一系列自我保护机制。它们相互协作，共同调节鮸鱼的生理状态，增强其抵抗鳗弧菌感染的能力。通过激活免疫应答相关基因，鮸鱼能够识别和清除病原体；调节代谢相关基因，为免疫反应提供充足的能量和物质支持；调控细胞凋亡相关基因，维持细胞和组织的正常功能和结构。深入研究这些基因的变化及其相互作用机制，有助于我们全面了解鮸鱼的免疫防御机制，为开发有效的病害防治策略提供理论依据。5.2miRNA在鮸鱼免疫应答中的调控作用miRNA在鮸鱼免疫应答中发挥着至关重要的调控作用，通过对免疫相关基因的精细调控，维持着鮸鱼免疫平衡。研究发现，miR-192对IL-1RI基因的表达具有显著的调控作用。在鳗弧菌感染的鮸鱼脾脏组织和LPS刺激的鮸鱼巨噬细胞中，IL-1RI的表达量显著升高，而miR-192的表达量则显著下降。这表明在免疫应答过程中，miR-192可能作为一种负调控因子，通过与IL-1RI3'UTR的相互作用，抑制IL-1RI的表达。IL-1RI是白细胞介素-1介导的信号传导通路中的关键受体，其表达的上调会激活下游的炎症反应相关信号通路。而miR-192的下调则解除了对IL-1RI的抑制作用，使得IL-1RI能够正常发挥其在免疫应答中的作用，促进炎症细胞因子的表达和释放，增强免疫细胞的活性，从而有效抵御鳗弧菌的感染。miR-217对TLR1基因的表达也具有重要的调控作用。在LPS刺激的脾脏组织和LPS刺激的鮸鱼巨噬细胞中，miR-217的表达量显著上调。通过双荧光素酶报告基因实验证实，miR-217可以直接靶作用于TLR13'UTR的“Seedregion”，抑制荧光素酶的表达，进而抑制TLR1的表达。TLR1是Toll/IL-1R家族的重要成员，在识别病原体相关分子模式、激活先天性免疫应答中发挥关键作用。当miR-217表达上调时，它会与TLR13'UTR结合，抑制TLR1的表达，从而在一定程度上调节先天性免疫应答的强度。这可能是鮸鱼机体的一种自我保护机制，避免免疫应答过度激活，导致炎症反应失控，对机体造成损伤。miR-192和miR-217等miRNA对免疫相关基因的调控，在调节鮸鱼免疫平衡中起着关键作用。它们通过对免疫相关基因表达的精确调控，使得免疫应答在适当的范围内进行。当病原体入侵时，miRNA能够及时调整免疫相关基因的表达，增强免疫应答，抵御病原体的感染。而在免疫应答后期，miRNA又能够抑制免疫相关基因的过度表达，避免免疫反应对机体造成损伤，从而维持免疫平衡。这些研究结果与其他鱼类免疫相关miRNA的研究具有一定的相似性。在鲈鱼中，研究发现某些miRNA能够通过调控免疫相关基因的表达，参与鲈鱼对细菌感染的免疫应答过程。在鲫鱼中，miRNA对免疫细胞的分化和功能调节也起着重要作用。这表明miRNA在鱼类免疫调控中具有普