乌蒙凤鸡FGF23基因多态位点筛查与生物信息学分析

研究旨在挖掘成纤维生长因子23(fibroblast growth factors,FGF23)基因的单核苷酸多态(SNP)位点,为进一步研究该基因遗传变异奠定基础。试验以乌蒙凤鸡为研究对象,采用PCR产物直接测序法对FGF23基因进行SNP位点筛查,并进行遗传特性、连锁不平衡、单倍型、双倍型和生物信息学分析。结果显示,仅在乌蒙凤鸡FGF23基因启动子区域检测到3个中度多态的SNPs位点,分别为:g.73424341 CA、g.73424417 AG和g.73424701 TA,每个SNP位点均产生3种基因型。χ~2检验结果发现,仅g.73424417 AG突变位点的基因型分布极显著偏离Hardy-Weinberg平衡(P0.01)。连锁不平衡、单倍型和双倍型分析结果显示,3个SNPs位点中仅g.73424417 AG和g.73424701 TA位点间存在强连锁不平衡,共发现4种单倍型和8种双倍型,双倍型H1H2频率最高,其次为H3H4,频率最低的为H2H2。生物信息学分析发现,FGF23基因的核心启动子区域很可能在-400～-300 bp处,本研究发现的3个SNPs位点不在核心启动子区;突变前g.73423055～g.73425055区域总转录因子数为293个,突变后为300个。g.73424341 CA位点突变前后转录因子不变,均为ICSBP;g.73424417 AG位点突变使转录因子由GR转变为C/EBPalp;g.73424701 TA位点突变前后均没有转录因子。推测SNP位点对调控启动子功能元件可能存在重要影响。     

湖羊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基因可作为湖羊生长性状选择的候选遗传标记。     

黔北麻羊MC4R基因多态性及其与生长性状的关联分析

【目的】探究黔北麻羊黑素皮质素受体-4（melanocortin receptor-4,MC4R）基因多态性及其与生长性状的关联性,为黔北麻羊的选种选育提供理论依据。【方法】选取196只6月龄黔北麻羊,利用DNA混合池结合Sanger测序筛选MC4R基因单核苷酸多态性（SNP）位点,并采用SPSS 22.0软件中的一般线性模型（GLM）对MC4R基因SNP位点与黔北麻羊生长性状进行关联分析。【结果】在黔北麻羊MC4R基因中发现4个SNPs位点：g.59345469 C>A和g.59346773 T>C位点分别位于5'-和3'-端非编码区（UTR）;g.59345871 A>C位点突变导致第37位天冬氨酸（Asp）变为丙氨酸（Ala）,引起RNA二级结构及蛋白质二级结构、三级结构发生明显改变;g.59346263 A>G位点突变导致第168位异亮氨酸（Ile）变为缬氨酸（Val）,对RNA二级结构和蛋白质二级结构有明显影响。关联分析结果显示,g.59345469 C>A位点对黔北麻羊体重、胸围和管围均有显著影响（P<0.05）;g.59345871 A>C位点对黔北麻羊体重、体高、管围和胸围均有显著影响（P<0.05）;g.59346263 A>G位点对黔北麻羊体高、体重和胸围均有显著影响（P<0.05）;g.59346773 T>C位点对黔北麻羊体重和胸围均有显著影响（P<0.05）。4个SNPs位点联合共检测到9种单倍型和21种双倍型,其联合对黔北麻羊体重、体高、胸围和管围均产生显著遗传效应（P<0.05）,纯合双倍型H5H5（AAAAGGCC）个体的生长性状优于其他双倍型。【结论】黔北麻羊MC4R基因存在4个SNPs位点,与黔北麻羊体重和胸围均存在显著关联,可作为黔北麻羊体重和胸围的遗传标记用于分子育种。     

贵州地方山羊BMPR-ⅠB基因启动子区多态性及生物信息学分析

为探究骨形态发生蛋白受体ⅠB（bone morphogenetic protein receptor ⅠB，BMPR-ⅠB）基因启动子区多态性及其与繁殖性状的相关性，本试验以贵州3个地方山羊品种（黔北麻羊、贵州黑山羊及贵州白山羊）为研究对象，通过构建混合DNA池PCR产物直接测序技术筛选SNPs位点，并采用多种生物信息学软件预测SNPs位点对核心启动子区、CpG岛和转录因子结合位点的影响。结果显示，BMPR-ⅠB基因启动子区存在4个SNPs位点：g.29893005 T > C、g.29893016 C > T、g.29893729 C > G和g.29894370 C > T。生物信息学软件预测得到BMPR-ⅠB基因核心启动子区和CpG岛，SNPs位点导致转录因子结合位点发生改变，g.29893005 T > C和g.29893016 C > T突变使原来的转录因子结合位点Sp1消失；g.29893729 C > G突变使原来转录因子结合位点AP-1消失，从而产生新的转录因子结合位点USF；g.29894370 C > T突变产生新的转录因子结合位点C/EBPalp，使原来的转录因子结合位点Sp1改变为C/EBPalp和Sp1。由此推测，SNPs位点对调控启动子功能元件可能存在重要影响。     

关岭牛TBC1D7基因多态性及其与生长性状的关联分析

为探究TBC1蛋白家族第7成员（TBC1 domain family member 7，TBC1D7）基因单核苷酸多态性（SNP）对关岭牛生长性状的影响，试验选取116头贵州关岭牛，测定其8项生长性状指标，提取血液基因组DNA，并利用DNA混池扩增TBC1D7基因全部外显子，通过直接测序法筛查SNP。利用不同独立样本DNA扩增SNP位点对应外显子并测序，以验证位点并进行基因分型。利用SPSS 19.0软件单因素方差分析TBC1D7基因不同基因型和关岭牛生长性状之间的相关性。结果显示，在TBC1D7基因第5、6和8外显子上共发现5个SNPs位点：g.12972 T>C、g.12984 A>G、g.20538 C>T、g.20577 T>C和g.24807 G>A，均存在3种基因型；χ²检验结果显示，5个SNPs位点均显著偏离Hardy-Weinberg平衡状态（P<0.05）；群体遗传参数分析发现，5个SNPs位点杂合度较高，多态信息含量均表现为中度多态（0.25<PIC<0.5）。关联性分析结果发现，关岭牛TBC1D7基因g.24807 G>A位点GA基因型个体体斜长、胸围和后臀围均显著高于AA基因型个体（P<0.05）。在本试验群体中，5个SNPs存在6种单倍型（H1~H6）和9种双倍型（H1H1、H2H2、H2H3、H2H6、H3H3、H3H4、H3H5、H3H6和H5H5），其中H3和H2为优势单倍型，H3H3、H2H2和H2H3为优势双倍型；双倍型关联分析结果显示，H3H6和H5H5个体体斜长和后臀围相比其他双倍型个体占有明显优势，可能是有利基因型。结果表明，TBC1D7基因g.24807 G>A位点为影响关岭牛生长发育的重要多态位点，本研究为筛选有助于关岭牛生长发育的分子标记提供了理论参考。     

长顺绿壳蛋鸡SFRP5基因外显子1的SNPs鉴定及其与常规肌肉品质的关联性分析

分泌型卷曲相关蛋白-5(SFRP5)对畜禽脂质代谢起着至关重要的调控作用。为了探讨鸡SFRP5基因SNPs与肌肉品质的相关性,采用PCR产物直接测序法检测SFRP5基因外显子1的SNPs位点,分析其对长顺绿壳蛋鸡胸肌肌肉品质的影响。结果显示:在SFRP5基因外显子1的CDS区中发现3个SNP沉默突变位点:g.23251706 C>T为低度多态位点,g.23251882 C>T和g.23252056 C>T为中度多态位点;3个位点的等位基因C和基因型CC均处于优势地位,且基因型分布均未偏离Hardy-Weinberg平衡状态。连锁不平衡分析表明:3个SNP位点之间不存在强连锁不平衡,发现4种单倍型和10种双倍型,单倍型H1和双倍型H1H1频率最高。关联性分析结果表明:g.23251706C>T位点对胸肌蒸煮损失的影响达到显著水平(P<0.05),TT基因型有利于降低蒸煮损失;g.23252056C>T位点则对失水率的影响达到显著水平(P<0.05),TT基因型有利于提高肌肉保水力;3个SNP位点联合组成的双倍型对胸肌的5个常规肉品质指标均有不同程度的影响。综合所有指标在双倍型之间的差异,双倍型H2H2是改善长顺绿壳蛋鸡肌肉品质最有利的双倍型,揭示在SFRP5基因外显子1中发现的3个SNP位点相互联合能显著影响肉质,可作为鸡肌肉品质选择的遗传标记。     

黔北麻羊UCP2基因组织表达及其多态性与生长性状的关联分析

本研究旨在探讨黔北麻羊解耦联蛋白2（uncoupling protein 2,UCP2）基因的遗传多态性及其与生长性状的关联性,以期为黔北麻羊的育种改良提供理论依据。以200头6月龄黔北麻羊为研究对象,采用直接测序法检测黔北麻羊UCP2基因的SNP位点,并将其与黔北麻羊生长性状进行关联分析;实时荧光定量PCR检测UCP2基因在黔北麻羊不同组织中的表达水平。结果显示,UCP2基因存在4个SNPs位点,分别为：g.29614172 A>G、g.29619320 G>T、g.29619721 C>T和g.29620037 A>G,均产生3种基因型,其中g.29619320 G>T和g.29620037 A>G为错义突变,g.29619320 G>T突变导致甘氨酸（Gly）变为半胱氨酸（Cys）,g.29620037 A>G突变导致苏氨酸（Thr）突变为丙氨酸（Ala）。多态信息含量显示,4个SNPs位点均表现为中度多态（0.25 < PIC < 0.5）。χ²检验结果表明,g.29614172 A>G、g.29619721 C>T和g.29620037 A>G突变位点在黔北麻羊群体中均处于Hardy-Weinberg平衡状态（P>0.05）,g.29619320 G>T突变位点极显著偏离了Hardy-Weinberg平衡状态（P<0.01）。连锁不平衡分析显示,g.29619320 G>T和g.29619721 C>T位点具有强连锁效应。实时荧光定量PCR结果表明,UCP2基因在黔北麻羊各组织中均有不同程度的表达,其在肺脏中的表达量最高,其次是脾脏、肝脏和心脏,其余组织表达量较低。关联性分析结果表明,g.29614172 A>G位点AG和GG基因型个体体重、体斜长、胸宽、胸深、胸围及管围均显著高于AA基因型,g.29619320 G>T位点TT基因型个体胸深、胸围及管围均显著高于GT和GG基因型,g.29619721 C>T位点CT基因型个体体高显著高于CC和TT基因型,g.29620037 A>G位点GG基因型个体体重、胸深及管围均显著高于AG和AA基因型个体（P<0.05）。本试验结果初步揭示UCP2基因对黔北麻羊部分生长性状有显著影响,发现的4个SNPs位点可作为黔北麻羊经济性状选择的遗传标记。     

贵州地方山羊BMPR-ⅠB基因启动子区多态性及生物信息学分析

为探究骨形态发生蛋白受体ⅠB(bone morphogenetic protein receptorⅠB,BMPR-ⅠB)基因启动子区多态性及其与繁殖性状的相关性,本试验以贵州3个地方山羊品种(黔北麻羊、贵州黑山羊及贵州白山羊)为研究对象,通过构建混合DNA池PCR产物直接测序技术筛选SNPs位点,并采用多种生物信息学软件预测SNPs位点对核心启动子区、CpG岛和转录因子结合位点的影响。结果显示,BMPR-ⅠB基因启动子区存在4个SNPs位点:g.29893005 TC、g.29893016 CT、g.29893729 CG和g.29894370 CT。生物信息学软件预测得到BMPR-ⅠB基因核心启动子区和CpG岛,SNPs位点导致转录因子结合位点发生改变,g.29893005 TC和g.29893016 CT突变使原来的转录因子结合位点Sp1消失;g.29893729 CG突变使原来转录因子结合位点AP-1消失,从而产生新的转录因子结合位点USF;g.29894370 CT突变产生新的转录因子结合位点C/EBPalp,使原来的转录因子结合位点Sp1改变为C/EBPalp和Sp1。由此推测,SNPs位点对调控启动子功能元件可能存在重要影响。     

水牛MFSD2A基因多态性及其与产奶性状的关联分析

【目的】研究水牛主要促进因子超家族成员2a(major facilitator superfamily domain containing 2a,MFSD2A)基因多态性,并筛选与水牛产奶性状相关联的SNP。【方法】采集383头健康水牛血液DNA,利用PCR扩增和Sequenom MassARRAY飞行时间质谱技术对SNP位点进行筛选及基因型检测,并进行遗传多样性、连锁不平衡(LD)和单倍型分析,以及对不同基因型和单倍型与水牛产奶性状进行关联分析。【结果】在水牛MFSD2A基因中检出7个SNPs位点,多态信息含量(PIC)均在0.25~0.40之间,属于中度多态,除g.8290 T＞C位点外,其余位点均处于Hardy-Weinberg平衡状态(P＞0.05);关联分析结果表明,这7个SNPs位点与305 d产奶量均显著关联(P＜0.05),且突变杂合子基因型的305 d产奶量是3种基因型中最低的,g.568 C＞G和g.701 A＞G位点与乳脂率显著关联(P＜0.05),乳脂率从大到小依次为突变纯合子＞野生纯合子＞突变杂合子,g.568 C＞G、g.3326 G＞A、g.8290 T＞C和g.8523 C＞G位点与乳蛋白率显著关联(P＜0.05),7个SNPs位点与乳尿素氮关联均不显著(P＞0.05);连锁不平衡和单倍型分析结果明,g.568 C＞G、g.701 A＞G和g.3203 T＞G位点可组成一个单倍型模块(Block 1),g.8523 C＞G和g.9138 G＞T位点可组成另一个单倍型模块(Block 2)。其中g.568 C＞G和g.701 A＞G、g.568 C＞G和g.3326 G＞A、g.8523 C＞G和g.9138 G＞T位点的r²分别为0.87、0.37和0.48,处于强连锁不平衡状态。Block 1中的H1(CAG)单倍型乳蛋白率显著高于H2(CAT)和H3(GCT)(P＜0.05),Block 2中的D3(GG)单倍型305 d产奶量极显著高于D1(CG)和D2(CT)(P＜0.01)。【结论】筛选出的MFSD2A基因7个SNPs位点与305 d产奶量均显著关联,g.568 C＞G和g.701 A＞G位点与乳脂率显著关联,g.568 C＞G、g.701 A＞G、g.3326 G＞A、g.8290 T＞C和g.8523 C＞G位点与乳蛋白率显著关联,这7个SNPs位点与乳尿素氮均无显著关联。纯合子的基因型在高305 d产奶量和乳脂率的选择上更有优势,H1(CAG)和D3(GG)单倍型是提高水牛305 d产奶量和乳蛋白率的有利单倍型。     

中国荷斯坦牛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位点可作为潜在分子标记,为中国荷斯坦牛的平衡育种提供理论依据。     

延边黄牛脂蛋白脂肪酶基因多态性及其与肉质性状的关联分析

为探究延边黄牛脂蛋白脂肪酶（lipoprteinlipase,LPL）基因遗传多态性及其对肉质性状的影响,本试验以80头健康的30月龄延边黄牛背最长肌作为试验材料,运用PCR直接测序法和高分辨率熔解曲线技术（HRM）对LPL基因的SNP位点进行检测,分析了试验群体的遗传效应,并与延边黄牛的肉质性状数据进行关联分析。PCR直接测序结果表明,延边黄牛LPL基因存在2个多态性位点（g.6215 A>G和g.18341 C>T）,其中g.6215 A>G位点表现为2种基因型：AA与AG,优势等位基因为A;g.18341 C>T位点表现为2种基因型：CC和CT,C为优势等位基因;对文献中2个位点的检测结果表明,g.355427 T>A位点表现为2种基因型：AA与AT,A为优势等位基因;c.322 G>A位点的HRM分型结果显示并未分型。这些多态性位点全部呈中度多态（0.25<PIC<0.5）,并处于Hardy-Weinberg平衡状态。HRM分型结果显示,这3个SNPs位点与延边黄牛肉质性状存在显著相关：g.6215 A>G位点与延边黄牛胴体重及脂肪、蛋白、水分含量有显著关联,表现为AG基因型个体胴体重、脂肪含量显著高于AA基因型个体（P<0.05）,AA基因型个体蛋白及水分含量显著高于AG基因型个体（P<0.05）;g.18341 C>T位点与延边黄牛宰前活重、胴体重及背膘厚有显著关联,表现为CT基因型个体这些性状均显著高于CC基因型个体（P<0.05）;g.355427 T>A位点AT基因型个体宰前活重、背膘厚和脂肪含量均高于AA基因型个体（P<0.05）。综上,LPL基因g.6215 A>G、g.18341 C>T和g.355427 T>A位点对延边黄牛的肉质有显著影响,LPL基因可以作为延边黄牛品种选育改良工作的优势候选基因和有效分子标记。     

ADSL基因启动子区遗传变异与鸡肉冻藏新鲜度的相关性分析

试验旨在探究鸡腺苷琥珀酸裂解酶（adenylosuccinate lyase,ADSL）基因启动子区多态性及其与鸡肉冻藏新鲜度的相关性。以161只青脚麻鸡为研究对象,通过DNA混池测序法检测ADSL基因启动子区遗传变异。采用等位基因特异性PCR完成个体基因型分析,并分析了各基因型与鸡胸肌肉色、体重及胸肌冷冻60 d后新鲜度K值和pH的相关性。结果显示,在ADSL基因启动子区ATG上游1 670 bp处发现1个SNP位点（c.-1670 C>A）,其AA基因型个体肌肉冷冻60 d后新鲜度K值（23.61%）极显著低于CA（33.00%）和CC（33.56%）基因型（P<0.01）,即冻藏60 d后AA基因型个体肌肉较CA和CC基因型新鲜,CA和CC基因型间新鲜度K值差异不显著（P>0.05）;各基因型间鸡体重、胸肌肉色及胸肌冷冻60 d后pH均无显著差异（P>0.05）。生物信息学预测发现,c.-1670 C>A导致ADSL基因启动子区转录因子结合位点的显著改变,该突变位点可能造成ADSL基因的表达差异。综上,ADSL基因启动子区c.-1670 C>A与青脚麻鸡胸肌冷冻60 d后新鲜度K值具有显著相关性,该位点可作为青脚麻鸡肉品新鲜度的候选分子遗传标记。     

QDPR基因克隆及其在乌蒙凤鸡不同生长阶段组织表达研究

旨在克隆乌蒙凤鸡醌型二氢生物喋呤还原酶（quinoid dihydropteridine reductase，QDPR）基因，并检测其在不同性别乌蒙凤鸡不同生长阶段各组织中的表达差异，以研究QDPR基因在乌蒙凤鸡生长发育过程中的调控作用。本试验对乌蒙凤鸡QDPR基因CDS区进行PCR扩增和克隆，并对得到的序列进行生物信息学分析，通过实时荧光定量PCR检测QDPR基因在4、12、20和28周龄乌蒙凤鸡公鸡及母鸡的心、肝、脾、肺、肾、胸肌和腿肌组织中的相对表达量。结果表明，QDPR基因开放阅读框长度为417 bp，共编码139个氨基酸，编码蛋白为酸性不稳定蛋白，乌蒙凤鸡的QDPR基因序列与原鸡、日本鹌鹑和珠鸡的QDPR基因同源性较高；QDPR基因在不同性别乌蒙凤鸡不同生长阶段各组织中均有表达，其中，公鸡4、12和20周龄肺组织中QDPR基因表达量最高，极显著高于同一周龄其他组织（P<0.01），母鸡4周龄肝、脾、肺中QDPR基因相对表达量为所有时期最高，极显著高于12、20周龄（P<0.01），公鸡大多数组织中QDPR基因表达量随时间增长总体呈现出先升后降的趋势，母鸡为先降后升。结果显示，QDPR基因在乌蒙凤鸡内脏组织中的表达均高于肌肉组织，且公母鸡不同时期表达情况有所差异，试验结果可为进一步研究QDPR基因在鸡生长发育中的调控作用提供参考。     

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基因与繁殖、泌乳性状的遗传关联，为中国荷斯坦牛高产高效选育提供了可用的遗传标记。     

山羊ZBP1基因多态性及与产羔性能的关联分析

【目的】试验旨在分析山羊Z-DNA结合蛋白(Z-DNA binding protein 1,ZBP1)基因3′-端非翻译区(3′-UTR)、内含子1和外显子7区域中单核苷酸多态性(single nucleotide polymorphism, SNP)位点与产羔性能之间的关联性,为高繁殖力山羊分子育种提供新的遗传标记。【方法】选取麻城黑山羊、黑头羊和波尔山羊作为研究样本,采集血样提取DNA,根据转录组数据所选SNP位点在山羊ZBP1基因序列中的位置信息设计引物,利用多重单碱基延伸SNP分型技术(SNaPshot)分析所选SNP位点的遗传多样性,采用一般线性模型(GLM)对ZBP1基因SNP位点与3个品种山羊的产羔性能进行关联分析。【结果】山羊ZBP1基因中存在12个潜在的SNPs位点,其中7个SNPs位点(g.9734 A>G、g.9772 T>C、g.352 C>T、g.955 C>T、g.1880 G>A、g.2566 T>C和g.7919 G>A)具有多态性,除g.9734 A>G在麻城黑山羊群体中仅有AA和GA 2种基因型外,其余...