Technical note: detection of the C allele of beta-casein (CSN2) in Czech Dairy goat breeds using LightCycler analysis

A protocol was developed for rapid genotyping of A and C variants at the CSN2 locus in goat species (White Shorthaired and Brown Shorthaired goat) by PCR and LightCycler analysis. The LightCycler technique combines rapid and efficient in vitro amplification of DNA in glass capillaries, with melting curve analysis based on fluorescence resonance energy transfer, for the sensitive detection of point mutation. Analysis of the CSN2 variability in the 2 goat breeds reared in the Czech Republic validated the genotyping test. Monitoring of CSN2 variability in the goat breeds indicates the predominance of the C allele. In both breeds, CSN2*A and CSN2*C showed almost similar frequencies. Variant CSN2*C occurred with a frequency of 0.699 in White Shorthaired goats and 0.570 in Brown Short-haired goats.     

Short communication: predominance of beta-casein (CSN2) C allele in goat breeds reared in Italy

A protocol for the rapid and simultaneous genotyping of A, C, and 0 'CSN2 alleles in goat was developed by single strand conformational polymorphism polymerase chain reaction (SSCP-PCR) technique. Screening the CSN2 variability in 7 goat breeds reared in Italy validated the genotyping test. The SSCP-PCR technique was also suitable for monitoring CSN2 polymorphism. In particular, the discrimination between CSN2*A and CSN2*C is important because the 2 corresponding protein variants cannot be separated by standard typing techniques. The monitoring of CSN2 variability in the goat breeds indicates the predominance of the C allele. In most breeds, CSN2*C occurred with the highest frequency, except in Saanen where CSN2*A and CSN2*C showed similar frequencies. Variant CSN2*C occurred with a frequency of 0.68 (Camosciata), 0.70 (Jonica), 0.71 (Garganica), 0.82 (Maltese), 0.87 (Cilentana), and 0.97 (Orobica). The alignment among the mature CSN2 sequences of different species suggests that CSN2*A is the ancestral allele compared with CSN2*C. Interestingly, the CSN2*A goat variant showed higher frequencies in selected breeds (Saanen and Camosciata).     

Technical note: simultaneous identification of CSN1S2 A, B, C, and E alleles in goats by polymerase chain reaction-single strand conformation polymorphism

Most variability in goat caseins originates from the high number of genetic polymorphisms often affecting the specific protein expression, with strong effects on milk composition traits and technological properties. At least 7 alleles have been found in the goat α_S2-CN gene (CSN1S2). Five of them (CSN1S2*A, CSN1S2*B, CSN1S2*C, CSN1S2*E, and CSN1S2*F) are widespread in most breeds, whereas the other 2 (CSN1S2*D and CSN1S2*0) are rarer alleles. Four different PCR-RFLP tests are needed to detect all of these variants at the DNA level. The objective of this study was to develop and validate a rapid method for typing 4 of the 5 most-common goat CSN1S2 alleles by means of PCR-single strand conformation polymorphism (SSCP). The method was validated by analyzing 37 goat samples at the protein and DNA level, respectively, by milk isoelectrofocusing and PCR-RFLP methods already described. The genotypes obtained using the PCR-SSCP approach were in full agreement with those obtained by the validation analyses. The newly developed PCR-SSCP approach provides an accurate and inexpensive assay highly suitable for genotyping goat CSN1S2.     

Genetic variation in the kappa-casein gene (CSN3) of Chinese yak (Bos grunniens) and phylogenetic analysis of CSN3 sequences in the genus Bos

Variants of κ-casein (CSN3) have been extensively studied in cattle and 13 alleles have been identified at the protein and DNA levels to date. Evolution of some of these alleles and a possible common ancestor remain unclear. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) analysis of CSN3 exon IV in domesticated yak revealed a 2-allele polymorphism showing migration patterns different from known cattle variants. The PCR products of both yak CSN3 SSCP alleles were sequenced. All yak had nucleotide sequences coding for Thr in AA position 136 (identical to bovine CSN3*A) and Ala in position 148 (identical to bovine CSN3*B). The sequencing results were confirmed by PCR-RFLP analysis using HindIII and TaqI. A 12-bp insertion in the coding region, representing a repeated nucleotide and AA motif, was found in 1 yak allele. The duplication corresponds to the codons for AA 147 to 150 (Glu-Ala-Ser-Pro) or 148 to 151 (Ala-Ser-Pro-Glu), which are repeated identically. In 21 yak samples genotyped by PCR-SSCP analysis, frequencies for the insertion variant and the short variant were about 68 and 32%, respectively. The loss of the insertion may have led to the ancestral CSN3 allele from which all currently known variants of CSN3 in the genus Bos evolved. This is the first report of polymorphisms in the yak CSN3 gene and may be helpful for future studies on genetic variation within and between yak populations or on associated traits.     

High polymorphism in the kappa-casein (CSN3) gene from wild and domestic caprine species revealed by DNA sequencing

We assessed polymorphisms in exon IV of the kappa-casein gene (CSN3) in ten different breeds of domestic goat (Capra hircus) from three continents and in three related wild caprine taxa (Capra ibex, Capra sibirica and Capra aegagrus). Thirty-five DNA samples were sequenced within a 558 bp fragment of exon IV. Nine polymorphic sites were identified in domestic goat, including four new polymorphisms. In addition to four previously described polymorphic positions, a total of 13 polymorphisms allowed the identification of 13 DNA variants, corresponding to 10 protein variants. Because of conflicting nomenclature of these variants, we propose a standardized allele designation. CSN3*A, CSN3*B, and CSN3*D were found as widely distributed alleles in European goat breeds. Within Capra ibex we identified three variants and showed that the sequence of Capra aegagrus is identical to the most common Capra hircus variant, consistent with Capra aegagrus being the wild progenitor of domestic goats. A dendrogram was drawn to represent the molecular network between the caprine CSN3 variants.     

Short communication: the beta-casein (CSN2) silent allele C1 is highly spread in goat breeds

Several single nucleotide polymorphisms have been identified in the goat milk casein genes, most of them modifying the amino acid sequence of the coded protein. At least 9 variants have been found in goat β-CN (CSN2); 6 of them were characterized at the DNA level (A, A1, C, E, 0, and 0′), whereas the other 3 variants were described only at the protein level. The recently identified silent A1 allele is characterized by a C→T transition at the 180th nucleotide of the ninth exon. In the present work, typing results from different breeds (3 Italian, 3 German, and a composite of African breeds for a total of 335 samples) demonstrated that the same mutation is carried by the CSN2*C allele. In addition, the T nucleotide at the 180th nucleotide of the ninth exon was always associated with CSN2*C in all the breeds analyzed. Thus, another silent allele occurs at goat CSN2 and can be named CSN2*C1. The much wider distribution of C1 with respect to the A1 allele indicates that the single nucleotide polymorphisms characterizing the silent mutation originated from CSN2*C. A method for the identification of this allele simultaneously with 5 of the 6 DNA-characterized alleles is also proposed. The mutation involved codifies for the same protein of the C allele; nevertheless, its location in the 3′ untranslated region of the gene might affect the specific casein expression.     

Molecular characterization of bovine CSN1S2*B and extensive distribution of zebu-specific milk protein alleles in European cattle

The B allele of the bovine α_S2-casein gene (CSN1S2) was characterized at the molecular level and the distribution of zebu-specific milk protein alleles was determined in 26 cattle breeds originating from 3 continents. The CSN1S2 *B allele is characterized by a C → T transition affecting nucleotide 17 of exon 3, which leads to a change in the eighth amino acid of the mature protein, from Ser to Phe (i.e., TCC →TCC). DNA-based methods were developed to identify carriers of CSN1S2*B and the other alleles (CSN1S2 *A, C, and D) at the same locus. CSN1S2*B and other zebu-specific milk protein alleles and casein haplotypes are widely distributed in European cattle breeds, particularly those of southeastern origin. Alleles CSN1S2 *B and CSN3*H are important in searching for zebu imprints in European cattle breeds. Diversity estimates at the milk protein loci were highest in the zebus followed by southeastern European taurines. Anatolian Black had the highest number of zebu alleles among European taurines. Common, group, and intergroup haplotypes occurred in the breeds and demonstrated relationships that concurred with developmental histories, genetic makeup, and, in particular, exposed the extent of zebu influence on southeastern European cattle.     

Caprine kappa-casein (CSN3) polymorphism: new developments in molecular knowledge

A high degree of polymorphism was recently found at the kappa-casein (CSN3) locus in the domesticated goat (Capra hircus). In the present study, 2 new patterns previously identified by PCR-single-strand conformation polymorphism analysis (SSCP) were characterized. The allele provisionally named "X" (GenBank Accession no. AY350425) differs from CSN3*C (AF485341) by a (silent) A-->G substitution at position 509 of the goat CSN3 reference sequence (X60763). As this newly identified sequence changes the amino acid sequence, and the already known CSN3*C allele (AF485341) has an additional silent mutation, we proposed a change in nomenclature to reflect these changes, indicating the silent mutation with the prime symbol (i.e.,'). The CSN3*M allele (provisionally named "Y") results in a new protein variant, differing by 2 nonsynonymous mutations from the CSN3*F allele. The new variant is characterized by a G-->A transition at nucleotide position 384, resulting in the amino acid exchange Asp90-->Asn90, and a C-->T transition at position 550, resulting in a Val145-->Ala145 substitution. Thus, the number of alleles identified in the domesticated goat has increased to 16, of which 13 are protein variants and 3 are silent mutations, involving a total of 15 polymorphic sites in CSN3 exon 4. Data on the distribution of the main alleles in 7 goat breeds of Europe, West Africa, and the Near East show differences in the occurrence and frequency of the alleles between breeds and geographic origin with the highest number of alleles found in goat breeds from the Near East.     

Short communication: new alleles of the bovine kappa-casein gene revealed by resequencing and haplotype inference analysis

We tested the hypothesis that extensive undiscovered genetic diversity exists in important functional genes from domestic and wild cattle species (Bos spp.). We resequenced 483 bp of a key exon (exon IV) from the kappa (κ)-casein gene (CSN3) for a panel of samples of domestic cattle from 8 countries and a close relative species, the gayal (Bos frontalis). Six single nucleotide polymorphisms were identified. Haplotype inference revealed 12 haplotypes, of which 8 were newly discovered. Among these 8 new haplotypes, 5 differed by one nonsynonymous mutation and 3 differed by one silent mutation from previously well-characterized CSN3 alleles. From those, one was shared by the gayal and Zebu, was different from CSN3*B at position Ile136Thr, and showed a close phylogenetic relationship with the banteng, gaur, and yak. The other 7 new haplotypes were detected in our panel of worldwide local cattle breeds but were absent from previously reported commercial breeds. These results support the hypothesis that genetic diversity at the coding region of CSN3 has been underestimated. This study also highlights how important it is to resequence functionally important genes in worldwide local cattle breeds, many of which are threatened by extinction or replacement by commercial breeds.     

Short communication: Polymorphism of casein cluster genes in Czech local goat breeds

The 4 casein loci were evaluated as haplotypes in 2 dairy goat breeds kept in the Czech Republic. Analysis of 41 White Shorthaired (WSH) trio families and 44 Brown Shorthaired (BSH) trio families revealed 14 and 20 haplotypes, respectively. Various genomic techniques were used to type the casein loci. Twenty-two different combinations of these alleles (casein haplotypes, in the order CSN1S1-CSN2-CSN1S2-CSN3) were found. Only 5 haplotypes in the WSH breed and 6 haplotypes in the BSH breed occurred at frequencies >0.05. For the WSH breed, the most common haplotype was FCFB (0.260), whereas for the BSH breed, the most common haplotype was FCFA (0.217). The information on the haplotype variability in both breeds could be used in breeding programs aimed at preserving biodiversity or selecting animals for specific protein production and cheesemaking.     

Genomic analysis of the origins of extant casein variation in goats

The variation in the casein genes has a major impact on the milk composition of goats. Even though many casein polymorphisms have been identified so far, we do not know yet whether they are evolutionarily ancient (i.e., they existed before domestication) or young (i.e., they emerged after domestication). Herewith, we identified casein polymorphisms in a data set of 106 caprine whole-genome sequences corresponding to bezoars (Capra aegagrus, the ancestor of domestic goats) and 4 domestic goat (Capra hircus) populations from Europe, Africa, the Far East, and the Near East. Domestic and wild goat populations shared a substantial number of casein SNP, from 36.1% (CSN2) to 55.1% (CSN1S2). The comparison of casein variation among bezoars and the 4 domestic goat populations demonstrated that more than 50% of the casein SNP are shared by 2 or more populations, and 18 to 44% are shared by all populations. Moreover, the majority of casein alleles reported in domestic goats also segregate in the bezoar, including several alleles displaying significant associations with milk composition (e.g., the A/B alleles of the CSN1S1 and CSN3 genes, the A allele of the CSN2 gene). We conclude that much of the current diversity of the caprine casein genes comes from ancient standing variation segregating in the ancestor of modern domestic goats.     

Characterization and genotyping of the caprine kappa-casein variants

Kappa-Casein (kappa-CN) is the milk protein that determines the size and specific function of milk micelles, and its cleavage by chymosin is responsible for milk coagulation. We have previously detected and characterized four variants of the goat kappa-CN in Spanish, French, and Italian breeds by screening the major part of the coding region in exon 4. Here we have sequenced and analyzed the full coding region of the kappa-CN gene which includes exons 3 and 4. No additional mutations were found, with exception of a single nucleotide substitution in exon 3, which had no amino acid change. However, the analysis of the association between the different mutations resulted in two new variants designated kappa-CN F and G. The novel variants are present in the Italian breeds Teramana, Girgentana, and Sarda (variant F). A protocol for rapid simultaneous genotyping of all known kappa-CN variants using the primer extension method was described, and a total of 210 animals from nine European breeds were genotyped. Alleles A and B are the most frequent variants occurring in the majority of breeds with highest prevalence of the B variant, except for the Canaria breed where the A allele is more frequent. Sequence data suggest that the F variant is the original type of caprine kappa-CN, other alleles being derived from this type following two different trunks by successive mutations.     

Effect of kappa-casein polymorphism on milk composition in the Orobica goat

The aim of this work was to study the effects of isoelectrofocusing (IEF) milk protein variants on milk composition in the Italian Orobica goat breed, which is characterized by a rather high frequency of the κ-casein (CSN3) B^IEF allele. Significant associations were found between the IEF phenotype and protein and casein percentages. A favorable effect of the CSN3 B^IEF variant was found for both protein and casein percentages, with a codominance trend for the 3 phenotypes: BB > AB > AA. Depending on the selection purpose, emphasis could be given to different κ-casein variants in breeding. The high frequency of B^IEF could be exploited in breeding strategies to improve the protein and casein percentages when cheese making is a selection objective.     

Characterization of the casein gene complex in West African goats and description of a new alpha(s1)-casein polymorphism

The analysis of casein polymorphisms was carried out in West Africa goat populations: Red Sokoto (n = 57), West African Dwarf Nigeria (n = 27), West African Dwarf Cameroon (n = 39), and Borno (n = 37). The 4 casein genes α_s1 (CSN1S1), β (CSN2), α_s2 (CSN1S2), and κ (CSN3) were typed at the DNA level. No null alleles were found in any of the genes analyzed. A PCR single-strand conformation polymorphism method was implemented for the identification of CSN1S1*F allele simultaneously with A/0₁, B/E, N and the new allele. The allele differed from CSN1S1*B by a synonymous transversion TCG→TCT in the codon corresponding to Ser₆₆ of the mature protein. The new allele, named CSN1S1*B′, occurred at a high frequency in all the populations, ranging from 0.295 (West African Dwarf Cameroon) to 0.405 (Borno). A greater frequency was found for alleles associated with high α_s1-casein quantity, as has already been observed in the goat populations from the Mediterranean area. The intermediate E allele occurred only in the Red Sokoto and at a low frequency. The faint F allele occurred in 3 populations at frequencies lower than 0.03. Linkage disequilibrium occurred in all the populations, with highly significant differences in Borno, Red Sokoto, and West Africa Dwarf Nigeria, and significant differences in West Africa Dwarf Cameroon. Only 10 haplotypes showed frequencies ≥0.05 in at least 1 of the 4 populations considered, and the overall frequency was >0.1 only for 4 haplotypes: BAAB, B′ACA, ACAB, and BACA (in the order CSN1S1-CSN2-CSN1S2-CSN3). Haplotype BAAB, postulated as an ancestral haplotype in previous studies, was the most common haplotype in all breeds except Borno, where B′ACA was predominant. The results obtained are of considerable significance given that very little information exists on the subject for African goats. The high frequency of strong alleles in the calcium-sensitive caseins as well as the high linkage disequilibrium found among the casein genes in the African breeds analyzed may suggest that specific casein haplotypes have already been selected due to their advantages for nutrition. Haplotypes providing greater protein and casein content would increase the energy content of milk, thus resulting in more favorable growth and survival of young goats and humans consuming the milk.     

Short communication: casein haplotype variability in sicilian dairy goat breeds

In the Mediterranean region, goat milk production is an important economic activity. In the present study, 4 casein genes were genotyped in 5 Sicilian goat breeds to 1) identify casein haplotypes present in the Argentata dell’Etna, Girgentana, Messinese, Derivata di Siria, and Maltese goat breeds; and 2) describe the structure of the Sicilian goat breeds based on casein haplotypes and allele frequencies. In a sample of 540 dairy goats, 67 different haplotypes with frequency ≥0.01 and 27 with frequency ≥0.03 were observed. The most common CSN1S1-CSN2-CSN1S2-CSN3 haplotype for Derivata di Siria and Maltese was FCFB (0.17 and 0.22, respectively), whereas for Argentata dell’Etna, Girgentana and Messinese was ACAB (0.06, 0.23, and 0.10, respectively). According to the haplotype reconstruction, Argentata dell’Etna, Girgentana, and Messinese breeds presented the most favorable haplotype for cheese production, because the casein concentration in milk of these breeds might be greater than that in Derivata di Siria and Maltese breeds. Based on a cluster analysis, the breeds formed 2 main groups: Derivata di Siria, and Maltese in one group, and Argentata dell’Etna and Messinese in the other; the Girgentana breed was between these groups but closer to the latter.     

New κ-Casein Alleles in Jakhrana Goat Affecting Milk Processing Properties

High genetic variability in the κ-Casein gene that can be exploited for the marker assisted selection milk production of desirable processing and nutritional properties. We resolved seven unique PCR-SSCP variants in a 425 bp segment of κ-Casein exon-4 (CSN*3) on PAGE gel and then sequenced the amplicons. All SSCP haplotypes revealed novel sequences (NCBI Genbank accessions EF053350- EF053356). Nucleotide substitutions at nucleotide positions 245, 247, 274, 345 and 471 caused single nucleotide polymorphism (SNP). We identified seven new alleles CSN3*N, CSN3*O, CSN3*P, CSN3*Q and CSN3*R being homozygous and CSN3*A^H, CSN3*K^H being heterozygous types. In silico analysis yielded an isoelectric point of 5.29 and 5.63 for the conceptual proteins based on CSN3*A^H and CSN3*K^H, CSN3*N, CSN3*O, CSN3*P, CSN3*Q, CSN3*R allele sequence, respectively, which were of IEF_A and IEF_B type κ-Casein (CSN3).     

Casein haplotype structure in five Italian goat breeds

The aim of this work was to investigate the genetic structure of the casein gene cluster in 5 Italian goat breeds and to evaluate the haplotype variability within and among populations. A total of 430 goats from Vallesana, Roccaverano, Jonica, Garganica, and Maltese breeds were genotyped at alphas1-casein (CSN1S1), alphas2-casein, (CSN1S2), beta-casein (CSN2), and kappa-casein (CSN3) loci using several genomic techniques and milk protein analysis. Casein haplotype frequencies were estimated for each breed. Principal component analysis was carried out to highlight the relationship among breeds. Allele and haplotype distributions indicated considerable differences among breeds. The haplotype CSN1S1*F- CSN1S2*F-CSN3*D occurred in all breeds with frequencies >0.100 and was the most common haplotype in the Southern breeds. A high frequency of CSN1S1*0-CSN1S2*C-CSN3*A haplotype was found in Vallesana population (0.162). Principal component analysis clearly separated the Northern and Southern breeds by the first component. The variability of the caprine casein loci and variety of resulting haplotypes should be exploited in the future using specific breeding programs aiming to preserve biodiversity and to select goat genetic lines for specific protein production.     

Focusing on the goat casein complex

The analysis of casein polymorphisms in goat species is rather difficult, because of a large number of mutations at each locus, and the tight linkage involving the 4 casein genes. Three goat breeds from Northern Italy, Orobica, Verzasca, and Frisa, were analyzed at the casein complex by milk isoelectrofocusing and analyses at the DNA level to identify the majority of all known polymorphisms. The casein gene structure of the 3 local breeds at α_S1-casein (CSN1S1), β-casein (CSN2), α_S2-casein (CSN1S2), and κ-casein (CSN3) was compared with that of Camosciata, a more widely distributed breed. A new allele was identified and characterized at CSN2 gene, which seemed to be specific to the Frisa breed. It was named CSN2*E, and was characterized by a transversion TCT → TAT responsible for the amino acid exchange Ser₁₆₆ → Tyr₁₆₆ in the mature protein. The casein haplotype structure is highly different among breeds. A total of 26 haplotypes showed a frequency higher than 0.01 in at least 1 of the 4 breeds considered, with 12, 3, 5, and 19 haplotypes in Frisa, Orobica, Verzasca, and Camosciata breeds, respectively. Only 13 haplotypes occurred at a frequency higher than 0.05 in at least 1 breed. With the molecular knowledge of each locus, the ancestral haplotype coding for CSN1S1*B, CSN2*A, CSN1S2*A, and CSN3*B protein variants can be postulated. A protein evolutionary model considering the whole casein haplotype is proposed.     

Single nucleotide polymorphisms in the ovine casein genes detected by polymerase chain reaction-single strand conformation polymorphism

Casein genetic polymorphisms are important and well known due to their effects on quantitative traits and technological properties of milk. At the DNA level, polymerase chain reaction (PCR)-single-strand conformation polymorphism (SSCP) allows for the simultaneous typing of several alleles at casein loci, as well as the detection of unknown polymorphisms. Here we describe the usefulness of the PCR-SSCP technique for casein typing in sheep. In particular, three single-nucleotide polymorphisms (SNP) are described at CSN1S1, CSN2, and CSN3, all resulting in amino acid exchanges. At CSN1S1, a transition T-->C was found, resulting in the deduced amino acid exchange Ile186-->Thr186. A transition A-->G resulting in the deduced amino acid exchange Met183-->Val183 was identified at CSN2. The 2 SNP showed a rather high frequency (ranging from 0.12 to 0.26) in 3 Italian breeds (Sarda, Comisana, Sopravissana). Another transition C-->T (Ser104-->Leu104) was found at CSN3 in one heterozygous animal.     

The occurrence of noncoagulating milk and the association of bovine milk coagulation properties with genetic variants of the caseins in 3 Scandinavian dairy breeds

Substantial variation in milk coagulation properties has been observed among dairy cows. Consequently, raw milk from individual cows and breeds exhibits distinct coagulation capacities that potentially affect the technological properties and milk processing into cheese. This variation is largely influenced by protein composition, which is in turn affected by underlying genetic polymorphisms in the major milk proteins. In this study, we conducted a large screening on 3 major Scandinavian breeds to resolve the variation in milk coagulation traits and the frequency of milk with impaired coagulation properties (noncoagulation). In total, individual coagulation properties were measured on morning milk collected from 1,299 Danish Holstein (DH), Danish Jersey (DJ), and Swedish Red (SR) cows. The 3 breeds demonstrated notable interbreed differences in coagulation properties, with DJ cows exhibiting superior coagulation compared with the other 2 breeds. In addition, milk samples from 2% of DH and 16% of SR cows were classified as noncoagulating. Furthermore, the cows were genotyped for major genetic variants in the α_S1- (CSN1S1), β- (CSN2), and κ-casein (CSN3) genes, revealing distinct differences in variant frequencies among breeds. Allele I of CSN2, which had not formerly been screened in such a high number of cows in these Scandinavian breeds, showed a frequency around 7% in DH and DJ, but was not detected in SR. Genetic polymorphisms were significantly associated with curd firming rate and rennet coagulation time. Thus, CSN1S1 C, CSN2 B, and CSN3 B positively affected milk coagulation, whereas CSN2 A², in particular, had a negative effect. In addition to the influence of individual casein genes, the effects of CSN1S1-CSN2-CSN3 composite genotypes were also examined, and revealed strong associations in all breeds, which more or less reflected the single gene results. Overall, milk coagulation is under the influence of additive genetic variation. Optimal milk for future cheese production can be ensured by monitoring the frequency of unfavorable variants and thus preventing an increase in the number of cows producing milk with impaired coagulation. Selective breeding for variants associated with superior milk coagulation can potentially increase raw milk quality and cheese yield in all 3 Scandinavian breeds.