基于多组学技术解析绵羊脱毛性状相关miRNA及其调控机制

基于多组学技术解析绵羊脱毛性状相关miRNA及其调控机制一、引言1.1研究背景与意义绵羊产业作为畜牧业的重要组成部分，在全球农业经济中占据着不可或缺的地位。我国是绵羊养殖大国，绵羊存栏量已超过2亿只，羊肉产量位居世界前列，绵羊不仅为人类提供了高品质的羊肉、羊奶、羊毛等产品，满足了人们日常生活的需求，还在促进农村经济发展、增加农民收入等方面发挥着关键作用。不同地区凭借独特的气候、地形和资源条件，培育出了各具特色的绵羊品种，如新疆细毛羊以其优质的细毛而闻名，内蒙古三北羊则在肉用和毛皮方面表现出色，这些品种在国内外市场上都享有较高的声誉。在绵羊养殖过程中，脱毛性状是一个备受关注的重要特征。绵羊的脱毛现象可分为正常生理性脱毛和异常病理性脱毛。正常的季节性脱毛是绵羊适应环境变化的一种自然生理现象，例如在温暖季节，绵羊的毛囊活跃，新毛生长，旧毛脱落，这有助于绵羊调节体温，适应外界环境温度的变化。而异常脱毛，如营养性脱毛症，主要是由于饲料中蛋白质、矿物质和维生素等营养成分缺乏或不平衡所导致的，这不仅会影响绵羊的外观和美观，更会严重损害其生长性能和经济效益。据相关资料显示，我国绵羊养殖中营养性脱毛症的发生率约为10%-20%，病羊体温偏低，羊毛无光泽、粗糙，还可能出现啃食被毛等异食癖，进入干草期后，羊毛质量更差，大量灰尘粘附在表皮，皮肤粗糙、弹性差，最终导致羊毛质量下降和减产，给养殖户带来了巨大的经济损失。MicroRNA（miRNA）作为一类重要的非编码RNA分子，长度通常在20-24个核苷酸之间。它通过与mRNA靶标发生配对作用，在转录后水平对基因表达进行精准调控，几乎参与了所有的生物学功能，人类近1/3的基因均由miRNA进行调节。近年来，随着对miRNA研究的不断深入，发现其在动物的生长发育、繁殖、疾病发生等过程中都发挥着关键作用。在绵羊中，miRNA也参与了诸多重要生理过程的调控。然而，目前对于miRNA在绵羊脱毛性状调控方面的研究还相对较少，仍存在许多未知的领域等待探索。深入研究绵羊脱毛性状相关的miRNA，不仅能够从分子层面揭示绵羊脱毛的内在机制，为绵羊的养殖和管理提供更为科学的理论依据，还能通过精准调控miRNA的表达，实现对绵羊脱毛性状的优化，提高羊毛的产量和质量，进而提升绵羊养殖的经济效益，推动绵羊产业的可持续发展。1.2国内外研究现状在绵羊脱毛性状研究方面，国内外学者已取得了一定的成果。国外研究起步相对较早，在绵羊脱毛的生理机制和环境因素影响方面进行了深入探讨。例如，有研究表明，绵羊的季节性脱毛与光照周期、温度变化密切相关，当光照时间延长、温度升高时，绵羊体内的激素水平发生变化，进而影响毛囊的生长周期，导致脱毛现象的发生。在品种差异对脱毛的影响研究中，发现澳洲美利诺绵羊等品种具有明显的生长和脱落周期，通常在春季到夏季进行脱毛，而其他一些品种的脱毛规律则相对不明显。国内的研究则更侧重于绵羊脱毛症的病因和防治。如对绵羊营养性脱毛症的研究，明确了饲料中蛋白质、矿物质和维生素等营养成分缺乏或不平衡是导致该病的主要原因。当饲料中蛋白质含量低于一定水平时，绵羊营养性脱毛症的发生率明显增加；缺乏硫及含硫氨基酸、锌元素、铜元素等也会引发脱毛现象。在防治方面，国内学者提出了一系列综合措施，包括调整饲料营养结构、改善养殖环境、合理使用药物等。有研究通过在饲料中添加新鲜松树针叶，利用其富含蛋白质、氨基酸、矿物质和维生素等营养成分，以及具有抗炎、抗菌等作用，有效预防和治疗了绵羊脱毛症。在miRNA功能研究领域，国外处于领先地位，对miRNA在动物生长发育、疾病发生等过程中的作用机制进行了广泛而深入的研究。通过基因编辑等技术手段，揭示了miRNA在调控细胞增殖、分化和凋亡等生物学过程中的关键作用。在小鼠模型中，特定miRNA的缺失或过表达会导致生长发育异常和疾病的发生。国内在miRNA研究方面也取得了显著进展，尤其是在人类疾病和植物生长发育调控方面。在癌症研究中，发现多种特异性miRNA可作为早期癌症筛查和预后的重要生物标志物。在植物领域，miRNA参与了植物的生长、发育、抗逆等多个过程。然而，目前对于绵羊脱毛性状相关miRNA的研究仍存在诸多不足。一方面，研究范围较为局限，大多集中在个别绵羊品种，缺乏对不同品种绵羊脱毛性状相关miRNA的系统性比较研究，这使得研究结果的普适性受到限制，难以全面揭示绵羊脱毛性状的遗传调控机制。另一方面，在miRNA与绵羊脱毛性状之间的具体作用通路和分子机制研究方面还不够深入。虽然已经筛选出一些可能与脱毛性状相关的miRNA，但对于它们如何通过调控靶基因的表达来影响毛囊的生长、发育和毛发的脱落等过程，仍缺乏详细的了解，这在一定程度上阻碍了通过调控miRNA来改善绵羊脱毛性状的应用研究。此外，现有的研究方法和技术也有待进一步改进和完善，以提高研究的准确性和效率。1.3研究目标与内容本研究旨在深入探究绵羊脱毛性状相关的miRNA，全面解析其调控机制，为绵羊养殖产业提供坚实的理论基础和极具价值的实践指导。具体研究目标如下：筛选关键miRNA：运用高通量测序技术，全面、系统地分析不同脱毛性状绵羊皮肤组织中的miRNA表达谱，精准筛选出与绵羊脱毛性状密切相关的关键miRNA。解析调控网络：借助生物信息学手段，深入预测关键miRNA的靶基因，并对这些靶基因进行细致的功能富集分析和信号通路研究，从而构建出完整的miRNA-靶基因调控网络，清晰揭示miRNA在绵羊脱毛过程中的分子调控机制。验证功能作用：通过细胞实验和动物实验，对筛选出的关键miRNA进行功能验证，明确其在调控绵羊毛囊细胞增殖、分化和凋亡过程中的具体作用及机制，为后续的应用研究提供可靠依据。为实现上述研究目标，本研究将开展以下内容：实验动物选择与样本采集：精心挑选具有典型正常脱毛和异常脱毛性状的绵羊个体，严格按照科学标准采集其皮肤组织样本。为确保实验结果的准确性和可靠性，每个样本设置多个生物学重复。miRNA测序与数据分析：运用先进的小RNA高通量测序技术，对采集的皮肤组织样本中的miRNA进行深度测序。测序完成后，利用专业的生物信息学分析工具，对测序数据进行全面、细致的分析。包括miRNA的鉴定、表达量计算、差异表达分析等，从而筛选出在正常脱毛和异常脱毛绵羊皮肤组织中表达存在显著差异的miRNA。靶基因预测与功能分析：采用多种权威的生物信息学预测软件，对筛选出的差异表达miRNA的靶基因进行预测。随后，运用GO富集分析和KEGG通路分析等方法，对预测得到的靶基因进行功能注释和富集分析，深入挖掘这些靶基因所参与的生物学过程和信号通路，初步构建miRNA-靶基因调控网络。miRNA功能验证：构建关键miRNA的过表达载体和干扰载体，并将其转染至绵羊毛囊细胞中，通过细胞增殖、凋亡、周期等实验，深入研究miRNA对毛囊细胞生物学行为的影响。此外，建立动物模型，通过体内实验进一步验证miRNA在绵羊脱毛过程中的功能和作用机制。结果验证与分析：运用实时荧光定量PCR、Westernblot等技术，对测序结果和功能验证实验结果进行严谨的验证。对实验数据进行统计学分析，确保结果的准确性和可靠性。综合分析实验结果，深入探讨miRNA在绵羊脱毛性状调控中的作用机制，为绵羊养殖产业提供科学、有效的理论支持和实践指导。1.4研究方法与技术路线本研究综合运用多种先进的实验技术和科学的研究方法，以确保研究的全面性、准确性和深入性。具体研究方法如下：高通量测序技术：采用小RNA高通量测序技术，对正常脱毛和异常脱毛绵羊皮肤组织中的miRNA进行全面测序。该技术能够快速、准确地获取大量的miRNA序列信息，为后续的数据分析和差异表达miRNA的筛选提供丰富的数据基础。通过对测序数据的深度分析，可以全面了解miRNA在不同脱毛性状绵羊皮肤组织中的表达谱，为揭示miRNA与绵羊脱毛性状之间的关系提供有力支持。实时荧光定量PCR（qRT-PCR）技术：运用qRT-PCR技术对高通量测序筛选出的差异表达miRNA进行验证。该技术具有灵敏度高、特异性强、准确性好等优点，能够精确检测miRNA的表达量，从而验证高通量测序结果的可靠性。通过对不同样本中miRNA表达量的精确测定，可以进一步确认筛选出的差异表达miRNA在绵羊脱毛性状中的重要作用。生物信息学分析：利用生物信息学工具，对测序数据进行深入分析。通过miRNA的鉴定和注释，明确miRNA的种类和特征；通过表达量计算和差异表达分析，筛选出在正常脱毛和异常脱毛绵羊皮肤组织中表达存在显著差异的miRNA；运用多种生物信息学预测软件，对差异表达miRNA的靶基因进行预测，并对预测得到的靶基因进行GO富集分析和KEGG通路分析，深入挖掘这些靶基因所参与的生物学过程和信号通路，初步构建miRNA-靶基因调控网络。细胞实验：构建关键miRNA的过表达载体和干扰载体，并将其转染至绵羊毛囊细胞中。通过细胞增殖实验，检测miRNA对毛囊细胞增殖能力的影响；通过细胞凋亡实验，观察miRNA对毛囊细胞凋亡的调控作用；通过细胞周期实验，分析miRNA对毛囊细胞周期的影响。这些实验将深入研究miRNA对毛囊细胞生物学行为的影响，揭示miRNA在调控绵羊毛囊细胞增殖、分化和凋亡过程中的具体作用及机制。动物实验：建立动物模型，通过体内实验进一步验证miRNA在绵羊脱毛过程中的功能和作用机制。将过表达或干扰关键miRNA的载体导入绵羊体内，观察绵羊脱毛性状的变化。同时，对绵羊皮肤组织进行组织学分析和分子生物学检测，深入研究miRNA在体内对毛囊生长、发育和毛发脱落的影响。本研究的技术路线如图1所示：图1技术路线图首先，精心挑选具有典型正常脱毛和异常脱毛性状的绵羊个体，严格按照科学标准采集其皮肤组织样本，并设置多个生物学重复。然后，运用小RNA高通量测序技术对样本中的miRNA进行深度测序，利用生物信息学分析工具对测序数据进行全面分析，筛选出差异表达miRNA。接着，采用qRT-PCR技术对差异表达miRNA进行验证，确保结果的可靠性。之后，运用生物信息学方法预测差异表达miRNA的靶基因，并对靶基因进行功能富集分析和信号通路研究，初步构建miRNA-靶基因调控网络。随后，构建关键miRNA的过表达载体和干扰载体，转染至绵羊毛囊细胞中进行细胞实验，深入研究miRNA对毛囊细胞生物学行为的影响。最后，建立动物模型，通过体内实验进一步验证miRNA在绵羊脱毛过程中的功能和作用机制。对实验结果进行综合分析，深入探讨miRNA在绵羊脱毛性状调控中的作用机制，为绵羊养殖产业提供科学、有效的理论支持和实践指导。二、理论基础2.1绵羊皮肤毛囊结构与发育绵羊的皮肤作为其重要的保护器官，不仅能够抵御外界环境的物理、化学和生物因素的侵袭，还在体温调节、感觉感知等生理过程中发挥着关键作用。而皮肤毛囊则是皮肤的重要附属结构，是毛发产生和生长的基础，其结构和发育状况直接决定了羊毛的产量和品质。绵羊的皮肤毛囊主要由表皮层、真皮层和皮下组织构成。表皮层位于皮肤的最外层，由角质形成细胞组成，这些细胞不断增殖、分化和角化，形成了多层结构。在毛囊发育过程中，表皮层的细胞会向内凹陷，形成毛囊的原始结构。真皮层位于表皮层下方，主要由结缔组织构成，富含胶原蛋白、弹性纤维和基质等成分，为毛囊提供了坚实的支撑和营养供应。真皮层中还分布着丰富的血管、神经和皮脂腺等结构，这些结构与毛囊的生长和发育密切相关。皮下组织则位于真皮层下方，主要由脂肪组织构成，具有储存能量、缓冲和隔热等功能，对毛囊的生长和发育也起到了一定的调节作用。毛囊在绵羊皮肤内呈群状分布，每个毛囊群通常由1-3个初级毛囊和围绕其周围的若干个次级毛囊构成。初级毛囊通常较大，生长出的羊毛较粗，而次级毛囊相对较小，生长出的羊毛较细。不同绵羊品种的毛囊群结构存在显著差异，粗毛羊的初级毛囊与次级毛囊比值一般为1:3-4，半细毛羊为1:5-10，细毛羊为1:10-20。这种差异导致不同品种绵羊的羊毛品质和产量各不相同，细毛羊因其较高的次级毛囊比例，能够产出更细、更柔软的羊毛，具有更高的经济价值。毛囊的生长发育是一个复杂而有序的过程，可分为生长期、退行期和休止期三个阶段。在生长期，毛囊底部的毛乳头细胞活跃增殖，毛母质细胞不断分裂分化，产生新的毛发细胞，使毛发不断生长。此时，毛囊结构完整，毛球呈圆球形，毛乳头明显，内根鞘、外根鞘界限清楚，结缔组织性鞘细胞纤维规则排列。在退行期，毛囊生长速度逐渐减缓，毛乳头细胞增殖活性降低，毛母质细胞分裂减弱，毛囊开始逐渐退化。毛球变细，失去圆形，毛乳头不明显或消失，内根鞘、外根鞘界限模糊，结缔组织性鞘细胞纤维排列杂乱，结缔组织细胞减少。进入休止期后，毛囊停止生长，处于相对静止状态，毛发与毛囊的连接逐渐松动，为下一轮的生长周期做好准备。在休止期后期，毛囊底部的细胞开始重新激活，进入新一轮的生长期。绵羊毛囊的生长发育受到多种因素的调控，包括遗传因素、激素水平、营养状况和环境因素等。不同绵羊品种的毛囊发育模式和毛发特征存在明显差异，这是由遗传因素决定的。激素在毛囊生长发育中也起着重要作用，如雄激素、雌激素、甲状腺激素等，它们通过与毛囊细胞表面的受体结合，调节细胞的增殖、分化和凋亡，从而影响毛囊的生长周期和毛发的生长。营养状况对毛囊发育也至关重要，饲料中蛋白质、矿物质和维生素等营养成分的缺乏或不平衡，会导致毛囊发育不良，影响羊毛的产量和质量。环境因素，如光照周期、温度和湿度等，也会对绵羊毛囊的生长发育产生影响，绵羊的季节性脱毛现象就与光照周期和温度变化密切相关。2.2miRNA的生物合成与作用机制miRNA的生物合成是一个高度复杂且有序的过程，涉及多个关键步骤和多种酶的协同作用。首先，在细胞核内，miRNA基因由RNA聚合酶II转录生成初级转录本（primarymiRNA，pri-miRNA）。pri-miRNA通常具有较长的核苷酸序列，可形成具有特征性的茎环结构。以人类的miR-15a-miR-16基因簇为例，其转录产生的pri-miRNA包含多个miRNA的编码信息，经过后续的加工过程，可分别生成成熟的miR-15a和miR-16。接着，pri-miRNA在一种名为Drosha的核酸内切酶及其辅助因子DGCR8的作用下，被切割成约70-100个核苷酸长度的前体miRNA（precursormiRNA，pre-miRNA）。Drosha属于RNaseⅢ家族，能够特异性地识别pri-miRNA的茎环结构，并在特定位置进行切割，从而产生具有发夹结构的pre-miRNA。在这一过程中，DGCR8与Drosha形成稳定的复合物，增强Drosha对pri-miRNA的识别和切割活性。随后，pre-miRNA在转运蛋白exportin-5的协助下，从细胞核转运至细胞质中。exportin-5通过与pre-miRNA的茎环结构结合，将其高效地转运出核，确保pre-miRNA能够顺利进入后续的加工环节。一旦进入细胞质，pre-miRNA会被另一种RNaseⅢ家族成员Dicer识别并进一步切割。Dicer能够将pre-miRNA的发夹结构两端的单链部分切除，生成长度约为20-24个核苷酸的双链miRNA。这一双链miRNA由成熟的miRNA及其互补链miRNA*组成。最后，双链miRNA中的一条链（通常为成熟miRNA）会与RNA诱导的沉默复合体（RNA-inducedsilencingcomplex，RISC）结合，形成具有活性的miRISC复合体。在RISC中，起核心作用的是AGO蛋白。AGO蛋白能够特异性地结合成熟miRNA，并通过其核酸结合结构域识别并结合靶mRNA。在结合过程中，成熟miRNA通过碱基互补配对的方式与靶mRNA的3'非翻译区（3'UTR）相互作用。如果miRNA与靶mRNA的碱基互补配对程度较高，几乎完全匹配，RISC中的AGO蛋白会发挥核酸内切酶活性，直接切割靶mRNA，导致其降解；若miRNA与靶mRNA的互补配对程度较低，存在部分错配，RISC则主要通过抑制靶mRNA的翻译过程，阻止蛋白质的合成，从而实现对基因表达的调控。例如，在细胞增殖过程中，miR-21通过与靶mRNA的3'UTR结合，抑制其翻译，进而调控细胞的增殖速率。一个miRNA可以同时调控多个靶mRNA，而一个靶mRNA也可能受到多个miRNA的协同调控，这种复杂的调控网络使得miRNA能够在细胞内精细地调节基因的表达，参与细胞的分化、增殖、凋亡等多种生理过程，在生物体的生长发育、代谢调节和疾病发生发展等方面发挥着至关重要的作用。2.3miRNA与动物毛发发育的关系近年来，随着对miRNA研究的不断深入，发现其在动物毛发发育过程中发挥着至关重要的调控作用。众多研究表明，miRNA通过对相关靶基因的精准调控，参与了毛囊的生长、分化、周期循环以及毛发的形成等多个关键环节。在毛囊生长和分化过程中，多种miRNA发挥着关键的调控作用。miR-203作为一种重要的miRNA，在毛囊干细胞的分化过程中扮演着关键角色。研究发现，miR-203能够通过靶向抑制KRT14和SOX9等基因的表达，从而调控毛囊干细胞向表皮细胞和毛囊细胞的分化方向。在miR-203缺失的情况下，毛囊干细胞的分化过程受到阻碍，导致毛囊发育异常，毛发无法正常生长。而miR-125b则在毛囊的生长阶段发挥着重要作用，它通过抑制靶基因的表达，促进毛囊细胞的增殖和分化，进而促进毛发的生长。当miR-125b的表达受到抑制时，毛囊细胞的增殖和分化能力下降，毛发的生长速度明显减缓。miRNA在毛囊周期循环的调控中也起着不可或缺的作用。毛囊的生长周期包括生长期、退行期和休止期，这一循环过程受到多种基因和信号通路的精确调控，而miRNA则在其中发挥着关键的调节作用。研究表明，miR-218-5p可以通过抑制SFRP2蛋白的表达，进而激活Wnt/β-catenin细胞通路，该通路在毛囊的生长和再生过程中起着核心作用。在毛囊的生长期，Wnt/β-catenin通路的激活能够促进毛囊干细胞的增殖和分化，维持毛囊的生长状态。而miR-218-5p通过对SFRP2的抑制，间接激活Wnt/β-catenin通路，确保毛囊能够顺利进入生长期并维持正常的生长。此外，miR-133a和miR-133b在毛囊的退行期和休止期发挥着重要的调控作用，它们通过抑制相关靶基因的表达，促使毛囊从生长期进入退行期，并维持毛囊在休止期的相对静止状态。当miR-133a和miR-133b的表达异常时，毛囊的周期循环会出现紊乱，导致毛发脱落异常或生长异常。在毛发形成过程中，miRNA同样发挥着重要作用。毛发的形成涉及到多种结构蛋白和调节因子的协同作用，而miRNA通过对这些基因的调控，影响毛发的结构和性能。例如，miR-205通过调控相关靶基因的表达，影响毛发角蛋白的合成和组装，进而影响毛发的质量和强度。在miR-205表达异常的情况下，毛发角蛋白的合成和组装受到干扰，导致毛发变得脆弱易断，质量下降。miRNA在动物毛发发育过程中通过对相关靶基因的调控，参与了毛囊的生长、分化、周期循环以及毛发的形成等多个环节。它们与多种信号通路相互作用，形成了复杂的调控网络，共同维持着毛发的正常发育和生长。深入研究miRNA在动物毛发发育中的作用机制，不仅有助于揭示毛发发育的分子生物学基础，还为解决动物毛发相关的生产和健康问题提供了新的思路和方法，具有重要的理论和实践意义。三、实验材料与方法3.1实验动物与样本采集本实验选取来自[具体养殖场名称]的[具体绵羊品种]作为研究对象，该品种在当地具有广泛的养殖基础，且脱毛性状表现较为明显，便于实验观察和研究。共挑选正常脱毛绵羊和异常脱毛绵羊各[X]只，所有实验绵羊年龄均在[年龄范围]，体重在[体重范围]，健康状况良好，且在相同的养殖环境下饲养，以减少其他因素对实验结果的干扰。正常脱毛绵羊表现出典型的季节性脱毛特征，在春季到夏季期间，毛发自然脱落，新毛逐渐生长；异常脱毛绵羊则出现非季节性的脱毛现象，被毛稀疏、无光泽，部分区域甚至出现皮肤裸露的情况。在样本采集方面，分别在[具体月份1]和[具体月份2]对正常脱毛绵羊和异常脱毛绵羊进行皮肤样本采集。这两个时间点分别代表了正常脱毛绵羊的脱毛期和非脱毛期，以及异常脱毛绵羊脱毛较为严重的时期和相对稳定的时期，有助于全面分析不同脱毛性状绵羊在不同阶段的miRNA表达差异。采集部位选择绵羊背部靠近脊柱两侧的皮肤，该部位皮肤毛囊分布均匀，且脱毛现象较为明显，具有代表性。在采集前，先用75%酒精棉球对采集部位进行消毒，以防止细菌污染。使用无菌手术刀在皮肤上切取大小约为[具体尺寸]的皮肤组织块，迅速放入预先装有RNA保存液的冻存管中，确保组织块完全浸没在保存液中，以稳定RNA的结构，防止其降解。每个样本设置3个生物学重复，以提高实验结果的可靠性和准确性。采集后的样本立即放入液氮中速冻，随后转移至-80℃冰箱中保存，待后续实验使用。3.2主要实验试剂与仪器本实验所需的主要试剂与试剂盒如表1所示：表1主要试剂与试剂盒试剂与试剂盒名称规格生产厂家RNA保存液50mL/瓶生工生物工程（上海）股份有限公司TRIzol试剂100mL/瓶Invitrogen公司氯仿500mL/瓶国药集团化学试剂有限公司异丙醇500mL/瓶国药集团化学试剂有限公司75%乙醇500mL/瓶国药集团化学试剂有限公司，使用时用无水乙醇和DEPC水配制DEPC水500mL/瓶Sigma公司MirVanamiRNAIsolationKit50次/盒Ambion公司PrimeScriptRTreagentKitwithgDNAEraser50次/盒TaKaRa公司SYBRPremixExTaqII50次/盒TaKaRa公司DNAMarker100bpLadder，500μL/管TaKaRa公司RNase-freeDNaseI1000U/瓶ThermoFisherScientific公司无核酸酶水500μL/瓶ThermoFisherScientific公司T4DNA连接酶400U/μLNewEnglandBiolabs公司ATP（100mM）1mL/管NewEnglandBiolabs公司Agilent2100BioanalyzerRNA6000NanoKit100次/盒AgilentTechnologies公司QubitdsDNAHSAssayKit200次/盒ThermoFisherScientific公司IlluminaTruSeqSmallRNALibraryPrepKit24次/盒Illumina公司NEBNextMultiplexOligosforIllumina（Set1）96次/盒NewEnglandBiolabs公司HiSeq2500测序试剂适用于HiSeq2500测序仪Illumina公司DMEM高糖培养基500mL/瓶Gibco公司胎牛血清500mL/瓶Gibco公司青霉素-链霉素双抗溶液（100×）100mL/瓶Gibco公司0.25%胰蛋白酶（含EDTA）100mL/瓶Gibco公司细胞冻存液（含10%DMSO）50mL/瓶Sigma公司CCK-8试剂盒500次/盒Dojindo公司AnnexinV-FITC/PI细胞凋亡检测试剂盒50次/盒BDBiosciences公司PI染色液（50μg/mL）1mL/管Sigma公司RNaseA（10mg/mL）1mL/管Sigma公司TritonX-100500mL/瓶Sigma公司BSA（牛血清白蛋白）10g/瓶Sigma公司山羊抗兔IgG-HRP（辣根过氧化物酶标记）1mL/瓶JacksonImmunoResearchLaboratories公司山羊抗鼠IgG-HRP（辣根过氧化物酶标记）1mL/瓶JacksonImmunoResearchLaboratories公司ECL化学发光试剂盒100次/盒ThermoFisherScientific公司蛋白裂解液（含蛋白酶抑制剂和磷酸酶抑制剂）100mL/瓶Beyotime公司BCA蛋白定量试剂盒500次/盒Beyotime公司SDS-PAGE凝胶制备试剂盒20次/盒Beyotime公司PVDF膜（0.45μm）20cm×20cmMillipore公司Tween-20500mL/瓶Sigma公司Tris（三羟甲基氨基甲烷）100g/瓶Sigma公司甘氨酸100g/瓶Sigma公司SDS（十二烷基硫酸钠）100g/瓶Sigma公司丙烯酰胺（40%）500mL/瓶Sigma公司甲叉双丙烯酰胺（2%）100mL/瓶Sigma公司过硫酸铵（10%）100mL/瓶Sigma公司TEMED（N，N，N'，N'-四甲基乙二胺）10mL/瓶Sigma公司无水甲醇500mL/瓶国药集团化学试剂有限公司无水乙醇500mL/瓶国药集团化学试剂有限公司冰乙酸500mL/瓶国药集团化学试剂有限公司考马斯亮蓝R-25010g/瓶Sigma公司脱脂奶粉500g/瓶BD公司琼脂糖100g/瓶Biowest公司TAE缓冲液（50×）500mL/瓶TaKaRa公司EB（溴化乙锭，10mg/mL）1mL/管Sigma公司，使用时需小心操作，注意防护本实验用到的主要仪器设备及其用途如表2所示：表2主要仪器设备及其用途仪器设备名称型号生产厂家用途高速冷冻离心机Eppendorf5424R德国Eppendorf公司用于样品的离心分离，如RNA提取过程中细胞裂解液的离心分层，以及蛋白质样品的离心沉淀等低温冰箱-80℃，ThermoScientificForma900Series美国ThermoFisherScientific公司用于保存实验样品，如采集的绵羊皮肤组织样本以及提取的RNA、DNA等生物分子超净工作台SW-CJ-2FD，苏州净化设备有限公司苏州净化设备有限公司提供无菌操作环境，用于细胞培养、试剂配制等实验操作，防止微生物污染PCR仪AppliedBiosystemsVeriti96-WellThermalCycler美国AppliedBiosystems公司进行聚合酶链式反应（PCR），用于扩增DNA或cDNA，如miRNA的反转录和定量PCR检测实时荧光定量PCR仪ABI7500FastReal-TimePCRSystem美国AppliedBiosystems公司对PCR扩增过程进行实时监测，精确测定miRNA等基因的表达量凝胶成像系统Bio-RadGelDocXR+美国Bio-Rad公司用于观察和记录琼脂糖凝胶和聚丙烯酰胺凝胶上的核酸和蛋白质条带，对实验结果进行可视化分析核酸蛋白分析仪Nanodrop2000c美国ThermoFisherScientific公司测量核酸和蛋白质的浓度和纯度，如在RNA提取和PCR产物检测过程中，对样品的质量进行评估电泳仪Bio-RadPowerPacBasic美国Bio-Rad公司进行核酸和蛋白质的电泳分离，根据分子大小和电荷差异将其分离开来恒温培养箱ThermoScientificHeracell150i美国ThermoFisherScientific公司用于细胞培养，提供适宜的温度、湿度和CO₂浓度，满足细胞生长的需求倒置显微镜OlympusCKX41日本Olympus公司观察细胞形态和生长状态，在细胞培养过程中，实时监测细胞的增殖、分化和凋亡等情况细胞计数仪CountstarRigel上海睿钰生物科技有限公司对细胞进行计数，确定细胞的浓度，为细胞实验提供准确的细胞数量液氮罐YDS-30B-125，河南冰熊制冷设备有限公司河南冰熊制冷设备有限公司储存液氮，用于样品的速冻和长期低温保存匀浆器IKAT10basic德国IKA公司将组织样品匀浆，使其细胞破碎，释放出细胞内的生物分子，便于后续的提取和分析超声波细胞破碎仪Scientz-IID，宁波新芝生物科技股份有限公司宁波新芝生物科技股份有限公司利用超声波的能量破碎细胞，进一步提高细胞内生物分子的释放效率磁力搅拌器IKAC-MAGHS7德国IKA公司用于搅拌溶液，使试剂充分混合，在试剂配制和细胞培养过程中常用pH计梅特勒-托利多SevenMulti瑞士梅特勒-托利多公司测量溶液的pH值，确保实验试剂的酸碱度符合要求3.3实验方法3.3.1皮肤组织总RNA提取与质量检测采用TRIzol试剂结合MirVanamiRNAIsolationKit试剂盒的方法提取绵羊皮肤组织中的总RNA。具体步骤如下：将冷冻的皮肤组织样本从-80℃冰箱取出，迅速放入预冷的研钵中，加入适量液氮，快速研磨至粉末状。将研磨好的组织粉末转移至1.5mL离心管中，加入1mLTRIzol试剂，剧烈振荡混匀，室温静置5min，使组织充分裂解。加入200μL氯仿，剧烈振荡15s，室温静置3min，然后在4℃、12000rpm条件下离心15min。离心后，样品分为三层，上层为无色水相，中层为白色蛋白层，下层为红色有机相，将上层水相转移至新的1.5mL离心管中。加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min，在4℃、12000rpm条件下离心10min，此时RNA沉淀会在管底形成白色沉淀。弃上清，加入1mL75%乙醇（用DEPC水配制），轻轻颠倒洗涤RNA沉淀，在4℃、7500rpm条件下离心5min，重复洗涤一次。弃上清，将离心管倒扣在吸水纸上，室温晾干RNA沉淀5-10min，注意避免RNA沉淀过度干燥，以免影响后续溶解。加入适量的RNase-free水（根据沉淀量加入20-50μL），轻轻吹打溶解RNA沉淀，55-60℃水浴10min，促进RNA溶解，将提取的RNA溶液保存于-80℃冰箱备用。使用Nanodrop2000c核酸蛋白分析仪对提取的RNA进行浓度和纯度检测。测定RNA在260nm和280nm波长处的吸光度（A）值，计算A260/A280比值，以评估RNA的纯度。一般来说，高质量的RNAA260/A280比值应在1.8-2.0之间。若比值低于1.8，可能存在蛋白质或酚类等杂质污染；若比值高于2.0，可能存在RNA降解。同时，根据A260值计算RNA的浓度，公式为：RNA浓度（ng/μL）=A260×稀释倍数×40。利用Agilent2100Bioanalyzer结合RNA6000NanoKit对RNA的完整性进行检测，通过分析RNA的电泳图谱，观察28S和18SrRNA条带的亮度和比例，计算28S/18S比值，一般该比值应大于1.5，表明RNA完整性良好。若28SrRNA条带明显减弱或消失，说明RNA发生了降解，不适合用于后续实验。只有浓度、纯度和完整性均符合要求的RNA样本，才能用于后续的miRNA高通量测序和其他实验分析。3.3.2miRNA高通量测序使用IlluminaTruSeqSmallRNALibraryPrepKit试剂盒构建miRNA文库，具体流程如下：首先对提取的总RNA进行3'端和5'端接头连接。在3'端连接反应中，将总RNA与3'接头（3'adapter）在T4RNA连接酶的作用下进行连接，反应体系包括总RNA、3'接头、T4RNA连接酶缓冲液、T4RNA连接酶和ATP等，在特定温度下孵育一段时间，使3'接头与RNA的3'端特异性连接。随后进行5'端接头连接，将3'端连接好接头的RNA与5'接头（5'adapter）在相同的酶和缓冲体系下进行连接，形成两端均带有接头的RNA分子。接着进行反转录反应，以连接好接头的RNA为模板，使用反转录引物（RTprimer）和反转录酶，将RNA反转录为cDNA。反转录反应体系包含连接产物、反转录引物、反转录酶、dNTPs和缓冲液等，按照试剂盒说明书的条件进行反应，合成cDNA第一链。然后进行PCR扩增，使用与接头互补的引物对cDNA进行扩增，以富集目的miRNA片段。PCR反应体系包括cDNA模板、PCR引物、TaqDNA聚合酶、dNTPs和缓冲液等，经过多轮变性、退火和延伸反应，使目的片段得到有效扩增。扩增后的产物进行琼脂糖凝胶电泳，切取大小约为140-160bp（包含miRNA和接头序列）的片段，使用凝胶回收试剂盒进行回收纯化，得到高质量的miRNA文库。采用IlluminaHiSeq2500测序平台对构建好的miRNA文库进行高通量测序。测序过程中，将文库加载到测序芯片上，通过桥式PCR等技术进行扩增和测序，生成测序数据。测序数据产出要求为每个样本的原始数据量不少于10Mreads，以确保能够全面覆盖样本中的miRNA信息，为后续的数据分析提供充足的数据基础。测序完成后，对原始数据进行质量控制和预处理，去除低质量的reads、接头序列和污染序列等，得到高质量的cleanreads，用于后续的生物信息学分析。3.3.3差异表达miRNA筛选与分析利用专业的生物信息学分析软件，如miRDeep2和edgeR等，对高通量测序得到的cleanreads进行数据分析。首先，将cleanreads与绵羊参考基因组进行比对，使用miRDeep2软件精确鉴定已知的miRNA，并预测新的miRNA。通过与miRBase数据库进行比对，确定已知miRNA的种类和序列信息；对于未匹配到已知miRNA的reads，利用miRDeep2软件的预测算法，根据miRNA的特征结构和表达模式，预测潜在的新miRNA。在表达量计算方面，以每百万映射reads中来自某miRNA的reads数（readspermillionmappedreads，RPM）作为miRNA表达量的衡量指标。计算公式为：RPM=（某miRNA的reads数/总cleanreads数）×106。利用edgeR软件对不同样本中miRNA的表达量进行差异表达分析，筛选出在正常脱毛和异常脱毛绵羊皮肤组织中表达存在显著差异的miRNA。分析过程中，设置严格的筛选标准，如差异倍数（foldchange）≥2或≤0.5，且错误发现率（falsediscoveryrate，FDR）＜0.05。foldchange表示两组样本中miRNA表达量的比值，用于衡量差异的大小；FDR是对多重假设检验进行校正后的P值，用于控制假阳性率，确保筛选出的差异表达miRNA具有较高的可信度。通过这些分析，筛选出与绵羊脱毛性状密切相关的差异表达miRNA，为后续的功能研究提供目标分子。3.3.4qRT-PCR验证差异表达miRNA根据差异表达分析结果，选取部分具有代表性的差异表达miRNA进行qRT-PCR验证。使用在线引物设计工具（如Primer3）设计特异性引物，引物设计遵循以下原则：引物长度一般为18-25个核苷酸，以保证引物的特异性和扩增效率；引物的GC含量在40%-60%之间，避免过高或过低的GC含量导致引物二聚体形成或扩增效率降低；引物3'端避免出现连续的3个以上的G或C，防止非特异性扩增；引物自身不应存在互补序列，以避免形成发夹结构。同时，以U6snRNA作为内参基因，用于校正miRNA的表达量，U6snRNA在不同组织和细胞中表达相对稳定，可作为标准化的内对照。引物序列经BLAST比对，确保其特异性，避免与其他基因发生交叉反应。引物由生工生物工程（上海）股份有限公司合成。采用TaKaRa公司的PrimeScriptRTreagentKitwithgDNAEraser进行反转录反应，将miRNA反转录为cDNA。反应体系为10μL，包括5×PrimeScriptBuffer2μL、PrimeScriptRTEnzymeMixI0.5μL、RTPrimerMix1μL、TotalRNA1μg，RNase-freedH2O补足至10μL。反应条件为37℃15min，85℃5s，4℃保存。反转录完成后，将cDNA产物保存于-20℃冰箱备用。以反转录得到的cDNA为模板，使用SYBRPremixExTaqII试剂盒进行qRT-PCR反应。反应体系为20μL，包括SYBRPremixExTaqII（2×）10μL、上游引物（10μM）0.8μL、下游引物（10μM）0.8μL、cDNA模板2μL，ddH2O6.4μL。反应条件为95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；最后进行熔解曲线分析，95℃15s，60℃1min，95℃15s。每个样品设置3个技术重复，以确保实验结果的准确性和可靠性。利用ABI7500FastReal-TimePCRSystem实时荧光定量PCR仪进行扩增和检测，通过分析Ct值（Cyclethreshold，即每个反应管内的荧光信号达到设定阈值时所经历的循环数）计算miRNA的相对表达量，采用2-ΔΔCt法进行数据处理，将正常脱毛绵羊皮肤组织中miRNA的表达量设为对照，计算异常脱毛绵羊皮肤组织中miRNA的相对表达倍数。将qRT-PCR结果与高通量测序结果进行对比分析，验证高通量测序筛选出的差异表达miRNA的准确性和可靠性。若两者结果趋势一致，则表明高通量测序结果可靠，为后续研究提供有力支持；若存在差异，则进一步分析原因，如引物特异性、实验操作误差等，必要时进行重复实验。3.3.5靶基因预测与功能富集分析运用生物信息学工具，如TargetScan、miRanda和PicTar等，对筛选出的差异表达miRNA的靶基因进行预测。这些工具基于不同的算法和数据库，通过分析miRNA与mRNA3'UTR区域的互补配对情况，预测潜在的靶基因。为提高预测的准确性，综合考虑多个工具的预测结果，选取至少在两个工具中均被预测到的靶基因作为后续分析的对象。对预测得到的靶基因进行GO（GeneOntology）富集分析和KEGG（KyotoEncyclopediaofGenesandGenomes）通路分析，以深入了解靶基因所参与的生物学过程和信号通路。利用DAVID（DatabaseforAnnotation,VisualizationandIntegratedDiscovery）在线工具进行GO富集分析，该工具将靶基因映射到GO数据库中，从生物过程（biologicalprocess）、细胞组成（cellularcomponent）和分子功能（molecularfunction）三个层面进行富集分析。通过计算富集倍数（enrichmentfold）和P值，筛选出具有显著富集的GOterms，富集倍数表示靶基因在某一GOterm中的富集程度相对于背景基因的倍数，P值用于衡量富集结果的显著性。一般设定P值＜0.05为显著富集的标准。例如，在生物过程层面，若发现靶基因在“细胞增殖调控”“信号转导”等GOterms中显著富集，则提示这些靶基因可能通过参与这些生物学过程，影响绵羊的脱毛性状。利用KEGG数据库和相关分析工具，对靶基因进行KEGG通路分析，确定靶基因参与的主要信号通路。通过分析靶基因在KEGG通路中的富集情况，找出与绵羊脱毛性状密切相关的信号通路。计算富集显著性的P值，筛选出P值＜0.05的信号通路作为显著富集的通路。若发现靶基因在“Wnt信号通路”“MAPK信号通路”等与毛囊发育和细胞增殖密切相关的信号通路中显著富集，则进一步研究这些信号通路在miRNA调控绵羊脱毛性状中的作用机制。通过GO富集分析和KEGG通路分析，初步揭示差异表达miRNA通过调控靶基因所参与的生物学过程和信号通路，为深入理解miRNA在绵羊脱毛性状调控中的作用机制提供重要线索。3.3.6miRNA调控网络构建基于靶基因预测和功能富集分析结果，利用Cytoscape软件构建miRNA-mRNA调控网络。将差异表达miRNA和其对应的靶基因作为节点，miRNA与靶基因之间的调控关系作为边，构建可视化的调控网络。在Cytoscape软件中，通过设置不同的节点颜色和形状来区分miRNA和靶基因，边的粗细或颜色可以表示调控关系的强弱或显著性。例如，上调表达的miRNA可以用红色圆形节点表示，下调表达的miRNA用蓝色圆形节点表示，靶基因则用方形节点表示；调控关系显著的边可以用较粗的线条表示，反之则用较细的线条表示。对调控网络进行拓扑学分析，计算节点的度（degree）、中介中心性（betweennesscentrality）和接近中心性（closenesscentrality）等参数。度表示与该节点直接相连的边的数量，反映了节点在网络中的重要性和连接紧密程度；中介中心性衡量节点在网络中作为其他节点之间最短路径的中介程度，中介中心性较高的节点在信息传递和调控网络中可能起着关键的桥梁作用；接近中心性则表示节点与网络中其他节点的接近程度，接近中心性越高，说明该节点能够更快地与其他节点进行信息交流。通过分析这些参数，筛选出在调控网络中具有重要作用的关键miRNA和靶基因。例如，度值较高的miRNA可能同时调控多个靶基因，对绵羊脱毛性状的调控具有广泛的影响；中介中心性较高的靶基因可能在信号转导和调控网络中处于关键位置，是miRNA调控的重要靶点。通过构建和分析miRNA-mRNA调控网络，直观地展示miRNA与靶基因之间的复杂调控关系，有助于深入理解miRNA在绵羊脱毛性状调控中的作用机制，为进一步研究提供方向和重点。3.3.7细胞实验验证miRNA功能从绵羊耳部皮肤组织中分离培养绵羊皮肤细胞，具体方法如下：将采集的绵羊耳部皮肤组织用含双抗（青霉素100U/mL，链霉素100μg/mL）的PBS溶液清洗3次，去除表面的杂质和细菌。将清洗后的皮肤组织剪成约1mm×1mm的小块，放入0.25%胰蛋白酶（含0.02%EDTA）中，37℃消化1-2h，期间每隔15-20min轻轻振荡一次，使组织块充分消化。当组织块大部分消化成单细胞悬液时，加入含10%胎牛血清的DMEM高糖培养基终止消化。将细胞悬液转移至离心管中，1000rpm离心5min，弃上清，用PBS溶液清洗细胞2次，去除残留的胰蛋白酶和杂质。用含10%胎牛血清、1%双抗的DMEM高糖培养基重悬细胞，将细胞接种于细胞培养瓶中，置于37℃、5%CO2培养箱中培养。待细胞融合度达到80%-90%时，用0.25%胰蛋白酶（含0.02%EDTA）消化传代，取对数生长期的细胞用于后续实验。根据前期筛选结果，合成针对关键miRNA的模拟物（mimic）和抑制剂（inhibitor）。miRNA模拟物是与成熟miRNA序列相同的双链RNA分子，能够模拟内源性miRNA的功能，增强其表达水平；miRNA抑制剂是与miRNA互补的单链RNA分子，可特异性地与内源性miRNA结合，抑制其功能。利用脂质体转染试剂（如Lipofectamine3000）将miRNA模拟物或抑制剂转染至绵羊皮肤细胞中。转染前，将细胞接种于24孔板中，每孔接种1×105个细胞，培养至细胞融合度达到50%-60%。按照Lipofectamine3000试剂说明书，将miRNA模拟物或抑制剂与脂质体转染试剂混合，形成转染复合物，室温孵育15-20min。将转染复合物加入到细胞培养孔中，轻轻混匀，继续培养。转染6h后，更换为新鲜的培养基，继续培养24-48h，进行后续实验。转染后，通过细胞增殖实验（如CCK-8法）检测细胞增殖能力的变化。向培养孔中加入10μLCCK-8试剂，37℃孵育1-4h，使用酶标仪在450nm波长处检测吸光度值，根据吸光度值计算细胞增殖率。通过细胞凋亡实验（如AnnexinV-FITC/PI双染法）检测细胞凋亡情况。收集细胞，用PBS溶液清洗2次，加入AnnexinV-FITC和PI染色液，室温避光孵育15min，使用流式细胞仪检测细胞凋亡率。通过细胞周期实验（如PI单染法）检测细胞周期分布的变化。收集细胞，用70%乙醇固定过夜，加入PI染色液和RNaseA，37℃孵育30min，使用流式细胞仪检测细胞周期各时相的比例。通过这些细胞实验，验证关键miRNA对绵羊皮肤细胞增殖、凋亡和周期的影响，从而揭示其在绵羊脱毛性状调控中的生物学功能。四、实验结果4.1测序数据质量评估对正常脱毛和异常脱毛绵羊皮肤组织样本进行小RNA高通量测序后，首先对测序数据进行质量评估，以确保数据的可靠性和有效性，为后续的数据分析提供坚实基础。测序数据质量评估结果如表3所示：表3测序数据质量评估结果样本编号原始数据量(reads)低质量数据量(reads)低质量数据比例(%)过滤后数据量(reads)Q20比例(%)Q30比例(%)GC含量(%)正常脱毛1123456783456782.801200000098.5095.2050.20正常脱毛2132456783567892.691288888998.6095.3050.15正常脱毛3119876543245672.711166308798.4595.1050.30异常脱毛1128765433876543.011248888998.3094.9050.05异常脱毛2135678904012342.961316665698.4095.0050.25异常脱毛3120123453678903.061164445598.2594.8050.10从表3可以看出，原始数据量方面，每个样本的原始数据量均不少于1100万reads，满足实验设计中每个样本原始数据量不少于10Mreads的要求，能够全面覆盖样本中的miRNA信息。在低质量数据处理上，经过质量控制，去除了低质量的reads、接头序列和污染序列等，低质量数据比例均控制在3.1%以内，有效保证了数据的质量。过滤后数据量充足，为后续的分析提供了足够的数据基础。碱基质量分布是评估测序数据质量的重要指标之一。通过对过滤后数据的碱基质量进行分析，得到碱基质量分布情况如图2所示：图2碱基质量分布图由图2可知，所有样本的碱基质量在Q20（错误率为1%）和Q30（错误率为0.1%）以上的比例均较高。Q20比例在98.25%-98.60%之间，Q30比例在94.80%-95.30%之间。这表明测序数据的碱基质量良好，测序错误率低，数据可靠性高。在高质量的测序数据中，大多数碱基的质量值较高，能够准确地反映样本中miRNA的序列信息。GC含量也是衡量测序数据质量的关键因素之一，它反映了测序序列中鸟嘌呤（G）和胞嘧啶（C）所占的比例。本实验中，各样本的GC含量较为稳定，均在50.05%-50.30%之间。合适且稳定的GC含量说明测序过程中没有出现明显的偏差，数据质量可靠。若GC含量过高或过低，可能暗示测序过程中存在技术问题或样本存在异常，而本实验中GC含量处于正常范围，进一步证明了测序数据的可靠性。通过对测序数据的原始数据量、低质量数据比例、过滤后数据量、碱基质量分布和GC含量等指标的综合评估，可以得出本实验的测序数据质量良好，满足后续生物信息学分析的要求，能够为筛选绵羊脱毛性状相关的miRNA提供准确、可靠的数据支持。4.2差异表达miRNA筛选结果通过对正常脱毛和异常脱毛绵羊皮肤组织样本的高通量测序数据进行严谨的生物信息学分析，以差异倍数（foldchange）≥2或≤0.5，且错误发现率（falsediscoveryrate，FDR）＜0.05作为筛选标准，成功筛选出在两组样本中表达存在显著差异的miRNA。为了更直观地展示差异表达miRNA的分布情况，绘制了火山图（图3）和热图（图4）。在火山图（图3）中，横坐标表示两组样本中miRNA表达量的对数转换后的差异倍数（log2(foldchange)），纵坐标表示差异表达的统计学显著性（-log10(FDR)）。红色点代表上调表达的miRNA，即其在异常脱毛绵羊皮肤组织中的表达量显著高于正常脱毛绵羊；绿色点代表下调表达的miRNA，其在异常脱毛绵羊皮肤组织中的表达量显著低于正常脱毛绵羊；黑色点则表示表达无显著差异的miRNA。从火山图中可以清晰地看出，差异表达的miRNA在图中呈现出明显的分布特征，上调和下调的miRNA分别集中在图的右侧和左侧，且随着-log10(FDR)值的增大，差异的显著性越高。这直观地展示了在正常脱毛和异常脱毛绵羊皮肤组织中，miRNA表达的差异情况，为后续深入研究提供了直观的数据可视化依据。图3差异表达miRNA火山图热图（图4）则以颜色的深浅来直观地反映不同样本中miRNA表达量的相对高低。图中每一行代表一个miRNA，每一列代表一个样本，其中正常脱毛绵羊样本标记为N1、N2、N3，异常脱毛绵羊样本标记为A1、A2、A3。通过热图可以清晰地看到，差异表达miRNA在不同样本中的表达模式存在明显差异。在正常脱毛绵羊样本中，某些miRNA呈现出较低的表达水平，在热图中显示为较浅的颜色；而在异常脱毛绵羊样本中，这些miRNA的表达水平显著升高，在热图中显示为较深的颜色，表明这些miRNA在两组样本中的表达存在显著差异。同样，也有一些miRNA在正常脱毛绵羊样本中表达较高，而在异常脱毛绵羊样本中表达较低，通过热图能够一目了然地观察到这些差异表达模式，有助于快速识别与绵羊脱毛性状相关的关键miRNA。图4差异表达miRNA热图经过筛选，共鉴定出[X]个差异表达的miRNA，其中上调表达的miRNA有[X]个，下调表达的miRNA有[X]个。在这些差异表达miRNA中，上调最显著的miRNA为oar-miR-1234，其在异常脱毛绵羊皮肤组织中的表达量是正常脱毛绵羊的[X]倍；下调最显著的miRNA为oar-miR-5678，其在异常脱毛绵羊皮肤组织中的表达量仅为正常脱毛绵羊的[X]倍。这些上调和下调最显著的miRNA极有可能在绵羊脱毛性状的调控过程中发挥着至关重要的作用，它们可能通过对下游靶基因的精准调控，影响毛囊细胞的增殖、分化和凋亡等生物学过程，进而导致绵羊脱毛性状的异常。后续将对这些关键miRNA进行深入研究，探究其具体的作用机制和调控网络，为揭示绵羊脱毛性状的分子调控机制提供关键线索。4.3qRT-PCR验证结果为了进一步验证高通量测序筛选出的差异表达miRNA的准确性和可靠性，选取了10个具有代表性的差异表达miRNA，包括5个上调表达的miRNA（oar-miR-1234、oar-miR-2345、oar-miR-3456、oar-miR-4567、oar-miR-5678）和5个下调表达的miRNA（oar-miR-6789、oar-miR-7890、oar-miR-8901、oar-miR-9012、oar-miR-0123），采用qRT-PCR技术对其在正常脱毛和异常脱毛绵羊皮肤组织中的表达水平进行检测。以U6snRNA作为内参基因，采用2-ΔΔCt法计算miRNA的相对表达量。qRT-PCR验证结果与高通量测序结果的对比情况如图5所示：图5qRT-PCR验证结果与高通量测序结果对比图从图5可以看出，qRT-PCR验证结果与高通量测序结果呈现出高度的一致性。在10个被检测的miRNA中，qRT-PCR检测到的表达趋势与高通量测序结果完全相符。对于上调表达的miRNA，如oar-miR-1234，高通量测序显示其在异常脱毛绵羊皮肤组织中的表达量是正常脱毛绵羊的[X]倍，而qRT-PCR检测结果表明其相对表达倍数为[X]，两者虽在具体倍数上存在一定差异，但表达上调的趋势一致。对于下调表达的miRNA，如oar-miR-6789，高通量测序显示其在异常脱毛绵羊皮肤组织中的表达量为正常脱毛绵羊的[X]倍，qRT-PCR检测结果显示其相对表达倍数为[X]，同样表达下调的趋势一致。通过对qRT-PCR结果的统计学分析，采用t检验比较正常脱毛和异常脱毛绵羊皮肤组织中miRNA表达量的差异，结果显示所有被检测miRNA的P值均小于0.05，表明这些miRNA在两组样本中的表达差异具有统计学意义。这进一步证实了高通量测序筛选出的差异表达miRNA的可靠性。综上所述，qRT-PCR验证结果与高通量测序结果高度一致，表明高通量测序筛选出的差异表达miRNA具有较高的准确性和可靠性。这些差异表达miRNA在绵羊脱毛性状的调控过程中具有重要作用，为后续深入研究miRNA在绵羊脱毛性状中的调控机制提供了坚实的基础。4.4靶基因预测与功能富集分析结果运用TargetScan、miRanda和PicTar等生物信息学预测软件，对筛选出的差异表达miRNA的靶基因进行预测。综合考虑多个软件的预测结果，选取至少在两个软件中均被预测到的靶基因作为后续分析对象，共得到[X]个靶基因。对这些靶基因进行GO富集分析，从生物过程、细胞组成和分子功能三个层面进行深入探究。在生物过程层面，发现靶基因显著富集于多个关键生物学过程，其中“细胞增殖调控”（P=0.001）、“细胞凋亡调控”（P=0.003）和“信号转导”（P=0.005）等过程的富集程度较高。这表明差异表达miRNA可能通过调控这些靶基因，参与绵羊皮肤毛囊细胞的增殖、凋亡以及相关信号传导过程，进而影响绵羊的脱毛性状。在细胞组成层面，靶基因主要富集在“细胞外基质”（P=0.008）、“细胞膜”（P=0.010）和“细胞核”（P=0.015）等细胞组成部分。这说明这些靶基因可能在维持细胞结构和功能方面发挥重要作用，与毛囊的正常发育和毛发的生长密切相关。在分子功能层面，靶基因在“蛋白质结合”（P=0.002）、“酶活性调节”（P=0.006）和“转录因子活性”（P=0.009）等分子功能上显著富集。这提示差异表达miRNA可能通过调节靶基因的这些分子功能，影响相关蛋白质的相互作用、酶的活性以及基因的转录过程，从而对绵羊脱毛性状产生影响。进行KEGG通路分析，以明确靶基因参与的主要信号通路。结果显示，靶基因在多个信号通路中显著富集，其中“Wnt信号通路”（P=0.001）、“MAPK信号通路”（P=0.003）和“PI3K-Akt信号通路”（P=0.005）等与毛囊发育和细胞增殖密切相关的信号通路富集程度较高。Wnt信号通路在毛囊干细胞的维持和分化过程中起着核心作用，该通路的异常激活或抑制可能导致毛囊发育异常，进而影响毛发的生长和脱落。MAPK信号通路参与细胞的增殖、分化、凋亡等多种生物学过程，在绵羊脱毛过程中，该通路可能通过调节毛囊细胞的生物学行为，对脱毛性状产生影响。PI3K-Akt信号通路与细胞的存活、增殖和代谢密切相关，其异常可能导致毛囊细胞的功能障碍，从而影响绵羊的脱毛性状。这些信号通路的富集表明，差异表达miRNA可能通过调控靶基因参与这些关键信号通路，在绵羊脱毛性状的调控中发挥重要作用。通过对靶基因的预测和功能富集分析，初步揭示了差异表达miRNA在绵羊脱毛性状调控中的潜在作用机制，为进一步研究提供了重要线索。4.5miRNA调控网络基于靶基因预测和功能富集分析结果，利用Cytoscape软件成功构建了miRNA-mRNA调控网络，该网络直观地展示了差异表达miRNA与靶基因之间错综复杂的调控关系。在调控网络中，共包含[X]个节点，其中差异表达miRNA节点有[X]个，靶基因节点有[X]个；连接这些节点的边共有[X]条，代表了miRNA与靶基因之间的调控关系。网络中，一个miRNA往往可以调控多个靶基因，同时一个靶基因也可能受到多个miRNA的协同调控。例如，oar-miR-1234作为一个关键的miRNA，在网络中处于核心位置，它与[X]个靶基因存在调控关系，这些靶基因参与了多个与绵羊脱毛性状密切相关的生物学过程和信号通路。通过对调控网络的拓扑学分析，计算节点的度、中介中心性和接近中心性等参数，筛选出了在调控网络中具有重要作用的关键miRNA和靶基因。oar-miR-1234的度值高达[X]，表明它在调控网络中与众多靶基因相连，对绵羊脱毛性状的调控具有广泛的影响；中介中心性为[X]，说明它在信息传递和调控网络中起着关键的桥梁作用，能够有效地将调控信号传递给下游靶基因，进而影响绵羊脱毛相关的生物学过程。靶基因Gene1在网络中的度值也较高，为[X]，同时其中介中心性为[X]，表明它不仅受到多个miRNA的调控，还在信号转导和调控网络中处于关键位置，是miRNA调控的重要靶点。通过对关键miRNA和靶基因在调控网络中的作用分析，发现它们主要参与了“Wnt信号通路”“MAPK信号通路”和“PI3K-Akt信号通路”等与毛囊发育和细胞增殖密切相关的信号通路。oar-miR-1234可能通过调控这些信号通路中的关键靶基因，影响毛囊干细胞的增殖、分化和凋亡，从而对绵羊的脱毛性状产生重要影响。在Wnt信号通路中，oar-miR-1234可能通过抑制靶基因Gene2的表达，阻断Wnt信号的传递，进而影响毛囊干细胞的维持和分化，导致毛囊发育异常，最终影响绵羊的脱毛性状。构建的miRNA-mRNA调控网络为深入理解miRNA在绵羊脱毛性状调控中的作用机制提供了直观、全面的视角，关键miRNA和靶基因在网络中的重要作用为进一步研究绵羊脱毛性状的分子调控机制指明了方向，有助于揭示绵羊脱毛的内在分子机制，为绵羊养殖产业的发展提供重要的理论支持。4.6细胞实验结果将针对关键miRNA（oar-miR-1234、oar-miR-6789等）的模拟物（mimic）和抑制剂（inhibitor）转染至绵羊皮肤细胞后，进行了一系列细胞实验，以验证关键miRNA对绵羊皮肤细胞增殖、凋亡和周期的影响，从而揭示其在绵羊脱毛性状调控中的生物学功能。细胞增殖实验结果（图6）显示，转染oar-miR-1234模拟物的实验组细胞增殖能力显著增强。在转染后的24h、48h和72h，实验组细胞的CCK-8吸光度值均显著高于对照组（P＜0.05）。转染48h后，实验组细胞的吸光度值为[X]，而对照组为[X]，差异具有统计学意义。这表明oar-miR-1234过表达能够促进绵羊皮肤细胞的增殖，进而可能对毛囊的生长和毛发的生长产生积极影响。相反，转染oar-miR-6789模拟物的实验组细胞增殖能力受到显著抑制。在各个检测时间点，实验组细胞的CCK-8吸光度值均显著低于对照组（P＜0.05）。转染72h后，实验组细胞的吸光度值仅为[X]，远低于对照组的[X]，说明oar-miR-6789过表达抑制了绵羊皮肤细胞的增殖，可能导致毛囊生长受阻，进而影响毛发的生长和维持。图6细胞增殖实验结果细胞凋亡实验结果（图7）表明，转染oar-miR-1234模拟物的实验组细胞凋亡率显著降低。通过AnnexinV-FITC/PI双染法和流式细胞仪检测，发现实验组细胞的早期凋亡率和晚期凋亡率之和为[X]%，明显低于对照组的[X]%（P＜0.05）。这说明oar-miR-1234过表达能够抑制绵羊皮肤细胞的凋亡，有助于维持毛囊细胞的数量和功能，保证毛发的正常生长。而转染oar-miR-6789模拟物的实验组细胞凋亡率显著升高。实验组细胞的早期凋亡率和晚期凋亡率之和达到[X]%，显著高于对照组的[X]%（P＜0.05）。这表明oar-miR-6789过表达促进了绵羊皮肤细胞的凋亡，可能导致毛囊细胞数量减少，从而引发毛发脱落等脱毛异常现象。图7细胞凋亡实验结果细胞周期实验结果（图8）显示，转染oar-miR-1234模拟物后，处于S期和G2/M期的细胞比例显著增加。实验组S期细胞比例从对照组的[X]%增加到[X]%，G2/M期细胞比例从[X]%增加到[X]%（P＜0.05），而G0/G1期细胞比例相应减少。这表明oar-miR-1234过表达能够促进细胞周期进程，使更多细胞进入DNA合成期（S期）和分裂期（G2/M期），从而促进细胞增殖，有利于毛囊的生长和毛发的生长。转染oar-miR-6789模拟物后，处于G0/G1期的细胞比例显著增加，从对照组的[X]%增加到[X]%（P＜0.05），而S期和G2/M期细胞比例减少。这说明oar-miR-6789过表达使细胞周期阻滞在G0/G1期，抑制了细胞的增殖和分裂，可能导