基于组学技术解析不同毒力PRRSVs感染猪肺组织的分子机制

基于组学技术解析不同毒力PRRSVs感染猪肺组织的分子机制一、引言1.1研究背景与意义猪繁殖与呼吸综合征（PorcineReproductiveandRespiratorySyndrome，PRRS），俗称蓝耳病，是由猪繁殖与呼吸综合征病毒（PorcineReproductiveandRespiratorySyndromeVirus，PRRSV）引起的一种对全球养猪业危害极为严重的传染病。自1987年首次在美国被发现以来，PRRSV迅速传播至世界各地，给养猪业带来了巨大的经济损失。PRRSV主要侵害猪的免疫系统和生殖系统，感染母猪常出现繁殖障碍，如流产、早产、产死胎、木乃伊胎和弱仔等；感染仔猪则表现出严重的呼吸道症状，如呼吸困难、咳嗽、体温升高，死亡率显著增加。同时，PRRSV还可导致猪群生长发育迟缓、饲料转化率降低，使养殖成本大幅上升。此外，PRRSV感染还会引起猪的免疫抑制，使猪对其他病原体的易感性增强，从而引发多种继发感染和混合感染，进一步加重病情，增加治疗难度和死亡率。在全球范围内，PRRSV的持续传播和不断变异给养猪业造成了沉重打击。美国作为养猪业发达的国家之一，每年因PRRSV感染导致的经济损失高达数亿美元。在中国，养猪业是农业的重要支柱产业，PRRSV的流行也给养猪业带来了巨大的冲击。据相关统计数据显示，中国每年因PRRS造成的直接经济损失可达数十亿元人民币。而且，随着养猪业规模化和集约化程度的不断提高，PRRSV的传播速度更快，危害范围更广，防控难度也越来越大。PRRSV具有高度的变异性，根据其基因组序列和抗原性的差异，可分为欧洲型（基因1型）和北美型（基因2型）两个基因型，两基因型之间的核苷酸同源性仅为60%-70%。不同基因型的PRRSV在毒力、致病性和免疫原性等方面存在显著差异，即使同一基因型内的不同毒株之间也存在较大的变异。这些变异使得PRRSV的致病机制变得极为复杂，给疫苗的研发和防控措施的制定带来了极大的挑战。目前，市场上已有的疫苗虽然在一定程度上对PRRSV的防控起到了积极作用，但由于病毒的不断变异，疫苗的保护效果往往不尽如人意，无法完全满足养猪业的实际需求。深入研究不同毒力PRRSV的致病机制，对于有效防控PRRS具有重要的理论和实际意义。从理论角度来看，通过探究不同毒力PRRSV感染宿主后引起的基因表达变化、信号通路激活以及蛋白质表达和修饰的改变等，能够揭示PRRSV感染的分子机制，丰富对病毒与宿主相互作用关系的认识，为病毒学和免疫学的发展提供新的理论依据。在实际应用方面，明确不同毒力PRRSV的致病机制，有助于开发更加有效的诊断方法，实现对PRRSV感染的早期准确检测；能够为研发新型、高效的疫苗和治疗药物提供关键靶点，提高疫苗的保护效力和药物的治疗效果；还可以为养猪业制定科学合理的防控策略提供有力支持，减少PRRSV的传播和感染，降低经济损失，保障养猪业的健康可持续发展。1.2PRRSV概述猪繁殖与呼吸综合征病毒（PRRSV）属于套式病毒目（Nidovirales）、动脉炎病毒科（Arteriviridae）、动脉炎病毒属（Arterivirus），是一种有囊膜的单股正链RNA病毒。其基因组长度约为15kb，包含9个开放阅读框（OpenReadingFrames，ORFs），分别编码病毒的结构蛋白和非结构蛋白。这些蛋白在病毒的复制、装配、感染和致病过程中发挥着关键作用。PRRSV具有高度的变异性，其变异主要源于基因突变和基因重组。不同地区、不同猪场流行的PRRSV毒株在基因组序列上存在差异，这种变异使得病毒的抗原性、毒力和致病性也发生改变，从而导致疫苗的免疫效果不稳定，给PRRS的防控带来了极大的困难。例如，高致病性PRRSV（HP-PRRSV）的出现，相较于经典毒株，其毒力更强，致病性更高，能引起猪更为严重的临床症状和更高的死亡率。PRRSV主要通过接触传播和空气传播两种途径在猪群中扩散。接触传播包括直接接触和间接接触，直接接触是指感染猪与健康猪之间的直接身体接触，如鼻对鼻接触、共同采食和饮水等；间接接触则是通过被病毒污染的饲料、饮水、器具、车辆以及人员等传播媒介，将病毒传播给易感猪。空气传播是PRRSV传播的重要方式之一，病毒可以随着空气中的飞沫和气溶胶远距离传播，尤其是在通风不良、猪群密度较大的养殖场，空气传播更容易导致病毒的快速扩散，使疫情迅速蔓延。此外，PRRSV还可通过精液传播，感染的公猪可通过精液将病毒传播给母猪，进而导致母猪繁殖障碍和仔猪感染。感染PRRSV的猪会表现出多种临床症状，母猪主要表现为繁殖障碍，如在妊娠后期发生流产，产出死胎、弱胎、木乃伊胎等；部分母猪还会出现返情、产后不发情等情况，严重影响母猪的繁殖性能。仔猪感染后，常出现呼吸困难、咳嗽、气喘等呼吸道症状，同时伴有体温升高、精神沉郁、食欲减退、消瘦等全身性症状，死亡率较高，尤其是在保育阶段，仔猪的死亡率可高达20%-50%。生长育肥猪感染PRRSV后，生长速度明显减缓，饲料转化率降低，料肉比增加，导致养殖成本上升，经济效益下降。此外，PRRSV感染还会引起猪的免疫抑制，使猪更容易感染其他病原体，如猪圆环病毒2型（Porcinecircovirustype2，PCV2）、猪流感病毒（Swineinfluenzavirus，SIV）、副猪嗜血杆菌（Haemophilusparasuis，HPS）等，引发混合感染和继发感染，加重病情，增加治疗难度和死亡率。PRRSV的感染给全球养猪业带来了沉重的经济负担。在美国，每年因PRRSV感染导致的经济损失高达数亿美元，包括直接损失和间接损失。直接损失主要来自病死猪的淘汰、治疗费用以及繁殖性能下降导致的仔猪数量减少等；间接损失则包括因生长育肥猪生长缓慢、饲料转化率降低而增加的养殖成本，以及为防控疫情而投入的大量人力、物力和财力等。在中国，养猪业是农业的重要支柱产业之一，PRRSV的流行对养猪业造成了巨大的冲击。据相关统计数据显示，中国每年因PRRS造成的直接经济损失可达数十亿元人民币，间接损失更是难以估量。而且，随着养猪业规模化和集约化程度的不断提高，PRRSV的传播风险进一步增加，经济损失也呈上升趋势。因此，有效防控PRRSV感染对于保障全球养猪业的健康发展和经济效益具有至关重要的意义。1.3转录组学与蛋白质组学技术转录组学是一门在整体水平上研究细胞中所有基因转录及转录调控规律的学科。转录组是指特定发育阶段或生理条件下，细胞内所有转录产物的集合，包括信使RNA（mRNA）、核糖体RNA（rRNA）、转运RNA（tRNA）及非编码RNA（ncRNA）等。通过转录组学技术，如RNA测序（RNA-Seq）、基因芯片等，可以全面、准确地分析转录组的组成和表达水平，揭示基因的转录模式、转录本结构以及不同样本间基因表达的差异。在病毒致病机制研究中，转录组学能够深入探究病毒感染宿主细胞后，宿主基因表达谱的变化情况。例如，在研究流感病毒感染宿主细胞时，通过转录组学分析发现，宿主细胞中与免疫应答、炎症反应相关的基因表达显著上调，而与细胞代谢、蛋白质合成相关的基因表达则受到抑制。这表明流感病毒感染会引起宿主细胞免疫反应的激活和细胞代谢的改变，从而影响细胞的正常功能，为深入理解流感病毒的致病机制提供了重要线索。蛋白质组学是以蛋白质组为研究对象，旨在从整体水平上研究蛋白质的组成、结构、功能及其相互作用的学科。蛋白质组是指一个基因组、一种生物或一个细胞/组织所表达的全部蛋白质。蛋白质作为生命活动的直接执行者，其表达水平、修饰状态和相互作用等变化直接反映了细胞的生理状态和功能变化。常用的蛋白质组学技术包括双向凝胶电泳（2-DE）、液相色谱-串联质谱（LC-MS/MS）、蛋白质芯片等。这些技术可以实现对蛋白质的分离、鉴定和定量分析，以及对蛋白质翻译后修饰、蛋白质-蛋白质相互作用等的研究。在病毒感染机制研究中，蛋白质组学发挥着重要作用。以乙肝病毒感染为例，利用蛋白质组学技术对感染乙肝病毒的肝细胞进行分析，发现了一系列与乙肝病毒感染相关的差异表达蛋白质，这些蛋白质涉及细胞信号转导、免疫调节、代谢等多个生物学过程。进一步研究这些蛋白质的功能和相互作用，有助于深入了解乙肝病毒的感染机制，为开发新的治疗方法和药物靶点提供了重要依据。转录组学和蛋白质组学技术在病毒致病机制研究中具有广泛的应用前景。在病毒感染宿主细胞的过程中，病毒基因的表达和宿主基因的表达都会发生复杂的变化，这些变化相互影响，共同决定了病毒的感染进程和宿主的免疫反应。通过转录组学和蛋白质组学技术，可以从基因转录和蛋白质表达两个层面全面解析病毒与宿主之间的相互作用机制，揭示病毒致病的分子基础。在研究人类免疫缺陷病毒（HIV）感染时，综合运用转录组学和蛋白质组学技术，不仅可以分析HIV感染后宿主细胞基因表达谱的变化，还能研究蛋白质表达水平和修饰状态的改变，从而深入了解HIV感染导致宿主免疫功能受损的机制。这对于开发有效的抗HIV药物和治疗策略具有重要意义。此外，转录组学和蛋白质组学技术还可以用于筛选和鉴定病毒感染相关的生物标志物，为病毒感染的早期诊断、病情监测和预后评估提供新的方法和指标。1.4研究目的与内容本研究旨在运用转录组学和蛋白质组学技术，深入探究不同毒力PRRSV感染后猪肺组织在基因转录和蛋白质表达层面的变化，从而初步揭示PRRSV的致病机制，为PRRS的防控提供理论依据和潜在的分子靶点。具体研究内容如下：不同毒力PRRSV感染猪模型的建立：选择合适的实验猪，分别用高毒力和低毒力PRRSV毒株进行感染，设置对照组。观察感染猪的临床症状，定期采集血液样本检测病毒血症水平，以确定感染模型的成功建立，并为后续的转录组和蛋白质组学分析提供实验材料。肺组织转录组学研究：在感染后的特定时间点采集感染组和对照组猪的肺组织样本，提取总RNA，利用RNA-Seq技术进行转录组测序。通过生物信息学分析，筛选出不同毒力PRRSV感染组与对照组之间的差异表达基因（DEGs），对这些DEGs进行功能注释和富集分析，明确其参与的生物学过程、细胞组成和分子功能。构建基因共表达网络，挖掘关键基因和核心调控模块，探究不同毒力PRRSV感染后宿主基因表达的调控机制。肺组织蛋白质组学研究：采用相同时间点采集的肺组织样本，提取总蛋白质，运用基于液相色谱-串联质谱（LC-MS/MS）的蛋白质组学技术，对蛋白质进行分离、鉴定和定量分析。筛选出不同毒力PRRSV感染组与对照组之间的差异表达蛋白质（DEPs），对DEPs进行功能注释和富集分析，确定其参与的生物学途径和分子功能。通过蛋白质-蛋白质相互作用网络分析，揭示蛋白质之间的相互关系和功能协同作用，寻找与PRRSV致病相关的关键蛋白质和信号通路。转录组与蛋白质组数据的整合分析：将转录组学和蛋白质组学数据进行整合，分析差异表达基因和差异表达蛋白质之间的关联，挖掘在转录和翻译水平均发生显著变化的基因和蛋白质，进一步验证其在PRRSV致病过程中的作用。综合运用生物信息学方法，构建病毒-宿主相互作用的分子调控网络，全面解析不同毒力PRRSV感染后肺组织的分子变化机制，为深入理解PRRSV的致病机制提供更全面、系统的信息。二、材料与方法2.1实验材料实验用猪：选取60头30日龄左右、体重相近、健康状况良好的SPF级仔猪，购自[具体猪场名称]。这些仔猪在实验前经过严格的检疫，确保未感染PRRSV及其他常见猪传染病，并在隔离的实验动物房中适应性饲养一周，期间给予常规饲料和清洁饮水。病毒毒株：高毒力PRRSV毒株[毒株具体名称1]，分离自[分离地点1]的发病猪群，经测序鉴定和毒力测定，具有典型的高致病性特征，能引起猪严重的临床症状和高死亡率；低毒力PRRSV毒株[毒株具体名称2]，分离自[分离地点2]，其毒力相对较弱，感染猪后临床症状较轻。两种毒株均由[病毒保存单位]提供，并在实验室中进行复苏和扩增，保存于-80℃冰箱备用。主要试剂：Trizol试剂（Invitrogen公司），用于提取组织总RNA；逆转录试剂盒（TaKaRa公司），将RNA逆转录为cDNA；SYBRGreen荧光定量PCR试剂盒（Roche公司），用于检测病毒血症水平和基因表达量；蛋白质提取试剂盒（ThermoFisherScientific公司），从肺组织中提取总蛋白质；胰蛋白酶（Promega公司），用于蛋白质酶解；乙腈、甲醇等色谱级有机溶剂（Merck公司），用于液相色谱-串联质谱分析；BCA蛋白定量试剂盒（Beyotime公司），测定蛋白质浓度；其他常规试剂如无水乙醇、异丙醇、***化钠、磷酸缓冲液等均为国产分析纯试剂。主要仪器设备：高速冷冻离心机（Eppendorf公司），用于样品离心分离；实时荧光定量PCR仪（ABI7500，ThermoFisherScientific公司），进行荧光定量PCR反应；核酸电泳仪（Bio-Rad公司），检测RNA和DNA的完整性和纯度；超低温冰箱（ThermoFisherScientific公司），保存病毒毒株和生物样品；超声波细胞破碎仪（SCIENTZ公司），破碎细胞以提取蛋白质；液相色谱-串联质谱仪（ThermoQExactiveHF，ThermoFisherScientific公司），用于蛋白质组学分析；恒温培养箱（ThermoFisherScientific公司），细胞培养和孵育；电子天平（Sartorius公司），称量试剂和样品。2.2实验方法2.2.1建立PRRSV感染模型将60头30日龄SPF级仔猪随机分为3组，每组20头。其中两组分别作为高毒力PRRSV感染组和低毒力PRRSV感染组，另一组作为对照组。感染前，对所有仔猪进行采血检测，确保其体内不存在PRRSV抗体及其他常见猪传染病病原体。高毒力PRRSV感染组仔猪经鼻腔接种高毒力PRRSV毒株[毒株具体名称1]，接种剂量为1×10^5TCID50（半数组织培养感染剂量）/头，接种时将病毒液用无菌PBS（磷酸盐缓冲液）稀释至合适浓度，缓慢滴入仔猪双侧鼻腔，每侧鼻腔接种量相同，接种过程中确保仔猪保持安静，避免病毒液流出。低毒力PRRSV感染组仔猪采用同样的接种方式，接种低毒力PRRSV毒株[毒株具体名称2]，接种剂量为1×10^5TCID50/头。对照组仔猪经鼻腔接种等量的无菌PBS。接种后，每天观察并记录仔猪的临床症状，包括体温、精神状态、食欲、呼吸状况、皮肤颜色等。定期采集仔猪的血液样本，采用实时荧光定量PCR方法检测血液中的病毒载量，以监测病毒血症水平。当感染组仔猪出现明显的临床症状，且病毒血症水平达到一定程度时，表明PRRSV感染模型成功建立。2.2.2样本采集与处理在PRRSV感染后的第3天、第7天和第14天，分别从感染组和对照组中随机选取5头仔猪进行安乐死。迅速采集猪的肺组织样本，每个样本采集约1g，分别取自左肺上叶、右肺上叶和右肺下叶等不同部位，以确保样本的代表性。采集后的肺组织样本立即放入预冷的无菌生理盐水中冲洗，去除表面的血液和杂质。然后将样本放入冻存管中，加入适量的RNA保护剂（如RNAlater），使样本完全浸没在保护剂中，以防止RNA降解。将冻存管迅速放入液氮中速冻5分钟，然后转移至-80℃冰箱中保存，直至进行后续的转录组学和蛋白质组学分析。对于用于蛋白质组学分析的肺组织样本，在采集和冲洗后，将其放入含有蛋白酶抑制剂的裂解缓冲液中，用组织匀浆器将组织充分匀浆，使细胞破碎，释放蛋白质。匀浆后的样本在4℃下以12000rpm离心15分钟，取上清液转移至新的离心管中，采用BCA蛋白定量试剂盒测定蛋白质浓度。将定量后的蛋白质样本分装成小份，加入适量的甘油（终浓度为10%-20%）以防止蛋白质在冻融过程中变性，然后放入液氮中速冻，再转移至-80℃冰箱中保存备用。2.2.3转录组学研究总RNA提取：从-80℃冰箱中取出保存的肺组织样本，在冰上解冻。使用Trizol试剂提取总RNA，具体操作按照Trizol试剂说明书进行。提取过程中，通过多次离心和洗涤步骤，去除蛋白质、DNA等杂质，以获得高纯度的RNA。提取后的RNA用超微量分光光度计测定其浓度和纯度，要求OD260/OD280比值在1.8-2.0之间，OD260/OD230比值大于2.0。同时，利用琼脂糖凝胶电泳检测RNA的完整性，观察28S和18SrRNA条带的亮度和清晰度，确保RNA无降解。文库构建与高通量测序：将提取的高质量总RNA送往专业的测序公司进行文库构建和高通量测序。首先，利用mRNA的poly-A尾巴特性，通过寡聚（dT）磁珠富集真核生物mRNA；对于原核生物RNA，则采用去除rRNA的方法富集mRNA。然后，将富集的mRNA进行片段化处理，以打断成适合测序的短片段。在逆转录酶的作用下，以mRNA为模板合成cDNA第一链，再合成cDNA第二链。对cDNA进行末端修复、加A尾和连接测序接头等一系列操作，构建成测序文库。文库构建完成后，通过Qubit荧光定量仪和Agilent2100生物分析仪对文库的浓度和质量进行检测。合格的文库在IlluminaHiSeq测序平台上进行双端测序，测序读长为150bp。数据分析：测序完成后，得到的原始数据首先进行质量控制，利用FastQC软件对原始测序数据进行评估，检查数据的质量分布、碱基组成、测序接头污染等情况。使用Trimmomatic软件去除低质量碱基、测序接头和含有过多N的序列，得到高质量的cleanreads。将cleanreads通过Hisat2软件比对到猪的参考基因组（如Susscrofa11.1）上，统计比对率。利用StringTie软件进行转录本组装和定量分析，计算每个基因的表达量，常用的表达量指标为FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）。差异基因筛选：使用DESeq2软件对感染组和对照组之间的基因表达数据进行差异分析，筛选出差异表达基因（DEGs）。设定筛选条件为|log2FC|>1且FDR（FalseDiscoveryRate）<0.05，其中log2FC表示两组基因表达量的对数2倍变化，FDR用于校正多重检验中的假阳性率。功能注释与富集分析：对筛选出的DEGs进行功能注释，通过BLAST软件将DEGs的核苷酸序列或氨基酸序列与多个数据库（如NCBI非冗余蛋白数据库（NR）、基因本体论数据库（GO）、京都基因与基因组百科全书数据库（KEGG）等）进行比对，获取基因的功能信息。利用clusterProfiler软件对DEGs进行GO功能富集分析和KEGG通路富集分析。GO富集分析从生物过程（BiologicalProcess）、细胞组成（CellularComponent）和分子功能（MolecularFunction）三个层面揭示DEGs参与的主要生物学过程和分子功能。KEGG通路富集分析则确定DEGs显著富集的代谢通路和信号转导通路，以了解病毒感染后宿主细胞内发生的关键生物学事件和信号调控网络。2.2.4蛋白质组学研究蛋白质提取与定量：从-80℃冰箱中取出保存的蛋白质样本，在冰上解冻。使用蛋白质提取试剂盒提取肺组织中的总蛋白质，提取过程中按照试剂盒说明书加入相应的试剂和进行操作，确保蛋白质的完整性和纯度。提取后的蛋白质采用BCA蛋白定量试剂盒测定浓度，具体操作按照试剂盒说明书进行。以牛血清白蛋白（BSA）为标准品，绘制标准曲线，根据样品的吸光度值从标准曲线中计算出蛋白质浓度。蛋白质酶解与肽段分离：将定量后的蛋白质样本进行酶解处理，采用胰蛋白酶在37℃下孵育过夜，将蛋白质酶解成肽段。酶解后的肽段通过固相萃取小柱（如C18小柱）进行脱盐和纯化处理，去除杂质和盐分。纯化后的肽段用真空浓缩仪浓缩至适当体积，用于后续的液相色谱-串联质谱分析。液相色谱-串联质谱分析（LC-MS/MS）：将浓缩后的肽段样品注入液相色谱系统（如ThermoScientificUltiMate3000RSLCnano系统）进行分离。采用C18反相色谱柱（如ThermoScientificAcclaimPepMapRSLCC18，2μm，100Å，75μm×50cm），以乙腈和0.1%甲酸水溶液为流动相进行梯度洗脱，使肽段在色谱柱上实现分离。分离后的肽段直接进入质谱仪（如ThermoQExactiveHF质谱仪）进行检测。质谱仪在数据依赖模式下运行，首先进行一级质谱扫描，采集肽段的母离子信息，然后选择强度较高的母离子进行二级质谱碎裂，采集子离子信息。通过对母离子和子离子的质量分析，获得肽段的氨基酸序列信息。蛋白质鉴定与定量：将质谱采集到的数据通过ProteomeDiscoverer软件进行分析。首先，将质谱数据与猪的蛋白质数据库（如UniProt猪数据库）进行比对，利用SequestHT算法进行蛋白质鉴定，确定肽段对应的蛋白质。采用Label-free定量方法对蛋白质进行定量分析，通过比较不同样本中肽段的峰面积或离子强度，计算蛋白质的相对表达量。差异蛋白筛选与分析：使用Perseus软件对感染组和对照组之间的蛋白质表达数据进行差异分析，筛选出差异表达蛋白质（DEPs）。设定筛选条件为|log2FC|>1且p-value<0.05，其中log2FC表示两组蛋白质表达量的对数2倍变化，p-value用于判断差异的显著性。对筛选出的DEPs进行功能注释和富集分析，注释数据库和分析方法与转录组学中的DEGs分析类似，通过与NR、GO、KEGG等数据库比对，揭示DEPs的功能和参与的生物学过程及信号通路。同时，利用STRING数据库构建蛋白质-蛋白质相互作用（PPI）网络，分析蛋白质之间的相互关系和功能协同作用，寻找网络中的关键节点蛋白质，这些关键蛋白质可能在PRRSV致病过程中发挥重要作用。2.2.5数据处理与分析转录组学和蛋白质组学数据的分析均在R语言环境下进行，结合多个生物信息学软件和工具包完成。对于转录组学数据，使用DESeq2进行差异基因分析，利用clusterProfiler进行GO和KEGG富集分析，通过ggplot2等绘图包绘制火山图、柱状图、热图等可视化图表，直观展示差异基因的分布、富集情况等信息。对于蛋白质组学数据，运用Perseus进行差异蛋白分析，同样使用clusterProfiler进行功能注释和富集分析，利用Cytoscape软件构建和分析PPI网络，并进行可视化展示。在整合分析转录组和蛋白质组数据时，将差异表达基因和差异表达蛋白质进行关联分析，寻找在转录和翻译水平均发生显著变化的基因和蛋白质，通过韦恩图等方式展示两者的重叠情况。利用生物信息学方法挖掘关键的基因-蛋白质调控关系，构建病毒-宿主相互作用的分子调控网络，深入解析不同毒力PRRSV感染后肺组织的分子变化机制。同时，对整合分析得到的关键基因和蛋白质进行验证，可采用实时荧光定量PCR、蛋白质免疫印迹（Westernblot）等实验技术，进一步确认其在PRRSV致病过程中的作用。三、不同毒力PRRSVs感染后肺组织的转录组学分析3.1转录组测序数据质量评估对不同毒力PRRSV感染组及对照组猪肺组织样本进行转录组测序后，首先对原始测序数据进行了全面的质量评估。利用FastQC软件对原始测序数据进行分析，结果显示各样本测序数据的碱基质量分布较为均匀，大部分碱基的质量值（Q值）均在30以上，表明测序数据的准确性较高，错误率较低。在碱基组成方面，A、T、C、G四种碱基的含量相对稳定，无明显的碱基偏好性，符合正常测序数据的特征。同时，对测序接头污染情况进行检查，发现各样本中测序接头的污染率均低于1%，说明在文库构建过程中接头去除较为彻底，不会对后续数据分析产生显著影响。使用Trimmomatic软件对原始数据进行过滤，去除低质量碱基、测序接头和含有过多N的序列后，得到高质量的cleanreads。统计各样本的cleanreads数量及与参考基因组的比对率，结果表明，各样本的cleanreads数量均达到了5000万以上，能够满足后续分析的需求。在比对率方面，高毒力PRRSV感染组样本的平均比对率为85.6%，低毒力PRRSV感染组样本的平均比对率为84.8%，对照组样本的平均比对率为86.2%，三组样本的比对率均处于较高水平，表明测序数据与猪参考基因组具有较好的匹配性，可用于后续的基因表达分析。通过对转录组测序数据的质量评估，证明了本次测序数据质量良好，能够为后续筛选差异表达基因、功能注释与富集分析等提供可靠的数据基础，确保了转录组学分析结果的准确性和可靠性。3.2差异表达基因筛选在完成转录组测序数据质量评估后，使用DESeq2软件对高毒力PRRSV感染组、低毒力PRRSV感染组与对照组的基因表达数据进行差异分析，筛选出差异表达基因（DEGs）。以|log2FC|>1且FDR<0.05作为筛选条件，得到了不同组间的差异表达基因数量及变化趋势。与对照组相比，高毒力PRRSV感染组在感染后第3天有1286个基因表达显著差异，其中上调基因892个，下调基因394个；第7天差异表达基因数量增加至2568个，上调基因1678个，下调基因890个；第14天差异表达基因数量为2154个，上调基因1320个，下调基因834个。这表明随着感染时间的延长，高毒力PRRSV感染组猪肺组织中基因表达的变化更为显著，且上调基因的数量始终多于下调基因，说明高毒力PRRSV感染可能导致猪肺组织中大量基因的表达被激活，从而引发一系列生物学过程的改变。低毒力PRRSV感染组与对照组相比，在感染后第3天有765个基因表达差异显著，其中上调基因486个，下调基因279个；第7天差异表达基因数量为1352个，上调基因895个，下调基因457个；第14天差异表达基因数量为1023个，上调基因650个，下调基因373个。低毒力PRRSV感染组的差异表达基因数量整体上少于高毒力PRRSV感染组，且在各时间点上调基因数量也相对较少，这提示低毒力PRRSV感染对猪肺组织基因表达的影响程度相对较弱，引发的生物学过程变化不如高毒力PRRSV感染明显。通过对不同毒力PRRSV感染组与对照组间差异表达基因的筛选和统计分析，直观地呈现了不同毒力PRRSV感染后猪肺组织基因表达的变化情况，为后续深入研究不同毒力PRRSV感染的分子机制奠定了基础，有助于进一步挖掘与PRRSV致病相关的关键基因和生物学通路。3.3差异表达基因的功能注释为深入探究不同毒力PRRSV感染后猪肺组织中差异表达基因（DEGs）的生物学功能，运用生物信息学方法，通过BLAST软件将筛选出的DEGs核苷酸序列或氨基酸序列与多个数据库进行比对，对其进行全面的功能注释，并利用clusterProfiler软件开展GO功能富集分析和KEGG通路富集分析。在GO功能富集分析中，从生物过程（BiologicalProcess）、细胞组成（CellularComponent）和分子功能（MolecularFunction）三个层面进行剖析。在生物过程方面，高毒力PRRSV感染组的DEGs显著富集于免疫应答、炎症反应、细胞因子介导的信号通路等过程。例如，大量与T细胞活化、B细胞增殖以及细胞因子分泌相关的基因表达上调，表明高毒力PRRSV感染强烈激活了猪肺组织的免疫反应和炎症进程。低毒力PRRSV感染组的DEGs也富集于免疫相关过程，但程度相对较弱，且还涉及一些细胞代谢调节过程，如碳水化合物代谢调控等，说明低毒力PRRSV感染在引发一定免疫反应的同时，对细胞代谢也产生了影响。从细胞组成角度来看，高毒力PRRSV感染组的DEGs主要富集于细胞膜、细胞外基质、溶酶体等细胞组成部分。其中，与细胞膜相关的基因表达变化可能影响病毒的入侵和感染过程，而细胞外基质相关基因的改变可能与肺组织的病理损伤和修复有关。低毒力PRRSV感染组的DEGs在细胞组成上的富集主要集中在细胞核、线粒体等细胞器，提示低毒力PRRSV感染对细胞内基本结构和功能的影响具有一定的特异性。在分子功能层面，高毒力PRRSV感染组的DEGs富集于细胞因子活性、趋化因子活性、蛋白激酶活性等分子功能。这些功能的改变与免疫应答和炎症反应的激活密切相关，例如细胞因子和趋化因子活性的增强有助于招募免疫细胞到感染部位，参与免疫防御。低毒力PRRSV感染组的DEGs则在核酸结合、氧化还原酶活性等分子功能上富集，反映了低毒力PRRSV感染后细胞内核酸代谢和氧化还原平衡的变化。通过KEGG通路富集分析，确定了DEGs显著富集的代谢通路和信号转导通路。高毒力PRRSV感染组的DEGs显著富集于Toll样受体信号通路、NF-κB信号通路、MAPK信号通路等。Toll样受体信号通路的激活可启动机体的天然免疫反应，识别PRRSV等病原体相关分子模式；NF-κB信号通路和MAPK信号通路的活化则进一步促进炎症因子的表达和释放，引发强烈的炎症反应，这与高毒力PRRSV感染后猪肺组织出现的严重炎症病变相吻合。低毒力PRRSV感染组的DEGs除了在部分免疫相关信号通路有一定富集外，还在一些代谢通路如糖酵解/糖异生通路、脂肪酸代谢通路等有富集。这表明低毒力PRRSV感染在影响免疫反应的同时，对猪肺组织的能量代谢和物质代谢也产生了较为明显的影响。通过GO和KEGG数据库注释及富集分析，全面揭示了不同毒力PRRSV感染后猪肺组织中DEGs在生物过程、细胞组成和分子功能方面的作用，以及它们参与的关键代谢通路和信号转导通路，为深入理解PRRSV的致病机制提供了重要线索。3.4差异表达基因的富集分析3.4.1GO富集分析利用clusterProfiler软件对不同毒力PRRSV感染组与对照组间的差异表达基因（DEGs）进行GO富集分析，从生物过程（BiologicalProcess）、细胞组成（CellularComponent）和分子功能（MolecularFunction）三个层面深入探究DEGs的生物学功能和潜在作用机制。在生物过程方面，高毒力PRRSV感染组在感染后第3天，DEGs显著富集于“免疫应答激活”（GO:0002250）、“炎症反应调节”（GO:0050727）、“细胞因子介导的信号通路”（GO:0019221）等条目。其中，“免疫应答激活”条目下包含众多与T细胞、B细胞活化相关的基因，如CD3D、CD3E、CD40LG等，这些基因的上调表明高毒力PRRSV感染早期就强烈刺激了机体的免疫反应，促使免疫细胞活化并参与免疫防御。“炎症反应调节”条目中，IL1B、IL6、TNF等炎症因子相关基因表达上调，显示炎症反应在感染初期就已被启动。随着感染时间延长至第7天，“抗病毒免疫反应”（GO:0051607）、“白细胞迁移调节”（GO:0030593）等条目显著富集，表明机体在持续对抗病毒感染的过程中，进一步加强了抗病毒免疫反应，同时通过调节白细胞迁移，增强免疫细胞对感染部位的浸润和清除病原体的能力。到第14天，“组织修复”（GO:0001501）、“细胞凋亡调控”（GO:0042981）等条目富集明显，这可能是由于感染后期肺组织损伤后启动自我修复过程，同时对受损细胞进行凋亡调控，以维持组织的正常结构和功能。低毒力PRRSV感染组在第3天，DEGs主要富集于“免疫应答启动”（GO:0002251）、“细胞代谢过程调节”（GO:0031323）等条目。“免疫应答启动”相关基因的表达变化相对较弱，说明低毒力PRRSV感染初期引发的免疫反应强度低于高毒力感染。而“细胞代谢过程调节”条目中涉及碳水化合物代谢、脂肪酸代谢等相关基因的表达改变，提示低毒力PRRSV感染早期就对细胞代谢产生影响。在第7天，“免疫细胞活化调节”（GO:0002694）、“氧化还原过程”（GO:0055114）等条目显著富集，表明随着感染的进行，机体在调节免疫细胞活化以对抗病毒的同时，细胞内氧化还原平衡也受到影响。第14天，“蛋白质折叠”（GO:0006457）、“RNA代谢过程”（GO:0016070）等条目富集，说明低毒力PRRSV感染后期对细胞内蛋白质和核酸的代谢过程产生了较为明显的作用。在细胞组成层面，高毒力PRRSV感染组在第3天，DEGs主要富集于“细胞膜”（GO:0005886）、“细胞外基质”（GO:0031012）、“溶酶体”（GO:0005764）等细胞组成部分。其中，细胞膜相关基因的变化可能影响病毒的吸附、侵入和释放过程；细胞外基质相关基因的改变可能与肺组织的结构完整性和病理损伤修复有关；溶酶体相关基因的表达变化可能参与病毒的降解和免疫调节。第7天，“线粒体”（GO:0005739）、“内质网”（GO:0005783）等细胞器相关条目富集，表明病毒感染对细胞内能量代谢和蛋白质合成加工等重要细胞器的功能产生影响。第14天，“细胞核”（GO:0005634）、“细胞连接”（GO:0030054）等条目显著富集，可能与基因转录调控和细胞间通讯在感染后期的变化有关。低毒力PRRSV感染组在第3天，DEGs主要富集于“细胞核”（GO:0005634）、“线粒体”（GO:0005739）等细胞组成，提示低毒力PRRSV感染早期就对细胞内重要的遗传物质储存和能量代谢场所产生影响。第7天，“细胞骨架”（GO:0005856）、“高尔基体”（GO:0005794）等条目富集，表明病毒感染对细胞的形态维持、物质运输和加工等功能产生作用。第14天，“质膜”（GO:0005887）、“细胞外囊泡”（GO:0070062）等条目显著富集，可能与细胞与外界环境的物质交换和信号传递在感染后期的改变有关。在分子功能方面，高毒力PRRSV感染组在第3天，DEGs富集于“细胞因子活性”（GO:0005125）、“趋化因子活性”（GO:0008009）、“蛋白激酶活性”（GO:0004672）等分子功能。细胞因子和趋化因子活性的增强有助于招募免疫细胞到感染部位，启动免疫应答；蛋白激酶活性的改变则可能参与细胞内信号转导通路的调控，影响细胞的增殖、分化和凋亡等过程。第7天，“受体活性”（GO:0004872）、“酶结合”（GO:0019899）等条目显著富集，表明病毒感染导致细胞表面受体和酶活性的变化，进而影响细胞对病毒的识别和感染后的代谢过程。第14天，“转录因子活性”（GO:0003700）、“核酸结合”（GO:0003676）等条目富集，说明感染后期基因转录调控和核酸代谢过程受到显著影响。低毒力PRRSV感染组在第3天，DEGs主要富集于“核酸结合”（GO:0003676）、“氧化还原酶活性”（GO:0016491）等分子功能，反映了低毒力PRRSV感染早期对细胞内核酸代谢和氧化还原平衡的影响。第7天，“转运蛋白活性”（GO:0005215）、“蛋白质结合”（GO:0005515）等条目显著富集，表明病毒感染对细胞内物质运输和蛋白质相互作用产生作用。第14天，“ATP结合”（GO:0005524）、“水解酶活性”（GO:0016787）等条目富集，说明低毒力PRRSV感染后期对细胞的能量代谢和物质分解过程产生影响。通过GO富集分析，全面揭示了不同毒力PRRSV感染后猪肺组织中DEGs在生物过程、细胞组成和分子功能方面的显著变化，为深入理解PRRSV的致病机制提供了丰富的生物学信息。3.4.2KEGG富集分析为进一步解析不同毒力PRRSV感染后猪肺组织中差异表达基因（DEGs）参与的重要生物学通路，运用clusterProfiler软件对DEGs进行KEGG通路富集分析，深入探究病毒感染引发的细胞内代谢和信号转导等过程的变化。高毒力PRRSV感染组在感染后第3天，DEGs显著富集于Toll样受体信号通路（ko04620）、NF-κB信号通路（ko04064）、MAPK信号通路（ko04010）等。Toll样受体信号通路是机体天然免疫的重要组成部分，通过识别病原体相关分子模式，激活下游信号转导，启动免疫应答。在高毒力PRRSV感染早期，该通路中关键基因如TLR3、TLR4、MyD88等表达上调，表明高毒力PRRSV能够被Toll样受体有效识别，从而激活天然免疫反应。NF-κB信号通路和MAPK信号通路在免疫调节和炎症反应中发挥关键作用。感染早期这两条通路的活化，促使炎症因子如IL1B、IL6、TNF等的表达和释放，引发强烈的炎症反应，这与高毒力PRRSV感染后猪肺组织出现的严重炎症病变相吻合。随着感染时间延长至第7天，除上述免疫相关信号通路持续富集外，细胞凋亡通路（ko04210）、自噬相关通路（ko04140）也显著富集。细胞凋亡是机体清除受损或感染细胞的重要机制，高毒力PRRSV感染后期诱导细胞凋亡相关基因如BAX、CASP3等表达上调，表明病毒感染可能导致肺组织细胞凋亡增加，以限制病毒复制和扩散。自噬是细胞内的一种自我降解和更新机制，在病毒感染过程中，自噬既可以协助病毒复制，也可能参与病毒的清除。高毒力PRRSV感染第7天自噬相关通路的富集，提示自噬在病毒感染后期的复杂作用，可能涉及病毒与宿主细胞之间的相互博弈。到第14天，代谢相关通路如糖酵解/糖异生通路（ko00010）、脂肪酸代谢通路（ko01212）等也出现显著富集。这表明高毒力PRRSV感染后期，肺组织的能量代谢和物质代谢发生明显改变。糖酵解/糖异生通路的变化可能是为了满足感染细胞对能量的需求增加；脂肪酸代谢通路的改变则可能影响细胞膜的合成和修复，以及炎症介质的产生，进一步影响肺组织的病理变化和免疫反应。低毒力PRRSV感染组在第3天，DEGs除了在部分免疫相关信号通路如Toll样受体信号通路有一定程度的富集外，还在一些代谢通路如糖酵解/糖异生通路、脂肪酸代谢通路等有明显富集。这表明低毒力PRRSV感染早期在引发一定免疫反应的同时，就对猪肺组织的能量代谢和物质代谢产生了较为明显的影响。相较于高毒力PRRSV感染，低毒力感染早期引发的免疫反应相对较弱，但对代谢通路的影响出现时间较早。在第7天，低毒力PRRSV感染组的DEGs在吞噬体通路（ko04145）、内吞作用通路（ko04144）等有显著富集。吞噬体通路和内吞作用是细胞摄取和处理病原体的重要途径，这表明低毒力PRRSV感染中期，机体通过加强吞噬和内吞作用来清除病毒，同时也可能影响病毒在细胞内的运输和复制过程。此外，代谢相关通路如氨基酸代谢通路（ko01230）也持续富集，说明低毒力PRRSV感染对细胞内物质代谢的影响在不断持续和深化。第14天，低毒力PRRSV感染组的DEGs在细胞粘附分子通路（ko04514）、细胞周期通路（ko04110）等有富集。细胞粘附分子通路的变化可能影响免疫细胞与肺组织细胞之间的相互作用，以及炎症细胞的浸润和迁移。细胞周期通路的富集则提示低毒力PRRSV感染后期对细胞的增殖和生长调控产生影响，可能与肺组织的修复和再生过程有关。通过KEGG富集分析，明确了不同毒力PRRSV感染后猪肺组织中DEGs显著富集的代谢通路和信号转导通路，揭示了病毒感染在免疫应答、炎症反应、细胞凋亡、自噬、能量代谢以及细胞增殖与生长调控等多个方面对肺组织的影响，为深入理解PRRSV的致病机制提供了关键线索。3.5转录因子调控网络分析转录因子（TranscriptionFactors，TFs）在基因表达调控中起着核心作用，它们通过与基因启动子区域的特定DNA序列结合，调控基因的转录起始和转录速率，从而影响细胞的生理功能和生物学过程。在病毒感染过程中，宿主细胞内的转录因子调控网络会发生复杂的变化，这些变化与病毒的致病机制密切相关。为深入探究不同毒力PRRSV感染后猪肺组织中基因表达调控的分子机制，本研究基于转录组学数据，对差异表达基因（DEGs）进行了转录因子调控网络分析。首先，利用生物信息学工具，从筛选出的DEGs中识别出潜在的转录因子。通过与已知的转录因子数据库（如AnimalTFDB等）进行比对，结合转录因子结合位点（TranscriptionFactorBindingSites，TFBSs）预测算法，确定了在不同毒力PRRSV感染组中差异表达的转录因子。结果显示，在高毒力PRRSV感染组中，发现了多个显著差异表达的转录因子，如NF-κB、AP-1、STAT1等。其中，NF-κB是一种重要的转录因子，在免疫应答和炎症反应中发挥关键作用。在高毒力PRRSV感染后，NF-κB的表达显著上调，提示其可能被激活并参与调控一系列与免疫和炎症相关基因的表达，这与KEGG富集分析中NF-κB信号通路的显著富集结果相呼应。AP-1也是一个在免疫和应激反应中起重要作用的转录因子，其表达变化可能影响细胞的增殖、分化和凋亡等过程。STAT1参与细胞因子介导的信号通路，在抗病毒免疫反应中具有重要功能，高毒力PRRSV感染后STAT1的差异表达表明其可能在宿主抗病毒免疫应答中发挥关键调控作用。在低毒力PRRSV感染组中，也鉴定出一些差异表达的转录因子，如SP1、E2F1等。SP1是一种广泛表达的转录因子，参与细胞的基本生理过程，如细胞增殖、分化和代谢等。低毒力PRRSV感染后SP1的表达变化可能影响细胞内相关基因的表达，进而对细胞的正常功能产生影响。E2F1在细胞周期调控中发挥重要作用，其表达差异可能与低毒力PRRSV感染后细胞周期的改变有关，这与KEGG富集分析中细胞周期通路的富集结果相一致。为了进一步揭示转录因子与靶基因之间的调控关系，构建了转录因子调控网络。利用Cytoscape软件，将识别出的转录因子和它们的靶基因（即差异表达基因）作为节点，转录因子与靶基因之间的调控关系作为边，构建了可视化的调控网络。在高毒力PRRSV感染组的转录因子调控网络中，NF-κB处于网络的核心位置，与众多免疫和炎症相关基因相连，如IL1B、IL6、TNF等。这表明NF-κB可能通过调控这些基因的表达，在高毒力PRRSV感染引发的强烈免疫应答和炎症反应中发挥关键的调控作用。AP-1和STAT1也与多个基因存在相互作用，它们可能协同NF-κB或独立地参与调控不同的生物学过程，共同影响高毒力PRRSV感染后的细胞命运和生理功能。在低毒力PRRSV感染组的转录因子调控网络中，SP1与一些参与细胞代谢和增殖相关的基因相互作用，如与糖酵解/糖异生通路中的关键酶基因相连，这进一步支持了低毒力PRRSV感染对细胞代谢产生影响的结论。E2F1则主要与细胞周期相关基因相互关联，表明其在低毒力PRRSV感染后细胞周期调控中的重要作用。通过对转录因子调控网络的拓扑结构分析，计算了网络的度（Degree）、中介中心性（BetweennessCentrality）和接近中心性（ClosenessCentrality）等指标。度表示节点与其他节点之间的连接数量，度值越高，说明该节点在网络中的连接越广泛，可能在调控网络中发挥更重要的作用。中介中心性反映了节点在网络中信息传递的能力，中介中心性高的节点往往处于网络的关键路径上，对网络中其他节点之间的信息交流和调控起着桥梁作用。接近中心性衡量节点与网络中其他所有节点的接近程度，接近中心性高的节点能够快速地将信息传递到整个网络。在高毒力PRRSV感染组的转录因子调控网络中，NF-κB的度、中介中心性和接近中心性均较高，表明其在网络中具有重要的地位，是调控网络的核心节点。AP-1和STAT1等转录因子也具有相对较高的这些指标值，说明它们在网络中也发挥着重要的调控作用。在低毒力PRRSV感染组的转录因子调控网络中，SP1和E2F1的这些指标相对较高，显示出它们在低毒力PRRSV感染相关的基因调控网络中的关键作用。通过转录因子调控网络分析，明确了不同毒力PRRSV感染后猪肺组织中差异表达的转录因子及其与靶基因之间的调控关系，揭示了转录因子在PRRSV致病过程中的重要调控作用。这些结果为深入理解PRRSV的致病机制提供了新的视角，有助于进一步挖掘潜在的治疗靶点和防控策略。四、不同毒力PRRSVs感染后肺组织的蛋白质组学分析4.1蛋白质鉴定与定量结果运用基于液相色谱-串联质谱（LC-MS/MS）的蛋白质组学技术，对不同毒力PRRSV感染组及对照组猪肺组织样本中的蛋白质进行分离、鉴定和定量分析。通过将质谱采集到的数据与猪的蛋白质数据库（如UniProt猪数据库）进行比对，利用SequestHT算法进行蛋白质鉴定，最终在所有样本中成功鉴定出[X]种蛋白质。在蛋白质定量方面，采用Label-free定量方法，通过比较不同样本中肽段的峰面积或离子强度，计算蛋白质的相对表达量。对鉴定出的蛋白质进行定量统计分析，结果显示各样本中蛋白质的定量重复性良好。以高毒力PRRSV感染组的三个生物学重复样本为例，计算各样本中相同蛋白质的定量相关性，Pearson相关系数均大于0.90，表明实验重复性高，定量结果可靠，能够准确反映不同样本中蛋白质的表达差异。对不同毒力PRRSV感染组与对照组之间的蛋白质表达数据进行初步比较分析，发现随着感染时间的延长，蛋白质表达的差异逐渐明显。在感染后第3天，高毒力PRRSV感染组与对照组相比，已有部分蛋白质的表达出现显著差异；到第7天，差异表达蛋白质的数量进一步增加；第14天，蛋白质表达的差异更为显著。低毒力PRRSV感染组也呈现出类似的趋势，但在相同感染时间点，其差异表达蛋白质的数量和变化幅度均低于高毒力PRRSV感染组。这些蛋白质鉴定与定量结果为后续筛选差异表达蛋白质、深入分析蛋白质功能和参与的生物学通路奠定了坚实基础，有助于全面揭示不同毒力PRRSV感染后猪肺组织在蛋白质水平的变化特征，为阐明PRRSV的致病机制提供重要的蛋白质组学数据支持。4.2差异表达蛋白筛选在明确蛋白质鉴定与定量结果后，运用Perseus软件对高毒力PRRSV感染组、低毒力PRRSV感染组与对照组的蛋白质表达数据展开差异分析，筛选出差异表达蛋白质（DEPs）。设定筛选条件为|log2FC|>1且p-value<0.05，以确保筛选出的差异表达蛋白具有统计学意义和显著的表达变化。与对照组相比，高毒力PRRSV感染组在感染后第3天有186个蛋白质表达显著差异，其中上调蛋白112个，下调蛋白74个；第7天差异表达蛋白数量上升至328个，上调蛋白205个，下调蛋白123个；第14天差异表达蛋白数量为276个，上调蛋白168个，下调蛋白108个。随着感染时间的推移，高毒力PRRSV感染组中差异表达蛋白的数量先增加后略有减少，但整体数量较多，且上调蛋白始终多于下调蛋白，这表明高毒力PRRSV感染对猪肺组织蛋白质表达的影响较为剧烈，在感染过程中持续激活或抑制大量蛋白质的表达，从而深刻改变肺组织的生理功能和生物学过程。低毒力PRRSV感染组与对照组相比，在感染后第3天有105个蛋白质表达差异显著，其中上调蛋白68个，下调蛋白37个；第7天差异表达蛋白数量为172个，上调蛋白108个，下调蛋白64个；第14天差异表达蛋白数量为134个，上调蛋白85个，下调蛋白49个。低毒力PRRSV感染组的差异表达蛋白数量明显少于高毒力PRRSV感染组，且在各时间点上调蛋白数量相对较少，这表明低毒力PRRSV感染对猪肺组织蛋白质表达的影响程度较轻，引发的蛋白质水平变化不如高毒力PRRSV感染显著。通过对不同毒力PRRSV感染组与对照组间差异表达蛋白的筛选和统计，清晰地呈现了不同毒力PRRSV感染后猪肺组织蛋白质表达的动态变化特征，为后续深入剖析差异表达蛋白的功能和参与的生物学通路奠定了基础，有助于进一步揭示PRRSV致病的分子机制，挖掘与PRRSV感染相关的潜在生物标志物和治疗靶点。4.3差异表达蛋白的功能注释为深入剖析不同毒力PRRSV感染后猪肺组织中差异表达蛋白（DEPs）的生物学功能，借助生物信息学手段，运用BLAST软件将筛选出的DEPs氨基酸序列与多个权威数据库进行比对，全面开展功能注释，并利用clusterProfiler软件进行GO功能富集分析和KEGG通路富集分析。在GO功能富集分析中，从生物过程、细胞组成和分子功能三个层面展开研究。在生物过程方面，高毒力PRRSV感染组的DEPs显著富集于免疫应答、炎症反应调节、细胞内物质转运等过程。其中，与免疫细胞活化、趋化因子分泌相关的蛋白，如CD4、CCL2等表达上调，表明高毒力PRRSV感染强烈激活了猪肺组织的免疫反应和炎症进程，通过吸引免疫细胞向感染部位聚集，增强免疫防御能力。参与细胞内物质转运的蛋白表达变化，可能影响病毒在细胞内的运输和装配过程。低毒力PRRSV感染组的DEPs也富集于免疫相关过程，但程度相对较弱，且还涉及细胞代谢调节、蛋白质折叠与修饰等过程。例如，一些参与碳水化合物代谢和脂肪酸代谢的酶类蛋白表达改变，提示低毒力PRRSV感染在引发一定免疫反应的同时，对细胞代谢产生了影响。同时，与蛋白质折叠和修饰相关的分子伴侣蛋白和修饰酶的表达变化，可能影响蛋白质的正常功能和稳定性。从细胞组成角度来看，高毒力PRRSV感染组的DEPs主要富集于细胞膜、线粒体、溶酶体等细胞组成部分。细胞膜相关蛋白的变化可能影响病毒的吸附、侵入和释放过程，线粒体相关蛋白的改变可能影响细胞的能量代谢和凋亡调控，溶酶体相关蛋白的表达变化则可能参与病毒的降解和免疫调节。低毒力PRRSV感染组的DEPs在细胞核、内质网、细胞骨架等细胞组成上有富集。细胞核相关蛋白的表达变化可能影响基因转录和调控，内质网相关蛋白的改变可能影响蛋白质合成和加工，细胞骨架相关蛋白的变化则可能影响细胞的形态和运动。在分子功能层面，高毒力PRRSV感染组的DEPs富集于细胞因子活性、受体活性、酶活性调节等分子功能。细胞因子活性的增强有助于招募免疫细胞到感染部位，启动免疫应答；受体活性的改变则可能影响细胞对病毒的识别和感染后的信号转导过程；酶活性调节相关蛋白的表达变化可能参与细胞内信号转导通路的调控，影响细胞的增殖、分化和凋亡等过程。低毒力PRRSV感染组的DEPs则在核酸结合、氧化还原酶活性、转运蛋白活性等分子功能上富集。核酸结合蛋白的表达变化可能影响基因表达调控，氧化还原酶活性的改变可能影响细胞内氧化还原平衡和代谢过程，转运蛋白活性的变化则可能影响细胞内物质的运输和交换。通过KEGG通路富集分析，确定了DEPs显著富集的代谢通路和信号转导通路。高毒力PRRSV感染组的DEPs显著富集于Toll样受体信号通路、NF-κB信号通路、MAPK信号通路等。这些通路在免疫调节和炎症反应中发挥关键作用，高毒力PRRSV感染后这些通路的活化，促使炎症因子如IL1B、IL6、TNF等的表达和释放，引发强烈的炎症反应，这与高毒力PRRSV感染后猪肺组织出现的严重炎症病变相吻合。此外，高毒力PRRSV感染组的DEPs还富集于细胞凋亡通路、自噬相关通路等。细胞凋亡通路的激活可能是机体清除受损或感染细胞的一种防御机制，自噬相关通路的变化则可能涉及病毒与宿主细胞之间的相互作用和博弈。低毒力PRRSV感染组的DEPs除了在部分免疫相关信号通路有一定富集外，还在一些代谢通路如糖酵解/糖异生通路、脂肪酸代谢通路等有明显富集。这表明低毒力PRRSV感染在影响免疫反应的同时，对猪肺组织的能量代谢和物质代谢也产生了较为明显的影响。此外，低毒力PRRSV感染组的DEPs在吞噬体通路、内吞作用通路等有显著富集。吞噬体通路和内吞作用是细胞摄取和处理病原体的重要途径，这表明低毒力PRRSV感染中期，机体通过加强吞噬和内吞作用来清除病毒，同时也可能影响病毒在细胞内的运输和复制过程。通过GO和KEGG数据库注释及富集分析，全面揭示了不同毒力PRRSV感染后猪肺组织中DEPs在生物过程、细胞组成和分子功能方面的作用，以及它们参与的关键代谢通路和信号转导通路，为深入理解PRRSV的致病机制提供了重要线索。4.4差异表达蛋白的富集分析4.4.1GO富集分析运用clusterProfiler软件对不同毒力PRRSV感染组与对照组间的差异表达蛋白（DEPs）展开GO富集分析，从生物过程（BiologicalProcess）、细胞组成（CellularComponent）和分子功能（MolecularFunction）三个维度深入剖析DEPs的生物学功能和潜在作用机制。在生物过程层面，高毒力PRRSV感染组在感染后第3天，DEPs显著富集于“免疫应答激活”（GO:0002250）、“炎症反应调节”（GO:0050727）、“细胞内蛋白质运输”（GO:0006886）等条目。“免疫应答激活”条目下，如CD4、CCL2等蛋白表达上调，表明高毒力PRRSV感染初期就强烈激发了机体的免疫反应，促进免疫细胞活化并向感染部位趋化。“炎症反应调节”条目中，IL1B、IL6等炎症因子相关蛋白表达上调，预示着炎症反应在感染早期就已启动。“细胞内蛋白质运输”相关蛋白的表达变化，可能影响病毒在细胞内的运输、装配和释放过程。到第7天，“抗病毒免疫反应”（GO:0051607）、“白细胞迁移调节”（GO:0030593）等条目显著富集，表明机体在持续对抗病毒感染的进程中，进一步强化了抗病毒免疫反应，同时通过调节白细胞迁移，增强免疫细胞对感染部位的浸润和病原体清除能力。第14天，“组织修复”（GO:0001501）、“细胞凋亡调控”（GO:0042981）等条目富集明显，这可能是由于感染后期肺组织损伤后启动自我修复过程，同时对受损细胞进行凋亡调控，以维持组织的正常结构和功能。低毒力PRRSV感染组在第3天，DEPs主要富集于“免疫应答启动”（GO:0002251）、“细胞代谢过程调节”（GO:0031323）等条目。“免疫应答启动”相关蛋白的表达变化相对较弱，说明低毒力PRRSV感染初期引发的免疫反应强度低于高毒力感染。而“细胞代谢过程调节”条目中涉及碳水化合物代谢、脂肪酸代谢等相关酶类蛋白的表达改变，提示低毒力PRRSV感染早期就对细胞代谢产生影响。在第7天，“免疫细胞活化调节”（GO:0002694）、“氧化还原过程”（GO:0055114）等条目显著富集，表明随着感染的进行，机体在调节免疫细胞活化以对抗病毒的同时，细胞内氧化还原平衡也受到影响。第14天，“蛋白质折叠”（GO:0006457）、“RNA代谢过程”（GO:0016070）等条目富集，说明低毒力PRRSV感染后期对细胞内蛋白质和核酸的代谢过程产生了较为明显的作用。从细胞组成角度来看，高毒力PRRSV感染组在第3天，DEPs主要富集于“细胞膜”（GO:0005886）、“线粒体”（GO:0005739）、“溶酶体”（GO:0005764）等细胞组成部分。细胞膜相关蛋白的变化可能影响病毒的吸附、侵入和释放过程；线粒体相关蛋白的改变可能影响细胞的能量代谢和凋亡调控，为病毒感染提供能量支持或影响细胞的死亡程序；溶酶体相关蛋白的表达变化则可能参与病毒的降解和免疫调节。第7天，“内质网”（GO:0005783）、“细胞骨架”（GO:0005856）等细胞器相关条目富集，表明病毒感染对细胞内蛋白质合成加工和细胞形态维持等重要细胞器的功能产生影响。第14天，“细胞核”（GO:0005634）、“细胞连接”（GO:0030054）等条目显著富集，可能与基因转录调控和细胞间通讯在感染后期的变化有关。低毒力PRRSV感染组在第3天，DEPs主要富集于“细胞核”（GO:0005634）、“线粒体”（GO:0005739）等细胞组成，提示低毒力PRRSV感染早期就对细胞内重要的遗传物质储存和能量代谢场所产生影响。第7天，“内质网”（GO:0005783）、“高尔基体”（GO:0005794）等条目富集，表明病毒感染对细胞的蛋白质合成加工和物质运输加工等功能产生作用。第14天，“质膜”（GO:0005887）、“细胞外囊泡”（GO:0070062）等条目显著富集，可能与细胞与外界环境的物质交换和信号传递在感染后期的改变有关。在分子功能方面，高毒力PRRSV感染组在第3天，DEPs富集于“细胞因子活性”（GO:0005125）、“趋化因子活性”（GO:0008009）、“蛋白激酶活性”（GO:0004672）等分子功能。细胞因子和趋化因子活性的增强有助于招募免疫细胞到感染部位，启动免疫应答；蛋白激酶活性的改变则可能参与细胞内信号转导通路的调控，影响细胞的增殖、分化和凋亡等过程。第7天，“受体活性”（GO:0004872）、“酶结合”（GO:0019899）等条目显著富集，表明病毒感染导致细胞表面受体和酶活性的变化，进而影响细胞对病毒的识别和感染后的代谢过程。第14天，“转录因子活性”（GO:0003700）、“核酸结合”（GO:0003676）等条目富集，说明感染后期基因转录调控和核酸代谢过程受到显著影响。低毒力PRRSV感染组在第3天，DEPs主要富集于“核酸结合”（GO:0003676）、“氧化还原酶活性”（GO:0016491）等分子功能，反映了低毒力PRRSV感染早期对细胞内核酸代谢和氧化还原平衡的影响。第7天，“转运蛋白活性”（GO:0005215）、“蛋白质结合”（GO:0005515）等条目显著富集，表明病毒感染对细胞内物质运输和蛋白质相互作用产生作用。第14天，“ATP结合”（GO:0005524）、“水解酶活性”（GO:0016787）等条目富集，说明低毒力PRRSV感染后期对细胞的能量代谢和物质分解过程产生影响。通过GO富集分析，全面揭示了不同毒力PRRSV感染后猪肺组织中DEPs在生物过程、细胞组成和分子功能方面的显著变化，为深入理解PRRSV的致病机制提供了丰富的生物学信息。4.4.2KEGG富集分析为进一步解析不同毒力PRRSV感染后猪肺组织中差异表达蛋白（DEPs）参与的重要生物学通路，运用clusterProfiler软件对DEPs进行KEGG通路富集分析，深入探究病毒感染引发的细胞内代谢和信号转导等过程的变化。高毒力PRRSV感染组在感染后第3天，DEPs显著富集于Toll样受体信号通路（ko04620）、NF-κB信号通路（ko04064）、MAPK信号通路（ko04010）等。Toll样受体信号通路是机体天然免疫的重要组成部分，通过识别病原体相关分子模式，激活下游信号转导，启动免疫应答。在高毒力PRRSV感染早期，该通路中关键蛋白如TLR3、TLR4、MyD88等表达上调，表明高毒力PRRSV能够被Toll样受体有效识别，从而激活天然免疫反应。NF-κB信号通路和MAPK信号通路在免疫调节和炎症反应中发挥关键作用。感染早期这两条通路的活化，促使炎症因子如IL1B、IL6、TNF等的表达和释放，引发强烈的炎症反应，这与高毒力PRRSV感染后猪肺组织出现的严重炎症病变相吻合。随着感染时间延长至第7天，除上述免疫相关信号通路持续富集外，细胞凋亡通路（ko04210）、自噬相关通路（ko04140）也显著富集。细胞凋亡是机体清除受损或感染细胞的重要机制，高毒力PRRSV感染后期诱导细胞凋亡相关蛋白如BAX、CASP3等表达上调，表明病毒感染可能导致肺组织细胞凋亡增加，以限制病毒复制和扩散。自噬是细胞内的一种自我降解和更新机制，在病毒感染过程中，自噬既可以协助病毒复制，也可能参与病毒的清除。高毒力PRRSV感染第7天自噬相关通路的富集，提示自噬在病毒感染后期的复杂作用，可能涉及病毒与宿主细胞之间的相互博弈。到第14天，代谢相关通路如糖酵解/糖异生通路（ko00010）、脂肪酸代谢通路（ko01212）等也出现显著富集。这表明高毒力PRRSV感染后期，肺组织的能量代谢和物质代谢发生明显改变。糖酵解/糖异生通路的变化可能是为了满足感染细胞对能量的需求增加；脂肪酸代谢通路的改变则可能影响细胞膜的合成和修复，以及炎症介质的产生，进一步影响肺组织的病理变化和免疫反应。低毒力PRRSV感染组在第3天，DEPs除了在部分免疫相关信号通路如Toll样受体信号通路有一定程度的富集外，还在一些代谢通路如糖酵解/糖异生通路、脂肪酸代谢通路等有明显富集。这表明低毒力PRRSV感染早期在引发一定免疫反应的同时，就对猪肺组织的能量代谢和物质代谢产生了较为明显的影响。相较于高毒力PRRSV感染，低毒力感染早期引发的免疫反应相对较弱，但对代谢通路的影响出现时间较早。在第7天，低毒力PRRSV感染组的DEPs在吞噬体通路（ko04145）、内吞作用通路（ko04144）等有显著富集。吞噬体通路和内吞作用是细胞摄取和处理病原体的重要途径，这表明低毒力PRRSV感染中期，机体通过加强吞噬和内吞作用来清除病毒，同时也可能影响病毒在细胞内的运输和复制过程。此外，代谢相关通路如氨基酸代谢通路（ko01230）也持续富集，说明低毒力PRRSV感染对细胞内物质代谢的影响在不断持续和深化。第14天，低毒力PRRSV感染组的DEPs在细胞粘附分子通路（ko04514）、细胞周期通路（ko04110）等有富集。细胞粘附分子通路的变化可能影响免疫细胞与肺组织细胞之间的相互作用，以及炎症细胞的浸润和迁移。细胞周期通路的富集则提示低毒力PRRSV感染后期对细胞的增殖和生长调控产生影响，可能与肺组织的修复和再生过程有关。通过KEGG富集分