Comparison between genetic parameters of cheese yield and nutrient recovery or whey loss traits measured from individual model cheese-making methods or predicted from unprocessed bovine milk samples using Fourier-transform infrared spectroscopy

Cheese yield is an important technological trait in the dairy industry. The aim of this study was to infer the genetic parameters of some cheese yield-related traits predicted using Fourier-transform infrared (FTIR) spectral analysis and compare the results with those obtained using an individual model cheese-producing procedure. A total of 1,264 model cheeses were produced using 1,500-mL milk samples collected from individual Brown Swiss cows, and individual measurements were taken for 10 traits: 3 cheese yield traits (fresh curd, curd total solids, and curd water as a percent of the weight of the processed milk), 4 milk nutrient recovery traits (fat, protein, total solids, and energy of the curd as a percent of the same nutrient in the processed milk), and 3 daily cheese production traits per cow (fresh curd, total solids, and water weight of the curd). Each unprocessed milk sample was analyzed using a MilkoScan FT6000 (Foss, Hillerød, Denmark) over the spectral range, from 5,000 to 900 wavenumber × cm⁻¹. The FTIR spectrum-based prediction models for the previously mentioned traits were developed using modified partial least-square regression. Cross-validation of the whole data set yielded coefficients of determination between the predicted and measured values in cross-validation of 0.65 to 0.95 for all traits, except for the recovery of fat (0.41). A 3-fold external validation was also used, in which the available data were partitioned into 2 subsets: a training set (one-third of the herds) and a testing set (two-thirds). The training set was used to develop calibration equations, whereas the testing subsets were used for external validation of the calibration equations and to estimate the heritabilities and genetic correlations of the measured and FTIR-predicted phenotypes. The coefficients of determination between the predicted and measured values in cross-validation results obtained from the training sets were very similar to those obtained from the whole data set, but the coefficient of determination of validation values for the external validation sets were much lower for all traits (0.30 to 0.73), and particularly for fat recovery (0.05 to 0.18), for the training sets compared with the full data set. For each testing subset, the (co)variance components for the measured and FTIR-predicted phenotypes were estimated using bivariate Bayesian analyses and linear models. The intraherd heritabilities for the predicted traits obtained from our internal cross-validation using the whole data set ranged from 0.085 for daily yield of curd solids to 0.576 for protein recovery, and were similar to those obtained from the measured traits (0.079 to 0.586, respectively). The heritabilities estimated from the testing data set used for external validation were more variable but similar (on average) to the corresponding values obtained from the whole data set. Moreover, the genetic correlations between the predicted and measured traits were high in general (0.791 to 0.996), and they were always higher than the corresponding phenotypic correlations (0.383 to 0.995), especially for the external validation subset. In conclusion, we herein report that application of the cross-validation technique to the whole data set tended to overestimate the predictive ability of FTIR spectra, give more precise phenotypic predictions than the calibrations obtained using smaller data sets, and yield genetic correlations similar to those obtained from the measured traits. Collectively, our findings indicate that FTIR predictions have the potential to be used as indicator traits for the rapid and inexpensive selection of dairy populations for improvement of cheese yield, milk nutrient recovery in curd, and daily cheese production per cow.     

Factors affecting variation of different measures of cheese yield and milk nutrient recovery from an individual model cheese-manufacturing process

Cheese yield (CY) is the most important technological trait of milk, because cheese-making uses a very high proportion of the milk produced worldwide. Few studies have been carried out at the level of individual milk-producing animals due to a scarcity of appropriate procedures for model-cheese production, the complexity of cheese-making, and the frequent use of the fat and protein (or casein) contents of milk as a proxy for cheese yield. Here, we report a high-throughput cheese manufacturing process that mimics all phases of cheese-making, uses 1.5-L samples of milk from individual animals, and allows the simultaneous processing of 15 samples per run. Milk samples were heated (35°C for 40 min), inoculated with starter culture (90 min), mixed with rennet (51.2 international milk-clotting units/L of milk), and recorded for gelation time. Curds were cut twice (10 and 15 min after gelation), separated from the whey, drained (for 30 min), pressed (3 times, 20 min each, with the wheel turned each time), salted in brine (for 60 min), weighed, and sampled. Whey was collected, weighed, and sampled. Milk, curd, and whey samples were analyzed for pH, total solids, fat content, and protein content, and energy content was estimated. Three measures of percentage cheese yield (%CY) were calculated: %CY_CURD, %CY_SOLIDS, and %CY_WATER, representing the ratios between the weight of fresh curd, the total solids of the curd, and the water content of the curd, respectively, and the weight of the milk processed. In addition, 3 measures of daily cheese yield (dCY, kg/d) were defined, considering the daily milk yield. Three measures of nutrient recovery (REC) were computed: REC_FAT, REC_PROTEIN, and REC_SOLIDS, which represented the ratio between the weights of the fat, protein, and total solids in the curd, respectively, and the corresponding components in the milk. Energy recovery, REC_ENERGY, represented the energy content of the cheese compared with that in the milk. This procedure was used to process individual milk samples obtained from 1,167 Brown Swiss cows reared in 85 herds of the province of Trento (Italy). The assessed traits exhibited almost normal distributions, with the exception of REC_FAT. The average values (± SD) were as follows: %CY_CURD = 14.97 ± 1.86, %CY_SOLIDS = 7.18 ± 0.92, %CY_WATER = 7.77 ± 1.27, dCY_CURD = 3.63 ± 1.17, dCY_SOLIDS = 1.74 ± 0.57, dCY_WATER = 1.88 ± 0.63, REC_FAT = 89.79 ± 3.55, REC_PROTEIN = 78.08 ± 2.43, REC_SOLIDS = 51.88 ± 3.52, and REC_ENERGY = 67.19 ± 3.29. All traits were highly influenced by herd-test-date and days in milk of the cow, moderately influenced by parity, and weakly influenced by the utilized vat. Both %CY_CURD and dCY_CURD depended not only on the fat and protein (casein) contents of the milk, but also on their proportions retained in the curd; the water trapped in curd presented an higher variability than that of %CY_SOLIDS. All REC traits were variable and affected by days in milk and parity of the cows. The described model cheese-making procedure and the results obtained provided new insight into the phenotypic variation of cheese yield and recovery traits at the individual level.     

Detection of milk powder and caseinates in Halloumi cheese

Halloumi cheese is traditionally manufactured from fresh milk. Nevertheless, dried dairy ingredients are sometimes illegally added to increase cheese yield. Lysinoalanine and furosine are newly formed molecules generated by heating and drying milk protein components. The levels of these molecular markers in the finished Halloumi have been investigated to verify their suitability to reveal the addition of skim milk powder and calcium caseinate to cheese milk. Because of the severe heating conditions applied in curd cooking, genuine Halloumi cheeses (n = 35), representative of the Cyprus production, were characterized by levels of lysinoalanine (mean value = 8.1 mg/100 g of protein) and furosine (mean value = 123 mg/100 g of protein) unusual for natural cheeses. Despite the variability of the values, a good correlation between the 2 parameters (R = 0.975) has been found in all cheeses, considering both the fresh and mature cheeses as well as those obtained from curd submitted to a prolonged cooking following a traditional practice adopted by a very small number of manufacturers. Experimental cheeses made by adding as low as 5% of skim milk powder, or calcium caseinate, or both, to cheese milk fell outside the prediction limits at ±2 standard deviation of the above-reported correlation regardless of curd cooking conditions or ripening length. This correlation may be adopted as a reliable index of Halloumi cheese genuineness.     

Potential influence of herd and animal factors on the yield of cheese and recovery of components from Sarda sheep milk,as determined by a laboratory bench-top model cheese-making

Individual milk samples from 169 Sarda sheep were collected to characterise the cheese-making potential through the use of a laboratory bench-top model cheese-manufacturing procedure. As the milk samples were not standardised before processing, the data collected at laboratory level fully reflected the great variability of milk from individual animals. The average cheese yield traits of fresh cheese, cheese dry matter and water retained in cheese (as percentages of the milk processed) were 20.6%, 10.1% and 10.6%. The average milk fat and protein recoveries in the curd were 94.0% and 76.7%, respectively. The values for daily production of curd and curd dry matter per sheep were 0.41 kg d⁻¹ and 0.20 kg d⁻¹, respectively. The cheese yield and cheese-related traits were mainly affected by the nutrient content of the milk and the individual effects of the stage of lactation and daily milk yield, respectively, but also by a large individual variation.     

Cheese yield,cheesemaking efficiency,and daily production of 6 breeds of goats

Little is known about the complex process of cheesemaking at the individual level of dairy goats because of the difficulties of producing a high number of model cheeses. The objectives of this work were (1) to study the cheesemaking ability of goat milk; (2) to investigate the variability of cheesemaking-related traits among different farms; (3) to assess the effects of stage of lactation and parity; and (4) to compare 6 breeds of goat (Saanen and Camosciata delle Alpi for the Alpine type; Murciano-Granadina, Maltese, Sarda and Sarda Primitiva for the Mediterranean type) for their cheesemaking ability. For each goat (n = 560) we studied (1) 8 milk quality traits (fat, protein, total solids, casein, lactose, pH, somatic cell score, and bacterial count); (2) 4 milk nutrient recovery traits (fat, protein, total solids, and energy) in curd; (3) 3 actual cheese yield traits (fresh cheese, cheese solids, and cheese water); (4) 2 theoretical cheese yield values (fresh cheese and cheese solids) and the related cheesemaking efficiencies; and (5) daily milk yield and 3 daily cheese yield traits (fresh cheese, cheese solids, and water retained in the curd). With respect to individual animal factors, farm was not particularly important for recovery traits or actual and theoretical cheese yield and estimates of efficiency, whereas it highly influenced daily productions. Parity of goats influenced daily cheese production, whereas DIM slightly affected recovery as well as percent and daily cheese yield traits. Breed was the most important source of variation for almost all cheesemaking traits. Compared with those of Alpine type, the 4 Mediterranean breeds had, on average, lower daily milk and cheese productions, greater actual and theoretical cheese yield, and higher recovery of nutrients in the curd. Among Alpine type, Camosciata delle Alpi was characterized by greater nutrients recovery than Saanen. Within the 4 Mediterranean types, the 3 Italians produced much less milk per day, with much more fat and protein and greater recovery traits than the Murciano-Granadina, resulting in greater actual cheese yield. Within the Italian breeds, milk from Sarda and Sarda Primitiva was characterized by lower daily yields, higher protein and fat content, and greater recoveries of nutrients than Maltese goats. These results confirmed the potential of goat milk for cheese production and could be useful to give new possibilities and direction in breeding programs.     

Effect of β-lactoglobulin A and B whey protein variants on cheese yield potential of a model milk system

Cheese yield mainly depends on the amount and proportion of milk constituents; however, genetic variants of the proteins present in milk may also have an important effect. The objective of this research was to study the effect of the variants A and B of β-lactoglobulin (LG) on cheese yield using a model system consisting of skim milk powder fortified with different levels of a mixture containing α-lactalbumin and β-LG genetic variants (A, B, or A-B) in a 1:2 ratio. Fortified milk samples were subjected to pasteurization at 65°C for 30 min. Miniature cheeses were made by acidifying (pH = 5.9) fortified milk and incubating with rennet for 1 h at 32°C. The clot formed was cut, centrifuged at 2,600 × g for 30 min at 20°C and drained for determining cheese yield. Cheese-yielding capacity was expressed as actual yield (grams of cheese curd per 100 g of milk) and dry weight yield (grams of dried cheese curd per 100 g of milk). Free-zone capillary electrophoresis was used for determining β-LG A or B recovery in the curd during rennet-induced coagulation. The presence of β-LG variant B resulted in a significantly higher actual and dried weight cheese yield than when A or A-B were present at levels ≤0.675% of whey protein (WP) addition. Results of free-zone capillary electrophoresis allowed us to infer that β-LG B associates with the casein micelles during renneting, as shown by an increase in the recovery of this variant in the curd when β-LG B was added up to a maximum at 0.45% (equivalent to 0.675% WP). In general, actual or dried weight cheese yield increased as WP addition was increased from 0.225 to 0.675%. However, when WP addition ranged from 0.675 to 0.90%, a drastic drop in cheese yield was observed. This behavior may be because an increase in the aggregation of casein micelles with a concomitant inclusion of whey protein in the gel occurs at low levels of WP addition, whereas once the association of WP with the casein micelles reach a saturation point at addition levels higher than 0.675%, rearrangements of the gel network result in larger whey expulsion and syneresis. This knowledge is expected to be useful to maximize cheese yield and optimize processing conditions during cheese and cheese analogs manufacturing.     

Genetic parameters of milk coagulation properties and their relationships with milk yield and quality traits in Italian Holstein cows

Milk coagulation properties (MCP) are an important aspect in assessing cheese-making ability. Several studies showed that favorable conditions of milk reactivity with rennet, curd formation rate, and curd strength, as well as curd syneresis, have a positive effect on the entire cheese-making process and subsequently on the ripening of cheese. Moreover, MCP were found to be heritable, but little scientific literature is available about their genetic aspects. The aims of this study were to estimate heritability of MCP and genetic correlations among MCP and milk production and quality traits. A total of 1,071 Italian Holstein cows (progeny of 54 sires) reared in 34 herds in Northern Italy were sampled from January to July 2004. Individual milk samples were collected during the morning milking and analyzed for coagulation time (RCT), curd firmness (a₃₀), pH, titratable acidity, fat, protein, and casein contents, and somatic cell count. About 10% of individual milk samples did not coagulate in 31 min, so they were removed from the analyses. Estimates of heritability for RCT and a₃₀ were 0.25 ± 0.04 and 0.15 ± 0.03, respectively. Estimates of genetic correlations between MCP traits and milk production traits were negligible except for a₃₀ with protein and casein contents (0.44 ± 0.10 and 0.53 ± 0.09, respectively). Estimates of genetic correlations between MCP traits and somatic cell score were strong and favorable, as well as those between MCP and pH and titratable acidity. Selecting for high casein content, milk acidity, and low somatic cell count might be an indirect way to improve MCP without reducing milk yield and quality traits.     

Effects of indirect indicators of udder health on nutrient recovery and cheese yield traits in goat milk

In dairy goats, very little is known about the effect of the 2 most important indirect indicators of udder health [somatic cell count (SCC) and total bacterial count (TBC)] on milk composition and cheese yield, and no information is available regarding the effects of lactose levels, pH, and NaCl content on the recovery of nutrients in the curd, cheese yield traits, and daily cheese yields. Because large differences exist among dairy species, conclusions from the most studied species (i.e., bovine) cannot be drawn for all types of dairy-producing animals. The aims of this study were to quantify, using milk samples from 560 dairy goats, the contemporary effects of a pool of udder health indirect indicators (lactose level, pH, SCC, TBC, and NaCl content) on the recovery of nutrients in the curd (%REC), cheese yield (%CY), and daily cheese yields (dCY). Cheese-making traits were analyzed using a mixed model, with parity, days in milk (DIM), lactose level, pH, SCC, TBC, and NaCl content as fixed effects, and farm, breed, glass tube, and animal as random effects. Results indicated that high levels of milk lactose were associated with reduced total solids recovery in the curd and lower cheese yields, because of the lower milk fat and protein contents in samples rich in lactose. Higher pH correlated with higher recovery of nutrients in the curd and higher cheese yield traits. These results may be explained by the positive correlation between pH and milk fat, protein, and casein in goat milk. High SCC were associated with higher recovery of solids and energy in the curd but lower recovery of protein. The higher cheese yield obtained from milk with high SCC was due to both increased recovery of lactose in the curd and water retention. Bacterial count proved to be the least important factor affecting cheese-making traits, but it decreased daily cheese yields, suggesting that, even if below the legal limits, TBC should be considered in order to monitor flock management and avoid economic losses. The effect of NaCl content on milk composition was linked with lower recovery of all nutrients in the curd during cheese-making. In addition, high milk NaCl content led to reductions in fresh cheese yield and cheese solids. The indirect indicators of the present study significantly affected the cheese-making process. Such information should be considered, to adjust the milk-to-cheese economic value and the milk payment system.     

Proteolysis, lipolysis, volatile compounds, texture, and flavor of Hispánico cheese made using frozen ewe milk curds pressed for different times

Hispánico cheese is manufactured in Spain from a mixture of cow and ewe milk. Production of ewe milk varies throughout the year, with a peak in spring and a valley in summer and autumn. To overcome this seasonal shortage, curd from spring ewe milk may be frozen and used for cheese manufacture some months later. In the present work, ewe milk curds pressed for 15, 60, or 120 min were held at −24°C for 4 mo, thawed, cut to 1-mm pieces, and mixed with fresh cow milk curd for the manufacture of experimental Hispánico cheeses. Control cheese was made from a mixture of pasteurized cow and ewe milk in the same (80:20) proportion. Cheeses, made in duplicate experiments, were analyzed throughout a 60-d ripening period. No significant differences between cheeses were found for lactic acid bacteria counts, dry matter content, hydrophilic peptides, 47 out of 68 vol.tile compounds, texture, and flavor characteristics. On the other hand, differences of minor practical significance between experimental and control cheeses, unrelated to the use of frozen ewe milk curd or the pressing time of ewe milk curd, were found for pH value, aminopeptidase activity, proteolysis, hydrophobic peptides, free amino acids, free fatty acids, and the remaining 21 vol.tile compounds. It may be concluded that the use of frozen ewe milk curd in the manufacture of Hispánico cheese does not alter its main characteristics.     

Preliminary evaluation of endogenous milk fluorophores as tracer molecules for curd syneresis

A front-face fluorescence spectroscopy probe was installed in the wall of a laboratory-scale cheese vat. Excitation and emission filters were chosen for the selective detection of vitamin A, tryptophan, and riboflavin fluorescence. The evolution of the fluorescence of each fluorophore during milk coagulation and syneresis was monitored to determine if they had the potential to act as intrinsic tracers of syneresis and also coagulation. The fluorescence profiles for 2 of the fluorophores during coagulation could be divided into 3 sections relating to enzymatic hydrolysis of κ-casein, aggregation of casein micelles, and crosslinking. A parameter relating to coagulation kinetics was derived from the tryptophan and riboflavin profiles but this was not possible for the vitamin A response. The study also indicated that tryptophan and riboflavin may act as tracer molecules for syneresis, but this was not shown for vitamin A. The evolution of tryptophan and riboflavin fluorescence during syneresis followed a first-order reaction and had strong relationships with curd moisture and whey total solids content (r = 0.86-0.96). Simple 1- and 2-parameter models were developed to predict curd moisture content, curd yield, and whey total solids using parameters derived from the sensor profiles (standard error of prediction = 0.0005-0.394%; R² = 0.963-0.999). The results of this study highlight the potential of tryptophan and riboflavin to act as intrinsic tracer molecules for noninvasive inline monitoring of milk coagulation and curd syneresis. Further work is required to validate these findings under a wider range of processing conditions.     

Genetic parameters of different measures of cheese yield and milk nutrient recovery from an individual model cheese-manufacturing process

Cheese yield (CY) is an important technological trait in the dairy industry, and the objective of this study was to estimate the genetic parameters of cheese yield in a dairy cattle population using an individual model-cheese production procedure. A total of 1,167 Brown Swiss cows belonging to 85 herds were sampled once (a maximum of 15 cows were sampled per herd on a single test day, 1 or 2 herds per week). From each cow, 1,500 mL of milk was processed according to the following steps: milk sampling and heating, culture addition, rennet addition, gelation-time recording, curd cutting, whey draining and sampling, wheel formation, pressing, salting in brine, weighing, and cheese sampling. The compositions of individual milk, whey, and curd samples were determined. Three measures of percentage cheese yield (%CY) were calculated: %CY_CURD, %CY_SOLIDS, and %CY_WATER, which represented the ratios between the weight of fresh curd, the total solids of the curd, and the water content of the curd, respectively, and the weight of the milk processed. In addition, 3 measures of daily cheese yield (dCY, kg/d) were defined, considering the daily milk yield. Three measures of nutrient recovery (REC) were computed: REC_FAT, REC_PROTEIN, and REC_SOLIDS, which represented the ratio between the weights of the fat, protein, and total solids in the curd, respectively, and the corresponding nutrient in the milk. Energy recovery, REC_ENERGY, represented the energy content of the cheese versus that in the milk. For statistical analysis, a Bayesian animal model was implemented via Gibbs sampling. The effects of parity (1 to ≥4), days in milk (6 classes), and laboratory vat (15 vats) were assigned flat priors; those of herd-test-date, animal, and residual were given Gaussian prior distributions. Intra-herd heritability estimates of %CY_CURD, _%CY_SOLIDS, and %CY_WATER ranged from 0.224 to 0.267; these were larger than the estimates obtained for milk yield (0.182) and milk fat content (0.122), and similar to that for protein content (0.275). Daily cheese yields showed heritability estimates similar to those of daily milk yield. The trait %CY_WATER showed a highly positive genetic correlation with %CY_SOLIDS (0.87), whereas their phenotypic correlation was moderate (0.37), and the fat and protein contents of milk showed high genetic correlations with %CY traits. The heritability estimates of REC_PROTEIN and REC_FAT were larger (0.490 and 0.208, respectively) than those obtained for the protein and fat contents of milk, and the genetic relationships between REC_PROTEIN and REC_FAT with milk protein and fat content were low or moderate; REC_PROTEIN and REC_FAT were moderately correlated with the %CY traits and highly correlated with REC_SOLIDS and REC_ENERGY. Both REC_SOLIDS and REC_ENERGY were heritable (0.274 and 0.232), and showed high correlations with each other (0.96) and with the %CY traits (0.83 to 0.97). Together, these findings demonstrate the existence of economically important, genetically determined variability in cheese yield that does not depend solely upon the fat and protein contents of milk, but also relies on the ability of the coagulum to retain the highest possible proportions of the available protein, fat, and water. Exploitation of this interesting genetic variation does not seem to be feasible through direct measurement of the phenotype in cows at the population level. Instead, further research is warranted to examine possible means for indirect prediction, such as through assessing the mid-infrared spectra of milk samples.     

Estimation of genetic parameters for cheese-making traits in Spanish Churra sheep

The global production of sheep milk is growing, and the main industrial use of sheep milk is cheese making. The Spanish Churra sheep breed is one of the most important native dairy breeds in Spain. The present study aimed to estimate genetic parameters for a wide range of traits influencing the cheese-making ability of Churra sheep milk. Using a total of 1,049 Churra ewes, we studied the following cheese-making traits: 4 traits related to milk coagulation properties (rennet coagulation time, curd-firming time, and curd firmness at 30 and 60 min after addition of rennet), 2 traits related to cheese yield (individual laboratory cheese yield and individual laboratory dried curd yield), and 3 traits measuring curd firmness over time (maximum curd firmness, time to attain maximum curd firmness, and syneresis). In addition, a list of milk traits, including the native pH of the milk and several milk production and composition traits (milk yield; the fat, protein, and dried extract percentages; and the somatic cell count), were also analyzed for the studied animals. After discarding the noncoagulating samples (only 3.7%), data of 1,010 ewes were analyzed with multiple-trait animal models by using the restricted maximum likelihood method to estimate (co)variance components, heritabilities, and genetic correlations. In general, the heritability estimates were low to moderate, ranging from 0.08 (for the individual laboratory dried curd yield trait) to 0.42 (for the fat percentage trait). High genetic correlations were found within pairs of related traits (i.e., 0.93 between fat and dried extract percentages, ?0.93 between the log of the curd-firming time and curd firmness at 30 min, 0.70 between individual laboratory cheese yield and individual laboratory dried curd yield, and ?0.94 between time to attain maximum curd firmness and syneresis). Considering all the information provided here, we suggest that in addition to the current consideration of the protein percentage trait for improving cheese yield traits, the inclusion of the pH of milk as a measured trait in the Churra dairy breeding program would represent an efficient strategy for improving the cheese-making ability of milk from this breed.     

Influence of curd cutting programme and stirring speed on the prediction of syneresis indices in cheese-making using NIR light backscatter

An NIR reflectance sensor, with a large field of view and a fibre-optic connection to a spectrometer for measuring light backscatter at 980 nm, was used to monitor the syneresis process online during cheese-making with the goal of predicting syneresis indices (curd moisture content, yield of whey and fat losses to whey) over a range of curd cutting programmes and stirring speeds. A series of trials were carried out in an 11 L cheese vat using recombined whole milk. A factorial experimental design consisting of three curd stirring speeds and three cutting programmes, was undertaken. Milk was coagulated under constant conditions and the casein gel was cut when the elastic modulus reached 35 Pa. Among the syneresis indices investigated, the most accurate and most parsimonious multivariate model developed was for predicting yield of whey involving three terms, namely light backscatter, milk fat content and cutting intensity (R² = 0.83, SEy = 6.13 g/100 g), while the best simple model also predicted this syneresis index using the light backscatter alone (R² = 0.80, SEy = 6.53 g/100 g). In this model the main predictor was the light backscatter response from the NIR light back scatter sensor. The sensor also predicted curd moisture with a similar accuracy.     

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. 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.     

Variations in milk protein fractions affect the efficiency of the cheese-making process

The aim of this study was to assess the influence of the amounts of the α_S1-, α_S2-, β-, and κ-casein (CN) and the α-lactalbumin and β-lactoglobulin protein fractions on the efficiency of the cheese-making process independently of their genetic polymorphisms. The study was carried out on milk samples from 1,271 Brown Swiss cows from 85 herds classified into 4 categories according to management, feeding, and housing characteristics (traditional and modern systems). To assess the efficiency of the cheese-making process, we processed the milk samples according to a laboratory cheese-making procedure (1,500 mL/sample) and obtained the following measures: (1) 3 percentage cheese yields (%CY_curd, %CY_solids, %CY_water), (2) 2 daily cheese yields obtained by multiplying %CY (curd and total solids) by daily milk yields (dCY_curd, dCY_solids), (3) 4 measures of nutrient recovery in the curd (REC_fat, REC_protein, REC_solids, REC_energy), and (4) 2 measures of cheese-making efficiency in terms of the ratio between the observed and theoretical %CY (Ef-%CY_curd, Ef-%CY_solids). All the aforementioned traits were analyzed by fitting 2 linear mixed models with protein fractions as fixed effects expressed as percentage in the milk (model M-%milk) and as percentage of the total casein content (model M-%cas) together with the effects of total casein content (only in model M-%cas), daily milk yield (only in model M-%milk; not for dCY traits), dairy system, herd (random effect), days in milk, parity, and vat. The efficiency of overall cheese yield (Ef-%CY_curd) was mostly positively associated with β-CN content in the milk, whereas Ef-%CY_solids was greater with higher amounts of κ-CN and α_S1-CN (M-%milk) due to the strong influence of both fractions on the recovery rate of milk components in the curd (fat and total solids, protein with α_S1-CN only) when expressed as percentage of milk and of total casein; only β-CN was more important for REC_protein. In contrast, we found β-lactoglobulin to be highly negatively related to all the traits related to the cheese-making process and to the daily cheese yield per cow, whereas α-lactalbumin was positively associated with the latter traits. Additional research on this topic is needed, with particular focus on the genetic and genomic aspects of the role of protein fractions in the cheese-making process and on the associations between genetic polymorphisms in milk protein and milk nutrient recovery in the curd.     

Validation of a curd-syneresis sensor over a range of milk composition and process parameters

An online visible-near-infrared sensor was used to monitor the course of syneresis during cheesemaking with the purpose of validating syneresis indices obtained using partial least squares, with cross-validation across a range of milk fat levels, gel firmness levels at cutting, curd cutting programs, stirring speeds, milk protein levels, and fat:protein ratio levels. Three series of trials were carried out in an 11-L cheese vat using recombined whole milk. Three factorial experimental designs were used, consisting of 1) 3 curd stirring speeds and 3 cutting programs; 2) 3 milk fat levels and 3 gel firmness levels at cutting; and 3) 2 milk protein levels and 3 fat:protein ratio levels, respectively. Milk was clotted under constant conditions in all experiments and the gel was cut according to the respective experimental design. Prediction models for production of whey and whey fat losses were developed in 2 of the experiments and validated in the other experiment. The best models gave standard error of prediction values of 6.6 g/100 g for yield of whey and 0.05 g/100 g for fat in whey, as compared with 4.4 and 0.013 g/100 g, respectively, for the calibration data sets. Robust models developed for predicting yield of whey and whey fat losses using a validation method have potential application in the cheese industry.     

Computer vision and color measurement techniques for inline monitoring of cheese curd syneresis

Optical characteristics of stirred curd were simultaneously monitored during syneresis in a 10-L cheese vat using computer vision and colorimetric measurements. Curd syneresis kinetic conditions were varied using 2 levels of milk pH (6.0 and 6.5) and 2 agitation speeds (12.1 and 27.2 rpm). Measured optical parameters were compared with gravimetric measurements of syneresis, taken simultaneously. The results showed that computer vision and colorimeter measurements have potential for monitoring syneresis. The 2 different phases, curd and whey, were distinguished by means of color differences. As syneresis progressed, the backscattered light became increasingly yellow in hue for circa 20 min for the higher stirring speed and circa 30 min for the lower stirring speed. Syneresis-related gravimetric measurements of importance to cheese making (e.g., curd moisture content, total solids in whey, and yield of whey) correlated significantly with computer vision and colorimetric measurements.     

Pathway-based genome-wide association analysis of milk coagulation properties,curd firmness,cheese yield,and curd nutrient recovery in dairy cattle

It is becoming common to complement genome-wide association studies (GWAS) with gene-set enrichment analysis to deepen the understanding of the biological pathways affecting quantitative traits. Our objective was to conduct a gene ontology and pathway-based analysis to identify possible biological mechanisms involved in the regulation of bovine milk technological traits: coagulation properties, curd firmness modeling, individual cheese yield (CY), and milk nutrient recovery into the curd (REC) or whey loss traits. Results from 2 previous GWAS studies using 1,011 cows genotyped for 50k single nucleotide polymorphisms were used. Overall, the phenotypes analyzed consisted of 3 traditional milk coagulation property measures [RCT: rennet coagulation time defined as the time (min) from addition of enzyme to the beginning of coagulation; k₂₀: the interval (min) from RCT to the time at which a curd firmness of 20 mm is attained; a₃₀: a measure of the extent of curd firmness (mm) 30 min after coagulant addition], 6 curd firmness modeling traits [RCT_eq: RCT estimated through the CF equation (min); CF_P: potential asymptotic curd firmness (mm); k_CF: curd-firming rate constant (% × min^?1); k_SR: syneresis rate constant (% × min^?1); CF_max: maximum curd firmness (mm); and t_max: time to CF_max (min)], 3 individual CY-related traits expressing the weight of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed and 4 milk nutrient and energy recoveries in the curd (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), milk pH, and protein percentage. Each trait was analyzed separately. In total, 13,269 annotated genes were used in the analysis. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases were queried for enrichment analyses. Overall, 21 Gene Ontology and 17 Kyoto Encyclopedia of Genes and Genomes categories were significantly associated (false discovery rate at 5%) with 7 traits (RCT, RCT_eq, k_CF, %CY_SOLIDS, REC_FAT, REC_SOLIDS, and REC_ENERGY), with some being in common between traits. The significantly enriched categories included calcium signaling pathway, salivary secretion, metabolic pathways, carbohydrate digestion and absorption, the tight junction and the phosphatidylinositol pathways, as well as pathways related to the bovine mammary gland health status, and contained a total of 150 genes spanning all chromosomes but 9, 20, and 27. This study provided new insights into the regulation of bovine milk coagulation and cheese ability that were not captured by the GWAS.     

Udder measurements and milk production in two Awassi sheep genotypes and their crosses

To improve the effectiveness of community-based breeding programs for increased milk production, the values of different udder measurements for predicting milk production traits during the milking period were assessed over 3 yr on 273 Awassi ewes. Machine milking of ewes began after weaning, 56 d after parturition, and continued until the milk yield of the ewes was <200 mL/d. Milk yield obtained by hand milking and milk composition were measured weekly, and days in milk, total milk yield, and total yields of protein, fat, and nonfat solids in milk were calculated for each ewe. On d 70 of milking, morphological traits of the whole udder (circumference, width, height, and length), udder cistern (height), and teats (length, width, and position score) were measured. On the same day, the milk yield of ewes was recorded by hand milking. Positive and moderate to strong correlations (r = 0.36 to 0.76) between udder circumference and width, teat width, and milk production traits of total milk yield, and total yields of protein, fat, and nonfat solids were found. However, a more accurate predictor of milk production traits was milk yield on d 70, as higher positive correlations between this variable and the milk production traits were found (r = 0.63 to 0.89). Nine farmers were invited to independently estimate the hand-milked milk yield performance of a sample of 169 ewes (d 15 to 45 of milking) by visually observing each ewe and making a subjective linear score (1 to 5). Their assessments were significantly correlated with milk yield on the day of the observation (r = 0.52), total milk yield (r = 0.50), and days in milk (r = 0.45). Considering the perception details provided by farmers concerning each of the subjective linear scores, it was found that most predictive linear udder measurements of udder circumference and width and teat width identified in this study were implicit in these scores. The predictive ability of the measurements studied have practical implications for community-based breeding programs involving improvement of milk production—not just in Syria, but in other countries in dry areas as well—because it is possible for experienced farmers to visually assess milk production of dairy ewes or take simple udder measurements with predictive value.     

Short communication: Occurrence of aflatoxin M1 in the Manchego cheese supply chain

The importance of ewe milk lies in the production of high quality cheeses, such as Manchego cheese with a Protected Designation of Origin, whose safety must be guaranteed. In a 2-yr study, 407 bulk tank milk samples from farms and 82 silo milk and curd samples from cheese factories were collected from southeast Spain and tested for aflatoxin M1 (AFM1) using 2 commercial ELISA tests. Of these, 99.3% of the bulk tank samples had AFM1 levels below the European Union (EU) legal limit for milk (50 ng/kg), and well below the limit adopted by the Codex Alimentarius (500 ng/kg). Moreover, 98.8% of the silo milk and curd samples from cheese factories had AFM1 levels below the EU limit for milk. When considering median AFM1 concentrations, an average 4-fold increase was found in the final curd in relation to the corresponding silo milk. Control of AFM1 in Manchega ewe milk would enhance dairy product safety by the possible detection of faults in the manufacture of Manchego cheese.