Identification of SNPs within the sheep <Emphasis Type="Italic">PROP1</Emphasis> gene and their effects on wool traits

Regarding mutations of PROP1 (Prophet of POU1F1) gene significantly associating with combined pituitary hormone deficiency (CPHD) in human patients and animals, PROP1 gene is a novel important candidate gene for detecting genetic variation and growth, reproduction, metabolism traits selection and breeding. The aim of this study was to detect PROP1 gene mutation of the exon 1–3 and its association with wool traits in 345 Chinese Merino sheep. In this study, on the basis of PCR-SSCP and DNA sequencing methods, ten novel SNPs within the sheep PROP1 gene, namely, AY533708: g.45A > G resulting in Glu15Glu, g.1198A > G, g.1341G > C resulting in Arg63Ser, g.1389G > A resulting in Ala79Ala, g.1402C > T resulting in Leu84Leu, g.1424A > G resulting in Asn91Ser, g.1522C > T, g.1556A > T, g.1574T > C, g.2430C > G were reported. In addition, association analysis showed that three genotypes of P4 fragment were significantly associated with fiber diameter in the analyzed population (P = 0.044). These results strongly suggested that polymorphisms of the PROP1 gene could be a useful molecular marker for sheep breeding and genetics through marker-assisted selection (MAS).     

Genetic variants in MYF5 affected growth traits and beef quality traits in Chinese Qinchuan cattle

Myogenic factor 5 plays actively roles in the regulation of myogenesis. The aims of this study are to identify the evolution information of MYF5 protein among 10 domestic and mammalian animals, to uncover the expression patterns of MYF5 gene in calves and adults of Qinchuan cattle, and to expose the genetic variants of the MYF5 gene and explore its effect on cattle growth traits and beef quality traits in Qinchuan cattle. The bioinformatics results showed that the MYF5 proteins highly conserved in different mammalian or domestic animals apart from chicken. The expression level of MYF5 gene in the heart, muscle, lung, large intestine and liver was greater than that of other tissues. PCR amplicons sequencing identified four novel SNPs at g.5738A>G, g.5785C>T and g.5816A>G in the 3rd exon region and g.6535A>G in the 3’ UTR. Genotypic frequencies of g.5785C>T was harshly deviated from the HWE (P < .05). Genetic diversity was low or intermediate for the four SNPs and those SNPs were in the weak linkage disequilibrium. Association analysis results indicated g.5785C>T, g.5816A>G and g.6535A>G significant effect on growth performance and beef quality traits of Qinchuan cattle. H1H3 diplotype had greater body size and better beef quality. All the results implicate that the MYF5 gene might be applied as a promising candidate gene in Qinchuan cattle breeding.     

Association of hormone-sensitive lipase (HSL) gene polymorphisms with the intramuscular fat content in two Chinese beef cattle breeds

Hormone-sensitive lipase (HSL) was considered as an essential enzyme in glucolipid metabolism. It has been proposed to be a lead candidate gene for genetic markers of lipid deposition in livestock. The aim of this study was to identify sequence variants (SVs) of the bovine HSL gene and evaluate the relations to intramuscular fat in two indigenous Chinese beef cattle breeds. Expression analysis by quantitative real-time polymerase chain reactions (qPCR) indicated that expression levels of bovine HSL gene were highest in the perirenal fat and heart within two different age stage (adult and calf), respectively. Five SVs were identified by direct DNA sequencing, which included four missense mutations (g.16563C>T, g.16734G>A, g.16896A>G, g.17388G>T) in exon 8 and a synonymous mutation (g.17402C>T) in exon 9. Population genetic analysis showed that except for g.16563C>T and g.17402C>T, all the other detected SVs strongly affected the bovine intramuscular fat content (P < 0.01 or P < 0.05). The individuals with Hap5/5 diplotypes (CC-GG-GG-GG-CC) was highly significantly associated with intramuscular fat content than the other diplotypes (P < 0.01). The above results suggested that the HSL gene can used as potential candidate markers gene for the beef breed improvement through marker assisted selection in Chinese cattle breeds.     

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.     

Identification of novel single nucleotide polymorphisms (SNPs) of the lipoprotein lipase (LPL) gene associated with fatty acid composition in Korean cattle

The lipoprotein lipase (LPL) gene can be considered a functional candidate gene that regulates fatty acid composition. In this study, genetic associations between fatty acid composition and exonic single nucleotide polymorphisms (SNPs) in the LPL gene were examined using 612 Korean cattle. We investigated the relationship between unsaturated fatty acids and five novel SNPs (c.322G>A, c.329A>T, c.527T>G, c.988C>T and c.1591G>A), and confirmed that three polymorphic SNPs (c.322G>A, c.329A>T and c.1591G>A) were associated with fatty acid composition. Korean cattle with an AA genotype of c.322G>A, c.329A>T, and GA genotype of c.1591G>A had higher levels of monounsaturated fatty acids and carcass traits (P < 0.05). Our findings confirmed that three novel SNPs we identified in the LPL gene can affect fatty acid composition and carcass traits. Therefore, selection for AA and GA genotypes should be recommended to genetically improve beef quality and flavor.     

Detection of SNPs in the TBC1D1 gene and their association with carcass traits in chicken

TBC1D1 plays an important role in numerous fundamental physiological processes including muscle metabolism, regulation of whole body energy homeostasis and lipid metabolism. The objective of the present study was to identify single nucleotide polymorphisms (SNPs) in chicken TBC1D1 using 128 Erlang mountainous chickens and to determine if these SNPs are associated with carcass traits. The approach consisted of sequencing TBC1D1 using a panel of DNA from different individuals, revealing twenty-two SNPs. Among these SNPs, two polymorphisms (g.69307744C>T and g.69307608T>G) of block 1, four polymorphisms (g.69322320C>T, g.69322314G>A, g.69317290A>G and g.69317276T>C) of block 2 and four polymorphisms of block 3 (g.69349746G>A, g.69349736C>G, g.69349727C>T and g.69349694C>T) exhibited a high degree of linkage disequilibrium in all test populations. An association analysis was performed between the twenty-two SNPs and seven performance traits. SNPs g.69307744C>T, g.69340192G>A and g.69355665T>C were demonstrated to have a strong effect on liveweight (BW), carcass weight (CW), semi-eviscerated weight (SEW) and eviscerated weight (EW) and g.69340070C>T polymorphism was related to BW, SEW and BMW in chicken populations. However, for the other SNPs, there were no significant correlations between different genotypes and carcass traits. Meanwhile, haplotype CT–TG of block 1 and combined genotype AG–TT–AC–CT of block 3 were significantly associated with BW, CW, SEW and EW. Overall, our results provide evidence that polymorphisms in TBC1D1 are associated with carcass traits and would be a useful candidate gene in selection programs for improving carcass traits.     

Nucleotide variants in prion-related protein (testis-specific) gene (PRNT) and effects on Chinese and Mongolian sheep phenotypes

ABSTRACT

Studies of the ovine prion-related protein (testis-specific) gene (PRNT), including studies of genetic diversity, have highlighted its potential relationship to scrapie infection and economically important ovine traits. PRNT was previously reported to be highly polymorphic in Portuguese sheep. To characterize genetic polymorphisms in this gene in Asian sheep, a direct sequencing method was used to detect polymorphic loci in PRNT in 285 individual sheep from four Chinese and one Mongolian breeds. Seven SNP variants in PRNT were identified, including three novel variants (g.93G>A, g.162G>T, and g.190A>G) and four previously reported variants (g.17 C>T, g.112G>C, g.129C>T, and g.144A>G). In the five breeds that we analyzed, the mutation frequencies of g.190A>G in Lanzhou Fat-tail sheep (LFTS) and g.129C>T in the other four varieties were high (F>0.5). Moreover, thirteen different haplotypes that had a comparable distribution in the tested breeds were also identified; ‘C-G-G-C-A-G-A’ occurred at the highest frequency in the five sheep breeds. Additionally, we previously explored the significance of relationships between polymorphisms in PRNP or PRND and ovine growth performance. Here, we also performed correlation analysis in all tested loci. These loci polymorphisms were significantly associated with ten different growth traits (P<0.05), except for g.93G>A. Meanwhile, in contrast to a previous study, there was no significant association between the seven SNP loci analyzed and our previously reported sheep PRND or PRNP insertion/deletion mutations. Our findings may provide new insights into polymorphic variation in ovine PRNT, which may contribute to genetic improvements in economic traits that are important for sheep breeding.     

Polymorphisms in the bovine ghrelin precursor (<Emphasis Type="Italic">GHRL</Emphasis>) and Syndecan-1 (<Emphasis Type="Italic">SDC1</Emphasis>) genes that are associated with growth traits in cattle

Transgenically expressed Syndecan-1 was found in the hypothalamic nuclei that control energy balance, and was associated with maturity-onset obesity, while ghrelin has been shown to play important roles in the control of food intake, gastric acid secretion, energy homeostasis, and glucose and lipid metabolism. However, the roles of genetic variations of Syndecan-1 and ghrelin on growth trait have few been reported in cattle. Herein, five Chinese cattle breeds were analyzed by PCR–SSCP and DNA sequencing methods. The bovine ghrelin gene showed eleven SNPs g.[267G>A, 271G>A, 290C>T, 326A>G, 327T>C, 420C>A, 569A>G, 945C>T, 993C>T, 4491A>G, 4644G>A] and three SNPs g.[420C>A, 569 A>G, 945C>T] were firstly detected in cattle. The bovine Syndecan-1 gene showed two SNPs. One SNP showed a transition C>G at position 21514, resulting in a synonymous mutation p.G(GGC)169G(GGG) and another showed a transversion C>T at position 22591, resulting in a synonymous mutation p.D(GAC)283D(GAT). In ghrelin gene, no significant associations were revealed between any variant sites and body weight, average daily gain, body sizes for different growth periods (6, 12, 18, and 24 months old), as well as for the milk yield at 305 days, milk protein rate and milk fat percentage. However, the polymorphism of Syndecan-1 gene was significantly associated with bovine birth weight and body length. Hence, we first suggested that Syndecan-1 gene could be regarded as molecular marker for superior birth weight and body length.     

Genetic association of <Emphasis Type="Italic">Phosphodiesterase 1B</Emphasis> (PDE1B) with carcass traits in Korean cattle

Quantitative trait loci for fat deposition and carcass traits have been identified in the vicinity of the gene encoding phosphodiesterase 1B (PDE1B) on bovine chromosome 5. Therefore, the PDE1B gene can be considered as a positional and functional candidate gene for carcass traits in beef cattle. This study aimed to identify single nucleotide polymorphisms (SNPs) in the PDE1B gene and to evaluate their associations with carcass traits in Korean cattle. Eight SNPs, g.440T>G, g.17122A>G, g.17507A>C, g.17575A>G, g.17607T>C, g.17609C>A, g.17692C>T, and g.17707C>G, were identified in the region ranged from exon 1 to intron 6. Five of them were used for association analysis because of their availability of restriction fragment length polymorphisms. As a result, g.17122A>G in intron 3 was significantly associated with backfat thickness (BFT), and g.17507A>C in exon 5 was associated with longissimus dorsi muscle area (LMA, P < 0.05). Animals with the AG genotype of g.17122 had thicker BFT than those with the AA genotype. Animals with the AA or AC genotype of g.17507A>C had larger LMA than those with the CC genotype. We suggested the PDE1B gene as a candidate gene for carcass traits of beef cattle. Fine mapping would be required for application to marker-assisted selection.     

Analysis of non‐synonymous SNPs of the porcine SERPINA6 gene as potential causal variants for a QTL affecting plasma cortisol levels on SSC7

Recently, the SERPINA6 gene encoding corticosteroid‐binding globulin (CBG) has been proposed as a candidate gene for a quantitative trait locus (QTL) affecting cortisol level on pig chromosome 7. The QTL was repeatedly detected in different lines, including a Piétrain × (German Landrace × German Large White) cross (PiF1) and purebred German Landrace (LR). In this study, we investigated whether the known non‐synonymous polymorphisms c.44G>T, c.622C>T, c.770C>T, c.793G>A, c.832G>A and c.919G>A of SERPINA6 are sufficient to explain the QTL in these two populations. Our investigations revealed that SNPs c.44G>T, c.622C>T, c.793G>A and c.919G>A are associated with cortisol level in PiF1 (P < 0.01). Haplotype analysis showed that these associations are largely attributable to differences between a major haplotype carrying SNPs c.793G>A and c.919G>A and a haplotype carrying SNPs c.44G>T and c.622C>T. Furthermore, some SNPs, particularly c.44G>T and c.622C>T and the carrier haplotype, showed association with meat quality traits including pH and conductivity (P < 0.05). In LR, the non‐synonymous SNPs segregate at very low frequency (<5%) and/or show only weak association with cortisol level (SNPs c.832G>A and c.919G>A; P < 0.05). These findings suggest that the non‐synonymous SNPs are not sufficient to explain the QTL across different breeds. Therefore, we examined whether the expression of SERPINA6 is affected by cis‐regulatory polymorphisms in liver, the major organ for CBG production. We found allelic expression imbalance of SERPINA6, which suggests that its expression is indeed affected by genetic variation in cis‐acting elements. This represents candidate causal variation for future studies of the molecular background of the QTL.     

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.     

Associations of FASN gene polymorphisms with economical traits in Nellore cattle (Bos primigenius indicus)

The aim of this study was to identify molecular markers to be applied to marker-assisted selection. Three SNPs of the FASN gene were studied. PCR–RFLP was used for genotyping. The SNPs g.17924A>G, g.17860C>T and g.15603A>G all in the FASN gene were genotyped using the enzymes MscI, DdeI and Hae III, respectively. The animals were raised in extensive systems and belong to three lines selected for growth as part of the Selection Program of Zebu and Caracu Breeds, S?o Paulo, Brazil. Allele and genotype frequencies were compared between selection lines using the Genepop 3.4. Associations between polymorphisms and the traits studied were evaluated using the PROC MIXED procedure of the SAS/STAT 9.1.3. The G and C alleles were the most frequent alleles of the g.15603A>G and g.17860C>T loci, respectively. The g.17924A>G locus showed no polymorphism in the population studied. Allele and genotype frequencies differed significantly between the NeT line and the NeC and NeS lines. The g.15603A>G polymorphism tended to exert an additive effect on rump fat thickness and male yearling height. For g.17860C>T, an additive effect on male yearling height was observed. Genotype combination analysis revealed a significant effect on loin eye area. Although this study provided evidence of an association between the FASN gene and some traits, more detailed analyses are needed to obtain more efficient molecular markers.     

Novel Mutation in the <Emphasis Type="Italic">MECP2</Emphasis> Gene Identified in a Group of Rett Syndrome Patients from Ukraine

Mutations in the MECP2 gene are known to cause Rett syndrome (RTT)—a neurodevelopmental disorder, one of the most common causes of intellectual disability in females, with an incidence of 1 in 10000–15000. We have investigated exons 3 and 4 of the MECP2 gene, that coding MBD and TRD domains of the MeCP2 protein, in 21 RTT patients from Ukraine by PCR-DGGE analysis followed by Sanger sequencing of PCR fragments with abnormal migration profiles. In 13 of 21 (61.9%) patients 7 different mutations were identified one nonsense mutation—c. NC_000023.11:g.154031326G>A (MECP2:c.502C>T) and 4 missense mutation NC_000023.11:g.154031409G>T (MECP2:c.419C>T), NC_000023.11:g.154031355G>A (MECP2:c.473C>T), NC_000023.11:g.154031354A>C (MECP2:c.472A>C), NC_000023.11:g.154031431G>A (MECP2:c.397C>T) located in exon 4, a rare RTT-causing splice site mutation NC_000023.10:g.153296903T>G (MECP2:c.378-2A>C) in intron 3 and deletion NC_000023.10:g.1532 96079_153296122del44 in exon 4. The novel mutation MECP2:c.472A>C identified in our study in patients withclassic RTT phenotype leds to T158P substitution. It is one more confirmation of crucial role that 158 codon in MECP2 protein function.     

Genetic variants of the <Emphasis Type="Italic">FABP4</Emphasis> gene are associated with marbling scores and meat quality grades in Hanwoo (Korean cattle)

This study was aimed to search new genetic variants in the bovine FABP4 gene as molecular markers for meat quality and carcass traits. PCR–RFLP analysis revealed that three SNPs located at nucleotide positions g.2834C>G, g.3533T>A, and g.3691G>A were identified based on a GenBank accession number (NC_007312.4). Sequence analysis revealed that SNPs were located in intron 1 (g.2834C>G) and 2 (g.3533T>A), and an exon 3 (g.3691G>A), showing allele frequencies as 0.592, 0.579, and 0.789, respectively. Genetic variabilities of heterozygosity (He) and polymorphic information contents (PIC) were estimated for g.2834C>G (0.608 and 0.531), g.3533T>A (0.615 and 0.539), and g.3691G>A (0.498 and 0.401) loci, respectively. A SNP located in the exon 3 of FABP4 was characterized and associated with desirable increases of MS (marbling scores) and MG (meat quality grades) in Hanwoo. The statistical analysis revealed that additive effects by GG genotypes in g.3691G>A SNP were significantly greater than AA genotypes in MS and MG traits. These findings suggest that the FABP4g.3691G>A SNP will be a useful candidate locus to maximize economic benefits for cattle populations.     

Polymorphism distribution of RYR1, PRKAG3, HFABP,MYF-5 and MC4R genes in crossbred pigs

This study was designed to screen the crossbred pigs for SNPs in five candidate genes, associated with pork quality traits and to differentiate their genotypes by PCR–RFLP. The results indicated that genotypes of crossbred pigs were NN (90%) and Nn (10%) for RYR1; RR (83%) and QR (17%) for PRKAG3; HH (98%), Hh (1%) and hh (1%) for HFABP; DD (99%) and CD (1%) for MYF-5; and AG (57%), GG (26%) and AA (17%) for MC4R SNPs, respectively. Allelic frequencies for five SNPs {RYR1 (1843C>T), PRKAG3 (c.599G>A), HFABP (c.1322C>T), MYF-5 (c.1205A>C) and MC4R (c.1426A>G)} were 0.95 and 0.05 (N/n), 0.08 and 0.92 (Q/R), 0.99 and 0.01 (H/h), 0.00 and 1.00 (C/D) and 0.45 and 0.55 (A/G), respectively. The effect of RYR1 (1843C>T) SNP was significant on pH₄₅ (P?<?0.05), pH₂₄ (P?<?0.05) and protein % (P?<?0.05). The PRKAG3 (c.599G>A) and MC4R (c.1426A>G) SNP had significant association with dressing percentages. The results revealed that RYR1, PRKAG3 and MC4R SNPs may be used in marker associated selection for pork quality traits in crossbred pigs.     

Associations of polymorphisms in the promoter I of bovine acetyl‐CoA carboxylase‐α gene with beef fatty acid composition

The objectives of this study were to identify single nucleotide polymorphisms (SNPs) in the promoter I (PI) region of the bovine acetyl‐CoA carboxylase‐α (ACACA) gene and to evaluate the extent to which they were associated with lipid‐related traits. Eight novel SNPs were identified, which were AJ276223:g.2064T>A (SNP1), g.2155C>T (SNP2), g.2203G>T (SNP3), g.2268T>C (SNP4), g.2274G>A (SNP5), g.2340A>G (SNP6), g.2350T>C (SNP7) and g.2370A>G (SNP8). Complete linkage disequilibrium was observed among SNP1, 2, 4, 5, 6 and 8. Phenotypic data were collected from 573 cross‐bred steers with six sire breeds, including Hereford, Angus, Brangus, Beefmaster, Bonsmara and Romosinuano. The genotypes of SNP1/2/4/5/6/8 were significantly associated with adjusted backfat thickness. The genotypes of SNP3 were significantly associated with triacylglycerol (TAG) content and fatty acid composition of longissimus dorsi muscle (LM) in Brangus‐, Romosinuano‐ and Bonsmara‐sired cattle. Cattle with g.2203GG genotype had greater concentrations of TAG, total lipid, total saturated fatty acid and total monounsaturated fatty acid than did cattle with g.2203GT genotype. The genotypes of SNP7 were significantly associated with fatty acid composition of LM. Cattle with genotype g.2350TC had greater amounts of several fatty acids in LM than did cattle with genotype g.2350CC. Our results suggested that the SNPs in the PI region of ACACA gene are associated with variations in the fatty acid contents in LM.     

Identification of novel SNPs in the Sarda breed goats POU1F1 gene and their association with milk productive performance

The aim of the study was to detect polymorphism in the POU1F1 gene in Sarda breed goat, as well as to establish if SNPs could be associated with milk productive traits. The research was conducted on 129 Sarda breed goats from 4 to 5 years old, multiparous, lactating and in their third to fifth lactation. We report nine exonic and seven non-coding regions SNPs within the Sarda goat POU1F1 gene, namely, Ex 1 61 G>C; Ex 1 108 G>A; Ex 3 C>T; Ex 3 92 C>T; Ex 4 110 A>G; Ex 5 34 G>A resulting in Arg213Lys change; IVS4 641 G>A, IVS4 643 A>C, IVS4 659 G>A, IVS4 677 A>C, IVS4 G699Del, IVS4 709 C>G, Ex 6 17 G>A resulting in Arg228Ser change, Ex 6 58 G>T, Ex 6 172 T>C, 3′UTR 110 T>C. A statistically significant association was found between genotype TT, in position 17 of the exon 6 (3.1 % of frequency), and increased milk yield (P < 0.01) while genotype GT (25.6 % of frequency) was associated with a higher fat content. Genotype TT in position 58 of the exon 6 (3.9 % of frequency) was found to be associated with a higher fat (P < 0.01) and protein content (P < 0.05). Twenty-eight haplotypes were detected, but no significant association between the haplotypes and the milk production traits have been found. Our data, as well as providing new SNPs extending the POU1F1 gene characterization, evidence a relationship between polymorphism and milk production traits in Sarda goat breed.     

Effect of FASN gene on milk yield and milk composition in the Chinese Holstein dairy population

Fatty acid synthase (FASN) is a multifunctional protein that catalyzes de novo synthesis of fatty acids in cells. It plays a key role in the lipid biosynthesis as well as in the general metabolism of all living animals. We herein investigated polymorphisms of FASN. As a result, six single nucleotide polymorphisms (SNPs) were found and then genotyped in 752 Chinese Holstein cows. It was found that g.17924A>G was non‐synonymous, g.13965 C>T, g.16907 T>C and g.18663T>C were synonymous mutations and two other two SNPs, g.8948 C>T (ss491228481) and g.14439T>C (rs133498277), were in intronic sequences of the gene. All such identified SNPs were found to be associated with milk yield and composition traits (P = 0.0441 to <0.0001). Significant additive and allele substitution effects were observed for three yield traits at all six loci as well (P < 0.05 to <0.01). Complete linkage disequilibrium among the five SNPs, with the exception of g.8948 C>T, was observed.