Cheesemaking aptitude of two Spanish dairy ewe breeds: Changes during lactation and relationship between physico-chemical and technological properties

Acid and rennet gelation properties of sheep, goat, and cow milks: Effects of processing and seasonal variation

Gelation is an important functional property of milk that enables the manufacture of various dairy products. This study investigated the acid (with glucono-δ-lactone) and rennet gelation properties of differently processed sheep, goat, and cow milks using small-amplitude oscillatory rheological tests. The impacts of ruminant species, milk processing (homogenization and heat treatments), seasonality, and their interactions were studied. Acid gelation properties were improved (higher gelation pH, shorter gelation time, and higher storage modulus (G') by intense heat treatment (95°C for 5 min) to comparable extents for sheep and cow milks, both better than those for goat milk. Goat milk produced weak acid gels with low G' (<100 Pa) despite improvements induced by heat treatments. Seasonality had a marked impact on the acid gelation properties of sheep milk. The acid gels of late-season sheep milk had a lower gelation pH, no maximum in tan δ following gel formation, and 70% lower G' values than those from other seasons. We propose the potential key role of a critical acid gelation pH that induces structural rearrangements in determining the viscoelastic properties of the final gels. For rennet-induced gelation, compared with cow milk, the processing treatments of the goat and sheep milks had much smaller impacts on their gelation properties. Intense heat treatment (95°C for 5 min) prolonged the rennet gelation time of homogenized cow milk by 8.6 min (74% increase) and reduced the G' of the rennet gels by 81 Pa (85% decrease). For sheep and goat milks, the same treatment altered the rennet gelation time by only less than 3 min and the G' of the rennet gels by less than 14 Pa. This difference may have been caused by the different physicochemical properties of the milks, such as differences in their colloidal stability, proportion of serum-phase caseins, and ionic calcium concentration. The seasonal variations in the gelation properties (both acid and rennet induced) of goat milk could be explained by the minor variation in its protein and fat contents. This study provides new perspectives and understandings of milk gelation by demonstrating the interactive effects among ruminant species, processing, and seasonality.

Genome-wide association studies (GWAS) and post-GWAS analyses for technological traits in Assaf and Churra dairy breeds

This study aimed to perform a GWAS to identify genomic regions associated with milk and cheese-making traits in Assaf and Churra dairy sheep breeds; second, it aimed to identify possible positional and functional candidate genes and their interactions through post-GWAS studies. For 2,020 dairy ewes from 2 breeds (1,039 Spanish Assaf and 981 Churra), milk samples were collected and analyzed to determine 6 milk production and composition traits and 6 traits related to milk coagulation properties and cheese yield. The genetic profiles of the ewes were obtained using a genotyping chip array that included 50,934 SNP markers. For both milk and cheese-making traits, separate single-breed GWAS were performed using GCTA software. The set of positional candidate genes identified via GWAS was subjected to guilt-by-association-based prioritization analysis with ToppGene software. Totals of 84 and 139 chromosome-wise significant associations for the 6 milk traits and the 6 cheese-making traits were identified in this study. No significant SNPs were found in common between the 2 studied breeds, possibly due to their genetic heterogeneity of the phenotypes under study. Additionally, 63 and 176 positional candidate genes were located in the genomic intervals defined as confidence regions in relation to the significant SNPs identified for the analyzed traits for Assaf and Churra breeds. After the functional prioritization analysis, 71 genes were identified as promising positional and functional candidate genes and proposed as targets of future research to identify putative causative variants in relation to the traits under examination. In addition, this multitrait study allowed us to identify variants that have a pleiotropic effect on both milk production and cheese-related traits. The incorporation of variants among the proposed functional and positional candidate genes into genomic selection strategies represent an interesting approach for achieving rapid genetic gains, specifically for those traits difficult to measure, such as cheese-making traits.

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 Composition, Antioxidant Status, Fatty Acid Profile, and Volatile Compounds of Milk and Fresh and Ripened Ewes' Cheese from a Sustainable Part-Time Grazing System

We conducted the first nutritional analysis of dairy products from the traditional Roja Mallorquina sheep breed. Samples of bulk raw milk were taken twice a month from December 2015 to March 2016 from sheep fed using a part-time grazing system, and fresh soft (FC, n = 8) and ripened (RC, n = 8) cheeses were made. The variability in vitamins, total phenolic compounds (TPC), total antioxidant capacity (TAC), and fatty acid (FA) content was influenced by the cheese-making process (differences between the cheese and the original milk) and by the type of cheese-making technology (mainly related to heating, the use of starter culture, and ripening). The most notable physicochemical characteristic of the cheeses was low fat content (24.1 and 29.6 g/100 g for FC and RC). Milk and RC were characterised by major concentrations of retinol (211.4 and 233.6 μg/100 g dry matter (DM), respectively) and TPC (18.7 and 54.6 μg/100 g DM, respectively), while FC was characterised by major concentrations of retinol (376.4 μg) and α-tocopherol (361.7 μg). The fat-soluble components of the FC generally exhibited better nutritional value for human health than those of the milk and RC, with a higher level of retinol and α-tocopherol; lower values for saturated FA, atherogenic, and thrombogenic indices; and higher levels of monounsaturated FA, polyunsaturated FA, n-3, and n-6. Acids, alcohols, and ketones comprised almost 95% of the volatile compounds detected. Acetoin and products of lactose and citrate metabolism played an important role in the development of the aromatic attributes of both kinds of cheese. This preliminary study can contribute to add value to these traditional products according to healthy nutritional criteria and supports the implementation of strategies to promote their commercialisation and obtain product labelling as “pasture-fed” or specific marks.

Gene Networks Driving Genetic Variation in Milk and Cheese-Making Traits of Spanish Assaf Sheep

Most of the milk produced by sheep is used for the production of high-quality cheese. Consequently, traits related to milk coagulation properties and cheese yield are economically important to the Spanish dairy industry. The present study aims to identify candidate genes and their regulators related to 14 milk and cheese-making traits and to develop a low-density panel of markers that could be used to predict an individual's genetic potential for cheese-making efficiency. In this study, we performed a combination of the classical genome-wide association study (GWAS) with a stepwise regression method and a pleiotropy analysis to determine the best combination of the variants located within the confidence intervals of the potential candidate genes that may explain the greatest genetic variance for milk and cheese-making traits. Two gene networks related to milk and cheese-making traits were created using the genomic relationship matrices built through a stepwise multiple regression approach. Several co-associated genes in these networks are involved in biological processes previously found to be associated with milk synthesis and cheese-making efficiency. The methodology applied in this study enabled the selection of a co-association network comprised of 374 variants located in the surrounding of genes showing a potential influence on milk synthesis and cheese-making efficiency.

Goat cheese yield and recovery of fat, protein, and total solids in curd are affected by milk coagulation properties

The aims of the present research were to quantify the effects of each coagulation trait, traditional milk coagulation properties [MCP: rennet coagulation time (RCT), curd-firming time (k₂₀), and curd firmness at 30 min (a₃₀)], and modeled curd-firming over time (CF_t) parameters [estimated rennet coagulation time (RCT_eq), curd-firming instant rate constant (k_CF), and potential curd firmness (CF_P)] directly on the following: (1) recovery of 3 milk components in the curd (%REC), (2) 3 measures of cheese yield (%CY), and (3) 3 daily cheese yield traits (dCY) from goat milk. Cheese-making traits were analyzed using 2 mixed different models, the first to test MCP and the second to test CF_t parameters. Pearson correlations were also calculated. Significant and favorable relationships (negative for time intervals and positive for CF measures) were found between the traditional MCP and the CF_t parameters and %REC and %CY traits. The effects of milk fat and protein contents were particularly important on all cheese-making traits, with the only exception being the effect of fat content on water retention in cheese (%CY_WATER). We found an optimum value of milk k₂₀, associated with the highest recovery of components and cheese yield in solids (%CY_SOLIDS). In addition, a lower level of curd water retention and an increased fresh curd yield (%CY_CURD) were associated with greater recovery of fat. The collection of all available information during the process of milk coagulation and curd-firming allowed us to discover the effect of RCT_eq on %REC traits and %CY_SOLIDS, which had not previously been revealed for traditional RCT. Moreover, higher k_CF values were associated with increased %CY_CURD and %CY_SOLIDS. Given that CF_t parameters showed a high level of independence from one another, these can also be easily used and characterized in future applications at the industry level. Information provided by traditional and modeled coagulation properties could efficiently support the goat dairy industry and lay the foundations for a quality payment scheme for goat milk.

Estimates of heritability and genetic correlations for milk coagulation properties and individual laboratory cheese yield in Sarda ewes

Objective of this study was to estimate genetic parameters of milk coagulation properties (MCPs) and individual laboratory cheese yield (ILCY) in a sample of 1018 Sarda breed ewes farmed in 47 flocks. Rennet coagulation time (RCT), curd-firming time (k 20) and curd firmness (a 30) were measured using Formagraph instrument, whereas ILCY were determined by a micromanufacturing protocol. About 10% of the milk samples did not coagulate within 30 min and 13% had zero value for k 20. The average ILCY was 36%. (Co)variance components of considered traits were estimated by fitting both single- and multiple-trait animal models. Flock-test date explained from 13% to 28% of the phenotypic variance for MCPs and 26% for ILCY, respectively. The largest value of heritability was estimated for RCT (0.23±0.10), whereas it was about 0.15 for the other traits. Negative genetic correlations between RCT and a 30 (-0.80±0.12), a 30 and k 20 (-0.91±0.09), and a 30 and ILCY (-0.67±0.08) were observed. Interesting genetic correlations between MCPs and milk composition (r G>0.40) were estimated for pH, NaCl and casein. Results of the present study suggest to use only one out of three MCPs to measure milk renneting ability, due to high genetic correlations among them. Moreover, negative correlations between ILCY and MCPs suggest that great care should be taken when using these methods to estimate cheese yield from small milk samples.

Quality traits and modeling of coagulation, curd firming, and syneresis of sheep milk of Alpine breeds fed diets supplemented with rumen-protected conjugated fatty acid

The aim of this study was to test the modeling of curd-firming (CF) measures and to compare the sheep milk of 3 Alpine breeds supplemented with or without rumen-protected conjugated linoleic acid (rpCLA). Twenty-four ewes of the Brogna, Foza, and Lamon breeds were allotted to 6 pens (2 pens/breed) and fed a diet composed of corn grain, corn silage, dried sugar beet pulp, soybean meal, wheat bran, wheat straw, and a vitamin-mineral mixture. The rpCLA supplement (12 g/d per ewe plus 4 g/d for each lamb older than 30 d) was mixed into the diet of 1 pen per sheep breed (3 pens/treatment) to provide an average of 0.945 and 0.915 g/d per ewe of the cis-9,trans-11 C18:2 and trans-10,cis-12 C18:2 conjugated linoleic acid isomers, respectively. The trial started at 38 ± 23 d after parturition, and individual morning milk samples were collected on d 16, 23, 37, 44, and 59 of the trial. Milk samples were analyzed for composition, and duplicate samples were assessed for milk coagulation properties (MCP). A total of 180 CF measures for each sample (1 every 15s) were recorded. Model parameters were the rennet coagulation time, the asymptotic potential CF, the CF instant rate constant, the syneresis instant rate constant, the maximum CF achieved within 45 min (CFmax), and the time at achievement of CFmax. The data were analyzed using a hierarchical model that considered the fixed effects of breed, diet, lamb birth, and initial days in milk, which were tested on individual ewe (random) variance; the fixed effect of sampling day, which was tested on the within-ewe sample (random) variance; and the fixed effect of instrument or cuvette position (only for MCP), which was tested on the residual (replicates within samples) variance. The local Alpine sheep breeds displayed similar milk compositions, traditional MCP, and CF modeling parameters. Supplementation with rpCLA triggered changes in milk composition and worsened MCP (e.g., delayed rennet coagulation time, slower CF instant rate constant, and a doubling of syneresis instant rate constant), but did not influence potential CF. Overall, our results indicate that rpCLA supplementation reduced the actual maximum CF (CFmax) but did not modify the interval between rennet addition and CFmax or time to CFmax.