Evidence of additive genetic variation for major milk proteins in dairy cows: A meta-analysis

In the past, there have been reports of genetic parameters for milk proteins in various dairy cattle populations. The high variability among genetic parameter estimates has been caused by this. This study aimed to use a random-effects meta-analysis model to compile published estimates of genetic parameter for major milk proteins of α-lactalbumin, β-lactoglobulin, sum of whey proteins, casein, αs1-casein, αs2-casein, β-casein, and κ-casein in dairy cows. The study used a total of 140 heritability and 256 genetic correlation estimates from 23 papers published between 2004 and 2022. The estimated range of milk protein heritability is from 0.284 (for α-lactalbumin in milk) to 0.596 (for sum of whey proteins). The genetic correlation estimates between casein and milk yield, milk fat and protein percentages were -0.461, 0.693, and 0.976, respectively (p < 0.05). The genetic correlation estimates between milk proteins expressed as a percentage of milk were significant and varied from 0.177 (between β-lactoglobulin and κ-casein) to 0.892 (between αs1-casein and αs2-casein). Moderate-to-high heritability estimates for milk proteins and their low genetic associations with milk yield and composition indicated the possibility for improving milk proteins in a genetic selection plan with negligible correlated effects on production traits in dairy cows.

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.

Predicting milk protein fractions using infrared spectroscopy and a gradient boosting machine for breeding purposes in Holstein cattle

In recent years, increasing attention has been focused on the genetic evaluation of protein fractions in cow milk with the aim of improving milk quality and technological characteristics. In this context, advances in high-throughput phenotyping by Fourier transform infrared (FTIR) spectroscopy offer the opportunity for large-scale, efficient measurement of novel traits that can be exploited in breeding programs as indicator traits. We took milk samples from 2,558 Holstein cows belonging to 38 herds in northern Italy, operating under different production systems. Fourier transform infrared spectra were collected on the same day as milk sampling and stored for subsequent analysis. Two sets of data (i.e., phenotypes and FTIR spectra) collected in 2 different years (2013 and 2019-2020) were compiled. The following traits were assessed using HPLC: true protein, major casein fractions [α_S1-casein (CN), α_S2-CN, β-CN, κ-CN, and glycosylated-κ-CN], and major whey proteins (β-lactoglobulin and α-lactalbumin), all of which were measured both in grams per liter (g/L) and proportion of total nitrogen (% N). The FTIR predictions were calculated using the gradient boosting machine technique and tested by 3 different cross-validation (CRV) methods. We used the following CRV scenarios: (1) random 10-fold, which randomly split the whole into 10-folds of equal size (9-folds for training and 1-fold for validation); (2) herd/date-out CRV, which assigned 80% of herd/date as the training set with independence of 20% of herd/date assigned as the validation set; (3) forward/backward CRV, which split the data set in training and validation set according with the year of milk sampling (FTIR and gold standard data assessed in 2013 or 2019-2020) using the "old" and "new" databases for training and validation, and vice-versa with independence among them; (4) the CRV for genetic parameters (CRV-gen), where animals without pedigree as assigned as a fixed training population and animals with pedigree information was split in 5-folds, in which 1-fold was assigned to the fixed training population, and 4-folds were assigned to the validation set (independent from the training set). The results (i.e., measures and predictions) of CRV-gen were used to infer the genetic parameters for gold standard laboratory measurements (i.e., proteins assessed with HPLC) and FTIR-based predictions considering the CRV-gen scenario from a bi-trait animal model using single-step genomic BLUP. We found that the prediction accuracies of the gradient boosting machine equations differed according to the way in which the proteins were expressed, achieving higher accuracy when expressed in g/L than when expressed as % N in all CRV scenarios. Concerning the reproducibility of the equations over the different years, the results showed no relevant differences in predictive ability between using "old" data as the training set and "new" data as the validation set and vice-versa. Comparing the additive genetic variance estimates for milk protein fractions between the FTIR predicted and HPLC measures, we found reductions of -19.7% for milk protein fractions expressed in g/L, and -21.19% expressed as % N. Although we found reductions in the heritability estimates, they were small, with values ranging from -1.9 to -7.25% for g/L, and -1.6 to -7.9% for % N. The posterior distributions of the additive genetic correlations (r_a) between the FTIR predictions and the laboratory measurements were generally high (>0.8), even when the milk protein fractions were expressed as % N. Our results show the potential of using FTIR predictions in breeding programs as indicator traits for the selection of animals to enhance milk protein fraction contents. We expect acceptable responses to selection due to the high genetic correlations between HPLC measurements and FTIR predictions.

Phenotypic Analysis of Fourier-Transform Infrared Milk Spectra in Dairy Goats

The infrared spectrum of bovine milk is used to predict many interesting traits, whereas there have been few studies on goat milk in this regard. The objective of this study was to characterize the major sources of variation in the absorbance of the infrared spectrum in caprine milk samples. A total of 657 goats belonging to 6 breeds and reared on 20 farms under traditional and modern dairy systems were milk-sampled once. Fourier-transform infrared (FTIR) spectra were taken (2 replicates per sample, 1314 spectra), and each spectrum contained absorbance values at 1060 different wavenumbers (5000 to 930 × cm^-1), which were treated as a response variable and analyzed one at a time (i.e., 1060 runs). A mixed model, including the random effects of sample/goat, breed, flock, parity, stage of lactation, and the residual, was used. The pattern and variability of the FTIR spectrum of caprine milk was similar to those of bovine milk. The major sources of variation in the entire spectrum were as follows: sample/goat (33% of the total variance); flock (21%); breed (15%); lactation stage (11%); parity (9%); and the residual unexplained variation (10%). The entire spectrum was segmented into five relatively homogeneous regions. Two of them exhibited very large variations, especially the residual variation. These regions are known to be affected by the absorbance of water, although they also exhibited wide variations in the other sources of variation. The average repeatability of these two regions were 45% and 75%, whereas for the other three regions it was about 99%. The FTIR spectrum of caprine milk could probably be used to predict several traits and to authenticate the origin of goat milk.

Spectroscopic technologies and data fusion: Applications for the dairy industry

Increasing consumer awareness, scale of manufacture, and demand to ensure safety, quality and sustainability have accelerated the need for rapid, reliable, and accurate analytical techniques for food products. Spectroscopy, coupled with Artificial Intelligence-enabled sensors and chemometric techniques, has led to the fusion of data sources for dairy analytical applications. This article provides an overview of the current spectroscopic technologies used in the dairy industry, with an introduction to data fusion and the associated methodologies used in spectroscopy-based data fusion. The relevance of data fusion in the dairy industry is considered, focusing on its potential to improve predictions for processing traits by chemometric techniques, such as principal component analysis (PCA), partial least squares regression (PLS), and other machine learning algorithms.

Role of CSN2, CSN3, and BLG genes and the polygenic background in the cattle milk protein profile

To devise better selection strategies in dairy cattle breeding programs, a deeper knowledge of the role of the major genes encoding for milk protein fractions is required. The aim of the present study was to assess the effect of the CSN2, CSN3, and BLG genotypes on individual protein fractions (α_S1-CN, α_S2-CN, β-CN, κ-CN, β-LG, α-LA) expressed qualitatively as percentages of total nitrogen content (% N), quantitatively as contents in milk (g/L), and as daily production levels (g/d). Individual milk samples were collected from 1,264 Brown Swiss cows reared in 85 commercial herds in Trento Province (northeast Italy). A total of 989 cows were successfully genotyped using the Illumina Bovine SNP50 v.2 BeadChip (Illumina Inc.), and a genomic relationship matrix was constructed using the 37,519 SNP markers obtained. Milk protein fractions were quantified and the β-CN, κ-CN, and β-LG genetic variants were identified by reversed-phase HPLC (RP-HPLC). All protein fractions were analyzed through a Bayesian multitrait animal model implemented via Gibbs sampling. The effects of days in milk, parity order, and the CSN2, CSN3, and BLG genotypes were assigned flat priors in this model, whereas the effects of herd and animal additive genetic were assigned Gaussian prior distributions, and inverse Wishart distributions were assumed for the respective co-variance matrices. Marginal posterior distributions of the parameters of interest were compared before and after the inclusion of the effects of the 3 major genes in the model. The results showed that a high portion of the genetic variance was controlled by the major genes. This was particularly apparent in the qualitative protein profile, which was found to have a higher heritability than the protein fraction contents in milk and their daily yields. When the genes were included individually in the model, CSN2 was the major gene controlling all the casein fractions except for κ-CN, which was controlled directly by the CSN3 gene. The BLG gene had the most influence on the 2 whey proteins. The genetic correlations showed the major genes had only a small effect on the relationships between the protein fractions, but through comparison of the correlation coefficients of the proteins expressed in different ways they revealed potential mechanisms of regulation and competitive synthesis in the mammary gland. The estimates for the effects of the CSN2 and CSN3 genes on protein profiles showed overexpression of protein synthesis in the presence of the B allele in the genotype. Conversely, the β-LG B variant was associated with a lower concentration of β-LG compared with the β-LG A variant, independently of how the protein fractions were expressed, and it was followed by downregulation (or upregulation in the case of the β-LG B) of all other protein fractions. These results should be borne in mind when seeking to design more efficient selection programs aimed at improving milk quality for the efficiency of dairy industry and the effect of dairy products on human health.

Breeding for improved protein fractions and free amino acids concentration in bovine milk

Considerable resources are required to routinely measure detailed milk compositional traits. Hence, an insufficient volume of phenotypic data can hinder genetic progress in these traits within dairy cow breeding programmes. The objective of the present study was to quantify the opportunities for breeding for improved milk protein and free amino acid (FAA) composition by exploiting mid‐infrared spectroscopy (MIRS) predictions routinely recorded from milk samples. Genetic parameters for protein fractions and FAA composition were estimated using 134,546 test‐day records from 16,166 lactations on 9,572 cows using linear mixed models. Heritability of MIRS‐predicted protein fractions ranged from 0.19 (α‐lactalbumin) to 0.55 (β‐lactoglobulin A), while heritability of MIRS‐predicted FAA ranged from 0.08 for glycine to 0.29 for glutamic acid. Genetic correlations among the MIRS‐predicted FAA were moderate to strong ranging from −0.44 (aspartic acid and lysine) to 0.97 (glutamic acid and total FAA). Adjustment of the genetic correlations for the genetic merit of 24‐h milk yield did not greatly affect the correlations. Results from the current study highlight the presence of exploitable genetic variation for both protein fractions and FAA in dairy cow milk. Besides, the direction of genetic correlations reveals that breeding programmes directly selecting for greater milk protein concentration carry with them favourable improvement in casein and whey fractions.

Polymorphisms of SORBS1 Gene and Their Correlation with Milk Fat Traits of Cattleyak

Simple Summary

Increasing milk fat rate has a good effect on the milk quality of cattleyak. SNPs can help us find potential molecular markers for the milk fat traits of cattleyak, and they can be screened according to molecular markers when they are young. It provides a reference for cultivating high milk fat cattle population in the future. The results of this study suggest that the SORBS1 gene polymorphism is closely related to the milk fat traits of cattleyak, which could be used as a candidate genetic marker for milk fat trait selection in cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak. It has a certain reference value for the research and improvement of milk quality.

Abstract

This study aimed to find the SNPs in the SORBS1 gene of cattleyak, analyze the relationship between its polymorphisms and the milk fat traits, and find potential molecular markers for the milk fat traits of cattleyak. The polymorphism of the SORBS1 gene in 350 cattleyak from Hongyuan County (Sichuan, China) were detected by PCR and DNA sequencing, and the correlation between these SNPs and the milk production traits of cattleyak was analyzed. The results showed that there were nine SNPs in the CDS and their adjacent non-coding regions of the SORBS1 gene, and all SNPs have three genotypes. The correlation analysis found that the genotypes with superior milk fat traits in the other eight alleles were homozygous genotypes with a high genotype frequency except the g.96284 G > A (c.3090 G > A) (p < 0.05). However, at locus g.96284 G > A, the milk fat percentage, monounsaturated fatty acids (MUFAs), polyunsaturated fatty acids (PUFAs) and saturated fatty acids (SFAs) of the GA genotype were significantly higher than that of GG and AA genotypes (p < 0.05). Among these SNPs, three SNPs (g.6256 C > T (c.298 C > T), g.24791 A > G (c.706 A > G) and g.29121 A > G (c.979 A > G)) caused the amino acids change. The genotypes of the three SNPs consist of three haplotypes and four diplotypes. The amino acid mutation degree of diplotype H1–H1 (CCAAAA) was the highest, and its milk fat percentage, MUFAs, PUFAs and SFAs were also the highest (p < 0.05). Taken together, we found nine SNPs in the SORBS1 gene that are closely related to the milk fat traits of cattleyak. Moreover, the mutation of amino acids caused by SNPs had positive effects on the milk fat traits of cattleyak. H1-H1 is the dominant diplotype which significantly related to the milk fat traits of cattleyak. This study provides a new molecular marker and theoretical basis for screening the milk fat traits of cattleyak.

Effects of breed, farm intensiveness, and cow productivity level on cheese-making ability predicted using infrared spectral data at the population level

Fourier-transform infrared (FTIR) spectra collected during milk recording schemes at population level can be used for predicting novel traits of interest for farm management, cows' genetic improvement, and milk payment systems. The aims of this study were as follows. (1) To predict cheese yield traits using FTIR spectra from routine milk recordings exploiting previously developed calibration equations. (2) To compare the predicted cheese-making abilities of different dairy and dual-purpose breeds. (3) To analyze the effects of herds' level of intensiveness (HL) and of the cow's level of productivity (CL). (4) To compare the patterns of predicted cheese yields with the patterns of milk composition in different breeds to discern the drivers of cheese-making efficiency. The major sources of variation of FTIR predictions of cheese yield ability (fresh cheese or cheese solids produced per unit milk) of individual milk samples were studied on 115,819 cows of 4 breeds (2 specialized dairy breeds, Holstein and Brown Swiss, and 2 dual-purpose breeds, Simmental and Alpine Grey) from 6,430 herds and exploiting 1,759,706 FTIR test-day spectra collected over 7 yr of milk sampling. Calibration equations used were previously developed on 1,264 individual laboratory model cheese procedures (cross-validation R² 0.85 and 0.95 for fresh and solids cheese yields, respectively). The linear model used for statistical analysis included the effects of parity, lactation stage, year of calving, month of sampling, HL, CL, breed of cow, and the interactions breed × HL and breed × CL. The HL and CL stratifications (5 classes each) were based on average daily secretion of milk net energy per cow. All effects were highly significant (P < 0.001). The major conclusions were as follows. (1) The FTIR-based prediction of cheese yield of milk goes beyond the knowledge of fat and protein content, partially explaining differences in cheese-making ability in different cows, breeds and herds. (2) Differences in cheese yields of different breeds are only partially explained by milk fat and protein composition, and less productive breeds are characterized by a higher milk nutrient content as well as a higher recovery of nutrients in the cheese. (3) High-intensive herds not only produce much more milk, but the milk has a higher nutrient content and a higher cheese yield, whereas within herds, compared with less productive cows, the more productive cows have a much greater milk yield, milk with a greater content of fat but not of protein, and a moderate improvement in cheese yield, differing little from expectations based on milk composition. Finally, (4) the effects of HL and CL on milk quality and cheese-making ability are similar but not identical in different breeds, the less productive ones having some advantage in terms of cheese-making ability. We can obtain FTIR-based prediction of cheese yield from individual milk samples retrospectively at population level, which seems to go beyond the simple knowledge of milk composition, incorporating information on nutrient retention ability in cheese, with possible advantages for management of farms, genetic improvement of dairy cows, and milk payment systems.

Sequence-based genome-wide association study of individual milk mid-infrared wavenumbers in mixed-breed dairy cattle

BACKGROUND

Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk.

RESULTS

Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk.

CONCLUSIONS

This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.

Structural equation modeling for unraveling the multivariate genomic architecture of milk proteins in dairy cattle

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, α_S1-, and α_S2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and α_S1-CN (-0.706) and between β-CN and κ-CN (-0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for α_S1-CN, and 14 SNP for α_S2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.

Investigation of Raman Spectroscopy (with Fiber Optic Probe) and Chemometric Data Analysis for the Determination of Mineral Content in Aqueous Infant Formula

This study investigated the use of Raman spectroscopy (RS) and chemometrics for the determination of eight mineral elements (i.e., Ca, Mg, K, Na, Cu, Mn, Fe, and Zn) in aqueous infant formula (INF). The samples were prepared using infant formula powder reconstituted to concentrations of 3%-13% w/w (powder: water) (n = 83). Raman spectral data acquisition was carried out using a non-contact fiber optic probe on the surface of aqueous samples in 50-3398 cm-1. ICP-AES was used as a reference method for the determination of the mineral contents in aqueous INF samples. Results showed that the best performing partial least squares regression (PLSR) models developed for the prediction of minerals using all samples for calibration achieved R2CV values of 0.51-0.95 with RMSECVs of 0.13-2.96 ppm. The PLSR models developed and validated using separate calibration (n = 42) and validation (n = 41) samples achieved R2CVs of 0.93, 0.94, 0.91, 0.90, 0.97, and 0.94, R2Ps of 0.75, 0.77, 0.31, 0.60, 0.84, and 0.80 with RMSEPs of 3.17, 0.29, 3.45, 1.51, 0.30, and 0.25 ppm for the prediction of Ca, Mg, K, Na, Fe, and Zn respectively. This study demonstrated that RS equipped with a non-contact fiber optic probe and combined with chemometrics has the potential for timely quantification of the mineral content of aqueous INF during manufacturing.

Genetic parameters for noncoagulating milk, milk coagulation properties, and detailed milk composition in Swedish Red Dairy Cattle

The rennet-induced coagulation ability of milk is important in cheese production. For Swedish Red Dairy Cattle (RDC), this ability is reduced because of a high prevalence of noncoagulating (NC) milk. In this study, we simultaneously combined genetic parameters for NC milk, milk coagulation properties, milk composition, physical traits, and milk protein composition. Our aim was to estimate heritability and genetic and phenotypic correlations for NC milk and 24 traits (milk coagulation properties, milk composition, physical traits, and milk protein composition). Phenotypes and ∼7,000 SNP genotypes were available for all 600 Swedish RDC. The genotypes were imputed from ∼7,000 SNP to 50,000 SNP. Variance components and genetic parameters were estimated with an animal model. In Swedish RDC, a moderate heritability estimate of 0.28 was found for NC milk. For the other 24 traits, heritability estimates ranged from 0.12 to 0.77 (standard errors from 0.08 to 0.18). A total of 300 phenotypic and genetic correlations were estimated. For phenotypic and genetic correlations, 172 and 95 were significant, respectively. In general, most traits showing significant genetic correlations also showed significant phenotypic correlations. In this study, phenotypic and genetic correlations with NC milk suggest that many correlations between traits exist, making it difficult to predict the real consequences on the composition of milk, if selective breeding is applied on NC milk. We speculate that some of these consequences may lead to changes in the composition of milk, most likely affecting its physical and organoleptic properties. However, our results suggest that κ-casein could be used as an indicator trait to predict the occurrence of NC milk at the herd level.

The evolving role of Fourier-transform mid-infrared spectroscopy in genetic improvement of dairy cattle

Over the last 100 years, significant advances have been made in the characterisation of milk composition for dairy cattle improvement programs. Technological progress has enabled a shift from labour intensive, on-farm collection and processing of samples that assess yield and fat levels in milk, to large-scale processing of samples through centralised laboratories, with the scope extended to include quantification of other traits. Fourier-transform mid-infrared (FT-MIR) spectroscopy has had a significant role in the transformation of milk composition phenotyping, with spectral-based predictions of major milk components already being widely used in milk payment and animal evaluation systems globally. Increasingly, there is interest in analysing the individual FT-MIR wavenumbers, and in utilising the FT-MIR data to predict other novel traits of importance to breeding programs. This includes traits related to the nutritional value of milk, the processability of milk into products such as cheese, and traits relevant to animal health and the environment. The ability to successfully incorporate these traits into breeding programs is dependent on the heritability of the FT-MIR predicted traits, and the genetic correlations between the FT-MIR predicted and actual trait values. Linking FT-MIR predicted traits to the underlying mutations responsible for their variation can be difficult because the phenotypic expression of these traits are a function of a diverse range of molecular and biological mechanisms that can obscure their genetic basis. The individual FT-MIR wavenumbers give insights into the chemical composition of milk and provide an additional layer of granularity that may assist with establishing causal links between the genome and observed phenotypes. Additionally, there are other molecular phenotypes such as those related to the metabolome, chromatin accessibility, and RNA editing that could improve our understanding of the underlying biological systems controlling traits of interest. Here we review topics of importance to phenotyping and genetic applications of FT-MIR spectra datasets, and discuss opportunities for consolidating FT-MIR datasets with other genomic and molecular data sources to improve future dairy cattle breeding programs.

SERPINA1 gene identified in RNA-Seq showed strong association with milk protein concentration in Chinese Holstein cows

The detection of candidate genes and mutations associated with phenotypic traits is important for livestock animals. A previous RNA-Seq study revealed that SERPINA1 gene was a functional candidate that may affect milk protein concentration in dairy cows. To further confirm the genetic effect of SERPINA1 on milk protein traits, genetic polymorphisms were identified and genotype-phenotype associations were performed in a large Chinese Holstein cattle population. The entire coding region and the 5′-regulatory region (5′-UTR) of SERPINA1 was sequenced using pooled DNA of 17 unrelated sires. Association studies for five milk production traits were performed using a mixed model with a population encompassing 1,027 Chinese Holstein cows. A total of four SNPs were identified in SERPINA1, among which rs210222822 and rs41257068 presented in exons, rs207601878 presented in an intron, and rs208607693 was in the 5′-UTR. Analyses of pairwise D′ measures of linkage disequilibrium (LD) showed strong linkage among these four SNPs (D′ = 0.99–1.00), and a 9 Kb haplotype block involving three main haplotypes with GTGT, CCCC and CCGT was inferred. An association study revealed that all four single SNPs and their haplotypes had significant genetic effects on milk protein percentage, milk protein yield and milk yield (P = 0.0458 −  < 0.0001). The phenotypic variance ratio for all 11 significant SNP-trait pairs ranged from 1.01% to 7.54%. The candidate gene of SERPINA1 revealed by our previous RNA-Seq study was confirmed to have pronounced effect on milk protein traits on a genome level. Two SNPs (rs208607693 and rs210222822) presented phenotypic variances of approximately 7% and may be used as key or potential markers to assist selection for new lines of cows with high protein concentration.

Association of UDP-galactose-4-epimerase with milk protein concentration in the Chinese Holstein population

Objective

An initial RNA-Sequencing study revealed that UDP-galactose-4-epimerase (GALE) was one of the most promising candidates for milk protein concentration in Chinese Holstein cattle. This enzyme catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. To further validate the genetic effect of GALE on milk protein traits, genetic variations were identified, and genotypes-phenotypes associations were performed.

Methods

The entire coding region and the 5′-regulatory region (5′-UTR) of GALE were re-sequenced using pooled DNA of 17 unrelated sires. Association studies for five milk production traits were performed using a mixed linear animal model with a population encompassing 1,027 Chinese Holstein cows.

Results

A total of three variants in GALE were identified, including two novel variants (g.2114 A>G and g.2037 G>A) in the 5′-UTR and one previously reported variant (g.3836 G>C) in an intron. All three single nucleotide polymorphisms (SNPs) were associated with milk yield (p<0.0001), fat yield (p = 0.0006 to <0.0001), protein yield (p = 0.0232 to <0.0001) and protein percentage (p<0.0001), while no significant associations were detected between the SNPs and fat percentage. A strong linkage disequilibrium (D’ = 0.96 to 1.00) was observed among all three SNPs, and a 5 Kb haplotype block involving three main haplotypes with GAG, AGC, and AGG was formed. The results of haplotype association analyses were consistent with the results of single locus association analysis (p<0.0001). The phenotypic variance ratio above 3.00% was observed for milk protein yield that was explained by SNP-g.3836G >C.

Conclusion

Overall, our findings provided new insights into the polymorphic variations in bovine GALE gene and their associations with milk protein concentration. The data indicate their potential uses for marker-assisted breeding or genetic selection schemes.

Genomic Analysis of Milk Protein Fractions in Brown Swiss Cattle

Simple Summary

Milk protein fractions are hugely important in the dairy industry because of the key role they play in milk technological properties. The selection of cows for milk protein fractions may, therefore, improve both the nutritional and technological characteristics of milk, and, consequently, the processing efficiency and value of the dairy product. This study estimated the genetic parameters of the major milk protein fractions (four caseins, and two whey proteins) determined variously as: (i) milk content (g/100g milk), (ii) percentage of milk nitrogen (%N) and (iii) daily yield (g/d) in Brown Swiss dairy cattle. The results showed that the (co)variances and genetic parameter estimates differed according to how the proteins were measured. These results provide useful information for developing selection strategies in dairy cattle breeding programs aimed at improving both the nutritional and technological properties of milk.

Abstract

Depending on whether milk protein fractions are evaluated qualitatively or quantitatively, different genetic outcomes may emerge. In this study, we compared the genetic parameters for the major milk protein fractions—caseins (α_S1-, α_S2-, β-, and к-CN), and whey proteins (β-lactoglobulin, β-LG; α-lactalbumin, α-LA)—estimated using the multi-trait genomic best linear unbiased prediction method and expressed variously as milk content (g/100g milk), percentage of milk nitrogen (%N) and daily yield per cow (g/d). The results showed that the genetic parameter estimates varied according to how the milk protein fractions were expressed. Heritability estimates for the caseins and whey protein fractions expressed as daily yields were lower than when they were expressed as proportions and contents, revealing important differences in genetic outcomes. The proportion and the content of β-CN were negatively correlated with the proportions and contents of α_S1-CN, α_S2-CN, and к-CN, while the daily yield of β–CN was negatively correlated with the daily yields of α_S1-CN and α_S2-CN. The Spearman’s rank correlations and the coincidence rates between the various predicted genomic breeding values (GEBV) for the milk protein fractions expressed in different ways indicated that these differences had a significant effect on the ranking of the animals. The results suggest that the way milk protein fractions are expressed has implications for breeding programs aimed at improving milk nutritional and technological characteristics.

A Novel SNPs in Alpha-Lactalbumin Gene Effects on Lactation Traits in Chinese Holstein Dairy Cows

Simple Summary

Alpha-lactalbumin (α-LA) is a major whey protein component in mammalian milk, such as human (approximately 36%), bovine (approximately 17%), and other species, etc., It is involved in the regulation of lactose synthesis and has high nutritional value, especially in infant formula. Previous studies have confirmed that bovine α-LA gene 5′-flanking region has single nucleotide polymorphisms (SNPs), but little is known about polymorphisms in other regions, especially sequence coding for amino acids in protein (CDS) and their adjacent non-coding regions, including Chinese Holstein dairy cows. This study focused on investigated SNPs in the CDS and their adjacent non-coding regions of the α-LA gene in Chinese Holstein dairy cows, and assessed the association between SNPs and lactation traits. Sequence alignment showed that a potential SNPs (562th, G/A) in CDS2 region affect protein spatial structure, suggesting that this SNPs might affect the lactation traits of cows (milk type (Holstein and Jersey), and non-milk type (Bos Taurus)) need more in-depth study. More importantly, a novel SNPs at 1847th (T/C) bp in non-coding region near CDS4 was significantly associated with milk lactose composition, and lactose contents were significantly correlated with milk protein content, indicating that the SNPs could be used as a novel potential molecular marker for lactation traits in Chinese Holstein dairy cows.

Abstract

Alpha-lactalbumin (α-LA) is a major whey protein in bovine and other mammalian milk, which regulates synthesis of lactose. Little is known about its genetic polymorphism and whether can be used as a potential marker for dairy ingredients, milk yield traits, and milk properties. To investigate its polymorphisms and their relationship with milk lactation traits in Chinese Holstein dairy cows, single-strand conformation polymorphism method (PCR-SSCP) and direct sequencing method were used to mark the α-LA gene SNPs. AA (0.7402) and AB (0.2598) genotypes were screened out by PCR-SSCP bands analysis in two independent populations. Direct sequencing revealed that there is one SNP at 1847th (T/C) bp in noncoding region of α-LA gene with highly polymorphic (0.5 < PIC = 0.5623 or 0.5822), of which T is in AA genotype while C in AB. Association analysis also showed that lactose content (p < 0.05) was negatively correlated with fat and protein contents within subgroup, indicating that the SNPs (1847th, T/C) in α-LA gene could be used as a novel potential molecular marker for lactation traits in Chinese Holstein dairy cows.

Sequence-based GWAS, network and pathway analyses reveal genes co-associated with milk cheese-making properties and milk composition in Montbéliarde cows

BACKGROUND

Milk quality in dairy cattle is routinely assessed via analysis of mid-infrared (MIR) spectra; this approach can also be used to predict the milk's cheese-making properties (CMP) and composition. When this method of high-throughput phenotyping is combined with efficient imputations of whole-genome sequence data from cows' genotyping data, it provides a unique and powerful framework with which to carry out genomic analyses. The goal of this study was to use this approach to identify genes and gene networks associated with milk CMP and composition in the Montbéliarde breed.

RESULTS

Milk cheese yields, coagulation traits, milk pH and contents of proteins, fatty acids, minerals, citrate, and lactose were predicted from MIR spectra. Thirty-six phenotypes from primiparous Montbéliarde cows (1,442,371 test-day records from 189,817 cows) were adjusted for non-genetic effects and averaged per cow. 50 K genotypes, which were available for a subset of 19,586 cows, were imputed at the sequence level using Run6 of the 1000 Bull Genomes Project (comprising 2333 animals). The individual effects of 8.5 million variants were evaluated in a genome-wide association study (GWAS) which led to the detection of 59 QTL regions, most of which had highly significant effects on CMP and milk composition. The results of the GWAS were further subjected to an association weight matrix and the partial correlation and information theory approach and we identified a set of 736 co-associated genes. Among these, the well-known caseins, PAEP and DGAT1, together with dozens of other genes such as SLC37A1, ALPL, MGST1, SEL1L3, GPT, BRI3BP, SCD, GPAT4, FASN, and ANKH, explained from 12 to 30% of the phenotypic variance of CMP traits. We were further able to identify metabolic pathways (e.g., phosphate and phospholipid metabolism and inorganic anion transport) and key regulator genes, such as PPARA, ASXL3, and bta-mir-200c that are functionally linked to milk composition.

CONCLUSIONS

By using an approach that integrated GWAS with network and pathway analyses at the whole-genome sequence level, we propose candidate variants that explain a substantial proportion of the phenotypic variance of CMP traits and could thus be included in genomic evaluation models to improve milk CMP in Montbéliarde cows.

Variation of Detailed Protein Composition of Cow Milk Predicted from a Large Database of Mid-Infrared Spectra

Simple Summary

Milk proteins are one of the most valuable milk components. The objective of the present study was to assess sources of variation of detailed protein composition predicted from infrared spectra in milk of dairy and dual-purpose cattle breeds. Results showed that protein fractions were primarily influenced by days in milk, and the relative proportion of each fraction through lactation was not constant. Protein fractions correlated with crude protein, total casein, fat and milk urea nitrogen. In perspective, mid-infrared predictions of milk fractions could be useful for the dairy sector to improve nutritional and technological properties of milk.

Abstract

This study aimed to investigate factors affecting protein fractions, namely α-casein (α-CN), β-casein (β-CN), κ-casein (κ-CN), β-lactoglobulin (β-LG) and α-lactalbumin (α-LA) predicted from milk infrared spectra in milk of dairy and dual-purpose cattle breeds. The dataset comprised 735,328 observations from 49,049 cows in 1782 herds. Results highlighted significant differences of protein fractions in milk of the studied breeds. Significant variations of protein fractions were found also through parities and lactation, with the latter thoroughly influencing protein fractions percentage. Interesting correlations (r) were estimated between β-CN, κ-CN and β-LG, expressed as percentage of crude protein, and milk urea nitrogen (r = 0.31, −0.20 and −0.26, respectively) and between α-LA and fat percentage (r = 0.41). The present study paves the way for future studies on the associations between protein fractions and milk technological properties, and for the estimation of genetic parameters of predicted protein composition.

Within-breed and multi-breed GWAS on imputed whole-genome sequence variants reveal candidate mutations affecting milk protein composition in dairy cattle

Background

Genome-wide association studies (GWAS) were performed at the sequence level to identify candidate mutations that affect the expression of six major milk proteins in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) dairy cattle. Whey protein (α-lactalbumin and β-lactoglobulin) and casein (αs1, αs2, β, and κ) contents were estimated by mid-infrared (MIR) spectrometry, with medium to high accuracy (0.59 ≤ R² ≤ 0.92), for 848,068 test-day milk samples from 156,660 cows in the first three lactations. Milk composition was evaluated as average test-day measurements adjusted for environmental effects. Next, we genotyped a subset of 8080 cows (2967 MON, 2737 NOR, and 2306 HOL) with the BovineSNP50 Beadchip. For each breed, genotypes were first imputed to high-density (HD) using HD single nucleotide polymorphisms (SNPs) genotypes of 522 MON, 546 NOR, and 776 HOL bulls. The resulting HD SNP genotypes were subsequently imputed to the sequence level using 27 million high-quality sequence variants selected from Run4 of the 1000 Bull Genomes consortium (1147 bulls). Within-breed, multi-breed, and conditional GWAS were performed.

Results

Thirty-four distinct genomic regions were identified. Three regions on chromosomes 6, 11, and 20 had very significant effects on milk composition and were shared across the three breeds. Other significant effects, which partially overlapped across breeds, were found on almost all the autosomes. Multi-breed analyses provided a larger number of significant genomic regions with smaller confidence intervals than within-breed analyses. Combinations of within-breed, multi-breed, and conditional analyses led to the identification of putative causative variants in several candidate genes that presented significant protein–protein interactions enrichment, including those with previously described effects on milk composition (SLC37A1, MGST1, ABCG2, CSN1S1, CSN2, CSN1S2, CSN3, PAEP, DGAT1, AGPAT6) and those with effects reported for the first time here (ALPL, ANKH, PICALM).

Conclusions

GWAS applied to fine-scale phenotypes, multiple breeds, and whole-genome sequences seems to be effective to identify candidate gene variants. However, although we identified functional links between some candidate genes and milk phenotypes, the causality between candidate variants and milk protein composition remains to be demonstrated. Nevertheless, the identification of potential causative mutations that underlie milk protein composition may have immediate applications for improvements in cheese-making.

Electronic supplementary material

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