Combined effects of four SNPs within goat PRLR gene on milk production traits

Xinong Saanen (SN, n = 323) and Guanzhong (GZ, n = 197) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions with their intron–exon boundaries of prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction (PIRA-PCR) and PCR-restriction fragment length polymorphism (PCR-RFLP). Four novel SNPs (g.40452T > C, g.40471G > A, g.61677G > A and g.61865G > A) were identified. The g.61677G > A and g.61865G > A SNPs caused amino acid variations p.Ser485Asn and p.Val548Met, respectively. Both g.40452T>C and g.40471G>A loci were closely linked in SN and GZ goat breeds (r² > 0.33). In addition, there was also a close linkage between g.61677G>A and g.61865G>A loci in both goat breeds. Statistical results indicated that the g.40452T > C, g.61677G > A and g.61865G > A SNPs were significantly associated with milk production traits in SN and GZ breeds. Further analysis revealed that combinative genotype C1 (TTAAGGGG) was better than the others for milk yield in SN and GZ goat breeds. These results are consistent with the regulatory function of PRLR in mammary gland development, milk secretion, and expression of milk protein genes, and extend the spectrum of genetic variation of the caprine PRLR gene, which might contribute to goat genetic resources and breeding.