Aberrant low expression level of bovine beta-lactoglobulin is associated with a C to A transversion in the BLG promoter region

Using QTL to Identify Genes and Pathways Underlying the Regulation and Production of Milk Components in Cattle

Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk composition, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focuses primarily on QTL discovered in cattle (Bos taurus) as a model species for the biology of lactation, and there are occasional references to sheep genetics. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype databases continue to grow and diversify, new QTL will continue to be discovered, and although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of the biology of lactation.

Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition

Background

Although bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes.

Results

Genome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein β-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and β-lactoglobulin protein levels. β-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield.

Conclusions

The genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00731-9.

Screening for phenotypic outliers identifies an unusually low concentration of a β-lactoglobulin B protein isoform in bovine milk caused by a synonymous SNP

Background

Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, β-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration.

Results

A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B’) on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B’ allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B’ allele which were homozygous for the BLG B protein isoform (BB’). The level of expression of the BLG B’ allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B’B’), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B’ downregulates milk BLG concentration remains to be determined.

Conclusions

High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B’), which reduced the milk concentration of β-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B’ allele is expected to have improved processing characteristics, particularly for cheese-making.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00711-z.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.

Genes encoding equine β-lactoglobulin (LGB1 and LGB2): Polymorphism, expression, and impact on milk composition

β-lactoglobulin is one of the most abundant milk whey proteins in many mammal species, including the domestic horse. The aim of this study was to screen for polymorphism in the equine LGB1 and LGB2 gene sequences (all exons, introns, and 5’-flanking region) and to assess potential relationship of particular genotypes with gene expression levels (measured in milk somatic cells) and milk composition traits (protein, fat, lactose, and total β-lactoglobulin content). Direct DNA sequencing analysis was performed for twelve horse breeds: Polish Primitive Horse (PPH), Polish Coldblood Horse (PCH), Polish Warmblood Horse (PWH), Silesian, Hucul, Fjording, Haflinger, Shetland Pony, Welsh Pony, Arabian, Thoroughbred, and Percheron—and revealed the presence of 83 polymorphic sites (47 and 36 for LGB1 and LGB2 genes, respectively), including eight that were previously unknown. Association analysis of the selected polymorphisms, gene expression, and milk composition traits (conducted for the PPH, PCH, and PWH breeds) showed several statistically significant relationships; for example, the two linked LGB1 SNPs (rs1143515669 and rs1144647991) were associated with total milk protein content (p < 0.01). Our study also confirmed that horse breed had significant impact on both gene transcript levels (p < 0.01) and on milk LGB content (p < 0.05), whereas an influence of lactation period was seen only for gene relative mRNA abundances (p < 0.01).

Whey protein polymorphisms in Sudanese goat breeds

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

Unravelling genetic variation underlying de novo-synthesis of bovine milk fatty acids

The relative abundance of specific fatty acids in milk can be important for consumer health and manufacturing properties of dairy products. Understanding of genes controlling milk fat synthesis may contribute to the development of dairy products with high quality and nutritional value. This study aims to identify key genes and genetic variants affecting de novo synthesis of the short- and medium-chained fatty acids C4:0 to C14:0. A genome-wide association study using 609,361 SNP markers and 1,811 animals was performed to detect genomic regions affecting fatty acid levels. These regions were further refined using sequencing data to impute millions of additional genetic variants. Results suggest associations of PAEP with the content of C4:0, AACS with the content of fatty acids C4:0-C6:0, NCOA6 or ACSS2 with the longer chain fatty acids C6:0-C14:0, and FASN mainly associated with content of C14:0. None of the top-ranking markers caused amino acid shifts but were mostly situated in putatively regulating regions and suggested a regulatory role of the QTLs. Sequencing mRNA from bovine milk confirmed the expression of all candidate genes which, combined with knowledge of their roles in fat biosynthesis, supports their potential role in de novo synthesis of bovine milk fatty acids.

Evaluation of the accuracy of imputed sequence variant genotypes and their utility for causal variant detection in cattle

Background

The availability of dense genotypes and whole-genome sequence variants from various sources offers the opportunity to compile large datasets consisting of tens of thousands of individuals with genotypes at millions of polymorphic sites that may enhance the power of genomic analyses. The imputation of missing genotypes ensures that all individuals have genotypes for a shared set of variants.

Results

We evaluated the accuracy of imputation from dense genotypes to whole-genome sequence variants in 249 Fleckvieh and 450 Holstein cattle using Minimac and FImpute. The sequence variants of a subset of the animals were reduced to the variants that were included on the Illumina BovineHD genotyping array and subsequently inferred in silico using either within- or multi-breed reference populations. The accuracy of imputation varied considerably across chromosomes and dropped at regions where the bovine genome contains segmental duplications. Depending on the imputation strategy, the correlation between imputed and true genotypes ranged from 0.898 to 0.952. The accuracy of imputation was higher with Minimac than FImpute particularly for variants with a low minor allele frequency. Using a multi-breed reference population increased the accuracy of imputation, particularly when FImpute was used to infer genotypes. When the sequence variants were imputed using Minimac, the true genotypes were more correlated to predicted allele dosages than best-guess genotypes. The computing costs to impute 23,256,743 sequence variants in 6958 animals were ten-fold higher with Minimac than FImpute. Association studies with imputed sequence variants revealed seven quantitative trait loci (QTL) for milk fat percentage. Two causal mutations in the DGAT1 and GHR genes were the most significantly associated variants at two QTL on chromosomes 14 and 20 when Minimac was used to infer genotypes.

Conclusions

The population-based imputation of millions of sequence variants in large cohorts is computationally feasible and provides accurate genotypes. However, the accuracy of imputation is low in regions where the genome contains large segmental duplications or the coverage with array-derived single nucleotide polymorphisms is poor. Using a reference population that includes individuals from many breeds increases the accuracy of imputation particularly at low-frequency variants. Considering allele dosages rather than best-guess genotypes as explanatory variables is advantageous to detect causal mutations in association studies with imputed sequence variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0301-x) contains supplementary material, which is available to authorized users.

Differential Expression of the Alpha S1 Casein and Beta-Lactoglobulin Genes in Different Physiological Stages of the Adani Goats Mammary Glands

Background

Milk proteins genes have been the focus of the researches as the candidate target genes that play a decisive role when animal breeding is desired.

Objectives

In the present study, the transcriptional levels of Beta-lactoglobulin (BLG) and Alpha S1 casein (CSN1S1) genes were investigated during prenatal, milking and drying times in mammary glands of the Adani goats which showed high and low breeding values.

Materials and Methods

The breeding values of the animals were estimated first by applying multi-trait random regression model. Using the biopsy gun, the mammary gland samples were taken and real-time PCR was applied to search the expression of the genes. Fixed factors of the model were the breeding value groups, sampling times and their interactions.

Results

The interactions were significant for both genes. At milking time, the high breeding value group exhibited more transcriptional levels for BLG and less transcriptional levels for CSN1S1 gene compared with the low breeding value group. The expression patterns of these genes were also different between the two breeding value groups. The maximum level of BLG and CSN1S1 transcriptions were found to occur at drying time.

Conclusions

A difference in the gene expression was observed between the two groups which indicate the change in the nucleotide sequence for transcription factor binding sites, or miRNA binding sites, otherwise in the coding regions. Therefore, the variations in the coding and promoter regions of this gene should be investigated in the further studies.

Exploiting biological priors and sequence variants enhances QTL discovery and genomic prediction of complex traits

BACKGROUND

Dense SNP genotypes are often combined with complex trait phenotypes to map causal variants, study genetic architecture and provide genomic predictions for individuals with genotypes but no phenotype. A single method of analysis that jointly fits all genotypes in a Bayesian mixture model (BayesR) has been shown to competitively address all 3 purposes simultaneously. However, BayesR and other similar methods ignore prior biological knowledge and assume all genotypes are equally likely to affect the trait. While this assumption is reasonable for SNP array genotypes, it is less sensible if genotypes are whole-genome sequence variants which should include causal variants.

RESULTS

We introduce a new method (BayesRC) based on BayesR that incorporates prior biological information in the analysis by defining classes of variants likely to be enriched for causal mutations. The information can be derived from a range of sources, including variant annotation, candidate gene lists and known causal variants. This information is then incorporated objectively in the analysis based on evidence of enrichment in the data. We demonstrate the increased power of BayesRC compared to BayesR using real dairy cattle genotypes with simulated phenotypes. The genotypes were imputed whole-genome sequence variants in coding regions combined with dense SNP markers. BayesRC increased the power to detect causal variants and increased the accuracy of genomic prediction. The relative improvement for genomic prediction was most apparent in validation populations that were not closely related to the reference population. We also applied BayesRC to real milk production phenotypes in dairy cattle using independent biological priors from gene expression analyses. Although current biological knowledge of which genes and variants affect milk production is still very incomplete, our results suggest that the new BayesRC method was equal to or more powerful than BayesR for detecting candidate causal variants and for genomic prediction of milk traits.

CONCLUSIONS

BayesRC provides a novel and flexible approach to simultaneously improving the accuracy of QTL discovery and genomic prediction by taking advantage of prior biological knowledge. Approaches such as BayesRC will become increasing useful as biological knowledge accumulates regarding functional regions of the genome for a range of traits and species.

Variability of the caprine whey protein genes and their association with milk yield, composition and renneting properties in the Sarda breed: 2. The BLG gene

The variability of the promoter region and the 3'UTR (exon-7) of the BLG gene, encoding the β-lactoglobulin, was investigated by sequencing in 263 lactating Sarda goats in order to assess its association with milk traits. Milk traits included: milk yield, fat, total protein and lactose content, pH, daily fat and protein yield (DFPY), freezing point, milk energy, somatic cell count, total microbial mesophilic count, rennet coagulation time (RCT), curd firming rate (k20) and curd firmness (a30). A total of 7 polymorphic sites were detected and the sequence analysed was given accession number KM817769. Only three SNPs (c.-381C>T, c.-323C>T and c.*420C>A) had minor allele frequency higher than 0.05. The effects of farm, stage of lactation and the interaction farm × stage of lactation significantly influenced all the milk traits (P T and c.*420C>A (P T (P < 0.001). The c.-381TT homozygous goats showed lower pH, RCT and k20 than c.-381CT (P < 0.05). In conclusion the polymorphism of the goat BLG gene did not affect the total protein content of the Sarda goat milk, and only weakly influenced RCT and k20. On the other hand, an interesting effect on milk yields and DFPY emerged in two SNPs. This information might be useful in dairy goat breeding programs.

Polymorphisms of beta-lactoglobulin promoter region in three Sicilian goat breeds

Several beta-lactoglobulin (BLG) polymorphisms have been described within the proximal promoter region and coding region of the caprine gene, although no genetic variants affecting the protein amino acid composition and/or expression level have been characterized so far. Binding sites for several transcription factors (TFs) are present in the BLG promoter region. The aims of this work were to sequence the full-length promoter region of three Sicilian goat breeds in order to identify polymorphisms, analyze the identified haplotypes, search for differences between breeds for the presence of polymorphisms in this gene region, search for putative TFs binding sites, and check if polymorphisms lay within the identified TFs binding sites. The promoter region of BLG gene in Sicilian goat breeds showed high level of polymorphism due to the presence of 36 single nucleotide polymorphisms (SNPs). Association between polymorphic sites was computed within the whole sample analyzed and 18 haplotypes were inferred. Binding sites for three milk protein binding factors (MPBFs) and four nuclear factor-I (NF-I) were found within BLG promoter region based on the ovine sequence. The identification of some SNPs within TFs binding sites allowed hypothesizing the loss of TFs. Further studies are in progress to evaluate the effect of these mutations on binding affinity of TFs, the functional interaction of the TFs with the goat BLG promoter, and the relationship of the polymorphisms with BLG gene expression and milk production and composition.

Study of polymorphisms in the promoter region of ovine β-lactoglobulin gene and phylogenetic analysis among the Valle del Belice breed and other sheep breeds considered as ancestors

The aim of this work was to sequence the promoter region of β-lactoglobulin (BLG) gene in four sheep breeds, in order to identify polymorphisms, infer and analyze haplotypes, and phylogenetic relationship among the Valle del Belice breed and the other three breeds considered as ancestors. Sequencing analysis and alignment of the obtained sequences showed the presence of 36 single nucleotide polymorphisms (SNPs) and one deletion. A total of 22 haplotypes found in "best" reconstruction were inferred considering the 37 polymorphic sites identified. Haplotypes were used for the reconstruction of a phylogenetic tree using the Neighbor-Joining algorithm. The number of polymorphisms identified showed high variability within breeds. Analysis of genetic diversity indexes showed that the Sarda breed presented the lowest nucleotide diversity, whereas the Comisana breed presented the highest one. Comparing the nucleotide diversity among breeds, the highest value was obtained between Valle del Belice and Pinzirita breeds, whereas the lowest one was between Valle del Belice and Sarda breeds. Considering that polymorphisms in the promoter region of BLG gene could have a functional role associated with milk composition, the lowest value of nucleotide diversity between Valle del Belice and Sarda breeds may be related to a higher similarity of milk composition of these two breeds compared to the others. Further analyses will be conducted in order to evaluate the possible correlation between the genetic diversity indexes and the BLG content in milk of our breeds.

Invited review: milk protein polymorphisms in cattle: effect on animal breeding and human nutrition

The 6 main milk proteins in cattle are encoded by highly polymorphic genes characterized by several nonsynonymous and synonymous mutations, with up to 47 protein variants identified. Such an extensive variation was used for linkage analysis with the description of the casein cluster more than 30 yr ago and has been applied to animal breeding for several years. Casein haplotype effects on productive traits have been investigated considering information on the whole casein complex. Moreover, mutations within the noncoding sequences have been shown to affect the specific protein expression and, as a consequence, milk composition and cheesemaking. Milk protein variants are also a useful tool for breed characterization, diversity, and phylogenetic studies. In addition, they are involved in various aspects of human nutrition. First, the occurrence of alleles associated with a reduced content of different caseins might be exploited for the production of milk with particular nutritional qualities; that is, hypoallergenic milk. On the other hand, the frequency of these alleles can be decreased by selection of sires using simple DNA tests, thereby increasing the casein content in milk used for cheesemaking. Furthermore, the biological activity of peptides released from milk protein digestion can be affected by amino acid exchanges or deletions resulting from gene mutations. Finally, the gene-culture coevolution between cattle milk protein genes and human lactase genes, which has been recently highlighted, is impressive proof of the nonrandom occurrence of milk protein genetic variation over the centuries.

Effect of exopolysaccharides on the hydrolysis of beta-lactoglobulin by Lactobacillus acidophilus CRL 636 in an in vitro gastric/pancreatic system

An analysis of the peptides generated by hydrolysis of BLG by nonproliferating cells of the strain Lactobacillus acidophilus CRL 636 was carried out. The effect of polysaccharides (pectin, and two EPS synthesized by two Streptococcus thermophilus strains, EPS1190 and EPS804) on BLG digestibility using an in vitro gastric/pancreatic system was analyzed. Polysaccharides are commonly used in the dairy industry to improve food texture; these hydrocolloids may interact with proteins, affecting their digestibility. Nonproliferating cells of Lb. acidophilus CRL 636 were able to hydrolyze 52% of BLG. Twenty-six resulting peptides with molecular masses in the range 544-4119 Da were identified by LC-MS/MS. These peptides resulted mostly from the hydrolysis of the more accessible N-terminal part of BLG. Degradation of BLG by pepsin was poor (8%). When BLG was previously hydrolyzed by Lb. acidophilus CRL 636, peptic hydrolysis was of 54.8%, while when pectin and EPS1190 were added, hydrolysis was higher (58.2 and 57.2%, respectively). Peptides crossing 8 kDa dialysis membranes after trypsin/chymotrypsin hydrolysis were analyzed by HPSEC. The produced peptides were smaller when BLG was hydrolyzed previously by the Lb. acidophilus strain. Moreover, in the presence of pectin, the amount of the larger peptide (3.5 kDa) observed in the size exclusion chromatograms was considerably decreased. Our studies showed that prehydrolysis of BLG by Lb. acidophilus CRL 636 had a positive influence on BLG digestibility and that polysaccharides may change the peptide profile yielded by trypsin/chymotrypsin hydrolysis, releasing smaller size peptides, which are known to be less immune-reactive. Moreover, Lb. acidophilus CRL 636 was able to hydrolyze the main epitopes (41-60, 102-124, and 149-162) of BLG, reducing its allergenic content.

Novel polymorphisms in the bovine beta-lactoglobulin gene and their effects on beta-lactoglobulin protein concentration in milk

The aim of our study was to detect new polymorphisms in the bovine beta-lactoglobulin (beta-LG) gene with significant effects on beta-LG protein concentration. Genomic DNA samples from 22 proven bulls were screened for polymorphisms in the coding and promoter regions of the beta-LG gene. In total, 50 polymorphisms were detected. Two single nucleotide polymorphisms (SNPs) (g.1772G>A and g.3054C>T) lead to amino acid changes and are the causal genetic polymorphisms of beta-LG protein variants A and B. Forty-two polymorphisms were in complete linkage disequilibrium (LD) with beta-LG protein variants A and B. Any of these 42 polymorphisms can be involved in the differential expression of the respective A and B alleles of the beta-LG gene. The eight polymorphisms not in complete LD with beta-LG protein variants A and B and the two polymorphisms causing the amino acid changes were genotyped in a set of 208 cows: 106 animals homozygous for beta-LG protein variant A and 102 animals homozygous for beta-LG protein variant B. Of these eight polymorphisms, six SNPs segregated only within the cows homozygous for beta-LG protein variant A and two SNPs segregated only within the cows homozygous for beta-LG protein variant B. One of the eight polymorphisms had a significant effect on beta-LG protein concentration. This SNP, g.-731G>A, segregated only within the 106 cows homozygous for beta-LG protein variant A. Within these cows, adjusted relative beta-LG protein concentration was reduced by 1.22% (w/w) in animals homozygous g.-731AA compared with animals homozygous g.-731GG.

Associations between 2 paternal casein haplotypes and milk yield traits of Swiss Fleckvieh cattle

Associations between casein haplotypes and milk yield traits of offspring from 5 Swiss Fleckvieh AI test bulls were investigated. The analysis was performed by using a daughter design, where each daughter inherited either paternal haplotype B-A1-A-A or B-A2-A-A for alleles of alpha s1-, beta-, alpha s2- and kappa-casein genes. The substitution effects of paternal CSN2 A1 versus A2 on protein yield deviations (YDs) were significant (P < 0.05), whereas their effects on milk and fat YDs were not. The paternal substitution effects of the CSN2 A1 versus the A2 allele on protein YDs within the 5 sires did not reach the significance level. This is due to the contrary allele substitution effect of a sire compared to the other 4 sires. The effects of maternal haplotypes on milk, protein and fat YDs were not significant. However, it is noteworthy that the effects of haplotypes with a low frequency in the population deviate largely from the most frequent haplotype B-A2-A-A. The effects of beta-lactoglobulin (BLG) genotypes were significant for protein YDs but not for milk and fat YDs. The association between the paternal CSN2 A1 and A2 alleles and milk protein YDs within sires but not milk and fat YDs indicate an interaction, which might be a consequence of CSN2 heterogeneity or a closely linked gene that is contributing to the estimated effects.