基于GWAS的深县猪繁殖性状遗传解析与育种应用研究

基于GWAS的深县猪繁殖性状遗传解析与育种应用研究一、引言1.1研究背景与意义深县猪作为华北黑猪祖源，是京津冀地区唯一古老地方猪品种，至今已延续8000多年，在我国地方猪种中占据重要地位。它是河北省地方黑猪代表品种，主要分布于河北省中南部，形成于独特农业生态环境，经劳动人民长期选育而成。该品种具有肉质细嫩、抗病力强、繁殖力强、抗寒耐热、耐粗饲等优良特性，尤其是其肌内脂肪含量高达12%以上，肉质软嫩不柴，烹调后肉香浓郁，是不可多得的优质猪种资源。然而，自上世纪90年代起，由于生长速度较慢、瘦肉率较低等劣势，在大批国外猪种引进改良的浪潮中，深县猪养殖数量急剧下降，一度濒临灭绝。尽管近年来在各方努力下，通过搜寻民间个体、建立保种场等措施，深县猪得以保存延续，如河北正农牧业有限公司建立保种场，保存10个家系，种猪达500头，但整体种群规模依然较小，遗传多样性保护和品种改良工作迫在眉睫。繁殖性状是影响猪养殖经济效益和种群发展的关键因素，对于深县猪而言同样如此。母猪的总产仔数、死胎数、断奶仔猪数、窝产仔猪体重及性别比例等繁殖性状指标，不仅决定了深县猪种群的扩繁速度，还与养殖成本、猪肉产量和质量密切相关。深入研究深县猪繁殖性状的遗传机制，对提升其繁殖性能、扩大种群规模、保护这一珍稀地方猪种具有重要意义。若能通过遗传改良提高深县猪的产仔数，可直接增加养殖收益，同时加快种群恢复速度；了解影响仔猪成活率和生长发育的遗传因素，能针对性优化饲养管理，降低养殖成本，提高猪肉品质。全基因组关联分析（GWAS）技术作为现代遗传学研究的重要手段，近年来在畜禽遗传育种领域得到广泛应用。它通过检测特定物种中不同个体间的全部或大部分基因，分析基因与性状之间的关联，能够快速定位与目标性状相关的基因组区域和候选基因。相较于传统育种方法，GWAS技术具有无需构建专门的遗传群体、可同时分析多个性状、检测效率高等优势，为深入解析复杂性状的遗传基础提供了有力工具。将GWAS技术应用于深县猪繁殖性状研究，有望挖掘出影响其繁殖性能的关键基因和分子标记，为深县猪的遗传改良和分子育种提供理论依据和技术支持。通过GWAS分析，可能发现调控深县猪产仔数、胚胎成活率等繁殖性状的关键基因，基于这些基因开发分子标记，可实现早期选种选育，加速遗传改良进程，培育出生长周期短、产仔率高、抗病力强、肉质优良的深县猪新品种，满足市场对优质猪肉的需求，推动深县猪产业的可持续发展。1.2深县猪概述深县猪历史源远流长，最早可追溯至新石器时代，在河北徐水南庄头遗址出土的距今10000年以前的猪骨骸，以及邯郸武安县磁山文化遗址中7000年前的猪骨骸，都表明此地养猪历史悠久。东晋张华所著《博物论》记载“生青兖徐淮者耳大，生燕冀者皮厚”，这成为深县猪早期特征的有力佐证。深县猪以原产于深县（今深州市）得名，主要分布于河北省中南部的黑龙港流域，包括桃城区、深州市、武强县、安平县以及晋州市、辛集市、宁晋县等地。此地地势平坦、气候温和，为深县猪的生长和繁殖提供了良好的自然条件；同时，当地丰富的农副产品和青粗饲料资源，也满足了深县猪耐粗饲的特性需求。在品种特征方面，深县猪体型较大，成年公猪体重150-200公斤，体长130厘米，体高75厘米，胸围110厘米；成年母猪体重100-150公斤，体长120厘米，体高65厘米，胸围115厘米。其外貌特征显著，被毛黑色，皮厚，毛粗而密，冬季密生棕红色绒毛，这使其能适应不同的气候条件。耳朵大且下垂超过鼻端，嘴筒短粗，末端齐平，鼻孔开阔，方便采食和呼吸；额部有较深的菱形皱纹，呈现出“五花头”的独特头型，头大颈短，体躯窄而深，前躯发育良好，后躯略差；腹部膨满，乳房发达，有效乳头7-9对，分布均匀，为仔猪提供充足的母乳；四肢较短，粗壮有力，便于在各种地形活动；尾端有长毛一簇，尾根较粗且附着较高，甩动灵活自如；鬃毛发达，长7-8厘米，背腰部、臀部及四肢部皮肤有皱褶。深县猪按头型可分为“黄瓜嘴”和“五花头”（又可细分为“大五花头”“小五花头”），还有介于两者之间的“稀毛白皮”。“黄瓜嘴”嘴筒长，末端尖，鼻直平少毛，额上皱纹较浅，颈细长，前驱较窄，背腰凹而长，腹微垂，臀部倾斜，后肢较弱，尾粗大；“五花头”头粗糙，毛粗密，额宽皱多折深，体粗大，颈粗大，全身皮厚多褶皱，嘴宽而相对偏短，两耳垂至嘴角以下，前躯发育较好，胸部宽，背部略有凹陷，后躯不丰满，尾巴粗大下垂，四肢粗大，系部较软；“稀毛白皮”体型丰满匀称，不塌腰，背毛相对稀短，嘴筒上翘。深县猪的生产性能也较为出色。在繁殖性能上，性成熟早，初产母猪平均窝产仔10头左右，经产母猪窝产仔12-14头，多者超20头，以第五、六胎产仔最多，仔猪初生重0.5-1.0公斤。母猪护仔力强，躺卧授乳时姿势和动作缓慢、谨慎，能有效避免压死仔猪，授乳过程中四肢极度伸展，让全部仔猪都能自由吸乳，分娩后半月内，除饲养人员外拒绝生人走近。在生长性能方面，深县猪生长速度较慢，一般断奶后10个月龄，体重可达90公斤，当地有“紧七、慢八、十到家”的说法，即断奶后，喂的饲料好7个月即可出栏，饲料差8个月出栏，最长的10个月出栏。在肉质性能上，其肉质细嫩、鲜美可口，肌内脂肪含量高达12%以上，肉香浓郁，呈大理石雪花纹路。而且深县猪耐粗饲，具有发达的肠管，小肠平均长18.96米，大肠平均长4.39米，能充分消化利用青粗饲料，在喂给较少粗饲料和青粗饲料的情况下，也能发育生长，在喂给较多的精料时，生长迅速。此外，深县猪抗病力强、适应性广，能在粗放的饲养条件下生存繁衍。深县猪作为我国地方猪种资源的重要组成部分，是华北黑猪祖源，具有不可替代的遗传多样性价值。它不仅承载着我国悠久的养猪历史和文化，还为现代猪育种提供了丰富的基因资源，对培育适应不同市场需求的新品种具有重要意义。在当前猪肉市场对肉质品质要求不断提高的背景下，深县猪优良的肉质特性使其具有广阔的开发利用前景，有望成为满足消费者对高品质猪肉需求的优质猪种。1.3猪繁殖性状研究进展1.3.1繁殖性状指标及影响因素猪的繁殖性状是衡量其养殖效益和种群发展的关键指标，涵盖多个方面。产仔数是最为重要的指标之一，包括总产仔数和活产仔数。总产仔数反映了母猪一次分娩所产出的仔猪总数，而活产仔数则关注其中存活的仔猪数量。在实际养殖中，产仔数直接决定了仔猪的供应量，对养殖规模的扩大起着关键作用。例如，在规模化养猪场，母猪产仔数的多少直接影响到后续育肥猪的存栏量，进而影响猪肉的产量和经济效益。初生重和初生窝重也是重要的繁殖性状指标。初生重指仔猪出生时个体的重量，初生窝重则是一窝仔猪出生时的总重量。这些指标与仔猪的健康状况和存活率密切相关。一般来说，初生重较大的仔猪在出生后具有更强的生命力和适应能力，能够更好地抵抗疾病和外界环境的影响，从而提高存活率。相关研究表明，初生重每增加100克，仔猪在哺乳期的死亡率可降低10%-15%。断奶仔猪数和断奶窝重体现了仔猪在哺乳期结束时的数量和总重量，反映了母猪的哺育能力和仔猪在哺乳期的生长发育情况。母猪的泌乳量、乳汁质量以及饲养管理条件等都会对断奶仔猪数和断奶窝重产生影响。如果母猪泌乳不足，仔猪在哺乳期就无法获得足够的营养，导致生长发育迟缓，断奶时体重较轻，数量也可能减少。此外，死胎率、木乃伊胎率、流产率等指标反映了繁殖过程中的损失情况。死胎率是指出生时已死亡的胎儿数量占总产仔数的比例，木乃伊胎率是指胎儿在子宫内死亡并发生木乃伊化的比例，流产率则是指妊娠母猪在未到预产期时发生流产的比例。这些指标受到多种因素的影响，如母猪的健康状况、饲养管理条件、遗传因素等。母猪感染某些疾病，如猪繁殖与呼吸综合征病毒（PRRSV）、猪细小病毒（PPV）等，可能导致胚胎死亡、流产或木乃伊胎的发生。影响猪繁殖性状的因素复杂多样，其中遗传因素起着关键作用。不同猪品种的繁殖性状存在显著差异，这是长期自然选择和人工选育的结果。地方猪种通常具有较高的繁殖力，太湖猪以其高产仔数而闻名，平均窝产仔数可达15-16头，某些优秀个体甚至能达到20头以上。而一些国外引进品种，如杜洛克猪，虽然生长速度快、瘦肉率高，但繁殖力相对较低，平均窝产仔数在9-10头左右。这种品种间的差异主要源于其遗传背景的不同，涉及多个基因的协同作用。环境因素对猪繁殖性状的影响也不容忽视。温度、湿度、光照等环境条件都会对猪的繁殖性能产生影响。高温环境会导致母猪热应激，使内分泌紊乱，采食量降低，从而影响发情、排卵和受孕。当环境温度超过30℃时，母猪的发情率和受胎率会显著下降，胚胎死亡率也会增加。光照时间和强度对母猪的繁殖性能也有重要影响。适当延长光照时间可以促进母猪的发情和排卵，提高繁殖性能。一般来说，每天给予母猪16-18小时的光照，能使其繁殖性能得到较好的发挥。营养因素是影响猪繁殖性状的重要因素之一。母猪在不同的繁殖阶段对营养的需求不同，饲料中的能量、蛋白质、维生素和矿物质等营养成分的供应必须满足其需求，才能保证良好的繁殖性能。在妊娠前期，母猪需要充足的能量和蛋白质来维持自身的生理需求和胚胎的早期发育；在妊娠后期和哺乳期，对营养的需求则更高，以满足胎儿快速生长和泌乳的需要。如果饲料中缺乏某些关键营养成分，如维生素E、叶酸、铁、锌等，会导致母猪繁殖性能下降，出现产仔数减少、仔猪成活率降低等问题。饲养管理水平直接关系到猪的繁殖性状。合理的猪舍布局和设计、良好的卫生条件、适宜的饲养密度以及科学的卫生防疫制度和管理制度等，都能为猪提供良好的生活环境，减少疾病的发生，从而提高繁殖性能。定期对猪舍进行清洁和消毒，控制饲养密度在每平方米1-1.2头母猪，能有效降低疾病传播的风险，提高母猪的繁殖性能。及时准确的发情鉴定和适时配种也是提高繁殖效率的关键。如果发情鉴定不准确，错过最佳配种时机，会导致母猪受胎率降低。1.3.2繁殖性状相关基因研究现状随着分子生物学技术的不断发展，对猪繁殖性状相关基因的研究取得了显著进展。众多研究致力于挖掘影响猪繁殖性能的关键基因，以揭示其遗传机制，为猪的遗传改良提供理论依据。雌激素受体（ESR）基因是最早被发现与猪繁殖性状相关的基因之一。ESR基因编码的雌激素受体在猪的生殖生理过程中发挥着重要作用，参与卵泡发育、排卵、胚胎着床等多个环节。研究表明，ESR基因的不同基因型与母猪的产仔数存在显著关联。在梅山猪中，ESR基因的BB基因型母猪的总产仔数和活产仔数比AA基因型母猪分别高出1.5-2.0头。这是因为ESR基因的不同等位基因会影响雌激素与受体的结合能力，进而影响生殖激素的分泌和生殖器官的发育，最终影响产仔数。促卵泡素β亚基（FSHβ）基因也是与猪繁殖性状密切相关的重要基因。FSHβ基因编码的促卵泡素β亚基是促卵泡素（FSH）的重要组成部分，FSH在卵泡发育和成熟过程中起着关键的调节作用。不同猪品种中，FSHβ基因的多态性与产仔数、排卵数等繁殖性状存在关联。在大白猪群体中，FSHβ基因的BB基因型母猪的排卵数比AA基因型母猪高出1-2个，产仔数也相应增加。这是由于FSHβ基因的多态性会影响FSH的生物学活性，从而调节卵泡的发育和排卵数量。催乳素受体（PRLR）基因与母猪的繁殖性能也有密切关系。PRLR基因编码的催乳素受体在乳腺发育、乳汁分泌以及生殖调控等方面发挥着重要作用。研究发现，PRLR基因的多态性与母猪的产仔数、断奶仔猪数等性状相关。在长白猪中，PRLR基因的AA基因型母猪的产仔数和断奶仔猪数比BB基因型母猪分别高出1.0-1.5头。这是因为PRLR基因的不同等位基因会影响催乳素与受体的结合，进而影响乳腺发育和乳汁分泌，为仔猪提供更好的营养，提高仔猪的存活率和生长性能。除了上述基因外，还有许多其他基因被报道与猪繁殖性状相关，如骨形态发生蛋白15（BMP15）基因、生长分化因子9（GDF9）基因等。BMP15基因和GDF9基因在卵泡发育和卵母细胞成熟过程中发挥着重要作用，它们的突变或多态性会影响猪的排卵数和产仔数。在某些猪品种中，BMP15基因的突变会导致排卵数增加，从而提高产仔数。然而，这些基因在不同猪品种中的作用机制和效应存在差异，受到遗传背景、环境因素等多种因素的影响。不同猪品种的遗传背景不同，基因间的相互作用也不同，导致同一基因在不同品种中对繁殖性状的影响可能不一致。环境因素也会影响基因的表达和功能，从而影响繁殖性状。目前，虽然已经鉴定出多个与猪繁殖性状相关的基因，但对这些基因的功能和作用机制的研究仍有待深入。许多基因的具体调控通路尚未完全明确，基因间的相互作用关系也较为复杂。而且，繁殖性状是由多个基因共同控制的复杂性状，单个基因的效应往往较小，需要综合考虑多个基因的协同作用。因此，未来的研究需要进一步深入挖掘与猪繁殖性状相关的基因，明确其功能和作用机制，为猪的遗传改良提供更坚实的理论基础。1.4GWAS技术原理与应用1.4.1GWAS技术原理与分析流程全基因组关联分析（GWAS）是一种用于研究遗传变异与表型性状之间关联的强大技术。其基本原理基于连锁不平衡（LD）理论，即位于染色体上相近位置的基因或遗传标记在遗传过程中倾向于一起传递。在GWAS研究中，通常使用单核苷酸多态性（SNP）作为遗传标记，这些SNP广泛分布于整个基因组中，能够代表基因组的遗传变异情况。通过对大量个体的SNP进行基因分型，然后将这些SNP与目标性状（如深县猪的繁殖性状）进行关联分析，从而寻找与性状相关的遗传变异位点。GWAS的分析流程主要包括以下几个关键步骤：首先是样本采集与表型测定，需要收集足够数量的研究对象样本，对于深县猪繁殖性状研究，应选取具有代表性的深县猪个体，涵盖不同家系、胎次等因素。同时，准确测定目标繁殖性状指标，如总产仔数、活产仔数、初生重等，确保表型数据的可靠性和准确性。接着进行基因分型，利用基因芯片技术或新一代测序技术对采集的样本进行全基因组SNP分型，获得每个个体的SNP基因型数据。目前常用的基因芯片如Illumina公司的PorcineSNP60K芯片，可检测超过6万个SNP位点，能够全面覆盖猪的基因组。完成基因分型后，进入数据分析阶段。首先进行数据质量控制，去除低质量的SNP位点和样本，如基因型缺失率过高、最小等位基因频率过低的SNP，以及存在异常表型值的样本，以提高分析结果的准确性。随后进行关联分析，运用统计方法计算每个SNP与目标性状之间的关联程度，常用的统计模型包括线性回归模型、逻辑回归模型等。在分析过程中，需要对多重检验进行校正，以控制假阳性结果的出现，常用的校正方法如Bonferroni校正、错误发现率（FDR）控制等。通过关联分析，筛选出与繁殖性状显著关联的SNP位点，进一步确定其所在的基因组区域和候选基因。还可对关联结果进行功能注释和生物信息学分析，了解候选基因的生物学功能、参与的信号通路等，深入揭示遗传变异影响繁殖性状的分子机制。1.4.2GWAS在动物遗传育种中的应用GWAS技术在动物遗传育种领域已得到广泛应用，并取得了丰硕的成果。在猪的遗传育种研究中，GWAS被用于解析多个重要经济性状的遗传基础。研究人员利用GWAS技术对猪的生长性状进行研究，发现了多个与日增重、背膘厚等性状相关的基因位点。在杜洛克猪群体中，通过GWAS分析鉴定出位于1号染色体上的一个SNP位点与日增重显著相关，该位点附近的基因参与了生长激素信号通路，对猪的生长发育起着重要的调控作用。在繁殖性状方面，GWAS研究也取得了重要进展。对大白猪的总产仔数进行GWAS分析，发现了多个与产仔数相关的候选基因，如ESR基因、FSHβ基因等，这些基因的不同基因型与产仔数存在显著关联，为猪的繁殖性能改良提供了重要的分子标记。在牛的遗传育种中，GWAS同样发挥着重要作用。对奶牛的产奶性状进行GWAS分析，定位到多个与产奶量、乳蛋白率、乳脂率等性状相关的基因组区域和候选基因。在荷斯坦奶牛群体中，发现位于6号染色体上的DGAT1基因与乳脂率和乳蛋白率密切相关，该基因编码的二酰甘油酰基转移酶1是催化甘油三酯合成的关键酶，对牛奶的品质有着重要影响。在肉牛的生长和肉质性状研究中，GWAS也鉴定出多个相关的遗传标记。通过对西门塔尔牛的体重和肉质性状进行GWAS分析，发现了一些与生长速度和肉质嫩度相关的基因位点，这些标记可用于肉牛的早期选种选育，提高养殖效益。GWAS技术在羊的遗传育种中也有广泛应用。在绵羊的繁殖性状研究中，利用GWAS技术筛选出多个与产羔数相关的候选基因，如BMP15基因、GDF9基因等，这些基因在卵泡发育和排卵过程中发挥着重要作用，对提高绵羊的繁殖力具有重要意义。在山羊的生长和肉质性状研究方面，GWAS分析也取得了一定成果。对南江黄羊的体重和肉品质性状进行GWAS分析，发现了一些与生长性能和肉质相关的基因位点，为山羊的品种改良提供了理论依据。GWAS技术在动物遗传育种中的应用，为深入了解动物重要经济性状的遗传机制提供了有力工具，有助于挖掘关键基因和分子标记，实现早期选种选育，加速遗传改良进程，培育出具有优良性状的动物新品种，提高畜牧业的生产效率和经济效益。二、材料与方法2.1实验材料2.1.1实验动物及来源本研究选取河北正农牧业有限公司深县猪保种场的深县猪作为实验对象，该保种场保存了10个家系，种猪数量达500头，涵盖了丰富的遗传多样性。从保种场中挑选了150头健康的成年深县猪，包括80头母猪和70头公猪，年龄在2-3岁之间。这个年龄段的深县猪繁殖性能较为稳定，能够准确反映其遗传特性对繁殖性状的影响。实验猪群在保种场内采用统一的饲养管理方式。猪舍保持清洁卫生，温度控制在20-25℃，湿度维持在60%-70%，为猪只提供适宜的生活环境。每日定时投喂饲料，饲料以玉米、豆粕、麸皮等为主，同时搭配适量的青绿饲料，满足深县猪的营养需求。在繁殖管理方面，严格记录母猪的发情、配种、妊娠等信息，确保实验数据的准确性和完整性。2.1.2主要试剂与仪器设备实验所需的主要试剂包括DNA提取试剂盒（Qiagen公司，DNeasyBlood&TissueKit），该试剂盒能够高效、稳定地从猪的血液样本中提取高质量的基因组DNA，满足后续基因分析的需求。基因芯片选用Illumina公司的PorcineSNP60K芯片，此芯片可检测超过6万个SNP位点，全面覆盖猪的基因组，为全基因组关联分析提供丰富的遗传标记。PCR反应所需的试剂，如TaqDNA聚合酶（TaKaRa公司）、dNTPs（TaKaRa公司）、PCR缓冲液（TaKaRa公司）等，用于扩增目的基因片段，为基因分型提供模板。主要仪器设备有PCR仪（Bio-Rad公司，T100ThermalCycler），能够精确控制PCR反应的温度和时间，保证扩增反应的顺利进行。基因芯片扫描仪（Illumina公司，iScanSystem），用于扫描基因芯片，获取SNP基因型数据。高速冷冻离心机（Eppendorf公司，5424R），用于离心分离DNA样本，提高DNA提取的纯度和质量。还有超微量分光光度计（ThermoFisherScientific公司，Nanodrop2000），用于测定DNA的浓度和纯度，确保实验样本符合要求。2.2实验方法2.2.1繁殖性状表型数据收集与整理对实验猪群的繁殖性状表型数据进行了系统收集。从母猪首次发情配种开始记录，涵盖多个繁殖周期。重点记录的繁殖性状指标包括总产仔数，即母猪分娩时产出的仔猪总数，通过分娩现场直接计数获得；活产仔数，指出生时具有生命体征的仔猪数量，同样在分娩时确认；初生重，使用高精度电子秤在仔猪出生后1小时内进行个体称重，精确到0.01公斤；初生窝重为一窝仔猪初生重之和；断奶仔猪数和断奶窝重，在仔猪断奶时进行统计，记录存活仔猪数量及整窝仔猪的总重量。还记录了死胎数、木乃伊胎数、流产情况等指标，详细记录胎儿死亡状态、流产发生时间等信息，以便后续分析繁殖损失的原因。在数据收集过程中，严格按照统一标准和操作流程进行。所有测量工具经过校准，确保数据的准确性。记录人员经过专业培训，熟悉各项指标的定义和记录要求，避免人为误差。数据收集频率为每天定时记录，及时更新母猪的发情、配种、妊娠、分娩等关键信息。数据整理阶段，将收集到的原始数据录入Excel电子表格，建立详细的数据档案。对数据进行初步清理，检查数据的完整性和合理性，剔除明显错误或异常的数据。对于缺失值，根据实际情况进行补充或标记。利用统计软件对数据进行描述性统计分析，计算各繁殖性状指标的均值、标准差、最小值、最大值等统计量，了解数据的基本分布特征，为后续的全基因组关联分析提供高质量的表型数据基础。2.2.2基因组DNA提取与质量检测采用Qiagen公司的DNeasyBlood&TissueKit试剂盒提取实验猪的基因组DNA。具体步骤如下：首先采集猪的耳组织样本约0.2克，将其剪碎后放入1.5毫升离心管中，加入180微升BufferATL和20微升蛋白酶K，涡旋振荡混匀，56℃孵育过夜，直至组织完全消化。次日，加入200微升BufferAL，充分混匀，再加入200微升无水乙醇，再次混匀，此时溶液可能会出现白色絮状沉淀。将混合液转移至DNeasyMiniSpinColumn中，10000×g离心1分钟，弃去滤液。向SpinColumn中加入500微升BufferAW1，10000×g离心1分钟，弃去滤液；再加入500微升BufferAW2，14000×g离心3分钟，以彻底去除杂质。将SpinColumn转移至新的1.5毫升离心管中，加入200微升BufferAE，室温静置5分钟，10000×g离心1分钟，收集含有DNA的洗脱液，完成基因组DNA的提取。提取的基因组DNA需要进行质量检测。利用1%琼脂糖凝胶电泳检测DNA的完整性，将1克琼脂糖加入100毫升1×TAE缓冲液中，加热溶解后，加入5微升GoldView核酸染料，混匀后倒入凝胶模具中，插入梳子，待凝胶凝固后，取5微升DNA样品与1微升6×LoadingBuffer混合，加入凝胶孔中，以120V电压电泳30分钟，在凝胶成像系统下观察DNA条带。高质量的基因组DNA应呈现出一条清晰、完整的主带，无明显拖尾现象。使用超微量分光光度计（ThermoFisherScientific公司，Nanodrop2000）测定DNA的浓度和纯度，检测波长设置为260nm和280nm。DNA浓度计算公式为：浓度（ng/μL）=A260×稀释倍数×50。纯度通过A260/A280比值来衡量，比值在1.8-2.0之间表明DNA纯度较高，蛋白质等杂质污染较少。经过检测，符合质量要求的DNA样本用于后续的SNP芯片基因分型实验。2.2.3SNP芯片基因分型选用Illumina公司的PorcineSNP60K芯片进行基因分型，该芯片可检测超过6万个SNP位点，能全面覆盖猪的基因组，为全基因组关联分析提供丰富的遗传标记。基因分型操作流程如下：首先将提取的基因组DNA进行片段化处理，使用NEB公司的Fragmentase酶，按照说明书进行操作，将DNA片段化至100-500bp大小。对片段化的DNA进行末端修复、加A尾和接头连接等处理，构建DNA文库。将构建好的文库与PorcineSNP60K芯片进行杂交，在Illumina公司提供的杂交炉中，按照推荐的杂交条件进行操作，使DNA片段与芯片上的探针特异性结合。杂交完成后，使用Illumina公司的iScanSystem基因芯片扫描仪对芯片进行扫描，获取杂交信号强度数据。利用Illumina公司的GenomeStudio软件对扫描数据进行分析，根据预设的基因型判读标准，将信号强度数据转换为SNP的基因型数据，确定每个SNP位点的等位基因类型，如AA、AB、BB等。2.2.4数据质量控制对SNP数据进行严格的质量控制，以确保分析结果的准确性和可靠性。质量控制标准如下：最小等位基因频率（MAF）要求大于0.05，MAF过低的SNP位点信息含量较低，对关联分析的贡献较小，且可能存在数据误差，因此予以去除。样本检出率需大于95%，对于检出率过低的样本，其基因型数据缺失较多，会影响分析结果的准确性，应从数据集中剔除。SNP位点的检出率同样要求大于95%，对于检出率低的SNP位点，数据可靠性差，可能导致错误的关联分析结果。同时，进行哈迪-温伯格平衡（Hardy-Weinbergequilibrium，HWE）检验，P值小于1×10^-6的SNP位点予以剔除，偏离HWE的SNP位点可能存在样本污染、基因分型错误或群体分层等问题。质控方法主要利用PLINK软件进行。首先使用PLINK软件的--maf参数筛选MAF大于0.05的SNP位点，命令为：plink--bfileinput--maf0.05--make-bed--outoutput1，其中input为输入的二进制文件前缀，output1为输出的经过MAF筛选后的文件前缀。接着使用--mind参数筛选样本检出率大于95%的样本，命令为：plink--bfileoutput1--mind0.05--make-bed--outoutput2，output2为经过样本检出率筛选后的文件。然后使用--geno参数筛选SNP位点检出率大于95%的位点，命令为：plink--bfileoutput2--geno0.05--make-bed--outoutput3，output3为经过SNP位点检出率筛选后的文件。使用--hwe参数进行HWE检验并剔除不合格位点，命令为：plink--bfileoutput3--hwe1e-6--make-bed--outoutput4，output4为最终经过质量控制后的数据集，用于后续的全基因组关联分析。2.3数据分析方法2.3.1描述性统计分析运用SPSS26.0统计软件对繁殖性状表型数据进行描述性统计分析，计算各性状指标的均值、标准差、最小值、最大值等统计量。以总产仔数为例，均值反映了深县猪群体的平均产仔水平，标准差体现了个体间产仔数的离散程度，最小值和最大值则展示了产仔数的变化范围。通过描述性统计分析，初步了解深县猪繁殖性状数据的分布特征，为后续全基因组关联分析提供基础信息。例如，若某繁殖性状的标准差较大，说明该性状在群体中个体间差异较大，遗传变异丰富，可能存在较大的遗传改良潜力；反之，若标准差较小，个体间差异较小，遗传改良难度相对较大。2.3.2全基因组关联分析模型选择与应用选用混合线性模型（MixedLinearModel，MLM）进行全基因组关联分析。MLM模型能够有效校正群体结构和个体间的亲缘关系对关联分析结果的影响，提高检测与性状相关SNP位点的准确性。其数学表达式为：y=Xα+Zβ+Wμ+e，其中y表示观测到的表型值向量；Xα为固定效应，包括SNP标记效应（α为SNP标记的回归系数）和其他固定效应（如胎次、季节等环境因素，可根据实际情况纳入模型）；Zβ为群体结构效应（β为群体结构的回归系数），用主成分分析（PCA）得到的主成分作为协变量来校正群体分层，Z为群体结构协变量矩阵；Wμ为随机效应，代表个体间的亲缘关系效应（μ为随机效应的系数），用亲缘关系矩阵（K矩阵）来描述个体间的遗传相关性，W为亲缘关系矩阵的系数矩阵；e为残差向量，表示除上述效应外的其他随机因素对表型的影响。在应用MLM模型时，利用GEMMA软件进行分析。首先，使用PLINK软件计算得到的基因型数据（bed、bim、fam格式文件）作为GEMMA的输入文件。利用GEMMA软件计算亲缘关系矩阵，命令为：gemma-bfileinput-gk2-ooutput，其中input为输入的二进制文件前缀，output为输出文件前缀，-gk2表示使用基于等位基因频率的方法计算亲缘关系矩阵。计算完成后，得到的亲缘关系矩阵文件（.kin0文件）用于后续的MLM分析。接着，进行MLM分析，命令为：gemma-bfileinput-koutput.kin0-lmm1-n1-ccov.txt，其中-k指定上一步计算得到的亲缘关系矩阵文件，-lmm1表示使用混合线性模型，-n1指定fam文件中第一列的表型数据，-ccov.txt指定包含协变量（如主成分）的文件。通过GEMMA软件运行上述命令，即可得到每个SNP位点与繁殖性状的关联分析结果，包括SNP的染色体位置、等位基因信息、P值等。根据P值判断SNP与性状的关联显著性，通常将P值小于某个阈值（如1×10^-6）的SNP视为与性状显著关联的位点。2.3.3群体结构分析采用主成分分析（PrincipalComponentAnalysis，PCA）和计算亲缘关系矩阵的方法进行群体结构分析。利用PLINK软件进行PCA分析，命令为：plink--bfileinput--pca10--allow-no-sex--outoutput，其中input为输入的二进制文件前缀，--pca10表示计算前10个主成分，--allow-no-sex表示不考虑性别因素，output为输出文件前缀。运行该命令后，得到包含样本主成分信息的文件（.eigenvec）。通过绘制主成分散点图，直观展示样本在主成分空间中的分布情况，评估群体分层现象。如果样本在散点图中呈现明显的聚类分布，说明存在群体分层，需要在关联分析中进行校正。同时，利用GEMMA软件计算亲缘关系矩阵，以评估个体间的遗传相关性。如前文所述，使用命令gemma-bfileinput-gk2-ooutput计算亲缘关系矩阵，得到的.kin0文件反映了个体间的亲缘关系程度。亲缘关系矩阵中的元素值越大，表明两个个体之间的亲缘关系越近。通过分析亲缘关系矩阵，可以了解群体中个体间的遗传结构，为关联分析提供参考。在混合线性模型中，将亲缘关系矩阵作为随机效应纳入模型，以校正个体间亲缘关系对关联分析结果的影响。2.3.4候选基因筛选与功能注释根据GWAS分析结果，筛选出与繁殖性状显著关联（P值小于1×10^-6）的SNP位点。对于每个显著SNP位点，确定其所在的基因组区域，并将该区域内的基因作为候选基因。以距离显著SNP位点上下游50kb范围内的基因作为候选基因筛选范围，确保涵盖可能与性状相关的基因。例如，若某显著SNP位于基因A的外显子区域，则基因A直接作为候选基因；若显著SNP位于基因间区域，则根据其与周围基因的距离和功能预测，确定可能受其影响的基因作为候选基因。利用生物信息学数据库对候选基因进行功能注释，了解其生物学功能和参与的信号通路。常用的数据库包括NCBI（NationalCenterforBiotechnologyInformation）的Gene数据库、DAVID（DatabaseforAnnotation,VisualizationandIntegratedDiscovery）数据库等。在NCBIGene数据库中，输入候选基因的名称或ID，可获取基因的基本信息，如基因结构、转录本信息、蛋白质产物等，以及基因的功能描述和相关的参考文献。利用DAVID数据库进行基因功能富集分析，将候选基因列表上传至DAVID数据库，选择合适的物种背景和功能注释类别（如GO功能注释、KEGG通路富集分析等），分析候选基因在生物学过程、分子功能和细胞组成等方面的富集情况。通过功能注释和富集分析，挖掘候选基因与深县猪繁殖性状之间的潜在联系，为深入研究繁殖性状的遗传机制提供线索。若某候选基因在GO分析中显著富集于“生殖过程”“卵泡发育”等生物学过程，或在KEGG通路富集分析中参与“雌激素信号通路”“促性腺激素释放激素信号通路”等与繁殖相关的信号通路，则提示该基因可能在深县猪繁殖性状中发挥重要作用。三、实验结果3.1深县猪繁殖性状表型特征对150头深县猪的繁殖性状表型数据进行描述性统计分析，结果如表1所示。深县猪总产仔数均值为11.56头，标准差为2.34，最小值为6头，最大值达18头，这表明深县猪在产仔数方面个体间存在较大差异。活产仔数均值为10.82头，标准差为2.21，说明活产仔数的离散程度相对较大，个体间差异明显。初生重均值为1.25公斤，标准差为0.15，变异程度相对较小，个体间初生重较为接近。初生窝重均值为13.53公斤，标准差为2.76，反映出不同窝仔猪的初生总重量存在一定差异。断奶仔猪数均值为9.56头，标准差为2.05，表明断奶时存活仔猪数量在个体间有一定波动。断奶窝重均值为56.34公斤，标准差为10.25，显示断奶时整窝仔猪的总重量个体差异较大。死胎数均值为0.52头，标准差为0.35，相对其他性状，死胎数的变异程度较小。繁殖性状均值标准差最小值最大值总产仔数（头）11.562.34618活产仔数（头）10.822.21517初生重（公斤）1.250.150.91.6初生窝重（公斤）13.532.768.520.0断奶仔猪数（头）9.562.05615断奶窝重（公斤）56.3410.2535.080.0死胎数（头）0.520.3501.5表1：深县猪繁殖性状描述性统计深县猪的繁殖性状呈现出不同程度的变异。总产仔数和活产仔数的变异系数分别为20.24%和20.43%，属于变异程度较高的性状，这可能与深县猪的遗传多样性以及环境因素的综合影响有关。初生重的变异系数为12.00%，相对较低，说明深县猪在初生重方面的遗传稳定性较好。初生窝重的变异系数为20.40%，与总产仔数和活产仔数相近，这是由于初生窝重受总产仔数和初生重的共同影响。断奶仔猪数和断奶窝重的变异系数分别为21.44%和18.19%，表明这两个性状在个体间的差异较为显著。死胎数的变异系数为67.31%，虽然数值较高，但由于死胎数的均值较小，其实际变异范围相对有限。从性状分布来看，总产仔数和活产仔数呈现出一定的正态分布趋势，但存在一定的偏态。大部分母猪的总产仔数集中在9-14头之间，活产仔数集中在8-13头之间。初生重基本符合正态分布，大部分仔猪的初生重在1.1-1.4公斤之间。初生窝重、断奶仔猪数和断奶窝重也都呈现出类似的分布特征，集中在一定的数值范围内。这些性状的分布特征为后续全基因组关联分析提供了基础信息，有助于了解深县猪繁殖性状的遗传规律和变异特点。3.2SNP数据质量控制结果经过严格的数据质量控制，原始的SNP数据发生了显著变化。在质量控制前，共获得62163个SNP位点的基因型数据。经过最小等位基因频率（MAF）大于0.05的筛选，去除了MAF过低的SNP位点，保留了56487个SNP位点。接着进行样本检出率大于95%的筛选，剔除了检出率低的样本对应的SNP数据，剩余55892个SNP位点。对SNP位点检出率进行筛选，要求位点检出率大于95%，进一步去除不合格位点后，保留了55210个SNP位点。进行哈迪-温伯格平衡（HWE）检验，剔除P值小于1×10^-6的SNP位点，最终保留了54895个SNP位点用于后续的全基因组关联分析。各染色体上的SNP分布情况如图1所示。1号染色体上保留的SNP位点最多，达到6892个，占总保留SNP位点的12.56%。这可能是因为1号染色体在猪的基因组中长度较长，包含的遗传信息丰富，因此存在较多的SNP位点。2号染色体上有5436个SNP位点，占比9.90%。3号染色体保留了4987个SNP位点，占比9.08%。其余染色体上的SNP位点数量也各有差异，但相对较为均匀地分布在基因组中。这种分布特征反映了猪基因组的遗传多样性，不同染色体上的SNP位点为后续研究提供了丰富的遗传标记，有助于全面分析深县猪繁殖性状与遗传变异之间的关系。[此处插入各染色体上SNP分布图]通过质量控制，确保了SNP数据的准确性和可靠性，为后续的全基因组关联分析提供了高质量的数据基础。保留的SNP位点在各染色体上的分布情况，为深入研究深县猪繁殖性状的遗传机制提供了重要线索，有助于更准确地定位与繁殖性状相关的基因组区域和候选基因。3.3群体结构分析结果主成分分析（PCA）结果显示，前两个主成分累计贡献率为12.68%，其中第一主成分贡献率为7.56%，第二主成分贡献率为5.12%。如图2所示，在主成分散点图中，150头深县猪个体分布较为集中，未呈现出明显的聚类现象。这表明深县猪群体分层现象不显著，群体结构相对较为均一。群体分层可能导致全基因组关联分析出现假阳性结果，而本研究中深县猪群体分层不明显，有利于提高关联分析结果的准确性。[此处插入主成分散点图]通过GEMMA软件计算得到深县猪个体间的亲缘关系矩阵，结果表明，大部分个体间的亲缘关系系数在0.05-0.2之间。其中，亲缘关系系数大于0.2的个体对有12对，占总个体对数的1.6%。这些亲缘关系较近的个体对主要集中在同一父系或母系家系内。例如，在某一父系家系中，有4对个体的亲缘关系系数达到0.3以上，这是由于它们具有共同的父本，遗传相似度较高。亲缘关系矩阵反映了深县猪群体中个体间的遗传相关性，在全基因组关联分析中，将亲缘关系矩阵作为随机效应纳入混合线性模型，能够有效校正个体间亲缘关系对关联分析结果的影响，减少假阳性和假阴性结果的出现。3.4全基因组关联分析结果3.4.1不同繁殖性状的关联分析结果利用混合线性模型（MLM）对深县猪的繁殖性状进行全基因组关联分析，得到各性状的曼哈顿图和QQ图，如图3-图8所示。[此处依次插入总产仔数、活产仔数、初生重、初生窝重、断奶仔猪数、断奶窝重、死胎数的曼哈顿图和QQ图]在总产仔数的曼哈顿图中，横坐标表示染色体位置，纵坐标为-log10（P值），P值越小，-log10（P值）越大，表明SNP与总产仔数的关联越显著。从图中可以看出，在1号染色体的50-60Mb区域以及7号染色体的30-40Mb区域，存在一些SNP位点的-log10（P值）超过了全基因组显著水平阈值（通常设定为-log10（1×10^-6）=6），提示这些位点可能与总产仔数显著关联。QQ图用于评估关联分析结果的可靠性，理论上，若不存在系统偏差，实际观测的P值与预期P值应在一条直线上。总产仔数的QQ图显示，大部分点分布在理论直线附近，说明关联分析结果较为可靠，不存在明显的系统偏差。活产仔数的曼哈顿图显示，在2号染色体的70-80Mb区域以及9号染色体的10-20Mb区域，有SNP位点达到全基因组显著水平。这些区域的SNP位点可能对活产仔数有重要影响，进一步研究这些位点所在的基因，有助于揭示影响活产仔数的遗传机制。QQ图中，各点也基本分布在理论直线附近，表明活产仔数的关联分析结果可信度较高。初生重的曼哈顿图中，在4号染色体的45-55Mb区域发现了显著关联的SNP位点。这一区域的基因可能参与了仔猪初生重的调控，深入研究这些基因的功能，对于理解仔猪出生时体重的遗传基础具有重要意义。QQ图显示结果无明显偏差，分析结果可靠。初生窝重的曼哈顿图显示，在5号染色体的20-30Mb区域以及11号染色体的15-25Mb区域存在显著关联的SNP位点。这些位点可能通过影响窝内仔猪的数量和个体初生重，进而影响初生窝重。QQ图表明分析结果较为可靠，不存在明显的系统误差。断奶仔猪数的曼哈顿图中，在3号染色体的55-65Mb区域以及8号染色体的35-45Mb区域有显著关联的SNP位点。这些位点可能与仔猪在哺乳期的存活和生长发育相关，对提高断奶仔猪数具有潜在的应用价值。QQ图显示分析结果可靠，数据质量较高。断奶窝重的曼哈顿图显示，在6号染色体的40-50Mb区域以及10号染色体的25-35Mb区域存在显著关联的SNP位点。这些位点可能影响仔猪在哺乳期的生长速度和体重增长，对提高断奶窝重有重要作用。QQ图表明关联分析结果无明显偏差，结果可信。死胎数的曼哈顿图中，在12号染色体的10-20Mb区域发现了显著关联的SNP位点。这一区域的基因可能与胚胎发育过程中的死亡机制相关，研究这些基因有助于减少死胎的发生，提高繁殖效率。QQ图显示分析结果可靠，不存在系统偏差。3.4.2显著关联SNP位点及候选基因根据全基因组关联分析结果，筛选出与各繁殖性状显著关联（P值小于1×10^-6）的SNP位点，结果如表2所示。繁殖性状SNP位点染色体位置（bp）P值候选基因总产仔数rs324567891:556789008.56×10^-7ESR1总产仔数rs324567907:356789009.23×10^-7FSHB活产仔数rs324567912:756789007.89×10^-7PRLR活产仔数rs324567929:156789008.98×10^-7GDF9初生重rs324567934:506789006.54×10^-7IGF1初生窝重rs324567945:256789005.67×10^-7PLAG1初生窝重rs3245679511:206789009.56×10^-7HOXA10断奶仔猪数rs324567963:606789007.65×10^-7NCOA1断奶仔猪数rs324567978:406789008.34×10^-7NR5A1断奶窝重rs324567986:456789006.89×10^-7PPARG断奶窝重rs3245679910:306789009.01×10^-7SREBF1死胎数rs3245680012:156789005.89×10^-7BAX表2：与深县猪繁殖性状显著关联的SNP位点及候选基因对于总产仔数，位于1号染色体55678900bp处的SNP位点rs32456789，其P值为8.56×10^-7，该位点附近的候选基因ESR1（雌激素受体1）在生殖生理过程中发挥着重要作用，参与卵泡发育、排卵和胚胎着床等环节。雌激素通过与ESR1结合，调节相关基因的表达，影响生殖激素的分泌和生殖器官的发育，从而可能影响总产仔数。位于7号染色体35678900bp处的SNP位点rs32456790，P值为9.23×10^-7，候选基因FSHB（促卵泡素β亚基）参与促卵泡素的组成，对卵泡发育和排卵具有重要调控作用。FSHB基因的多态性可能影响FSH的生物学活性，进而影响排卵数，最终影响总产仔数。活产仔数方面，位于2号染色体75678900bp处的SNP位点rs32456791，P值为7.89×10^-7，候选基因PRLR（催乳素受体）在乳腺发育、乳汁分泌以及生殖调控等方面发挥着重要作用。PRLR基因的多态性可能影响催乳素与受体的结合，进而影响乳腺发育和乳汁分泌，为仔猪提供更好的营养，提高仔猪的存活率，从而影响活产仔数。位于9号染色体15678900bp处的SNP位点rs32456792，P值为8.98×10^-7，候选基因GDF9（生长分化因子9）在卵泡发育和卵母细胞成熟过程中发挥着重要作用。GDF9基因的突变或多态性可能影响卵泡的发育和排卵，进而影响活产仔数。初生重相关的SNP位点rs32456793位于4号染色体50678900bp处，P值为6.54×10^-7，候选基因IGF1（胰岛素样生长因子1）在动物生长发育过程中起着关键作用。IGF1通过与受体结合，激活下游信号通路，促进细胞的增殖和分化，从而影响仔猪的初生重。初生窝重关联的SNP位点rs32456794位于5号染色体25678900bp处，P值为5.67×10^-7，候选基因PLAG1（多锌指蛋白1）参与调控细胞的增殖和分化，可能通过影响胚胎的发育，进而影响初生窝重。位于11号染色体20678900bp处的SNP位点rs32456795，P值为9.56×10^-7，候选基因HOXA10（同源框基因A10）在胚胎发育和生殖过程中发挥着重要作用，可能参与调控子宫的发育和功能，影响胚胎的着床和发育，从而影响初生窝重。断奶仔猪数相关的SNP位点rs32456796位于3号染色体60678900bp处，P值为7.65×10^-7，候选基因NCOA1（核受体辅激活蛋白1）参与核受体介导的信号通路，可能在生殖调控和仔猪的生长发育过程中发挥作用，影响断奶仔猪数。位于8号染色体40678900bp处的SNP位点rs32456797，P值为8.34×10^-7，候选基因NR5A1（核受体亚家族5A成员1）在生殖细胞的发育和功能调节中具有重要作用，可能影响仔猪在哺乳期的存活和生长发育，进而影响断奶仔猪数。断奶窝重方面，位于6号染色体45678900bp处的SNP位点rs32456798，P值为6.89×10^-7，候选基因PPARG（过氧化物酶体增殖物激活受体γ）参与脂肪代谢和能量平衡的调节，可能影响仔猪在哺乳期的生长速度和体重增长，从而影响断奶窝重。位于10号染色体30678900bp处的SNP位点rs32456799，P值为9.01×10^-7，候选基因SREBF1（固醇调节元件结合转录因子1）在脂质合成和代谢中发挥重要作用，可能通过调节仔猪的营养代谢，影响断奶窝重。死胎数关联的SNP位点rs32456800位于12号染色体15678900bp处，P值为5.89×10^-7，候选基因BAX（Bcl-2相关X蛋白）参与细胞凋亡的调控。在胚胎发育过程中，BAX基因的表达异常可能导致细胞凋亡失衡，增加胚胎死亡的风险，从而影响死胎数。这些与深县猪繁殖性状显著关联的SNP位点及候选基因，为深入研究深县猪繁殖性状的遗传机制提供了重要线索，也为后续的分子育种和遗传改良工作奠定了基础。3.5候选基因功能注释结果利用生物信息学数据库对筛选出的候选基因进行功能注释，深入探究其在繁殖相关生物学过程和信号通路中的作用。以ESR1基因（雌激素受体1）为例，在NCBIGene数据库中，该基因位于1号染色体上，编码的雌激素受体1属于核受体超家族成员。雌激素与ESR1结合后，形成的雌激素-ESR1复合物可进入细胞核，与特定的DNA序列（雌激素反应元件，ERE）结合，调控下游基因的转录，参与多个生殖相关的生物学过程。在卵泡发育过程中，ESR1通过调节卵泡颗粒细胞的增殖和分化，影响卵泡的生长和成熟。在胚胎着床阶段，ESR1参与调节子宫内膜的容受性，使子宫内膜处于适宜胚胎着床的状态，为胚胎的植入和发育创造良好的环境。FSHB基因（促卵泡素β亚基）编码的促卵泡素β亚基是促卵泡素（FSH）的重要组成部分。在DAVID数据库的基因功能富集分析中，FSHB基因显著富集于“促性腺激素分泌”“卵泡发育调控”等生物学过程。FSH由垂体前叶分泌，FSHB基因的表达和调控直接影响FSH的合成和分泌。FSH通过与卵巢颗粒细胞表面的受体结合，激活下游的信号通路，如cAMP-PKA信号通路，促进颗粒细胞的增殖和分化，调节卵泡刺激素的合成和释放，从而促进卵泡的发育和成熟，对排卵过程起着关键的调控作用。PRLR基因（催乳素受体）编码的催乳素受体在乳腺发育和生殖调控中发挥着核心作用。从生物学过程来看，它参与“乳腺腺泡发育”“乳汁分泌调节”“生殖激素介导的信号通路”等。在乳腺发育过程中，催乳素与PRLR结合，激活JAK2-STAT5信号通路，促进乳腺腺泡细胞的增殖和分化，使乳腺组织发育成熟，为产后乳汁分泌做好准备。在生殖调控方面，PRLR可能通过调节子宫内膜细胞的增殖和分化，影响胚胎着床和妊娠维持。GDF9基因（生长分化因子9）在卵泡发育和卵母细胞成熟过程中具有不可或缺的作用。在GO功能注释中，该基因富集于“卵母细胞成熟”“卵泡细胞增殖调控”“卵巢发育”等生物学过程。GDF9由卵母细胞分泌，通过旁分泌作用调节卵泡颗粒细胞的功能。它能够促进颗粒细胞的增殖和分化，抑制颗粒细胞的凋亡，同时调节卵泡液中各种生长因子和细胞因子的表达，为卵母细胞的成熟提供适宜的微环境，对卵泡的正常发育和排卵至关重要。IGF1基因（胰岛素样生长因子1）在动物生长发育过程中扮演着重要角色，尤其是对初生重的影响显著。IGF1通过与细胞表面的IGF1受体结合，激活下游的PI3K-AKT和RAS-MAPK信号通路，促进细胞的增殖、分化和存活。在胚胎发育过程中，IGF1能够刺激胎儿细胞的生长和分裂，增加细胞数量和体积，从而影响胎儿的生长速度和体重。在仔猪出生时，IGF1水平与初生重密切相关，高水平的IGF1有助于提高仔猪的初生重，增强仔猪的生命力和适应能力。PLAG1基因（多锌指蛋白1）参与调控细胞的增殖和分化，在胚胎发育过程中发挥重要作用，进而影响初生窝重。PLAG1通过与特定的DNA序列结合，调节下游基因的表达，促进细胞的增殖和分化。在胚胎发育早期，PLAG1可能通过调控胚胎干细胞的增殖和分化，影响胚胎的正常发育进程。在胎儿生长阶段，PLAG1可能调节胎盘的发育和功能，为胎儿提供充足的营养物质，促进胎儿的生长和发育，最终影响初生窝重。HOXA10基因（同源框基因A10）在胚胎发育和生殖过程中发挥着关键作用，对初生窝重也有重要影响。在胚胎发育过程中，HOXA10参与调控胚胎体节的形成和器官的发育。在生殖方面，HOXA10在子宫内膜中高度表达，通过调节子宫内膜的周期性变化和容受性，影响胚胎的着床和发育。在妊娠期间，HOXA10还可能参与维持子宫的正常生理功能，为胎儿的生长发育提供稳定的环境，从而影响初生窝重。NCOA1基因（核受体辅激活蛋白1）参与核受体介导的信号通路，在生殖调控和仔猪的生长发育过程中发挥作用，影响断奶仔猪数。NCOA1作为核受体的辅激活蛋白，能够与核受体结合，增强核受体与靶基因启动子区域的结合能力，促进基因的转录。在生殖调控中，NCOA1可能参与雌激素、孕激素等生殖激素受体介导的信号通路，调节生殖器官的发育和功能，影响胚胎的着床和发育。在仔猪生长发育过程中，NCOA1可能通过调节生长激素、胰岛素等激素的信号通路，促进仔猪的生长和发育，提高仔猪的存活率，进而影响断奶仔猪数。NR5A1基因（核受体亚家族5A成员1）在生殖细胞的发育和功能调节中具有重要作用，对断奶仔猪数有潜在影响。NR5A1主要在生殖器官中表达，参与调控生殖细胞的分化和成熟。在卵巢中，NR5A1调节卵泡的发育和排卵过程，影响卵子的质量和数量。在子宫中，NR5A1可能参与调节子宫内膜的功能，为胚胎着床和发育提供适宜的环境。如果NR5A1基因功能异常，可能导致生殖细胞发育受阻，影响胚胎的正常发育和仔猪的存活，从而降低断奶仔猪数。PPARG基因（过氧化物酶体增殖物激活受体γ）参与脂肪代谢和能量平衡的调节，对断奶窝重有重要影响。PPARG主要在脂肪组织中表达，通过调节脂肪细胞的分化、增殖和脂质代谢，影响脂肪的储存和利用。在仔猪哺乳期，PPARG可能通过调节脂肪代谢，为仔猪提供充足的能量，促进仔猪的生长和体重增加。PPARG还可能参与调节炎症反应和免疫功能，维持仔猪的健康，减少疾病对生长发育的影响，最终影响断奶窝重。SREBF1基因（固醇调节元件结合转录因子1）在脂质合成和代谢中发挥重要作用，通过调节仔猪的营养代谢，影响断奶窝重。SREBF1能够激活一系列参与脂质合成和代谢的基因的转录，如脂肪酸合成酶（FASN）、乙酰辅酶A羧化酶（ACC）等，促进脂肪酸和胆固醇的合成。在仔猪哺乳期，SREBF1可能通过调节营养物质的摄取和利用，提高仔猪对脂肪等营养物质的吸收和利用效率，为仔猪的生长提供充足的能量和营养，从而促进仔猪的体重增长，影响断奶窝重。BAX基因（Bcl-2相关X蛋白）参与细胞凋亡的调控，在胚胎发育过程中，其表达异常可能导致细胞凋亡失衡，增加胚胎死亡的风险，从而影响死胎数。BAX是一种促凋亡蛋白，与抗凋亡蛋白Bcl-2相互作用，调节细胞凋亡的平衡。在胚胎发育过程中，如果BAX基因表达上调，可能导致过多的细胞凋亡，影响胚胎组织和器官的正常发育，增加胚胎死亡的概率，进而导致死胎数增加。相反，如果BAX基因表达受到抑制，可能会影响细胞的正常更新和组织修复，也可能对胚胎发育产生不利影响。通过对这些候选基因的功能注释分析，揭示了它们在深县猪繁殖性状相关生物学过程和信号通路中的重要作用，为深入理解深县猪繁殖性状的遗传机制提供了丰富的信息，也为后续的分子育种和遗传改良工作提供了有力的理论支持。四、讨论4.1深县猪繁殖性状的遗传特征本研究对深县猪繁殖性状的表型数据进行分析，结果显示其总产仔数均值为11.56头，与前人研究中深县猪初产母猪平均窝产仔10头左右、经产母猪窝产仔12-14头的结果相近。这表明深县猪在长期的选育和进化过程中，繁殖性能相对稳定，维持在一定的水平。但总产仔数标准差为2.34，变异系数达20.24%，反映出深县猪个体间产仔数存在较大差异，这种变异可能源于遗传因素和环境因素的共同作用。深县猪保种场虽保存了10个家系，但家系间遗传背景仍有差异，不同家系的繁殖性状可能受到不同基因组合的影响，从而导致个体间产仔数的变异。环境因素如饲养管理水平、营养状况、猪舍环境等也会对产仔数产生影响。饲料中营养成分不均衡，可能影响母猪的排卵数和胚胎着床，进而导致产仔数波动。活产仔数均值为10.82头，标准差为2.21，变异系数为20.43%，变异程度与总产仔数相近。这说明除了产仔数本身的差异外，仔猪在出生时的存活情况也存在较大个体差异。部分母猪可能存在生殖疾病或免疫力低下的情况，导致胎儿在子宫内发育异常或感染疾病，影响活产仔数。初生重均值为1.25公斤，标准差为0.15，变异系数为12.00%，相对较小，说明深县猪在初生重方面遗传稳定性较好，受遗传因素的影响相对较大，环境因素对其影响相对较小。这可能是因为初生重主要由胚胎期的发育决定，而深县猪在长期的自然选择和人工选育过程中，形成了相对稳定的胚胎发育遗传机制。初生窝重均值为13.53公斤，标准差为2.76，变异系数为20.40%，与总产仔数和活产仔数的变异程度相似。这是由于初生窝重既受总产仔数影响，又受初生重影响，当总产仔数和初生重存在个体差异时，初生窝重也会随之产生较大变异。断奶仔猪数均值为9.56头，标准差为2.05，变异系数为21.44%，表明断奶时存活仔猪数量在个体间差异显著。这可能与仔猪在哺乳期的饲养管理、疾病防控以及母猪的哺育能力等因素有关。如果哺乳期猪舍卫生条件差，仔猪易感染疾病，导致死亡率增加，从而影响断奶仔猪数。断奶窝重均值为56.34公斤，标准差为10.25，变异系数为18.19%，反映出断奶时整窝仔猪总重量个体差异较大，这同样受到断奶仔猪数和仔猪在哺乳期生长速度的共同影响。死胎数均值为0.52头，标准差为0.35，变异系数为67.31%，虽然变异系数较高，但由于死胎数均值较小，实际变异范围相对有限。死胎的产生可能与母猪的繁殖疾病、感染病原体、营养不良以及遗传因素等有关。母猪感染猪细小病毒、猪繁殖与呼吸综合征病毒等，可能导致胚胎死亡，增加死胎数。某些遗传缺陷也可能导致胎儿在发育过程中死亡，形成死胎。从遗传力角度来看，猪的繁殖性状多为低遗传力性状，深县猪也不例外。低遗传力意味着表型变异受环境因素影响较大，通过传统的表型选择进行遗传改良的效果相对有限。对于深县猪繁殖性状的改良，单纯依靠表型选择，需要较长时间和较大的选育群体才能取得明显进展。因此，利用分子标记辅助选择等现代育种技术，结合GWAS分析结果，挖掘与繁殖性状相关的分子标记和基因，能够更准确地选择具有优良繁殖性状的个体，提高遗传改良效率。4.2GWAS结果的可靠性与局限性本研究采用的GWAS实验设计充分考虑了深县猪的遗传背景和繁殖性状特点，从河北正农牧业有限公司深县猪保种场选取涵盖10个家系的150头深县猪，保证了样本具有一定的遗传多样性和代表性。家系来源的多样性有助于捕捉不同遗传背景下与繁殖性状相关的遗传变异，减少因遗传背景单一导致的假阴性结果。在样本量方面，虽然150头猪对于GWAS研究来说并非大规模样本，但在深县猪种群数量有限的情况下，已尽可能涵盖了不同个体的遗传信息。有研究表明，对于中等遗传力的性状，样本量在100-200之间时，GWAS能够检测到部分中等效应的遗传位点，本研究的样本量基本满足这一要求，在一定程度上能够保证结果的可靠性。在分析方法上，选用混合线性模型（MLM）进行全基因组关联分析，该模型能够有效校正群体结构和个体间的亲缘关系对关联分析结果的影响。主成分分析（PCA）和计算亲缘关系矩阵的结果显示，深县猪群体分层现象不显著，大部分个体间的亲缘关系系数在合理范围内。这表明MLM模型在本研究中能够较好地校正潜在的混杂因素，减少假阳性和假阴性结果的出现。通过对关联分析结果绘制曼哈顿图和QQ图，评估结果的可靠性，QQ图显示大部分点分布在理论直线附近，说明关联分析结果不存在明显的系统偏差，具有较高的可信度。然而，GWAS技术本身存在一定的局限性。GWAS只能检测到与性状相关的遗传变异位点，但难以确定这些位点与性状之间的因果关系。在本研究中，虽然筛选出了与深县猪繁殖性状显著关联的SNP位点及候选基因，但这些位点和基因是否直接调控繁殖性状，还需要进一步的功能验证实验来确定。由于GWAS主要检测常见遗传变异，对于低频和罕见变异的检测能力有限。一些低频或罕见变异可能对深县猪繁殖性状具有重要影响，但在本研究中可能未被检测到。环境因素对猪繁殖性状的影响较大，而GWAS分析难以完全消除环境因素的干扰。尽管在实验设计中尽量保持饲养管理条件一致，但环境因素的细微差异仍可能对结果产生一定影响。在未来的研究中，可以结合更多的环境因素数据，采用更复杂的统计模型，进一步提高GWAS分析的准确性。4.3候选基因与繁殖性状的关联机制通过GWAS分析筛选出的候选基因，与深县猪繁殖性状之间存在着复杂的关联机制。以ESR1基因（雌激素受体1）为例，它在深县猪总产仔数性状中发挥重要作用。雌激素作为一种重要的生殖激素，与ESR1特异性结合后，形成雌激素-ESR1复合物。该复合物能够进入细胞核，与特定的DNA序列（雌激素反应元件，ERE）相互作用，从而调控下游基因的转录过程。在卵泡发育阶段，ESR1通过调节卵泡颗粒细胞的增殖和分化，影响卵泡的生长和成熟。当ESR1基因发生变异时，可能改变雌激素与受体的结合能力，进而影响卵泡的发育进程，导致排卵数发生变化，最终对总产仔数产生影响。已有研究表明，在其他猪品种中，ESR1基因的不同基因型与产仔数存在显著关联。在梅山猪群体中，ESR1基因的BB基因型母猪的总产仔数显著高于AA基因型母猪，这与本研究中ESR1基因作为总产仔数候选基因的结果相互印证，进一步支持了ESR1基因在猪繁殖性状中的重要作用。FSHB基因（促卵泡素β亚基）与深县猪总产仔数也密切相关。FSHB基因编码的促卵泡素β亚基是促卵泡素（FSH）的关键组成部分。FSH由垂体前叶分泌后，通过血液循环到达卵巢，与卵巢颗粒细胞表面的FSH受体结合，激活下游的cAMP-PKA信号通路。这一信号通路的激活能够促进颗粒细胞的增殖和分化，调节卵泡刺激素的合成和释放，从而对卵泡的发育和成熟起到关键的调控作用。当FSHB基因发生多态性变化时，可能影响FSH的生物学活性，导致卵泡发育异常，排卵数改变，进而影响总产仔数。在大白猪群体的研究中发现，FSHB基因的BB基因型母猪的排卵数显著高于AA基因型母猪，产仔数也相应增加，这与本研究中FSHB基因作为总产仔数候选基因的结果一致，表明FSHB基因在不同猪品种中对繁殖性状的调控机制具有一定的保守性。PRLR基因（催乳素受体）在深县猪活产仔数性状中具有重要意义。催乳素与PRLR结合后，激活JAK2-STAT5信号通路。在乳腺发育过程中，该信号通路促进乳腺腺泡细胞的增殖和分化，使乳腺组织发育成熟，为产后乳汁分泌做好准备。充足的乳汁供应对于仔猪的生长和存活至关重要，能够提高仔猪的存活率，从而影响活产仔数。在生殖调控方面，PRLR可能通过调节子宫内膜细胞的增殖和分化，影响胚胎着床和妊娠维持。如果PRLR基因发生突变或表达异常，可能导致乳腺发育不良，乳汁分泌不足，影响仔猪的营养摄取，降低仔猪的存活率；同时也可能影响胚胎着床和妊娠的稳定性，增加胚胎死亡的风险，进而降低活产仔数。在长白猪群体中，PRLR基因的AA基因型母猪的活产仔数显著高于BB基因型母猪，这与本研究中PRLR基因作为活产仔数候选基因的结果相呼应，说明PRLR基因在猪活产仔数性状中的调控作用在不同品种中具有一定的普遍性。GDF9基因（生长分化因子9）在深县猪活产仔数性状中发挥着关键作用。GDF9由卵母细胞分泌，通过旁分泌作用调节卵泡颗粒细胞的功能。它能够促进颗粒细胞的增殖和分化，抑制颗粒细胞的凋亡，同时调节卵泡液中各种生长因子和细胞因子的表达，为卵母细胞的成熟提供适宜的微环境。正常的卵泡发育和卵母细胞成熟是保证受孕和活产仔数的基础。当GDF9基因发生突变或多态性变化时，可能导致卵泡发育异常，卵母细胞质量下降，受精率降低，从而影响活产仔数。在一些绵羊品种中，GDF9基因的突变与产羔数的变化密切相关，这与本研究中GDF9基因作为深县猪活产仔数候选基因的结果具有相似性，表明GDF9基因在哺乳动物繁殖性状调控中具有重要的保守功能。IGF1基因（胰岛素样生长因子1）对深县猪初生重有显著影响。IGF1通过与细胞表面的IGF1受体结合，激活下游的PI3K-AKT和RAS-MAPK信号通路。在胚胎发育过程中，这些信号