基于染色质可及性剖析大白猪与民猪胚胎期骨骼肌发育的分子分野

基于染色质可及性剖析大白猪与民猪胚胎期骨骼肌发育的分子分野一、引言1.1研究背景猪作为全球最重要的肉类来源之一，其骨骼肌的生长发育直接决定了猪肉的产量和品质，是动物遗传育种领域的重要研究课题。骨骼肌约占猪体重的40%，不仅是猪肉的主要组成部分，还在猪的运动、代谢等生理过程中发挥着关键作用。在畜牧业生产中，提高猪骨骼肌的生长性能和改善肉质品质，对于满足不断增长的肉类需求、提升养猪业的经济效益具有重要意义。猪胚胎肌肉发育主要集中在两个生长波，胚胎期35-60天为初级纤维形成阶段，胚胎期54-90天则是次级纤维生成阶段。在这一复杂的发育过程中，众多基因参与其中，它们的表达受到精细且严格的调控，从而确保骨骼肌正常发育。然而，目前对于猪胚胎期骨骼肌发育相关基因的调控机制，我们的了解仍十分有限，这在很大程度上制约了猪经济性状分子机理解析以及基因组育种技术的创新发展。染色质可及性在基因表达调控中起着核心作用，是表观遗传学领域的重要研究方向。在细胞核内，染色质并非均匀分布，而是存在着一些特定区域，这些区域的染色质结构相对松散，DNA序列能够与转录因子、染色质修饰酶等调控因子相互作用，从而启动或调节基因的转录过程，这些区域被称为染色质开放区，其具备的这种特性即为染色质可及性。研究表明，染色质开放区约占基因组DNA总序列的2%-3%，而90%以上的开放区域与转录因子的结合密切相关。通过改变染色质的可及性，如利用能量依赖的染色质重塑酶改变核小体的位置和动态，能够选择性地关闭或打开细胞中的基因，进而对细胞的发育和功能产生影响。近年来，随着高通量测序技术的迅猛发展，基于高通量测序技术的染色质可及性分析方法取得了显著进展，为深入研究染色质可及性提供了强大的技术支持。目前，常见的染色质可及性分析技术包括脱氧核糖核酸酶I超敏位点测序（DNase-seq）、微球菌核酸酶测序（MNase-seq）、甲醛辅助分离调控元件测序（FAIRE-seq）以及转座酶可及性测序（ATAC-seq）等。这些技术各自具有独特的原理和优势，能够从不同角度揭示染色质的可及性状态及其在基因调控中的作用机制。例如，DNase-seq主要利用DNase酶1对基因组上具有DNase敏感型的位置进行切割，通过对消化后的片段进行扩增和测序，分析峰值来确定有蛋白质保护的区域，这些区域通常是转录因子以及核小体结合的位置；ATAC-seq则使用高度活跃的Tn5转座酶将Illumina测序适配器插入可达染色质区域，通过特异的标签引物进行PCR扩增和后续测序，从而高效地检测染色质的开放区域。大白猪作为全球广泛养殖的瘦肉型猪品种，具有生长速度快、饲料转化率高、瘦肉率高等优良特性，在现代养猪业中占据着重要地位。民猪是我国宝贵的地方猪种，具有抗寒、耐粗饲、肉质鲜美、繁殖性能好等独特优点，但生长速度相对较慢，瘦肉率较低。这两个猪种在骨骼肌生长发育和肉质性状方面存在明显差异，为研究猪胚胎期骨骼肌发育的分子机制提供了理想的实验材料。通过对大白猪和民猪胚胎期骨骼肌染色质可及性进行研究，有望揭示不同猪种在骨骼肌发育过程中的基因调控差异，挖掘影响骨骼肌生长发育和肉质性状的关键调控元件和基因，为猪的分子育种提供重要的理论依据和潜在的分子标记，对于推动猪遗传育种技术的创新和发展，提高猪肉的产量和品质具有重要的科学意义和应用价值。1.2研究目的和意义本研究旨在通过对大白猪和民猪胚胎期骨骼肌染色质可及性进行深入研究，全面绘制两个猪种胚胎期骨骼肌染色质开放区域图谱，分析染色质可及性在不同发育阶段的动态变化规律，鉴定与骨骼肌发育相关的关键调控元件和基因，进而揭示大白猪和民猪在胚胎期骨骼肌发育过程中的分子机制差异。在猪的遗传育种领域，深入了解骨骼肌发育的分子机制对于培育优良猪种具有至关重要的意义。目前，虽然对猪骨骼肌发育的研究已经取得了一些进展，但对于染色质可及性在猪胚胎期骨骼肌发育过程中的调控作用，以及不同猪种之间在这方面的差异，我们的认识还非常有限。大白猪和民猪在骨骼肌生长发育和肉质性状上的显著差异，为我们研究胚胎期骨骼肌发育的分子机制提供了独特的素材。通过本研究，有望填补这一领域的空白，为猪的分子育种提供新的理论依据和潜在的分子标记。从实践应用角度来看，本研究的成果对于提高猪肉产量和品质具有重要的指导意义。随着人们生活水平的提高，对猪肉的品质要求也越来越高。通过揭示大白猪和民猪胚胎期骨骼肌发育的分子机制差异，可以为猪的遗传改良提供更加精准的方向。例如，针对民猪生长速度慢、瘦肉率低的问题，可以通过分子育种技术，引入大白猪中与骨骼肌生长相关的优势基因或调控元件，同时保留民猪肉质鲜美、繁殖性能好等优良特性，从而培育出既具有良好生长性能，又具备高品质肉质的猪种。这不仅可以满足市场对高品质猪肉的需求，还能提高养猪业的经济效益和竞争力。在基础研究方面，本研究有助于深化我们对骨骼肌发育生物学的理解。骨骼肌发育是一个复杂的生物学过程，涉及多个基因和信号通路的协同作用。染色质可及性作为基因表达调控的重要环节，其在骨骼肌发育过程中的动态变化对于理解基因调控网络具有关键作用。通过对大白猪和民猪胚胎期骨骼肌染色质可及性的研究，可以揭示不同猪种在骨骼肌发育过程中的基因调控模式，为进一步研究其他物种的骨骼肌发育提供参考，推动整个生命科学领域对肌肉发育机制的认识。二、材料与方法2.1实验动物与样本采集本研究选取健康的大白猪和民猪作为实验动物，大白猪购自[具体猪场名称1]，民猪来自[具体猪场名称2]，两个猪场均具备完善的养殖管理体系和良好的养殖环境，确保实验动物的健康和质量。在母猪配种成功后，通过B超技术准确监测怀孕情况，以配种当天为胚胎发育的第0天。分别在胚胎期65天（E65）、75天（E75）和90天（E90）这三个关键发育阶段获取胚胎。选择这三个阶段是因为胚胎期65天处于初级纤维形成阶段后期，此时初级肌纤维已经基本形成，开始向次级纤维形成阶段过渡；胚胎期75天处于次级纤维生成阶段的关键时期，次级肌纤维大量生成；胚胎期90天则是次级纤维生成阶段接近尾声，肌肉组织进一步发育成熟，这三个阶段对于研究猪胚胎期骨骼肌发育具有重要代表性。在获取胚胎时，首先对怀孕母猪进行全身麻醉，采用[具体麻醉药物及剂量]进行静脉注射，确保母猪在手术过程中无痛觉且生命体征稳定。然后通过无菌手术操作，打开母猪的腹腔，小心取出含有胚胎的子宫角。将胚胎从子宫中轻轻分离出来，迅速放入预冷的生理盐水中，清洗掉表面的血液和组织液。在体视显微镜下，仔细辨别并分离出胚胎的骨骼肌组织，包括背最长肌、股四头肌等主要部位，这些部位的骨骼肌在猪的生长发育过程中具有重要功能，且在不同猪种间可能存在明显的发育差异。每个发育阶段，从3头不同母猪所产的胚胎中分别采集骨骼肌样本，每头母猪选取1-2个胚胎，每个胚胎采集适量的骨骼肌组织，以保证样本的生物学重复性和代表性。采集到的骨骼肌样本立即放入液氮中速冻，以迅速停止细胞的代谢活动，防止基因表达和染色质状态发生变化。随后将样本转移至-80℃冰箱中保存，等待后续实验分析。在整个样本采集过程中，严格遵循无菌操作原则，避免样本受到微生物污染，同时确保样本的完整性和质量，以保证实验结果的准确性和可靠性。2.2染色质可及性检测技术2.2.1ATAC-seq原理与实验流程转座酶可及性测序（ATAC-seq）是一种在全基因组水平上研究染色质可及性的重要技术，自2013年由斯坦福大学的HowardChang和WilliamGreenleaf实验室的首席研究员JasonBuenrostro首次发表以来，因其高效、灵敏等优势，在表观遗传学研究领域得到了广泛应用。ATAC-seq技术的核心原理基于DNA转座现象。在自然界中，DNA转座是指DNA转座子在DNA转座酶的作用下，从染色体的一个区域转移到另一个区域。而在ATAC-seq中，利用了高度活跃的Tn5转座酶突变体。这种突变体的转座酶能够将携带已知DNA序列标签的测序接头，以“剪切和粘贴”的方式插入到基因组的开放染色质区域。由于开放染色质区域的核小体分布较少，染色质结构相对松散，使得转座酶能够更容易地接近并切割这些区域的DNA。具体来说，在“tagmentation”过程中，Tn5转座酶会识别并结合到开放染色质区域，然后切割双链DNA，并同时将测序接头标记在切割后的DNA片段两端。这样，原本难以检测的开放染色质区域就被加上了可用于后续分析的测序标签。在本研究中，进行ATAC-seq实验时，首先需要制备高质量的单细胞悬液。从-80℃冰箱中取出保存的大白猪和民猪胚胎期骨骼肌样本，迅速放入预冷的含有蛋白酶抑制剂的裂解缓冲液中，使用组织匀浆器将组织充分匀浆，使细胞分散。随后通过30-70μm的细胞筛网过滤，去除未破碎的组织块和细胞团，获得单细胞悬液，并使用台盼蓝染色法对细胞进行计数和活力检测，确保细胞活力在80%以上。接着，将单细胞悬液在4℃、500-1000g条件下离心5-10分钟，弃去上清液，用预冷的无核酸酶的PBS缓冲液重悬细胞，再次离心洗涤2-3次，以去除杂质和残留的裂解缓冲液。之后，向细胞沉淀中加入适量的含有非离子型去污剂的细胞核裂解缓冲液，轻轻吹打混匀，冰上孵育5-10分钟，使细胞膜裂解，释放出细胞核。将裂解后的细胞悬液通过30μm的细胞筛网过滤，去除未裂解的细胞和细胞碎片，得到纯净的细胞核悬液。将纯化后的细胞核与含有Tn5转座酶和测序接头的转座反应混合物充分混合，在37℃条件下孵育30-60分钟，使Tn5转座酶能够高效地切割开放染色质区域，并将测序接头插入到DNA片段中。转座反应结束后，使用磁珠法对反应产物进行纯化，去除未反应的转座酶、接头和其他杂质。向纯化后的DNA片段中加入特异性的PCR引物，进行12-15个循环的PCR扩增，以富集带有测序接头的DNA片段，构建成高质量的测序文库。使用Agilent2100生物分析仪对文库的片段大小和浓度进行检测，确保文库的质量符合要求。将合格的测序文库采用IlluminaHiSeqXTen等高通量测序平台进行测序，测序模式为双端测序，测序读长为150bp，以获得高质量的测序数据，用于后续的生物信息学分析。2.2.2其他可选用技术对比（如DNase-seq、MNase-seq）脱氧核糖核酸酶I超敏位点测序（DNase-seq）是最早用于研究染色质可及性的技术之一。其原理基于脱氧核糖核酸酶I（DNaseI）能够特异性地切割暴露的、无核小体保护的DNA序列这一特性。在细胞核内，染色质的开放区域由于缺乏核小体的紧密包裹，DNA序列相对暴露，这些区域对DNaseI的消化作用更为敏感，容易被切割成小片段。通过对这些被切割的DNA片段进行高通量测序，并与参考基因组进行比对，就可以确定染色质的开放区域，即DNaseI超敏位点（DHSs）。这些位点通常与基因的调控区域相关，如启动子、增强子等，因为这些区域需要与转录因子等调控蛋白结合，以启动基因的转录过程，所以染色质结构较为开放，对DNaseI的敏感性较高。微球菌核酸酶测序（MNase-seq）则是利用微球菌核酸酶（MNase）对染色质进行消化。MNase能够优先切割核小体之间连接的DNA序列，而对于缠绕在核小体上的DNA具有相对的抗性。通过控制MNase的消化时间和酶量，可以使染色质被消化成不同长度的片段。其中，长度约为147bp的片段通常对应于一个核小体所保护的DNA序列，而更短或更长的片段则反映了核小体的分布情况和染色质的开放程度。对这些消化后的DNA片段进行测序和分析，可以绘制出核小体在基因组上的分布图谱，进而推断染色质的可及性。例如，在染色质开放区域，核小体分布较为稀疏，MNase消化后会产生更多的短片段；而在染色质紧密区域，核小体排列紧密，消化后主要产生147bp左右的片段。与ATAC-seq相比，DNase-seq和MNase-seq虽然在染色质可及性研究中都发挥了重要作用，但也存在一些明显的局限性。从实验操作的复杂性来看，DNase-seq需要精确控制DNaseI的用量和消化时间，因为酶量过多或消化时间过长会导致过度消化，使染色质上的所有DNA都被切割，无法准确区分开放区域和封闭区域；而酶量过少或消化时间过短，则可能无法充分切割开放区域的DNA，导致检测灵敏度降低。MNase-seq同样对酶的用量和消化条件要求严格，且消化过程较为繁琐，需要进行多次优化实验才能获得理想的结果。相比之下，ATAC-seq的实验流程相对简单，只需将细胞核与转座酶混合孵育即可完成染色质的切割和接头标记，大大缩短了实验时间，提高了实验效率。在样本需求量方面，DNase-seq和MNase-seq通常需要大量的起始细胞，一般要求达到10^7数量级，这对于一些珍贵的样本或难以获取大量细胞的实验材料来说，是一个较大的限制。而ATAC-seq具有更高的灵敏度，仅需500-50000个细胞即可完成实验，这使得它能够应用于更多类型的样本，包括临床穿刺样本、少量的胚胎组织等。此外，ATAC-seq在检测染色质可及性时，能够同时获得核小体定位和转录因子结合位点的信息，而DNase-seq和MNase-seq在这方面的能力相对较弱，通常需要结合其他实验技术才能实现对这些信息的全面分析。2.3数据分析方法在本研究中，对ATAC-seq实验获得的原始测序数据进行了一系列严谨且全面的分析，以深入探究大白猪和民猪胚胎期骨骼肌染色质可及性的特征和差异。原始测序数据首先进行质量控制处理，使用FastQC软件对原始测序数据进行评估，以检查数据的质量，包括测序读段的碱基质量分布、GC含量、测序接头污染情况等指标。利用Trimmomatic软件去除测序数据中的低质量碱基（质量值低于20）、接头序列以及含N碱基比例过高的读段，确保后续分析数据的准确性和可靠性。经质量控制后的数据，运用Bowtie2软件将其比对到猪的参考基因组（如Sscrofa11.1版本）上。在比对过程中，设置合适的参数，以保证比对的准确性和效率。例如，调整最大错配数、种子长度等参数，使测序读段能够尽可能准确地定位到基因组上。比对完成后，使用SAMtools软件对结果进行处理，将比对结果从SAM格式转换为BAM格式，并进行排序和索引，以便后续分析。同时，通过计算比对率来评估比对效果，一般要求比对率达到70%以上，若比对率过低，则需检查实验和分析过程，找出问题并进行优化。利用MACS2软件识别染色质开放区域，该软件基于模型的分析算法，能够准确地从比对后的BAM文件中识别出染色质开放区域，即峰（peak）。在分析过程中，设置合适的参数，如设置P值阈值为1e-5，FDR（错误发现率）阈值为0.05，以确保识别出的峰具有较高的可信度。同时，考虑到不同样本间的差异，采用合适的对照组进行比较分析，以排除非特异性信号的干扰。通过MACS2软件的分析，得到每个样本的染色质开放区域的位置信息和相关统计数据，为后续的分析提供基础。在差异分析方面，为了深入了解大白猪和民猪胚胎期骨骼肌染色质可及性的差异，使用DiffBind软件对不同猪种、不同发育阶段的样本进行差异分析。该软件通过对不同样本的染色质开放区域进行比较，能够准确地识别出在不同条件下具有显著差异的开放区域。在分析过程中，设置严格的筛选标准，如倍数变化（foldchange）大于2，FDR小于0.05，以确保筛选出的差异开放区域具有生物学意义。通过差异分析，得到不同猪种、不同发育阶段间染色质开放区域的差异信息，为进一步探究骨骼肌发育的分子机制提供关键线索。对于识别出的染色质开放区域，使用HOMER软件进行注释，以确定这些区域在基因组中的位置和功能。HOMER软件能够将开放区域注释到基因的启动子、增强子、外显子、内含子等不同区域，并提供相关的功能信息。例如，通过注释可以了解开放区域与哪些基因相关，这些基因参与了哪些生物学过程和信号通路，从而为深入研究染色质可及性与基因表达调控的关系提供重要依据。利用MEME-Suite软件进行转录因子结合位点（TFBS）分析，该软件可以在染色质开放区域中搜索潜在的转录因子结合位点，预测可能参与调控的转录因子。通过分析转录因子结合位点的序列特征和保守性，结合已知的转录因子数据库，如JASPAR、TRANSFAC等，确定可能与染色质开放区域相互作用的转录因子，进而深入探讨转录因子在骨骼肌发育过程中的调控作用。三、大白猪与民猪胚胎期骨骼肌染色质可及性特征3.1大白猪胚胎期骨骼肌染色质可及性图谱构建利用ATAC-seq技术，对大白猪胚胎期65天（E65）、75天（E75）和90天（E90）的骨骼肌样本进行分析，成功构建了大白猪胚胎期骨骼肌染色质可及性图谱。在这三个发育阶段，分别鉴定出了数量可观的染色质开放区域，其中E65阶段鉴定出[X1]个染色质开放区域，E75阶段为[X2]个，E90阶段则达到了[X3]个。这些染色质开放区域在基因组上的分布并非随机，而是呈现出一定的特征。从染色体分布来看，染色质开放区域在各条染色体上均有分布，但分布密度存在差异。部分染色体上的开放区域相对集中，如1号、3号和7号染色体，这些染色体上的开放区域可能与骨骼肌发育的关键基因相关，对大白猪胚胎期骨骼肌的生长和分化起着重要作用。通过对开放区域在基因结构上的分布分析发现，约[X4]%的开放区域位于基因间区，这些区域可能包含了大量尚未被揭示的调控元件，对基因的表达调控具有潜在影响；约[X5]%位于内含子区域，虽然内含子不直接编码蛋白质，但它们在基因转录后的加工、调控等过程中发挥着重要作用；而位于启动子区域的开放区域约占[X6]%，启动子区域的染色质开放状态直接关系到基因转录的起始，这些开放区域可能与转录因子的结合密切相关，从而启动基因的表达。对大白猪胚胎期骨骼肌染色质开放区域相关基因进行基因本体（GO）功能富集分析，结果显示，这些基因在多个生物学过程中显著富集。在生物过程方面，主要富集于肌肉器官发育、肌肉系统过程、细胞分化调控等与骨骼肌发育密切相关的过程。其中，肌肉器官发育过程涉及到肌肉组织的形成、发育和成熟，如肌纤维的分化、肌肉组织的构建等，相关基因的表达变化直接影响着骨骼肌的发育进程；肌肉系统过程则涵盖了肌肉的收缩、舒张以及肌肉与神经系统的相互作用等方面，这些过程对于维持骨骼肌的正常功能至关重要；细胞分化调控过程则参与了肌细胞从干细胞向成熟肌纤维的分化过程，调控着细胞命运的决定和分化方向。在分子功能方面，富集于DNA结合、转录因子活性、蛋白质结合等功能类别。DNA结合和转录因子活性使得相关蛋白能够与染色质开放区域的DNA序列特异性结合，从而调节基因的转录；蛋白质结合功能则参与了蛋白质-蛋白质相互作用网络的构建，在信号传导、基因表达调控等过程中发挥重要作用。在细胞组成方面，主要富集于细胞核、染色质等细胞结构相关的类别，表明这些染色质开放区域主要在细胞核内发挥作用，参与染色质的结构重塑和基因表达调控。京都基因与基因组百科全书（KEGG）通路富集分析结果表明，大白猪胚胎期骨骼肌染色质开放区域相关基因显著富集于多条与骨骼肌发育相关的信号通路。其中，MAPK信号通路在细胞增殖、分化、凋亡等过程中发挥着关键作用，在大白猪胚胎期骨骼肌发育过程中，该信号通路的激活可能通过调节相关基因的表达，促进肌细胞的增殖和分化；PI3K-Akt信号通路则与细胞的生长、存活和代谢密切相关，在骨骼肌发育中，该通路的异常激活或抑制可能影响肌细胞的生长和发育；Wnt信号通路在胚胎发育过程中起着重要的调控作用，其在大白猪胚胎期骨骼肌发育中的富集，提示该通路可能参与了肌细胞的命运决定和分化过程。此外，还富集于钙信号通路、肌动蛋白细胞骨架调节等信号通路，这些通路在肌肉收缩、细胞运动和细胞骨架的构建等方面发挥重要作用，对于维持骨骼肌的正常结构和功能具有重要意义。进一步分析大白猪胚胎期骨骼肌染色质开放区域中的调控元件，发现了一些重要的转录因子结合位点。其中，MyoD、Myf5等生肌调节因子的结合位点在染色质开放区域中显著富集。MyoD和Myf5是骨骼肌发育过程中的关键转录因子，它们能够与染色质开放区域的特定DNA序列结合，激活下游与骨骼肌发育相关基因的表达，从而促进肌细胞的分化和肌纤维的形成。此外，还发现了一些与细胞周期调控、代谢调节相关的转录因子结合位点，如E2F家族、PPARγ等。E2F家族转录因子在细胞周期调控中发挥重要作用，其结合位点的富集表明在大白猪胚胎期骨骼肌发育过程中，细胞周期的调控可能与染色质可及性密切相关；PPARγ则参与了脂肪代谢和能量平衡的调节，其在染色质开放区域的结合位点富集，提示脂肪代谢和能量平衡可能在骨骼肌发育过程中也起到一定的作用。这些转录因子结合位点的发现，为深入研究大白猪胚胎期骨骼肌发育的分子调控机制提供了重要线索。3.2民猪胚胎期骨骼肌染色质可及性图谱构建运用ATAC-seq技术对民猪胚胎期65天（E65）、75天（E75）和90天（E90）的骨骼肌样本展开分析，成功构建了民猪胚胎期骨骼肌染色质可及性图谱。在不同发育阶段，民猪胚胎期骨骼肌染色质开放区域的数量和分布特征呈现出动态变化。E65阶段共鉴定出[X7]个染色质开放区域，E75阶段增加至[X8]个，E90阶段达到[X9]个，整体呈现出随着胚胎发育进程，染色质开放区域逐渐增多的趋势。从染色体分布情况来看，染色质开放区域在民猪的各条染色体上均有分布，但分布并不均匀。例如，2号、5号和10号染色体上的染色质开放区域相对较为密集，这些区域可能包含了较多与民猪胚胎期骨骼肌发育密切相关的基因和调控元件。进一步对染色质开放区域在基因结构上的分布进行分析发现，在民猪胚胎期骨骼肌中，约[X10]%的开放区域位于基因间区，这些基因间区的开放区域可能包含了尚未被发现的长链非编码RNA（lncRNA）基因、增强子、沉默子等调控元件，它们通过远程调控等方式影响基因的表达；约[X11]%的开放区域位于内含子区域，内含子不仅参与基因转录后的剪接过程，还可能包含一些顺式作用元件，对基因表达起到调控作用；而位于启动子区域的开放区域约占[X12]%，启动子区域的染色质开放是基因转录起始的重要前提，这些开放区域能够与转录起始复合物等相互作用，启动基因的转录过程。对民猪胚胎期骨骼肌染色质开放区域相关基因进行GO功能富集分析，结果显示，在生物过程方面，主要富集于肌肉结构发育、肌肉收缩调节、细胞增殖与分化调控等生物学过程。肌肉结构发育过程涉及到肌纤维的形成、排列以及肌肉组织的构建等多个环节，对民猪骨骼肌的正常发育至关重要；肌肉收缩调节过程则与肌肉的运动功能密切相关，通过调控相关基因的表达，影响肌肉的收缩和舒张；细胞增殖与分化调控过程参与了肌细胞从干细胞向成熟肌纤维的转化过程，决定了肌细胞的数量和类型。在分子功能方面，富集于蛋白质结合、酶活性调节、转录因子结合等功能类别。蛋白质结合功能参与了蛋白质-蛋白质相互作用网络的形成，在信号传导、基因表达调控等过程中发挥重要作用；酶活性调节功能则对细胞内的各种代谢反应起到调控作用，维持细胞的正常生理功能；转录因子结合功能使得转录因子能够与染色质开放区域的特定DNA序列结合，从而启动或抑制基因的转录。在细胞组成方面，主要富集于肌原纤维、线粒体、细胞核等细胞结构相关的类别。肌原纤维是肌肉收缩的基本单位，其结构和功能的正常维持与染色质可及性密切相关；线粒体是细胞的能量工厂，在肌肉细胞中，线粒体的功能对于提供肌肉运动所需的能量至关重要，染色质开放区域相关基因可能通过调控线粒体的生物发生和功能，影响肌肉的能量代谢；细胞核作为遗传物质的储存场所，染色质开放区域在细胞核内的动态变化直接关系到基因的表达调控。KEGG通路富集分析表明，民猪胚胎期骨骼肌染色质开放区域相关基因显著富集于多条与骨骼肌发育、能量代谢和信号传导相关的信号通路。其中，钙信号通路在肌肉收缩过程中起着关键作用，通过调节细胞内钙离子浓度，影响肌肉的收缩和舒张；AMPK信号通路参与了细胞的能量代谢调节，在民猪胚胎期骨骼肌发育过程中，该通路的激活可能通过调节能量代谢相关基因的表达，维持肌肉细胞的能量平衡；Notch信号通路在细胞分化和发育过程中发挥重要作用，其在民猪胚胎期骨骼肌发育中的富集，提示该通路可能参与了肌细胞的命运决定和分化过程。此外，还富集于PI3K-Akt信号通路、MAPK信号通路等，这些通路在细胞的生长、增殖、存活和凋亡等过程中发挥重要作用，共同调控着民猪胚胎期骨骼肌的发育进程。深入分析民猪胚胎期骨骼肌染色质开放区域中的调控元件，发现了一些关键转录因子的结合位点。其中，Myogenin、Pax7等转录因子的结合位点在染色质开放区域中显著富集。Myogenin是生肌调节因子家族的重要成员，在肌细胞分化过程中发挥着核心作用，它能够与染色质开放区域的特定DNA序列结合，激活下游与肌肉分化相关基因的表达，促进肌细胞向成熟肌纤维的分化；Pax7则主要在卫星细胞中表达，对卫星细胞的维持、增殖和分化具有重要调控作用，其结合位点的富集表明Pax7可能在民猪胚胎期骨骼肌卫星细胞的功能调控中发挥关键作用。此外，还发现了一些与代谢调节、细胞周期调控相关的转录因子结合位点，如PPARα、E2F1等。PPARα参与了脂肪酸代谢和能量平衡的调节，其在染色质开放区域的结合位点富集，提示脂肪酸代谢和能量平衡在民猪胚胎期骨骼肌发育过程中可能起到重要作用；E2F1在细胞周期调控中发挥重要作用，其结合位点的存在表明细胞周期的调控与民猪胚胎期骨骼肌发育密切相关。这些转录因子结合位点的鉴定，为深入研究民猪胚胎期骨骼肌发育的分子调控机制提供了重要线索。3.3两者染色质可及性图谱的共性分析通过对大白猪和民猪胚胎期骨骼肌染色质可及性图谱的深入分析，发现两个猪种在胚胎期骨骼肌染色质可及性方面存在诸多共性。在染色质开放区域的分布特征上，大白猪和民猪呈现出相似的模式。从染色体分布来看，染色质开放区域在各条染色体上均有分布，且在一些染色体上的分布密度表现出相似的偏好。例如，大白猪的1号、3号和7号染色体与民猪的2号、5号和10号染色体，虽然染色体编号不同，但这些染色体上染色质开放区域相对集中，这可能暗示着在不同猪种中，这些染色体区域对胚胎期骨骼肌发育具有重要的调控作用，可能包含了许多保守的调控元件和关键基因。在基因结构上的分布，两个猪种也具有高度一致性。约[X4]%（大白猪）与[X10]%（民猪）的开放区域位于基因间区，[X5]%（大白猪）与[X11]%（民猪）位于内含子区域，[X6]%（大白猪）与[X12]%（民猪）位于启动子区域。基因间区的开放区域可能包含大量尚未被揭示的调控元件，通过长距离调控影响基因表达；内含子区域的开放区域可能参与基因转录后的加工和调控；启动子区域的开放则直接关系到基因转录的起始，这些相似的分布模式表明在猪胚胎期骨骼肌发育过程中，染色质开放区域在基因结构上的功能分布具有一定的保守性。对大白猪和民猪胚胎期骨骼肌染色质开放区域相关基因进行基因本体（GO）功能富集分析，结果显示在多个生物学过程、分子功能和细胞组成方面存在显著的共性。在生物过程方面，均显著富集于肌肉器官发育、肌肉系统过程、细胞分化调控等与骨骼肌发育密切相关的过程。肌肉器官发育过程涵盖了肌纤维的形成、肌肉组织的构建和成熟等环节，是骨骼肌发育的基础；肌肉系统过程涉及肌肉的收缩、舒张以及与神经系统的相互作用，对于维持骨骼肌的正常功能至关重要；细胞分化调控过程则决定了肌细胞从干细胞向成熟肌纤维的分化方向和进程，两个猪种在这些关键生物过程上的富集，表明它们在胚胎期骨骼肌发育的基本生物学过程调控机制上具有相似性。在分子功能方面，都富集于DNA结合、转录因子活性、蛋白质结合等功能类别。DNA结合和转录因子活性使得转录因子能够与染色质开放区域的DNA序列特异性结合，启动或调节基因的转录；蛋白质结合功能参与了蛋白质-蛋白质相互作用网络的构建，在信号传导、基因表达调控等过程中发挥重要作用，这些共性反映了两个猪种在基因表达调控的分子机制上具有一定的保守性。在细胞组成方面，主要富集于细胞核、染色质、肌原纤维、线粒体等细胞结构相关的类别。细胞核是遗传物质的储存和基因转录的场所，染色质的动态变化直接影响基因表达；肌原纤维是肌肉收缩的基本单位，其结构和功能的正常维持与染色质可及性密切相关；线粒体是细胞的能量工厂，为肌肉细胞的活动提供能量，两个猪种在这些细胞组成相关类别的富集，表明它们在胚胎期骨骼肌细胞的基本结构和功能维持方面具有相似的调控机制。京都基因与基因组百科全书（KEGG）通路富集分析结果表明，大白猪和民猪胚胎期骨骼肌染色质开放区域相关基因在多条信号通路上存在共性。其中，MAPK信号通路、PI3K-Akt信号通路、Wnt信号通路、钙信号通路等在两个猪种中均显著富集。MAPK信号通路在细胞增殖、分化、凋亡等过程中发挥着关键作用，通过激活一系列下游激酶，调节基因表达，影响细胞的生物学行为；PI3K-Akt信号通路与细胞的生长、存活和代谢密切相关，在胚胎期骨骼肌发育中，该通路的激活可以促进肌细胞的增殖和存活，调节肌肉的生长和发育；Wnt信号通路在胚胎发育过程中起着重要的调控作用，参与细胞命运决定、分化和组织器官的形成，在猪胚胎期骨骼肌发育中，Wnt信号通路可能通过调控相关基因的表达，影响肌细胞的分化和肌纤维的形成；钙信号通路在肌肉收缩过程中起着关键作用，通过调节细胞内钙离子浓度，控制肌肉的收缩和舒张，两个猪种在这些信号通路上的共同富集，表明这些信号通路在猪胚胎期骨骼肌发育过程中具有保守的调控作用，可能是维持骨骼肌正常发育的关键信号传导途径。四、大白猪与民猪胚胎期骨骼肌染色质可及性差异分析4.1差异染色质开放区域的识别在对大白猪和民猪胚胎期骨骼肌染色质可及性图谱构建及共性分析的基础上，进一步深入挖掘两个猪种之间染色质可及性的差异。利用DiffBind软件对不同猪种、不同发育阶段的样本进行差异分析，以倍数变化（foldchange）大于2且错误发现率（FDR）小于0.05作为严格的筛选标准，旨在识别出在大白猪和民猪胚胎期骨骼肌中具有显著差异的染色质开放区域，为后续探究骨骼肌发育的分子机制差异提供关键线索。在胚胎期65天（E65），共识别出[X13]个差异染色质开放区域，其中大白猪特异性开放区域有[X14]个，民猪特异性开放区域为[X15]个。这些差异开放区域在基因组上的分布呈现出一定的特征，部分区域集中在特定的染色体上，如大白猪特异性开放区域在5号染色体上较为集中，而民猪特异性开放区域在8号染色体上相对较多。在基因结构上，大白猪特异性开放区域主要分布在基因间区（约[X16]%）和内含子区域（约[X17]%），民猪特异性开放区域则在启动子区域（约[X18]%）和基因间区（约[X19]%）分布较为显著。这表明在胚胎发育的早期阶段，两个猪种在染色质开放区域的分布上就存在明显差异，这些差异可能与基因的表达调控以及早期骨骼肌发育的不同进程密切相关。随着胚胎发育至75天（E75），差异染色质开放区域的数量增加至[X20]个，其中大白猪特异性开放区域为[X21]个，民猪特异性开放区域有[X22]个。此时，差异开放区域在染色体上的分布发生了一些变化，大白猪特异性开放区域在3号和7号染色体上的分布密度增加，而民猪特异性开放区域在6号和10号染色体上更为集中。在基因结构上，大白猪特异性开放区域在启动子区域的比例有所上升（约[X23]%），而民猪特异性开放区域在基因间区的比例进一步提高（约[X24]%）。这一时期差异染色质开放区域的变化，可能反映了在胚胎期骨骼肌发育的关键时期，两个猪种在基因表达调控和细胞分化等过程中的差异进一步加剧。到胚胎期90天（E90），差异染色质开放区域的数量达到[X25]个，其中大白猪特异性开放区域有[X26]个，民猪特异性开放区域为[X27]个。在染色体分布上，大白猪特异性开放区域在1号和11号染色体上表现出较高的密度，民猪特异性开放区域则在2号和9号染色体上更为集中。在基因结构上，大白猪特异性开放区域在基因间区和启动子区域的分布相对均衡（分别约为[X28]%和[X29]%），民猪特异性开放区域则主要集中在基因间区（约[X30]%）。这表明在胚胎发育的后期阶段，两个猪种的染色质可及性差异依然存在，且可能对骨骼肌的成熟和功能完善产生重要影响。通过对不同发育阶段差异染色质开放区域的识别和分析，发现随着胚胎发育进程，大白猪和民猪之间差异染色质开放区域的数量和分布呈现出动态变化。这些差异可能与两个猪种在骨骼肌生长速度、肉质性状等方面的差异密切相关，为进一步研究猪胚胎期骨骼肌发育的分子机制提供了重要的研究方向。4.2差异区域相关基因功能注释与富集分析对识别出的差异染色质开放区域相关基因进行全面而深入的功能注释与富集分析，是揭示大白猪和民猪胚胎期骨骼肌发育分子机制差异的关键步骤。本研究运用多种生物信息学工具和数据库，从生物学过程、分子功能和细胞组成等多个维度，对这些基因的功能进行了系统解析。通过DAVID、Metascape等功能注释工具，将差异区域相关基因映射到基因本体（GO）数据库和京都基因与基因组百科全书（KEGG）数据库中，以确定它们在不同生物学过程、分子功能和细胞组成中的作用。在生物学过程方面，发现差异基因显著富集于肌肉发育调控、能量代谢调节、细胞增殖与凋亡调控等关键过程。其中，在肌肉发育调控过程中，大白猪中与肌纤维生长和分化相关的基因如MyoD1、Myf5等，其染色质开放区域在特定发育阶段表现出较高的可及性，这些基因通过激活下游靶基因的表达，促进肌细胞的增殖和分化，从而推动大白猪骨骼肌的快速生长；而民猪中与肌肉结构稳定性和收缩功能相关的基因如TNNT2、MYH3等，其染色质开放区域的可及性变化更为显著，这些基因可能通过调节肌肉的结构和功能，影响民猪骨骼肌的发育进程。在能量代谢调节过程中，大白猪中与脂肪酸氧化和糖代谢相关的基因如CPT1B、PFKFB3等，其染色质开放区域的可及性较高，这表明大白猪在胚胎期骨骼肌发育过程中，可能更倾向于利用脂肪酸和糖类作为能量来源，以满足其快速生长的能量需求；而民猪中与线粒体呼吸链和氧化磷酸化相关的基因如MT-CO1、MT-ATP6等，其染色质开放区域的可及性变化更为明显，这可能反映了民猪在能量代谢方面的独特机制，通过增强线粒体的功能，提高能量利用效率，以适应其生长环境和生长速度相对较慢的特点。在细胞增殖与凋亡调控过程中，大白猪中与细胞周期调控和凋亡抑制相关的基因如CCND1、BCL2等，其染色质开放区域的可及性较高，这有助于促进大白猪肌细胞的增殖和存活，加快骨骼肌的生长；而民猪中与细胞凋亡诱导和细胞周期阻滞相关的基因如BAX、CDKN1A等，其染色质开放区域的可及性变化更为显著，这可能导致民猪肌细胞的增殖速度相对较慢，但同时也有助于维持细胞的稳定性和正常发育。在分子功能方面，差异基因主要富集于转录因子活性、蛋白激酶活性、离子结合等功能类别。其中，转录因子活性在基因表达调控中起着核心作用，大白猪和民猪中均存在一些与骨骼肌发育相关的关键转录因子，如MyoD、Myogenin等，它们通过与染色质开放区域的特定DNA序列结合，激活或抑制下游靶基因的表达，从而调控骨骼肌的发育进程。蛋白激酶活性则参与了细胞内信号传导通路的激活和调节，通过磷酸化修饰下游蛋白，影响细胞的生物学行为。离子结合功能对于维持细胞内离子平衡和信号传导具有重要意义，如钙离子结合蛋白在肌肉收缩过程中起着关键作用，通过调节细胞内钙离子浓度，控制肌肉的收缩和舒张。在细胞组成方面，差异基因主要富集于细胞核、线粒体、肌原纤维等细胞结构相关的类别。细胞核是遗传物质的储存和基因转录的场所，染色质开放区域在细胞核内的动态变化直接关系到基因的表达调控；线粒体是细胞的能量工厂，在肌肉细胞中，线粒体的功能对于提供肌肉运动所需的能量至关重要，染色质开放区域相关基因可能通过调控线粒体的生物发生和功能，影响肌肉的能量代谢；肌原纤维是肌肉收缩的基本单位，其结构和功能的正常维持与染色质可及性密切相关，染色质开放区域相关基因可能通过调节肌原纤维的组成和结构，影响肌肉的收缩功能。KEGG通路富集分析结果显示，差异基因显著富集于多条与骨骼肌发育、能量代谢和信号传导相关的信号通路。其中，MAPK信号通路在细胞增殖、分化、凋亡等过程中发挥着关键作用，大白猪中该信号通路的激活可能通过调节相关基因的表达，促进肌细胞的增殖和分化，而民猪中该信号通路的调控可能更为精细，通过不同的信号分子和下游靶基因，影响肌细胞的命运决定和分化方向；PI3K-Akt信号通路与细胞的生长、存活和代谢密切相关，在大白猪胚胎期骨骼肌发育过程中，该通路的激活可能促进肌细胞的生长和增殖，而民猪中该通路可能更多地参与了细胞的存活和代谢调节；Wnt信号通路在胚胎发育过程中起着重要的调控作用，大白猪和民猪中该通路的差异可能导致肌细胞的分化和肌纤维的形成过程存在差异，进而影响骨骼肌的生长和发育；此外，还富集于钙信号通路、AMPK信号通路等，这些通路在肌肉收缩、能量代谢调节等方面发挥重要作用，其在大白猪和民猪中的差异调控可能与两个猪种骨骼肌的生长速度、肉质性状等差异密切相关。4.3关键差异基因及调控网络分析在差异染色质开放区域相关基因功能注释与富集分析的基础上，进一步筛选出在大白猪和民猪胚胎期骨骼肌发育过程中起关键作用的差异基因，并构建它们之间的调控网络，深入分析这些基因在调控网络中的作用和相互关系，有助于全面揭示两个猪种胚胎期骨骼肌发育分子机制的差异。通过严格的筛选标准，结合基因功能注释和富集分析结果，筛选出了一批与骨骼肌发育、能量代谢和信号传导密切相关的关键差异基因。在骨骼肌发育相关基因中，MyoD1、Myf5、Myogenin等生肌调节因子在大白猪和民猪中的表达模式和染色质可及性存在显著差异。MyoD1在大白猪胚胎期骨骼肌中的染色质开放区域可及性较高，其表达水平也相对较高，这可能促进了大白猪肌细胞的增殖和分化，使其骨骼肌生长速度更快；而在民猪中，Myogenin的染色质开放区域可及性变化更为显著，其表达可能在民猪骨骼肌的分化和成熟过程中发挥更为关键的作用，有助于维持民猪骨骼肌的正常结构和功能。在能量代谢相关基因方面，CPT1B、PFKFB3、MT-CO1等基因在两个猪种中的差异明显。CPT1B和PFKFB3在大白猪胚胎期骨骼肌中染色质开放区域的可及性较高，它们分别参与脂肪酸氧化和糖代谢过程，这表明大白猪在胚胎期骨骼肌发育过程中可能更依赖脂肪酸和糖类作为能量来源，以满足其快速生长的能量需求；而MT-CO1在民猪中的染色质开放区域可及性变化更为显著，该基因参与线粒体呼吸链和氧化磷酸化过程，可能反映了民猪在能量代谢方面的独特机制，通过增强线粒体的功能，提高能量利用效率，以适应其生长环境和生长速度相对较慢的特点。利用Cytoscape等软件，基于基因之间的相互作用关系、信号传导通路以及转录因子调控等信息，构建了关键差异基因的调控网络。在这个调控网络中，MyoD1、Myf5等生肌调节因子处于核心地位，它们通过与其他基因之间的相互作用，形成了复杂的调控网络。MyoD1可以激活下游与肌细胞增殖和分化相关的基因，如Myf5、Myogenin等，同时也受到其他转录因子和信号通路的调控。例如，Wnt信号通路中的关键分子β-catenin可以与MyoD1的启动子区域结合，促进其表达，从而进一步调控骨骼肌的发育进程。能量代谢相关基因在调控网络中也与其他基因存在密切的相互作用。CPT1B和PFKFB3可以通过调节能量代谢产物的水平，影响细胞内的信号传导通路，进而调控与骨骼肌发育相关基因的表达。MT-CO1等线粒体相关基因则通过维持线粒体的正常功能，为整个调控网络提供能量支持，确保基因表达和信号传导过程的正常进行。通过对调控网络的拓扑结构分析，发现一些基因具有较高的节点度和中介中心性，这些基因在调控网络中起着关键的桥梁和枢纽作用。例如，MyoD1不仅与多个生肌调节因子和骨骼肌发育相关基因相互作用，还参与了多个信号通路的调控，其在调控网络中的重要性不言而喻。这些关键基因的鉴定，为进一步深入研究大白猪和民猪胚胎期骨骼肌发育的分子机制提供了重要的靶点。五、差异染色质可及性对骨骼肌发育的影响5.1与肌肉发育相关信号通路的关联差异染色质可及性区域相关基因在肌肉发育相关信号通路中发挥着关键作用，这些信号通路的异常调控可能导致骨骼肌发育异常，进而影响猪的生长性能和肉质性状。通过对大白猪和民猪胚胎期骨骼肌差异染色质可及性区域相关基因的深入分析，发现它们与多条重要的肌肉发育信号通路密切相关，其中Wnt信号通路和Notch信号通路在两个猪种的骨骼肌发育过程中表现出显著的差异调控，对骨骼肌的生长、分化和成熟产生重要影响。Wnt信号通路在胚胎发育和组织稳态维持中起着至关重要的作用，在骨骼肌发育过程中，该通路参与了肌细胞的增殖、分化和肌纤维的形成等关键过程。在大白猪胚胎期骨骼肌中，差异染色质可及性区域相关基因在Wnt信号通路中的富集程度较高，其中一些关键基因如Wnt3a、β-catenin等的染色质可及性发生显著变化。研究表明，Wnt3a作为Wnt信号通路的重要配体，能够与细胞膜上的Frizzled受体和LRP5/6共受体结合，激活下游的β-catenin信号。当Wnt3a与受体结合后，会抑制β-catenin的降解，使其在细胞质中积累并进入细胞核，与TCF/LEF转录因子家族结合，从而激活下游靶基因的表达，如c-Myc、CyclinD1等，这些基因参与细胞周期调控和细胞增殖，促进大白猪肌细胞的增殖和生长，使其骨骼肌生长速度更快。而在民猪胚胎期骨骼肌中，虽然Wnt信号通路也有一定程度的富集，但相关基因的染色质可及性变化模式与大白猪存在差异。例如，民猪中Wnt5a基因的染色质可及性相对较高，Wnt5a属于非经典Wnt信号通路的配体，它与Frizzled受体结合后，通过激活下游的PKC、JNK等信号分子，调节细胞骨架的重组和细胞迁移，在民猪骨骼肌发育过程中，可能更侧重于调节肌细胞的迁移和组织形态的构建，对维持骨骼肌的正常结构和功能具有重要作用，而非像大白猪那样主要促进细胞增殖。这种差异可能与两个猪种的生长特性和肉质性状差异密切相关，大白猪生长速度快，需要更强的细胞增殖能力来支持骨骼肌的快速生长；而民猪生长速度相对较慢，但肉质鲜美，其骨骼肌发育可能更注重结构和功能的稳定性，因此Wnt信号通路的调控方式也有所不同。Notch信号通路同样在骨骼肌发育过程中发挥着关键作用，它参与了肌细胞的命运决定、分化和增殖等过程。在大白猪和民猪胚胎期骨骼肌中，差异染色质可及性区域相关基因在Notch信号通路中的富集情况和调控机制也存在差异。在大白猪中，Notch1基因的染色质可及性较高，Notch1是Notch信号通路的重要受体，当Notch1与配体Delta-like或Jagged结合后，会发生蛋白水解切割，释放出Notch胞内结构域（NICD），NICD进入细胞核与CSL转录因子结合，激活下游靶基因如Hes1、Hey1等的表达，这些基因抑制肌细胞的分化，促进肌细胞的增殖和自我更新，从而有助于大白猪骨骼肌的快速生长。在民猪胚胎期骨骼肌中，Notch信号通路相关基因的染色质可及性变化模式与大白猪有所不同。例如，民猪中Delta-like1基因的染色质可及性相对较高，Delta-like1作为Notch信号通路的配体，与Notch受体结合后，激活的信号通路可能更侧重于调节肌细胞的分化和成熟。研究表明，Delta-like1与Notch受体结合后，通过激活下游的信号分子，促进肌细胞向成熟肌纤维的分化，抑制肌细胞的增殖，这可能使得民猪骨骼肌在发育过程中，肌纤维的分化和成熟更为充分，有助于形成更稳定的肌肉结构，进而影响民猪的肉质性状，使其肉质更加鲜美。这种Notch信号通路调控机制的差异，进一步体现了大白猪和民猪在胚胎期骨骼肌发育过程中的分子机制差异，为深入理解两个猪种的生长性能和肉质性状差异提供了重要线索。5.2对肌肉细胞增殖、分化和凋亡的影响差异染色质可及性对肌肉细胞的增殖、分化和凋亡过程具有显著影响，这种影响在大白猪和民猪胚胎期骨骼肌发育过程中表现出明显的差异，进一步揭示了两个猪种在骨骼肌生长和发育机制上的不同。在肌肉细胞增殖方面，大白猪胚胎期骨骼肌中，一些与细胞增殖相关的基因，如PCNA（增殖细胞核抗原）、CCND1（细胞周期蛋白D1）等，其染色质开放区域的可及性较高。PCNA是DNA聚合酶的辅助蛋白，在DNA复制过程中发挥关键作用，其染色质可及性的增加，使得相关转录因子更容易结合到其调控区域，促进基因转录，从而提高细胞内PCNA的表达水平，加速DNA复制进程，促进肌肉细胞的增殖。CCND1则是细胞周期G1期向S期转变的关键调节因子，其染色质可及性的升高，有助于激活细胞周期相关信号通路，推动肌肉细胞进入增殖周期，增加细胞数量，这与大白猪生长速度快的特性相符合。相比之下，民猪胚胎期骨骼肌中与细胞增殖相关基因的染色质可及性变化相对较小。这可能导致民猪肌肉细胞的增殖速度相对较慢，但也使得细胞增殖过程更为稳定和有序。研究表明，民猪中一些抑制细胞增殖的基因，如p21、p53等，其染色质开放区域的可及性相对较高。p21是一种细胞周期蛋白依赖性激酶抑制剂，能够与细胞周期蛋白-细胞周期蛋白依赖性激酶复合物结合，抑制其活性，从而阻止细胞周期的进展，抑制细胞增殖。p53作为一种重要的肿瘤抑制因子，在细胞受到应激或损伤时，能够激活下游一系列基因的表达，诱导细胞周期阻滞、凋亡或衰老，以维持细胞基因组的稳定性。民猪中这些抑制基因染色质可及性的增加，可能使得民猪肌肉细胞在增殖过程中更加注重质量控制，减少异常增殖的发生，有助于维持骨骼肌的正常发育和结构稳定性。在肌肉细胞分化方面，大白猪和民猪也表现出明显的差异。大白猪胚胎期骨骼肌中，生肌调节因子MyoD1和Myf5等基因的染色质可及性在早期较高，这有利于促进肌细胞向成肌细胞的分化。MyoD1和Myf5能够与染色质开放区域的特定DNA序列结合，招募相关的转录辅助因子，形成转录起始复合物，激活下游一系列与肌肉分化相关基因的表达，如Myogenin、MEF2等，这些基因进一步促进成肌细胞的分化和肌纤维的形成，使得大白猪在胚胎期能够快速形成大量的肌纤维，为其后期的快速生长奠定基础。而民猪胚胎期骨骼肌中，Myogenin基因的染色质可及性在分化后期相对较高。Myogenin是肌肉分化的关键调节因子，在成肌细胞融合和肌纤维成熟过程中发挥着重要作用。民猪中Myogenin染色质可及性在后期的增加，表明民猪在胚胎期骨骼肌发育过程中，更注重肌纤维的成熟和功能完善，通过增强Myogenin的表达，促进成肌细胞的融合和肌纤维的进一步分化，形成更为稳定和成熟的肌肉结构，这与民猪肉质鲜美、肌肉品质优良的特性密切相关。对于肌肉细胞凋亡，大白猪和民猪胚胎期骨骼肌中相关基因的染色质可及性也存在差异。在大白猪中，一些抗凋亡基因，如BCL2（B细胞淋巴瘤-2）等，其染色质开放区域的可及性较高。BCL2是一种重要的抗凋亡蛋白，能够抑制细胞色素C从线粒体释放，阻止caspase级联反应的激活，从而抑制细胞凋亡，促进肌肉细胞的存活，这有助于维持大白猪骨骼肌细胞在快速生长过程中的数量和功能。民猪胚胎期骨骼肌中，促凋亡基因，如BAX（Bcl-2相关X蛋白）等，其染色质可及性相对较高。BAX是一种促凋亡蛋白，能够在线粒体外膜上形成孔道，促进细胞色素C的释放，激活caspase级联反应，诱导细胞凋亡。民猪中BAX染色质可及性的增加，可能使得民猪在胚胎期骨骼肌发育过程中，对异常或受损细胞进行及时清除，维持骨骼肌细胞群体的质量和稳定性，这对于民猪骨骼肌的正常发育和肉质品质的形成具有重要意义。5.3对猪肉品质相关性状的潜在影响差异染色质可及性对猪肉品质相关性状具有潜在的重要影响，这种影响通过多种途径得以体现，与肌肉纤维类型、肌内脂肪含量以及肉质风味等关键品质性状密切相关。在肌肉纤维类型方面，不同类型的肌肉纤维在代谢特性、收缩速度和能量消耗等方面存在显著差异，进而对猪肉的嫩度、多汁性和风味产生重要影响。研究表明，I型肌纤维（慢肌纤维）含量较高的猪肉，通常具有更好的嫩度和多汁性，口感更为鲜美。大白猪和民猪胚胎期骨骼肌染色质可及性的差异，可能导致其肌肉纤维类型组成的不同。例如，大白猪中与快肌纤维发育相关的基因，如MYH1、MYH2等，其染色质开放区域的可及性相对较高，这可能促进了快肌纤维的生成，使得大白猪骨骼肌在生长过程中，快肌纤维的比例相对较高。快肌纤维具有较高的收缩速度和爆发力，但在肉质方面，其嫩度和多汁性相对较差。相比之下，民猪中与慢肌纤维发育相关的基因，如MYH7等，其染色质开放区域的可及性变化更为显著，这可能有利于慢肌纤维的发育，使得民猪骨骼肌中慢肌纤维的比例相对较高，从而在肉质上表现出更好的嫩度和多汁性。肌内脂肪含量是影响猪肉品质的另一个重要因素，它与猪肉的风味、嫩度和多汁性密切相关。适量的肌内脂肪能够赋予猪肉独特的风味和良好的口感。大白猪和民猪在胚胎期骨骼肌染色质可及性上的差异，可能对肌内脂肪的沉积和代谢产生不同的影响。研究发现，民猪中一些与脂肪代谢相关的基因，如FASN（脂肪酸合酶）、LPL（脂蛋白脂肪酶）等，其染色质开放区域的可及性相对较高，这可能促进了脂肪的合成和沉积，使得民猪的肌内脂肪含量相对较高。FASN是脂肪酸合成的关键酶，其染色质可及性的增加，有利于提高FASN的表达水平，从而促进脂肪酸的合成，增加肌内脂肪的含量。LPL则在脂蛋白代谢中发挥重要作用，能够促进脂肪的摄取和沉积。而大白猪中与脂肪分解相关的基因，如ATGL（脂肪甘油三酯脂肪酶）、HSL（激素敏感性脂肪酶）等，其染色质开放区域的可及性相对较高，这可能导致脂肪分解代谢增强，肌内脂肪含量相对较低。ATGL和HSL是脂肪分解的关键酶，它们的染色质可及性增加，使得这些酶的表达水平升高，促进脂肪的分解，减少肌内脂肪的沉积。肉质风味是猪肉品质的重要体现，它受到多种因素的影响，包括挥发性风味物质的组成和含量等。大白猪和民猪胚胎期骨骼肌染色质可及性的差异，可能通过影响相关基因的表达，进而影响挥发性风味物质的合成和代谢，最终导致肉质风味的不同。例如，一些与挥发性风味物质合成相关的基因，如醛脱氢酶、醇脱氢酶等，其染色质开放区域的可及性在两个猪种中存在差异。在民猪中，这些基因的染色质可及性相对较高，可能促进了挥发性风味物质的合成，使得民猪的肉质风味更加浓郁。醛脱氢酶能够催化醛类物质的氧化，生成相应的酸，这些酸是挥发性风味物质的重要组成部分；醇脱氢酶则参与醇类物质的代谢，与挥发性风味物质的形成密切相关。而在大白猪中，这些基因的染色质可及性相对较低，可能导致挥发性风味物质的合成减少，肉质风味相对较淡。六、结论与展望6.1研究主要结论本研究通过对大白猪和民猪胚胎期骨骼肌染色质可及性进行深入研究，成功构建了两个猪种胚胎期不同发育阶段的染色质可及性图谱，揭示了其染色质可及性的特征、差异以及这些差异对骨骼肌发育和猪肉品质的影响，取得了一系列重要研究成果。在染色质可及性图谱构建方面，利用ATAC-seq技术，分别在大白猪和民猪胚胎期65天（E65）、75天（E75）和90天（E90）的骨骼肌样本中，鉴定出大量染色质开放区域。这些开放区域在基因组上呈现出非随机分布，且在不同发育阶段存在动态变化。在染色体分布上，部分染色体上的开放区域相对集中，如大白猪的1号、3号和7号染色体，民猪的2号、5号和10号染色体；在基因结构上，开放区域主要分布在基因间区、内含子区域和启动子区域，且不同猪种在各区域的分布比例存在一定差异。对染色质开放区域相关基因的功能分析发现，大白猪和民猪在多个生物学过程、分子功能和细胞组成方面存在共性。在生物过程中，均显著富集于肌肉器官发育、肌肉系统过程、细胞分化调控等与骨骼肌发育密切相关的过程；在分子功能上，都富集于DNA结合、转录因子活性、蛋白质结合等功能类别；在细胞组成方面，主要富集于细胞核、染色质、肌原纤维、线粒体等细胞结构相关的类别。KEGG通路富集分析显示，两个猪种在多条信号通路上存在共性，如MAPK信号通路、PI3K-Akt信号通路、Wnt信号通路、钙信号通路等，这些通路在猪胚胎期骨骼肌发育过程中具有保守的调控作用。然而，大白猪和民猪在胚胎期骨骼肌染色质可及性方面也存在显著差异。通过严格的差异分析，在不同发育阶段识别出大量差异染色质开放区域，这些差异区域的数量和分布随着胚胎发育进程呈现动态变化。对差异区域相关基因的功能注释与富集分析表明，差异基因在肌肉发育调控、能量代谢调节、细胞增殖与凋亡调控等关键生物学过程中发挥重要作用。例如，在肌肉发育调控方面，大白猪中与肌纤维生长和分化相关的基因如MyoD1、Myf5等，其染色质开放区域的可及性较高，促进了肌细胞的增殖和分化，使其骨骼肌生长速度更快；而民猪中与肌肉结构稳定性和收缩功能相关的基因如TNNT2、MYH3等，其染色质开放区域的可及性变化更为显著，有助于维持民猪骨骼肌的正常结构和功能。在能量代谢调节方面，大白猪中与脂肪酸氧化和糖代谢相关的基因如CPT1B、PFKFB3等，其染色质开放区域的可及性较高，更依赖脂肪酸和糖类作为能量来源，以满足快速生长的能量需求；民猪中与线粒体呼吸链和氧化磷酸化相关的基因如MT-CO1、MT-ATP6等，其染色质开放区域的可及性变化更为明显，通过增强线粒体的功能，提高能量利用效率，以适应其生长环境和生长速度相对较慢的特点。在细胞增殖与凋亡调控方面，大白猪中与细胞周期调控和凋亡抑制相关的基因如CCND1、BCL2等，其染色质开放区域的可及性较高，促进肌细胞的增殖和存活；民猪中与细胞凋亡诱导和细胞周期阻滞相关的基因如BAX、CDKN1A等，其染色质开放区域的可及性变化更为显著，导致肌细胞的增殖速度相对较慢，但有助于维持细胞的稳定性和正常发育。进一步筛选出关键差异基因，并构建其调控网络。在骨骼肌发育相关基因中，MyoD1、Myf5、Myogenin等生肌调节因子在两个猪种中的表达模式和染色质可及性存在显著差异；在能量代谢相关基因方面，CPT1B、PFKFB3、MT-CO1等基因在两个猪种中的差异明显。这些关键差异基因通过相互作用，形成复杂的调控网络，在猪胚胎期骨骼肌发育过程中发挥着核心作用。差异染色质可及性对骨骼肌发育具有重要影响，与肌肉发育相关信号通路密切关联。Wnt信号通路和Notch信号通路在大白猪和民猪胚胎期骨骼肌发育过程中表现出显著的差异调控。在Wnt信号通路中，大白猪中Wnt3a、β-catenin等基因的染色质可及性变化促进肌细胞的增殖，而民猪中Wnt5a基因的染色质可及性相对较高，更侧重于调节肌细胞的迁移和组织形态的构建；在Notch信号通路中，大白猪中Notch1基因的染色质可及性较高，抑制肌细胞的分化，促进肌细胞的增殖和自我更新，民猪中Delta-like1基因的染色质可及性相对较高，促进肌细胞向成熟肌纤维的分化，抑制肌细胞的增殖。差异染色质可及性还对肌肉细胞的增殖、分化和凋亡过程产生显著影响。在肌肉细胞增殖方面，大白猪中与细胞增殖相关基因的染色质可及性较高，促进肌肉细胞的增殖，民猪中与细胞增殖相关基因的染色质可及性变化相对较小，且一些抑制细胞增殖的基因染色质可及性相对较高，使得民猪肌肉细胞的增殖速度相对较慢，但更注重质量控制；在肌肉细胞分化方面，大白猪中生肌调节因子MyoD1和Myf5等基因的染色质可及性在早期较高，促进肌细胞向成肌细胞的分化，民猪中Myogenin基因的染色质可及性在分化后期相对较高，更注重肌纤维的成熟和功能完善；在肌肉细胞凋亡方面，大白猪中抗凋亡基因的染色质可及性较高，促进肌肉细胞的存活，民猪中促凋亡基因的染色质可及性相对较高，对异常或受损细胞进行及时清除，维持骨骼肌细胞群体的质量和稳定性。此外，差异染色质可及性对猪肉品质相关性状具有潜在影响。在肌肉纤维类型方面，大白猪中与快肌纤维发育相关的基因染色质可及性相对较高，快肌纤维比例相对较高，民猪中与慢肌纤维发育相关的基因染色质可及性变化更为显著，慢肌纤维比例相对较高，导致两者在肉质嫩度和多汁性上存在差异；在肌内脂肪含量方面，民猪中与脂肪代谢相关的基因染色质可及性相对较高，促进脂肪的合成和沉积，肌内脂肪含量相对较高，大白猪中与脂肪分解相关的基因染色质可及性相对较高，导致肌内脂肪含量相对较低；在肉质风味方面，一些与挥发性风味物质合成相关的基因染色质可及性在两个猪种中存在差异，民猪中这些基因的染色质可及性相对较高，可能促进了挥发性风味物质的合成，使得民猪的肉质风味更加浓郁。6.2研究的创新点与不足本研究在猪胚胎期骨骼肌染色质可及性研究领域具有多个创新点。在研究方法上，创新性地运用ATAC-seq技术对大白猪和民猪胚胎期不同发育阶段的骨骼肌样本进行分析。相较于传统的染色质可及性分析技术，如DNase-seq和MNase-seq，ATAC-seq具有操作简便、样本需求量少、灵敏度高等显著优势。其独特的原理基于高度活跃的Tn5转座酶，能够高效地将测序接头插入到开放染色质区域，大大缩短了实验周期，且仅需少量细胞即可完成实验，这使得对珍贵的胚胎期骨骼肌样本的研究成为可能，为该领域的研究提供了更为高效和精准的技术手段。在研究内容方面，本研究首次全面系