Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

A novel SNP of α-lactalbumin gene in Chinese dairy goats

The α-lactalbumin (α-LA) plays a key role in lactose synthesis in mammary glands of domestic animals. Mutations in the α-LA gene are associated with the milk traits in dairy cattle. In our study, a novel SNP: NO_X06366: g.875 C > T was detected in 708 dairy goat individuals—268 of the Xinong Saanen breed and 440 of Guanzhong breed, which revealed a synonymous mutation in the exon 1 of α-LA gene. The Polymerase Chain Reaction-Single Strand Conformation Polymorphism (PCR-SSCP) and sequencing techniques showed that there were three genotypes: CC, CT and TT. Moreover, the χ²-test showed that the genotype frequencies of the two breeds were in good agreement with the Hardy-Weinberg equilibrium (P > 0.05). The relationship of the polymorphism of dairy goat α-LA gene with the milk trait and the body size trait was revealed. Individuals with the CC genotype were significantly smaller at chest circumference than those with CT (P < 0.05) in both breeds. But the milk trait and other body size traits of the two dairy goat breeds had no significant association with genotypes studied.     

Effects of genetic variability of the dairy goat growth hormone releasing hormone receptor (<Emphasis Type="Italic">GHRHR</Emphasis>) gene on growth traits

Growth hormone-releasing hormone receptor (GHRHR) plays a critical role in growth hormone (GH) synthesis, release and regulation of pituitary somatotroph expansion in vertebrates. The objective of this study was to investigate variations in goat GHRHR gene and their associations with growth traits in 668 dairy goats. The results showed four novel single nucleotide polymorphisms (SNPs): NC_007302:g.5203C>T, 7307C>G, 9583G>A and 9668A>C. In detail, the novel SNP C>T in the 5203rd nucleotide identified a missense mutation: CCC (Pro)>TCC (Phe) at position 116aa of the goat GHRHR (423aa). Besides, 9583G>A and 9668A>C polymorphism were in complete linkage disequilibrium. The genetic diversity analysis revealed that the Guanzhong dairy goat possessed intermediate genetic diversity in P3 and P7 loci, and the Xinong Sannen dairy goat belonged to poor genetic diversity in P4 locus. Significant associations between the genotypes of P3 locus and body length, body height and chest circumference was observed in Guanzhong goat (P < 0.05). However, in Xinong saanen population, significant statistical difference was only found in body height and body length (P < 0.05). In P4 and P7 loci, no significant associations were detected between any variant sites and body length, body height and chest circumference, as well as for the milk traits (P > 0.05). These results strongly suggested that the goat GHRHR gene is a candidate gene that influences growth traits in dairy goat.     

Two mutations in the 5′‐flanking region of the KITLG gene are associated with litter size of dairy goats

In this study, Xinong Saanen (SN) and Guanzhong (GZ) dairy goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the 5′‐flanking region of the KITLG gene by DNA sequencing and primer‐introduced restriction analysis–polymerase chain reaction. Two novel SNPs (g.13090G>T and g.13664C>A) were identified (GenBank Accession no. KM658964). Furthermore, g.13090G>T and g.13664C>A loci were closely linked in SN and GZ breeds (r² > 0.33). Association analysis results showed that g.13090G>T and g.13664C>A SNPs significantly affected litter size (P < 0.05). The litter size of individuals with the combined genotype GG/CC from both dairy goat breeds was greater than that of individuals with TT/AA in average parity (P < 0.05). Known biochemical and physiological functions, along with our results, indicated that GG/CC could be used in marker‐assisted selection to choose individuals with greater litter size from both breeds. These results extend the spectrum of genetic variation in the caprine KITLG gene and may contribute to genetic resources and breeding of goats.     

Polymorphisms of myostatin gene (<Emphasis Type="Italic">MSTN</Emphasis>) in four goat breeds and their effects on Boer goat growth performance

Polymorphisms of myostatin (MSTN) gene were investigated as a candidate marker for goat growth in 687 individuals by gene sequencing and polymerase chain reaction restriction fragment length polymorphism methods. Three potential genotypes (AA, AB and BB) of 5 bp indel (1,256 TTTTA/−) in 5′UTR were detected in four breeds. The polymorphism (CC, CD and DD) of substitution (1,388 T/A) in exon 1 was only segregating in Boer. Genotype AB resulted in significant increases in body weights at birth (BW0), 90 (BW90) and 300 (BW300) days of age, and birth body length compared to BB (P < 0.05). Whilst genotype CD contributed to heavier BW0 and BW90, and larger birth body height (BH) compared to CC (P < 0.05). Individuals either with AB or CD genotype also had greater values in other body sizes, although no significant differences appeared (P > 0.05). When in combination, the combined genotype ABCD displayed positive impacts on better growth traits in BW0, BW90, BW300, BH and BCG (P < 0.05). These results indicate that these markers in myostatin (MSTN) are associated with Boer growth and may be useful for marker assisted selection.     

GnRHR gene polymorphisms and their effects on reproductive performance in Chinese goats

In this study, polymorphisms in the goat GnRHR gene exon 1 were detected by PCR-SSCP and DNA sequencing methods in 786 individuals from two different goat breeds. Two haplotypes (A and B), two observed genotypes (AA and AB), and two single nucleotide polymorphisms (SNPs) were detected, which resulted in five amino acid substitutions. The frequencies of haplotypes A and B in the two goat breeds were 0.78–0.83 and 0.17–0.22, respectively. The SNP locus was in Hardy–Weinberg disequilibrium in the two goat breeds (P < 0.05). Polymorphisms of the GnRHR gene were shown to be associated with litter size in the two goat breeds. The SNPs in the goat GnRHR gene had significant effects on litter size (P < 0.05). Therefore, these results suggest that the GnRHR gene is a strong candidate gene that affects litter size in goat.     

Single nucleotide polymorphisms of the prolactin receptor (PRLR) gene and its association with growth traits in chinese cattle

Genetic polymorphism of the prolactin receptor (PRLR) gene was detected by PCR-SSCP and DNA sequencing methods in 665 individuals from five Chinese cattle breeds. The results showed that at the P1 locus, three observed genotypes (AA, AB and BB), two linked SNPs (G1267A and T1268C), and one missense mutation (S18N) within a putative signal peptide were determined. The frequencies of haplotypes A and B in the five breeds were 0.596–0.802 and 0.198–0.404, respectively. Polymorphism of the PRLR gene was shown to be significantly associated with growth traits in the Nanyang breed. Individuals with genotype BB had greater hucklebone width, body weight and average daily gain than those with genotype AA at 6 months old (P < 0.01), as well as better body height, body length and heart girth when 6 months (P < 0.05). This study revealed for the first time that the PRLR gene is a promising candidate gene that affects growth traits in cattle.     

Single-nucleotide polymorphisms g.151435C>T and g.173057T>C in PRLR gene regulated by bta-miR-302a are associated with litter size in goats

Single-nucleotide polymorphisms (SNPs) located at microRNA-binding sites (miR-SNPs) can affect the expression of genes. This study aimed to identify the miR-SNPs associated with litter size. Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect SNPs in the caprine prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction, and polymerase chain reaction-restriction fragment length polymorphism. Three novel SNPs (g.151435C>T, g.151454A>G, and g.173057T>C) were identified in the caprine PRLR gene. Statistical results indicated that the g.151435C>T and g.173057T>C SNPs were significantly associated with litter size in Guanzhong and Boer goat breeds. Further analysis revealed that combinative genotype C6 (TTAACC) was better than the others for litter size in both goat breeds. Furthermore, the PRLR g.173057T>C polymorphism was predicted to regulate the binding activity of bta-miR-302a. Luciferase reporter gene assay confirmed that 173057C to T substitution disrupted the binding site for bta-miR-302a, resulting in the reduced levels of luciferase. Taken together, these findings suggested that bta-miR-302a can influence the expression of PRLR protein by binding with 3′untranslated region, resulting in that the g.173057T>C SNP had significant effects on litter size.     

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

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 ² = 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.     

Analysis of genetic variability at codon 42 within caprine <Emphasis Type="Italic">prion protein</Emphasis> gene in relation to production traits in Chinese domestic breeds

In this study, genetic variability at codon 42 within prion protein (PRNP) gene and its associations with production traits were investigated in 2002 goats from four Chinese domestic breeds. The frequencies of allele “A” ranged from 0.353 to 0.562 in analyzed goat breeds with Hardy–Weinberg equilibrium (P > 0.05) except Xinong Sannen (XNSN) dairy breed. The establishment of relationships between different genotypes and growth traits was performed in Inner Mongolia white Cashmere (IMWC) breed and revealed an association of the polymorphism with body weight at 7-year-old goats (P = 0.033). The individuals with genotype GG showed heavier body weight than those with genotype AA. Moreover, association analysis detected two significant associations between different genotypes and cashmere yield and fiber length in IMWC breed (P = 0.009, P = 0.048, respectively). In addition, three significant associations of different genotypes with density of milk (a.m. and p.m.), solids-not-fat of milk (P = 0.013, P = 0.009 and P = 0.002), respectively, were found in XNSN breed. Genotype GG had better milk quality than others. These findings suggested that the polymorphism of codon 42 within PRNP was a useful DNA marker for eliminating or selecting excellent individuals in relation to production traits in marker-assist selection breeding of goat.     

Polymorphisms in BMP-2 gene and their associations with growth traits in goats

BMPs plays an important role in skeletal development. In this study, PCR-SSCP and DNA sequencing methods were employed to screen the genetic variation of caprine BMP-2 gene in 299 goats from three breeds (Boer goat, Chinese Xuhuai white goat and Chinese Haimen goat). Three fragments of BMP-2 gene were investigated, only 3′ flanking region of BMP-2 gene showed polymorphism. The alignment between nucleotide sequences of AY714781.1 in GenBank and the sequencing results of three PCR products with different patterns revealed that there was one mutation (AY714781.1: g.A172G) in 3′ flanking region of BMP-2 gene, which constructed three genotypes (AA, AB, BB). The frequencies of allele A and genotype AA were dominant in all three breeds. Frequency of allele B in Haimen goat breed was significantly lower than that of the other two breeds. The genotype distributions were in good agreement with Hardy-Weinberg equilibrium (P>0.05) in each breed. The PIC (Polymorphism Information Content) values in three populations were 0.3261, 0.2558, 0.1663 for Boer goat, Xuhuai White goat and Haimen goat respectively. It was indicated that individuals of Boer goats with genotype BB were significantly higher than those of individuals with genotype AA in body trunk index (P<0.01). No polymorphism was detected in the exon3.     

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.     

Association of goat (<Emphasis Type="Italic">Capra hircus</Emphasis>) CD4 gene exon 6 polymorphisms with ability of sperm internalizing exogenous DNA

As one of the transport systems on the sperm plasma membrane, CD4 molecule plays a distinct role in the process of sperm/DNA interaction. This makes it possible to explain the mechanism of sperm-mediated gene transfer (SMGT), which at present is still a mystery in this area. In this study, seminal samples of 60 individuals from seven breed bucks were collected to detect the ability of sperm in internalizing exogenous DNA, and genomic DNA from 147 individual blood samples (including 60 bucks referred ahead) were extracted to test the polymorphisms of CD4 genes by using PCR-SSCP technique. Then the correlation between them was evaluated. The results showed that: (1) it was a novel finding that breed-dependence of exogenous DNA binding to goat spermatozoa. There was the most significant difference among the buck breeds of sperm in binding exogenous DNA (F_{(6, 53)} = 4.811, P = 0.001) and in internalizing them into nuclei (F_{(6, 53)} = 4.587, P = 0.001). The ability of Lezhi Black goat was the lowest (P < 0.01) among the seven breeds. (2) There was no significant correlation between the ability of sperm in internalizing exogenous DNA and each semen quality parameter (P > 0.05). (3) In particular, three single-nucleotide polymorphisms (SNP) were described and there was one SNP (G/A⁷⁰⁰) of CD4 gene that made G234R substitution in the amino acid sequence of CD4 molecule. Nanjiang Yellow goat and Lezhi Black goat had higher hereditary variation compared with other breeds. (4) CD4 polymorphisms were highly associated with the ability of sperm in internalizing exogenous DNA. The SNP of Caprine CD4 gene exon 6 might be an important molecular marker of the ability to internalize exogenous DNA into sperm.     

Identification of a functional SNP in the 3′‐UTR of caprine MTHFR gene that is associated with milk protein levels

Xinong Saanen (n = 305) and Guanzhong (n = 317) dairy goats were used to detect SNPs in the caprine MTHFR 3′‐UTR by DNA sequencing. One novel SNP (c.*2494G>A) was identified in the said region. Individuals with the AA genotype had greater milk protein levels than did those with the GG genotype at the c.*2494 G>A locus in both dairy goat breeds (P < 0.05). Functional assays indicated that the MTHFR:c.2494G>A substitution could increase the binding activity of bta‐miR‐370 with the MTHFR 3′‐UTR. In addition, we observed a significant increase in the MTHFR protein level of AA carriers relative to that of GG carriers. These altered levels of MTHFR protein may account for the association of the SNP with milk protein level.     

Effects of genetic variability of the carpine homeobox transcription factor <Emphasis Type="Italic">HESX1</Emphasis> gene on performance traits

HESX1 plays a key role in the development of the forebrain and pituitary gland and produces potential effects on performance traits. The objective of this study was to detect and assess the associations of the possible polymorphisms of six loci within HESX1 gene with performance traits in Chinese 1,119 goats. Only one novel SNP (NM_001494116:g.307049A > G) locating on IVS1 + 348A > G was identified and detected by HaeIII forced-RFLP-PCR. The frequencies of allele “G” varied from 0.025 to 0.245 in analyzed populations with the Hardy–Weinberg equilibrium (P > 0.05). Genotypic and allelic frequencies were found to be significantly different in four breeds (χ² = 147.674, df = 6, P < 0.001; χ² = 157.250, df = 3, P < 0.001, respectively), implying that the distribution of genotypic and allelic frequencies of goat HESX1 gene was significantly associated with different goat utilities (cashmere, meat and dairy). Association analysis results revealed no significant effects of caprine HESX1 gene on body sizes in XNSN population (P > 0.05) and cashmere traits in IMWC population (P > 0.05). Significant statistical of HESX1 gene with body weight was found (*P < 0.05). The genotype AA showed significantly higher body weight than those of AG in 2-year-old age (*P < 0.05), while the AA genotype was senior to AG genotype in 4-year-old body weight trait (*P < 0.05). These suggestions indicated that the HESX1 gene has significant effect on goat body weight depending on ages, which is accordance with the function repressor of the HESX1.     

Polymorphism of <Emphasis Type="Italic">DLK1</Emphasis> and <Emphasis Type="Italic">CLPG</Emphasis> gene and their association with phenotypic traits in Chinese cattle

DLK1 and CLPG were located in DLK1-GTL2 imprinted cluster. They all affected muscle growth and meat tenderness. The functional importance of DLK1 and CLPG imply that the variation of the genes could affect the growth traits of animal. PCR-SSCP and sequencing were used to analyze the four loci of DLK1 gene and CLPG gene in 1109 individuals, which belong to eight breeds/species of bovidae, including cattle, buffalo and yak. A synonymous mutation (C451T) was detected in exon 5 of DLK1 in Qinchuan cattle, but didn’t change significantly with phenotypic traits. Three genotypes AA, AB and AC of CLPG were identified in Jiaxian cattle. The associations analyst of different genotypes showed that the individuals with genotypes AA and AC had a greater body weight and longer body length than those with genotype AB (P < 0.05 and P < 0.01, respectively); the AA individuals were different from those AB (P < 0.05) in the circumference of cannon bone. No polymorphism was observed in the other populations at other loci. These results were in agreement with the homology analysis: DLK1 and CLPG genes were in a highly conserved.     

A post‐GWAS confirming the SCD gene associated with milk medium‐ and long‐chain unsaturated fatty acids in Chinese Holstein population

The stearoyl‐CoA desaturase (delta‐9‐desaturase) gene encodes a key enzyme in the cellular biosynthesis of monounsaturated fatty acids. In our initial genome‐wide association study (GWAS) of Chinese Holstein cows, 19 SNPs fell in a 1.8‐Mb region (20.3–22.1 Mb) on chromosome 26 underlying the SCD gene and were highly significantly associated with C14:1 or C14 index. The aims of this study were to verify whether the SCD gene has significant genetic effects on milk fatty acid composition in dairy cattle. By resequencing the entire coding region of the bovine SCD gene, a total of six variations were identified, including three coding variations (g.10153G>A, g.10213T>C and g.10329C>T) and three intronic variations (g.6926A>G, g.8646G>A and g.16158G>C). The SNP in exon 3, g.10329C>T, was predicted to result in an amino acid replacement from alanine (GCG) to valine (GTG) in the SCD protein. An association study for 16 milk fatty acids using 346 Chinese Holstein cows with accurate phenotypes and genotypes was performed using the mixed animal model with the proc mixed procedure in sas 9.2. All six detected SNPs were revealed to be associated with six medium‐ and long‐chain unsaturated fatty acids (P = 0.0457 to P < 0.0001), specifically for C14:1 and C14 index (P = 0.0005 to P < 0.0001). Subsequently, strong linkage disequilibrium (D′ = 0.88–1.00) was observed among all six SNPs in SCD and the five SNPs (rs41623887, rs109923480, rs42090224, rs42092174 and rs42091426) within the 1.8‐Mb region identified in our previous GWAS, indicating that the significant association of the SCD gene with milk fatty acid content traits reduced the observed significant 1.8‐Mb chromosome region in GWAS. Haplotype‐based analysis revealed significant associations of the haplotypes encompassing the six SCD SNPs and one SNP (rs109923480) in a GWAS with C14:1, C14 index, C16:1 and C16 index (P = 0.0011 to P < 0.0001). In summary, our findings provide replicate evidence for our previous GWAS and demonstrate that variants in the SCD gene are significantly associated with milk fatty acid composition in dairy cattle, which provides clear evidence for an increased understanding of milk fatty acid synthesis and enhances opportunities to improve milk‐fat composition in dairy cattle.     

A novel SNP of the GHRL gene in goat and its association with growth traits

Ghrelin, an endogenous gland for growth hormone secretagogue receptor, has been shown to stimulate food intake and control energy homeostasis and lipid metabolism. So, ghrelin precursor (GHRL) gene is a potential candidate gene for caprine growth traits. In this study, we detected the polymorphism of the caprine GHRL gene by PCR–SSCP and DNA sequencing methods in 459 individuals from four goat breeds. A novel single nucleotide polymorphism (SNP) (IVS2 + 147G > A) was detected. Frequencies of IVS2 + 147G allele varied from 0.842 to 1.000. The association of IVS2 + 147G > A with growth traits was analyzed and IVS2 + 147G > A was shown to be associated with growth traits. Individuals with genotype AG were significantly higher than those of individuals with genotype GG in circumference of chest and cannon and trunk index (P < 0.05 or P < 0.01). Moreover, there was a tendency that genotype AG individuals had better performance in other aspects such as body height and body length than genotype GG individuals although no significant differences appeared (P > 0.05). We suggested that IVS2 + 147G > A could be a perfect molecular marker in marker-assisted selection (MAS).     

A coding SNP of <Emphasis Type="Italic">LHX4</Emphasis> gene is associated with body weight and body length in bovine

Heterozygous mutations in LHX4 are associated with combined pituitary hormone deficiency. In this study, the polymorphism of LHX4-HaeIII locus was revealed in 822 individuals from four Chinese cattle breeds. The PCR–RFLP analysis showed that there were three genotypes: GG, GA, AA. The frequencies of genotype GG ranged from 0.6620 to 0.9789 in analyzed populations. The genotypic frequencies of LHX4 locus in the four populations all agreed with Hardy–Weinberg equilibrium (P > 0.05). Distributions of genotypic frequencies of different breeds (QC, NY, JX, CH) at this locus were found to be significantly different based on a χ ² test (P < 0.001). The genetic diversity analysis revealed the JX cattle possessed intermediate genetic diversity, and the other three Chinese cattle breeds belonged to poor genetic diversity. Correlation analysis with growth traits in the NY breed indicated that: the animals with genotype GA had greater body weight than those with genotype GG (P < 0.05); the animals with GA genotype owned significantly longer body length than the ones with GG genotype (P < 0.05) at 18 and 24 months.