miR-29a在酉州乌羊不同体组织中表达特征及其在B16细胞中的动态表达分析

为探究miR-29a在酉州乌羊不同体组织表达特征及在B16细胞中的表达特性,本研究采用RT-qPCR方法检测miR-29a在酉州乌羊不同体组织以及B16细胞不同发育时期的表达水平。结果表明,miR-29a在酉州乌羊各体组织均有表达,其中在大脑中的相对表达量最高,肺和皮肤中的相对表达量最低;同时在渝东白羊皮肤中的相对表达量是酉州乌羊皮肤的3.44倍,且差异极显著(P<0.01);在细胞增殖过程中,miR-29a相对表达量极显著增加(P<0.01);在分化初期miR-29a相对表达量较高,之后趋于稳定且差异不显著(P>0.05)。在加入分化培养基后,细胞的形态发生变化,黑色素分泌增多。综上所述,miR-29a在酉州乌羊各组织中广泛表达,miR-29a的相对表达量随着B16细胞的增殖而递增,但随着B16细胞的分化而降低,推测miR-29a与B16细胞的增殖与分化密切相关。本研究结果为明确miR-29a在酉州乌山羊各组织中的表达特性及其是否参与了黑色素细胞行为学的调控提供了一定的理论证据。     

陇东绒山羊KRTAP2-1基因鉴定及其对羊绒性状影响的研究

角蛋白关联蛋白(KAPs)是羊绒纤维的主要结构成分,但尚未在山羊基因组中鉴定出KRTAP2-1基因,本研究以人类KRTAP2-1基因为模板,在山羊基因组中进行同源性搜索,在19号染色体上发现了一个399 bp的开放阅读框。采用聚合酶链反应-单链构象多态性(PCR-SSCP)方法,在249只陇东山羊上检测到7条不同的序列(A-G)。进化树结果表明,这7条序列与人类和绵羊的KAP2-1序列有较高的同源性,表明该序列是山羊KRTAP2-1的7个等位基因。测序结果表明,KRTAP2-1基因中存在12个单核苷酸多态性(SNPs)(5个SNPs位于编码区,7个SNPs位于非编码区),以及c.-61_-63delCTC和c.93_95delCCG两个碱基缺失,其中c.93_95delCCG引起了第32位精氨酸的移码缺失。山羊KAP2-1蛋白含有丰富的半胱氨酸、脯氨酸、苏氨酸和丝氨酸。相关性分析表明,等位基因C的存在与较小的羊绒纤维直径相关(存在:13.4±0.05 μm;缺失:13.6±0.03 μm;P=0.010)。本研究表明,山羊KRTAP2-1基因有丰富的多态性,可作为陇东绒山羊羊绒纤维直径选育的分子标记用于生产实践。     

LncRNA XLOC_015132对山羊黑色素沉积的影响及其lncRNA-miRNA-mRNA轴的预测

LncRNA在哺乳动物的基因组中被广泛转录,能够通过相应的miRNA和靶基因促进黑色素瘤的侵袭和转移。为探究LncRNA XLOC_015132的组织表达谱及其与黑色素沉积之间的关系,以酉州乌羊和酉州本地白山羊为试验对象,检测LncRNA XLOC_015132在各组织中以及黑色素瘤细胞中的表达情况;同时,通过靶miRNA和靶基因预测,构建lncRNA-miRNA-mRNA调控网络。qRT-PCR结果显示,LncRNA XLOC_015132在酉州乌羊脾脏、肺、肾脏、大脑、皮肤等组织中的表达量显著高于本地白山羊(P<0.05);LncRNA XLOC_015132在B16-F10黑色素瘤细胞中的定量分析结果显示,随着B16-F10细胞的增殖,LncRNA XLOC_015132表达量呈递增趋势;LncRNA XLOC_015132与互作分子的生物信息学分析结果显示,LncRNA XLOC_015132/miR-150-5p/CCND2轴可能通过激活PI3K-AKT、Wnt等信号通路,在黑色素沉积中发挥重要作用。综合上述结果,LncRNA XLOC_015132的表达与黑色素沉积可能存在正相关关系。本研究结果为深入探究山羊皮肤黑色素沉积的lncRNA调控理论研究提供了基础。     

大肠杆菌对藏羊子宫内膜上皮细胞相关炎性因子表达的影响

为了探究大肠杆菌(Escherichia coli)对藏羊子宫内膜上皮细胞(EECs)相关炎性因子表达的影响,本研究采用组织块贴壁法和酶消化法培养藏羊EECs;用不同感染复数(MOI)进行E. coli感染试验,通过实时荧光定量PCR(qRT-PCR)及酶联免疫(ELISA)技术,分析白细胞介素1β(IL-1β)、白细胞介素6(IL-6)、白细胞介素8(IL-8)和肿瘤坏死因子α(TNF-α)等炎性因子的RNA和蛋白表达水平。结果表明,培养液中加入50 ng·mL^-1表皮生长因子(EGF)能成功培养出藏羊EECs,纯度达到98%以上;E. coli感染试验中,当MOI为1:1时,IL-1β、IL-6和IL-8基因和蛋白的表达量与对照组相比差异不显著,而TNF-α表达量显著升高(P<0.05);当MOI为20:1时,IL-8基因和蛋白的表达量显著升高(P<0.05);当MOI为50:1时,IL-1β、IL-6基因和蛋白的表达量显著升高(P<0.05),而TNF-α表达量与对照组差异不显著(P>0.05)。结果表明,酶消化法可成功培养藏羊EECs;不同MOI E.coli感染后,藏羊EECs中IL-1β、IL-6、IL-8和TNF-α等炎性因子的表达量较对照组均显著升高(P<0.05),它们可作为识别子宫内膜炎病理变化程度的指标。本研究结果为预防子宫内膜炎、提高藏羊的生产性能以及生态养殖提供了理论依据。     

海沃德猕猴桃货架期预测模型的建立

为建立海沃德猕猴桃货架期品质变化动力学模型,以实现对海沃德猕猴桃货架期的预测。本试验将经0.5 μL·L^-1 1-甲基环丙烯(1-MCP)处理(处理组)与不处理(对照组)的猕猴桃于0、4、15、25℃条件下贮藏,测定其硬度、可溶性固形物含量、可滴定酸含量、Vc含量;并基于Arrhenius理论建立货架期预测模型,选取20℃样本对该模型进行验证。结果表明,1-MCP对试验设置温度条件下贮藏的海沃德猕猴桃硬度、可滴定酸含量和Vc含量的降低以及可溶性固形物含量的升高均有一定的延缓作用;同时对其硬度和Vc含量变化情况进行拟合,比较拟合后的反应速率常数k,发现温度越高Vc含量和硬度变化越快。比较零级反应和一级反应线性回归决定系数R²可知,对照组和处理组Vc含量变化动力学方程遵循一级反应动力学,对照组R²≥0.97,处理组R²≥0.94;对照组和处理组硬度含量变化动力学方程遵循零级反应动力学,对照组R²≥0.96,处理组R²≥0.92。对照组硬度下降及Vc含量下降的lnk和1/T的R²>0.98,处理组硬度下降及Vc含量下降的lnk和1/T的R²>0.89。将预测值和实测值进行比较,结果显示,对照组猕猴桃货架期预测值和实测值相对误差均在±10%以内。综上表明,用硬度变化预测猕猴桃货架期,操作简单,相对误差较低。本研究结果在猕猴桃贮藏、销售中具有较好的应用价值。     

高羊茅硝态氮转运蛋白基因NRT1.1的克隆及表达模式分析

为探究高羊茅(Festuca arundinacea)硝态氮转运蛋白基因(NRT1.1)的表达模式,本研究以黔草1号高羊茅为试验材料,采用RACE和RT-qPCR技术对高羊茅NRT1.1基因的cDNA全长序列进行扩增,并对其不同胁迫处理下的表达情况进行分析。生物信息学分析发现,高羊茅NRT1.1的理论等电点为4.81,平均亲小性为0.919,含有约32.63% α-螺旋、7.63% β-转角和53.73%不规则卷曲。结果表明,NRT1.1基因的cDNA序列全长为2 328 bp,编码606个氨基酸,预测蛋白质分子量为193.9 kDa,且高羊茅NRT1.1与黑麦草NRT1.1氨基酸序列的相似性最高。RT-qPCR表达分析发现,高羊茅叶片NRT1.1受低氮处理0.5~1 h时表达量达到峰值,显著(P< 0.05)高于对照组;在干旱和热处理下,NRT1.1表达量分别在6 h和12 h时达到峰值,且显著(P< 0.05)高于对照组;在盐处理下,仅在6 h时NRT1.1表达量高于对照组,其余时间均受显著(P< 0.05)抑制。本研究结果为解析高羊茅NRT1.1基因的表达模式提供了分子生物学基础。     

茶树CsDREB-A2转录因子基因的克隆与非生物胁迫响应分析

干旱应答元件结合蛋白(DREB)类转录因子在植物逆境信号转导途径中具有重要的调控作用。为了解AP2/ERF转录因子在茶树逆境胁迫的分子调控机理,本研究从茶树龙井43叶片的cDNA中克隆得到一个编码CsDREB-A2转录因子的基因;对CsDREB-A2基因及其编码蛋白序列特征进行分析,并利用实时荧光定量PCR法检测该基因在茶树不同非生物胁迫处理下的表达水平。结果表明,CsDREB-A2基因开放阅读框为1 056 bp,编码351个氨基酸,其编码氨基酸序列具有AP2保守结构域,包含典型的YRG元件和WLG基序。AP2结构域第14、第19位氨基酸分别为缬氨酸和谷氨酸。拟南芥AP2/ERF家族转录因子的进化分析表明,该转录因子属于DREB亚族的A2组。CsDREB-A2相对分子质量为39 080 Da,理论等电点为5.32,属于亲水性蛋白,主要由α-螺旋和随机卷曲组成,无序化特征明显且存在一个LM无序区域;可能定位于细胞核,不存在信号肽和跨膜结构,属于非分泌蛋白。CsDREB-A2基因在高温(38℃)和干旱(200 g·L^-1 PEG)胁迫下均能快速诱导表达,并显著高于对照,分别在处理4、2 h达到最大值,为对照的20.70和42.90倍,植物在渗透胁迫下,可能通过ABA信号途径调节该基因对干旱的耐受性。高盐(200 mmol·L^-1 NaCl)胁迫下CsDREB-A2基因的表达受抑制,推测其可能存在负调控结构域降低该基因在盐胁迫下的表达量。本试验结果为研究DREB类转录因子在茶树抗逆胁迫的分子调控机制提供了一定的理论参考。     

竹类植物中CYP85A1基因的克隆及表达分析

油菜素内酯(BRs)是一种促进植物茎秆伸长的甾类植物激素,CYP85A1基因参与了BRs合成中重要的催化反应步骤。为阐明CYP85A1基因在竹类植物节间伸长中的作用,本研究以20个竹类植物为试验材料,利用生物信息学分析了CYP85A1基因的特点,并利用实时荧光定量PCR技术分析了该基因在孝顺竹及其变种凤尾竹的不同组织和不同发育阶段的相对表达量。结果表明, 20个竹种中CYP85A1基因序列在DNA水平上的同源性为92.23%,cDNA序列的同源性为98.79%,氨基酸序列的同源性为97.86%;这些基因均含有9个外显子和8个内含子,编码465~480个氨基酸,蛋白分子量约为54 kDa,理论等电点均约为9.2。氨基酸序列中均含有细胞色素P450家族保守的血红素结合域、富含脯氨酸和氧的结构域及Glu-X-X-Arg结构域;CYP85A1基因孝顺竹和凤尾竹母竹的叶中表达量均最高,茎中次之,在根中的表达量最低。孝顺竹和凤尾竹幼竹中,该基因在顶部的表达量明显高于第2、第4、第5节间中的表达量;CYP85A1基因在快速伸长阶段的表达量明显高于生长缓慢时期,节间伸长速率越快该基因的表达量越高。综上,CYP85A1基因在竹类植物茎秆及节间伸长发育过程中具有促进作用。本研究结果为阐明CYP85A1在竹类植物中的生物学功能奠定了一定的理论基础。     

桃PpMADS1的克隆及对类胡萝卜素合成基因的调控研究

为探究桃MADS转录因子基因的序列特征、表达模式及其对类胡萝卜素合成的调控作用,本研究从桃果实中克隆得到一个MADS-box家族基因,命名为PpMADS1。结果表明,PpMADS1的开放阅读框(ORF)为732 bp,编码243个氨基酸,蛋白具有亲水性,稳定性较差,其二级结构以α-螺旋为主要构成元件,且与三级结构相似。亚细胞定位预测显示PpMADS1蛋白定位于细胞核。氨基酸序列分析及系统进化分析结果表明,PpMADS1属于MADS-box转录因子AG亚族,与乌梅亲缘关系最近。实时荧光定量PCR结果表明,PpMADS1基因在芽和叶中几乎不表达,在花和果实中大量表达,黄肉桃品种金丽果实成熟后期PpMADS1表达量显著高于白肉桃品种湖景果实。双荧光素酶试验结果表明,PpMADS1转录因子对PpPSY和PpCHYB启动子具有转录激活作用。因此推测PpMADS1可能通过调控PpPSY和PpCHYB基因表达促进桃果实类胡萝卜素合成。本研究为桃MADS-box转录因子的进一步功能分析奠定了理论基础。     

谷子SBP蛋白基因家族的鉴定及表达分析

为了深入分析SBP基因家族在谷子中的功能作用,本研究在谷子全基因组水平,通过NCBI blast及HMMER search等生物信息学方法进行了谷子SBP基因家族鉴定,分析了SBP蛋白理化性质及系统发育,利用RT-qPCR分析了谷子SBP基因在PEG和ABA胁迫下的表达情况。结果表明,谷子SBP基因家族有19个成员,位于9条染色体上,编码的氨基酸序列长度为199~1 118 aa。谷子SBP蛋白系统发育树结果显示,该蛋白家族可分为两大类:第一类分为六个亚类,第二类仅有一个SBP蛋白。谷子SBP蛋白家族成员的外显子数从2~11个不等,且均为亲水性蛋白。亚细胞定位预测结果显示,除SBP7外,所有的SBP蛋白在细胞核中均有信号。谷子SBP蛋白与拟南芥SBP蛋白在进化关系上较远,而与高粱、水稻SBP蛋白进化关系较近。谷子SBP基因家族启动子区域的顺式作用元件主要为植物激素、干旱、光照及病原菌、创伤类的响应元件。RT-qPCR分析结果显示,PEG和ABA胁迫诱导下,SBP基因在叶和茎中的表达均上调,表达量增幅表现为:叶>茎。PEG胁迫后SBP基因的表达增长量显著高于外源ABA诱导;其中,SiSBP13的表达量增幅最大,PEG胁迫4 h叶和茎中的表达量分别达到未胁迫时的180倍和15倍,ABA诱导4 h茎中的表达量最高达到未胁迫时的8倍,而叶中的表达量持续升高,24 h时达到未胁迫时的28倍。本研究为下一步深入研究SiSBPs基因在谷子生长发育阶段中的功能奠定了基础,也为谷子功能基因的挖掘利用提供了候选基因。     

贵州本地山羊STAT5A基因多态性及其与生长性状关联分析

本研究通过探讨STAT5A基因SNP与山羊生长性状关联性,旨在为山羊选种选育提供更好的科学依据。实验以贵州白山羊和贵州黑山羊为研究对象,采用DNA池法及PCR-SSCP技术检测STAT5A基因单核苷酸多态性。结果发现,在贵州白山羊和贵州黑山羊STAT5A基因内含子10均检测到1个SNP位点G127A,表现为3种基因型,分别命名为GG、GA和AA。基因型与生长性状关联分析显示,贵州黑山羊基因型GA个体的胸围和体重指标显著高于基因型GG个体（p〈0.05）,而其余3个指标均差异不显著（p〉0.05）;贵州白山羊基因型GA个体的体斜长、胸围和管围指标显著高于基因型GG个体（p〈0.05）,而其余指标均差异不显著（p〉0.05）。研究结果提示：STAT5A基因可能是影响山羊体重、体斜长、管围和胸围的主效基因或与主效基因连锁,G127A位点可望作为提高山羊个体生长性状的分子遗传标记。     

山羊FSHR基因第10外显子的PCR-SSCP检测及其序列分析

利用PCR-SSCP技术检测山羊(Caprahircus)包括西农萨能奶山羊、关中奶山羊、陕南白山羊、安哥拉山羊和波尔山羊173个个体FSHR基因第10外显子的单核苷酸多态性(SNP)。结果未发现SNP位点。测序后获得山羊FSHR基因第10外显子的核苷酸序列,并在NCBI数据库中获得GenBank登录No.DQ069909和DQ069910。通过DNA序列分析发现,FSHR基因第10外显子第120位碱基不存在C→T的转换,也不存在颠换等其它遗传变化。山羊、绵羊(OvisariesL.)和普通牛(Bostaurus)FSHR基因第10外显子序列同源性比较和聚类分析结果表明,山羊、普通牛和绵羊该部分序列的相似性最高为99.3%;在物种间比较中,绵羊和普通牛纯合子该基因外显子序列的不相似性最高为3.4%;据FSHR基因外显子序列构建的分子系统树结果显示,山羊、绵羊和普通牛物种内的个体各自聚为一类;山羊和绵羊先聚为一类,然后再与普通牛聚为一类。提示FSHR基因第10外显子的核苷酸序列适合于物种间的动物分子树的构建。     

of Goat Follicle-stimulating Hormone Receptor (FSHR) Gene

利用PCR-SSCP技术检测山羊(Capra hircus)包括西农萨能奶山羊、关中奶山羊、陕南白山羊、安哥拉山羊和波尔山羊173个个体FSHR基因第10外显子的单核苷酸多态性 (SNP)。结果未发现SNP位点。测序后获得山羊FSHR基因第10外显子的核苷酸序列， 并在NCBI数据库中获得GenBank登录No.DQ069909和DQ069910。通过DNA序列分析发现， FSHR基因第10外显子第120位碱基不存在C→T的转换，也不存在颠换等其它遗传变化。山羊、绵羊 (Ovisaries L.) 和普通牛( Bos taurus) FSHR基因第10外显子序列同源性比较和聚类分析结果表明，山羊、普通牛和绵羊该部分序列的相似性最高为99.3%；在物种间比较中，绵羊和普通牛纯合子该基因外显子序列的不相似性最高为3.4%；据FSHR基因外显子序列构建的分子系统树结果显示，山羊、绵羊和普通牛物种内的个体各自聚为一类；山羊和绵羊先聚为一类，然后再与普通牛聚为一类。提示FSHR基因第10外显子的核苷酸序列适合于物种间的动物分子树的构建。     

神经肽Y-Y1受体基因多态性及其与山羊繁殖性能关系的研究

设计5对引物(P1～P5),用限制性片段长度多态性（RFLP）和单链构象多态性（SSCP）方法检测神经肽Y1受体（neuropeptide Y Y1 receptor,NPY-Y1R）基因在性早熟、高繁殖力品种（济宁青山羊）和性晚熟、中等繁殖力（波尔山羊）以及性晚熟、低繁殖力品种（安哥拉山羊和内蒙古绒山羊）中的单核苷酸多态性,同时研究该基因对济宁青山羊性早熟和高繁殖力的影响。结果发现仅引物P2、P4和P5的扩增片段具有多态性。引物P2在济宁青山羊、安哥拉山羊和内蒙古绒山羊中检测到AA、AB和BB型,在波尔山羊中只检测到AA和AB型;测序发现BB型与AA型相比存在1个突变（C625G）,没有引起氨基酸变化;济宁青山羊AA、AB与BB基因型之间产羔数差异均不显著（P>0.05）。引物P4在4个山羊品种中均检测到CC和CD基因型;测序发现CD型与CC型相比发生了G1141A突变,没有引起氨基酸变化;济宁青山羊CC和CD基因型之间产羔数差异不显著（P>0.05）。引物P5在济宁青山羊中检测到EE和EF基因型,在其余3个山羊品种中只检测到EE型;测序发现EF型与EE型相比发生了C1330T突变,没有引起氨基酸变化;济宁青山羊与波尔山羊、安哥拉山羊和内蒙古绒山羊之间NPY-Y1R基因型分布差异均达到极显著（P<0.01）水平,后3个品种之间均无显著差异（P>0.05）;EF型济宁青山羊产羔数均值比EE型多0.58只（P<0.01）。本研究结果初步表明NPY-Y1R基因的F等位基因与济宁青山羊性早熟和高繁殖力相关。     

催乳素受体基因多态性及其与济宁青山羊高繁殖力关系的研究

设计8对引物，采用PCR-SSCP技术检测催乳素受体(prolactin receptor，PRLR)基因外显子3至外显子9和内含子9在高繁殖力山羊品种（济宁青山羊）和低繁殖力山羊品种（辽宁绒山羊、安哥拉山羊）中的单核苷酸多态性，同时研究该基因对济宁青山羊高繁殖力的影响。结果表明：8对引物中，只有扩增内含子9的P8引物的扩增区域具有多态性，其余7对引物的扩增片段都不存在多态性。对于P8扩增片段，在济宁青山羊、辽宁绒山羊和安哥拉山羊中都只检测到AA型和AB型。 测序表明AB型与AA型相比有1处突变（209A→G），该突变位于PRLR的第9内含子。济宁青山羊2 种基因型之间产羔数的最小二乘均值差异不显著（P>0.05）。这些结果初步表明检测的PRLR基因位点对济宁青山羊高繁殖力没有显著影响。     

山羊美臀突变基因（CLPG）的多态性鉴定及其与生产性能的关系

Callipyge绵羊的臀部肌肉明显肥大,呈极性超显性遗传。研究克隆了山羊（Capra hircas）美臀突变基因(CLPG) 250 bp的序列(GenBank登录号：EU753362), 对其多态性及其与山羊生产性能的关系进行了分析。结果表明,所扩增片段与绵羊和猪相应片段的同源性分别为96.04%和88.65%。PCR-SSCP分析表明,在该DNA区段存在多态性。测序结果显示,在对应于绵羊（Ovis aries）的CLPG基因没有发现A→G的突变,但在该位点下游147 bp处检测到一个A→C的颠换,记为SNP216。对波尔山羊(n　=63)、莱芜黑山羊(n　=70)、鲁波杂交山羊(n　=40)、鲁北白山羊(n　=29)和内蒙古阿拉善白绒山羊(n　=115)共5个群体的多态性进行了检测,结果表明：除鲁北白山羊以外,A等位基因均为优势等位基因;波尔山羊、莱芜黑山羊、鲁北白山羊和鲁波杂交山羊群体均处于Hardy-Weinberg平衡状态（P　＞0.05）,而内蒙古阿拉善白绒山羊群体显著偏离Hardy-Weinberg平衡状态（P　＜0.01）;不同基因型的内蒙古阿拉善白绒山羊的初生重、从出生到断奶时的增重及产绒量的最小二乘均值差异不显著（P　＞0.5）。     

促性腺激素释放激素受体(GnRHR)基因多态性及其与济宁青山羊高繁殖力关系

设计8对引物, 应用聚合酶链式反应-单链构象多态性（PCR-SSCP）技术检测促性腺激素释放激素受体（gonadotropin releasing hormone receptor, GnRHR）基因外显子1、外显子2和外显子3在高繁殖力山羊(Caprar hircus)品种（济宁青山羊）和低繁殖力山羊品种（安哥拉、波尔和内蒙古绒山羊）中的单核苷酸多态性,同时研究该基因对济宁青山羊高繁殖力的影响。结果发现引物P4扩增片段具有多态性,其余7对引物的扩增片段都不存在多态性。对于引物P4扩增片段,在济宁青山羊、波尔山羊和内蒙古绒山羊中均检测到AA、AB和BB基因型,而在安哥拉山羊中只检测到AA和AB基因型,未检测到BB基因型;测序分析发现BB型与AA型相比在外显子1有1个突变（757G→A）,但未引起氨基酸改变;济宁青山羊AA、AB和BB基因型频率分别为0.623、0.300和0.077,BB型济宁青山羊平均产羔数比AB型和AA型分别多0.69只（P　<　0.05）和0.82只（P　<　0.05）,AA型和AB型济宁青山羊产羔数的最小二乘均值差异不显著（P　>　0.05）。     

转化生长因子β1基因(TGF-β1)外显子2的单核苷酸多态性及其与两个山羊品种产羔数的相关分析

为寻找与山羊(Capra hircus)产羔数密切相关的分子标记,加快山羊育种进程,研究转化生长因子β1基因(TGF-β1)遗传多态性与山羊产羔数的相关性,根据牛的TGF-β1基因序列设计3对引物,采用PCR-SSCP.和DNA测序技术检测山羊TGF-β1基因外显子2和内含子3的多态性,同时用最小二乘法研究了其多态性与产羔数的关系。结果表明,在波尔山羊和陕南白山羊中,仅第二外显子148bp位点碱基发生A→G的突变,导致苏氨酸(Thr)变为丙氨酸(Ala);多态位点均以A为优势等位基因,基因频率分别为0.530和0.648。在波尔山羊的1～3胎产羔数中,AB型个体的产羔数显著高于BB和AA型个体(P<0.05);在第四胎产羔数和平均产羔数中,AB型个体的产羔数显著高于AA(P<0.05),显著高于BB(P<0.05);BB型个体的平均产羔数显著高于AA型(P<0.05)。在陕南白山羊中,AB型个体的2～4胎产羔数和平均产羔数显著高于AA和BB型个体(P<0.05);在第一胎中,AB型个体显著高于AA型个体(P<0.05)。研究结果表明TGF-β1基因多态位点与产羔数显著相关,可以用于山羊分子遗传育种的候选基因。