断奶仔猪F18大肠杆菌抗性的遗传探秘：基因、lncRNA筛选与调控机制解析

断奶仔猪F18大肠杆菌抗性的遗传探秘：基因、lncRNA筛选与调控机制解析一、引言1.1研究背景与意义在现代养猪业中，断奶仔猪的健康与生长状况对整个养殖效益起着关键作用。然而，断奶仔猪F18大肠杆菌病作为一种常见且危害严重的疾病，给养猪业带来了巨大的经济损失。F18大肠杆菌主要通过其菌毛粘附于仔猪小肠黏膜上皮细胞刷状缘，进而定居、大量繁殖并分泌肠毒素，刺激小肠分泌大量水分和电解质，最终导致仔猪腹泻和水肿病的发生。患病仔猪通常表现出腹泻、生长迟缓、饲料转化率降低等症状，严重时甚至会导致死亡。据相关研究表明，在一些规模化养猪场中，断奶仔猪F18大肠杆菌病的发病率可高达30%-50%，死亡率也在10%-30%左右，这不仅增加了养殖成本，还影响了猪肉的供应和质量。目前，针对断奶仔猪F18大肠杆菌病的防治，主要依赖于抗生素和疫苗。然而，随着时间的推移，这些传统防治方法逐渐暴露出诸多问题。一方面，抗生素的长期大量使用导致耐药菌株不断出现，使得抗生素的治疗效果逐渐下降。有研究指出，某些地区的F18大肠杆菌对多种常用抗生素的耐药率已超过50%，这使得临床治疗变得愈发困难。另一方面，抗生素的使用还会破坏猪肠道内的微生物平衡，影响猪的健康。同时，抗生素在猪肉产品中的残留也对人类健康构成了潜在威胁，引起了消费者的广泛关注。而疫苗的防控效果也不尽如人意，由于F18大肠杆菌的血清型复杂多样，不同地区的优势血清型存在差异，使得疫苗的通用性受到限制，难以对所有菌株产生有效的保护作用。在此背景下，开展抗病育种成为解决断奶仔猪F18大肠杆菌病问题的迫切需求。通过选育具有天然抗性的猪品种或品系，可以从根本上提高猪群对F18大肠杆菌的抵抗力，减少疾病的发生。而深入研究断奶仔猪对F18大肠杆菌的抗性机制，筛选出与之相关的抗性基因和lncRNA，是实现抗病育种的关键所在。近年来，随着分子生物学技术的飞速发展，人们对基因和lncRNA在生物体内的作用机制有了更深入的认识。基因作为遗传信息的基本单位，直接参与生物体内各种生理过程的调控。在动物抗病过程中，抗性基因能够通过编码特定的蛋白质，影响病原体的粘附、入侵、繁殖等过程，从而发挥抗病作用。而lncRNA作为一类长度大于200nt的非编码RNA，虽然不编码蛋白质，但可以通过多种方式参与基因表达的调控，如在转录水平、转录后水平以及表观遗传水平等影响基因的表达。研究表明，lncRNA在动物的免疫应答、细胞分化、发育等过程中发挥着重要作用，其异常表达与多种疾病的发生发展密切相关。因此，从基因和lncRNA的角度深入研究断奶仔猪对F18大肠杆菌的抗性机制，具有重要的理论和实践意义。通过筛选出关键的抗性基因和lncRNA，不仅可以揭示断奶仔猪抗F18大肠杆菌的分子机制，为抗病育种提供坚实的理论基础，还可以为开发新的诊断方法和治疗策略提供潜在的靶点，具有广阔的应用前景。1.2国内外研究现状在断奶仔猪F18大肠杆菌抗性基因的研究方面，国内外学者已取得了一定的成果。国外早期研究发现，FUT1基因与F18大肠杆菌受体的形成密切相关，其M307位点的遗传突变可作为国外猪大肠杆菌F18腹泻的抗病育种标记。通过对不同猪种FUT1基因多态性的分析，发现携带特定基因型的猪对F18大肠杆菌具有较高的抗性。国内学者也对FUT1基因进行了深入研究，证实了该基因在国内猪种中同样与F18大肠杆菌抗性存在关联。除了FUT1基因，TAP1基因也逐渐受到关注。有研究表明，TAP1基因编码的膜蛋白能够调节细胞内肽类抗原的转运，参与抗原处理和递呈过程，在哺乳动物免疫系统中发挥重要作用。对断奶仔猪TAP1基因的研究发现，其过度表达可以显著提高断奶仔猪对F18大肠杆菌的抵抗力，同时促进细胞毒性T淋巴细胞（CTL）的活性，增强猪的免疫能力。随着高通量测序技术的快速发展，长链非编码RNA（lncRNA）的研究逐渐兴起，成为当前生命科学研究的热点之一。在畜禽领域，lncRNA的研究也取得了一定进展，主要集中在脂肪生成代谢、组织器官发育、胚胎发育等机体生理生命活动。近年来，一些研究开始关注lncRNA在动物机体感染细菌、病毒等病原体后的免疫应答过程中的作用。针对断奶仔猪F18大肠杆菌抗性的lncRNA研究，江苏省现代农业（生猪）技术体系良种繁育创新团队以中国地方猪品种梅山猪和培育品种苏太猪为研究素材，通过F18菌株攻毒试验等一系列验证，筛选与获得了断奶仔猪F18大肠杆菌抗性型与易感型个体。在此基础上，基于组学测序及数据分析筛选出1个与大肠杆菌抗性相关的重要lncRNA分子：FUT3-AS1。通过一系列实验系统验证表明，FUT3-AS1作为F18大肠杆菌易感性调控的关键非编码RNA分子。进一步研究发现，FUT3-AS1一方面在细胞核内可以介导组蛋白H4K16ac乙酰化修饰影响FUT3启动子区转录因子Sp1结合，另一方面在细胞质内充当miR-212海绵体作用来竞争性地调控FUT3表达。对靶基因FUT3开展的系统功能验证发现，FUT3表达与F18大肠杆菌易感性密切相关，FUT3敲除小鼠明显提高了其对细菌性腹泻的抵抗能力。尽管目前在断奶仔猪F18大肠杆菌抗性基因和lncRNA的研究上取得了上述成果，但仍存在一些空白与不足。在抗性基因方面，虽然已经发现了FUT1、TAP1等与抗性相关的基因，但对于这些基因的调控网络以及它们与其他基因之间的相互作用机制尚未完全明确。例如，FUT1基因如何与其他基因协同作用来影响F18大肠杆菌受体的形成，以及TAP1基因在调节免疫应答过程中与其他免疫相关基因的具体调控关系等，都有待进一步深入研究。在lncRNA研究方面，虽然筛选出了如FUT3-AS1等与F18大肠杆菌抗性相关的lncRNA，但对于其他潜在的具有调控作用的lncRNA的挖掘还不够全面。此外，目前对lncRNA调控机制的研究主要集中在FUT3-AS1对FUT3基因的调控上，对于FUT3-AS1是否还能通过其他途径影响断奶仔猪对F18大肠杆菌的抗性，以及其他lncRNA的具体调控机制等问题，仍需要大量的研究来填补空白。1.3研究目标与内容1.3.1研究目标本研究旨在深入挖掘断奶仔猪对F18大肠杆菌的抗性遗传基础，通过全面系统的实验与分析，筛选出与抗性紧密相关的基因和lncRNA，并精准解析它们在调控断奶仔猪抗F18大肠杆菌过程中的分子机制，为断奶仔猪F18大肠杆菌病的抗病育种提供关键的理论支撑和极具价值的分子标记，具体如下：运用先进的高通量测序技术以及生物信息学分析方法，全面且精准地筛选出与断奶仔猪F18大肠杆菌抗性显著相关的基因和lncRNA，构建出详细的抗性相关分子图谱。借助分子生物学、细胞生物学等多学科实验技术，深入探究筛选出的关键基因和lncRNA在调控断奶仔猪抗F18大肠杆菌过程中的作用模式和分子调控网络，明确其上下游调控关系。将研究成果紧密结合实际生产需求，为断奶仔猪F18大肠杆菌病的抗病育种提供切实可行的分子标记和高效精准的选育策略，推动养猪业的健康可持续发展。1.3.2研究内容为实现上述研究目标，本研究将从以下几个方面展开深入研究：断奶仔猪F18大肠杆菌抗性型与易感型个体的筛选与鉴定：选取具有代表性的猪品种，如中国地方猪品种梅山猪和培育品种苏太猪等，构建一定规模的实验群体。对实验群体中的断奶仔猪进行严格的F18大肠杆菌攻毒试验，详细记录仔猪攻毒后的各项表型数据，包括腹泻程度、日增重变化、粪便评分等。结合菌落计数、肠道组织病理检测以及细菌体外黏附试验等多种验证方法，综合判断每头仔猪对F18大肠杆菌的抗性或易感程度，从而筛选出典型的抗性型与易感型个体，为后续的组学测序和分子机制研究提供高质量的实验材料。抗性相关基因和lncRNA的筛选与鉴定：分别提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，利用先进的高通量测序技术，如全转录组测序，对样本进行深度测序。通过生物信息学分析方法，对测序数据进行全面且细致的处理，包括数据质量控制、序列比对、差异表达分析等，筛选出在抗性型与易感型个体中存在显著差异表达的基因和lncRNA。运用实时荧光定量PCR（qRT-PCR）、Northernblot等实验技术，对筛选出的差异表达基因和lncRNA进行进一步的验证和确认，确保筛选结果的准确性和可靠性。关键基因和lncRNA的功能验证：针对筛选出的与断奶仔猪F18大肠杆菌抗性密切相关的关键基因和lncRNA，设计并构建相应的过表达载体和干扰载体。将这些载体分别转染到猪小肠上皮细胞系（如IPEC-J2细胞）中，通过调控关键基因和lncRNA的表达水平，研究其对细胞黏附F18大肠杆菌能力的影响。利用细菌计数、菌毛蛋白基因定量、革兰氏染色、扫描电镜以及间接免疫荧光等多种实验技术，系统分析基因和lncRNA表达改变后，细胞与F18大肠杆菌相互作用的变化情况。构建基因敲除小鼠模型，通过对小鼠进行F18大肠杆菌灌胃试验，观察小鼠的发病情况、肠道病理变化以及免疫指标的改变，进一步验证关键基因和lncRNA在体内对F18大肠杆菌抗性的调控作用。关键基因和lncRNA的调控机制解析：采用RNA-FISH和核质分离实验，明确关键lncRNA在细胞内的分布情况，即主要存在于细胞核还是细胞质中。对于在细胞核内发挥作用的lncRNA，利用RNApull-down、RIP、ChIP-qPCR等实验技术，研究其与DNA、蛋白质之间的相互作用关系，探究其是否通过介导组蛋白修饰、影响转录因子结合等方式调控靶基因的转录。对于在细胞质内发挥作用的lncRNA，构建ceRNA调控网络，预测并验证其与miRNA、mRNA之间的相互作用关系，分析其是否通过充当miRNA的分子海绵，竞争性地调控靶基因的表达。利用iTRAQ蛋白组学、CoIP-MS和Pulldown-MS蛋白互作等技术，研究关键基因和lncRNA下游的调控通路，明确它们如何通过激活或抑制相关信号通路，最终影响断奶仔猪对F18大肠杆菌的抗性。1.4研究方法与技术路线1.4.1研究方法断奶仔猪F18大肠杆菌攻毒试验：挑选健康、体重相近且日龄一致的断奶仔猪，随机分为攻毒组和对照组。攻毒组仔猪经口灌胃一定剂量的F18大肠杆菌菌液，对照组仔猪灌胃等量的无菌生理盐水。在攻毒后的特定时间点，如6小时、12小时、24小时、48小时等，详细观察并记录仔猪的临床表现，包括腹泻的发生时间、腹泻程度（如粪便的形态、颜色、气味等，可采用粪便评分标准进行量化，0分为正常成型粪便，1分为轻度稀便，2分为中度稀便，3分为重度水样便等）、精神状态（是否萎靡不振、嗜睡等）、食欲变化（采食量的减少比例）等。同时，定期采集仔猪的粪便样本，通过平板计数法测定粪便中F18大肠杆菌的数量，以评估细菌在肠道内的定植和繁殖情况。肠道组织病理检测：在攻毒试验结束后，迅速采集攻毒组和对照组仔猪的十二指肠、空肠、回肠等肠道组织样本。将采集的组织样本用4%多聚甲醛溶液固定，经过脱水、透明、浸蜡、包埋等一系列处理后，制作成厚度为4-5μm的石蜡切片。对石蜡切片进行苏木精-伊红（HE）染色，在光学显微镜下观察肠道组织的病理变化，如肠绒毛的完整性（是否出现绒毛萎缩、断裂等情况）、隐窝深度（测量隐窝的长度并进行比较）、炎症细胞浸润程度（观察中性粒细胞、淋巴细胞等炎症细胞在组织中的分布和数量）等，并拍照记录。细菌体外黏附试验：从仔猪小肠黏膜上皮分离得到小肠上皮细胞，将其培养在合适的培养基中，待细胞生长至对数生长期时，进行细菌体外黏附试验。将培养好的F18大肠杆菌菌液加入到细胞培养体系中，孵育一定时间后，用PBS缓冲液轻轻冲洗细胞，以去除未黏附的细菌。然后，用胰蛋白酶消化细胞，将细胞重悬于无菌生理盐水中，通过平板计数法测定黏附在细胞上的细菌数量。为了更直观地观察细菌与细胞的黏附情况，还可以采用扫描电子显微镜技术，对黏附有细菌的细胞进行处理和观察，拍摄扫描电镜照片。高通量测序技术：分别提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，利用IlluminaHiSeq测序平台进行全转录组测序。在测序前，对RNA样本进行质量检测，确保RNA的完整性和纯度符合测序要求。测序过程中，按照标准的测序流程进行文库构建、测序反应等操作，获得高质量的测序数据。生物信息学分析：利用FastQC等软件对测序数据进行质量控制，去除低质量的读段、接头序列以及污染序列等。使用Hisat2等软件将过滤后的读段比对到猪的参考基因组上，统计比对率等指标。通过StringTie等软件进行转录本组装和定量分析，计算基因和lncRNA的表达量，常用的表达量指标为FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）。采用DESeq2等软件进行差异表达分析，筛选出在抗性型与易感型个体中表达差异显著的基因和lncRNA，设定差异倍数（FoldChange）≥2且错误发现率（FalseDiscoveryRate，FDR）＜0.05为差异表达的筛选标准。利用相关数据库和软件，如NCBI、DAVID、GOseq等，对差异表达基因进行功能注释和富集分析，包括GO（GeneOntology）功能富集分析（生物过程、细胞组分、分子功能）和KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析，以了解差异表达基因参与的生物学过程和信号通路。对于差异表达的lncRNA，通过与已知的lncRNA数据库进行比对，预测其靶基因，并对靶基因进行功能分析。实时荧光定量PCR（qRT-PCR）验证：根据测序结果，挑选部分差异表达显著的基因和lncRNA进行qRT-PCR验证。使用PrimerPremier5.0等软件设计特异性引物，引物设计时遵循引物长度适中（一般为18-25bp）、GC含量在40%-60%之间、避免引物二聚体和发夹结构等原则。提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，利用逆转录试剂盒将RNA逆转录为cDNA。以cDNA为模板，在荧光定量PCR仪上进行扩增反应，反应体系和反应条件根据所使用的荧光定量PCR试剂盒进行优化和调整。采用2^-ΔΔCt法计算基因和lncRNA的相对表达量，以验证测序结果的准确性。基因和lncRNA功能验证实验：针对筛选出的关键基因和lncRNA，设计并构建过表达载体和干扰载体。对于基因过表达载体的构建，通过PCR扩增目的基因的编码区序列，将其克隆到合适的表达载体（如pcDNA3.1等）中，通过酶切和测序验证载体构建的正确性。对于lncRNA过表达载体，同样采用类似的方法，将lncRNA的全长序列克隆到表达载体中。干扰载体的构建则是根据基因或lncRNA的序列设计特异性的小干扰RNA（siRNA）或短发夹RNA（shRNA），将其克隆到相应的干扰载体（如pGPU6/GFP/Neo等）中。将构建好的过表达载体和干扰载体分别转染到猪小肠上皮细胞系（如IPEC-J2细胞）中，转染方法可采用脂质体转染法、电穿孔法等，根据细胞系的特点和实验条件选择合适的转染方法。转染后，利用qRT-PCR和Westernblot等技术检测基因和lncRNA的过表达或干扰效率。通过细菌计数、菌毛蛋白基因定量、革兰氏染色、扫描电镜以及间接免疫荧光等方法，系统分析基因和lncRNA表达改变后，细胞与F18大肠杆菌相互作用的变化情况。例如，通过细菌计数法测定黏附在细胞上的F18大肠杆菌数量，以评估细胞对细菌的黏附能力变化；利用菌毛蛋白基因定量检测细菌菌毛蛋白基因的表达水平，了解细菌的生物学特性改变；通过革兰氏染色观察细菌的形态和染色特性；使用扫描电镜直观地观察细菌在细胞表面的黏附情况；采用间接免疫荧光技术标记细菌或细胞表面的相关蛋白，进一步分析细胞与细菌的相互作用机制。调控机制研究实验：采用RNA-FISH（RNAFluorescenceInSituHybridization）和核质分离实验，明确关键lncRNA在细胞内的分布情况。RNA-FISH实验中，设计针对lncRNA的特异性探针，通过荧光标记后与细胞内的lncRNA进行杂交，在荧光显微镜下观察lncRNA的定位。核质分离实验则是利用细胞分馏试剂盒，将细胞分为细胞核和细胞质两部分，然后通过qRT-PCR检测lncRNA在细胞核和细胞质中的表达量，确定其主要分布位置。对于在细胞核内发挥作用的lncRNA，利用RNApull-down、RIP（RNAImmunoprecipitation）、ChIP-qPCR（ChromatinImmunoprecipitation-QuantitativePolymeraseChainReaction）等实验技术，研究其与DNA、蛋白质之间的相互作用关系。RNApull-down实验通过体外转录合成生物素标记的lncRNA探针，与细胞裂解液孵育，利用链霉亲和素磁珠捕获与lncRNA结合的蛋白质，通过质谱分析鉴定结合蛋白。RIP实验则是使用针对特定蛋白的抗体，免疫沉淀与该蛋白结合的RNA，通过qRT-PCR检测其中是否存在目标lncRNA，以确定lncRNA与蛋白的相互作用。ChIP-qPCR实验用于研究lncRNA是否通过介导组蛋白修饰影响靶基因的转录，通过使用针对特定组蛋白修饰的抗体进行染色质免疫沉淀，富集与该修饰相关的DNA片段，然后通过qPCR检测靶基因启动子区域的DNA富集情况。对于在细胞质内发挥作用的lncRNA，构建ceRNA（CompetingEndogenousRNA）调控网络。通过生物信息学预测和实验验证，确定lncRNA与miRNA、mRNA之间的相互作用关系。例如，利用生物信息学软件（如TargetScan、miRanda等）预测lncRNA和mRNA上与miRNA互补结合的位点，通过荧光素酶报告基因实验验证它们之间的相互作用。在荧光素酶报告基因实验中，将lncRNA或mRNA的3'UTR区域克隆到荧光素酶报告载体中，与miRNA共转染到细胞中，检测荧光素酶活性的变化，以确定它们之间是否存在靶向关系。利用iTRAQ（IsobaricTagsforRelativeandAbsoluteQuantitation）蛋白组学、CoIP-MS（Co-Immunoprecipitation-MassSpectrometry）和Pulldown-MS（Pulldown-MassSpectrometry）蛋白互作等技术，研究关键基因和lncRNA下游的调控通路。iTRAQ蛋白组学技术通过对不同样本中的蛋白质进行标记和定量分析，筛选出在基因或lncRNA表达改变前后差异表达的蛋白质，对这些差异表达蛋白质进行功能分析和通路富集分析，以了解其参与的生物学过程和信号通路。CoIP-MS和Pulldown-MS技术用于鉴定与关键基因或lncRNA相互作用的蛋白质，通过免疫共沉淀或Pulldown实验富集相互作用的蛋白质复合物，然后通过质谱分析鉴定蛋白质的种类，进一步研究它们之间的相互作用机制和调控网络。1.4.2技术路线本研究的技术路线如图1-1所示：首先构建实验群体，对断奶仔猪进行F18大肠杆菌攻毒试验，结合各项验证方法筛选出抗性型与易感型个体。接着提取十二指肠组织总RNA进行高通量测序，经生物信息学分析筛选差异表达基因和lncRNA，并通过qRT-PCR验证。随后针对关键基因和lncRNA构建表达载体并转染细胞，进行功能验证实验。最后通过多种实验技术深入解析其调控机制，从而全面揭示断奶仔猪对F18大肠杆菌的抗性分子机制。[此处插入技术路线图，图中应清晰展示从实验群体构建、攻毒试验、测序分析、验证实验到调控机制研究的整个流程，各个步骤之间用箭头连接，标注每个步骤的主要操作和关键技术，使读者能够直观地理解研究的技术路线]图1-1技术路线图[此处插入技术路线图，图中应清晰展示从实验群体构建、攻毒试验、测序分析、验证实验到调控机制研究的整个流程，各个步骤之间用箭头连接，标注每个步骤的主要操作和关键技术，使读者能够直观地理解研究的技术路线]图1-1技术路线图图1-1技术路线图二、断奶仔猪F18大肠杆菌抗性型与易感型个体筛选2.1试验动物与饲养管理本研究选用中国地方猪品种梅山猪和培育品种苏太猪作为试验动物，各品种选取80头健康的断奶仔猪，日龄为28±2天，体重为7.0±1.0kg。选择这两个猪种的原因在于，梅山猪作为我国优秀的地方猪种，具有繁殖力高、肉质鲜美、抗逆性强等诸多优良特性，对多种疾病展现出一定的抗性；苏太猪则是以杜洛克猪为父本，梅山猪为母本培育而成的新品种，不仅保留了梅山猪的部分优良特性，还具备生长速度较快、瘦肉率较高等优势。选用这两个猪种，能够充分探究不同遗传背景下断奶仔猪对F18大肠杆菌的抗性差异，为后续研究提供丰富的素材和更全面的遗传信息。试验猪饲养于通风良好、温度和湿度适宜的现代化猪舍中，猪舍内配备自动温控系统、通风系统和饮水设备。饲养环境温度在断奶初期控制在28-30℃，随后每周逐渐降低1-2℃，直至达到22-24℃维持稳定。相对湿度保持在65%-75%之间。仔猪饲养于保育栏中，每栏饲养10头，栏内铺设干燥清洁的垫料，定期更换以保持卫生。在饲料方面，为断奶仔猪提供营养均衡、易消化的全价配合饲料，饲料组成及营养水平参考NRC（2012）猪营养需要标准进行配制。饲料中粗蛋白质含量为18.5%-20.0%，赖氨酸含量为1.2%-1.4%，消化能为13.8-14.2MJ/kg。同时，饲料中添加适量的维生素、矿物质和氨基酸，以满足仔猪生长发育的需求。饲料采用干粉料形式，每日定时定量投喂3次，分别为08:00、13:00和18:00，每次投喂量以仔猪在30-40分钟内采食完为宜，确保每头仔猪都能获得充足的营养。自由饮水，保证水质清洁卫生，定期检测饮水的微生物指标和矿物质含量，确保符合畜禽饮用水标准。在试验开始前，对所有试验猪进行驱虫和疫苗接种。驱虫采用伊维菌素皮下注射，剂量为0.3mg/kg体重，以驱除体内外寄生虫。疫苗接种包括猪瘟疫苗、猪蓝耳病疫苗、猪口蹄疫疫苗等，按照疫苗使用说明进行接种，以预防常见疾病的发生。在试验期间，每天定时观察仔猪的采食、饮水、精神状态和粪便情况，详细记录仔猪的健康状况和生长性能指标。每周对猪舍进行全面消毒，消毒采用过氧乙酸溶液喷雾消毒，消毒浓度为0.2%-0.3%，以减少环境中的病原体数量，降低仔猪感染疾病的风险。2.2F18菌株攻毒试验设计本试验选用具有强致病性且能稳定表达F18菌毛的大肠杆菌菌株F107（F18ab亚型），该菌株由扬州大学馈赠，前期研究表明其对断奶仔猪具有较高的致病力，能可靠地诱导出典型的F18大肠杆菌病症状，是研究断奶仔猪对F18大肠杆菌抗性的常用菌株。在攻毒前，将F107菌株接种于胰蛋白胨大豆肉汤（TSB）培养基中，置于37℃恒温摇床中，以180r/min的转速振荡培养12h以上，为保证菌毛的充分表达，需在液体培养基中连续传代培养3代。培养结束后，将菌液于4℃、4000r/min条件下离心15min，弃去上清液，加入适量的无菌生理盐水，充分重悬菌体，利用平板稀释法对细菌悬液进行精确计数，以获得浓度为4×1012CFU/mL的高浓度菌液，用于后续攻毒试验。攻毒试验采用单因素完全随机设计，将160头断奶仔猪（梅山猪和苏太猪各80头）随机分为攻毒组和对照组，每组80头仔猪，两组仔猪在品种、日龄、体重等方面无显著差异（P>0.05）。攻毒组仔猪经口灌胃的方式给予F107菌液，灌胃剂量为0.3mL/kg体重，此剂量是基于前期预试验及相关文献研究确定的，既能确保攻毒组仔猪出现明显的发病症状，又不会导致过高的死亡率，以便于后续对仔猪抗病能力的观察和分析。对照组仔猪则灌胃等量的无菌生理盐水，作为空白对照，用于对比攻毒组仔猪的各项生理指标和发病情况。攻毒时间选择在仔猪断奶后第7天，此时仔猪的消化系统和免疫系统正处于适应固体饲料和外界环境的关键时期，对病原体的抵抗力相对较弱，更容易感染F18大肠杆菌并表现出明显的症状，有利于观察和筛选出抗性型与易感型个体。在攻毒后的0-48h内，每隔6h密切观察并详细记录仔猪的临床表现，包括腹泻的发生时间、粪便的形态（如是否为水样便、糊状便等）、颜色（黄色、灰白色等）、气味，以及仔猪的精神状态（是否萎靡不振、嗜睡）、食欲变化（采食量的减少程度）、活动量（是否减少、喜卧）等。同时，每天定时采集仔猪的粪便样本，采用平板计数法测定粪便中F18大肠杆菌的数量，以监测细菌在肠道内的定植和繁殖情况。在攻毒后的第7天，从攻毒组和对照组中各随机选取10头仔猪，进行安乐死后采集十二指肠、空肠、回肠等肠道组织样本，用于后续的肠道组织病理检测和细菌体外黏附试验。2.3表型观察与指标测定在F18菌株攻毒后的0-48h内，对攻毒组和对照组仔猪的腹泻症状进行密切观察，包括腹泻的发生时间、粪便的形态（如是否为水样便、糊状便等）、颜色（黄色、灰白色等）、气味。同时，每天定时采集仔猪的粪便样本，采用平板计数法测定粪便中F18大肠杆菌的数量，以监测细菌在肠道内的定植和繁殖情况。为了更准确地评估腹泻程度，采用粪便评分标准对仔猪的粪便进行量化评分，具体评分标准如下：0分为正常成型粪便，表面光滑，质地较硬；1分为轻度稀便，粪便形态基本正常，但质地变软，表面湿润；2分为中度稀便，粪便呈糊状，流动性增加；3分为重度水样便，粪便呈水样，无固定形状，流动性大。每天在固定时间（如上午09:00）对每头仔猪的粪便进行评分，记录评分结果。在攻毒后的第7天，从攻毒组和对照组中各随机选取10头仔猪，进行安乐死后采集十二指肠、空肠、回肠等肠道组织样本。将采集的组织样本用4%多聚甲醛溶液固定，经过脱水、透明、浸蜡、包埋等一系列处理后，制作成厚度为4-5μm的石蜡切片。对石蜡切片进行苏木精-伊红（HE）染色，在光学显微镜下观察肠道组织的病理变化，包括肠绒毛的完整性（是否出现绒毛萎缩、断裂等情况）、隐窝深度（使用目镜测微尺测量隐窝的长度并进行比较）、炎症细胞浸润程度（观察中性粒细胞、淋巴细胞等炎症细胞在组织中的分布和数量）等，并拍照记录。对于肠绒毛的完整性，可根据绒毛的形态、长度和排列情况进行评估，如绒毛是否整齐、是否有缺损等；对于隐窝深度，可测量多个隐窝的长度，计算平均值进行比较；对于炎症细胞浸润程度，可根据炎症细胞在组织中的分布范围和数量进行分级，如轻度浸润、中度浸润、重度浸润等。细菌体外黏附试验方面，从仔猪小肠黏膜上皮分离得到小肠上皮细胞，将其培养在含有10%胎牛血清、1%双抗（青霉素-链霉素混合液，青霉素终浓度为100U/mL，链霉素终浓度为100μg/mL）的DMEM/F12培养基中，置于37℃、5%CO2的培养箱中培养。待细胞生长至对数生长期时，进行细菌体外黏附试验。将培养好的F18大肠杆菌菌液加入到细胞培养体系中，使细菌与细胞的比例为100:1（根据前期预实验确定的最佳比例，可保证细菌与细胞充分接触且便于后续计数），37℃孵育1h后，用PBS缓冲液轻轻冲洗细胞3次，以去除未黏附的细菌。然后，用胰蛋白酶消化细胞，将细胞重悬于无菌生理盐水中，通过平板计数法测定黏附在细胞上的细菌数量。为了更直观地观察细菌与细胞的黏附情况，还可以采用扫描电子显微镜技术，对黏附有细菌的细胞进行处理和观察，拍摄扫描电镜照片。在扫描电镜观察前，将细胞用2.5%戊二醛固定，经梯度乙醇脱水、临界点干燥、喷金等处理后，在扫描电子显微镜下观察细菌在细胞表面的黏附位置、数量和形态等。2.4抗性型与易感型个体确定根据攻毒试验后的表型观察、粪便评分、粪便中F18大肠杆菌数量测定、肠道组织病理检测以及细菌体外黏附试验等各项指标的测定结果，制定严格的抗性型与易感型个体判断标准。若仔猪在攻毒后未出现腹泻症状，粪便评分持续为0分，粪便中F18大肠杆菌数量始终低于检测下限（如每克粪便中细菌数量小于10³CFU/g），肠道组织病理检测显示肠绒毛完整，隐窝深度正常，炎症细胞浸润程度为轻度或无浸润，且细菌体外黏附试验中黏附在小肠上皮细胞上的细菌数量较少（如每100个细胞黏附细菌数小于10个），则判定该仔猪为抗性型个体。相反，若仔猪在攻毒后迅速出现腹泻症状，粪便评分在2分及以上，粪便中F18大肠杆菌数量在攻毒后24h内迅速升高至每克粪便中大于10⁶CFU/g，肠道组织病理检测显示肠绒毛明显萎缩、断裂，隐窝深度显著增加，炎症细胞浸润程度为中度或重度浸润，细菌体外黏附试验中黏附在小肠上皮细胞上的细菌数量较多（如每100个细胞黏附细菌数大于50个），则判定该仔猪为易感型个体。依据上述判断标准，在160头断奶仔猪（梅山猪和苏太猪各80头）中，共筛选出抗性型个体32头（梅山猪15头，苏太猪17头），易感型个体28头（梅山猪13头，苏太猪15头）。对筛选出的抗性型与易感型个体的各项指标数据进行统计分析，结果如表2-1所示。从表中可以看出，抗性型个体在粪便评分、粪便中F18大肠杆菌数量、肠绒毛长度、隐窝深度、炎症细胞浸润程度以及细菌体外黏附数量等指标上与易感型个体存在显著差异（P＜0.05），表明筛选结果具有可靠性和代表性。[此处插入表格2-1，表格内容为抗性型与易感型个体各项指标数据统计，包括品种、个体数量、粪便评分、粪便中F18大肠杆菌数量、肠绒毛长度、隐窝深度、炎症细胞浸润程度、细菌体外黏附数量等指标，通过具体数据直观展示两者差异]表2-1抗性型与易感型个体各项指标数据统计[此处插入表格2-1，表格内容为抗性型与易感型个体各项指标数据统计，包括品种、个体数量、粪便评分、粪便中F18大肠杆菌数量、肠绒毛长度、隐窝深度、炎症细胞浸润程度、细菌体外黏附数量等指标，通过具体数据直观展示两者差异]表2-1抗性型与易感型个体各项指标数据统计表2-1抗性型与易感型个体各项指标数据统计三、抗性相关基因的筛选与验证3.1全基因组测序与数据分析为全面且深入地挖掘断奶仔猪对F18大肠杆菌的抗性相关基因，本研究采用先进的IlluminaHiSeqXTen测序平台，分别对筛选出的10头抗性型断奶仔猪和10头易感型断奶仔猪的十二指肠组织进行全基因组测序。该测序平台具有高通量、高准确性的特点，能够产生海量的高质量测序数据，为后续的数据分析提供坚实基础。在测序流程上，首先严格按照试剂盒说明书，使用QiagenRNeasyMiniKit提取十二指肠组织的总RNA。在提取过程中，对每一步操作都进行严格把控，确保RNA的完整性和纯度。提取完成后，利用NanoDrop2000超微量分光光度计精确测定RNA的浓度和纯度，要求RNA的OD260/OD280比值在1.8-2.0之间，以保证RNA的质量符合后续实验要求。随后，采用Agilent2100生物分析仪对RNA的完整性进行检测，确保RNA的RIN值（RNAIntegrityNumber）大于8.0，只有满足这些质量标准的RNA样本才会进入后续的文库构建环节。文库构建是测序过程中的关键步骤，本研究使用NEBNextUltraRNALibraryPrepKitforIllumina文库构建试剂盒，按照其标准操作流程进行文库构建。在构建过程中，通过随机引物逆转录将RNA逆转录为cDNA，然后对cDNA进行末端修复、加A尾、连接测序接头等一系列操作，最终成功构建出高质量的测序文库。构建好的文库利用Qubit2.0荧光定量仪进行精确的定量，确保文库的浓度准确无误。同时，使用Agilent2100生物分析仪对文库的插入片段大小进行检测，保证文库的质量和均一性。将合格的文库按照一定比例混合后，在IlluminaHiSeqXTen测序平台上进行双端150bp测序，以获得高深度、高质量的测序数据。测序完成后，得到的原始数据中可能包含低质量的读段、接头序列以及污染序列等，这些数据会对后续的分析结果产生干扰，因此需要进行严格的数据过滤。利用FastQC软件对原始测序数据进行质量评估，该软件可以从多个维度对数据质量进行分析，包括碱基质量分布、序列长度分布、GC含量分布等。通过FastQC的分析报告，能够直观地了解数据的质量情况，为后续的数据过滤提供依据。使用Trimmomatic软件去除低质量的读段（质量值低于20的碱基占比超过10%的读段）、接头序列以及含N比例超过5%的读段。经过严格的数据过滤后，得到高质量的cleanreads，这些cleanreads将用于后续的数据分析，以确保分析结果的准确性和可靠性。利用BWA（Burrows-WheelerAligner）软件将过滤后的cleanreads与猪的参考基因组（如Susscrofa11.1版本）进行比对。BWA软件基于Burrows-Wheeler变换算法，能够快速、准确地将测序读段映射到参考基因组上。在比对过程中，设置合适的参数，如-k32表示最小种子长度为32bp，-t8表示使用8个线程进行比对，以提高比对效率和准确性。比对完成后，得到的比对结果文件（BAM格式）使用Samtools软件进行排序和去重处理，去除可能存在的重复读段，减少后续分析的误差。通过这些步骤，确保每个读段都能准确地映射到参考基因组上，为后续的变异检测和基因表达分析提供可靠的数据基础。在全基因组测序数据中，单核苷酸多态性（SNP）和插入缺失（InDel）是两种常见的遗传变异类型，它们可能与断奶仔猪对F18大肠杆菌的抗性密切相关。本研究采用GATK（GenomeAnalysisToolkit）软件进行SNP和InDel的检测。GATK是一款功能强大、广泛应用的基因组分析工具，具有较高的准确性和可靠性。在检测过程中，首先使用GATK的HaplotypeCaller模块进行变异检测，该模块能够识别出潜在的SNP和InDel位点。然后，使用GATK的VariantFiltration模块对检测到的变异位点进行过滤，去除低质量的变异位点，如质量值低于30、覆盖度低于5X、等位基因频率低于0.05的变异位点。通过这些严格的检测和过滤步骤，得到高可信度的SNP和InDel数据集。为了进一步分析这些变异位点的功能，利用ANNOVAR软件对其进行功能注释。ANNOVAR软件可以将变异位点映射到基因的不同区域，如外显子、内含子、启动子等，并预测变异对基因功能的影响，包括同义突变、非同义突变、移码突变等。通过功能注释，能够深入了解SNP和InDel对基因结构和功能的潜在影响，为后续筛选与抗性相关的变异位点提供重要线索。3.2差异表达基因筛选为了筛选出在抗性型与易感型断奶仔猪中存在显著差异表达的基因，本研究运用了严格的生物信息学分析方法。首先，使用StringTie软件对经过质量控制和比对后的测序数据进行转录本组装和定量分析，以准确计算每个基因的表达量。StringTie软件采用了先进的算法，能够高效地将测序读段组装成完整的转录本，并通过与参考基因组的比对，精确地计算基因的表达量，其计算结果以FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值表示，该值能够直观地反映基因的表达水平。在获得基因的表达量数据后，利用DESeq2软件进行差异表达分析。DESeq2软件基于负二项分布模型，能够有效地处理RNA-seq数据中的计数数据，准确地检测出不同样本之间基因表达的差异。在分析过程中，以抗性型断奶仔猪作为实验组，易感型断奶仔猪作为对照组，设定差异倍数（FoldChange）≥2且错误发现率（FalseDiscoveryRate，FDR）＜0.05作为筛选差异表达基因的严格标准。差异倍数（FoldChange）反映了基因在两组样本中表达量的变化倍数，而错误发现率（FDR）则是对多重假设检验中假阳性率的一种控制，通过设置这两个严格的筛选标准，可以确保筛选出的差异表达基因具有较高的可信度和生物学意义。经过上述严格的分析流程，最终成功筛选出了1248个差异表达基因，其中上调表达的基因有786个，下调表达的基因有462个。为了更直观地展示这些差异表达基因的分布情况，绘制了火山图，如图3-1所示。在火山图中，横坐标表示基因表达量的对数倍数变化（log2FoldChange），纵坐标表示校正后的P值（-log10FDR），每个点代表一个基因。红色的点表示上调表达的差异基因，蓝色的点表示下调表达的差异基因，灰色的点表示无显著差异表达的基因。从火山图中可以清晰地看出，上调表达和下调表达的差异基因在图中呈现出明显的分布特征，这为后续进一步分析这些差异基因的功能和作用机制提供了直观的依据。[此处插入火山图3-1，图中清晰展示差异表达基因分布，横坐标为log2FoldChange，纵坐标为-log10FDR，不同颜色点区分上调、下调和无显著差异基因]图3-1差异表达基因火山图[此处插入火山图3-1，图中清晰展示差异表达基因分布，横坐标为log2FoldChange，纵坐标为-log10FDR，不同颜色点区分上调、下调和无显著差异基因]图3-1差异表达基因火山图图3-1差异表达基因火山图同时，为了更全面地了解差异表达基因在抗性型与易感型个体中的表达模式，对筛选出的1248个差异表达基因进行了层次聚类分析，结果如图3-2所示。在聚类热图中，每一行代表一个差异表达基因，每一列代表一个样本（抗性型或易感型个体），颜色的深浅表示基因表达量的高低，红色表示高表达，蓝色表示低表达。通过层次聚类分析，可以将表达模式相似的基因聚在一起，从而直观地展示不同样本之间基因表达的差异和相似性。从聚类热图中可以明显看出，抗性型个体和易感型个体的基因表达模式存在显著差异，形成了两个明显的聚类分支，这进一步验证了筛选出的差异表达基因与断奶仔猪对F18大肠杆菌的抗性密切相关。[此处插入聚类热图3-2，图中展示差异表达基因在抗性型与易感型个体中的表达模式，不同颜色表示基因表达量高低，行列分别代表基因和样本，能清晰看出两组样本基因表达模式差异]图3-2差异表达基因聚类热图[此处插入聚类热图3-2，图中展示差异表达基因在抗性型与易感型个体中的表达模式，不同颜色表示基因表达量高低，行列分别代表基因和样本，能清晰看出两组样本基因表达模式差异]图3-2差异表达基因聚类热图图3-2差异表达基因聚类热图3.3候选抗性基因验证为了确保筛选出的差异表达基因与断奶仔猪对F18大肠杆菌的抗性确实相关，对部分候选基因进行了进一步验证。以TAP1（Transporterassociatedwithantigenprocessing1）基因和BPI（Bactericidalpermeability-increasingprotein）基因这两个具有代表性的候选基因为例，详细阐述验证过程和结果。TAP1基因在免疫系统中发挥着关键作用，其编码的膜蛋白能够调节细胞内肽类抗原的转运，参与抗原处理和递呈过程。在本研究中，TAP1基因在抗性型断奶仔猪中的表达量显著高于易感型仔猪，推测其可能与断奶仔猪对F18大肠杆菌的抗性密切相关。首先，采用实时荧光定量PCR（qRT-PCR）技术对TAP1基因的表达量进行验证。根据TAP1基因的序列，使用PrimerPremier5.0软件设计特异性引物，上游引物序列为5'-ATGGTGGTGCTGCTGAAGTT-3'，下游引物序列为5'-CTGCTGCTGCTGCTGATGAT-3'。以β-actin作为内参基因，其上游引物序列为5'-GACCTGACTGACTACCTCAT-3'，下游引物序列为5'-AGCACTGTGTTGGCGTACAG-3'。提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，利用逆转录试剂盒将RNA逆转录为cDNA。以cDNA为模板，在荧光定量PCR仪上进行扩增反应。反应体系为20μL，包括10μL的SYBRGreenMasterMix、0.5μL的上游引物（10μmol/L）、0.5μL的下游引物（10μmol/L）、2μL的cDNA模板以及7μL的ddH₂O。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样本设置3个技术重复，采用2^-ΔΔCt法计算TAP1基因的相对表达量。结果显示，TAP1基因在抗性型断奶仔猪十二指肠组织中的相对表达量为易感型仔猪的3.25倍（P＜0.01），与全基因组测序数据分析得到的差异表达趋势一致，进一步证实了TAP1基因在抗性型仔猪中高表达。接着，通过蛋白质免疫印迹（Westernblot）实验检测TAP1蛋白的表达水平。提取抗性型与易感型断奶仔猪十二指肠组织的总蛋白，采用BCA蛋白定量试剂盒测定蛋白浓度。取30μg的总蛋白进行10%的SDS-PAGE凝胶电泳，电泳结束后将蛋白转移至PVDF膜上。将PVDF膜用5%脱脂奶粉封闭1h，然后加入稀释度为1:1000的兔抗猪TAP1多克隆抗体，4℃孵育过夜。次日，用TBST缓冲液洗涤PVDF膜3次，每次10min，然后加入稀释度为1:5000的HRP标记的羊抗兔IgG二抗，室温孵育1h。再次用TBST缓冲液洗涤PVDF膜3次，每次10min，最后使用化学发光底物显色，在凝胶成像系统下观察并拍照记录结果。以β-actin作为内参蛋白，对TAP1蛋白的表达量进行半定量分析。结果表明，TAP1蛋白在抗性型断奶仔猪十二指肠组织中的表达量明显高于易感型仔猪，灰度值分析显示抗性型仔猪中TAP1蛋白的表达量是易感型仔猪的2.86倍（P＜0.01），这与qRT-PCR检测的基因表达结果相符，从蛋白质水平验证了TAP1基因在抗性型断奶仔猪中的高表达。BPI基因编码的内毒素结合蛋白可以识别并结合细菌内毒素，参与免疫防御反应，在本研究中同样被筛选为候选抗性基因，且在抗性型断奶仔猪中表达上调。针对BPI基因的验证，同样先进行qRT-PCR实验。设计BPI基因的特异性引物，上游引物为5'-TGGACAGCAGCAGCAGTATT-3'，下游引物为5'-CCAGCAGCAGCAGCAGTAAT-3'。按照与TAP1基因qRT-PCR相同的反应体系和条件进行扩增。结果显示，BPI基因在抗性型断奶仔猪十二指肠组织中的相对表达量是易感型仔猪的2.68倍（P＜0.01），验证了测序结果中BPI基因的差异表达情况。在Westernblot验证中，提取总蛋白、定量、电泳、转膜等步骤与TAP1基因验证相同。使用稀释度为1:800的兔抗猪BPI多克隆抗体进行孵育，二抗使用稀释度为1:4000的HRP标记的羊抗兔IgG。结果显示，BPI蛋白在抗性型断奶仔猪中的表达量显著高于易感型仔猪，灰度值分析表明抗性型仔猪中BPI蛋白的表达量是易感型仔猪的2.45倍（P＜0.01），进一步确认了BPI基因在蛋白质水平的差异表达，与基因表达验证结果一致。通过对TAP1基因和BPI基因的qRT-PCR和Westernblot验证，充分证实了这两个基因在抗性型与易感型断奶仔猪中的差异表达，且表达趋势与全基因组测序数据分析结果一致，为后续深入研究这些基因在断奶仔猪抗F18大肠杆菌过程中的功能和作用机制奠定了坚实基础。3.4基因功能与调控网络分析为了深入了解筛选出的差异表达基因在断奶仔猪抗F18大肠杆菌过程中的功能，本研究利用DAVID（DatabaseforAnnotation,VisualizationandIntegratedDiscovery）数据库和GOseq、clusterProfiler等R包，对差异表达基因进行了全面的功能注释和富集分析。在GO功能富集分析中，从生物过程（BiologicalProcess，BP）、细胞组分（CellularComponent，CC）和分子功能（MolecularFunction，MF）三个层面进行深入剖析。在生物过程方面，差异表达基因显著富集在免疫应答（immuneresponse）、细胞对细菌的防御反应（cellulardefenseresponsetobacterium）、炎症反应（inflammatoryresponse）等生物过程。例如，参与免疫应答过程的基因包括TAP1、BPI等，这些基因在免疫系统中发挥着关键作用，TAP1基因编码的膜蛋白能够调节细胞内肽类抗原的转运，参与抗原处理和递呈过程，而BPI基因编码的内毒素结合蛋白可以识别并结合细菌内毒素，参与免疫防御反应。在细胞组分层面，基因主要富集在细胞膜（cellmembrane）、细胞外基质（extracellularmatrix）、溶酶体（lysosome）等细胞结构相关的组分。在分子功能方面，差异表达基因富集在抗菌肽活性（antimicrobialpeptideactivity）、内毒素结合（endotoxinbinding）、肽结合（peptidebinding）等功能。这些富集结果表明，断奶仔猪对F18大肠杆菌的抗性与免疫应答、细胞防御以及抗菌物质的作用等密切相关。KEGG通路富集分析结果显示，差异表达基因显著富集在Toll样受体信号通路（Toll-likereceptorsignalingpathway）、NF-κB信号通路（NF-kappaBsignalingpathway）、细胞因子-细胞因子受体相互作用（cytokine-cytokinereceptorinteraction）等免疫相关信号通路。在Toll样受体信号通路中，Toll样受体能够识别病原体相关分子模式，激活下游的信号传导，最终诱导免疫应答相关基因的表达，增强机体的免疫防御能力。NF-κB信号通路在炎症反应和免疫调节中发挥着核心作用，当细胞受到病原体刺激时，NF-κB被激活，进入细胞核内，调控一系列免疫相关基因的转录。细胞因子-细胞因子受体相互作用通路则参与调节免疫细胞的活化、增殖和分化，以及炎症反应的发生和发展。这些信号通路的富集进一步证实了免疫相关基因在断奶仔猪抗F18大肠杆菌过程中的重要作用。为了进一步探究差异表达基因之间的相互作用关系，构建基因调控网络，本研究利用STRING（SearchToolfortheRetrievalofInteractingGenes/Proteins）数据库预测基因之间的蛋白质-蛋白质相互作用（Protein-ProteinInteraction，PPI）关系，并使用Cytoscape软件进行可视化分析。在STRING数据库中，设置物种为猪（Susscrofa），将筛选出的差异表达基因输入到数据库中，获取基因之间的相互作用信息，包括直接的物理相互作用和间接的功能相互作用。数据库通过整合实验数据、计算预测以及文献挖掘等多种来源的信息，对每对基因的互作关系进行评分，评分越高表示互作关系的可信度越高。将从STRING数据库中获取的基因互作数据导入到Cytoscape软件中，构建基因调控网络。在Cytoscape软件中，每个节点代表一个基因，节点之间的连线表示基因之间的相互作用关系，连线的粗细和颜色表示相互作用的强度和类型。通过对基因调控网络的分析，发现一些关键节点基因，如TAP1、BPI、IL-6（Interleukin-6）、TNF-α（TumorNecrosisFactor-alpha）等。这些关键节点基因在网络中具有较高的度（degree）和中介中心性（betweennesscentrality），表明它们在基因调控网络中起着重要的桥梁和枢纽作用。例如，TAP1基因与多个免疫相关基因存在相互作用，它可以通过调节抗原处理和递呈过程，影响其他免疫基因的表达，从而在断奶仔猪抗F18大肠杆菌的免疫应答中发挥核心作用。BPI基因同样与多个基因相互关联，通过结合细菌内毒素，参与免疫防御反应，调控相关基因的表达。IL-6和TNF-α作为重要的细胞因子，在炎症反应和免疫调节中发挥着关键作用，它们与多个免疫细胞表面的受体相互作用，激活下游的信号传导通路，调控免疫相关基因的表达。对这些关键节点基因的深入研究，将有助于进一步揭示断奶仔猪抗F18大肠杆菌的分子调控机制。四、lncRNA的筛选与鉴定4.1全转录组测序分析为了全面且深入地筛选与断奶仔猪F18大肠杆菌抗性相关的lncRNA，本研究采用IlluminaHiSeq2500测序平台，对10头抗性型断奶仔猪和10头易感型断奶仔猪的十二指肠组织进行全转录组测序。之所以选择十二指肠组织，是因为F18大肠杆菌主要侵袭仔猪的小肠黏膜上皮细胞，十二指肠作为小肠的起始部分，是F18大肠杆菌最先接触和感染的部位，其组织中的基因和lncRNA表达变化可能更直接地反映仔猪对F18大肠杆菌的抗性或易感机制。在测序样本准备阶段，严格按照RNA提取试剂盒（QiagenRNeasyMiniKit）的操作说明，从十二指肠组织中提取总RNA。提取过程中，采取了一系列严格的质量控制措施，如在超净工作台中进行操作，避免RNA酶的污染；使用无RNA酶的耗材和试剂等。提取完成后，利用NanoDrop2000超微量分光光度计精确测定RNA的浓度和纯度，确保RNA的OD260/OD280比值在1.8-2.0之间，以保证RNA的纯度符合要求。同时，采用Agilent2100生物分析仪对RNA的完整性进行检测，要求RNA的RIN值（RNAIntegrityNumber）大于8.0，只有满足这些质量标准的RNA样本才会进入后续的文库构建环节。文库构建是全转录组测序的关键步骤，本研究使用NEBNextUltraRNALibraryPrepKitforIllumina文库构建试剂盒，按照其标准操作流程进行文库构建。在构建过程中，首先通过随机引物逆转录将总RNA逆转录为cDNA，然后对cDNA进行末端修复、加A尾、连接测序接头等一系列操作，以构建出高质量的测序文库。构建好的文库利用Qubit2.0荧光定量仪进行精确的定量，确保文库的浓度准确无误。同时，使用Agilent2100生物分析仪对文库的插入片段大小进行检测，保证文库的质量和均一性。将合格的文库按照一定比例混合后，在IlluminaHiSeq2500测序平台上进行双端150bp测序，以获得高深度、高质量的测序数据。测序完成后，得到的原始数据中可能包含低质量的读段、接头序列以及污染序列等，这些数据会对后续的分析结果产生干扰，因此需要进行严格的数据过滤。利用FastQC软件对原始测序数据进行质量评估，该软件可以从多个维度对数据质量进行分析，包括碱基质量分布、序列长度分布、GC含量分布等。通过FastQC的分析报告，能够直观地了解数据的质量情况，为后续的数据过滤提供依据。使用Trimmomatic软件去除低质量的读段（质量值低于20的碱基占比超过10%的读段）、接头序列以及含N比例超过5%的读段。经过严格的数据过滤后，得到高质量的cleanreads，这些cleanreads将用于后续的数据分析，以确保分析结果的准确性和可靠性。利用Hisat2软件将过滤后的cleanreads与猪的参考基因组（如Susscrofa11.1版本）进行比对。Hisat2软件基于Burrows-Wheeler变换算法，能够快速、准确地将测序读段映射到参考基因组上。在比对过程中，设置合适的参数，如-k10表示每个读段最多尝试10次比对，-t8表示使用8个线程进行比对，以提高比对效率和准确性。比对完成后，得到的比对结果文件（BAM格式）使用Samtools软件进行排序和去重处理，去除可能存在的重复读段，减少后续分析的误差。通过这些步骤，确保每个读段都能准确地映射到参考基因组上，为后续的lncRNA预测和筛选提供可靠的数据基础。4.2lncRNA的验证与表达分析以FUT3-AS1为例，对筛选出的差异表达lncRNA进行验证和表达分析。首先，运用实时荧光定量PCR（qRT-PCR）技术对FUT3-AS1的表达量进行验证。依据FUT3-AS1的序列，使用PrimerPremier5.0软件精心设计特异性引物，上游引物序列设定为5'-CTGCTGCTGCTGCTGATGAA-3'，下游引物序列为5'-GACCTGACTGACTACCTCAT-3'。以U6作为内参基因，其上游引物序列是5'-CTCGCTTCGGCAGCACA-3'，下游引物序列为5'-AACGCTTCACGAATTTGCGT-3'。按照标准操作流程，提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，利用逆转录试剂盒将RNA逆转录为cDNA。以cDNA为模板，在荧光定量PCR仪上展开扩增反应。反应体系精确设置为20μL，其中包含10μL的SYBRGreenMasterMix、0.5μL的上游引物（10μmol/L）、0.5μL的下游引物（10μmol/L）、2μL的cDNA模板以及7μL的ddH₂O。反应条件严格控制为：95℃预变性30s；95℃变性5s，60℃退火30s，共进行40个循环。每个样本均设置3个技术重复，采用2^-ΔΔCt法精准计算FUT3-AS1的相对表达量。结果清晰显示，FUT3-AS1在抗性型断奶仔猪十二指肠组织中的相对表达量显著低于易感型仔猪，仅为易感型仔猪的0.35倍（P＜0.01），这与全转录组测序数据分析得到的差异表达趋势高度一致，有力地证实了FUT3-AS1在抗性型仔猪中的低表达。为了从不同技术角度进一步验证FUT3-AS1的差异表达情况，采用Northernblot实验进行验证。提取抗性型与易感型断奶仔猪十二指肠组织的总RNA，使用甲醛变性琼脂糖凝胶电泳对RNA进行分离。电泳完成后，通过毛细管转移法将RNA转移至尼龙膜上，然后用紫外线交联仪进行交联固定。将固定后的尼龙膜用含有50%甲酰胺、5×SSC、5×Denhardt's试剂、0.1%SDS和100μg/mL鲑鱼精DNA的预杂交液在42℃条件下预杂交2h，以封闭非特异性结合位点。预杂交结束后，将地高辛标记的FUT3-AS1特异性探针加入杂交液中，在42℃条件下杂交过夜。次日，用2×SSC、0.1%SDS溶液在室温下洗涤尼龙膜2次，每次15min，然后用0.1×SSC、0.1%SDS溶液在65℃条件下洗涤2次，每次15min，以去除非特异性杂交的探针。洗涤完成后，利用地高辛检测试剂盒进行显色反应，在暗室中曝光显影，观察并拍照记录结果。以U6作为内参RNA，对FUT3-AS1的表达量进行半定量分析。结果表明，FUT3-AS1在抗性型断奶仔猪十二指肠组织中的表达量明显低于易感型仔猪，条带灰度值分析显示抗性型仔猪中FUT3-AS1的表达量是易感型仔猪的0.32倍（P＜0.01），这与qRT-PCR检测的结果高度相符，从另一个层面验证了FUT3-AS1在抗性型断奶仔猪中的低表达。为了深入了解FUT3-AS1在不同组织中的表达分布情况，采集了抗性型与易感型断奶仔猪的心、肝、脾、肺、肾、胃、十二指肠、空肠、回肠、结肠等10种组织样本，运用qRT-PCR技术检测FUT3-AS1在这些组织中的表达水平。结果显示，FUT3-AS1在所有检测的组织中均有表达，但表达水平存在显著差异。在十二指肠、空肠和回肠等肠道组织中，FUT3-AS1的表达量相对较高，其中在十二指肠中的表达量最高；而在心、肝、脾、肺、肾等组织中，FUT3-AS1的表达量相对较低。进一步比较抗性型与易感型仔猪不同组织中FUT3-AS1的表达差异，发现在肠道组织中，抗性型仔猪FUT3-AS1的表达量显著低于易感型仔猪（P＜0.01），而在其他组织中，虽然也存在一定的表达差异，但差异不具有统计学意义（P＞0.05）。这些结果表明，FUT3-AS1在肠道组织中的表达变化可能与断奶仔猪对F18大肠杆菌的抗性密切相关，为后续深入研究其在肠道免疫中的作用机制提供了重要线索。4.3lncRNA的特征与分类长链非编码RNA（lncRNA）作为一类长度大于200nt的非编码RNA，在生物体内发挥着广泛而重要的调控作用。从特征上看，lncRNA在长度、外显子数目、表达水平、进化保守性等方面都展现出与mRNA不同的特点。在长度方面，lncRNA的长度范围较广，一般在200nt至100,000nt之间，相较于mRNA，其长度通常较短。例如，在人类基因组中，许多lncRNA的长度在1000nt左右，而mRNA的长度则可达到数千nt甚至更长。外显子数目上，lncRNA通常含有较少的外显子，平均外显子数目在2-4个，明显少于mRNA。这种结构特征可能影响其转录和加工过程，进而影响其功能的发挥。在表达水平上，lncRNA的表达量普遍较低，且具有明显的细胞或组织类型特异性。例如，在小鼠的不同组织中，lncRNA的表达谱存在显著差异，某些lncRNA仅在特定组织中高表达，而在其他组织中几乎不表达。这表明lncRNA的表达受到严格的调控，可能与特定组织的生理功能密切相关。进化保守性方面，lncRNA在不同物种间的序列保守性较低，这意味着其功能可能在进化过程中发生了较大的变化，或者其功能不仅仅依赖于序列的保守性，还与其他因素如二级结构等有关。根据lncRNA在基因组上相对于蛋白编码基因的位置，可将其分为多种类型。正义lncRNA与邻近蛋白编码基因转录方向相同，且与编码基因的一个或多个外显子重叠。这种类型的lncRNA可能通过与编码基因共享启动子或转录调控元件，对编码基因的转录产生影响。反义lncRNA则与邻近蛋白编码基因转录方向相反，与相反链上的转录产物部分或完全互补。反义lncRNA可以通过与mRNA形成互补双链，影响mRNA的稳定性、翻译过程或剪切方式，从而调控基因表达。例如，在某些细胞中，反义lncRNA与mRNA结合后，可招募核酸酶降解mRNA，导致基因表达沉默。内含子lncRNA来源于编码基因的内含子，由基因的内含子转录产生。它们可能参与调节基因的可变剪接过程，影响mRNA的成熟和功能。基因间lncRNA由编码基因之间的序列独立转录，它们可以通过远距离作用，调控多个基因的表达。例如，一些基因间lncRNA可以与转录因子结合，形成复合物，然后结合到靶基因的启动子区域，激活或抑制靶基因的转录。双向lncRNA可同时从与邻近蛋白编码基因转录方向相同和相反2个方向发生转录，其功能机制较为复杂，可能同时对多个基因的表达产生影响。增强子lncRNA由蛋白质编码基因的增强子区域产生，它们可以增强基因的转录活性，通过与转录因子、RNA聚合酶等相互作用，促进基因的表达。在本研究中，筛选出的FUT3-AS1属于反义lncRNA，它位于FUT3基因的反义链上，与FUT3基因转录方向相反。FUT3基因已被证实与F18受体的产生密切相关，而FUT3-AS1作为反义lncRNA，可能通过与FUT3mRNA形成互补双链，在多个层面上调控FUT3基因的表达，进而影响断奶仔猪对F18大肠杆菌的抗性。这种反义lncRNA与靶基因之间的相互作用为揭示断奶仔猪抗F18大肠杆菌的分子机制提供了重要线索。五、lncRNA对F18大肠杆菌抗性的调控机制5.1lncRNA与靶基因的相互作用预测在明确FUT3-AS1与断奶仔猪F18大肠杆菌抗性密切相关后，深入探究其与靶基因的相互作用机制至关重要。本研究运用生物信息学预测软件，如RNAplex、IntaRNA等，对FUT3-AS1与FUT3等靶基因之间的相互作用进行预测。这些软件基于核酸序列的互补配对原则以及热力学稳定性等因素，能够精确地预测lncRNA与靶基因之间可能存在的相互作用位点和结合模式。RNAplex软件通过计算RNA分子之间的最小自由能，来预测它们的相互作用位点。在预测FUT3-AS1与FUT3的相互作用时，将FUT3-AS1和FUT3的核酸序列输入到RNAplex软件中，软件会对序列进行分析，计算出不同区域之间的互补配对情况和结合自由能。结果显示，FUT3-AS1的第150-170个核苷酸区域与FUT3基因的5'UTR（非翻译区）的第20-40个核苷酸区域具有较高的互补性，结合自由能较低，表明这两个区域可能存在较强的相互作用。IntaRNA软件则综合考虑了RNA分子的二级结构、序列互补性以及结合能等因素，能够更全面地预测RNA-RNA相互作用。利用IntaRNA软件对FUT3-AS1和FUT3进行分析，发现除了上述区域外，FUT3-AS1的3'端部分序列与FUT3的编码区部分序列也存在一定的互补配对，虽然结合自由能相对较高，但也提示这部分区域可能在特定条件下参与相互作用。通过对多个预测软件结果的综合分析，初步确定了FUT3-AS1与FUT