探寻牛流行热病毒感染宿主细胞中差异表达miRNA：筛选、鉴定与功能解析

探寻牛流行热病毒感染宿主细胞中差异表达miRNA：筛选、鉴定与功能解析一、引言1.1研究背景与意义1.1.1牛流行热病毒概述牛流行热病毒（BovineEphemeralFeverVirus，BEFV），隶属弹状病毒科暂时热病毒属，是引发牛流行热的病原体，该病毒呈子弹状或圆锥形，成熟病毒粒子长130-220nm、宽60-70nm，具有囊膜结构，基因组为单股RNA。病毒粒子由核衣壳和外层的脂质囊膜构成，囊膜表面存在纤突，这些纤突在病毒感染宿主细胞过程中发挥关键作用，它们能够与宿主细胞表面的特异性受体结合，从而介导病毒进入细胞内。牛流行热病毒主要感染黄牛、奶牛和水牛等偶蹄类动物，其中3-5岁的牛易感性较高。病牛是主要传染源，处于高热期的病牛，其血液、呼吸道分泌物以及粪便中均含有大量病毒。病毒的传播途径主要包括吸血昆虫叮咬和呼吸道感染。在自然条件下，库蚊、伊蚊、牛虻等吸血昆虫是重要的传播媒介，它们在叮咬病牛后，病毒会在其体内繁殖，当再次叮咬健康牛时，就会将病毒传播给健康牛。此外，健康牛吸入含有病毒的飞沫，也可能感染牛流行热病毒。牛流行热在全球范围内广泛分布，热带和亚热带地区是其高发区域，温带地区也时有发生。该病具有明显的季节性，多在蚊虫繁殖活跃的夏季和秋季流行。牛流行热的爆发往往具有周期性，通常3-5年流行一次。每次流行都会给养牛业带来沉重打击，造成巨大的经济损失。例如，病牛会出现高热、呼吸急促、跛行、食欲不振等症状，导致奶牛产奶量大幅下降，乳品质变差；役用牛的工作能力降低；怀孕母牛可能发生流产，新生犊牛的死亡率增加。据相关统计，在牛流行热流行期间，牛群的发病率可达14%-39.89%，虽无继发疾病时病死率为0.73%-1.80%，但由于发病牛数量众多，总体经济损失依然十分可观。1.1.2miRNA在病毒感染中的作用miRNA（microRNA）是一类内源性非编码小分子RNA，长度约为18-24个核苷酸。它在生物体内参与多种重要的生理和病理过程，尤其是在病毒感染宿主的过程中，miRNA扮演着至关重要的角色。在病毒感染宿主细胞时，miRNA可以通过多种方式调控病毒的复制和传播。一方面，宿主细胞可以产生特异性的miRNA来识别并结合病毒的核酸序列，抑制病毒基因的表达，从而限制病毒的复制。例如，在乙型肝炎病毒（HBV）感染中，宿主细胞产生的某些miRNA能够与HBV的mRNA互补配对，通过RNA诱导的沉默复合体（RISC）降解病毒mRNA，或者抑制其翻译过程，进而减少病毒的产生。另一方面，病毒也可以利用宿主细胞的miRNA来促进自身的复制和生存。一些病毒会编码产生与宿主miRNA相似的序列，干扰宿主细胞正常的miRNA调控网络，为病毒的感染和繁殖创造有利条件。比如，Epstein-Barr病毒（EBV）编码的miRNA可以靶向宿主细胞的多个基因，抑制宿主的免疫应答，有利于病毒的潜伏感染。此外，miRNA还能够调节宿主细胞在病毒感染后的免疫反应。它可以通过调控免疫细胞的分化、活化以及细胞因子的分泌，影响宿主的抗病毒免疫能力。例如，miR-155在病毒感染时能够调节巨噬细胞和T细胞的功能，增强机体的抗病毒免疫应答。然而，miRNA的异常表达也可能导致免疫失衡，引发过度的炎症反应，对宿主造成损伤。1.1.3研究目的和意义本研究旨在筛选和鉴定牛流行热病毒感染宿主细胞后差异表达的miRNA，通过深入分析这些miRNA的功能，揭示它们在牛流行热病毒感染过程中的作用机制。这一研究具有重要的理论和实际意义。从理论层面来看，目前对于牛流行热病毒与宿主细胞之间的相互作用机制，尤其是miRNA在其中的调控作用，了解还十分有限。本研究通过系统地筛选和鉴定差异表达的miRNA，能够填补这一领域的空白，丰富我们对牛流行热病毒致病机制的认识，为进一步研究病毒与宿主的相互作用提供新的视角和理论基础。在实际应用方面，牛流行热给养牛业带来了严重的经济损失，迫切需要有效的防控措施。明确差异表达miRNA的功能，有望为牛流行热的诊断、治疗和预防提供新的靶点和策略。例如，可以基于差异表达的miRNA开发新型的诊断方法，提高诊断的准确性和时效性；通过调控这些miRNA的表达，研发新的抗病毒药物，为牛流行热的治疗提供新的手段；此外，还可以将相关miRNA作为疫苗设计的靶点，增强疫苗的免疫效果，提高牛群对牛流行热病毒的抵抗力。1.2国内外研究现状1.2.1牛流行热病毒感染机制研究进展牛流行热病毒入侵宿主细胞是感染的起始关键步骤。病毒表面的纤突蛋白在这一过程中发挥着核心作用，它能够特异性地识别并结合宿主细胞表面的受体，目前研究发现，牛流行热病毒可能利用宿主细胞表面的硫酸乙酰肝素等作为受体结合位点，通过受体介导的内吞作用进入细胞。当病毒与受体结合后，病毒囊膜与细胞膜发生融合，将病毒基因组释放到宿主细胞内，从而开启病毒的感染进程。进入细胞后，牛流行热病毒以自身的单股RNA为模板，在病毒编码的RNA依赖性RNA聚合酶的作用下，进行基因组的复制和转录。病毒首先转录出多个mRNA，这些mRNA分别编码病毒的结构蛋白和非结构蛋白。结构蛋白如核蛋白（N）、磷蛋白（P）、基质蛋白（M）、糖蛋白（G）等，它们参与病毒粒子的组装；非结构蛋白则在病毒的复制、转录以及逃避宿主免疫监视等过程中发挥重要作用。在病毒复制过程中，会利用宿主细胞的物质和能量代谢系统，如宿主细胞的核糖体、ATP等，来合成病毒自身的蛋白质和核酸。随着病毒蛋白和核酸的大量合成，它们在宿主细胞内组装成新的病毒粒子，通过出芽的方式从宿主细胞释放，继续感染其他细胞。牛流行热病毒感染对宿主细胞产生多方面的影响。在细胞形态上，感染细胞会出现肿胀、变形等变化，这是由于病毒感染导致细胞内的代谢紊乱和结构损伤。在细胞生理功能方面，病毒感染会抑制宿主细胞的正常基因表达，干扰细胞的信号传导通路。例如，病毒感染可能抑制宿主细胞中与免疫应答相关基因的表达，从而削弱宿主细胞的抗病毒免疫能力。此外，病毒感染还会引发细胞的应激反应，如内质网应激等，若应激反应持续过度，可能导致细胞凋亡。细胞凋亡是宿主细胞抵御病毒感染的一种防御机制，但同时也可能导致组织器官的损伤，进而引发牛流行热的各种临床症状。1.2.2miRNA在病毒感染中的研究现状在众多病毒感染的研究中，miRNA展现出复杂且关键的调控作用。以乙型肝炎病毒（HBV）感染为例，宿主细胞内的miR-122对HBV的复制具有重要调控作用。miR-122能够与HBV的mRNA结合，促进其稳定性，从而有利于HBV的复制。研究发现，通过抑制miR-122的表达，可以显著降低HBV的复制水平。在丙型肝炎病毒（HCV）感染中，miR-122同样发挥着重要作用，它不仅可以促进HCV的复制，还参与调控HCV感染相关的脂质代谢紊乱。此外，在艾滋病病毒（HIV）感染过程中，宿主细胞的miRNA也参与调节HIV的潜伏和激活。例如，miR-28、miR-125b等可以靶向HIV的mRNA，抑制其表达，从而影响HIV的感染进程。这些研究成果与牛流行热病毒感染研究存在一定的关联和借鉴意义。一方面，其他病毒感染中miRNA的作用机制提示我们，牛流行热病毒感染宿主细胞后，宿主细胞可能也会产生特异性的miRNA来应对病毒感染，这些miRNA可能通过与牛流行热病毒的核酸或相关蛋白相互作用，影响病毒的复制、转录和组装等过程。另一方面，病毒也可能利用宿主细胞的miRNA来促进自身的感染和生存，这为研究牛流行热病毒与宿主细胞之间的相互作用提供了新的思路。通过对其他病毒感染中miRNA研究的深入分析，可以为筛选和鉴定牛流行热病毒感染宿主细胞后差异表达的miRNA提供方法和技术上的参考，有助于更高效地开展相关研究。1.2.3研究存在的问题与不足当前在牛流行热病毒与miRNA关联研究中存在诸多空白和待解决问题。从研究广度来看，目前对于牛流行热病毒感染宿主细胞后，哪些miRNA会发生差异表达，以及这些miRNA在病毒感染的不同阶段（如入侵、复制、释放等）的动态变化情况，缺乏全面系统的研究。这使得我们难以从整体上把握miRNA在牛流行热病毒感染过程中的调控网络。在研究深度方面，虽然已经知晓miRNA在病毒感染中具有重要作用，但对于牛流行热病毒感染宿主细胞后差异表达的miRNA具体的作用靶点和分子机制，研究还十分有限。例如，不清楚这些miRNA是如何与牛流行热病毒的基因组或宿主细胞的相关基因相互作用，进而影响病毒的感染和宿主细胞的免疫应答。此外，目前对于miRNA在牛流行热病毒感染的临床应用研究也相对滞后，如何将miRNA的研究成果转化为有效的诊断、治疗和预防手段，还需要进一步探索。例如，如何基于miRNA开发高灵敏度和特异性的诊断方法，以及如何利用miRNA设计安全有效的抗病毒药物和疫苗等，都是亟待解决的问题。二、材料与方法2.1实验材料2.1.1病毒株与宿主细胞本实验选用的牛流行热病毒株为[具体毒株名称]，其分离自[分离地点]的发病牛群。该毒株经过多次传代和鉴定，具有典型的牛流行热病毒特性，如病毒粒子呈子弹状，基因组为单股RNA，能够在宿主细胞内高效复制并引发细胞病变效应。宿主细胞为[宿主细胞类型]，该细胞系来源于牛的[组织名称]，具有良好的生长特性和对牛流行热病毒的易感性。在实验前，将宿主细胞置于含有10%胎牛血清（FetalBovineSerum，FBS）、100U/mL青霉素和100μg/mL链霉素的[细胞培养液名称]中，在37℃、5%CO₂的培养箱中进行培养。每隔2-3天进行一次传代，以确保细胞处于对数生长期，用于后续的病毒感染实验。2.1.2主要试剂与仪器实验所需的主要试剂包括：RNA提取试剂TRIzol（Invitrogen公司），该试剂能够有效裂解细胞，提取高质量的总RNA；反转录试剂PrimeScriptRTMasterMix（TaKaRa公司），可将提取的RNA反转录为cDNA，用于后续的PCR扩增和miRNA分析；SYBRGreenPCRMasterMix（Roche公司），用于实时荧光定量PCR（qRT-PCR），以检测差异表达miRNA的相对表达量；RNA酶抑制剂（RNaseInhibitor），可防止RNA在操作过程中被降解。此外，还包括氯仿、异丙醇、无水乙醇等常规试剂，用于RNA提取过程中的纯化步骤。主要实验仪器有：高速冷冻离心机（Eppendorf公司），用于细胞和病毒的离心分离；PCR仪（Bio-Rad公司），进行PCR扩增反应；实时荧光定量PCR仪（ABI公司），实现对miRNA表达量的精确检测；核酸蛋白测定仪（NanoDrop公司），用于测定RNA和cDNA的浓度和纯度；凝胶成像系统（Bio-Rad公司），观察和分析PCR产物的电泳结果。2.2实验方法2.2.1牛流行热病毒感染宿主细胞模型的建立将处于对数生长期的[宿主细胞类型]以[细胞接种密度]接种于6孔细胞培养板中，每孔加入2mL细胞培养液。将细胞置于37℃、5%CO₂的培养箱中培养24h，待细胞贴壁且融合度达到70%-80%时，进行病毒感染实验。用无血清的细胞培养液将牛流行热病毒[具体毒株名称]稀释至不同的感染复数（MOI），分别为0.1、0.5、1、5、10。弃去6孔板中的细胞培养液，用PBS轻轻洗涤细胞2-3次，以去除残留的血清和杂质。向每孔中加入1mL不同MOI的病毒稀释液，确保病毒均匀覆盖细胞表面。将细胞培养板置于37℃、5%CO₂的培养箱中孵育1-2h，期间每隔15-20min轻轻摇晃培养板，使病毒与细胞充分接触。孵育结束后，弃去病毒液，用PBS再次洗涤细胞2-3次，以去除未吸附的病毒。然后向每孔中加入2mL含有2%胎牛血清的细胞培养液，继续在37℃、5%CO₂的培养箱中培养。在感染后的不同时间点（6h、12h、24h、36h、48h），通过倒置显微镜观察细胞病变效应（CPE），记录细胞的形态变化，如细胞肿胀、变圆、脱落等情况。同时，收集细胞培养上清液，采用实时荧光定量PCR（qRT-PCR）方法检测病毒的核酸拷贝数，以确定病毒在细胞内的复制情况。通过对不同MOI和感染时间下病毒复制和CPE的观察分析，确定最佳的病毒接种量和感染时间。最终确定以MOI=1的病毒接种量感染细胞24h，可建立稳定且有效的牛流行热病毒感染宿主细胞模型。2.2.2miRNA的提取与高通量测序分别收集牛流行热病毒感染24h的宿主细胞和未感染的正常宿主细胞，用预冷的PBS洗涤细胞3次，以去除细胞表面的杂质和培养液。向每孔细胞中加入1mLTRIzol试剂，充分吹打混匀，使细胞裂解。将裂解液转移至无RNA酶的离心管中，室温静置5min，使核酸蛋白复合物充分解离。向离心管中加入0.2mL氯仿，剧烈振荡15s，室温静置3min。然后在4℃、12000g条件下离心15min，此时溶液会分为三层，上层为无色透明的水相，含有RNA；中间层为白色的蛋白层；下层为红色的有机相。小心吸取上层水相转移至新的离心管中，避免吸取到中间层的蛋白和下层的有机相。向水相中加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min，使RNA沉淀。在4℃、12000g条件下离心10min，可见管底出现白色的RNA沉淀。弃去上清液，加入1mL75%的乙醇（用DEPC水配制），轻轻洗涤RNA沉淀，以去除残留的杂质和盐分。在4℃、7500g条件下离心5min，弃去上清液。将离心管倒置在滤纸上，室温晾干RNA沉淀5-10min，注意不要让RNA沉淀过于干燥，以免影响后续的溶解。向晾干的RNA沉淀中加入适量的DEPC水，轻轻吹打溶解RNA。使用核酸蛋白测定仪测定RNA的浓度和纯度，要求RNA的A₂₆₀/A₂₈₀比值在1.8-2.0之间，A₂₆₀/A₂₃₀比值大于2.0，以确保提取的RNA质量良好。将提取的RNA保存于-80℃冰箱中备用。将提取的RNA样品送至专业的测序公司进行高通量测序。测序公司首先对RNA样品进行质量检测，确保RNA的完整性和纯度符合测序要求。然后，利用Illumina测序平台进行小RNA文库的构建和测序。具体流程如下：首先，在RNA的3'端和5'端分别连接特异性接头，然后进行逆转录反应，将RNA反转录为cDNA。接着，通过PCR扩增富集cDNA片段，构建小RNA文库。对文库进行质量检测和定量后，在Illumina测序仪上进行测序，得到原始的测序数据。测序公司对原始测序数据进行初步的质量控制，去除低质量的reads、接头序列和污染序列。然后，将高质量的reads与参考基因组进行比对，确定miRNA的表达谱。通过生物信息学分析，识别出已知的miRNA和潜在的新miRNA，并计算它们在感染组和未感染组细胞中的表达量。2.2.3差异表达miRNA的筛选与生物信息学分析采用DESeq2软件对高通量测序得到的miRNA表达数据进行分析，筛选出在牛流行热病毒感染组和未感染组宿主细胞中差异表达的miRNA。筛选标准为：|log₂(fold-change)|≥1且调整后的P值（adj.P.Val）＜0.05。其中，log₂(fold-change)表示感染组与未感染组miRNA表达量的对数比值，反映了miRNA表达的变化倍数；adj.P.Val是通过多重检验校正后的P值，用于控制假阳性率。利用miRanda、TargetScan和PicTar等生物信息学工具，对筛选出的差异表达miRNA进行靶基因预测。这些工具基于miRNA与靶基因mRNA之间的互补配对原则，通过算法预测miRNA可能结合的靶基因。为了提高预测的准确性，取这三个工具预测结果的交集作为最终的靶基因集合。对预测得到的靶基因进行基因本体论（GO）功能富集分析和京都基因与基因组百科全书（KEGG）通路富集分析。GO功能富集分析从生物过程（BiologicalProcess）、细胞组分（CellularComponent）和分子功能（MolecularFunction）三个层面，对靶基因的功能进行注释和富集分析，以了解差异表达miRNA可能参与调控的生物学过程和细胞功能。KEGG通路富集分析则是将靶基因映射到KEGG数据库中的各种信号通路，找出显著富集的信号通路，从而揭示差异表达miRNA在细胞信号传导和代谢途径中的作用机制。使用clusterProfiler等R软件包进行富集分析，设定富集的显著性阈值为P＜0.05。2.2.4miRNA表达验证采用实时荧光定量PCR（qRT-PCR）方法对筛选出的差异表达miRNA进行表达验证。首先，将提取的总RNA按照逆转录试剂盒（PrimeScriptRTMasterMix）的说明书进行逆转录反应，将RNA反转录为cDNA。反应体系为：5×PrimeScriptRTMasterMix2μL，总RNA1μg，RNaseFreedH₂O补足至10μL。反应条件为：37℃15min，85℃5s，4℃保存。根据miRNA的成熟序列，设计特异性的茎环引物用于qRT-PCR扩增。以U6snRNA作为内参基因，其引物序列为已知通用序列。qRT-PCR反应体系为：2×SYBRGreenPCRMasterMix10μL，上下游引物（10μmol/L）各0.5μL，cDNA模板1μL，ddH₂O8μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；最后进行熔解曲线分析，以验证扩增产物的特异性。每个样品设置3个生物学重复和3个技术重复，采用2⁻ΔΔCt法计算miRNA的相对表达量。其中，ΔCt=Ct(miRNA)-Ct(U6)，ΔΔCt=ΔCt(感染组)-ΔCt(未感染组)。通过比较qRT-PCR结果与高通量测序结果中miRNA的表达趋势，验证差异表达miRNA的可靠性。除了qRT-PCR外，还可采用Northernblot方法对部分差异表达miRNA进行验证。将提取的总RNA进行变性聚丙烯酰胺凝胶电泳（PAGE），使RNA按照大小分离。然后将凝胶中的RNA转移至尼龙膜上，通过紫外交联固定RNA。用放射性同位素或地高辛标记的miRNA探针与尼龙膜上的RNA进行杂交，洗去未结合的探针后，通过放射自显影或化学发光检测杂交信号。根据杂交信号的强度判断miRNA的表达水平，与高通量测序和qRT-PCR结果进行对比分析。2.2.5miRNA功能验证为了验证差异表达miRNA的功能，构建miRNA过表达载体和miRNA抑制剂。对于需要过表达的miRNA，合成其前体序列（pre-miRNA），并将其克隆到真核表达载体（如pcDNA3.1）中，构建成miRNA过表达载体。对于需要抑制表达的miRNA，设计并合成相应的反义寡核苷酸（antagomir），即miRNA抑制剂。将宿主细胞以[细胞接种密度]接种于24孔细胞培养板中，每孔加入500μL细胞培养液。培养24h后，待细胞贴壁且融合度达到50%-60%时，进行转染实验。按照脂质体转染试剂（如Lipofectamine3000）的说明书，将miRNA过表达载体或miRNA抑制剂与脂质体混合，形成转染复合物。将转染复合物加入到细胞培养液中，轻轻混匀，继续在37℃、5%CO₂的培养箱中培养。转染6-8h后，更换为含有正常血清的细胞培养液，继续培养。分别设置过表达组（转染miRNA过表达载体）、抑制组（转染miRNA抑制剂）、阴性对照组（转染空载体或阴性对照寡核苷酸）和空白对照组（未转染任何物质）。在转染后的24h、48h和72h，收集细胞，采用qRT-PCR方法检测miRNA的表达水平，以验证转染效果。过表达或敲低miRNA后，检测细胞表型和病毒感染相关指标的变化。观察细胞的形态变化，通过细胞计数法检测细胞增殖能力，采用流式细胞术检测细胞周期和凋亡情况。对于病毒感染相关指标，在转染后的细胞中接种牛流行热病毒，按照建立的感染模型条件进行感染。感染一定时间后，收集细胞培养上清液，采用qRT-PCR方法检测病毒的核酸拷贝数，用ELISA方法检测病毒蛋白的表达水平，以评估miRNA对病毒感染和复制的影响。三、牛流行热病毒感染宿主细胞差异表达miRNA的筛选结果3.1高通量测序数据质量评估3.1.1测序数据的基本统计经过高通量测序，本研究共获得了[X]条原始读长（reads）。在对原始数据进行严格的质量控制后，去除了低质量的reads、接头序列以及污染序列，最终得到了[X]条高质量的有效读长，有效读长占原始读长的比例为[X]%，这一比例表明测序数据的质量较高，能够满足后续分析的需求。对有效读长的长度分布进行分析，结果显示，读长主要集中在18-24nt之间，这与miRNA的长度范围相符，进一步验证了测序数据的可靠性。其中，长度为22nt的读长数量最多，占总有效读长的[X]%。在测序数据的碱基组成方面，对GC含量进行了统计分析。结果表明，GC含量为[X]%，处于正常范围之内。一般来说，GC含量过高或过低都可能影响测序数据的质量和后续分析的准确性，本研究中合适的GC含量说明测序过程较为稳定，数据质量可靠。此外，对测序数据的碱基质量值进行评估，Q30（碱基错误率低于0.1%）的碱基比例达到了[X]%，这意味着大部分碱基的测序质量较高，能够为后续的分析提供准确的基础。3.1.2数据比对与注释将高质量的有效读长与牛的参考基因组进行比对，使用了Bowtie等比对软件，以确保比对结果的准确性。比对结果显示，共有[X]条读长成功比对到参考基因组上，比对率为[X]%。这一比对率表明，大部分测序数据能够与参考基因组进行有效匹配，为后续的miRNA鉴定和分析提供了有力支持。在miRNA的注释方面，通过与miRBase数据库进行比对，鉴定出了已知的miRNA。共注释到[X]种已知的miRNA，其中表达量较高的miRNA包括miR-[X1]、miR-[X2]等。这些已知miRNA在细胞的生长、发育、分化以及免疫应答等过程中发挥着重要作用，它们在牛流行热病毒感染宿主细胞后的表达变化，可能与病毒感染的机制密切相关。除了已知的miRNA，还通过生物信息学分析预测了潜在的新miRNA。利用miRDeep2等软件，根据miRNA的前体结构、成熟序列以及表达特征等信息，预测出了[X]个潜在的新miRNA。对这些新miRNA的二级结构进行分析，发现它们具有典型的发夹结构，这是miRNA的重要特征之一。虽然这些新miRNA的功能尚未明确，但它们在牛流行热病毒感染宿主细胞过程中的差异表达，暗示着它们可能参与了病毒与宿主细胞之间的相互作用，为进一步研究病毒感染机制提供了新的线索。3.2差异表达miRNA的筛选3.2.1差异表达分析结果通过对牛流行热病毒感染组和未感染组宿主细胞的miRNA表达数据进行DESeq2分析，筛选出了差异表达的miRNA。以|log₂(fold-change)|≥1且adj.P.Val＜0.05为筛选标准，共得到[X]个差异表达的miRNA，其中上调表达的miRNA有[X]个，下调表达的miRNA有[X]个。为了直观地展示差异表达miRNA的分布情况，绘制了火山图（图1）。在火山图中，横坐标表示log₂(fold-change)，反映了miRNA在感染组和未感染组之间的表达倍数变化；纵坐标表示-log₁₀(adj.P.Val)，代表差异表达的显著性水平。图中的每个点代表一个miRNA，红色的点表示上调表达显著的miRNA，绿色的点表示下调表达显著的miRNA，黑色的点表示差异表达不显著的miRNA。从火山图中可以清晰地看出，上调和下调表达显著的miRNA在图中呈现出明显的聚集分布，表明牛流行热病毒感染宿主细胞后，确实导致了部分miRNA的表达发生显著变化。[此处插入火山图，图题：牛流行热病毒感染宿主细胞差异表达miRNA火山图，图注：红色点表示上调表达显著的miRNA，绿色点表示下调表达显著的miRNA，黑色点表示差异表达不显著的miRNA]进一步绘制了差异表达miRNA的热图（图2），以展示这些miRNA在感染组和未感染组中的表达模式。热图中，每一行代表一个差异表达的miRNA，每一列代表一个样本（感染组或未感染组）。颜色的深浅表示miRNA表达量的高低，红色表示高表达，蓝色表示低表达。通过热图可以直观地看到，上调表达的miRNA在感染组中的表达量明显高于未感染组，而下调表达的miRNA在感染组中的表达量则明显低于未感染组。此外，还可以发现不同的miRNA在感染组和未感染组中的表达模式存在一定的差异，这暗示着它们可能在牛流行热病毒感染过程中发挥着不同的作用。[此处插入热图，图题：牛流行热病毒感染宿主细胞差异表达miRNA热图，图注：红色表示高表达，蓝色表示低表达]表1列出了上调和下调表达显著的前10个miRNA及其log₂(fold-change)和adj.P.Val值。其中，上调表达最显著的miRNA为miR-[X1]，其log₂(fold-change)达到了[X]，adj.P.Val＜0.001，表明在牛流行热病毒感染宿主细胞后，miR-[X1]的表达量显著增加。下调表达最显著的miRNA为miR-[X2]，log₂(fold-change)为-[X]，adj.P.Val＜0.001，说明miR-[X2]在感染后表达量急剧下降。这些差异表达显著的miRNA是后续研究的重点对象，它们可能在牛流行热病毒感染的过程中发挥着关键的调控作用。表1：上调和下调表达显著的前10个miRNAmiRNA名称log₂(fold-change)adj.P.Val表达趋势miR-[X1][X]＜0.001上调miR-[X3][X]＜0.001上调miR-[X4][X]＜0.001上调miR-[X5][X]＜0.001上调miR-[X6][X]＜0.001上调miR-[X7][X]＜0.001上调miR-[X8][X]＜0.001上调miR-[X9][X]＜0.001上调miR-[X10][X]＜0.001上调miR-[X2]-[X]＜0.001下调miR-[X11]-[X]＜0.001下调miR-[X12]-[X]＜0.001下调miR-[X13]-[X]＜0.001下调miR-[X14]-[X]＜0.001下调miR-[X15]-[X]＜0.001下调miR-[X16]-[X]＜0.001下调miR-[X17]-[X]＜0.001下调miR-[X18]-[X]＜0.001下调miR-[X19]-[X]＜0.001下调3.2.2差异表达miRNA的聚类分析为了深入了解差异表达miRNA的表达模式和相互关系，对筛选出的[X]个差异表达miRNA进行了聚类分析。聚类分析采用层次聚类算法，通过计算miRNA之间的表达相似性，将表达模式相近的miRNA聚为一类。聚类分析结果以聚类树状图（图3）的形式展示。在聚类树状图中，横坐标表示不同的样本（感染组和未感染组），纵坐标表示差异表达的miRNA。树状图中的分支表示miRNA的聚类关系，分支越短，说明miRNA之间的表达模式越相似。从聚类树状图中可以看出，差异表达的miRNA被分为了多个不同的簇，每个簇内的miRNA具有相似的表达模式。其中，一些簇中的miRNA在感染组中表达上调，而另一些簇中的miRNA在感染组中表达下调。这表明，牛流行热病毒感染宿主细胞后，不同的miRNA可能通过不同的表达模式参与病毒感染的调控过程。[此处插入聚类树状图，图题：牛流行热病毒感染宿主细胞差异表达miRNA聚类树状图，图注：横坐标表示样本，纵坐标表示差异表达的miRNA，分支表示miRNA的聚类关系]进一步对每个簇内的miRNA进行分析，发现同一簇内的miRNA可能参与相似的生物学过程或信号通路。例如，在一个上调表达的簇中，包含了miR-[X1]、miR-[X3]和miR-[X4]等miRNA。通过生物信息学分析预测，这些miRNA的靶基因富集在细胞增殖、免疫应答等生物学过程相关的信号通路中，提示它们可能在牛流行热病毒感染后，通过调控这些生物学过程来影响病毒的感染和宿主细胞的反应。而在一个下调表达的簇中，miR-[X2]、miR-[X11]等miRNA的靶基因则主要富集在细胞凋亡、代谢调节等相关的信号通路中，表明这些miRNA可能通过抑制相关信号通路，来改变细胞的生理状态，从而适应牛流行热病毒的感染。聚类分析结果为进一步研究差异表达miRNA的功能和作用机制提供了重要线索，有助于揭示牛流行热病毒感染宿主细胞的分子调控网络。四、差异表达miRNA的鉴定与生物信息学分析4.1差异表达miRNA的验证4.1.1qRT-PCR验证结果为了验证高通量测序筛选出的差异表达miRNA的可靠性，采用qRT-PCR方法对部分差异表达miRNA进行了表达验证。选取了上调表达显著的miR-[X1]、miR-[X3]和下调表达显著的miR-[X2]、miR-[X11]等共8个miRNA进行qRT-PCR检测。以U6snRNA作为内参基因，每个样品设置3个生物学重复和3个技术重复，采用2⁻ΔΔCt法计算miRNA的相对表达量。qRT-PCR结果显示，所选的8个miRNA的表达趋势与高通量测序结果基本一致（图4）。其中，miR-[X1]在牛流行热病毒感染组中的相对表达量为[X]，相较于未感染组显著上调，与高通量测序中log₂(fold-change)值所反映的上调趋势相符；miR-[X2]在感染组中的相对表达量为[X]，明显低于未感染组，呈现下调趋势，也与高通量测序结果一致。对qRT-PCR和高通量测序结果进行相关性分析，结果显示两者具有显著的正相关关系，相关系数r=[X]（P＜0.01）。这表明qRT-PCR验证结果与高通量测序结果具有高度的一致性，进一步证明了高通量测序筛选出的差异表达miRNA的可靠性，为后续对这些miRNA的功能研究奠定了坚实的基础。[此处插入图4，图题：差异表达miRNA的qRT-PCR验证结果，图注：*表示P＜0.05，**表示P＜0.01，与未感染组相比；横坐标为miRNA名称，纵坐标为相对表达量，感染组和未感染组的数据以柱状图表示]4.1.2Northernblot验证结果除了qRT-PCR验证外，本研究还采用Northernblot方法对miR-[X1]和miR-[X2]这两个具有代表性的差异表达miRNA进行了进一步验证。将提取的总RNA进行变性聚丙烯酰胺凝胶电泳（PAGE），成功使RNA按照大小分离。随后，将凝胶中的RNA转移至尼龙膜上，并通过紫外交联固定RNA。用放射性同位素标记的miR-[X1]和miR-[X2]探针与尼龙膜上的RNA进行杂交，洗去未结合的探针后，通过放射自显影检测杂交信号。Northernblot结果显示，miR-[X1]在牛流行热病毒感染组中的杂交信号强度明显强于未感染组，表明其表达量上调；而miR-[X2]在感染组中的杂交信号强度较弱，表达量明显低于未感染组，呈现下调趋势（图5）。这一结果与qRT-PCR和高通量测序结果高度一致，进一步证实了miR-[X1]和miR-[X2]在牛流行热病毒感染宿主细胞后的差异表达情况。通过多种验证方法的相互印证，不仅增强了研究结果的可靠性和说服力，也为深入探究这些差异表达miRNA在牛流行热病毒感染过程中的功能和作用机制提供了有力的证据支持。[此处插入图5，图题：miR-[X1]和miR-[X2]的Northernblot验证结果，图注：M为Marker，1为未感染组，2为感染组；图片展示了杂交后的条带，上方条带为miR-[X1]，下方条带为miR-[X2]]4.2差异表达miRNA的生物信息学分析4.2.1靶基因预测为深入探究差异表达miRNA在牛流行热病毒感染过程中的作用机制，利用miRanda、TargetScan和PicTar这三种生物信息学工具对筛选出的差异表达miRNA进行靶基因预测。这三种工具的预测原理均基于miRNA与靶基因mRNA之间的互补配对原则，但在具体算法和参数设置上存在差异。miRanda是一款经典的靶基因预测软件，它通过计算miRNA与靶基因mRNA之间的最小自由能，评估两者的结合稳定性。该软件在预测过程中，考虑了miRNA种子区域（miRNA的第2-8位核苷酸）与靶基因mRNA的互补配对情况，以及miRNA与mRNA结合形成的双链结构的热力学稳定性。TargetScan则主要依据miRNA种子区域与靶基因mRNA3'UTR的互补配对来预测靶基因。它通过分析不同物种间靶位点的保守性，提高预测的准确性。在预测过程中，TargetScan还会考虑一些其他因素，如靶位点周围的序列特征、mRNA的二级结构等，以更全面地评估miRNA与靶基因的相互作用。PicTar的预测算法综合考虑了miRNA与靶基因mRNA的互补配对、结合位点的保守性以及结合位点在mRNA二级结构中的可及性等因素。它通过构建一个综合的评分系统，对每个潜在的靶基因进行评分，从而筛选出最有可能的靶基因。为确保预测结果的可靠性，取这三个工具预测结果的交集作为最终的靶基因集合。经过分析，共预测得到[X]个靶基因，这些靶基因涉及细胞代谢、信号传导、免疫应答等多个生物学过程。部分差异表达miRNA及其预测的靶基因见表2。表2：部分差异表达miRNA及其预测的靶基因miRNA名称预测靶基因数量部分预测靶基因名称miR-[X1][X]Gene1、Gene2、Gene3miR-[X2][X]Gene4、Gene5、Gene6miR-[X3][X]Gene7、Gene8、Gene94.2.2靶基因功能富集分析对预测得到的[X]个靶基因进行GO功能富集分析和KEGG通路富集分析，以揭示差异表达miRNA在牛流行热病毒感染过程中可能参与的生物学过程和信号通路。GO功能富集分析从生物过程、细胞组分和分子功能三个层面进行注释和富集分析。在生物过程方面，靶基因显著富集于细胞增殖、凋亡调控、免疫应答等过程（图6）。其中，在细胞增殖相关的GOterms中，如“细胞周期调控”（GO:0051726）、“DNA复制”（GO:0006260）等，有大量靶基因富集，这表明差异表达miRNA可能通过调控这些基因，影响宿主细胞的增殖能力，进而对牛流行热病毒的感染和复制产生影响。在凋亡调控相关的GOterms中，“细胞凋亡过程的正调控”（GO:0043065）、“细胞凋亡的内在调控途径”（GO:0097199）等也有显著富集的靶基因，暗示miRNA可能参与调节宿主细胞在病毒感染后的凋亡过程，这对于维持细胞内环境稳定以及病毒的生存和传播具有重要意义。在免疫应答方面，“抗病毒免疫应答”（GO:0051607）、“T细胞活化”（GO:0042110）等GOterms富集了众多靶基因，说明差异表达miRNA可能在宿主抗病毒免疫反应中发挥关键作用。[此处插入图6，图题：靶基因GO功能富集分析（生物过程），图注：横坐标为富集的GOterms，纵坐标为富集的靶基因数量，不同颜色表示不同的富集程度，颜色越深表示富集越显著]在细胞组分层面，靶基因主要富集在细胞核、细胞质和细胞膜等部位（图7）。其中，在细胞核相关的GOterms中，“细胞核”（GO:0005634）、“染色体”（GO:0005694）等富集了大量靶基因，这表明差异表达miRNA的靶基因在细胞核内参与了基因转录、染色质结构调控等重要过程。在细胞质和细胞膜相关的GOterms中，也有一定数量的靶基因富集，如“细胞质”（GO:0005737）、“细胞膜”（GO:0005886）等，说明这些靶基因在细胞的物质运输、信号传递等方面发挥作用。[此处插入图7，图题：靶基因GO功能富集分析（细胞组分），图注：横坐标为富集的GOterms，纵坐标为富集的靶基因数量，不同颜色表示不同的富集程度，颜色越深表示富集越显著]从分子功能角度来看，靶基因主要富集在核酸结合、蛋白质结合和酶活性等功能（图8）。在核酸结合相关的GOterms中，“DNA结合”（GO:0003677）、“RNA结合”（GO:0003723）等有较多靶基因富集，这表明这些靶基因可能通过与核酸相互作用，参与基因表达的调控。在蛋白质结合相关的GOterms中，“蛋白质异二聚化活性”（GO:0046982）、“蛋白质同二聚化活性”（GO:0042803）等也有显著富集，说明靶基因可能通过与其他蛋白质形成复合物，参与细胞内的各种生物学过程。此外，在酶活性相关的GOterms中，“蛋白激酶活性”（GO:0004672）、“磷酸酶活性”（GO:0016791）等富集了一定数量的靶基因，暗示这些靶基因可能通过调节酶的活性，影响细胞的代谢和信号传导。[此处插入图8，图题：靶基因GO功能富集分析（分子功能），图注：横坐标为富集的GOterms，纵坐标为富集的靶基因数量，不同颜色表示不同的富集程度，颜色越深表示富集越显著]KEGG通路富集分析结果显示，靶基因显著富集在多条信号通路中（图9）。其中，“Toll样受体信号通路”（ko04620）是重要的免疫相关信号通路，该通路中富集了多个靶基因。Toll样受体在识别病原体相关分子模式（PAMP）后，会激活下游的信号传导，诱导细胞因子和趋化因子的表达，从而启动免疫应答。差异表达miRNA的靶基因富集在该通路中，表明miRNA可能通过调控Toll样受体信号通路，影响宿主细胞对牛流行热病毒的免疫识别和应答。“MAPK信号通路”（ko04010）也是一条关键的信号通路，它参与细胞增殖、分化、凋亡等多种生物学过程。在牛流行热病毒感染过程中，MAPK信号通路的激活与病毒的复制和宿主细胞的应激反应密切相关。本研究中，该通路富集了大量靶基因，提示差异表达miRNA可能通过调节MAPK信号通路，影响病毒感染后的细胞生理状态。此外，“细胞凋亡信号通路”（ko04210）、“PI3K-Akt信号通路”（ko04151）等也有显著富集的靶基因，这些通路在细胞生存、增殖和凋亡调控中发挥重要作用，进一步说明差异表达miRNA可能通过调控这些信号通路，参与牛流行热病毒感染过程中宿主细胞的生理病理变化。[此处插入图9，图题：靶基因KEGG通路富集分析，图注：横坐标为富集的KEGG通路，纵坐标为富集的靶基因数量，不同颜色表示不同的富集程度，颜色越深表示富集越显著]4.2.3miRNA-mRNA调控网络构建基于差异表达miRNA及其预测的靶基因，利用Cytoscape软件构建miRNA-mRNA调控网络。在构建的调控网络中，每个节点代表一个miRNA或mRNA，边表示miRNA与mRNA之间的调控关系，即miRNA对mRNA的靶向作用。通过对调控网络的分析，可以直观地了解差异表达miRNA在网络中的调控作用和地位。调控网络分析结果显示，一些miRNA在网络中处于核心调控地位，与多个靶基因存在相互作用。例如，miR-[X1]与[X]个靶基因存在调控关系，这些靶基因涉及细胞增殖、免疫应答等多个重要生物学过程。miR-[X1]可能通过同时调控多个靶基因，在牛流行热病毒感染过程中发挥关键的调控作用。进一步分析发现，部分靶基因也受到多个miRNA的共同调控。例如，Gene1同时受到miR-[X1]、miR-[X3]和miR-[X5]等多个miRNA的靶向作用。这种复杂的调控关系表明，在牛流行热病毒感染宿主细胞后，miRNA-mRNA调控网络通过多种方式协同作用，精细地调节宿主细胞的生理功能，以应对病毒的感染。为了深入分析关键miRNA在调控网络中的作用，对网络中的节点进行拓扑学分析，计算节点的度（degree）、介数中心性（betweennesscentrality）和接近中心性（closenesscentrality）等参数。度表示与节点直接相连的边的数量，反映了节点在网络中的连接紧密程度；介数中心性衡量节点在网络中最短路径上的出现频率，体现了节点在信息传递中的重要性；接近中心性则表示节点到其他所有节点的最短路径的平均值，反映了节点在网络中的信息传播效率。通过分析这些参数，确定了miR-[X1]、miR-[X2]等几个关键miRNA。这些关键miRNA具有较高的度、介数中心性和接近中心性，表明它们在调控网络中与众多靶基因相互作用，在信息传递和调控过程中发挥着核心作用。对关键miRNA的进一步研究，将有助于深入揭示牛流行热病毒感染宿主细胞的分子机制。五、差异表达miRNA的功能验证与机制探讨5.1miRNA功能验证实验结果5.1.1miRNA过表达对病毒感染的影响为深入探究差异表达miRNA在牛流行热病毒感染过程中的功能，构建了miR-[X1]过表达载体，并将其转染至宿主细胞中。通过qRT-PCR检测发现，转染miR-[X1]过表达载体后，细胞内miR-[X1]的表达量显著上调，与阴性对照组相比，上调倍数达到[X]倍（图10A），表明miR-[X1]过表达载体成功转染并发挥作用。在转染后的细胞中接种牛流行热病毒，检测病毒感染相关指标的变化。结果显示，与阴性对照组相比，过表达miR-[X1]组的病毒滴度显著降低（图10B）。在感染后24h，阴性对照组的病毒滴度为[X]PFU/mL，而过表达miR-[X1]组的病毒滴度仅为[X]PFU/mL，降低了约[X]倍。同时，采用qRT-PCR检测病毒的核酸拷贝数，发现过表达miR-[X1]组的病毒核酸拷贝数也明显低于阴性对照组（图10C）。在感染后12h，阴性对照组的病毒核酸拷贝数为[X]copies/μL，而过表达miR-[X1]组为[X]copies/μL，减少了[X]%。此外，通过Westernblot检测病毒蛋白的表达水平，结果表明过表达miR-[X1]能够显著抑制病毒蛋白的表达（图10D）。以病毒糖蛋白G为例，过表达miR-[X1]组中病毒糖蛋白G的表达量相较于阴性对照组降低了[X]%。上述结果表明，过表达miR-[X1]能够有效抑制牛流行热病毒的感染和复制，在病毒感染过程中发挥负向调控作用。[此处插入图10，图题：miR-[X1]过表达对牛流行热病毒感染的影响，图注：A为miR-[X1]过表达载体转染后miR-[X1]的表达量检测；B为病毒滴度检测；C为病毒核酸拷贝数检测；D为病毒蛋白表达水平的Westernblot检测；*表示P＜0.05，**表示P＜0.01，与阴性对照组相比]5.1.2miRNA敲低对病毒感染的影响为进一步验证miRNA的功能，设计并合成了miR-[X2]的抑制剂（antagomir），将其转染至宿主细胞中，以敲低miR-[X2]的表达。qRT-PCR检测结果显示，转染miR-[X2]抑制剂后，细胞内miR-[X2]的表达量显著下降，与阴性对照组相比，降低了[X]%（图11A），表明miR-[X2]抑制剂成功发挥作用。在敲低miR-[X2]表达的细胞中接种牛流行热病毒，检测病毒感染相关指标。结果发现，与阴性对照组相比，敲低miR-[X2]组的病毒滴度显著升高（图11B）。在感染后36h，阴性对照组的病毒滴度为[X]PFU/mL，而敲低miR-[X2]组的病毒滴度达到[X]PFU/mL，升高了约[X]倍。通过qRT-PCR检测病毒核酸拷贝数，结果显示敲低miR-[X2]组的病毒核酸拷贝数明显高于阴性对照组（图11C）。在感染后24h，阴性对照组的病毒核酸拷贝数为[X]copies/μL，敲低miR-[X2]组为[X]copies/μL，增加了[X]%。此外，Westernblot检测结果表明，敲低miR-[X2]能够显著促进病毒蛋白的表达（图11D）。以病毒核蛋白N为例，敲低miR-[X2]组中病毒核蛋白N的表达量相较于阴性对照组增加了[X]%。综合以上结果，敲低miR-[X2]能够促进牛流行热病毒的感染和复制，说明miR-[X2]在病毒感染过程中发挥正向调控作用，即对病毒感染具有抑制作用，当miR-[X2]表达被敲低时，病毒感染能力增强。[此处插入图11，图题：miR-[X2]敲低对牛流行热病毒感染的影响，图注：A为miR-[X2]抑制剂转染后miR-[X2]的表达量检测；B为病毒滴度检测；C为病毒核酸拷贝数检测；D为病毒蛋白表达水平的Westernblot检测；*表示P＜0.05，**表示P＜0.01，与阴性对照组相比]5.2miRNA调控病毒感染的机制探讨5.2.1miRNA与靶基因的相互作用验证为进一步探究miR-[X1]和miR-[X2]在牛流行热病毒感染过程中的作用机制，采用荧光素酶报告实验验证它们与靶基因的直接相互作用。根据生物信息学预测结果，选择miR-[X1]的潜在靶基因Gene1和miR-[X2]的潜在靶基因Gene4进行验证。首先构建荧光素酶报告载体，将Gene1和Gene4的3'UTR区域分别插入到荧光素酶报告基因载体（如pGL3-Basic）的下游，得到pGL3-Gene1-3'UTR和pGL3-Gene4-3'UTR重组质粒。同时，针对miR-[X1]和miR-[X2]，分别合成其模拟物（mimics）和阴性对照（negativecontrol）。将宿主细胞以[细胞接种密度]接种于24孔细胞培养板中，每孔加入500μL细胞培养液。培养24h后，待细胞贴壁且融合度达到50%-60%时，进行转染实验。将pGL3-Gene1-3'UTR或pGL3-Gene4-3'UTR重组质粒分别与miR-[X1]mimics、miR-[X2]mimics或相应的阴性对照共转染至细胞中。转染采用脂质体转染试剂（如Lipofectamine3000），按照试剂说明书进行操作。转染48h后，收集细胞，采用双荧光素酶报告基因检测试剂盒（Dual-LuciferaseReporterAssayKit）检测荧光素酶活性。实验结果显示，与阴性对照相比，共转染miR-[X1]mimics和pGL3-Gene1-3'UTR重组质粒的细胞中，荧光素酶活性显著降低（图12A），表明miR-[X1]能够与Gene1的3'UTR区域结合，抑制荧光素酶报告基因的表达，从而验证了miR-[X1]与Gene1之间存在直接的相互作用。对于miR-[X2]，共转染miR-[X2]mimics和pGL3-Gene4-3'UTR重组质粒的细胞中，荧光素酶活性也明显下降（图12B），说明miR-[X2]可以与Gene4的3'UTR区域相互作用，抑制其表达。为了进一步确定miRNA与靶基因的结合位点，对Gene1和Gene4的3'UTR区域中预测的miRNA结合位点进行突变。构建含有突变结合位点的荧光素酶报告载体pGL3-Gene1-3'UTR-mut和pGL3-Gene4-3'UTR-mut。重复上述转染和检测实验，结果发现，当结合位点突变后，miR-[X1]mimics与pGL3-Gene1-3'UTR-mut共转染以及miR-[X2]mimics与pGL3-Gene4-3'UTR-mut共转染的细胞中，荧光素酶活性与阴性对照相比无显著差异（图12A、B）。这表明miR-[X1]和miR-[X2]确实是通过与靶基因3'UTR区域的特定结合位点相互作用，来调控靶基因的表达。[此处插入图12，图题：miR-[X1]和miR-[X2]与靶基因相互作用的荧光素酶报告实验结果，图注：A为miR-[X1]与Gene1的相互作用验证；B为miR-[X2]与Gene4的相互作用验证；*表示P＜0.05，**表示P＜0.01，与阴性对照组相比；NC为阴性对照，WT为野生型3'UTR，Mut为突变型3'UTR]5.2.2miRNA对宿主细胞信号通路的影响基于生物信息学分析中靶基因在信号通路的富集结果，深入研究miR-[X1]和miR-[X2]通过调控靶基因对宿主细胞内信号通路的激活或抑制作用，以及对病毒感染的间接影响。研究发现，miR-[X1]通过抑制靶基因Gene1的表达，影响Toll样受体信号通路。Toll样受体信号通路在宿主的免疫应答中发挥关键作用，当病原体入侵时，Toll样受体识别病原体相关分子模式，激活下游的信号传导，诱导细胞因子和趋化因子的表达，从而启动免疫应答。在牛流行热病毒感染宿主细胞后，miR-[X1]表达上调，抑制Gene1的表达。Gene1是Toll样受体信号通路中的关键分子，它的表达受到抑制后，导致Toll样受体信号通路的激活受到阻碍。通过Westernblot检测通路中关键蛋白的磷酸化水平，发现与未感染组相比，感染组中Toll样受体信号通路下游关键蛋白（如MyD88、NF-κB等）的磷酸化水平显著降低（图13A）。这表明miR-[X1]通过抑制靶基因Gene1，负向调控Toll样受体信号通路，进而削弱了宿主细胞的免疫应答能力，为病毒的感染和复制创造了有利条件。对于miR-[X2]，它通过调控靶基因Gene4，影响MAPK信号通路。MAPK信号通路参与细胞增殖、分化、凋亡等多种生物学过程。在牛流行热病毒感染过程中，miR-[X2]表达下调，导致Gene4的表达上调。Gene4是MAPK信号通路的正调控因子，它的表达增加会激活MAPK信号通路。通过检测通路中关键蛋白的磷酸化水平，发现感染组中MAPK信号通路下游关键蛋白（如ERK1/2、JNK等）的磷酸化水平明显高于未感染组（图13B）。激活的MAPK信号通路促进了细胞的增殖和存活，这可能有助于宿主细胞在病毒感染时维持一定的生理功能，同时也可能为病毒的复制提供了更多的物质和能量。然而，过度激活的MAPK信号通路也可能导致细胞的应激反应和炎症反应加剧，对宿主细胞造成损伤。综上所述，miR-[X1]和miR-[X2]通过调控不同的靶基因，分别对Toll样受体信号通路和MAPK信号通路产生影响，进而间接影响牛流行热病毒的感染过程。这些结果揭示了miRNA在牛流行热病毒感染宿主细胞过程中，通过调节宿主细胞信号通路来参与病毒与宿主相互作用的分子机制。[此处插入图13，图题：miR-[X1]和miR-[X2]对宿主细胞信号通路关键蛋白磷酸化水平的影响，图注：A为Toll样受体信号通路；B为MAPK信号通路；*表示P＜0.05，**表示P＜0.01，与未感染组相比；Con为未感染组，Inf为感染组]六、结论与展望6.1研究主要结论本研究成功筛选并鉴定出牛流行热病毒感染宿主细胞后差异表达的miRNA，通过一系列实验深入分析了它们的功能和调控机制，取得了以下主要研究成果：差异表达miRNA的筛选与鉴定：利用高通量测序技术，对牛流行热病毒感染组和未感染组宿主细胞的miRNA表达谱进行分析，筛选出[X]个差异表达的miRNA，其中上调表达的miRNA有[X]个，下调表达的miRNA有[X]个。通过qRT-PCR和Northernblot验证，证实了这些差异表达miRNA的可靠性，为后续研究奠定了坚实基础。生物信息学分析：对差异表达miRNA进行靶基因预测，共得到[X]个靶基因。GO功能富集分析和KEGG通路富集分析表明，这些靶基因涉及细胞增殖、凋亡、免疫应答等多个生物学过程，以及Toll样受体信号通路、MAPK信号通路等多条关键信号通路。进一步构建miRNA-mRNA调控网络，确定了miR-[X1]、miR-[X2]等关键miRNA，它们在调控网络中与众多靶基因相互作用，可能在牛流行热病毒感染过程中发挥核心调控作用。功能验证与机制探讨：通过构建miRNA过表达载体和miRNA抑制剂，对miR-[X1]和miR-[X2]的功能进行验证。结果表明，过表达miR-[X1]能够显著抑制牛流行热病毒的感染和复制，而敲低miR-[X2]则促进病毒感染。荧光素酶报告实验证实了miR-[X1]与Gene1、miR-[X2]与Gene4之间存在直接的相互作用。进一步研究发现，miR-[X1]通过抑制靶基因Gene1，负向调控Toll样受体信号通路，削弱宿主细胞的免疫应答能力；miR-[X2]通过调控靶基因Gene4，激活MAPK信号通路，影响细胞的增殖和存活，进而间接影响牛流行热病毒的感染过程。6.2研究的创新点与不足之处本研究在牛流行热病毒感染宿主细胞差异表达miRNA的筛选与鉴定方面具有一定的创新点。首次系统地对牛流行热病毒感染宿主细胞后的miRNA表达谱进行全面分析，相较于以往对牛流行热病毒感染机制的研究，本研究从miRNA层面入手，为揭示病毒与宿主细胞的相互作用提供了全新的视角。通过高通量测序技术，不仅鉴定出已知miRNA的差异表达情况，还预测并发现了潜在的新miRNA