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

设计5对引物, 采用PCR-SSCP技术检测催乳素受体(prolactin receptor, PRLR)基因外显子10及部分3′非翻译区在高繁殖力山羊品种(济宁青山羊)和低繁殖力山羊品种(辽宁绒山羊、波尔山羊和安哥拉山羊)中的单核苷酸多态性, 同时研究该基因对济宁青山羊高繁殖力的影响。结果表明: 首次拼接出的山羊PRLR基因外显子10及部分3′非翻译区的核苷酸序列长度为1,118 bp, 与已公布的绵羊、牛、人PRLR基因mRNA相应序列的同源性分别为98.33%、93.92%、74.52%, 与已公布的羊驼PRLR基因部分序列的同源性为78.29%。引物P1、P2与P4扩增片段具有多态性, 其余2对引物的扩增片段不存在多态性。对于P1扩增片段, 在济宁青山羊和辽宁绒山羊中检测到AA型和AB型, 在安哥拉山羊中检测到AA型和AC型, 在波尔山羊中只检测到AA型; 克隆测序表明AB型与AA型相比有两处突变(186G→A和220T→C), 分别导致氨基酸由天冬氨酸变为天冬酰胺、亮氨酸变为脯氨酸; AC型与AA型相比有1处突变(140A→G), 该突变没有导致氨基酸变化; 济宁青山羊AA和AB基因型之间产羔数的最小二乘均值差异不显著(P>0.05)。对于P2扩增片段, 在济宁青山羊、辽宁绒山羊和波尔山羊中都检测到DD型和DE型, 而在安哥拉山羊中只检测到DD型; 克隆测序表明DE型和DD型相比有两处突变(52G→A和122G→A), 其中122 bp处的突变导致氨基酸由精氨酸变为甘氨酸; 济宁青山羊DD和DE基因型之间产羔数的最小二乘均值差异不显著(P>0.05)。对于P4扩增片段, 在济宁青山羊中检测到FF型和FG型, 在辽宁绒山羊中检测到FF型和GG型, 在波尔山羊中只检测到FF型, 在安哥拉山羊中检测到FF型、FG型和GG型; 克隆测序表明GG型和FF型相比在扩增片段的143 bp处发生1处碱基突变(A→G), 并导致氨基酸由蛋氨酸变为缬氨酸; FG基因型济宁青山羊产羔数最小二乘均值比FF基因型的多0.76只(P<0.05)。研究结果初步表明: PRLR基因可能是控制济宁青山羊多胎性能的一个主效基因或是与之存在紧密遗传连锁的一个标记。     

FSHβ 基因PCR-SSCP多态性及其与济宁青山羊高繁殖力关系的研究

采用PCR-SSCP技术检测促卵泡素b亚基(follicle-stimulating hormone β, FSHβ)基因5′调控区、外显子1和外显子2在高繁殖力山羊品种(济宁青山羊)和低繁殖力山羊品种(辽宁绒山羊、波尔山羊、安哥拉山羊)中的单核苷酸多态性, 同时研究该基因对济宁青山羊高繁殖力的影响。结果表明: 山羊与绵羊的FSHβ 基因该段核苷酸序列同源性为98%。9对引物中, 只有P9的扩增片段存在多态性。P9的扩增片段在济宁青山羊和辽宁绒山羊中检测到AA、AB和AC 3种基因型; 在波尔山羊中检测到AA、CC和AC 3种基因型; 在安哥拉山羊中检测到AA、BB、CC、AB、AC和BC共6种基因型。测序分析发现BB型与AA型相比在外显子2的第94 bp处有G→A突变, 并引起氨基酸改变(丙氨酸→苏氨酸); CC型与AA型相比在外显子2的第174 bp有一处C→T沉默突变。济宁青山羊AA、AB和AC基因型频率分别为0.686、0.137和0.177。AA基因型济宁青山羊产羔数最小二乘均值比AB基因型的多0.78只(P<0.05), 比AC基因型的多0.64只(P<0.05)。     

三种绵羊BMP15基因第二外显子多态性及与产羔数的相关性分析

绵羊存在影响多胎性状的不同主效基因,选择影响Romney Hanna绵羊和Cambridge绵羊高繁殖力的骨形态发生蛋白15 (bone morphogenetic protein 15, BMP15)为候选基因,采用PCR-SSCP的方法检测BMP15基因外显子Ⅱ第747位点(T747→C)和755位点(T755→C)在蒙古羊、甘肃高山细毛羊、小尾寒羊三种绵羊母羊中的多态性,同时还研究了上述两处突变对三种绵羊产羔数的影响。表明:(1)一共检测到野生纯合型AA、突变杂合型AB (T747→C)、AC (T755→C)三种不同的基因型,AA为优势基因型,A为优势等位基因;(2)三种基因型在甘肃高山细毛羊中均被检测到,而蒙古羊和小尾寒羊中未检测出AB基因型;(3)突变杂合型蒙古羊(AC)比野生纯合型(AA)的平均产羔数多0.27只(p<0.05)。(4)AC的基因型频率,双羔母羊和多羔母羊均高于单羔母羊。根据以上实验推测,BMP15第755位点发生的T→C突变(AC型)对蒙古羊一胎产双羔影响十分显著,甘肃高山细毛羊中AC基因型的绵羊其产羔数有比AA基因型和AB基因型多的趋势,因此该位点可能是一个影响绵羊高繁殖力潜在的DNA标记。     

生长激素基因多态性与山羊体重性状的关系

以鲁北白山羊、波尔山羊以及波尔山羊与鲁北白山羊的杂交一代、回交一代共计224只山羊为材料,根据山羊生长激素基因5′调控区的序列(GenBank登录号:D00476)设计两对引物,用PCRSSCP法进行多态性分析。多态性片段的纯合基因型经克隆测序,发现共有5处突变。对突变位点产生的不同基因型与体重、体尺性状进行分析表明:第一对引物扩增片段波尔山羊AA型个体的初生重、周岁重显著高于BB型和AB型(P<0.05);杂交一代中断奶重、周岁重也以AA型偏高,但不同基因型间的差异均不显著(P>0.05);而鲁北白山羊BB基因型的体重相对于另外两种基因型的偏低,且断奶重的差异达到显著水平(P<0.05)。第二对引物扩增片段不同基因型对体重、体尺没有显著影响。由此初步推断生长激素基因可能是影响山羊体重性状的主基因或与主基因相连锁,可以用该位点对山羊体重性状进行标记辅助选择。     

小尾寒羊高繁殖力候选基因BMP15和GDF 9的研究

以控制Belclare和Cambridge绵羊高繁殖力的骨形态发生蛋白 15 (bonemorphogeneticprotein 15 ,BMP15 )基因和生长分化因子 9(growthdifferentiationfactor 9,GDF9)基因为候选基因 ,采用PCR RFLP技术检测BMP15基因和GDF9基因在高繁殖力绵羊品种 (小尾寒羊、湖羊 )以及低繁殖力绵羊品种 (多赛特羊、特克塞尔羊、德国肉用美利奴羊 )中的单核苷酸多态性 ,同时研究这两个基因对小尾寒羊高繁殖力的影响。结果表明 :在 5个绵羊品种中都没有检测到GDF9基因的G8突变 (C→T) ,也没有检测到BMP15基因的B4突变 (G→T)。高繁殖力的小尾寒羊在BMP15基因编码序列第 718位碱基处发生了与Belclare绵羊和Cambridge绵羊相同的B2突变 (C→T) ,而其余 4个绵羊品种则没有发生这种突变。对于BMP15基因的B2突变 ,在小尾寒羊中检测到AA、AB两种基因型 ,A等位基因频率为 0 734,B等位基因频率为 0 2 6 6。小尾寒羊与其余 4个绵羊品种间B2突变基因型分布差异极显著 (P <0 0 0 1)。突变杂合基因型 (AB)小尾寒羊平均产羔数比野生纯合基因型 (AA)多 0 6 2只 (P <0 0 1)。研究结果表明 ,BMP15B2突变对小尾寒羊高繁殖力影响作用十分明显 ,同时排除了GDF9G8突变和BMP15B4突变影响小尾寒羊高繁殖力的可能性     

山羊生长激素基因5′调控区的多态性分析

以鲁北白山羊、引进波尔山羊、纯繁波尔山羊以及鲁北白山羊与波尔山羊的杂交一代、回交一代共计274个个体为研究材料,用两对引物分别扩增山羊生长激素(GH)基因5′区的26～239bp以及225～429bp片段,扩增产物经SSCP分析发现均存在多态性.在26～239bp片段上,波尔山羊及杂交后代以AA型个体占多数,而鲁北白山羊则BB型个体较多;在225～429bp片段上,所有种群均以CC型个体较多.对两个片段的纯合型(AA,BB;CC,DD)分别克隆测序发现:(1)26～239bp片段上AA型在第60位发生了C→T的突变,第211位发生碱基C的丢失,(2)225～429bp片段上,DD型存在3处突变,分别为264位由T→C,292位由T→A,372位由C→T.上述结果为首次实验证实山羊生长激素5′调控区存在序列多态性.     

山羊生长激素基因5调控区的多态性分析

以鲁北白山羊、引进波尔山羊、纯繁波尔山羊以及鲁北白山羊与波尔山羊的杂交一代、回交一代共计274个个体为研究材料,用两对引物分别扩增山羊生长激素(GH)基因5^'区的26－239bp以及225－429bp片段,扩增产物经SSCP分析发现均存在多态性。在26－239bp片段上,波尔山羊及杂交后代以 AA型个体占多数,而鲁北白山羊则BB型个体较多;在225－429bp片段上,所有种群均以 CC型个体较多。对两个片段的纯合型（AA,BB;CC,DD）分别克隆测序发现：（1）26－239bp片段上AA型在第60位发生了C→T的突变,第211位发生碱基C的丢失,（2）225－429bp片段上,DD型存在3处突变,分别为264位由T→C,292位由T→A,372位由C→T。上述结果为首次实验证实山羊生长激素5^'调控区存在序列多态性。     

乐至黑山羊PRLR基因外显子10多态性与产羔数的关系研究

设计2对特异性引物对乐至黑山羊PRLR基因第10外显子进行了PCR-SSCP检测,并研究该基因与产子性能的相关性。结果表明,P1引物扩增片段不存在多态性;P2引物扩增片段存在多态性,表现为AA,AB,AD和CD 4种基因型,测序结果表明,4种基因型都在该片段第89、94、146和157位存在C→T、A→C、C→G、G→C的突变;此外AA型还在61位发生C→T的突变;AD型还在175位发生A→G的突变;CD型还在24位发生T→C的突变,96位发生C→T的突变,通过统计分析发现AD型平均产羔数优于其他3种基因型,并且与AB型差异达到显著水平(P<0.05)。因此认为PRLR基因对于乐至黑山羊产子性能有一定的影响。     

山羊抑制素α亚基基因HaeⅡ酶切多态性及其与产羔数性状的关联分析

以323 只马头山羊、努比山羊、波尔山羊和海门山羊为试验材料, 利用PCR-SSCP、PCR-RFLP和克隆测序等方法对山羊抑制素α亚基基因(INHA)编码区一个片段进行多态性分析, 发现在该基因第284位(登录号: L28815)存在G→A突变, 该突变引起HaeⅡ酶切位点改变; HaeⅡ-RFLP分型结果显示G等位基因为优势等位基因, INHA不同基因型平均产羔数表现为GG>AG>AA。结果显示: INHA可能是影响山羊产羔数性状的一个主效基因, G等位基因可能与高产性状呈正相关。     

山羊GOLA-DQB1基因外显子2多态性与免疫性状的相关分析

利用PCR-RFLP技术, 对莱芜黑山羊、鲁波山羊和波尔山羊3个山羊种群共 175 只个体的GOLA-DQB1基因外显子2进行遗传多态性研究, 并对山羊种群的血液免疫指标的效应进行了分析。结果表明: 3个山羊种群共检测到(AA、BB、CC、AB、AC、BC、DD)7种基因型, GOLA-DQB1基因外显子2的第24、151位的碱基表现出多态性。多数指标品种效应是主要效应。莱芜黑山羊中, BC基因型的淋巴细胞百分比(W-SCR)显著高于AC 、CC基因型(P<0.05), 中性球比例(W-LCR)显著低于CC基因型(P<0.05), 大型白细胞数(W-LCC)低于AC、CC基因型, 但差异不显著(P>0.05)。波尔山羊中, BC基因型的W-LCC低于AA 、AB 、BB基因型, 但差异不显著(P>0.05)。鲁波山羊中, BC 、AC基因型的W-LCC显著低于AA基因型(P<0.05)。揭示GOLA-DQB1基因与血液免疫性状有一定的相关性。     

Analysis on DNA sequence of goat RFRP gene and its possible association with average daily sunshine duration

Goat RFRP gene was cloned and its mutations were detected in thirteen goat breeds whose reproductive seasonality and litter size were different. Then sequence characteristics were analyzed and association analyses were performed to reveal the relationships between mutations of RFRP gene and average daily sunshine duration, reproductive seasonality as well as litter size in goats. A 4,862 bp DNA fragment of goat RFRP gene was obtained and the complete CDS of 591 bp encodes 196 amino acids, having high homology with that of other mammals. The protein was predicted to be a secreted protein with a signal peptide of 21 amino acids. Moreover, two mutations (A712G, T1493C) in 5' regulatory region and one mutation (A3438T) in exon 2 were detected. The test of genotype distribution in six selective goat breeds showed that there was no uniform significant association between the three polymorphisms and seasonal reproduction. The association just existed in some goat breeds for each locus. Interestingly, however, there was a strong positive correlation (r = 0.830, P = 0.003) between the G allele frequency of the A712G locus and average daily sunshine duration in ten local goat breeds, suggesting that RFRP gene has undergone a selective pressure in sunshine duration and may have indirect relationship with reproductive seasonality in goats. Additionally, no significant difference was found in litter size between genotypes in prolific Jining Grey goats.     

Nucleotide sequencing and DNA polymorphism studies of BMP 15 gene in Corriedale and local Kashmir valley sheep (Ovis aries)

The families of TGF-β proteins are the most important growth factors in the ovary for growth and differentiation of early ovarian follicles. Three related oocyte-derived members of the transforming growth factor-β superfamily, namely GDF9, BMP15 and BMPR-IB have been shown to be essential for follicular growth and ovulation. The objective of the present study was to detect the incidence of mutation in intronic portion of BMP 15 gene in Corriedale and local Kashmir valley sheep breeds. Blood samples were collected from 85 ewes and genomic DNA was extracted using the modified phenol chloroform method. The quantity and quality of extracted DNA was examined using spectrophotometry and gel electrophoresis, respectively. A fragment with the size of 356 bp was amplified using polymerase chain reaction (PCR) with a pair of specific primers. The amplified PCR products were digested with Mph11031 restriction enzyme. In the presence of mutation at this locus, the Mph11031 enzyme cannot recognize the restriction site. However, here in the absence of mutations, the enzyme recognizes one restriction site and divides the amplified fragment into two fragments of 152 and 204 bp. The 356 bp fragment was also analyzed for polymorphism by SSCP technique. The results indicated two different banding patterns AA and AB for this fragment. Later on two different allelic forms A and B were confirmed by nucleotide sequencing. The 356 bp nucleotide sequence was subjected to alignment analysis and it was observed that sequence similarity of this fragment with that of other sheep and Jining grey goat was more than 97.8%. Phylogenetic analysis revealed that both designated A and B alleles as well as published sequence of sheep form a common cluster indicating their evolutionary closeness. The origin of Jining grey goat was located some distance away from the sheep. The overall frequencies of AA and AB genotypes were 0.79 and 0.21. The breed wise frequencies were 0.78 and 0.22 in Corriedale sheep and the frequencies in Kashmir valley sheep were 0.80 and 0.20 for AA and AB genotypes, respectively. The overall allelic frequencies of A and B alleles were 0.89 and 0.11 whereas allelic frequencies Corriedale sheep was 0.89 and 0.11 and that of Kashmir valley sheep were 0.90 and 0.10.     

Polymorphism of bone morphogenetic protein 4 gene and its relationship with litter size of Jining Grey goats

Two pairs of primers (P1 and P2) were designed to detect single nucleotide polymorphisms of exon 2 and intron 2 of bone morphogenetic protein 4 (BMP4) gene in both high fecundity breed (Jining Grey goat) and low fecundity breeds (Boer, Angora and Inner Mongolia Cashmere goats) by single strand conformation polymorphism. Results showed that no polymorphism was detected for exon 2 (primer P1) of BMP4 gene in four goat breeds. For intron 2 (primer P2), three genotypes (AA, AB and BB) were detected in Jining Grey and Inner Mongolia Cashmere goats, two genotypes (AB and BB) in Angora goats, and only one genotype (AA) in Boer goats. Sequencing revealed one mutation (2203G>A) of BMP4 gene in the genotype BB in comparison to the genotype AA. The differences of litter size between AA, AB and BB genotypes were not significant (P > 0.05) in Jining Grey goats. A pair of primer (P3) was designed to detect polymorphism in the 3' flanking region of BMP4 gene that contained dinucleotide repeated sequence (CA) in the four goat breeds by microsatellite analysis. For primer P3, three genotypes (CC, CD and DD) were detected in four goat breeds. Sequencing revealed one more CA dinucleotide in genotype DD than in genotype CC. The Jining Grey does with genotype CC had 0.55 (P < 0.05) or 0.72 (P < 0.05) kids more than those with genotype CD or DD. These results preliminarily indicated that allele C of BMP4 gene is a potential DNA marker for improving litter size in goats.     

None of polymorphism of ovine fecundity major genes FecB and FecX was tested in goat

The polymorphism of mutation Q249R in BMPR-IB gene (FecB) and loci FecX^I, FecX^H, FecX^G, FecX^B in BMP15 gene was analyzed by forced PCR-RFLP method in 550 individuals from 6 flocks or breeds of goats with litter size varied from 1.4 to 2.7 including Boer (209), Haimen (128), second generation of Boer goat crossed with Huanghuai goat (82), Huanghuai (71), Nubi (37) and Matou (23) goat. None of mutations was detected in these goat breeds and their crossbreed. These results suggest that fecundity of goat is not linked to the same loci in BMPR-IB and BMP15 as sheep. Therefore, it is necessary to seek for other genes or loci in order to develop marker assistance selection techniques and study the prolific mechanism of the goat.     

Polymorphisms of caprine <Emphasis Type="Italic">POU1F1</Emphasis> gene and their association with litter size in Jining Grey goats

Seven pairs of primers were designed to amplify 5′ promoter region, six exons and partial introns and to detect the polymorphisms of POU1F1 gene in five goat breeds with different prolificacy. The results showed that six mutations were identified in caprine POU1F1 gene including C256T in exon 3, C53T and T123G in intron 3, and G682T (A228S), T723G and C837T in exon 6. The former four mutations were novel SNPs in goat POU1F1 gene. The 53 and 123 loci were in complete linkage disequilibrium in five caprine breeds. Regarding the 256 locus, the Jining Grey goat does with genotype CT had 0.66 kids more than those with genotype CC (P < 0.05), while does with genotype GT had 0.63 (P < 0.05) kids more than those with genotype GG at the 682 locus. The present study preliminarily showed an association between allele T at the 256 and 682 loci of POU1F1 gene and high litter size in Jining Grey goats. Totally 16 haplotypes and 50 genotypes were identified at the above six loci in POU1F1 gene of five goat breeds. Three common haplotypes (hap2, hap3 and hap4) were identified in five goat breeds joined. Four specific haplotypes (hap7, hap9, hap11 and hap13) were detected in Jining Grey goats. The predominant haplotype was hap1 (35.29% and 48.25%) in both Jining Grey and Guizhou White goats, while hap4 (50%) in Boer goats, and hap2 (42.86% and 38.75%) in both Wendeng Dairy and Liaoning Cashmere goats. The most frequent genotypes at six loci in the above five goat breeds were hap1/hap2 (14.38%) and hap1/hap4 (14.38%), hap1/hap2 (38.60%), hap4/hap4 (40.91%), hap2/hap4 (26.53%), hap2/hap5 (20.00%), respectively. The Jining Grey goat does with nine genotypes analyzed of POU1F1 gene showed no obvious difference in litter size.     

体细胞克隆山羊微卫星DNA分析

用10对山羊微卫星DNA多态性引物对2只体细胞克隆济宁青山羊、青山羊供体细胞、受体奶山羊母羊以及具有亲缘关系的3只对照济宁青山羊进行微卫星DNA分析.结果表明有5对山羊微卫星DNA多态性引物,即SR-CRSP1, SR-CRSP5, SR-CRSP6, SR-CRSP7和SR-CRSP24,扩增产物有明显的多态性.扩增产物经过6%变性聚丙烯酰胺凝胶电泳后银染,结果2只体细胞克隆山羊的微卫星DNA指纹与供体细胞完全相同,而且不同于其受体母亲也不同于其他所有同品种不同个体的对照青山羊.证明体细胞克隆山羊基因组来源于供体细胞.     

Analysis on DNA sequence of <Emphasis Type="Italic">GPR54</Emphasis> gene and its association with litter size in goats

The kisspeptin/GPR54 pathway is crucial in the process of puberty onset. Six pairs of primers were designed to clone goat GPR54 and scan polymorphisms and one pair of primers to detect polymorphisms of GPR54 in sexual precocious and sexual late-maturing goat breeds. A DNA fragment of 4258 bp of goat GPR54 was obtained, which contains an open reading frame (ORF) of 1137 bp and encodes 378 amino acids, having a high homology with other mammals. The protein was predicted to have seven transmembrane regions. There were no base pair variation in exons 1–4 and three base changes (G4014A, G4136A and C4152T) in exon 5 by sequencing and the three mutations may have some correlation with sexual precocity in goats. For the 4152 locus, the Jining Grey goat does with genotype TT and CT had 1.02 and 0.84 (P < 0.01) kids more than those with genotype CC, respectively. No significant difference (P > 0.05) was found in litter size between TT and CT genotypes in Jining Grey goat. For the other two loci, no significant difference (P > 0.05) was found in litter size between different genotypes in Jining Grey goats. The present study preliminarily indicated an association between allele T of the 4152 locus in GPR54 and high litter size in Jining Grey goats.     

Polymorphisms of BMPR-IB gene and their relationship with litter size in goats

The bone morphogenetic protein receptor IB (BMPR-IB) gene was studied as a candidate gene for the prolificacy of goats. According to mRNA sequence of ovine BMPR-IB gene, ten pairs of primers were designed to detect single nucleotide polymorphisms (SNPs) of exon 1, exon 2, exon 6 to exon 10 and 3′ untranslated region (UTR) of the BMPR-IB gene in both high prolificacy breed (Jining Grey goat) and low prolificacy breeds (Wendeng Dairy and Inner Mongolia Cashmere goats) by polymerase chain reaction (PCR)-single strand conformation polymorphism (SSCP) method. Only the products amplified by primers P8 and P10 of the 3′UTR displayed polymorphisms. For primer P8, three genotypes (AA, AB and BB) were detected in Jining Grey and Wendeng Dairy goats, two genotypes (AA and AB) were in Inner Mongolia Cashmere goats. Sequencing revealed one mutation (71C→T) of the BMPR-IB gene in genotype BB compared with AA. The differences of least squares mean (LSM) for litter size between genotypes AA, AB and BB were non-significant (P > 0.05) in Jining Grey goats. For primer P10, three genotypes (CC, CD and DD) were detected in Jining Grey and Wendeng Dairy goats and one genotype (CC) in Inner Mongolia Cashmere goats. Sequencing revealed one mutation (130T→C) of the BMPR-IB gene in genotype DD compared with CC. The differences of LSM for litter size between genotypes CC, CD and DD were non-significant (P > 0.05) in Jining Grey goats. These results preliminarily showed that the detected loci of the BMPR-IB gene had no significant effect on prolificacy of Jining Grey goats.     

Polymorphism of fecundity genes (BMPR1B, BMP15 and GDF9) in the Indian prolific Black Bengal goat

The Black Bengal is a prolific goat breed in India. Natural mutations in prolific sheep breeds have shown that the transforming growth factor beta (TGF-β) super family ligands such as growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and their type I receptor (bone morphogenetic protein receptor, BMPR1B) are crucial for ovulation and as well as for increasing litter size. Mutations in any of these genes increased prolificacy in sheep. Based on the known mutation information in sheep PCR primers were designed to screen known polymorphism in 88 random Black Bengal goats. Only the BMPR1B gene was polymorphic. Three genotypes of animals were detected in tested animals with mutant (FecB^B) and wild type (FecB⁺) alleles were 0.57 and 0.43, respectively. Non-carrier, heterozygous carrier and homozygous carrier Black Bengal does had 2.7, 3.04 and 3.11 kids, respectively. All known point mutations of BMP15 and GDF9 genes were monomorphic in the animals tested. These results preliminarily showed that the BMPR1B gene might be a major gene that influences prolificacy of Black Bengal goats.