Polymorphism of the b-Lactoglobulin Gene and Its Association with Milk Production Traits in Iranian Najdi Cattle

Document Type: Research Paper


Department of Animal Sciences, Faculty of Agricultural Sciences, Ramin Agricultural University, P.O. Box 62233-10544, Molsani, Ahvaz, I.R. Iran


The objective of this study was to determine the polymorphism of bovine β-lactoglobulin (β-LG) gene and its association with milk production traits in Iranian Najdi cattle. Blood samples were collected from 80 Najdi cattle from the Najdi cattle breeding station located in Shoshtar, Khuzestan. Polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) techniqus were used for identification of various genotypes. The results indicated that allele frequency of β-LG with regard to the B-allele (0.9125) was higher than that of the A-allele (0.0875). Frequencies of AB and BB genotypes were 0.175 and 0.825, respectively, while the AA genotype was found to be absent. Genotype frequencies were in accordance with the Hardy-Weinberg equilibrium. The effect of β-LG genotypes on milk production traits were analysed using a general linear model (GLM). There was no significant association between different genotypes of β-LG and milk production traits of the analysed cows.



The use of polymorphic genes as detectable molecular markers is a promising alternative to the current methods of trait selection once these genes are proven to be associated with traits of interest in animals. Selection efficiency, however, depends on allelic frequencies in the breeds and on the effect of these polymorphisms on selected traits (e.g. dairy traits and technological properties of milk). Among specific genes that may affect economically important traits in cattle, the b-lactoglobulin (b-LG) locus has been extensively studied (Tsiaras et al., 2005). β-LG is amphiphatic and an extremely acid stable protein which exists at the normal pH of bovine milk as a dimer with a molecular weight of 36,000 Daltons. It is a single chain polypeptide of 18 kDa comprising of 162 amino acid residues. The complete amino acid sequence of β-LG has been reported and genetic variation in amino acids sequence has been identified (Rachagani et al., 2006; Creamer et al., 1983). The biological functions of this protein are still not known. It could have a role in metabolism of phosphate in the mammary gland and the transport of retinol and fatty acids in the gut (Formaggioni et al., 1999; Hill, 1993).
The b-LG gene is situated on bovine chromosome 11 and encodes the main protein of whey. β-LG loci affect the milk production parameters and quality of milk protein. Their polymorphisms partly explain the genetic variance and improve the estimation of breeding value. Such loci can be taken into account as suitable supplements to conventional breeding procedures.  Polymorphism of this gene was discovered in 1955 and a total of 15 alleles are known so far (Matejicek et al., 2007). Common alleles are A, B, C and D, with alleles A and B being the most frequent. The point mutations in exon IV of the bovine b-LG gene determine two allelic variants A and B. The bovine b-LG A variant differs from the B variant by two amino acids only, i.e., aspartate-64 and valine-118. These amino acids are substituted by glycine and alanine, respectively in the B variant (Rachagani et al., 2006).
Effects of milk protein polymorphisms on milk production traits have been investigated during the past decades and, in some cases, results are still conflicting (Kucerova et al., 2006; Prinzenberg et al., 2003; Bovenhuis et al., 1992; Aleandri et al., 1990). Results of studies on the effect of b-LG genotypes on milk production traits have been rather more consistent. The AA genotype of b-LG has been shown to have a favorable effect on protein yield (Lunden et al., 1997; Sabour et al., 1996), whereas positive effects of the BB genotype on fat content have also been reported (Tsiaras et al., 2005; Ron et al., 1994; Hill, 1993). The b-LG locus affects mainly milk composition and milk quality and the B allele has been especially recognized as superior for milk quality in European cattle breeds, whereas allele A is associated rather with yield parameters (Strzalkowska et al., 2002). However, significant association between production traits and different genotypes of milk protein has not been reported in other researches (Curi et al., 2005; Ng-Kwai-Hang., 1998). Often alkaline and acidic polyacrylamide gel electrophoresis and PCR-based markers, especially PCR-RFLP, have been used for determination of different genotypes in such researches.
The Najdi breed of cattle (Bos indicus) represent native and low productive dairy cattle that are mainly bred in Khuzestan located in southern Iran, Resistance to diseases, adaptation to tropical climate and their role in the economy of rural families are important reasons for consideration of breeding and genetic improvement of this animal. The aim of this study was to identify genotypes of the b-LG gene and its effects on milk production traits in the Najdi cattle.


Blood samples were collected from 80 Najdi cattle at the Najdi cattle breeding station located in Shoshtar, Khuzestan. DNA was extracted from whole blood using the guanidin thiocyanate-silica gel method, as described by Boom et al. (1989). The gel monitoring and the spectrophotometeric methods were used for determination of DNA quality and quantity. The primers used for amplification of a 247 bp fragment of the b-LG gene were those described by Strzalkowska et al. (2002), with the following nucleotide sequence: 5´-TGTGCTGGACACCGACTACAAAAA-3´ (forward); 5´-GCTCCCGGTATATGACCACCCTCT-3´ (reverse).
Amplification reactions were carried out in a final volume of 25 ml, containing 100 ng of DNA, 0.5 mM of each primer, 1X PCR buffer, 1.5 mM MgCl2, 0.2 mM dNTPs and 1U of Taq DNA polymerase. The following cycles were applied: initial denaturation step at 95ºC for 5 min followed by 30 cycles : denaturation at 94ºC for 45 sec, primer annealing at 60ºC for 60 sec, PCR product synthesis at 72ºC for 60 sec and final synthesis step at 72ºC for 5 min. PCR products were recognized by electrophoresis on 1.2% agarose gel (0.48 g agarose was dissolved in 40 ml TBE 1X buffer), stained with ethidium bromide. The restriction digestion of the PCR products was carried out with the HaeIII enzyme. The digested products were loaded and visualized on 2.5% agarose gel (1 g agarose was dissolved in 40 ml TBE 1X buffer) agarose gel after staining with ethidium bromide. The allele and genotype frequencies were estimated by direct counting.
Records of milk yield (kg/day), fat and protein percentages were collected from the Najdi cattle station. Milk samples were analyzed monthly for fat and protein percentages. The time frame of all data collections is from 1997 to 2007. Detailed information was available for some reproductive events and was used to determine the traits of reproductive performance for each cow. Statistical analysis was carried out by means of the programme package SAS (SAS Inst., Inc., Cary, NC., U.S.A). The effect of the observed b-lactoglobulin genotypes on milk production traits were analysed using a general linear model (GLM). The following statistical model was used:
Yijklm = ì + Gi + Lj + Sk + Yl + eijklm

Yijklnm = records of analysed traits (milk yield, fat and protein percentages).
m = overall mean of population  
Gi = fixed effect of â-lactoglobulin genotypes ( i = 1, 2).
Lj = fixed effect of lactation number (j = 1, 2, …, 5; lactation numbers >5 were pooled with lactation numbers of 5).
Sk = fixed effect of calving season (k = 1, 2 for hot and cool seasons, respectively).
Yl =  fixed effect of calving year (l = 1, 2, …., 10 for 1997, 1998, ……., 2006, respectively).
eijklm =  random residual error.


The restriction digestion analysis of the 247 bp PCR product of the β-lactoglobulin gene indicated the presence of two types of restriction pattern; two fragments of 99, 74 bp (BB-genotype) and three fragments of 148, 99, 74 bp (AB-genotype), were observed. The restriction pattern that includes two fragments of 148, 99 bp (AA-genotype) was not found. In the analysed population of the Najdi cattle, the BB genotype was most frequent. The genotype frequencies of AA, AB and BB were 0, 0.175, and 0.825, respectively. Gene frequencies of the A and B alleles were 0.0875 and 0.9125, respectively. Expected frequency of the AA genotype was less than 5, therefore AA and AB genotypes were placed in the same category. Then chi-square (c2) test was performed on the basis of the Yates’ correction (Stansfield, 2004). Genotype frequencies were in accordance with the Hardy-Weinberg equilibrium (Table 1). The above results reveal that there was polymorphism in the β-LG gene, but that the B allele frequency was very much higher than the A allele thus indicating that genetic variability in this breed is low.
Milk from cows with the b-LG BB genotype had a greater fat percentage, while the milk yield and protein percentages of the b-LG AB genotype cows were higher, but this difference was not significant. There was no significant association between different genotypes of b-LG and milk production traits of the analysed cows. The means of the analysed traits for the observed b-LG genotypes are shown in Table 2.


The frequency of the A allele was found to be lower than that of the B allele in Najdi cattle. The frequencies of b-LG alleles obtained in this study are similar to those reported for Gyr, Nelore, Sindi (Del Lama and Zago, 1996), Sahiwal and Tharparkar breeds (Rachagani et al., 2006). These results mainly show that the frequency of A allele in the native breeds of cattle is lower than the European breeds of cattle. This study did not confirm the association between b-LG genotypes with milk production traits in Najdi cattle. In other studies of the b-LG genotypic effects on milk yield, several authors have reported no significant associations (Ghasemzadeh et al., 2007; Kucerova et al., 2006; Lunden et al., 1997; Ojala et al., 1997; Van Eenennaam and Medrano, 1991). However, reports exist where the b-LG genotypes AA (Kaminiski et al., 2002; Ikonen et al., 1999) or AB (Tsiaras et al., 2005) have been positively associated with milk yield. Several authors have reported significant effects of b-LG on the milk fat percentage and in all of these studies, the b-LG BB genotype was associated whit higher milk fat content (Bobe et al., 1999; Van der Berg et al., 1992). This inconsistency is also found in sheep studies. According to Dario et al. (2005) the breed differences, population size, frequency distributions of genetic variants, the structure of data analysed, the model used for statistical analysis and a failure to consider the relationship among animals are possible reasons for this inconsistency. Population size and structure of data analysed may further affect our findings. Hence it is essential that we improve methods of the standard recording of data. The limited number of animals and restricted methods applied so far did not help us to draw firm conclusions. Hence it should be verified in future studies with the more balanced populations of dairy cattle.

Aleandri R, Buttazzoni LG, Schneider JC, Caroli A, Davoli R (1990). The effects of milk protein polymorphisms on milk components and cheese-producing ability. J Dairy Sci. 73: 241-255.
Bobe G , Beitz DC , Freeman AE, Lindberg GL (1999). Effect of milk protein genotypes on milk protein composition and its genetic parameter estimates. J Dairy Sci. 82: 2797-2804.
Boom R, Sol CJ , Salimans MM (1989). Rapid and simple method for purification of nucleic acids. J Clin Microbiol. 28: 495-503.
Bovenhuis H, Van Arendonk AM, Korver S (1992). Associations between milk protein polymorphisms and milk production traits. J Dairy Sci. 75: 2549-2559.
Creamer L, Parry D, Malcolm G (1983). Secondary structure of b-lactoglobulin B. Arch Biochem Biophys. 227: 98-105.
Curi RA, Oliveria ND, Gimenes MA, Silveira AC (2005). Effect of CSN3 and LGB gene polymorphism on production traits in beef cattle. Genet Mol Biol. 28: 262-266.
Dario C, Carnicella D, Bufano G (2005). Effect of b-lactoglobulin genoypes on ovine milk composition in Altamurana breed. Arch Zootec. 54: 105-108.
Del Lama SN, Zago MA (1996). Identification for the kappa-caesin and beta-actoglobulin genotypes in Brazilian Bos indicus and Bubalis bubalis population. Brazilian J Genet. 19: 73-77.
Formaggioni P, Summer A, Malacarne M, Mariani P (1999). Milk protein polymorphism: Detection and diffusion of the genetic variants in Bos Genus. Ann Fac Med Vet Univ Parma. 19: 127-165.
Ghasemzadeh M, Rahimi GH, Khanmohammadi AR (2007). Detection of b lactoglobulin gene polymorphism in Montbeliard cattle and their association with milk yield using PCR-RFLP marker. In Proceedings of the 2th Iranian Conference of animal science, Iran, Tehran, Animal science society of Iran. 1252-1255.
Hill JP (1993). The relationship between b-lactoglobulin phenotype and milk composition in New Zealand dairy cattle. J Dairy Sci. 76: 281-286.
Kaminski S, Rymkiewicz-Schymczyk J, Wojcik E, Rusc A (2002). Associations between bovine milk protein genotypes and haplotypes and the breeding value of Polish Black-and-White bulls. J Anim Feed Sci. 11: 205-221.
Kucerova J, Matejicek A, Jandurova OM, Sqrensen P, Kott T, Frelich J (2006). Milk protein gene CSN1S1, CSN2, CSN3, LGB and their relation to genetic value of milk production parameters in Czech Fleckvieh. Czech J Anim Sci. 51: 241-247.
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.
Lunden A, Nilsson M, Janson L (1997). Marked effect of beta-Lactoglobulin polymorphism on the ratio of casein to total protein in milk. J Dairy Sci. 80 : 2996-3005.
Matejicek A, Matejickova J, Nemcova E, Frelich J (2007). Joint effect of CSN3 and LGB genotypes and their relation to breeding values of milk production parameters in Czech Fleckvieh. Czech J Anim Sci. 52:83-87.
Ng-Kwai-Hang KF (1998). Genetic polymorphism of milk proteins: Relationships with production traits, milk composition and technological properties. Can J Anim Sci. 78: 131-147.
Ojala M, Famula TR, Medrano JF (1997). Effects of milk protein genotypes on the variation for milk production traits of Holstein and Jersey cows in California. J Dairy Sci. 80: 1776-1785.
Pribyl J (1995). A way of using markers for farm animal selection. Czech J Anim Sci. 40: 375-382.
Prinzenberg EM, Gutscher K, Chessa S, Caroli A, Erhardt G (2003). Caprine k-Casein (CSN3) Polymorphism: New Developments in Molecular Knowledge. J Dairy Sci. 88: 1490-1498.
Rachagani S, Gupta ID, Gupta N, Gupta SC (2006). Genotyping of beta-lactoglobulin gene by PCR-RFLP in Sahiwal and Tharparkar cattle breeds. BMC Genet. 7: 31-34 .
Ron M, Yoffe O, Ezra E, Meldrano JF (1994). Determination of effect milk protein genotype on production traits of Israeli Holsteins. J Dairy Sci. 77: 1106-1113.
Sabour MP, Lin CY, Lee AJ, Mcallister AJ (1996). Association between milk protein genetics variants bulls for milk yield traits. J Dairy Sci. 79: 1050-1056.
Stansfield WD (2004). Theory and Problems of Genetic. MacGraw-Hill Book Company. NewYork. pp. 521
Strzalkowska N, Krzyzewski J, Zwierzchowski L, Ryniewicz Z (2002). Effect of kappa-casein and beta-lactoglobulin polymorphism cows age, stage of lactation and somatic cell count on daily milk yield and milk composition in Polish Black and white cattle. Anim Sci Pap Rep. 20: 21-35.
Tsiaras AM, Bargouli GG, Banos G, Boscos CM (2005). Effect of kappaCasein and betaLactoglobulin loci on milk production traits and reproductive performance. J Dairy Sci. 88: 324-334.
Van Eenennaam A, Medrano JF (1991). Milk protein polymorphisms in California dairy cattle. J Dairy Sci. 74: 1730-1742.
Van der Berg G, Escher JTM, Konning PJ, Bovenhuis H (1992). Genetic polymorphism of k-casein and b-lactoglobulin in relation to milk composition and processing. Netherl Milk Dairy J. 46: 145-168.