Open access

Milk Protein Genotype Associations with Milk Coagulation and Quality Traits

Written By

Elli Pärna, Tanel Kaart, Heli Kiiman, Tanel Bulitko and Haldja Viinalass

Submitted: August 2nd, 2011 Published: September 26th, 2012

DOI: 10.5772/29123

Chapter metrics overview

2,220 Chapter Downloads

View Full Metrics

1. Introduction

Cheese production is of substantial economic importance in most European countries, where an increasing amount of the produced milk is used for manufacturing cheese (Eurostat 2010a). In Estonia 60% (Statistics Estonia, 2010), in Italy more than 75% (De Marchi et al., 2008), and in Scandinavian countries 33% (Wedholm et al., 2006) of milk is used for cheese production. Recent trends indicate that per capita consumption of cheese is also increasing (Eurostat 2010b).

Milk quality is an essential factor to the dairy industry due to its economic impact. Milk coagulation ability is one of the most important factors affecting cheese yield and quality and has been reviewed (Jakob & Puhan, 1992; Johnson et al., 2001), and therefore is becoming more important as an increasing percentage of milk is used for cheese manufacturing. Milk coagulation properties (MCP) are commonly defined by milk coagulation time (RCT) and curd firmness (A30). It is feasible to design raw milk according to its specific technological use.

Improving cheese yield and quality, through the direct selection of breeding animals on the basis of milk coagulation property traits, is an option due to genetic variation of MCP traits (Ojala et al., 2005). MCP are heritable, quantitative traits; up to 40% of the variation among animals is caused by genetic factors (Ikonen et al., 2004). Estimates of heritability for MCP traits are from 0.30 to 0.40 (Bittante et al., 2002; Ikonen et al., 1999a), and from 0.25 for RCT, and 0.15 for A30 (Cassandro et al., 2008) to 0.28 for RCT to 0.41 for A30 (Vallas et al., 2010). Predictions of MCP provided by mid-infrared spectroscopy (MIR) techniques have been proposed as indicator traits for the genetic enhancement of MCP (Cecchinato et al., 2009; De Marchi et al., 2009). The expected response of RCT and A30 ensured by the selection using MIR predictions as indicator traits was equal to, or slightly less than, the response achievable through a single measurement of these traits. Accordingly, breeding strategies for the enhancement of MCP based on MIR predictions as indicator traits could be easily and immediately implemented for dairy cattle populations where the routine acquisition of spectra from individual milk samples is already measured (Cecchinato et al., 2009). Nevertheless, MCP traits analyzed with different methodologies have significantly different values, due to the diversity of the instruments used and the coagulant activity (Pretto et al., 2011). The type of coagulant could also have an effect, since different coagulants have been used. The method proposed for the prediction of non-coagulation probability of milk samples showed that non coagulating samples from one methodology were highly predictable based on the rennet coagulation time measured with another methodology (Pretto et al., 2011). A standard definition of MCP traits analysis is needed to enable reliable comparisons between MCP traits recorded in different laboratories, and in different animal populations and breeds.

More than 95% of the proteins contained in ruminant milk are coded by six structural genes (Martin et al., 2002). The four casein genes (CSN1S1, CSN2, CSN1S2, and CSN3) are linked in a cluster, referred to as the CN locus, mapped on chromosome 6 and encode αs1-CN, β-CN, αs2-CN, and κ-CN, respectively, as previously reviewed (Caroli et al., 2009). κ-CN is a key element in renneting, but interactions with the other milk protein systems have to be taken into account, in particular, β-CN and β-LG. Monitoring milk protein variation in different breeds of cattle avoids an increase of alleles with unfavourable effects on cheesemaking (Caroli et al., 2000; Comin et al., 2008; Erhardt et al., 1997; Ikonen et al., 1994, 1999a; Lodes et al., 1996). Therefore, another option for enchancing cheese yield and quality is indirect selection against, or for, some milk protein alleles. Selection against κ-CN E-allele would be a good means to indirectly improve milk quality for cheese production because the E-allele is unfavourably associated with non-coagulating milk, which is common (10%) in Finnish Ayrshire cows (Elo et al., 2007). Likewise the κ-CN G-allele in the Pinzgauer breed (Erhardt et al., 1997) and the κ-CN E-allele in the Italian Friesian breed (Caroli et al., 2000) were both associated with unfavourable coagulation properties. As for a positive association, it is well known that the B variants of β- and k-casein and β -lactoglobulin (β-LG) are favourable for milk coagulation and cheese-making (Dovc & Buchberger, 2000; Losi et al., 1973; Patil et al., 2003; Schaar et al., 1985; Walsh et al., 1998). k-CN allelic variants have been associated with variation in total casein and k-CN concentrations in milk (Hallén et al., 2008; Van den Berg et al., 1992), variation in casein micelle size (Di Stasio & Mariani, 2000; Walsh et al., 1998) and differences in coagulating properties and the cheesemaking quality of milk (Di Stasio & Mariani, 2000). Genetic variants of β -LG have been shown to have an indirect effect on cheese yield through their effect on the ratio of casein to total protein (Coulon et al., 1998; Lundén et al., 1997; McLean, 1986).

In the Dutch Holstein-Friesian population selection for CSN2-CSN3 haplotype A2-B, together with LGB B, would result in cows that produce milk more suitable for cheesemaking (Heck et al., 2009). It has also been argued that haplotypes have similar effects in the different breeds and the CN genes themselves were responsible for the haplotype effects observed, rather than genes physically linked to the CN complex (Boettcher et al., 2004).

The objective of this study was to estimate the contribution of the aggregate β-κ-CN and β-LG genotypes on first lactation milk coagulation and quality traits of Estonian Holstein cows. A parallel objective was to identify the variation in genetic polymorphism of milk proteins with the aim to improve the protein composition in milk by selecting for variants of specific genes.


2. Material and methods

2.1. Performance of Estonian cattle populations

In Estonia there are three breeds of dairy cattle – Estonian Red (ER), Estonian Holstein (EHF) and Estonian Native (EN). Distribution of breeds has come to favour EHF (Fig.1). In Estonia, 93.0% of cows are enrolled in an official milk recording programme (Results of Animal Recording in Estonia 2010, 2011) (Fig. 2). Since 1995, average milk yield in Estonia has risen 3947 kg (48%, Fig. 3). The 2010 average actual production for Estonian Holstein herd that were enrolled in producing-testing programs and eligible for genetic evaluation was 7778 kg milk, 317 kg of fat and 260 kg of protein per year. Holstein dairy cattle dominate in Estonian milk production industry because of their excellent production and greater income.

Figure 1.

Changes over time in number of cows of each indicated breed on a milk recording programme

Figure 2.

Changes over time in number of Estonian dairy cows and the proportion on a milk recording programme

In the evaluation programme for young bulls of Estonian Holstein breed, ca 25 bulls are tested each year, parallel testing is carried out on 10-12 foreign bulls. The selection of bulls is made from imported American and Canadian embryos, the Estonian Holstein best bull dams and imported young bulls. Often the sons of imported cows are used whose sires are world-known top bulls. The bulls come mainly from the USA, Canada, Germany and the Netherlands.

Currently, estimated breeding values (EBV) for production, conformation and udder health traits for bulls and cows in Estonia are computed by the Animal Recording Centre four times per year (Pentjärv & Uba, 2004). Breeding value estimation is carried out separately for the EHF and the Estonian Red breed (ER), using the best linear unbiased predictor (BLUP) test day animal model for production and udder health traits and the BLUP animal model for conformation traits. The EBV for each production trait – milk (kg), fat (kg) and protein (kg) – is the mean breeding value of the first, second and third lactations, adjusted by the mean average breeding value of the cows born in a defined base year (currently, 1995).

The milk production index (SPAV) is expressed as relative breeding value (RBV) with a mean of 100 and a standard deviation of 12 points for base animals, combining breeding values for milk, fat and protein yield weighted by relative economic values of 0:1:4 for EHF and 0:1:6 for ER (Pentjärv & Uba, 2004).

Figure 3.

Changes over time in annual milk yield per cow of each indicated breed

The information source for breeding value estimation of udder health traits is somatic cell count (SCC) in one millilitre of milk, transformed into a somatic cell score (SCS) using the formula SCS = log2 (SCC/100000) + 3 (Pentjärv & Uba, 2004).

The udder health index SSAV is calculated as the sum of EBVs of the first, second and third lactations with index weights 0.26, 0.37 and 0.37, respectively, and is expressed as RBV for genetic evaluation of conformation traits. Data from first lactation cows are used to compute RBVs for 16 linear traits for EHF and 14 linear traits for ER, as well as for three general traits. The conformation index SVAV is expressed as RBV, combining relative breeding values for type, udder and feet by relative economic weights of 0.3:0.5:0.2 for ER and 0.3:0.4:0.3 for EHF.

2.2. Data collection and laboratory milk analysis

First lactation milk samples were collected during routine milk recording as part of a development project for the Bio-Competence Centre of Healthy Dairy Products in Estonia during the period April 2005 – June 2008. The herds had twice-a-day or thrice-a-day milkings. The individual milk samples collected from the cows were either bulked test-day milkings or separate samples from each of the milkings on the test-day. Milk samples were immediately preserved after collection with Bronopol (2-bromo-2-nitropropane-1,3-diol, Knoll Pharmaceuticals, Nottingham, UK) and stored at 4˚C during transportation and analyzing periods. Milk samples with a pH lower than 6.5, indicative of colostrum, and non-coagulated milk samples (n=33) were excluded from the analysis. Furthermore, farms with less than 10 cows, and cows with fewer than three test-day records, were removed. The final dataset used for analyses consist of 11,437 test-day records from 2,769 Estonian Holstein cows which were located in 66 herds across the country and were daughters of 169 sires. The number of daughters per sire ranged from 1 to 267. Each cow had from 3 – 6 measurements collected during the different stages (7 – 305 DIM) of the first lactation. Information about the cows, herds and pedigree was obtained from the Estonian Animal Recording Centre (EARC), the Animal Breeders’ Association of Estonia and a database, COAGEN®, was produced. The test-day milk yield was recorded and individual milk samples were analyzed for fat percentage, protein percentage and urea content using the MilkoScan 4000 and MilcoScan FT6000, and for SCC using the Fossomatic 4000 and Fossomatic 5000 cell counter at the Milk Analysis Laboratory of the EARC, using methods suggested by the International Committee for Animal Recording (2009). Values of SCC were log-transformed to SCS.

The pH and milk coagulation properties were determined at the Laboratory of Milk Quality of the Estonian University of Life Sciences, usually three days after sampling. The proportion of milk samples with a maximum storage age of seven days was very small, less than 1%. The pH level of the milk was determined using a pH meter (Seven Multi; Mettler Toledo GmbH, Greifensee, Switzerland) at a temperature of 20˚C before analyzing the milk coagulation properties. The latter were milk coagulation time in minutes and firmness of curd in volts. Prior to the assessment of the milk coagulation properties, milk samples were heated to the renneting temperature (35˚C). The rennet (Milase MRS 750 IMCU/ml; CSK Food Enrichment B.V., The Netherlands) used in the analyses was diluted 1:100 (v/v) with distilled water and 0.2 ml of the solution was added to 10 ml milk. The milk coagulation properties were determined using the Optigraph (Ysebaert, Frepillon, France), which was developed by YDD (Ysebaert Dairy Division) in partnership with the INRA (LGMPA, lab. G. CORRIEU) to define coagulation characteristics in the laboratory, specifically to answer the needs of cheese makers (Ysebaert Dairy Division, 2009).

Measurements made with the Optigraph are not based on a rheological method but on an optical signal in the near-infrared spectrum. During a coagulation test, the light transmitted through the milk gradually weakens because of changes in the micellar structure of casein. The Optigraph then calculates the coagulation parameters (coagulation time, curd firmness, speed of aggregation) by means of particular feature points extracted from the optical information acquired in real time (Optigraph User’s Manual).

2.3. Data collection and laboratory blood analysis

Blood samples were collected as part of a development project for the Bio-Competence Centre of Healthy Dairy Products in Estonia during the period of June 2005 to December 2007. Blood samples (n=2,959) were stored in tubes containing K3EDTA. DNA was extracted from whole blood according to the method described by (Miller et al., 1998) or by using a commercial Puregene Gentra Blood kit (Minneapolis, USA). The quantity of template DNA was approximately 40 to 100 ng for Allele-specific oligonucleotide (ASO) PCR and PCR-RFLP, respectively. Polymorphisms of five milk protein genes were analyzed, four from the casein cluster (CSN1S1, csn1s2, CSN2, CSN3) and LGB. The list of single nucleotide polymorphisms (SNP) previously reported by (Chessa et al., 2007) was considered to distinguish genetic variants of milk proteins of the sampled cows. ASO primers were designed for the detection of polymorphisms in the CSN2 (primer sequences in (Värv et al., 2009)) and CSN1S1 (present study). CSN1S1 genotyping included amplification of a 155-bp sequence of the gene at exon 17. In accordance with the SNP a26181g in CSN1S1, two specific forward oligos were designed to distinguish B and C alleles paired with one reverse primer. An extra mismatch was added to both forward primers at position 2 at the 3’ end. Restriction analysis was carried out to genotype CSN1S2 (Ibeagha-Awemu et al., 2007), CSN3 (Velmala et al., 1993) and LGB (Medrano & Aquilar-Cordova, 1990). In this study, amplified regions of CSN1S2 with a digestion site to discriminate the protein variants from A-allele, were 330 bp (114 and 216 bp after restriction with MboII to detect B-allele), 354 bp (NlaIV restriction fragments 211 and 143bp to detect C-allele) and 356 bp (MnlI restriction fragments 160 and 196 bp for non-D allele). Products of allele-specific PCR and digestion fragments were separated on agarose gel. Sequencing to verify the studied DNA regions was performed with a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA) and analyzed using an ABI Prism 3130 Genetic Analyzer (Applied Biosystems, USA).

2.4. Statistical analyses

Preliminary analyses for testing the significance of fixed effects and single genotype effects were carried out on the SAS System® (SAS, Cary, NC, USA) using the MIXED procedure. Aggregate β-κ-CN genotypes were formed for further analysis. The genotypes with relative frequencies of less than 1% were grouped together into rare β-κ-CN genotype (A1A1-BB, A1A1-BE, A1A1-EE, A1B-BB, A1B-BE, A2A2-AE, A2A2-BE, A2A3-AA, A2B-AA, A2B-BB, BB-AB, BB-BB). Further statistical analysis was carried out using ASReml (VSN International Ltd., Hemel Hempstead, UK), using the following univariate repeatability animal model:


where y – vector of observations of the dependent variable (log-transformed RCT, A30, milk yield, milk protein percentage, fat percentage, SCS); b – vector of fixed effects (quadratic polynomial of day in milk, calving age, sample age, sampling year-season, calving year-season); g – fixed effects of the β-κ-CN genotypes or β-LG genotypes; a – vector of random effects (herdy=Xb+Cg+Za+e,, additive genetic N(0,Iσh2) and permanent environmental effectN(0,Aσa2)); e – vector of residual random error effectsN(0,Iσpe2); X, C, Z – known incidence matrices for fixed, genotype and random effects, respectively.

Sample age was included as a covariate in the model only for milk coagulation traits. Sampling year-season and calving year-season were grouped into 3-month classes, 14 classes from April 2005 to June 2008 and 11 classes from December 2004 to August 2007, respectively. Three generations of ancestors with a total number of 17,185 animals in the relationship matrix were included in the analysis.

The class with the largest number of observations, genotype A2A2 -AA was used as the class of comparison. It is also homozygous for both loci. Accordingly, the standard errors of the genotype effects are standard errors of the differences between each genotype and the most frequent A2A2 -AA genotype (Comin et al., 2008).


3. Results and discussion

3.1. Genetic structure of milk proteins

Allelic variants of casein and LGB loci and genotypes are presented in Tables 1 and 2. Four allelic variants of CSN2, three allelic variants of CSN3 and two variants of the LGB gene in Estonian Holstein population were detected. CSN2 occured at a significantly higher frequency for the A2 allele. The B-allele and A3 allele at CSN2 were rare. The further scan of β-CN, to discriminate protein variant I from variant A2 performed in the Dutch Holstein Friesian population revealed frequency of I allele of 0.14 for proven bulls, 0.27 in young bulls and 0.192 in cows and demonstrated that it is actually one of the common variants for the Holstein population (Visker et al., 2011). This frequency of the β-CN protein variant in Dutch Holstein Friesian population is relatively high compared with the frequencies of other cattle breeds. The Italian Holstein’s frequency of the I allele was 0.12 (Jann et al., 2002), while a survey of 30 cattle breeds yielded frequencies up to 0.14 (Jann et al., 2004). As pointed out (Visker et al., 2011), the associations of β-CN protein variant A2 should still interpreted with care, because they may have been combined associations of protein variants A2 and I, which are very different in some traits. Therefore further investigation is needed for Estonian Holstein β-CN I allele. The mean frequency of the β-CN A2 allele in Estonian Holstein was 0.647 (Table 1), ranging in farms in Estonia from 0.35 to 0.73. The mean frequency of the β-CN A1 allele was 0.320, ranging in farms in of Estonia from 0.19 to 0.50. At the CSN2 locus the A2 allele was highly predominant over other alleles. The heterozygous genotype A1A2 and the homozygous genotype A2A2 at β-CN were represented at almost the same frequencies (0.434 and 0.411, respectively). The most common β-CN genotypes A1A2 and A2A2 frequencies comprised jointly 84%. The occurrence of milk protein genetic variants revealed somewhat higher frequencies of the β-CN A2 allele and the homozygotic genotype A2A2 (Table 1 and 2) in the Estonian Holstein population in comparison with the common European cattle breeds (European Food Safety Authority, 2004). In some large Estonian farms the β-CN A2 allele mean frequency was about 0.70 and homozygotic genotype mean frequency about 0.50.

κ-CNA4,3560.737 αs1-CNB5,4460.983
E3910.066 αs2-CNA5,3290.962

Table 1.

Allele frequencies of κ-CN, β-CN, αs1-CN, αs2-CN and β-LG in cows of the Estonian Holstein breed

This finding, of a high frequency of the A2 allele, confirms the advantage of the Estonian Holstein breed that their milk naturally might lower health risks associated with the occurence of the β-CN A1 allele. The other advantage of the β-CN A2 allele is its positive association with protein yield (Olenski et al., 2010). The positive effect of the rare A2A2-BB genotype, that is the A2-B haplotype, on milk and protein yields has been reported in previous studies on Californian Holstein (Ojala et al., 1997), Finnish Ayrshire (Ikonen et al., 2001) and Dutch Holstein-Friesian cows (Heck et al., 2009).

AB 8500.288BC920.033
AE 2940.100 αs2-CNAA2,7640.926
BB 1160.039AB1590.057
β-CNA1A1 2700.092BD30.001
A1A21,2820.434 β-LGAA6310.214
A2B 1070.036

Table 2.

Genotype frequencies of κ-CN, β-CN, αs1-CN, αs2-CN and β-LG in cows of the Estonian Holstein breed

CSN3 had A, B and E allelic variants (Table 1 and 2). κ-CN shows a prevalence of the A allele at a mean frequency of 0.737, followed by the B allele at frequency of 0.197 and the E allele at frequency of 0.066. The most frequent κ-CN genotype of all genotyped Estonian Holstein cows was AA, which was found in slightly more than half of the cows (54.4%), followed by AB (28.8%) and AE (10%). A favourable genetic marker for protein yield, MCP and cheese production, κ-CN B, was rare (3.9% overall) in the homozygous state. The unfavourable κ-CN E allele was also very rare, only eight of the 2,954 sampled cows had the EE genotype.

CSN1S1 had B and C allelic variants (Table 1 and 2), showing the prevalence of the B allele (0.983) and CSN1S2 revealed three allelic variants: A, B, and D, showing the prevalence of the A allele (0.962). Because αs1-CN and αs2-CN were almost monomorphic (Table 1), they were excluded from the aggregate casein genotypes.

β-LG was represented with two allelic variants A (0.498) and B (0.502). Comparing the results of the genetic structure of milk proteins of the Estonian Holstein breed of Jõudu, 2008 and Värv et al., 2009 to those of this study, the quantity of detected alleles, and their frequencies, are somewhat different. The reason for this could be attributed to the sampling capacity (n=42) in (Värv et al., 2009) and the sampling procedure. Sampling was carried out in this study on a large population (n=2,954, 42 farms) across the whole breed. With a small sampling size the rare allele was not exposed and in one previous paper (Jõudu, 2008) the sampling was not performed across the whole breed, causing some different conclusions about allele frequencies of Estonian cattle breeds in their investigations compared to those of this study.

Expected frequencies of the β-κ-CN genotypes were calculated by multiplying the expected frequencies of the β-CN and κ-CN genotypes. Some alleles at one locus were associated only with certain alleles at the other locus, causing distinct differences between observed and expected frequencies of certain β-κ-CN genotypes (Table 3).

κ-CN β-CN genotype
AA4.0 (117)23.1 (683)0.027.1 (801)0.2 (5)0.1 (2)0.0
AB1.6 (48)10.2 (301)1.5 (45)12.5 (369)2.8 (82)0.1 (3)
AE2.5 (73)7.2 (214)0.00.2 (7)0.00.0
BB0.1 (3)1.2 (34)0.4 (12)1.2 (36)0.8 (24)0.2 (6)
BE0.7 (20)1.7 (51)0.3 (9)0.0 (1)0.0
EE0.3 (8)0.00.0

Table 3.

Observed and expected frequencies (upper and lower line respectively, each given as percentage) of the aggregate β-κ-CN genotypes (numbers of cows in the brackets) in 2,954 Estonian Holstein cows

Some genotypes were observed two to fourfold more frequently than expected (A1A1-BE, A1B-BE, A1B-BB etc) and A2B-BB at eightfold more frequently than expected. All cows carrying the κ-CN EE genotype had association only with the β-CN A1A1 genotype as has also been reported for Finnish Ayrshire cows (Ikonen et al., 1999a). Some genotypes were less frequent than expected (A1A1-BB a quarter and A2A2-AE one-twentieth of the expected frequency). Linkage disequilibrium in the casein loci has been observed in dairy cattle population differing in breed and geographical location, leading to unbalanced data (Bovenhuis et al., 1992; Ikonen et al., 1999b; Van Eenennaam & Medrano, 1991). The probable reason for unbalanced data in the Estonian Holstein population could be also the frequent use of few sires carrying, and transmitting, specific casein haplotypes to their offspring. Disequilibrium in the CSN2 and CSN3 loci can also be produced and maintained by selection favouring one combination of alleles over another (Falconer & Mackay, 1996).

3.2. Associations between milk protein genotypes and milk coagulation and quality traits

The associations of β-CN and κ-CN genotypes with milk coagulation (RCT, A30), quality traits (SCS, fat and protein contents), and milk yield was investigated (Table 4).

3.2.1. Milk coagulation traits

Milk coagulation traits (RCT and A30) were affected by aggregate β-κ-CN genotypes (p<0.001, Table 4). The most favourable β-κ-CN genotypes for RCT included the B allele at both loci, as has also been reported elsewhere (Comin et al., 2008) for Italian Holstein cows. Favourable aggregate genotypes for RCT were A1B-AB and A2B-AB. The best aggregate genotypes for A30 had two B alleles κ-CN, A1A2-BB, and the second best had the genotype A2A2-BB. κ-CN B was the most favourable for MCP, in every combination with β-CN, as has also been reported by (Comin et al., 2008).

RCT* (min)A30 (V)MILK (kg)PROTEIN ( %)FAT (%)SCS***
A1A1 -AA110-0.0370.014-0.5060.255-1.2370.4240.0080.0190.0410.0470.1410.139
A1A1 -AB42-0.0390.0212.1640.388-0.7030.6450.0730.0290.0080.071-0.0950.213
A1A1 -AE70-0.0330.016-0.8940.305-1.1320.5080.0530.0230.0350.0560.3110.168
A1A2 -AA633-0.0190.007-0.2370.130-0.6390.2170.0190.0100.0460.024-0.0490.072
A1A2 -AB271-0.0560.0092.4010.173-0.7080.2880.0540.0130.0640.0320.0720.095
A1A2 -AE207-0.0090.010-0.6830.1900.0280.317-0.0040.014-0.0350.0350.0510.105
A1A2 - BB34-0.0730.0234.3570.422-2.0130.7050.1300.0310.2460.077-0.5600.235
A1A2 - BE49-0.0520.0192.2140.355-0.9630.5950.0390.026-0.0080.0650.0840.197
A1B -AB42-0.1370.0212.3900.3900.0570.6480.0300.0290.0090.071-0.0810.214
A2A2 -AA768000000
A2A2 -AB337-0.0170.0092.1570.161-0.3830.2680.0570.0120.0360.029-0.1440.088
A2A2 - BB34-0.0720.0233.8650.418-0.4600.6990.0760.031-0.0420.077-0.2370.233
A2B -AB75-0.0920.0161.9660.295-0.4430.4930.0060.022-0.1210.054-0.0560.163

Table 4.

Statistical significance of milk protein genotypes (p), the number of Estonian Holstein cows (n) per β-κ-CN aggregate genotype and β-LG genotype, estimated genotype effects (Est.) with their standard errors (SE) on milk coagulation time (RCT), curd firmness (A30) and milk production and composition

As for the impact of CSN2 locus, the β-CN A2 allele was more favourable for A30 than A1. The A1 allele did not show, in this investigation, the most favourable effect (superior to β-CN A2) which was described by (Comin et al., 2008). The composite β-κ-CN genotypes, including the κ-CN B allele, were also associated with the best MCP in Finnish Ayrshire cattle (Ikonen et al., 1999a; Ojala et al., 2005), but the CSN2 locus in that sampled population did not include the β-CN B allele. Comparing both casein loci, it seems that CSN3 affected milk coagulation traits more than CSN2. Since the discovery of the micelle-stabilizing protein κ-casein, in 1956, it became evident that κ-CN had an important effect on the stabilization of casein micelles (Waugh & Von Hippel, 1956). In our study, the most frequent β-κ-CN genotype, A2A2-AA, and genotype A1A2-AE were associated with poor milk coagulation time, which is consistent with (Comin et al., 2008). Also, the rare E allele of the κ-CN in the β-κ-CN aggregate genotype had an unfavourable effect on milk coagulation properties. The association of this rare allele of κ-CN with poor milk coagulation has been previously reported (Ikonen et al., 1999a) and (Comin et al., 2008).

3.2.2. Milk yield and protein and fat percentage, SCS

Milk yield and protein and fat contents were affected by aggregate β-κ-CN genotype (p<0.05, Table 4). β-κ-CN genotype A1B-AB is favourable for milk yield. Similarly to A30, the most favourable β-κ-CN genotypes for milk protein content were homozygous for the B allele for κ-CN, A1A2-BB and A2A2-BB. The most favourable aggregate genotype for fat content was also A1A2-BB and unfavourable genotype for fat percentage, containing E allele in κ-CN locus, A1A2-AE, but also genotype A2B-AB.

The most favourable for protein content was BB for κ-CN and A1A2 for β-CN (the second best was A2A2, where the A2A2 genotype of β-CN had a slight advantage over the A1A1 genotype). These results were in agreement with those previously reported (Heck, 2009), that the κ-CN genotype was associated with protein content (B>A). Milk with the aggregate genotype A1A2-BB had the best firmness of curd and also the best protein and fat contents. This is in agreement with another investigation (Vallas et al., 2010), where curd firmness had the highest genetic correlation with milk protein percentage (0.48), suggesting that a high protein percentage results in a favourable curd firmness. It has been reported (Cassandro et al., 2008) that there is a correlation coefficient of 0.44 between curd firmness and protein percentage, which is in agreement with the results found in this experiment. The genetic correlations of −0.24 and −0.07 reported (Ikonen et al., 1999a, 2004) for the same traits, however, are different. These inconsistencies indicate that different methodologies used for the investigations may influence the results (Pretto et al., 2011). Curd firmness showed a weak positive genetic correlation (Vallas et al., 2010) with milk fat percentage (0.25) and a weak negative genetic correlation with milk yield (−0.29). Therefore, selection for improved curd firmness may be associated with a somewhat higher protein and fat percentage and slightly reduced milk yield. Genetic correlations for curd firmness with milk yield and fat percentage were negligible in previous studies (Cassandro et al., 2008; Ikonen et al., 1999a). As for the impact of the CSN2 A1 and A2 alleles on milk production, the A2 allele seems to have slight advantage over A1 in the aggregate β-κ-CN genotype. Genotypes of β-LG were associated with both milk coagulation traits (p<0.001), but had no a significant effect on either milk yield (p=0.462), protein percentage (p=0.648), nor fat percentage (p=0.356) and SCS (p=0.571).


4. Conclusion

The β-κ-CN locus had a strong effect on protein and fat content and milk coagulation properties. Milk with the β-κ-CN aggregate genotype A1A2-BB had the best firmness of curd and also the best protein and fat contents. The aggregate genotype A2A2-BB, haplotype A2-B, was also favourable for milk coagulation property traits and protein content. The β-LG locus had no impact on SCS, milk production nor protein and fat contents. The β-LG BB genotype had better curd firmness and AA better milk coagulation time. Linkage disequilibrium in the CSN2 and CSN3 loci, which most probably led to unbalanced data, provided justification for the use of aggregate β-κ-CN and β-LG in selection for better milk technological and quality traits.


The research leading to these results is co-financed by the European Community`s Regional Development Fund in the framework of the Competence Centre Programme of the Enterprise Estonia under project No EU22868; EU27789; EU28662; EU30002 of Bio-Competence Centre Of Healthy Dairy Products (Tervisliku Piima Biotehnoloogiate Arenduskeskus OÜ) and by the Targeted Finance Project 1080045s07.


  1. 1. BittanteG.MarusiM.CesariniF.PovinelliM.CassandroM.2002Genetic analysis on milk rennet-coagulation ability in Italian Holstein cows, Proceedings of the 7th World Congress on Genetics Applied to Livestock Production (WCGALP), CD-ROM Communication 09-030903978-2-73801-052-0Montpellier, France, August 19-23, 2002
  2. 2. BoettcherP. J.CaroliA.StellaA.ChessaS.BudelliE.CanavesiF.GhiroldiS.PagnaccoG.2004Effects of casein haplotypes on production traits in Italian Holstein and Brown cattle,Journal of Dairy Science8712December 2004), 431143170022-0302
  3. 3. BovenhuisH.Van ArendonkJ. A. M.KorverS.1992Associations between milk protein polymorphisms and milk production traits,Journal of Dairy Science739September 1992), 254925590022-0302
  4. 4. CaroliA.BollaP.BudelliE.BarbieriG.LeoneP.2000Effect of k-casein E allele on clotting aptitude of Italian Friesian milk, Zootecnica e Nutrizione Animale, 263June 2000), 1271300390-0487
  5. 5. CaroliA. M.ChessaS.ErhardtG. J.2009Invited review: Milk protein polymorphisms in cattle: Effect on animal breeding and human nutritionJournal of Dairy Science9211November 2009), 533553520022-0302
  6. 6. CassandroM.CominA.OjalaM.DalZotto. R.De MarchiM.GalloL.CarnierP.BittanteG.2008Genetic parameters of milk coagulation properties and their relationships with milk yield and quality traits in Italian Holstein cowsJournal of Dairy Science911January 2008), 3713760022-0302
  7. 7. CecchinatoA.De MarchiM.GalloL.BittanteG.CarnierP.2009Mid-infrared spectroscopy predictions as indicator traits in breeding programs for enhanced coagulation properties of milkJournal of Dairy Science9210October 2009), 530453130022-0302
  8. 8. ChessaS.ChiattiF.CeriottiG.CaroliA.ConsolandiC.PagnaccoG.CastiglioniB.2007Development of a Single Nucleotide Polymorphism Genotyping Microarray Platform for the Identification of Bovine Milk Protein Genetic Polymorphisms,Journal of Dairy Science901January 2007), 4514640022-0302
  9. 9. CominA.CassandroM.ChessaS.OjalaM.DalZotto. R.De MarchiM.CarnierP.GalloL.PagnaccoG.BittanteG.2008Effects of composite β- and κ-casein genotypes on milk coagulation, quality, and yield traits in Italian Holstein cows, Journal of Dairy Science, 9110October 2008), 402240270022-0302
  10. 10. CoulonJ.B.VerdierI.PradelP.AlmenaM.1998Effect of lactation stage on the cheesemaking properties of milk and the quality of Saint-Nectaire-type cheese,Journal of Dairy Research, 652May 1998), 2953050022-0299
  11. 11. De MarchiM.BittanteG.DalZotto. R.DalvitC.CassandroM.2008Effect of Holstein Friesian and Brown Swiss Breeds on Quality of Milk and CheeseJournal of Dairy Science9110October 2008), 409241020022-0302
  12. 12. De MarchiM.FaganC. C.O’DonnellC. P.CecchinatoA.DalZotto. R.CassandroM.PenasaM.BittanteG.2009Prediction of coagulation properties, titratable acidity, and pH of bovine milk using mid-infrared spectroscopyJournal of Dairy Science921January 2009), 4234320022-0302
  13. 13. Di StasioL.MarianiP.2000The role of protein polymorphism in the genetic improvement of milk productionZootecnica e Nutrizione Animale263June 2000), 69900390-0487
  14. 14. DovcP.BuchbergerJ.2000Lactoprotein genetic variants in cattle and cheese-making abilityFood Technology and Biotechnology382April-June 2000), 91981330-9862
  15. 15. EloK.TyriseväA.M.AnttilaP.VilvaV.OjalaM.2007Genomic mapping of non-coagulation of milk in the Finnish Ayrshire, Journal of Animal and Feed Sciences, 16No.: suppl. 1, 1951991230-1388
  16. 16. ErhardtG.PrinzenbergE. M.BuchbergerJ.Krick-SaleckH.KrauseI.MillerM.1997Bovine k-casein G detection, occurrence, molecular genetic characterization, genotyping and coagulation properties. Proceedings of the IDF Milk Protein Polymorphism Seminar II, International Dairy Federation, Brussels, 3283299-29098-026-9North, New Zealand, February 1997
  17. 17. European Food Safety Authority [EFSA].2004Opinion of the Scientific Panel on Dietetic Products, Nutrition and Allergies on a request from the Commission relating to the evaluation of goats’ milk protein as a protein source for infant formulae and follow-on formulae, TheEFSA Journal, 301151831-4732
  18. 18. Eurostat2010aProduction and external trade of foodstuffs: Dairy products and eggs, 27.06.2011, Available from
  19. 19. Eurostat2010bGross human apparent consumption of main food items per capita. Cheese, 27.06.2011, Available from
  20. 20. FalconerD. S.MackayT. F. C.1996Introduction to Quantitative Genetics,th ed. Longman Group Ltd., Essex, UK, 978-0-58224-302-6
  21. 21. HallénE.WedholmA.AndrénA.LundénA.2008Effect of β-casein, κ-casein and β-lactoglobulin genotypes on concentration of milk protein variants, Journal of Animal Breeding and Genetics, 1252April 2008), 1191290931-2668
  22. 22. HeckJ. M. L.2009Milk genomics, opportunities to improve the protein and fatty acid composition in raw milk,PhD thesis, Wageningen University, The Netherlands, 32 pp., 978-9-08585-332-9
  23. 23. HeckJ. M. L.SchenninkA.Van ValenbergH. J. F.BovenhuisH.ViskerM. H. P. W.Van ArendonkJ. A. M.Van HooijdonkA. C. M.2009Effects of milk protein variants on the protein composition of bovine milkJournal of Dairy Science923March 2009), 119212020022-0302
  24. 24. Ibeagha-AwemuE.PrinzenbergE.JannO.LuhkenG.IbeaghaA.ZhaoX.ErhardtG.2007Molecular characterization of bovine CSN1S2*B and extensive distribution of zebu-specific milk protein alleles in European cattle, Journal of Dairy Science907July 2007), 352235290022-0302
  25. 25. IkonenT.AhlforsK.KempeR.OjalaM.RuottinenO.1999aGenetic parameters for the milk coagulation properties and prevalence of noncoagulating milk in Finnish dairy cows,Journal of Dairy Science821January 1999), 2052140022-0302
  26. 26. IkonenT.BovenhuisH.OjalaM.RuottinenO.GeorgesM.2001Associations Between Casein Haplotypes and First Lactation Milk Production Traits in Finnish Ayrshire Cows,Journal of Dairy Science842February 2001), 5075140022-0302
  27. 27. IkonenT.MorriS.TyriseväA.M.RuottinenO.OjalaM.2004Genetic and phenotypic correlations between milk coagulation properties, milk production traits, somatic cell count, casein content, and pH of milk,Journal of Dairy Science872February 2004), 4584670022-0302
  28. 28. IkonenT.OjalaM.RuottinenO.1999bAssociations between milk protein polymorphism and first lactation milk production traits in Finnish Ayrshire cows,Journal of Dairy Science825May 1999), 102610330022-0302
  29. 29. IkonenT.OjalaM.SyväojaE. L.1997Effects of composite casein and β-lactoglobulin genotypes on renneting properties and composition of bovine milk by assuming an animal model, Agricultural and Food Science in Finland, 642832941239-0992
  30. 30. International Committee for Animal Recording [ICAR]. (2009). International agreement of recording practices. Guidelines approved by the General Assembly held June 18, 2008, 13.05.2010, Available from
  31. 31. JakobE.PuhanZ.1992Technological properties of milk as influenced by genetic polymorphism of milk proteins- A review, International Dairy Journal, 231571780958-6946
  32. 32. JannO.CeriottiG.CaroliA.ErhardtG.2002A new variant in exon VII of bovine β-CN gene (CSN2) and its distribution among European cattle breeds, Journal of Animal Breeding and Genetics, 1191February 2002), 65680931-2668
  33. 33. JannO. C.Ibeagha-AwemuE. M.ÖzbeyazC.ZaragozaP.WilliamsJ. L.Ajmone-MarsanP.LenstraJ. A.Moazami-GoudarziK.ErhardtG.2004Geographic distribution of haplotype diversity at the bovine casein locusGenetics Selection Evolution, 362March-April 2004), 2432570099-9193X
  34. 34. JohnsonM. E.ChenC. M.JaeggiJ. J.2001Effect of Rennet Coagulation Time on Composition, Yield, and Quality of Reduced-Fat Cheddar Cheese,Journal of Dairy Science845May 2001), 102710330022-0302
  35. 35. JõuduI.2008Effect of milk protein composition and genetic polymorphism on milk rennet coagulation properties, PhD thesis, Estonian University of Life Sciences, Tartu, 114 pp., 978-9-94942-654-6
  36. 36. LodesA.BuchbergerJ.KrauseI.AumannJ.KlostermeyerH.1996The influence of genetic variants of milk proteins on the compositional and technological properties of milk. 2. Rennet coagulation time and firmness of the rennet curdMilchwissenschaft5110December 1996), 5435480026-3788
  37. 37. LosiG.CastagnettiG. B.GraciaL.ZambonettiC.MarianiP.RussoV.1973Influenza delle varianti genetiche della caseina κ sulla formazione e sulle caratteristiche della cagliata, Scienza e Tecnologia degli Alimenti, 33733740304-0410
  38. 38. LundénA.NilssonM.JansonL.1997Marked effect of β-lactoglobulin polymorphism on the ratio of casein to total protein in milk, Journal of Dairy Science, 8011November 1997), 299630050022-0302
  39. 39. MartinP.SzymanowskaM.ZwierzchowskiL.LerouxC.2002The impact of genetic polymorphisms on the protein composition of ruminants milksReproductionNutrition Development, 425September-October 2002), 4334590926-5287
  40. 40. Mc LeanD. M.1986Influence of milk protein genetic variants on milk composition, yield and cheese making properties, Animal Genetics, 18suppl. 1, 1001020268-9146
  41. 41. MedranoJ.Aquilar-CordovaE.1990Polymerase chain reaction amplification of bovine β-lactoglobulin genomic sequences and identification of genetic variants by RFLP analysis, Animal Biotechnology, 1173771049-5398
  42. 42. MillerS. A.DykesD. D.PoleskyH. F.1998A simple salting procedure for extracting DNA from human nucleated cells, Nucleic Acids Research, 163February 1988), 12150305-1048
  43. 43. OjalaM.FamulaT. R.MedranoJ. F.1997Effects of milk protein genotypes on the variation for milk production traits of Holstein and Jersey cows in California,Journal of Dairy Science808August 1997), 177617850022-0302
  44. 44. OjalaM.TyriseväA.M.IkonenT.2005Genetic improvement of milk quality traits for cheese productionIn Indicators of Milk and Beef Quality, Ed. Hocquette, J.F. & Gigli, S.: Wageningen Academic Publishers, Wageningen, The Netherlands, 3073110071-2477
  45. 45. OlenskiK.KamińskiS.SzydaJ.CieslinskaA.2010Polymorphism of the beta-casein gene and its associations with breeding value for production traits of Holstein-Friesian bullsLivestock Science1311June 2010), 1371401871-1413
  46. 46. PatilM. R.BorkhatriyaV. N.BoghraV. R.SharmaR. S.2003Effect of bovine milk k-casein genetic polymorphs on curd characteristics during Cheddar cheese manufacture, Journal of Food Science and Technology406November 2003), 5825860022-1155
  47. 47. PentjärvA.UbaM.2004years of animal recording, In: Viinalass, H. (ed) Animal breeding in Estonia, 1416OÜ Paar, 9-94910-631-1
  48. 48. PrettoD.KaartT.VallasM.JõuduI.HennoM.AncilottoL.CassandroM.PärnaE.2011Relationships between milk coagulation property traits analyzed with different methodologies, Journal of Dairy Science, 949September 2011), 433643460022-0302
  49. 49. Results of Animal Recording in Estonia 2010.2011Estonian Animal Recording Centre, 52 pp., Elmatar, 0140-67341406734X, Tartu
  50. 50. SchaarJ.HanssonB.PetterssonH.E.1985Effects of genetic variants of k-casein and β-lactoglobulin on cheesemaking, Journal of Dairy Research, 5234294370022-0299
  51. 51. Statistics Estonia.2010Agriculture statistical database, Production of milk products (months), 27.06.2011, Available from
  52. 52. VallasM.BovenhuisH.KaartT.PärnaK.KiimanH.PärnaE.2010Genetic parameters for milk coagulation properties in Estonian Holstein cowsJournal of Dairy Science938August 2010), 378937960022-0302
  53. 53. Van denBerg. G.EscherJ. T. KoningP. J.BovenhuisH.1992Genetic polymorphism of k-casein and b-lactoglobulin in relation to milk composition and processing properties, Netherland Milk and Dairy Journal, 463-41451680002-8209X
  54. 54. Van EenennaamA. L.MedranoJ. F.1991Milk protein polymorphisms in California dairy cattle, Journal of Dairy Science, 745May 1991), 173017420022-0302
  55. 55. VärvS.BelousovaA.SildE.ViinalassH.2009Genetic diversity in milk proteins among Estonian dairy cattle, Veterinarija ir Zootechnika4870December 2009), 93981392-2130
  56. 56. VelmalaR.MäntysaariE. A.Mäki-TanilaA.1993Molecular genetic polymorphism at the k-casein and b-lactoglobulin loci in Finnish dairy bulls, Journal of Agricultural Sciences in Finland, 254314350782-4386
  57. 57. ViskerM. H. P. W.DibbitsB. W.KindersS. M.Van ValenbergH. J. F.Van ArendonkJ. A. M.BovenhuisH.2011Association of bovine β-casein protein variant I with milk production and milk protein composition, Animal Genetics, 422April 2011), 2122180268-9146
  58. 58. WalshC. D.GuineeT. P.RevilleW. D.HarringtonD.MurphyJ. J.O’KennedyB. T.FitzGerald. R. J.(1998Influence of k-casein genetic variant on rennet gel microstructure, Cheddar cheesemaking properties and casein micelle size,88August 1998), 7077140958-6946
  59. 59. WaughD. F.VonHippel. P. H.1956κ-Casein and stabilization of casein micelles, Journal of American Chemical Society, 7818September 1956), 457645820002-7863
  60. 60. WedholmA.LarsenL. B.Lindmark-MånssonH.KarlssonA. H.AndrénA.2006Effect of protein composition on the cheese making properties of milk from individual dairy cows, Journal of Dairy Science, 899September 2006), 329633050022-0302
  61. 61. Ysebaert Dairy Division, Optigraph,2009Available from

Written By

Elli Pärna, Tanel Kaart, Heli Kiiman, Tanel Bulitko and Haldja Viinalass

Submitted: August 2nd, 2011 Published: September 26th, 2012