Genome-wide transcriptome analysis in the ovaries of two goats identifies differentially expressed genes related to fecundity

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.

Growth differentiation factor 9 (GDF9) was studied as a candidate gene for high prolificacy in goats. The polymorphism of exon 1 and flanking of GDF9 gene was detected by PCR-SSCP in five goat breeds with different prolificacy. Three genotypes (AA, AB and BB) were detected in goat breeds joined and two silent mutations (c.183A>C and c.336C>T) were identified in comparison genotype AA with genotype BB. Heterozygous genotype AB and wild type BB were detected in all five goat breeds and homozygous genotype AA was only detected in Jining Grey goats. The frequencies of genotypes AA, AB and BB were 0.18, 0.42 and 0.40 in Jining Grey goats, respectively. The genotype distribution was different (P < 0.01) between high prolificacy breed (Jining Grey goat) and low prolificacy breeds (Boer, Wendeng Dairy, Liaoning Cashmere and Beijing native goats). The Jining Grey goat does with genotype AA and AB had 0.72 (P < 0.01) and 0.56 (P < 0.01) kids more than those with genotype BB, respectively. The does with genotype AA had 0.16 (P > 0.05) kids more than those with genotype AB. These indicated that the allele A may have certain correlation with prolificacy in Jining Grey goats.

Abstract. Lin28a and Lin28b serve as crucial regulators in gametogenesis. However, it remains unclear whether the two genes are expressed in goat ovarian tissues. In the present study, we examined the distribution of Lin28a proteins in the ovarian tissues of juvenile and pubertal Jining Grey (JG) goats. Additionally, we investigated the mutations in Lin28a and Lin28b and their association with litter size in JG goats. Lin28a proteins were predominantly expressed in the oocytes and granulosa cells of the follicles at different stages of development. Lin28a was also expressed in the corpus luteum, and its expression level declined as the corpus luteum formed and degenerated into the corpus albicans. A mutation (12298T&gt;C) was identified in intron 2 of Lin28a, and there was no significant difference (P&gt;0.05) in terms of the kidding number among the three genotypes of JG goats. For Lin28b, a mutation (718A&gt;T) in exon 4 was found, presenting three genotypes. The frequency of genotypes AA, AT, and TT was 0.466, 0.504, and 0.03, respectively. The A allele was the dominant allele in JG goats, with an allele frequency of 0.718. At this locus, JG female goats with genotypes TT and AT had 0.83 (P&lt;0.01) and 0.48 (P&lt;0.05) more kids than the goats with genotype AA, and no significant difference (P&gt;0.05) was observed in terms of litter size between the TT and AT genotypes. The present study preliminarily indicated an association between the T allele at the 718 locus in Lin28b and a larger litter size in JG goats. Moreover, Lin28a was widely expressed in follicles. These results suggest a potential DNA marker for improving the kidding number in goats and hold significance for goat breeding.

The exons 1, 2 and flanking region of growth differentiation factor 9 (GDF9) gene in five randomly selected does of Jining Grey, Boer and Liaoning Cashmere goats were amplified and analyzed. Thirteen nucleotide differences were identified in GDF9 gene between sheep (AF078545) and goats. Four SNPs (G3288A in intron 1, G423A, A959C [Gln320Pro] and G1189A [Val397Ile] in exon 2) were detected in four goat breeds with different prolificacy, in which G3288A was a new SNP in goats. The results showed that loci 3288, 423 and 1189 in Boer goats, loci 3288 and 423 in Guizhou White goats, loci 423 and 1189 in Liaoning Cashmere goats were all in complete linkage disequilibrium (D' = 1, r (2) = 1), respectively. In moderate (Boer goat) and low prolificacy (Liaoning Cashmere goat) breeds, linkage analysis indicated that there were more fervent linkage disequilibrium among loci 3288, 423 and 1189 than high prolificacy (Jining Grey and Guizhou White goats) breeds. For the 959 locus, the genotype distribution showed obvious difference between high prolificacy breeds and moderate or low prolificacy breeds (P < 0.05 or P < 0.01). The Jining Grey goat does with genotype CC or AC had 0.81 (P < 0.01) or 0.63 (P < 0.01) kids more than those with genotype AA, respectively. The present study preliminarily showed an association between allele C at 959 locus of GDF9 gene and high litter size in Jining Grey goats. These results provide further evidence that the GDF9 gene may be significantly correlated with high prolificacy in goats.

Gene cloning, homology comparison and analysis of the main functional structure domains of beta estrogen receptor in Jining Gray goat

The Jining Grey is a prolific local goat breed in P.R. China. Bone morphogenetic protein 15 (BMP15) gene that controls high fecundity of Inverdale, Hanna, Lacaune, Belclare, Cambridge, and Small Tailed Han ewes was studied as a candidate gene for the prolificacy of Jining Grey goats. According to the sequence of ovine BMP15 gene, six pairs of primers were designed to detect single nucleotide polymorphisms in exon 1 and exon 2 of the BMP15 gene in both high fecundity breed (Jining Grey goats) and low fecundity breeds (Boer, Liaoning Cashmere, and Inner Mongolia Cashmere goats) by single strand conformation polymorphism (SSCP). Two pairs of primers (F1/R1 and F2/R2) were used to amplify the exon 1. Four pairs of primers (F3/R3, F4/R4, F5/R5, and F6/R6) were used to amplify the exon 2. Only the products amplified by primer F5/R5 displayed polymorphism. Results indicated that two genotypes (AA and AB) were detected in prolific Jining Grey goats and only one genotype (AA) was detected in low fecundity goat breeds. In Jining Grey goats frequencies of genotypes AA and AB were 0.10 and 0.90, respectively. Sequencing revealed two point mutations (G963A and G1050C) of BMP15 gene in the AB genotype in comparison to the AA genotype. In Jining Grey goats the heterozygous AB does had 1.13 (p < 0.01) kids more than the homozygous AA does. These results preliminarily showed that the BMP15 gene is either a major gene that influences the prolificacy of Jining Grey goats or a molecular genetic marker in close linkage with such a gene.

Prolificacy of most local goat breeds in China is low. Jining Grey goat is one of the most prolific goat breeds in China, it is an important goat breed for the rural economy. ASMT (acetylserotonin O‐methyltransferase) and ADAMTS1 (ADAM metallopeptidase with thrombospondin type 1 motif) are essential for animal reproduction. Single nucleotide polymorphisms (SNPs) of ASMT and ADAMTS1 genes in the highly prolific breed (Jining Grey goats), medium prolific breed (Boer goats and Guizhou White goats) and low prolific breeds (Angora goats, Liaoning Cashmere goats and Inner Mongolia Cashmere goats) were detected by polymerase chain reaction‐restriction fragment length polymorphism and sequencing. Two SNPs (g.158122T>C, g.158700G>A) of ASMT gene and two SNPs (g.7979798A>G, g.7979477C>T) of ADAMTS1 gene were identified. For g.158122T>C of ASMT gene, further analysis revealed that genotype TC or CC had 0.66 (p < 0.05) or 0.75 (p < 0.05) kids more than those with genotype TT in Jining Grey goats. No significant difference (p > 0.05) was found in litter size between TC and CC genotypes. The SNP (g.158122T>C) caused a p.Tyr298His change and this SNP mutation resulted in changes in protein binding sites and macromolecule‐binding sites. The improvement in reproductive performance may be due to changes in the structure of ASMT protein. For g.7979477C>T of ADAMTS1 gene, Jining Grey does with genotype CT or TT had 0.82 (p < 0.05) or 0.86 (p < 0.05) more kids than those with genotype CC. No significant difference (p > 0.05) was found in litter size between CT or TT genotypes. These results preliminarily indicated that C allele (g.158122T>C) of ASMT gene and T allele (g.7979477C>T) of ADAMTS1 gene are potential molecular markers which could improve litter size of Jining Grey goats and be used in goat breeding.

Genome-wide detection of selective signals for fecundity traits in goats (Capra hircus).

In our early work, one mutation (A898G, which caused Ser32Gly (S32G) change of mature peptide) was identified in the goat BMP15 gene and the allele G was associated with high litter size in Jining Grey goats (the goat breed with the highest fecundity in China). The aim of this research was to investigate the genetic structure of the prolific BMP15 gene in 18 goat breeds reared in China, including 12 native breeds, 3 introduced breeds, and 2 cultivated breeds. Genomic DNA of 1011 goats was screened for the S32G mutation. The results showed that this mutation existed in all 3 cultivated and 8 native goat breeds but in none of the 3 introduced goat breeds, which hinted that this mutation may originate from native goat breeds of China. Besides the Jining Grey goats, the BMP15 gene was also a potential prolificacy gene in Matou goats. Moreover, the structure prediction indicated that the S32G mutation might participate in the binding of BMP15 with receptors.

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.

Single nucleotide polymorphisms of the 5' regulatory region of the FSHR gene were detected in 3 goat breeds (Jining Grey, Boer, and Inner Mongolia Cashmere) by the polymerase chain reaction-single-strand conformation polymorphism method, and its effects on litter size in Jining Grey goats were also evaluated. Concerning primer P1, 2 genotypes (AA and BB) were detected in Boer goats, while only the AA genotype was identified in Jining Grey and Inner Mongolia Cashmere goats. Sequencing revealed a T->A transversion at the 26th position and an A->C transversion at the 61st position of the amplified region in genotype BB in comparison to genotype AA. Regarding primer P2, 3 genotypes (CC, CD, and DD) were identified in the 3 goat breeds mentioned above. Sequencing revealed a T->A transversion at the 70th position and a G->C transversion at the 130th position of the amplified region in genotype DD compared with CC. Does with genotype CC had 0.46 (P < 0.05) or 1.03 (P < 0.05) kids more than those with genotype CD or DD, respectively, while does with CD had 0.57 (P < 0.05) kids more than those with DD in Jining Grey goats. These results preliminarily showed that allele C of the FSHR gene was a potential marker for improving litter size in goats.

Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair-wise genetic differentiation (FST ) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within-breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.

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.

Three pairs of primers were designed to clone the goat KiSS-1 and scan polymorphisms and four pairs to detect polymorphisms in sexual precocious and sexual late-maturing goat breeds. A 4118 bp DNA fragment was obtained, which contains an ORF of 408 bp and encodes 135 amino acids, having a high homology with other mammals. The protein was predicted containing a signal peptide of 17 amino acids. There are two mutations (G3433A [A86T] and C3688A) in exon 3, three mutations (G296C, G454T and T505A) in intron 1 and a 18 bp deletion (-)/insertion (+) (1960-1977) in intron 2 and no mutations in exon 2. The genotype distribution didn't show obvious difference between sexual precocious and sexual late-maturing goat breeds and no consistency within the sexual late-maturing breeds. For the 296 locus, the Jining Grey goats with genotype CC had 0.80 (P < 0.01) or 0.77 (P < 0.01) kids more than those with genotype GG or GC, respectively. No significant difference (P > 0.05) was found in litter size between GG and GC. For the 1960-1977 locus, the Jining Grey goat does with genotype -/- had 0.77 (P < 0.01) or 0.73 (P < 0.01) kids more than those with +/+ or +/-, respectively. No significant difference (P > 0.05) was found in litter size between +/+ and +/- genotypes. For the other four 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 C of the 296 locus and allele (-) of the 1960-1977 locus in KiSS-1 and high litter size in Jining Grey goats.

Litter size is a favorable economic trait for the goat industry, but remains a complex trait controlled by multiple genes in multiple organs. Several genes have been identified that may affect embryo survival, follicular development, and the health of fetuses during pregnancy. Jining Grey goats demonstrate the largest litter size among goat breeds indigenous to China. In order to better understand the genetic basis of this trait, six suppression subtractive hybridization (SSH) cDNA libraries were constructed using pooled mRNAs from hypothalamuses, pituitaries, and ovaries of sexually mature and adult polytocous Jining Grey goats, as testers, versus the pooled corresponding mRNAs of monotocous Liaoning Cashmere goats, as drivers. A total of 1,458 true-positive clones--including 955 known genes and 481 known and 22 unknown expressed sequence tags--were obtained from the SSH libraries by sequencing and alignment. The known genes were categorized into cellular processes and signaling information storage and processing, and metabolism. Three genes (FTH1, GH, and SAA) were selected to validate the SSH results by quantitative real-time PCR; all three were up-regulated in the corresponding tissues in the tester group indicating that these are candidate genes associated with the large litter size of Jining Grey goats. Several other identified genes may affect embryo survival, follicular development, and health during pregnancy. This study provides insights into the mechanistic basis by which the caprine hypothalamic-pituitary-gonadal axis affects reproductive traits and provides a theoretical basis for goat production and breeding.