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.

Composition, coagulation properties, and predicted cheesemaking traits of bulk goat milk from different farming systems, breeds, and stages of production

At the global level, the quantity of goat milk produced and its gross production value have increased considerably over the last 2 decades. Although many scientific papers on this topic have been published, few studies have been carried out on bulk goat milk samples. The aim of the present study was to investigate in the field the effects of farming system, breed type, individual flock, and stage of production on the composition, coagulation properties (MCP), curd firming over time parameters (CF_t), predicted cheese yield (CY%), and nutrient recovery traits (REC) of 432 bulk milk samples from 161 commercial goat farms in Sardinia, Italy. We found that the variance due to individual flock was of the same order as the residual variance for almost all composition and cheesemaking traits. With regard to the fixed effects, the effect of farming system on bulk milk variability was not highly significant for the majority of traits (it was lower than individual flock), whereas the effects of breed type and stage of production were much higher. More specifically, the intensive farms produced milk with the best concentrations of almost all constituents, whereas extensive farms exhibited faster rennet coagulation times, a slower rate of curd firming, lower potential curd firmness, and lower percentages of fat and energy recoveries in the fresh curd. Farms rearing the local breed, Sarda, alone or together with the Maltese breed, produced milk with the best concentrations of fat and protein, superior curd firmness, and better predicted percentage of fresh curd (CY_CURD) and recovery traits. The results show the potential of both types of breed, either for their quantitative (specialized breeds) or their qualitative (local breeds) attributes. As expected, the concentrations of fat, protein fractions, and lactose were influenced by the stage of production, with samples collected in the early stage of production (in February and March) having a greater quantity of the main constituents. Somatic cells reached the highest levels in the late stage of production, which corresponds to the goats' advanced stage of lactation (June-July), although no differences were present in the logarithmic bacterial counts between the early and late stages. Regarding cheesemaking potential, bulk milk samples of the late stage were characterized by delayed rennet coagulation and curd firming times, the lowest values of curd firmness, and a general reduction in CY%, and REC traits. In conclusion, we highlight several issues regarding the effects of the most important sources of variation on bulk goat milk, and point to some critical factors relevant for improving dairy goat farming and milk production.

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.

Composition and aptitude for cheese-making of milk from cows, buffaloes, goats, sheep, dromedary camels, and donkeys

Bovines account for about 83% of the milk and dairy products consumed by humans worldwide, the rest represented by bubaline, caprine, ovine, camelid, and equine species, which are particularly important in areas of extensive pastoralism. Although milk is increasingly used for cheese production, the cheese-making efficiency of milk from the different species is not well known. This study compares the cheese-making ability of milk sampled from lactating females of the 6 dairy species in terms of milk composition, coagulation properties (using lactodynamography), curd-firming modeling, nutrients recovered in the curd, and cheese yield (through laboratory model-cheese production). Equine (donkey) milk had the lowest fat and protein content and did not coagulate after rennet addition. Buffalo and ewe milk yielded more fresh cheese (25.5 and 22.9%, respectively) than cow, goat, and dromedary milk (15.4, 11.9, and 13.8%, respectively). This was due to the greater fat and protein contents of the former species with respect to the latter, but also to the greater recovery of fat in the curd of bubaline (88.2%) than in the curd of camelid milk (55.0%) and consequent differences in the recoveries of milk total solids and energy in the curd; protein recovery, however, was much more similar across species (from 74.7% in dromedaries to 83.7% in bovine milk). Compared with bovine milk, the milk from the other Artiodactyla species coagulated more rapidly, reached curd firmness more quickly (especially ovine milk), had a more pronounced syneresis (especially caprine milk), had a greater potential asymptotical curd firmness (except dromedary and goat milk), and reached earlier maximum curd firmness (especially caprine and ovine milk). The maximum measured curd firmness was greater for bubaline and ovine milk, intermediate for bovine and caprine milk, and lower for camelid milk. The milk of all ruminant species can be used to make cheese, but, to improve efficiency, cheese-making procedures need to be optimized to take into account the large differences in their coagulation, curd-firming, and syneresis properties.

Associations between differential somatic cell count and milk yield, quality, and technological characteristics in Holstein cows

The aim of this study was to investigate the associations between differential somatic cell count (DSCC) and milk quality and udder health traits, and for the first time, between DSCC and milk coagulation properties and cheesemaking traits in a population of 1,264 Holstein cows reared in northern Italy. Differential somatic cell count represents the combined proportions of polymorphonuclear neutrophils plus lymphocytes (PMN-LYM) in the total somatic cell count (SCC), with macrophages (MAC) making up the remaining proportion. The milk traits investigated in this study were milk yield (MY), 8 traits related to milk composition and quality (fat, protein, casein, casein index, lactose, urea, pH, and milk conductivity), 9 milk coagulation traits [3 milk coagulation properties (MCP) and 6 curd firming (CF) traits], 7 cheesemaking traits, 3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits. A linear mixed model was fitted to explore the associations between SCS combined with DSCC and the aforementioned milk traits. An additional model was run, which included DSCC expressed as the PMN-LYM and MAC counts, obtained by multiplying the percentage of PMN-LYM and MAC by SCC in the milk for each cow in the data set. The unfavorable association between SCS and milk quality and technological traits was confirmed. Increased DSCC was instead associated with a linear increase in MY, casein index, and lactose proportion and a linear decrease in milk fat and milk conductivity. Accordingly, DSCC was favorably associated with all MCP and CF traits (with the exception of the time needed to achieve maximum, CF), particularly with rennet coagulation time, and it always displayed linear relationships. Differential somatic cell count was also positively associated with the recovery of milk nutrients in the curd (protein, fat, and energy), which increased linearly with increasing DSCC. The PMN-LYM count was rarely associated with milk traits, even though the pattern observed confirmed the results obtained when both SCS and DSCC were included in the model. The MAC count, however, showed the opposite pattern: MY, casein index, and lactose percentage decreased and milk conductivity increased with an increasing MAC count. No significant association was found between PMN-LYM count and MCP, CF, CY, and REC traits, whereas MAC count was unfavorably associated with MCP, CF traits, some CY traits, and all REC traits. Our results showed that the combined information derived from SCS and DSCC might be useful to monitor milk quality and cheesemaking-related traits.

Integration of Wet-Lab Measures, Milk Infrared Spectra, and Genomics to Improve Difficult-to-Measure Traits in Dairy Cattle Populations

The objective of this study was to evaluate the contribution of Fourier-transformed infrared spectroscopy (FTIR) data for dairy cattle breeding through two different approaches: (i) estimating the genetic parameters for 30 measured milk traits and their FTIR predictions and investigating the additive genetic correlation between them and (ii) evaluating the effectiveness of FTIR-derived phenotyping to replicate a candidate bull’s progeny testing or breeding value prediction at birth. Records were available from 1,123 cows phenotyped using gold standard laboratory methodologies (LAB data). This included phenotypes related to fine milk composition and milk technological characteristics, milk acidity, and milk protein fractions. The dataset used to generate FTIR predictions comprised 729,202 test-day records from 51,059 Brown Swiss cows (FIELD data). A first approach consisted of estimating genetic parameters for phenotypes available from LAB and FIELD datasets. To do so, a set of bivariate animal models were run, and genetic correlations between LAB and FIELD phenotypes were estimated using FIELD information obtained at the population level. Heritability estimates were generally higher for FIELD predictions than for the corresponding LAB measures. The additive genetic correlations (r_a) between LAB and FIELD phenotypes had different magnitudes across traits but were generally strong. Overall, these results demonstrated the potential of using FIELD information as indicator traits for the indirect genetic improvement of LAB measures. In the second approach, we included genotype information for 1,011 cows from the LAB dataset, 1,493 cows from the FIELD dataset, 181 sires with daughters in both LAB and FIELD datasets, and 540 sires with daughters in the FIELD dataset only. Predictions were obtained using the single-step GBLUP method. A four fold cross-validation was used to assess the predictive ability of the different models, assessed as the ability to predict masked LAB records from daughters of progeny testing bulls. The correlation between observed and predicted LAB measures in validation was averaged over the four training-validation sets. Different sets of phenotypic information were used sequentially in cross-validation schemes: (i) LAB cows from the training set; (ii) FIELD cows from the training set; and (iii) FIELD cows from the validation set. Models that included FIELD records showed an improvement for the majority of traits. This study suggests that breeding programs for difficult-to-measure traits could be implemented using FTIR information. While these programs should use progeny testing, acceptable values of accuracy can be achieved also for bulls without phenotyped progeny. Robust calibration equations are, deemed as essential.

Effect of Feeding Adaptation of Italian Simmental Cows before Summer Grazing on Animal Behavior and Milk Characteristics

Simple Summary

The traditional transhumant system of rearing dairy cows in mountain areas expects animals to remain indoors in the valley during the cold season, whereas during the summer they are moved to pastures at progressively higher altitudes. The animals transferred from the valley farm to the alpine pasture must adapt to various management changes. This study aimed to evaluate whether a gradual inclusion of fresh grass in the diet of dairy cows in the valley farm can improve the performance and milk characteristics during summer grazing. Three groups of six animals each were considered: one group was kept in the stable, one was transferred from the valley to the summer farm without adaptation, and the other was progressively adapted to grazing with a feeding adaptation period. Compared to animals kept indoors, grazing animals had similar performance and milk characteristics, higher rumination time and, with respect to volatile compounds in milk, higher concentrations of alcohols, aldehydes, hydrocarbons, and ketones but lower concentrations of organic acids, phenolic compounds, and dimethyl sulfone, regardless of the feeding adaptation. In conclusion, the gradual inclusion of fresh grass in the diet in the valley farm did not improve the performance and milk characteristics during summer grazing.

Abstract

According to the alpine transhumance system, dairy cows are moved from indoor feeding with conserved forage to fresh herbage feeding on pasture. The aim of this study was to assess, as a feeding adaptation technique, the effect of a gradual inclusion of fresh herbage in the diet of Italian Simmental dairy cows before their transfer to alpine pasture on performance, behavior, and milk characteristics. Eighteen cows were assigned to three groups: animals transferred to alpine pasture with a 10-d feeding adaptation period consisting in gradual access to a pasture close to the valley farm (GT), animals transferred to alpine pasture without a feeding adaptation period (AT), and animals kept in the valley farm (IND). During the first two weeks of summer grazing, GT and AT showed higher rumination time and different concentrations of ketones, hydrocarbons, organic acids, toluene, alcohols, phenols, and dimethyl sulfone in milk as compared to IND, whereas no differences were found in milk yield, composition, or coagulation properties. No differences between GT and AT were evident for the studied variables. The feeding adaptation technique used in this study did not influence the performance and milk characteristics of Italian Simmental dairy cows grazing on alpine pasture.

Effects of Summer Transhumance of Dairy Cows to Alpine Pastures on Body Condition, Milk Yield and Composition, and Cheese Making Efficiency

Simple Summary

Summer transhumance of dairy cows is a seasonal pastoral system practiced in many European countries from ancient times. This practice provides additional forage supply for mountain dairy farms and plays a role in the preservation of landscape, biodiversity, and natural habitats and conservation of local traditional dairy products, but it may affect cows’ physiological and nutritional status. This study aimed to investigate the effects of transhumance of Brown Swiss cows to summer pastures on the yield, composition, and coagulation properties of milk, and on cheese yield. For this study, twelve multiparous cows from a mountain lowland permanent farm were divided into two groups of six cows: One group stayed at the permanent farm while the other moved to the alpine pasture (1860 m above sea level). Cows at the alpine pasture had reduced milk yield and body condition, and greater fat and lower protein contents in milk compared to cows at the permanent farm. Conversely, neither milk coagulation properties nor cheese yield were affected by summer transhumance. In conclusion, summer transhumance did not affect the cheese making efficiency of milk compared to permanent farm, but the negative effect on milk yield depressed daily cheese yield, which was 2 kg/d lower in cows moved to Alpine pasture.

Abstract

Summer transhumance to alpine pastures (ALP) is widespread in dairy systems of alpine regions. This study aimed to investigate the effects of transhumance of Brown Swiss cows to ALP on the yield, composition, and coagulation properties of milk (MCP), and on cheese yield (CY). The study involved 12 multiparous cows kept at a mountain lowland permanent farm (PF), which were divided into two equal groups: One remained at the PF, the other was moved to the ALP (1860 m above sea level) from July to September. Every month (June to October), daily milk yield (MY) and body condition score (BCS) were recorded, and individual milk samples (n = 60, 2000 mL each) were collected to assess milk composition, MCP, and CY. Compared with PF, ALP cows had a reduced MY and BCS, which was maintained on return to the PF, greater fat and lower protein contents of milk. Neither MCP nor CY were affected by summer transhumance. In conclusion, summer transhumance did not affect the cheese making efficiency of milk but depressed MY and consequently daily cheese yield, which was nearly 2 kg/d lower for the ALP than the PF cows and was only partially recovered after returning to the PF in autumn.

Inferring individual cow effects, dairy system effects and feeding effects on latent variables underlying milk protein composition and cheese-making traits in dairy cattle

We examined the latent structure of 26 cheese related phenotypes in dairy cattle. Traits related to milk yield and quality (8 traits), milk protein fractions (8 traits), coagulation and curd firmness indicators (CF, 5 traits) and cheese-making phenotypes (cheese yields (%CY) and nutrient recoveries in the curd (REC), 5 traits) were analysed through multivariate factor analysis (MFA) using a varimax rotation. All phenotypes were measured in 1264 Brown Swiss cows. Ten mutual orthogonal, latent variables (factors; Fs) were obtained explaining 74% of the original variability. These Fs captured basic concepts of the cheese-making process. More precisely, the first 4 Fs, sorted by variance explained, were able to capture the underlying structure of the CY percentage (F1: %CY), the CF process with time (F2: CFt), the milk and solids yield (F3: Yield) and the presence of nitrogen (N) in the cheese (F4: Cheese N). Moreover, 4 Fs (F5: as1-β-CN, F7: κ-β-CN, F8: as2-CN and F9: as1-CN-Ph) were related to the basic milk caseins and 1 factor was associated with the α-LA whey protein (F10: α-LA). A factor describing udder health status (F6: Udder health), mainly loaded on lactose, other nitrogen compounds in the milk and SCS, was also obtained. Further, we inferred the effects of some potential sources of variation (e.g. stage of lactation and parity) including feeding and management systems. Stage of lactation had a significant effect for 7 of the 10 Fs, followed by parity of the cow (3 Fs), dairy system and feeding (3 Fs). Our work demonstrates the usefulness of MFA in reducing a large number of variables to a few latent factors with biological meaning and representing groups of traits that describe a complex process like cheese-making. Such an approach would be a valuable tool for studying the influence of different production environments and individual animal factors on protein composition and cheese-making related traits.

Genome-wide association study for cheese yield and curd nutrient recovery in dairy cows

Cheese production and consumption are increasing in many countries worldwide. As a result, interest has increased in strategies for genetic selection of individuals for technological traits of milk related to cheese yield (CY) in dairy cattle breeding. However, little is known about the genetic background of a cow's ability to produce cheese. Recently, a relatively large panel (1,264 cows) of different measures of individual cow CY and milk nutrient and energy recoveries in the cheese (REC) became available. Genetic analyses showed considerable variation for CY and for aptitude to retain high proportions of fat, protein, and water in the coagulum. For the dairy industry, these characteristics are of major economic importance. Nevertheless, use of this knowledge in dairy breeding is hampered by high costs, intense labor requirement, and lack of appropriate technology. However, in the era of genomics, new possibilities are available for animal breeding and genetic improvement. For example, identification of genomic regions involved in cow CY might provide potential for marker-assisted selection. The objective of this study was to perform genome-wide association studies on different CY and REC measures. Milk and DNA samples from 1,152 Italian Brown Swiss cows were used. Three CY traits expressing the weight (wt) of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed, and 4 REC (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY, calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk) were analyzed. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2. Single marker regressions were fitted using the GenABEL R package (genome-wide association using mixed model and regression-genomic control). In total, 103 significant associations (88 single nucleotide polymorphisms) were identified in 10 chromosomes (2, 6, 9, 11, 12, 14, 18, 19, 27, 28). For REC_FAT and REC_PROTEIN, high significance peaks were identified in Bos taurus autosome (BTA) 6 and BTA11, respectively. Marker ARS-BFGL-NGS-104610 (∼104.3 Mbp) was highly associated with REC_PROTEIN and Hapmap52348-rs29024684 (∼87.4 Mbp), closely located to the casein genes on BTA6, with REC_FAT. Genomic regions identified may enhance marker-assisted selection in bovine cheese breeding beyond the use of protein (casein) and fat contents, whereas new knowledge will help to unravel the genomic background of a cow's ability for cheese production.

From cow to cheese: genetic parameters of the flavour fingerprint of cheese investigated by direct-injection mass spectrometry (PTR-ToF-MS)

BACKGROUND

Volatile organic compounds determine important quality traits in cheese. The aim of this work was to infer genetic parameters of the profile of volatile compounds in cheese as revealed by direct-injection mass spectrometry of the headspace gas from model cheeses that were produced from milk samples from individual cows.

METHODS

A total of 1075 model cheeses were produced using raw whole-milk samples that were collected from individual Brown Swiss cows. Single spectrometry peaks and a combination of these peaks obtained by principal component analysis (PCA) were analysed. Using a Bayesian approach, we estimated genetic parameters for 240 individual spectrometry peaks and for the first ten principal components (PC) extracted from them.

RESULTS

Our results show that there is some genetic variability in the volatile compound fingerprint of these model cheeses. Most peaks were characterized by a substantial heritability and for about one quarter of the peaks, heritability (up to 21.6%) was higher than that of the best PC. Intra-herd heritability of the PC ranged from 3.6 to 10.2% and was similar to heritabilities estimated for milk fat, specific fatty acids, somatic cell count and some coagulation parameters in the same population. We also calculated phenotypic correlations between PC (around zero as expected), the corresponding genetic correlations (from -0.79 to 0.86) and correlations between herds and sampling-processing dates (from -0.88 to 0.66), which confirmed that there is a relationship between cheese flavour and the dairy system in which cows are reared.

CONCLUSIONS

This work reveals the existence of a link between the cow's genetic background and the profile of volatile compounds in cheese. Analysis of the relationships between the volatile organic compound (VOC) content and the sensory characteristics of cheese as perceived by the consumer, and of the genetic basis of these relationships could generate new knowledge that would open up the possibility of controlling and improving the sensory properties of cheese through genetic selection of cows. More detailed investigations are necessary to connect VOC with the sensory properties of cheese and gain a better understanding of the significance of these new phenotypes.

Breed of cow and herd productivity affect milk composition and modeling of coagulation, curd firming, and syneresis

Milk coagulation properties (MCP) have been widely investigated in the past using milk collected from different cattle breeds and herds. However, to our knowledge, no previous studies have assessed MCP in individual milk samples from several multi-breed herds characterized by either high or low milk productivity, thereby allowing the effects of herd and cow breed to be evaluated independently. Multi-breed herds (n=41) were classified into 2 categories based on milk productivity (high vs. low), defined according to the average milk net energy yielded daily by lactating cows. Milk samples were taken from 1,508 cows of 6 different breeds: 3 specialized dairy (Holstein-Friesian, Brown Swiss, Jersey) and 3 dual-purpose (Simmental, Rendena, Alpine Grey) breeds, and analyzed in duplicate (3,016 tests) using 2 lactodynamographs to obtain 240 curd firming (CF) measurements over 60min (1 every 15 s) for each duplicate. The 5 traditional single-point MCP (RCT, k₂₀, a₃₀, a₄₅, and a₆₀) were yielded directly by the instrument from the available CF measures. All 240 CF measures of each replicate were also used to estimate 4 individual equation parameters: RCT estimated according to curd firm change over time modeling (RCT_eq), asymptotic potential curd firmness (CF_P), curd firming instant rate constant (k_CF), and syneresis instant rate constant (k_SR) and 2 derived traits: maximum curd firmness achieved within 45min (CF_max) and time at achievement of CF_max (t_max) by curvilinear regression using a nonlinear procedure. Results showed that the effect of herd-date on traditional and modeled MCP was modest, ranging from 6.1% of total variance for k₂₀ to 10.7% for RCT, whereas individual animal variance was the highest, ranging from 32.0% for t_max to 82.5% for RCT_eq. The repeatability of MCP was high (>80%) for all traits except those associated with the last part of the lactodynamographic curve (i.e., a₆₀, k_SR, k_CF, and t_max: 57 to 71%). Reproducibility, taking into account the effect of instrument, was equal to or slightly lower than repeatability. Milk samples collected in farms characterized by high productivity exhibited delayed coagulation (RCT_eq: 18.6 vs. 16.3min) but greater potential curd firmness (CF_P: 76.8 vs. 71.9mm) compared with milk samples collected from low-productivity herds. Parity and days in milk influenced almost all MCP. Large differences in all MCP traits were observed among breeds, both between specialized and dual-purpose breeds and within these 2 groups of breeds, even after adjusting for milk quality and yield. Milk quality and MCP of samples from Jersey cows, and coagulation time of samples from Rendena cows were better than in milk from Holstein-Friesian cows, and intermediate results were found with the other breeds of Alpine origin. The results of this study, taking into account the intrinsic limitation of this technique, show that the effects of breed on traditional and modeled MCP are much greater than the effects of herd productivity class, parity, and DIM. Moreover, the variance in individual animals is much greater than the variance in individual herds within herd productivity class. It seems that improvement in MCP depends more on genetics (e.g., breed, selection) than on environmental and management factors.

Cheesemaking in highland pastures: Milk technological properties, cream, cheese and ricotta yields, milk nutrients recovery, and products composition

Summer transhumance of dairy cows to high Alpine pastures is still practiced in many mountainous areas. It is important for many permanent dairy farms because the use of highland pastures increases milk production and high-priced typical local dairy products often boost farm income. As traditional cheese- and ricotta-making procedures in Alpine pastures are central to this dairy system, the objective of this study was to characterize the quality and efficiency of products and their relationships with the quality and availability of grass during the grazing season. The milk from 148 cows from 12 permanent farms reared on a temporary farm located in Alpine pastures was processed every 2wk during the summer (7 cheesemakings from late June to early September). During each processing, 11 dairy products (4 types of milk, 2 by-products, 3 fresh products, and 2 ripened cheeses) were sampled and analyzed. In addition, 8 samples of fresh forage from the pasture used by the cows were collected and analyzed. At the beginning of the pasture season the cows were at 233±90d in milk, 2.4±1.7 parities, and produced 23.6±5.7kg/d of milk. The milk yield decreased with the move from permanent to temporary farms and during the entire summer transhumance, but partly recovered after the cows returned to the permanent farms. Similar trends were observed for the daily yields of fat, protein, casein, lactose, and energy, as we found no large variations in the quality of the milk, with the exception of the first period of Alpine pasture. The somatic cell counts of milk increased during transhumance, but this resulted from a concentration of cells in a lower quantity of milk rather than an increase in the total number of cells ejected daily from the udder. We noted a quadratic trend in availability of forage (fresh and dry matter weight per hectare), with a maximum in late July. The quality of forage also varied during the summer with a worsening of chemical composition. The evening milk (before and after natural creaming), the whole morning milk, and the mixed vat milk had different chemical compositions, traditional coagulation properties, and curd-firming modeling parameters. These variations over the pasture season were similar to the residual variations with respect to chemical composition, and much lower with respect to coagulation and curd-firming traits. Much larger variations were noted in cream, cheese, and ricotta yields, as well as in nutrient recoveries in curd during the pasture season. The protein content of forage was correlated with some of the coagulation and curd-firming traits, the ether extract of forage was positively correlated with milk fat content and cheese yields, and fiber fractions of forage were unfavorably correlated with some of the chemical and technological traits. Traditional cheese- and ricotta-making procedures showed average cream, cheese, and ricotta yields of 6.3, 14.2, and 4.9%, respectively, and an overall recovery of almost 100% of milk fat, 88% of milk protein, and 60% of total milk solids.

The nonlinear effect of somatic cell count on milk composition, coagulation properties, curd firmness modeling, cheese yield, and curd nutrient recovery

The aim of this study was to investigate the relationships between somatic cell count (SCC) in milk and several milk technological traits at the individual cow level. In particular, we determined the effects of very low to very high SCC on traits related to (1) milk yield and composition; (2) coagulation properties, including the traditional milk coagulation properties (MCP) and the new curd firming model parameters; and (3) cheese yield and recovery of milk nutrients in the curd (or loss in the whey). Milk samples from 1,271 Brown Swiss cows from 85 herds were used. Nine coagulation traits were measured: 3 traditional MCP [rennet coagulation time (RCT, min), curd firming rate (k20, min), and curd firmness after 30 min (a30, mm)] and 6 new curd firming and syneresis traits [potential asymptotic curd firmness at infinite time (CFP, mm), curd firming instant rate constant (kCF, % × min(-1)), syneresis instant rate constant (kSR, % × min(-1)), rennet coagulation time estimated using the equation (RCTeq, min), maximum curd firmness achieved within 45 min (CFmax, mm), and time at achievement of CFmax (tmax, min)]. The observed cheese-making traits included 3 cheese yield traits (%CYCURD, %CYSOLIDS, and %CYWATER, which represented the weights of curd, total solids, and water, respectively, as a percentage of the weight of the processed milk) and 4 nutrient recoveries in the curd (RECFAT, RECPROTEIN, RECSOLIDS, and RECENERGY, which each represented the percentage ratio between the nutrient in the curd and milk). Data were analyzed using a linear mixed model with the fixed effects of days in milk, parity, and somatic cell score (SCS), and the random effect of herd-date. Somatic cell score had strong influences on casein number and lactose, and also affected pH; these were traits characterized by a quadratic pattern of the data. The results also showed a negative linear relationship between SCS and milk yield. Somatic cell score influenced almost all of the tested coagulation traits (both traditional and modeled), with the exceptions of k20, CFP, and kSR. Gelation was delayed when the SCS decreased (slightly) and when it increased (strongly) with respect to a value of 2, as confirmed by the quadratic patterns observed for both RCT and RCTeq. The SCS effect on a30 showed a quadratic pattern almost opposite to that observed for RCT. With respect to the CFt parameters, kCF decreased linearly as SCS increased, resulting in a linear decrease of CFmax and a quadratic pattern for tmax. Milk SCS attained significance for %CYCURD, %CYWATER, and RECPROTEIN. As the SCS increased beyond 3, we observed a progressive quadratic decrease of the water retained in the curd (%CYWATER), which caused a parallel decrease in %CYCURD. With respect to RECPROTEIN, the negative effect of SCS was almost linear. Recovery of fat and (consequently) RECENERGY was characterized by a more evident quadratic trend, with the most favorable values associated with an intermediate SCS. Together, our results confirmed that high SCS has a negative effect on milk composition and technological traits, highlighting the nonlinear trends of some traits across the different classes of SCS. Moreover, we report that a very low SCS has a negative effect on some technological traits of milk.

Proton transfer reaction time-of-flight mass spectrometry: A high-throughput and innovative method to study the influence of dairy system and cow characteristics on the volatile compound fingerprint of cheeses

The aim of this work was to study the effect of dairy system and individual cow-related factors on the volatile fingerprint of a large number of individual model cheeses analyzed by proton transfer reaction time-of-flight mass spectrometry (PTR-ToF-MS). A total of 1,075 model cheeses were produced using milk samples collected from individual Brown Swiss cows reared in 72 herds located in mountainous areas of Trento province (Italy). The herds belonged to 5 main dairy systems ranging from traditional to modern and the cows presented different daily milk yields (24.6±7.9kg × d(-1)), stages of lactation (199±138 d in milk), and parities (2.7±1.8). The PTR-ToF-MS revealed 619 peaks, of which the 240 most intense were analyzed, and 61 of these were tentatively attributed to relevant volatile organic compounds on the basis of their fragmentation patterns and data from the literature. Principal component analysis was used to convert the multiple responses characterizing the PTR-ToF-MS spectra into 5 synthetic variables representing 62% of the total information. These principal components were related to groups of volatile compounds tentatively attributed to different peaks and used to investigate the relationship of the volatile compound profile obtained by PTR-ToF-MS to animal and farm characteristics. Lactation stage is related to 4 principal components which brought together 52.9% of the total variance and 57.9% of the area of analyzed peaks. In particular, 2 principal components were positively related to peaks tentatively attributed to aldehydes and ketones and negatively related to alcohols, esters, and acids, which displayed a linear increase during lactation. The second principal component was affected by dairy system; it was higher in the modern system in which cows received total mixed rations. The third principal component was positively related to daily milk production. In summary, we report the first application of this innovative, high-throughput technique to study the effects of dairy system and individual animal factors on volatile organic compounds of model cheeses. Individual cheesemaking procedures together with this spectrometric technique open new avenues for genetic selection of dairy species with respect to both milk and cheese quality.

Bayesian regression models outperform partial least squares methods for predicting milk components and technological properties using infrared spectral data

The aim of this study was to assess the performance of Bayesian models commonly used for genomic selection to predict "difficult-to-predict" dairy traits, such as milk fatty acid (FA) expressed as percentage of total fatty acids, and technological properties, such as fresh cheese yield and protein recovery, using Fourier-transform infrared (FTIR) spectral data. Our main hypothesis was that Bayesian models that can estimate shrinkage and perform variable selection may improve our ability to predict FA traits and technological traits above and beyond what can be achieved using the current calibration models (e.g., partial least squares, PLS). To this end, we assessed a series of Bayesian methods and compared their prediction performance with that of PLS. The comparison between models was done using the same sets of data (i.e., same samples, same variability, same spectral treatment) for each trait. Data consisted of 1,264 individual milk samples collected from Brown Swiss cows for which gas chromatographic FA composition, milk coagulation properties, and cheese-yield traits were available. For each sample, 2 spectra in the infrared region from 5,011 to 925 cm(-1) were available and averaged before data analysis. Three Bayesian models: Bayesian ridge regression (Bayes RR), Bayes A, and Bayes B, and 2 reference models: PLS and modified PLS (MPLS) procedures, were used to calibrate equations for each of the traits. The Bayesian models used were implemented in the R package BGLR (http://cran.r-project.org/web/packages/BGLR/index.html), whereas the PLS and MPLS were those implemented in the WinISI II software (Infrasoft International LLC, State College, PA). Prediction accuracy was estimated for each trait and model using 25 replicates of a training-testing validation procedure. Compared with PLS, which is currently the most widely used calibration method, MPLS and the 3 Bayesian methods showed significantly greater prediction accuracy. Accuracy increased in moving from calibration to external validation methods, and in moving from PLS and MPLS to Bayesian methods, particularly Bayes A and Bayes B. The maximum R(2) value of validation was obtained with Bayes B and Bayes A. For the FA, C10:0 (% of each FA on total FA basis) had the highest R(2) (0.75, achieved with Bayes A and Bayes B), and among the technological traits, fresh cheese yield R(2) of 0.82 (achieved with Bayes B). These 2 methods have proven to be useful instruments in shrinking and selecting very informative wavelengths and inferring the structure and functions of the analyzed traits. We conclude that Bayesian models are powerful tools for deriving calibration equations, and, importantly, these equations can be easily developed using existing open-source software. As part of our study, we provide scripts based on the open source R software BGLR, which can be used to train customized prediction equations for other traits or populations.

Selection for milk coagulation properties predicted by Fourier transform infrared spectroscopy in the Italian Holstein-Friesian breed

Milk coagulation is based on a series of physicochemical changes at the casein micelle level, resulting in formation of a gel. Milk coagulation properties (MCP) are relevant for cheese quality and yield, important factors for the dairy industry. They are also evaluated in herd bulk milk to reward or penalize producers of Protected Designation of Origin cheeses. The economic importance of improving MCP justifies the need to account for this trait in the selection process. A pilot study was carried out to determine the feasibility of including MCP in the selection schemes of the Italian Holstein. The MCP were predicted in 1,055 individual milk samples collected in 16 herds (66 ± 24 cows per herd) located in Brescia province (northeastern Italy) by means of Fourier transform infrared (FTIR) spectroscopy. The coefficient of determination of prediction models indicated moderate predictions for milk rennet coagulation time (RCT=0.65) and curd firmness (a₃₀=0.68), and poor predictions for curd-firming time (k₂₀=0.49), whereas the range error ratio (8.9, 6.9, and 9.5 for RCT, k₂₀, and a₃₀, respectively) indicated good practical utility of the predictive models for all parameters. Milk proteins were genotyped and casein haplotypes (αS₁-, β-, αS₂-, and κ-casein) were reconstructed. Data from 51 half-sib families (19.9 ± 16.4 daughters per sire) were analyzed by an animal model to estimate (1) the genetic parameters of predicted RCT, k₂₀, and a₃₀; (2) the breeding values for these predicted clotting variables; and (3) the effect of milk protein genotypes and casein haplotypes on predicted MCP (pMCP). This is the first study to estimate both genetic parameters and breeding values of pMCP, together with the effects of milk protein genotypes and casein haplotypes, that also considered k₂₀, probably the most important parameter for the dairy industry (because it indicates the time for the beginning of curd-cutting). Heritability of predicted RCT (0.26) and k₂₀ (0.31) were close to the average heritability described in literature, whereas the heritability of a₃₀ was higher (0.52 vs. 0.27). The effects of milk proteins were statistically significant and similar to those obtained on measured MCP. In particular, haplotypes including uncommon variants showed positive (B-I-A-B) or negative (B-A(1)-A-E) effects. Based on these findings, FTIR spectroscopy-pMCP is proposed as a potential selection criterion for the Italian Holstein.