Combined Metabolomics and Biochemical Analyses of Serum and Milk Revealed Parity-Related Metabolic Differences in Sanhe Dairy Cattle

The production performance of dairy cattle is closely related to their metabolic state. This study aims to provide a comprehensive understanding of the production performance and metabolic features of Sanhe dairy cattle across different parities, with a specific focus on evaluating variations in milk traits and metabolites in both milk and serum. Sanhe dairy cattle from parities 1 to 4 (S1, n = 10; S2, n = 9; S3, n = 10; and S4, n = 10) at mid-lactation were maintained under the same feeding and management conditions. The milk traits, hydrolyzed milk amino acid levels, serum biochemical parameters, and serum free amino acid levels of the Sanhe dairy cattle were determined. Multiparous Sanhe dairy cattle (S2, S3, and S4) had a greater milk protein content, lower milk lactose content, and lower solids-not-fat content than primiparous Sanhe dairy cattle (S1). Moreover, S1 had a higher ratio of essential to total amino acids (EAAs/TAAs) in both the serum and milk. The serum biochemical results showed the lower glucose and total protein levels in S1 cattle were associated with milk quality. Furthermore, ultra-high-resolution high-performance liquid chromatography with tandem MS analysis (UPLC-MS/MS) identified 86 and 105 differential metabolites in the serum and milk, respectively, and these were mainly involved in amino acid, carbohydrate, and lipid metabolism. S1 and S2/S3/S4 had significantly different metabolic patterns in the serum and milk, and more vitamin B-related metabolites were significantly higher identified in S1 than in multiparous cattle. Among 36 shared differential metabolites in the serum and milk, 10 and 7 metabolites were significantly and strongly correlated with differential physiological indices, respectively. The differential metabolites identified were enriched in key metabolic pathways, illustrating the metabolic characteristics of the serum and milk from Sanhe dairy cattle of different parities. L-phenylalanine, dehydroepiandrosterone, and linoleic acid in the milk and N-acetylornithine in the serum could be used as potential marker metabolites to distinguish between Sanhe dairy cattle with parities of 1-4. In addition, a metabolic map of the serum and milk from the three aspects of carbohydrates, amino acids, and lipids was created for the further analysis and exploration of their relationships. These results reveal significant variations in milk traits and metabolites across different parities of Sanhe dairy cattle, highlighting the influence of parity on the metabolic profiles and production performance. Tailored nutritional strategies based on parity-specific metabolic profiles are recommended to optimize milk production and quality in Sanhe cattle.

Observational study on the associations between milk yield, composition, and coagulation properties with blood biomarkers of health in Holstein cows

The considerable increase in the production capacity of individual cows owing to both selective breeding and innovations in the dairy sector has posed challenges to management practices in terms of maintaining the nutritional and metabolic health status of dairy cows. In this observational study, we investigated the associations between milk yield, composition, and technological traits and a set of 21 blood biomarkers related to energy metabolism, liver function or hepatic damage, oxidative stress, and inflammation or innate immunity in a population of 1,369 high-yielding Holstein-Friesian dairy cows. The milk traits investigated in this study included 4 production traits (milk yield, fat yield, protein yield, daily milk energy output), 5 traits related to milk composition (fat, protein, casein, and lactose percentages and urea), 11 milk technological traits (5 milk coagulation properties and 6 curd-firming traits). All milk traits (i.e., production, composition, and technological traits) were analyzed according to a linear mixed model that included the days in milk, the parity order, and the blood metabolites (tested one at a time) as fixed effects and the herd and date of sampling as random effects. Our findings revealed that milk yield and daily milk energy output were positively and linearly associated with total cholesterol, nonesterified fatty acids, urea, aspartate aminotransferase, γ-glutamyl transferase, total bilirubin, albumin, and ferric-reducing antioxidant power, whereas they were negatively associated with glucose, creatinine, alkaline phosphatase, total reactive oxygen metabolites, and proinflammatory proteins (ceruloplasmin, haptoglobin, and myeloperoxidase). Regarding composition traits, the protein percentage was negatively associated with nonesterified fatty acids and β-hydroxybutyrate (BHB), while the fat percentage was positively associated with BHB, and negatively associated with paraoxonase. Moreover, we found that the lactose percentage increased with increasing cholesterol and albumin and decreased with increasing ceruloplasmin, haptoglobin, and myeloperoxidase. Milk urea increased with an increase in cholesterol, blood urea, nonesterified fatty acids, and BHB, and decreased with an increase in proinflammatory proteins. Finally, no association was found between the blood metabolites and milk coagulation properties and curd-firming traits. In conclusion, this study showed that variations in blood metabolites had strong associations with milk productivity traits, the lactose percentage, and milk urea, but no relationships with technological traits of milk. Specifically, increasing levels of proinflammatory and oxidative stress metabolites, such as ceruloplasmin, haptoglobin, myeloperoxidase, and total reactive oxygen metabolites, were shown to be associated with reductions in milk yield, daily milk energy output, lactose percentage, and milk urea. These results highlight the close connection between the metabolic and innate immunity status and production performance. This connection is not limited to specific clinical diseases or to the transition phase but manifests throughout the entire lactation. These outcomes emphasize the importance of identifying cows with subacute inflammatory and oxidative stress as a means of reducing metabolic impairments and avoiding milk fluctuations.

Exploring Interrelationships between Colour, Composition, and Coagulation Traits of Milk from Cows, Goats, and Sheep

This study explores the interrelationships between the composition, coagulation, and colour of sheep, goat, and cow milk to identify their similarities and differences and to assess whether the relationships between the variables are common to all species or whether they emerge from species-specific relationships. For this purpose, 2400 individual milk samples were analysed. The differences and similarities between the species were determined using discriminant analysis and cluster analysis. The results show a clear differentiation between species. Sheep milk stands out for its cheesemaking capacity and shows similarities with goat milk in composition and coagulation. Nonetheless, colorimetry highlights a greater similarity between sheep and cow milk. Composition and colorimetry were more discriminating than coagulation, and the variables that differed the most were fat, protein, curd yield, lightness, and red-green balance. Using canonical correlation, the interrelationships between the different sets of variables were explored, revealing patterns of common variation and species-specific relationships. Colorimetric variables were closely related to milk solids in all species, while in sheep milk, an inverse relationship with lactose was also identified. Furthermore, a strong relationship was revealed for all species between colour and curd yield. This could be modelled and applied to estimate the technological value of milk, proving colorimetry as a useful tool for the dairy industry.

Metabolic Changes in Serum and Milk of Holstein Cows in Their First to Fourth Parity Revealed by Biochemical Analysis and Untargeted Metabolomics

The performance of dairy cows is closely tied to the metabolic state, and this performance varies depending on the number of times the cows have given birth. However, there is still a lack of research on the relationship between the metabolic state of Holstein cows and the performance of lactation across multiple parities. In this study, biochemical analyses and metabolomics studies were performed on the serum and milk from Holstein cows of parities 1-4 (H1, N = 10; H2, N = 7; H3, N = 9; H4, N = 9) in mid-lactation (DIM of 141 ± 4 days) to investigate the link between performance and metabolic changes. The results of the milk quality analysis showed that the lactose levels were highest in H1 (p = 0.036). The total protein content in the serum increased with increasing parity (p = 0.013). Additionally, the lipase activity was found to be lowest in H1 (p = 0.022). There was no difference in the composition of the hydrolyzed amino acids in the milk among H1 to H4. However, the free amino acids histidine and glutamate in the serum were lowest in H1 and highest in H3 (p < 0.001), while glycine was higher in H4 (p = 0.031). The metabolomics analysis revealed that 53 and 118 differential metabolites were identified in the milk and serum, respectively. The differential metabolites in the cows' milk were classified into seven categories based on KEGG. Most of the differential metabolites in the cows' milk were found to be more abundant in H1, and these metabolites were enriched in two impact pathways. The differential metabolites in the serum could be classified into nine categories and enriched in six metabolic pathways. A total of six shared metabolites were identified in the serum and milk, among which cholesterol and citric acid were closely related to amino acid metabolism in the serum. These findings indicate a significant influence of blood metabolites on the energy and amino acid metabolism during the milk production process in the Holstein cows across 1-4 lactations, and that an in-depth understanding of the metabolic changes that occur in Holstein cows during different lactations is essential for precision farming, and that it is worthwhile to further investigate these key metabolites that have an impact through controlled experiments.

Metabolic and physiological adaptations to first and second lactation in Holstein dairy cows

Huge differences exist between cow yields and body sizes during their first and second lactations. The transition period is the most critical and investigated phase of the lactation cycle. We compared metabolic and endocrine responses between cows at different parities during the transition period and early lactation. Eight Holstein dairy cows were monitored at their first and second calving during which they were reared under the same conditions. Milk yield, dry matter intake (DMI), and body weight (BW) were regularly measured, and energy balance, efficiency, and lactation curves were calculated. Blood samples were collected on scheduled days from -21 d relative to calving (DRC) to 120 DRC for the assessment of metabolic and hormonal profiles (biomarkers of metabolism, mineral status, inflammation, and liver function). Large variations in the period in question for almost all variables investigated were observed. Compared with their first lactation, cows during their second lactation had higher DMI (+15%) and BW (+13%), their milk yield was greater (+26%), lactation peak was higher and earlier (36.6 kg/d at 48.8 DRC vs. 45.0 kg/d at 62.9 DRC), but persistency was reduced. Milk fat, protein, and lactose contents were higher during the first lactation and coagulation properties were better (higher titratable acidity, faster and firmer curd formation). Postpartum negative energy balance was more severe the during the second lactation (1.4-fold at 7 DRC) and plasma glucose was lower. Circulating insulin and insulin-like growth factor-1 were lower in second-calving cows during the transition period. At the same time, markers of body reserve mobilization (β-hydroxybutyrate and urea) increased. Moreover, albumin, cholesterol, and γ-glutamyl transferase were higher during second lactation, whereas bilirubin and alkaline phosphatase were lower. The inflammatory response after calving was not different, as suggested by the similar haptoglobin concentrations and only transient differences in ceruloplasmin. Blood growth hormone did not differ during the transition period but was lower during the second lactation at 90 DRC, whereas circulating glucagon was higher. These results agree with the differences in milk yield and confirmed the hypothesis of a different metabolic and hormonal status between the first and second lactation partly related to different degrees of maturity.

Efficiency of Manchega Sheep Milk Intended for Cheesemaking and Determination of Factors Causing Inefficiency

Understanding the factors that determine and regulate cheese yield would allow, through deterministic parametric efficiency models, the determination of the most appropriate milk characteristics for the industry, and the estimation of a technological value for these characteristics. The present study aims to evaluate coagulation performance of Manchega sheep milk intended for cheesemaking and explores two models to determine milk technological efficiency. For this purpose, 1200 Manchega sheep milk samples were collected, and analyses were performed for composition, milk coagulation properties (MCP), somatic cell count (SCC), and milk color values. A first model was built based on curd yield (CE) and a second one based on dry curd yield (DCE). GLM and MANCOVA analyses were used to identify the factors that determine curd yield efficiency, which mainly depended on pH, casein, and lactose content and, to a lesser extent, on the speed of coagulation and curd firmness. When comparing both models, differences were linked to the water retention capacity of the curd. Based on this, the DCE model was considered much more accurate for prediction of coagulation efficiency in a wider variety of cheeses, as it does not seem to be affected by moisture loss.

Phenotypic and genetic characterization of the occurrence of noncoagulating milk in dairy sheep

Milk coagulation ability is of central importance for the sheep dairy industry because almost all sheep milk is destined for cheese processing. The occurrence of milk with impaired coagulation properties is an obstacle to cheese processing and, in turn, to the profitability of the dairy companies. In this work, we investigated the causes of noncoagulation of sheep milk; specifically, we studied the effect of milk physicochemical properties on milk coagulation status [coagulating and noncoagulating (NC) milk samples, which do or do not coagulate within 30 min, respectively], and whether mid-infrared spectroscopy (MIR) could be used to assess variability in coagulation status. We also investigated the genetic background of milk coagulation ability. Individual milk samples were collected from 996 Sarda ewes farmed in 47 flocks located in Sardinia (Italy). Considered traits were daily milk yield, milk composition traits, and milk coagulation properties (rennet coagulation time, curd firming time, and curd firmness), and MIR spectra were acquired. About 9% of samples did not coagulate within 30 min. A logistic regression approach was used to test the effect of milk-related traits on milk coagulation status. A principal component (PC) analysis was carried out on the milk MIR spectra, and PC scores were then used as covariates in a logistic regression model to assess their relationship with milk coagulation status. Results of the present work demonstrated that the probability of having NC samples increases as milk contents of proteins and chlorides and somatic cell score increase. The analysis of PC extracted from milk spectra that influenced coagulation status highlighted key regions associated with lactose and protein concentrations, and others not associated with routinely collected milk composition traits. These results suggest that the occurrence of NC is mostly related to damage of the epithelium secretory mammary cells, which occurs with the advancement of a lactation or due to unhealthy mammary gland status. Genetic analysis of milk coagulation status and of the extracted PC confirmed the genetic background of the milk coagulability of sheep milk.

Physicochemical Characteristics, Fatty Acids Profile and Lipid Oxidation during Ripening of Graviera Cheese Produced with Raw and Pasteurized Milk

The production of cheese can be made from either pasteurized or non-pasteurized milk, depending on the country or dietary habits. In this work, the effect of pasteurization of milk on the progress of the physicochemical properties, fatty acids profile and lipid oxidation of cheese throughout a maturation period of 90 days is presented. This research was carried out on two types of Graviera cheese produced in Greece, one made from raw milk and the other from pasteurized milk. The proximal composition of each sample was evaluated, the fatty acids profile was analyzed by Gas Chromatography, whereas lipid oxidation was determined on the basis of the formation of malondialdehyde (MDA). Significant differences (p < 0.05) in the values of pH, fat and density between raw and pasteurized milk were observed. The physicochemical parameters during the ripening of the cheeses showed significant differences according to the type and the stage of maturation. Specifically, the two types of cheese differed significantly (p < 0.05) in terms of pH, protein, fat in dry matter (FDM), and water-soluble nitrogen/total nitrogen (WSN/TN). Although the fatty acids profile was similar for the two types of cheese, differences were observed during the ripening stages as well as between the milk and the final product. The lipid oxidation levels increased during maturation, whereas they seemed to be lower in the pasteurized cheeses. Therefore, it can be concluded that the use of raw or pasteurized milk affects the physicochemical characteristics, fatty acids profile and lipid oxidation of Graviera cheese during ripening.

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.

Chemical Composition, Fatty Acid Profile, and Lipid Quality Indices in Commercial Ripening of Cow Cheeses from Different Seasons

The aim of the study was to compare and demonstrate whether commercial rennet ripening cheeses available on the market in summer and winter differ in their chemical composition, fatty acid profile, content of cis9trans11 C18:2 (CLA) acid and other trans isomers of C18:1 and C18:2 acid and whether they are characterized by different values of lipid quality assessment indices. The experimental material consisted of rennet ripening of cheeses produced from cow's milk available in the Polish market. The first batch contained cheeses produced in winter and purchased from the market between May and June. The second batch contained cheeses produced in summer and purchased between November and December. Chemical composition was analyzed by FoodScan apparatus. The gas chromatography (GC) method was used to determine the content of fatty acids. Results obtained in the presented study indicate that the chemical composition, content of fatty acids trans isomers, and lipid quality indices varied between summer and winter cheeses. The summer cheeses were richer sources of MUFA and PUFA compared to winter cheeses. Summer cheeses were also characterized by lower content of SFA, higher content n - 3, lower n - 6/n - 3 ratio, and higher content of DFA. Higher contents of CLA and trans C18:1 and C18:2 were found in summer cheeses.

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.

Multiple Breeds and Countries' Predictions of Mineral Contents in Milk from Milk Mid-Infrared Spectrometry

Measuring the mineral composition of milk is of major interest in the dairy sector. This study aims to develop and validate robust multi-breed and multi-country models predicting the major minerals through milk mid-infrared spectrometry using partial least square regressions. A total of 1281 samples coming from five countries were analyzed to obtain spectra and in ICP-AES to measure the mineral reference contents. Models were built from records coming from four countries (n = 1181) and validated using records from the fifth country, Austria (n = 100). The importance of including local samples was tested by integrating 30 Austrian samples in the model while validating with the remaining 70 samples. The best performances were achieved using this second set of models, confirming the need to cover the spectral variability of a country before making a prediction. Validation root mean square errors were 54.56, 63.60, 7.30, 59.87, and 152.89 mg/kg for Na, Ca, Mg, P, and K, respectively. The built models were applied on the Walloon milk recording large-scale spectral database, including 3,510,077. The large-scale predictions on this dairy herd improvement database provide new insight regarding the minerals' variability in the population, as well as the effect of parity, stage of lactation, breeds, and seasons.

Graduate Student Literature Review: Current understanding of the influence of on-farm factors on bovine raw milk and its suitability for cheesemaking

Relationships between dairy farm practices, the composition and properties of raw milk, and the quality of the resulting cheese are complex. In this review, we assess the effect of farm factors on the quality of bovine raw milk intended for cheesemaking. The literature reports several prominent farm-related factors that are closely associated with milk quality characteristics. We describe their effects on the composition and technological properties of raw milk and on the quality of the resulting cheese. Cow breed, composite genotype, and protein polymorphism all have noticeable effects on milk coagulation, cheese yield, and cheese composition. Feed and feeding strategy, dietary supplementation, housing and milking system, and seasonality of milk production also influence the composition and properties of raw milk, and the resulting cheese. The microbiota in raw milk is influenced by on-farm factors and by the production environment, and may influence the technological properties of the milk and the sensory profile of certain cheese types. Advances in research dealing with the technological properties of raw milk have undoubtedly improved understanding of how on-farm factors affect milk quality attributes, and have refuted the concept of one milk for all purposes. The specific conditions for milk production should be considered when the milk is intended for the production of cheese with unique characteristics. The scientific identification of these conditions would improve the current understanding of the complex associations between raw milk quality and farm and management factors. Future research that considers dairy landscapes within broader perspectives and develops multidimensional approaches to control the quality of raw milk intended for long-ripening cheese production is recommended.

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.

Derivation of multivariate indices of milk composition, coagulation properties, and curd yield in Manchega dairy sheep

This study approaches the interrelation patterns between composition of milk and whey, curd yield, chromaticity, syneresis, and technological quality of Manchega sheep milk using multivariate factor analysis. In addition, the effect of the main husbandry components (flock, prolificacy, season of the year, stage of lactation, and parity) on the common latent factors that define the pattern of variation of Manchega milk was assessed. For this purpose, 1,200 individual Manchega ewe milk samples from 4 different flocks registered under the Protected Designation of Origin Queso Manchego were analyzed (50 ewes/flock). Samples were collected in 2 different seasons of the year (spring and autumn) and at 3 time points per season: early, mid-, and late lactation. The obtained results suggested that curd yield mainly depends on milk composition, and the retention of water in the curd is related to coagulation traits. Thus, composition and moisture content could be useful indicators to assess the efficiency and quality of milk intended for cheesemaking, regardless of the analysis of coagulation properties. Finally, in terms of husbandry, a direct effect of flock and stage of lactation was observed on all analyzed factors, with a lower influence of season and parity.

Coagulation process in Manchega sheep milk from Spain: A path analysis approach

Characteristics of sheep milk are of great interest for the dairy industry, as almost the totality of production is intended for cheesemaking. However, the existing relationships between these variables are complex. This study assessed composition, hygienic quality, coagulation properties, and curd yield of 1,200 individual Manchega sheep milk samples. The aim was to compare the effect of composition and hygienic quality on coagulation and curdling, and to evaluate the relationship between curd yields and the coagulation process and the effect of other features by using path analysis methodologies. Outcomes proved path analysis to be a useful and effective tool to assess these relationships through direct and indirect paths within the same model. Results showed that the factors that had a direct influence on milk coagulation were lactose concentration, casein content, and initial pH of milk. Contrastingly, somatic cells did not seem to have any effect (direct or indirect) on the coagulation process. Factors that directly affected curd yield were fat content, lactose concentration, casein content, and curd moisture. However, technological parameters showed little effect over curd yield.

Physicochemical, Microbiological and Technological Properties of Red Deer (Cervus elaphus) Milk during Lactation

Simple Summary

Milk from red deer is richer in fat and proteins than that of cow or other ruminants. The semi-captive breeding of this species has traditionally focused on meat, velvet or hunting purposes, but recent studies suggested that the high level of nutrients, the promising content of bioactive peptides and the better digestibility than that of milk from other species could open innovative alternatives for the dairy industry. As for other non-commercial mammalian species that are gaining technological interest for the elaboration of dairy products, it is necessary to understand the aptitude and performance of milk from red deer to be used for the production of cheese, fermented milks or other products. Our study aims to assess some chemical, physical, microbiological and technological properties of red deer milk during a lactation period of 18 weeks. The results show that milk from this species is similar to that of other ruminant species whose milk is commercialized. In addition, our results indicate the best period to industrialize the milk during lactation. To the best of our knowledge, this is the first study to explore the benefits of using red deer milk with a technological approach.

Abstract

This study describes chemical, physical, microbiological and technological characteristics of red deer milk and the effect of lactation on these parameters in order to know their potential aptitude to elaborate dairy products. During 18 weeks, milk from five hinds was monitored for composition, bacteriology, somatic cell count (SCC), physical properties and rennet coagulation. Mean values (g/100 g) for fat, protein, lactose and dry matter were 10.4, 7.1, 4.3 and 24.2, respectively, and for urea, 265 mg/100 mL. Except for lactose, a significant increase in these components was observed (p < 0.01) as lactation progressed. The average values for bacteriology and SCC were 5.3 log cfu/mL and 4.7 log cells/mL, respectively. Regarding physical properties, conductivity (mean: 2.8 ms/cm), viscosity (3.1 Cp), coordinates L* (89.9) and a* (−3.1) and milk fat globule diameter (D_4,3: 6.1 µm) increased along with lactation while density (1.038 g/mL) decreased (p < 0.01). The pH (6.7), acidity (22.9° Dornic), coordinate b* (8.4) and ethanol stability (66.6% v/v) were stable during the study period. The stage of lactation also has a significant impact on milk coagulation properties and mean curd yield was 3.29 g/10 mL. These results suggest that red deer milk could be a potential innovative source of milk for the dairy industry.

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.

Goat farm variability affects milk Fourier-transform infrared spectra used for predicting coagulation properties

Driven by the large amount of goat milk destined for cheese production, and to pioneer the goat cheese industry, the objective of this study was to assess the effect of farm in predicting goat milk-coagulation and curd-firmness traits via Fourier-transform infrared spectroscopy. Spectra from 452 Sarda goats belonging to 14 farms in central and southeast Sardinia (Italy) were collected. A Bayesian linear regression model was used, estimating all spectral wavelengths' effects simultaneously. Three traditional milk-coagulation properties [rennet coagulation time (min), time to curd firmness of 20 mm (min), and curd firmness 30 min after rennet addition (mm)] and 3 curd-firmness measures modeled over time [rennet coagulation time estimated according to curd firmness change over time (RCT_eq), instant curd-firming rate constant, and asymptotical curd firmness] were considered. A stratified cross validation (SCV) was assigned, evaluating each farm separately (validation set; VAL) and keeping the remaining farms to train (calibration set) the statistical model. Moreover, a SCV, where 20% of the goats randomly taken (10 replicates per farm) from the VAL farm entered the calibration set, was also considered (SCV₈₀). To assess model performance, coefficient of determination (R²_VAL) and the root mean squared error of validation were recorded. The R²_VAL varied between 0.14 and 0.45 (instant curd-firming rate constant and RCT_eq, respectively), albeit the standard deviation was approximating half of the mean for all the traits. Although average results of the 2 SCV procedures were similar, in SCV₈₀, the maximum R²_VAL increased at about 15% across traits, with the highest observed for time to curd firmness of 20 mm (20%) and the lowest for RCT_eq (6%). Further investigation evidenced important variability among farms, with R²_VAL for some of them being close to 0. Our work outlined the importance of considering the effect of farm when developing Fourier-transform infrared spectroscopy prediction equations for coagulation and curd-firmness traits in goats.

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.

Modeling weight loss of cheese during ripening and the influence of dairy system, parity, stage of lactation, and composition of processed milk

The yield, flavor, and texture of ripened cheese result from numerous interrelated microbiological, biochemical, and physical reactions that take place during ripening. The aims of the present study were to propose a 2-compartment first-order kinetic model of cheese weight loss over the ripening period; to test the variation in new informative phenotypes describing this process; and to assess the effects on these traits of dairy farming system, individual farms within dairy system, animal factors, and milk composition. A total of 1,211 model cheeses were produced in the laboratory using individual 1.5-L milk samples from Brown Swiss cows reared on 83 farms located in Trento Province. During ripening (60 d; temperature 15°C, relative humidity 85%), the weight of all model cheeses was measured, and cheese yield (cheese weight/processed milk weight, %CY) was calculated at 7 intervals from cheese-making (0, 1, 7, 14, 28, 42, and 60 d). Using these measures, a 2-compartment first-order kinetic model (3-parameter equation) was developed for modeling %CY during the ripening period, as follows: [Formula: see text] , where %CY_t is the %CY at ripening time t; %CY_i and %CY_f are the modeled %CY traits at time 0 d (%CY_i = initial %CY) and at the end of a ripening period sufficient to reach a constant wheel weight (%CY_f = final %CY after 60 d ripening in the case of small model cheeses); k_CY is the instant rate constant for cheese weight loss (%/d). Cheese weight and protein and fat losses were calculated as the % difference between the model cheeses at 0 and after 60 d of ripening. The variation in cheese pH was calculated as the % difference between pH at 0 and after 60 d. Dairy system, individual herd within dairy system, and the cow's parity and lactation stage (tested with a linear mixed model) strongly affected almost all the traits collected during model cheese ripening. Milk fat, protein, lactose, pH, and somatic cell score also greatly affected almost all the traits, although k_CY was affected only by milk protein. After including milk composition in the linear mixed model, the importance of all the herd and animal sources of variation was greatly reduced for all traits. The proposed model and novel traits could be tested, first, with the aim of establishing new monitoring procedures enabling the dairy industry to improve milk quality-based payment systems at the herd level and, second, with a view to exploring possible genetic improvements to dairy cow populations.

Can We Observe Expected Behaviors at Large and Individual Scales for Feed Efficiency-Related Traits Predicted Partly from Milk Mid-Infrared Spectra?

Phenotypes related to feed efficiency were predicted from records easily acquired by breeding organizations. A total of 461,036 and 354,148 records were collected from the first and second parity Holstein cows. Equations were applied to the milk mid-infrared spectra to predict the main milk components and coupled with animal characteristics to predict the body weight (pBW). Dry matter intake (pDMI) was predicted from pBW using the National Research Council (NRC) equation. The consumption index (pIC) was estimated from pDMI and fat, and protein corrected milk. All traits were modeled using single trait test-day models. Descriptive statistics were within the expected range. Milk yield, pDMI, and pBW were phenotypically positively related (r ranged from 0.08 to 0.64). As expected, pIC was phenotypically negatively correlated with milk yield (-0.77 and -0.80 for the first and second lactation) and slightly positively correlated with pBW (0.16 and 0.07 for the first and second lactation). Later, parity cows seemed to have a better feed efficiency as they had a lower pIC. Although the prediction accuracy was moderate, the observed behaviors of studied traits by year, stage of lactation, and parity were in agreement with the literature. Moreover, as a genetic component was highlighted (heritability around 0.18), it would be interesting to realize a genetic evaluation of these traits and compare the obtained breeding values with the ones estimated for sires having daughters with reference feed efficiency records.

Structural equation modeling for investigating multi-trait genetic architecture of udder health in dairy cattle

Mastitis is one of the most prevalent and costly diseases in dairy cattle. It results in changes in milk composition and quality which are indicators of udder inflammation in absence of clinical signs. We applied structural equation modeling (SEM) - GWAS aiming to explore interrelated dependency relationships among phenotypes related to udder health, including milk yield (MY), somatic cell score (SCS), lactose (%, LACT), pH and non-casein N (NCN, % of total milk N), in a cohort of 1,158 Brown Swiss cows. The phenotypic network inferred via the Hill-Climbing algorithm was used to estimate SEM parameters. Integration of multi-trait models-GWAS and SEM-GWAS identified six significant SNPs for SCS, and quantified the contribution of MY and LACT acting as mediator traits to total SNP effects. Functional analyses revealed that overrepresented pathways were often shared among traits and were consistent with biological knowledge (e.g., membrane transport activity for pH and MY or Wnt signaling for SCS and NCN). In summary, SEM-GWAS offered new insights on the relationships among udder health phenotypes and on the path of SNP effects, providing useful information for genetic improvement and management strategies in dairy cattle.

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.

Changes in whey proteome with lactation stage and parity in dairy cows using a label-free proteomics approach

Milk yield and several components of milk that are affected by physiological factors have been widely investigated. However, the effects of lactation stage and parity on bovine milk whey proteins have not been well elucidated. To aid in unraveling the proteome profile and exploring the protein biosynthesis of mammary glands, a label-free proteomic approach was used to characterize whey proteomes depending on the lactation stage and parity of dairy cows. The results of this study show that the abundances of several proteins, such as early lactation protein, syntenin, and heparanase, were associated with specific stages of the lactation cycle; this was evidenced by a principal component analysis. In addition, several proteins, such as hemoglobin subunits beta and alpha, β-lactoglobulin, CD320, and apolipoprotein E, corresponded to the parity of the dairy cows and were herein considered as useful biomarkers to distinguish different parities. Most of the differentially expressed proteins from specific lactation stages and parity milk groups were annotated in the response to stimulus and protein metabolic processes. The findings reveal that developmental changes in whey proteomes correspond to lactation stages and parities, which in turn provides new insight into the underlying implications of the production of specific proteins to meet the health benefits of offspring and host, and allow us to explore the mechanisms of protein biosynthesis in mammary glands associated with physiological changes in dairy cows.

Chemometric authentication of farming systems of origin of food (milk and ripened cheese) using infrared spectra, fatty acid profiles, flavor fingerprints, and sensory descriptions

Milk samples from 1264 cows in 85 farms were authenticated for different farming-systems using a 10-fold cross-validated linear-discriminant-analysis using Fourier-transform infrared spectra (FTIRS) and gas-chromatographic fatty-acid (FA) profiles. FTIRS gave correct classification greater than FAs (97.4% vs. 81.1%) during calibration, but slightly worse in validation (73.5% vs 77.3%) and their combination improved the results. All milk samples were processed into ripened model-cheeses, and analyzed by near-infrared-spectrometry (NIRS), by proton-transfer-reaction time-of-flight mass-spectrometry for their volatile organic compound (VOCs) fingerprint and by panel sensory profiling (SENS). Farming-system authentication on cheese samples was less efficient than on milk, but still possible. The instrumental methods yielded similar validation results, better than SENS, and their combination improved the correct classification rate. The efficiency of the different technics was affected by specific farming systems. In conclusion, dairy products could be discriminated for farming-systems with acceptable accuracy, but the methods tested differ in sampling procedure, rapidity and costs.

Detailed macro- and micromineral profile of milk: Effects of herd productivity, parity, and stage of lactation of cows of 6 dairy and dual-purpose breeds

The aim of this study was to quantify the major sources of variation in the levels of 15 minerals in individual milk samples collected from cows raised in multibreed dairy herds. The herds (n = 27) were classified into 2 categories, according to milk productivity. Milk productivity was based on the net energy of lactating cows' average daily milk yield. Milk samples were collected from 240 cows of 6 different breeds: 3 specialized dairy (Holstein-Friesian, Brown Swiss, and Jersey) and 3 dual-purpose (Simmental, Rendena, and Alpine Grey). The samples were analyzed for macro-elements (Na, Mg, P, S, K, and Ca), essential micro-elements (Mn, Fe, Cu, Zn, and Se), and environmental micro-elements (B, Si, Sr, and Sn), using inductively coupled plasma-optical emission spectrometry. Data were analyzed using a linear mixed model that included fixed effects of days in milk (DIM), parity, breed, and herd productivity, and a random effect of herd-date within productivity level. Results showed that the effect of herd-date varied across minerals. It was especially large for environmental minerals (ranging from 47 to 91% of total variance) and ranged from 11 to 61% for macrominerals and essential microminerals. Milk samples collected from farms with a high level of herd productivity had a richer mineral profile than samples from low-productivity herds. Parity only influenced macrominerals, with the exception of S and Ca, while DIM influenced almost all minerals, with a few exceptions among the environmental elements. Differences in mineral profile were small between specialized and dual-purpose breeds, but they were large within the group of the specialized cows. These breed differences were reduced after adjusting for milk quality and yield, particularly in the case of milk Mg, S, Ca, Mn, and Zn levels. Milk samples from the Jersey and Brown Swiss cows had higher mineral levels (Sn excluded) than milk from the Holstein-Friesian cows; the other breeds of Alpine origin produced milk of intermediate quality. Our findings suggest that breed has a stronger effect on macrominerals and some of the essential microminerals than herd productivity, parity, and DIM. The modification of the mineral profile in milk seems possible for many minerals, but it likely depends on genetics (e.g., breed, selection) and on environmental and management factors in variable proportions according to the mineral considered.

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.

A Study on the Effects of Rumen Acidity on Rumination Time and Yield, Composition, and Technological Properties of Milk from Early Lactating Holstein Cows

Simple Summary

The increase in milk yield achieved in recent decades by the dairy sector has been sustained by feeding dairy cows with more concentrates and less forage. This leads to increasing rumen acidity, a status widespread in high-producing dairy cows that may affect feed intake, impair ruminal digestion, and cause diarrhea, laminitis, inflammation, and liver abscesses. The effects of rumen acidity on milk yield and composition are controversial, while those on milk coagulation properties and cheese yield have not yet been explored. This study investigated whether the rumen acidity status affects rumination time, and the production, composition, coagulation properties and cheese yield of milk obtained by 100 early-lactating Holstein cows. The variation in rumen acidity was associated with changes in the cows’ rumen fluid composition and circadian pattern of rumination time. Moreover, daily milk yield linearly decreased as the rumen acidity increased. Conversely, the composition and technological properties of milk were unaffected, even when there were differences in rumen acidity, suggesting that variation in rumen acidity has little impact on cheese-making traits.

Abstract

The use of high grain rations in dairy cows is related to an increase in rumen acidity. This study investigated whether the rumen acidity status affects rumination time (RT), and the production, composition, coagulation properties (MCPs) and cheese yield (CY) of milk. One hundred early-lactating Holstein cows with no clinical signs of disease and fed total mixed rations were used. Rumen fluid was collected once from each cow by rumenocentesis to determine pH and volatile fatty acid (VFA) content. The cows were classified according to the quartile of rumen acidity (QRA), a factor defined by multivariate analysis and associated with VFA and pH. Rumen fluid pH averaged 5.61 in the first quartile and 6.42 in the fourth, and total VFA content increased linearly with increasing rumen acidity. In addition, RT increased as rumen acidity increased, but only in the daily time interval from 08:00 to 12:00. Milk yield linearly decreased as rumen acidity increased, whereas QRA did not affect pH, fat or protein contents of milk. Furthermore, the MCPs, assessed by lactodynamograph, and CY were unaffected by QRA. It is suggested that differences in rumen acidity have little influence on the nutrient content, coagulation properties and CY of milk.

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.

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.

Variation of milk coagulation properties, cheese yield, and nutrients recovery in curd of cows of different breeds before, during and after transhumance to highland summer pastures

This paper aimed at evaluating the effect of summer transhumance to mountain pastures of dairy cows of different breeds on cheese-making ability of milk. Data were from 649 dairy cows of specialized (Holstein Friesian and Brown Swiss) dual purpose (Simmental) and local (mostly Rendena and Alpine Grey) breeds. The Fourier-Transform Infra-Red Spectra (FTIRS) of their milk samples were collected before and after transhumance in 109 permanent dairy farms, and during transhumance in 14 summer farms (with multi-breeds herds) of the Trento Province, north-eastern Italy. A variety of 18 traits describing milk coagulation, curd firming, cheese yield and nutrients recovery in curd/loss in whey were predicted on the basis of FTIRS collected at the individual cow level. Moving the cows to summer farms improved curd firming traits but reduced cheese yields because of an increase of water and fat lost in the whey. During summer grazing, most of cheese-making traits improved, often non-linearly. The milk from summer farms supplementing cows with more concentrates showed better curd firming and cheese yield, because of lower fat lost in the whey. The breed of cows affected almost all the traits with a worst cheese-making ability for milk samples of Holsteins through all the trial, and interacted with concentrate supplementation because increasing compound feed tended to improve cheese-making traits for all breed, with the exception of local breeds for coagulation time and of Brown Swiss for curd firming time. In general, summer transhumance caused a favourable effect on cheese-making aptitude of milk, even though with some difference according to parity, initial days in milk, breed and concentrate supplementation of cows.

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.