Association of Genetic Variants ofB-Lactoglobulin Gene with Milk Production in a Herd and a Superior Family of Holstein Cattle

Document Type: Brief Report


1 Department of Animal Sciences, Faculty of Agriculture, Gorgan University of Agricultural Sciences and Natural Resources, P.O. Box 49189-43464 , Gorgan, I.R. Iran

2 Department of Animal Sciences, Faculty of Agriculture, P.O. Box 163, Gonbad, I.R. Ira


Polymorphism of the b-lactoglobulin (b-LG) gene in 101 cows belonging to the Holstein herd and a superior cow, producing more than 150 Kg milk/day, together with four offsprings was investigated by the Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) method. In the Holstein herd, the alleles A and B of the b-LG gene had frequencies of 0.53 and 0.47, respectively. The genotypes AA, AB and BB of the b-LG gene were estimated to have frequencies of 0.257, 0.544 and 0.198, respectively. Genotypes were distributed according to the Hardy-Weinberg equilibrium. Results indicated that the b-LG genotypes significantly affected (P< 0.01) milk yield (genotype AA being more effective than genotype BB). The superior cow and her progenies were all heterozygotes (AB).


b-lactoglobulin (b-LG) is the major whey protein of ruminant species and is also  present in the milk of many, but not all mammalian species (Kontopidis et al., 2004). b-LG is located on chromosome 11 in the cow (Daniela et al., 2008) and has 12 known variants in which A and B variants are the most frequent ones (Rachagani et al., 2006). The b-LG binding affinity for retinol (Frapin et al., 1993) and fatty acids (Fugate and Song, 1980) has been demonstrated, which suggests a possible role of b-LG in the transport and metabolism of these components (Lum et al., 1997; Eigel et al., 1984). There are many contradictory reports on the association between genetic variants of b-LG and milk production and composition (Tsiaras et al., 2005; Aleandri et al., 1990; Ng-Kawi-Hang et al., 1986; Ng-Kawi-Hang et al., 1984). Several researches have shown that the BB genotype of b-LG is associated with higher fat and  increased cheese yields (Wedholm et al., 2006; Aleandri et el., 1990; Ng-Kawi-Hang et al., 1986). An association of β-LG BB with higher protein yield has also been reported (Ng-Kawi-Hang et al., 1984). On the contrary, Bovenhuis et al. (1992) reported lower protein yield for this genotype. In various studies, the AA genotype of the b-LG gene has been associated with higher milk production (Ikonen et al., 1999; Bovenhuis et al., 1992; Ng-Kwai-Hang et al., 1990) and higher milk protein (Strazalkowska et al., 2002; Ng-Kwai-Hang et al., 1990; Ng-Kwai-Hang et al., 1984). In other reports, the AB genotype has been associated with higher milk production (Tsiaras et al., 2005; Kaygisiz and Douan, 1999), lactose and protein yield (Tsiaras et al., 2005).
The current study was carried out from July 2008 to November 2008 in the Molecular Genetics laboratory at the Department of Animal Sciences at Gorgan University of Agricultural Sciences and Natural Resources. The aim was to investigate b-LG polymorphism and its association with milk yield in a herd of Holstein cows. The blood samples were collected randomly from 101 Holstein cows (Behin Talise farm in Gorgan, Iran) having at least one lactation record. In the meantime, blood samples of a superior cow (producing more than 150 Kg milk/day) and her four female progenies from a farm near the region of the investigated herd were also collected. DNA was extracted, using the salting out extraction protocol (Miller et al., 1988). The 247 bp fragment, comprising a part of the IV exon and intron of the genomic DNA was amplified by using primers as suggested by Strazalkowska et al. (2002). The sequences of the primers were as follows:
forward: 5´-TGT GCT GGA CAC CGA CTA CAA AAA G-3´; and reverse: 3´-GCT CCC GGT ATA TGA CCA CCC TCT-5´.
Polymerase Chain Reaction (PCR) was carried out, using the Personal Cycler™ thermocycler (Biometra, Germany) and the PCR Master Kit (CinnaGen Inc., Iran). The kit master mix consisted of 0.04 U/ml of Taq DNA polymerase, 10X PCR buffer, 3mM MgCl2 and 0.04 mM dNTPs (each). Each reaction mixture consisted of 12.5 ml of the master mix, 1 ml of the DNA solution (50 to 100 ng/ml), 1 ml of each primer (5 pmol/µl) and some deionized water making up a final volume of 25 ml. The amplification program was as follows: initial denaturation at 94ºC for 5 min followed by 35 cycles of denaturation at 94°C for 1 min, annealing at 60ºC for 1 min, extension at 72°C for 1 min and a final extension at 72°C for 5 min. Digestion of PCR products was carried out, using the HaeIII restriction endonuclease (Fermentas, European Union) at 37°C for 16 h and then analyzed by electrophoresis on a 2.5% (w/v) agarose gel (Fig.1). Determination of gene and genotype frequencies and the c² test were carried out using the POPGene 32 software (Yeh et al., 1997).
The b-LG allele and genotype frequencies of the herd and those of other studies are summarized in Table 1. In most of the studies, the frequency of the b-LG B allele was more than the A allele. The frequencies of the genetic variants A and B varied from 0.23 to 0.52 and 0.47 to 0.77, respectively in the Holstein breed. The results of the current study are in agreement with other studies, particularly with those of Ron et al. (1994) and Tsiaras et al. (2005). The maximum differences between A and B alleles’ frequencies have been observed in the Ayrshire (Lin and Mcallister, 1986) and Sahiwal (Rachagani et al., 2006) breeds. It must be noted that differences between allele frequencies of the gene depend on economic strategies and animal breeding programs for every breed and herd. Genotype of the superior cow and her progenies were determined as AB by using the same method.
A single-trait fixed model was used to examine the effect of b-LG genotypes on milk production. The mature equivalent milk yield records, corrected for missing data, were used. Statistical analysis was performed using the Generalized Least squares Means (GLM) procedure of the SAS software (2002) and comparison of the least squares means was carried out using the Tukey-Kramer test (SAS software, 2002).
The following fixed model was used:
y = Xb + e,
with E(y) = Xb and var(y) = V = var(e)
X = known incidence matrix relating observation in y to the fixed effects in b,
b = vector of levels of fixed effects,          e = vector of random residual terms,
E = the expectation of y and        V = the variance of y
The fixed effects are represented by the calving years (2006, 2007 and 2008), seasons with 4 classes, parity with 8 classes and b-LG genotypes (AA, AB, BB).
Comparisons of the lactation milk yields of the b-LG genotypes showed that cows with the AA genotype produced more milk than animals with the BB genotype (P< 0.006) (Table 2). Comparison of expected and observed genotypic frequencies yielded a c² value of 0.78, suggesting that the b-LG locus was in Hardy-Weinberg equilibrium (P< 0.05).
The significantly higher milk production yield of cows with the AA genotype was in agreement with the results of other researches (Ikonen et al., 1999; Bovenhuis et al., 1992; Aleandri et el., 1990; Ng-Kwai-Hang et al., 1990).
Although there was a difference between AA and AB genotypes, but it was statistically not significant (Table 2). The superior cow, because of having both the A and B alleles and the AB genotypes, apparently was superior in milk production and fat yields, as well as the quality of cheese. Some reports have confirmed significant association between the AB genotype and higher milk yield (Tsiaras et al., 2005; Kaygisiz and Douan, 1999). Therefore, these genotypes (AA, AB) seem to be suitable candidates for selection aiming at improving milk production.
It must be pointed out that other factors might also affect milk yield of the cow, such as physiological status and other major genes. Hence, sequencing the genome of this cow will be an attractive theme for future studies to understand the main factors affecting this high production level.


The authors would like to thank Mr. Kabiri, owner of the Behin Talise farm, and Mr. Khari, owner of superior cow, for their kind help in providing blood samples. This research was financially supported by the University of Agricultural Sciences and Natural Resources, Gorgan, I.R. Iran.

Aleandri R, Buttazzoni LG, Schnerder JC (1990). The Effects of milk protein polymorphisms on milk components and cheese-producing ability. J Dairy Sci. 73: 241-255.
Bovenhuis H, Van Arendonk JAM, Korver S (1992). Associations between milk protein polymorphisms and milk production traits. J Dairy Sci. 75: 2549-2559.
Celik S (2003). b-lactoglobulin genetic variants in Brown Swiss breed and its association with compositional properties and rennet clotting time of milk. Int Dairy J. 13: 727-731.
Daniela I, Aurelia S, Anuta M, Claudia S, Vintila I (2008). Genetic polymorphism at the b-lactoglobulin locus in a dairy herd of Romanian Spotted and Brown of Maramures breeds. Zootehnie si Biotehnologii. 41: 104-107. 
Eenennaam AV, Medrano JF (1991). Milk protein polymorphisms in California dairy cattle. J Dairy Sci. 74: 1730-1742.
Eigel WN, Butler JE, Ernstrom CA, Farrell HM, Harwalker VR, Jenness R, Whitney RM (1984). Nomenclature of proteins of cow’s milk: fifth edition. J Dairy Sci. 76: 1599.
Frapin D, Dufour E, Haertle T (1993). Probing the fatty acid binding site of b-lactoglobulins. J Protein Chem. 12: 443-449.
Fugate RD, Song PS (1980). Spectrocopic characterization of b-lactoglobulin retinol complex. Biochem Biophys Acta. 625: 28-42.
Ikonen T, Ojala M, Ruottinen O (1999). Associations between milk protein polymorphism and first lactation milk production traits in Finnish Ayrshire cows. J Dairy Sci. 82: 1026-1033.
Kaygisiz A, Douan M (1999). Genetics of milk protein polymorphism and its relation to milk yield traits in Holstein cows. Tr J Vet Anim Sci. 23: 447-454.
Kontopidis G, Holt C, Sawyer L (2004). Invited review: b-lactoglobulin: binding properties, structure, and function. J Dairy Sci. 87: 785-796.
Lin CY, Mcallister AJ (1986). Effects of milk protein loci on first lactation production in dairy cattle. J Dairy Sci. 69: 704-712.
Lum LS, Dovc P, Medrano JF (1997). Polymorphisms of bovine b-lactoglobulin promoter and differences in the binding affinity of activator protein-2 transcription factor. J Dairy Sci. 80: 1389-1397.
Lunden A, Nilsson M, Janson L (1997). Marked effect of b-lactoglobulin polymorphism on the ratio of casein to total protein in milk. J Dairy Sci. 80: 2996- 3005.
Meignanalakshmi A, Nainar AM, Nachimuthu K (2001). Identification of genetic polymorphism of b-lactoglobulin gene locus in Red-Sindhi cow by PCR-RFLP analysis. Int J Anim Sci. 6: 223-226.
Miller SA, Dykes DD, Polesky HF (1988). A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 16: 1215.
Ng-Kwai-Hang KF, Hayes JF, Moxley JE, Monardes HG (1984). Association of genetic variants of casein and milk serum proteins with milk, fat, and protein production by dairy cattle. J Dairy Sci. 67: 835-840.
Ng-Kwai-Hang  KF, Hayes JF, Moxley JE, Monardes HG (1986). Relationships between milk protein polymorphisms and major milk constituents in Holstein- Friesian cows. J Dairy Sci. 69: 22-26.
Ng-Kwai-Hang KF, Monardes HG, Hayes JF (1990). Association between genetic of milk proteins and production traits during three lactations. J Dairy Sci. 73: 3414-3420. 
Rachagani S, Dayal Gupta I, Gupta N, Gupta SC (2006). Genotyping of b-lactoglobulin gene by PCR-RFLP in Sahiwal and Tharparkar cattle breeds. BMC Genet. 7: 31.
Ron M, Yoffe O, Ezra E, Medrano JF, Welle JI (1994). Determination of effects of milk protein genotype on production traits of Israeli Holsteins. J Dairy Sci. 77: 1106-1113. 
Sabour MP, Lin CY, Keough A, Mechanda SM, Lee AJ (1993). Effects of selection practiced on the frequencies of k-casein and β-lactoglobulin genotypes in Canadian artificial insemination bulls. J Dairy Sci. 76: 274-280.
SAS® User’s Guide: Statistics, Version 8.2 Edition (2002). SAS Inst., Inc., Cary, NC.
Strazalkowska N, Kvzewski J, Ryniewiez Z, (2002). Effect of Kappa-casein and beta-lactoglobulin polymorphism on cows age, stage of lactation and somatic cell count on dairy milk composition in Polish Black and White cattle. Anim Sci. Papers and Reports 20: 21-35.
Tsiaras AM, Bargouli GG, Banos G, Boscos CM (2005). Effect of Kappa-casein and b-lactoglobulin loci on milk production traits and reproductive performance of Holstein cows. J Dairy Sci. 88: 327-334.
Wedholm A, Larsen LB, Lindmark-Mansson H, Karlsson AH, Andern A (2006). Effect of protein composition on the cheese-making properties of milk from individual dairy cows. J Dairy Sci. 89:3296-3305.
Yeh FC, Yang RC, Timothy BJ, Ye Z, Judy M (1997). Pop Gene, the user-friendly shareware for population genetic analysis. Molecular Biology and Biotechnology Center. Univ. Alberta.