北京黑猪NR6A1、VSX2、VRTN、LTBP2基因多态性与脊椎数及胴体性状的关联分析

旨在研究NR6A1、VSX2、VRTN、LTBP2基因在北京黑猪群体中基因的多态性与脊椎数和胴体性状(胴体长、胴体直长、胴体斜长及胴体重)的关联。本研究提取410头(210±7)日龄健康北京黑猪耳组织DNA,通过PCR技术和Sanger测序技术对4个基因外显子区进行基因分型,计算变异位点基因型频率以及等位基因频率的分布,并与脊椎数和胴体性状进行关联分析。结果表明,在VSX2、VRTN、LTBP2三个基因中共筛选到15个SNPs,其中包含1个CDS区的同义突变和14个3'UTR突变。使用Duncan's多重检验统计分析发现,VSX2基因中CDS区的1个突变c.536 G＞A与脊椎数关联显著(P＜0.05)。VRTN中3'UTR区共4个突变与性状关联显著(P＜0.05),其中c.2598 A＞G和c.2607 G＞T与脊椎数、胴体长、胴体直长及胴体斜长关联显著,c.3087 G＞C只与胴体斜长关联显著,c.3094 A＞G与脊椎数及胴体斜长关联显著(P＜0.05)。LTBP2上6个3'UTR突变位点(c.7158 A＞G、c.6869 C＞T、c.6319 G＞C、c.6246 G＞T、c.6170 A＞G、c.6079 G＞T)与脊椎数关联显著(P＜0.05),但与胴体性状关联不显著。综上所述,VSX2、VRTN、LTBP2基因中共有10个SNPs与脊椎数性状关联显著(P＜0.05),仅VRTN基因中4个SNPs与胴体性状关联显著(P＜0.05)。这些关联显著的SNPs可以作为北京黑猪脊椎数及胴体性状变异的候选基因功能位点。     

拜城油鸡FSHR基因多态性及其与生产性能的关联分析

【目的】探究促卵泡激素受体(follicle stimulating hormone receptor,FSHR)基因多态性与拜城油鸡产蛋性状、蛋品质及孵化性能的关系,为选育高生产性能的种鸡提供新的分子标记。【方法】以129只拜城油鸡为研究对象,利用DNA混池重测序技术筛选SNPs位点,采取Sequenom飞行时间质谱技术进行相关位点分型,通过SPSS 26.0软件的一般线性模型进行FSHR基因多态性与拜城油鸡生产性能的关联分析。【结果】FSHR基因外显子1上存在2个SNPs位点:SNP1(g.7640519 C＞A)、SNP2(g.7640639 T＞C);外显子4和5上分别存在1个SNP位点:SNP3(g.7692270 C＞T)、SNP4(g.7698478 A＞G);外显子10上存在3个SNPs位点:SNP5(g.7716725 G＞C)、SNP6(g.7715900 A＞G)、SNP7(g.7716470 T＞C),除SNP6外,检出率均在90%以上。SNP1位点仅检测到2个基因型:CC和CA,CC为优势基因型;SNP2、SNP3和SNP7位点均存在3种基因型:TT、TC和CC;SNP4位点存在3种基因型:AA、AG和GG;SNP5位点存在3种基因型:CC、CG和GG。相关分析结果表明,SNP2位点CC基因型个体开产日龄极显著高于TT基因型(P＜0.01),显著高于CT基因型(P＜0.05);TT基因型个体蛋重显著高于CT基因型(P＜0.05),入孵蛋死胚率显著高于CC基因型(P＜0.05)。SNP3位点TT基因型个体蛋壳厚度显著高于CC基因型(P＜0.05)。SNP5位点CG基因型个体300日龄总产蛋量极显著高于CC和GG基因型(P＜0.01),300日龄平均蛋重和开产蛋重均显著高于GG基因型(P＜0.05),开产体重显著高于CC基因型(P＜0.05);GG基因型个体蛋壳厚度显著低于CC和CG基因型(P＜0.05)。SNP7位点TT基因型个体开产蛋重和个体蛋重均显著高于CC和CT基因型(P＜0.05),蛋白高度和蛋黄重均显著高于CT基因型(P＜0.05),哈氏单位值极显著高于CC和CT基因型(P＜0.01)。【结论】FSHR基因SNP5(g.7716725 G＞C)位点CG基因型可作为影响拜城油鸡产蛋性能的标记基因型,SNP7(g.7716470 T＞C)位点TT基因型可作为影响拜城油鸡蛋品质的优势基因型。     

中国荷斯坦牛PIK3CB基因多态性及其与繁殖和产奶性状的关联分析

【目的】 探究磷脂酰肌醇-4,5-二磷酸3-激酶催化亚单位β（PIK3CB）基因多态性及其与中国荷斯坦牛繁殖和产奶性状的关系。【方法】 通过混池测序对中国荷斯坦牛PIK3CB基因进行单核苷酸多态性（SNP）位点筛选,采用竞争性等位基因特异性PCR （KASP）技术在1 160头健康泌乳中国荷斯坦牛中进行SNP分型并进行群体遗传学分析,采用线性模型进行SNP与11个繁殖和产奶性状基于单位点和单倍型组合的关联分析。【结果】 在PIK3CB基因中共检测到了17个SNPs,筛选出7个SNPs用于后续分析。关联分析发现,7个SNPs与多个目标性状存在显著或极显著的关联（P<0.05;P<0.01）;位于外显子区域的g.130433743 A>G位点AA基因型个体和位于可变剪接区域的g.130448069 G>A位点GG基因型个体,其经产牛首末次配种间隔、产奶量、乳蛋白量和乳脂量最低,体细胞评分最高,上述基因型个体具有较短的首末次配种间隔,而产奶性能相对较差;g.130387717 G>A位点AA基因型个体,其初配日龄和青年牛首末次配种间隔最低,产奶量、乳蛋白量和乳脂量最高,该基因型个体的繁殖和产奶性能均较好,上述3个SNPs位点可作为中国荷斯坦牛繁殖和产奶性状的候选位点重点关注。单倍型分析发现,PIK3CB基因的g.130387717 G>A、g.130430832 A>－、g.130433743 A>G、g.130433982 C>T、g.130446073 C>T和g.130448069 G>A 6个SNPs紧密连锁形成一个单倍型块,且与多个目标性状存在显著或极显著关联（P<0.05;P<0.01）,其中H2H3和H2H4单倍型组合个体的繁殖和产奶性能较好,为优势单倍型组合。【结论】 中国荷斯坦牛PIK3CB基因存在丰富的遗传变异,其多态性与繁殖和产奶性状存在关联,g.130433743 A>G、g.130448069 G>A和g.130387717 G>A位点可作为潜在分子标记,为中国荷斯坦牛的平衡育种提供理论依据。     

藏羊BMPRⅡ基因多态性及其单倍型与产羔性状相关分析

【目的】 研究骨形态发生蛋白受体Ⅱ(bone morphogenetic protein receptor Ⅱ, BMPRⅡ)基因多态性及其单倍型与藏羊产羔性状的相关性。【方法】 以433只藏羊母羊为研究对象, 采用改良多重高温连接酶检测反应技术(improved multiple ligase detection reaction, iMLDR)对BMPRⅡ基因9个单核苷酸多态性(single nucleotide polymorphism, SNP)位点在藏羊群体中的多态性进行检测, 使用Haploview 4.2软件分析其连锁不平衡性并构建单倍型, 应用基因关联分析探究BMPRⅡ基因多态性及单倍型与藏羊产羔性状的关联。【结果】 藏羊BMPRⅡ基因外显子12中存在6个SNPs, 分别为: g.90192 T>C、g.142532 C>T、g.142614 A>G、g.142751 T>C、g.143138 A>G和g.143189 G>A, 外显子3、9和13中各存在1个SNP, 分别为: g.143570 C>G、g.124843 A>C和g.145233 A>G, 这9个SNPs均为同义突变。9个SNPs中, 除g.90192 T>C外, 均存在3种基因型。多态信息含量分析显示, 群体在g.90192 T>C、g.142614 A>G和g.145233 A>G位点处于低度多态(PIC<0.25), 在g.124843 A>C、g.142532 C>T、g.142751 T>C、g.143138 A>G、g.143189 G>A和g.143570 C>G位点均处于中度多态(0.25<PIC<0.5)。χ²适合性检验结果显示, 所有突变位点均未偏离哈代-温伯格平衡状态。相关性分析显示, 不同位点不同基因型与藏羊产羔性状均无显著相关(P>0.05), 但g.142532 C>T位点CT基因型平均产羔数高于CC和TT基因型, g.142751 T>C位点CC基因型平均产羔数高于TT和TC基因型, g.143189 G>A位点GA基因型平均产羔数高于GG和AA基因型, g.143570 C>G位点CG基因型平均产羔数高于CC和GG基因型, g.145233 A>G位点GG基因型平均产羔数高于AA和AG基因型, 差异均不显著(P>0.05)。BMPRⅡ基因中除g.90192 T>C位点外的8个SNPs位点在藏羊群体中共形成14种单倍型(H1~H14), 单倍型与藏羊产羔数间均无显著相关(P>0.05)。【结论】 BMPRⅡ基因g.142532 C>T、g.142751 T>C、g.143189 G>A、g.143570 C>G和g.145233 A>G位点对藏羊产羔数有一定潜在的影响。     

中国荷斯坦奶牛RPL23A、ACACB基因3个SNPs检测及其与产奶性状的关联分析

旨在开展奶牛群体高乳成分功能基因的验证与筛选,利用微流控芯片自主选育技术分析北京地区母牛群体高乳蛋白、高乳脂基因频率,同时分析基因多态性及其与产奶性状的相关性。本研究对北京地区8个大型奶牛场1 596头中国荷斯坦奶牛RPL23A、ACACB基因的多态性进行了检测,所有个体均为3胎以内的泌乳牛,收集每头牛所有测定日的产奶量、乳脂率、乳蛋白率等数据,同时对3个多态位点不同基因型与产奶性状进行了关联分析。RPL23A基因的SNP位点g.20146771G>A,在第一泌乳期,与产奶量、乳脂量、乳脂率和乳蛋白量均达到极显著关联（P<0.01）,第二泌乳期,g.20146771G>A与5个产奶性状均呈极显著关联（P<0.01）。ACACB基因的SNP位点g.63962768C>T,在第一泌乳期,与产奶量、乳脂量、乳脂率和乳蛋白量均呈极显著关联（P<0.01）,在第二泌乳期,与产奶量、乳脂量、乳脂率和乳蛋白率呈极显著关联（P<0.01）。结果表明,RPL23A、ACACB基因可以作为影响中国荷斯坦奶牛产奶性状的候选基因用于标记辅助选择,以上基因位点可能通过直接或间接的途径影响奶牛的乳脂或乳蛋白性状,对产奶性状起到重要调控作用。本研究为荷斯坦奶牛后续的标记辅助选择奠定了良好基础。     

黔东南小香鸡MEF2C基因外显子多态性及其与体重、体尺性状的关联分析

【目的】探究肌细胞增强因子2C（myocyte enhancer factor 2C,MEF2C）基因多态性对黔东南小香鸡体重和体尺性状的影响。【方法】试验以251羽300日龄黔东南小香鸡为研究对象,利用PCR扩增和直接测序技术对MEF2C基因全部（15个）外显子进行单核苷酸多态性（single nucleotide polymorphism,SNP）位点筛选,并运用SPSS 18.0统计分析软件对不同SNPs位点所对应的基因型与黔东南小香鸡的体重、体尺性状进行相关性分析。【结果】黔东南小香鸡群体中仅在MEF2C基因外显子12上发现3个SNPs位点：g.1819 T＞C、g.2426 T＞C和g.2543 C＞T,在其他外显子中未发现多态位点。χ²检验结果表明,3个SNPs在黔东南小香鸡群体中均偏离Hardy-Weinberg平衡状态（P＞0.05）。关联分析结果显示,g.1819 T＞C位点TC基因型个体体斜长、胸深和胫围均极显著高于TT和CC基因型,体重极显著高于TT基因型（P＜0.01）;TC基因型个体胫长显著高于TT和CC基因,龙骨长显著高于TT基因型（P＜0.05）。g.2426 T＞C和g.2543 C＞T位点互为连锁平衡位点,对胫围有极显著影响（P＜0.01）。【结论】黔东南小香鸡MEF2C基因外显子具有较强的保守性,g.1819 T>C位点可以作为黔东南小香鸡分子标记辅助选育参考位点。     

梅花鹿成纤维细胞因子受体2基因多态性及其与茸重性状的关联分析

【目的】探索梅花鹿成纤维细胞因子受体2(fibroblast growth factor receptor 2,FGFR2)基因多态性及其对茸重性状的影响。【方法】应用直接测序法对梅花鹿FGFR2基因的全部外显子进行测序分析,通过MassARRAY^® SNP分型技术对314头24月龄梅花鹿进行基因分型和单倍型分析,分析FGFR2基因不同基因型和单倍型与茸重的关联性。【结果】在梅花鹿FGFR2基因中共发现12个多态性位点,其中5个位于外显子区域,且突变均未引起氨基酸改变,属于同义突变,其余7个位点均存在于内含子区域。分型结果显示,g.80975864 T>G位点未分型成功,后续对其余11个位点进行了分析,g.80943673 T>C、g.80943683 C>A及g.80938352 C>T 3个位点属于中度多态位点(0.25＜P＜0.5),其余位点均属于低度多态位点(P＜0.25)。χ²检验结果表明,g.80998742 G>A和g.80987708 G>A 2个位点偏离Hardy-Weinberg平衡(P＜0.05),其他9个位点均处于Hardy-Weinberg平衡(P＞0.05)。关联分析结果表明,11个多态性位点各基因型之间的茸重差异均不显著(P＞0.05)。单倍型结果显示,FGFR2基因存在5种单倍型,不同单倍型间梅花鹿茸重差异均不显著(P＞0.05)。【结论】FGFR2基因的11个突变位点可能不是影响梅花鹿茸重性状的关键位点。     

梅花鹿MSRB3基因多态性及其与茸重性状的关联分析

【目的】 本研究旨在探索梅花鹿甲硫氨酸亚砜还原酶B3（methionine sulfoxide reductase B3，MSRB3）基因多态性及其与茸重性状的关联性。【方法】 选取高、低产茸量的梅花鹿各8只，应用DNA直接测序法检测基因变异位点。采用MassARRAY^® SNP分型技术对314头24月龄梅花鹿MSRB3基因进行基因分型，并结合梅花鹿茸重性状数据进行了关联分析和单倍型分析。【结果】 在梅花鹿MSRB3基因中共发现6个SNPs位点，其中2个SNPs位点位于外显子区域，且突变均未引起氨基酸改变，属于同义突变，其余4个SNPs位点均存在于内含子区域。分型结果显示4个样本未分型成功，其余310个样本进行后续分析。各位点观测杂合度和期望杂合度基本一致，g.44455582 T>C、g.44455759 C>T、g.44414424 T>C、g.44350306 T>C及g.44340836 G>A等5个位点的杂合度较高，均属于中度多态（0.25<PIC<0.5），g.44340734 G>C位点杂合度较低，属于低度多态（PIC<0.25）。卡方检验结果表明，g.44455759 C>T位点偏离Hardy-Weinberg平衡状态（P<0.05），其他5个位点均处于Hardy-Weinberg平衡状态（P>0.05）。对6个SNPs位点的不同基因型与梅花鹿茸重的关联分析结果表明，g.44455582 T>C位点的TT基因型个体茸重极显著高于CC和TC基因型（P<0.01）。单倍型分析结果显示，g.44340734 G>C和g.44350306 T>C、g.44455582 T>C和g.44455759 C>T位点存在强连锁，各产生3种单倍型，在Block 2区块中单倍型TC茸重显著高于单倍型CT （P<0.05）。【结论】 MSRB3基因与梅花鹿茸重性状密切相关，g.44455582 T>C位点和单倍型TC可作为筛选高产梅花鹿的分子标记。     

湖羊PLAG1基因5'调控区多态性及其与早期体重的关联分析

旨在克隆湖羊PLAG1基因5'调控区序列,明确PLAG1基因5'调控区多态性与湖羊早期体重的关系,寻找用于湖羊生长性状辅助选择的分子标记。本研究利用5'RACE技术鉴定PLAG1基因转录起始位点,以456只断奶湖羊为对象,利用测序法筛选PLAG1 5'调控区SNP位点,使用SPSS 18.0软件分析不同基因型对湖羊初生体重和断奶体重的影响。结果表明,PLAG1基因转录起始位点均位于g.30535位点,测序发现,PLAG1基因5'调控序列存在15个SNPs位点,在湖羊群体中均有3种基因型,均处于Hardy-Weinberg平衡状态。关联分析发现,g.28888A>T位点AA型、g.28901G>T位点GG型、g.30802C>T位点CC型、g.30825A>G位点AA型的初生体重均显著大于相应的杂合基因型（P<0.05）,g.30737T>C位点CC型断奶体重显著小于TC型（P<0.05）。连锁分析发现,g.28888A>T、g.28901G>T、g.28913T>A、g.29717G>C、g.29724A>C、g.29725G>A、g.29768T>C、g.29771A>G、g.30141T>A、g.30144A>G、g.30691C>A、g.30694A>T、g.30737T>C、g.30802C>T、g.30825A>G,15个位点连锁且有18种单倍型。将样本量较多的5种单倍型与湖羊早期体重进行关联分析发现,AAGGTAGGAAGGTTAATTAACCAATCCCAA单倍型初生体重显著大于ATGTTAGCACGATCAGTAAGCAATTCCTAG单倍型（P<0.05）,断奶体重显著大于AAGGAAGGAAGGTTAATTAACCAACCCCAA单倍型（P<0.05）。上述结果表明,湖羊群体中PLAG1基因5'调控区存在15个连锁的SNPs位点,其中g.28888A>T、g.28901G>T、g.30802C>T、g.30825A>G位点与初生体重显著相关（P<0.05）,g.30737T>C位点与断奶体重显著相关（P<0.05）,AAGGTAGGAAGGTTAATTAACCAATCCCAA单倍型与初生体重和断奶体重显著相关（P<0.05）,说明PLAG1基因可作为湖羊生长性状选择的候选遗传标记。     

F3代波杂山羊VNN1基因多态性及其与产羔数的关联分析

【目的】 研究泛酰巯基乙胺酶-1(panthenoyl mercaptoethamine-1,VNN1)基因多态性与F3代波杂山羊产羔数的关系,从而从分子水平指导F3代波杂山羊的育种工作。【方法】 选取107只F3代波杂山羊,采集血样提取DNA,根据VNN1基因(登录号:NC_030816.1)外显子序列设计7对引物,进行PCR扩增及测序,运用DNAStar软件分析测序结果,对可能存在的单核苷酸多态性(SNP)位点外显子全群进行PCR扩增测序,并进行Hardy-Weinberg平衡状态及遗传多样性分析;应用SPSS 16.0统计学软件对VNN1基因不同基因型与F3代波杂山羊产羔数进行关联分析。【结果】 在VNN1基因在第5和7外显子中共发现7个SNPs位点:g.10956 A→G、g.11010 C→T、g.11103 C→T、g.11284 C→A、g.11313 T→C、g.18094 C→T和g.18158 G→C,且每个位点均存在3种不同基因型。χ²检验显示,g.10956 A→G、g.11010 C→T、g.11103 C→T、g.11313 T→C和g.18094 C→T位点均处于Hardy-Weinberg平衡状态,g.11284 C→A和g.18158 G→C位点均处于Hardy-Weinberg不平衡状态。多态信息含量均为中度多态((0.25<PIC<0.5)。关联分析发现,g.11010 C→T和g.11313 T→C位点与F3代波杂山羊产羔数显著相关,其中g.11010 C→T位点CC基因型产羔数显著低于CT和TT基因型,g.11313 T→C位点TT基因型产羔数显著低于TC和CC基因型(P<0.05);g.10956 A→G、g.11103 C→T和g.18158 G→C位点对F3代波杂山羊产羔数均无显著影响(P>0.05)。【结论】 VNN1基因在第5和7外显子中共发现7个SNPs位点,其中g.11010 C→T和g.11313 T→C位点对F3代波杂山羊产羔数有显著影响,VNN1基因可作为F3代波杂山羊的候选基因。     

绵羊FUT8基因多态性及其与生长性状的关联分析

【目的】研究岩藻糖基转移酶8(fucosyltransferase 8,FUT8)基因g.74522417 C＞T位点多态性及其对绵羊生长性状的影响,为绵羊分子育种筛选新的遗传标记。【方法】采用Sequenom MassARRAY? SNP分型技术、Illumina OvineSNP 600K BeadChip及Illumina OvineSNP 50K BeadChip芯片检测技术对湖羊(n=992)和乌珠穆沁羊群体(n=333)FUT8基因g.74522417 C＞T位点进行分型,并将其多态性与湖羊和乌珠穆沁羊的生长性状进行关联分析。【结果】FUT8基因g.74522417 C＞T位点在乌珠穆沁羊和湖羊群体中均检测出3种基因型:CC、TC和TT。群体遗传分析表明,g.74522417 C＞T位点在乌珠穆沁羊和湖羊群体中均处于中度多态(0.25＜PIC＜0.50)。卡方适应性检验结果表明,g.74522417 C＞T位点在乌珠穆沁羊和湖羊群体中均处于哈代-温伯格平衡状态(P＞0.05)。关联分析结果表明,FUT8基因g.74522417 C＞T位点与乌珠穆沁羊4月龄体长和6月龄体高均呈显著相关(P＜0.05);与湖羊3月龄管围、5月龄腰角宽和管围、6月龄眼肌面积和管围均呈显著或极显著相关(P＜0.05;P＜0.01)。【结论】FUT8基因g.74522417 C＞T位点可作为绵羊生长发育潜在分子标记,为FUT8基因在绵羊分子育种中的应用提供参考。     

摩拉水牛PRKAA2基因多态性与生长性状的关联分析

试验旨在探讨PRKAA2基因在摩拉水牛群体中的遗传多态性及其与生长性状的关联性,进而得到显著相关的遗传标记,为摩拉水牛的分子标记辅助选择提供理论依据。利用DNA测序法等技术进行候选基因PRKAA2外显子4及部分内含子3单核苷酸多态性（SNPs）检测,并采用生物信息学方法结合统计软件分析PRKAA2基因与摩拉水牛生长性状的关联性。结果显示,摩拉水牛PRKAA2基因中共检测到4个SNPs位点：c.462 G>A、IVS3.557 T>C、IVS3.560 C>T和IVS3.565 G>A,均包含3种基因型,且符合Hardy-Weinberg平衡。4个SNPs位点在水牛群体中均为中度多态,可以构建3种单倍型,T-C-G-G为优势单倍型,其中IVS3.560 C>T与c.462 G>A位点之间存在完全连锁不平衡,关联分析结果显示,单倍型H2与摩拉水牛体斜长和腰角宽呈显著相关（P<0.05）。c.462 G>A位点与摩拉水牛体高、十字部高、尻长、坐骨端宽和腰角宽呈显著相关,其中GG和AA基因型个体的体高和十字部高均显著高于GA基因型个体,GG基因型个体的尻长显著高于GA和AA基因型个体,AA基因型个体的坐骨端宽和腰角宽显著高于GA基因型个体（P<0.05）;IVS3.557 T>C位点与摩拉水牛的生长性状指标未达到显著相关（P>0.05）;IVS3.560 C>T位点与摩拉水牛体高、十字部高和尻长呈显著相关,其中CC和TT基因型个体的体高和十字部高均显著高于CT基因型个体,CC基因型个体的尻长显著高于CT基因型个体（P<0.05）;IVS3.565 G>A位点与体斜长呈显著相关,其中GG和GA基因型个体的体斜长显著高于AA基因型个体（P<0.05）。综上,PRKAA2基因的c.462 G>A、IVS3.560 C>T和IVS3.565 G>A位点对摩拉水牛部分生长性状均有显著影响,可作为摩拉水牛品种早期选育的候选基因和分子辅助标记。     

灌服黄芪多糖对奶牛产奶量、免疫性能及抗氧化性能的影响

【目的】研究黄芪多糖(APS)对奶牛产奶量、免疫功能及抗氧化性能的影响。【方法】选取体况良好,体重、产奶量、泌乳天数、胎次接近的中国荷斯坦奶牛44头,随机分为对照组(CK)及APS20、APS55和APS90组,每组11头。所有试验牛均饲喂基础饲粮。对照组灌服温水200 mL,试验组分别灌服20、55、和90 g/d的APS溶液200 mL,预试期10 d,正试期60 d。正试期内,每10 d统计1次产奶量;于正试期的第0、30、60天采集乳样,测定乳中体细胞数(SCC);于正试期第0、30、60天晨饲前,采集奶牛尾根静脉血样,测定血清中相关免疫指标、抗氧化指标。【结果】在试验全期,各组产奶量均差异不显著(P＞0.05)。在体细胞数方面,与CK组相比,第30和60天时,APS55组体细胞评分(SCS)均显著降低(P＜0.05),APS20和APS90组SCS在各时间点均差异不显著(P＞0.05),且APS90组极显著高于APS55组(P＜0.01)。在免疫指标方面,与CK组相比,第30天时,APS20组血清中白细胞介素2(IL-2)含量显著增加(P＜0.05),APS55与APS90组血清中IL-2含量极显著增加(P＜0.01);第60天时,APS20、APS55和APS90组血清中IL-2含量均极显著增加(P＜0.01),且APS55与APS90组IL-2含量极显著高于APS20组(P＜0.01);第60天时,APS55组血清中IL-6含量显著增加(P＜0.05)。在抗氧化指标方面,与CK组相比,第30、60天时,APS20、APS55和APS90组血清中过氧化氢酶(CAT)活性均极显著增加(P＜0.01);第60天时,ASP90组血清中丙二醛(MDA)浓度显著降低(P＜0.05),APS20、APS55和APS90组血清中谷胱甘肽过氧化物酶(GSH-Px)活性均显著增加(P＜0.05),APS55和APS90组血清中总抗氧化能力(T-AOC)极显著增加(P＜0.01),且显著高于APS20组(P＜0.05)。【结论】55 g/d APS可以提高奶牛免疫功能、抗氧化能力,并且降低乳中体细胞数,结果为APS在奶牛生产应用提供参考。     

北京黑猪DKK3和CCR1基因多态性检测及其与背膘厚的关联分析

旨在探讨DKK3和CCR1基因多态性及其对北京黑猪背膘厚的影响，以期能够筛选出可用于提高北京黑猪背膘厚性状的分子标记，为北京黑猪优良性状定向选育提供依据。本研究以413头健康的北京黑猪为试验材料，对其进行DNA提取、PCR扩增、基因测序筛选候选基因的SNPs，同时利用SAS 9.2软件对突变位点的不同基因型分别与北京黑猪的六七肋背膘厚和腰荐部背膘厚进行关联分析。结果显示，在北京黑猪DKK3基因的外显子区共检测到10个SNPs，包含1个可变剪切突变、1个错义突变和8个3'UTR突变，其中有4个3'UTR突变位点均与腰荐部背膘厚显著关联(P<0.05)，分别为g.47443779 T>C、g.47444258 C>T、g.47444912 G>T和g.47445304 T>C，其余位点与六七肋背膘厚和腰荐部背膘厚均不显著关联(P>0.05)。在CCR1基因的外显子区共检测到3个SNPs，包含2个同义突变和1个错义突变，其中1个同义突变位点g.29229037 G>A与六七肋背膘厚显著关联(P<0.05)。对DKK3的10个SNPs进行连锁不平衡分析，结果显示，g.47443779 T>C和g.47443783 C>T处于完全连锁状态(D’=1)，g.47443783 C>T和g.47443858 C>T处于完全连锁状态(D’=1)，g.47444258 C>T和g.47444275 C>T呈现高度连锁状态(D’=0.98)。综上所述，DKK3和CCR1基因的SNPs分别与北京黑猪的腰荐部背膘厚和六七肋背膘厚有显著关联性，可为北京黑猪背膘厚的早期分子选育提供参考。     

MET基因多态性与中国荷斯坦牛繁殖和泌乳性状的关联分析

为筛选与中国荷斯坦牛繁殖和泌乳性状相关的遗传标记，本研究从分子水平探究了MET基因单核苷酸多态（single nucleotide polymorphisms，SNPs）与繁殖和泌乳性状的相关性。通过混池测序法在70头无血缘关系的健康中国荷斯坦公牛中扫描目标基因的SNPs，采用竞争性等位基因特异性PCR （kompetitive allele-specific PCR，KASP）对其后代1 160头中国荷斯坦泌乳牛进行部分SNP的基因分型，基于连锁不平衡分析获得双倍型信息，并采用线性模型对SNP及其双倍型与5个繁殖和6个泌乳性状的估计育种值（estimated breeding values，EBV）进行关联分析。随后，针对位于功能区域的SNPs进行生物信息学分析和基因表达量相关性分析，初步探究其潜在的调控机制。本研究在MET基因中共检测到19个SNPs，说明该基因存在丰富的遗传变异；就其中7个SNPs进行了基因分型，分型结果与目标性状的关联分析结果显示，7个SNPs与多个性状存在显著（P<0.05）或极显著（P<0.01）的关联；位于上游调控区的g.51737889T>C位点GG基因型个体的初产日龄、青年牛首末配种间隔、经产牛首末配种间隔、体细胞评分的EBV最低，说明该基因型个体同时具有较好的繁殖性能和泌乳性能；另一个位于上游调控区的g.51736640A>C位点GG基因型个体的初产日龄、青年牛首末配种间隔、经产牛首末配种间隔、体细胞评分的EBV最低，但乳蛋白率、乳脂率的EBV最低，说明该基因型个体繁殖性能较好，但泌乳性能相对较差；此外，位于第6外显子的g.51660569G>A位点TC基因型个体初产日龄的EBV较低，青年牛首末配种间隔、经产牛首末配种间隔、体细胞评分的EBV最低，说明该基因型个体具有较好的繁殖性能。连锁不平衡分析发现，该基因的7个SNPs形成2个单倍型块，其双倍型与目标性状的关联分析进一步表明，单倍型块1中H1H3个体的繁殖性能较好且泌乳性能中等，单倍型块2中H4H4个体的繁殖性能和泌乳性能均较优秀，为两个优势双倍型；此外，H1H3包含g.51660569G>A位点的TC基因型，H4H4为g.51737889T>C、g.51736640A>C两位点GG基因型的组合，其关联分析结果与单个位点的结果一致。基序分析结果表明，g.51737889T>C与g.51736640A>C位点的优势等位基因G可富集到与基因激活、细胞增殖、分化相关的转录因子，且两位点GG基因型的MET基因表达量均最高，其双倍型H4H4表达量也较高，进一步验证了关联分析结果。综上，本研究获得了中国荷斯坦牛MET基因的多态性图谱，并通过关联分析、生物信息学预测及基因表达量分析发现并验证了MET基因与繁殖、泌乳性状的遗传关联，为中国荷斯坦牛高产高效选育提供了可用的遗传标记。     

牦牛EPO、PPARα基因SNP与高原低氧适应的相关性分析

试验旨在探究促红细胞生成素基因（EPO）、过氧化物酶体增殖剂激活受体α（PPARα）的遗传多样性,并分析其多态性与牦牛高原低氧适应性的相关性。采集不同海拔高度的6个牦牛类群（中甸牦牛、麦洼牦牛、斯布牦牛、类乌齐牦牛、帕里牦牛、申扎牦牛）以及三江黄牛共375头耳样,提取DNA并分别构建DNA池,采用直接测序法结合PCR-RFLP检测分析EPO、PPARα基因的多态性,最后应用SHEsis软件统计分析候选基因SNPs与牦牛高原适应性的相关性。结果表明,EPO基因存在3个SNPs位点：rs527G→A、rs1031A→T、rs1192T→C;PPARα基因存在3个SNPs位点：rs77363C→T、rs77471C→A和rs77534C→T。χ²适合性检验结果显示,EPO基因3个SNPs位点均符合Hardy-Weinberg平衡状态;PPARα基因的rs77363C→T位点上,6个牦牛类群都处于平衡状态（P>0.05）,在PPARα基因的rs77471C→A和rs77534C→T位点,麦洼牦牛偏离Hardy-Weinberg平衡（P<0.05）。单倍型分析得出,EPO基因的ATC单倍型在高海拔地区牦牛中的分布频率随海拔升高而升高;PPARα基因的TAC单倍型在6个牦牛类群中分布频率显著高于其他单倍型。研究表明,EPO、PPARα基因可作为牦牛适应高原环境的分子标记,为进一步探讨牦牛高原低氧适应性机制提供一定理论帮助。