Article

# Associations between Milk Coagulation Properties and Microbiological Quality in Sheep Bulk Tank Milk

Ramón Arias<fnref n="1" />,*, Lorena Jiménez <fnref n="1" />,*, Ana Garzón <fnref n="2" />, Javier Caballero-Villalobos <fnref n="2" />, Bonastre Oliete <fnref n="3" />, Nicolò Amalfitano <fnref n="4" />, Alessio Cecchinato <fnref n="4" /> and José M. Perea <fnref n="2" />

**Abstract:** This study conducted a seasonal analysis of bulk tank milk from 77 sheep farms to establish relationships between the concentration of major microbial groups and milk coagulation properties. The investigated milk traits included composition (pH, fat, casein, lactose), coagulation properties (curd firmness: A<sub>60</sub>-, rennet clotting time: RCT-, curd firming time: k<sub>20</sub>-, curd yield: CY-), and somatic cell score (SCS). The main microbial groups analyzed were total mesophilic bacteria (SPC), thermodurics (THERMO), psychrotrophs (PSYCHRO), *Pseudomonas* spp. (PSEUDO), lactic acid bacteria (LAB), catalase-negative gram-positive cocci (GPCNC), *Escherichia coli* (ECOLI), coliforms other than *Escherichia coli* (COLI), coagulase-positive staphylococci (CPS), coagulase-negative staphylococci (CNS), and spores of lactate-fermenting *Clostridium* (BAB). Mixed linear models were used to explore associations between coagulation properties and the aforementioned variables. Results demonstrated that incorporating microbial loads into the models improves their fit and the relative quality of the outcomes. An important seasonality is demonstrated by an increase in CY and A<sub>60</sub>, along with a decrease in RCT and k<sub>20</sub> during autumn and winter, contrasting with spring and summer. BAB concentration resulted in a reduction of A<sub>60</sub> and an increase in RCT, whereas SPC concentration led to an enhancement of A<sub>60</sub> and a reduction in RCT. An increase in GPCNC concentration was associated with an increase in k<sub>20</sub> and a decrease in CY.

**Keywords:** coagulation traits; dairy industry; microbiological quality; milk sheep quality; tank milk

## 1. Introduction

Sheep milk production holds significant importance in countries across the Mediterranean basin. Traditionally, native breeds have been raised, with many now facing the threat of extinction. These breeds have played a crucial role in sustaining livestock activity and associated industries in vast rural areas. This not only helps to maintain local populations but also contributes to environmental protection [1]. Sheep milk production is primarily directed towards cheese and dairy product processing, with many of these products holding protected designation of origin (PDO) status [2]. Many of these cheeses are known for their unique sensory characteristics, which result from traditional manufacturing techniques. The quality of sheep milk is a key factor in creating exceptional cheeses and is subject to rigorous quality control systems in the industry [3]. Traditionally, milk quality assessment has focused on physicochemical aspects, mainly fat and protein content. Recently, hygienic parameters have been integrated into these quality systems, including

<figure><img src="image_1.png" /></figure>

check for
updates

**Citation:** Arias, R.; Jiménez, L.; Garzón, A.; Caballero-Villalobos, J.; Oliete, B.; Amalfitano, N.; Cecchinato, A.; Perea, J.M. Associations between Milk Coagulation Properties and Microbiological Quality in Sheep Bulk Tank Milk. *Foods* **2024**, *13*, 886.
https://doi.org/10.3390/foods13060886

Received: 17 February 2024
Revised: 7 March 2024
Accepted: 12 March 2024
Published: 14 March 2024

<figure><img src="image_2.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Instituto Regional de Investigación y Desarrollo Agroalimentario y Forestal de Castilla-La Mancha (IRIAF), CERSYRA de Valdepeñas, 13300 Ciudad Real, Spain; rarias@jccm.es</footnote>

<footnote n="2">² Departamento de Producción Animal, Universidad de Córdoba, 14071 Córdoba, Spain; pa1gasia@uco.es (A.G.); javier.caballero@uco.es (J.C.-V.); jmperea@uco.es (J.M.P.)</footnote>

<footnote n="3">³ Institut Agro, Université Bourgogne Franche-Comté, PAM UMR A 02.102, 21000 Dijon, France; bonastre.oliete@u-bourgogne.fr</footnote>

<footnote n="4">⁴ Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, Legnaro, 35020 Padova, Italy; nicolo.amalfitano@unipd.it (N.A.); alessio.cecchinato@unipd.it (A.C.)</footnote>

<footnote>* Correspondence: ljimenez@quesomanchego.es</footnote>
analyses like total mesophilic microorganism and somatic cell counts (SCC). However, practices considering aspects such as differential milk microbiology are still limited [3].

Research in dairy cattle has provided valuable insights into the microbial composition of milk and its relationship with several sources of contamination, emphasizing the importance of this aspect in milk quality [4]. These contaminants introduce a variety of microorganisms associated with farm environments, such as lactic acid bacteria, hygiene practices on farms leading to the presence of coliforms, and intramammary infections caused by staphylococci [5]. Other researchers [6,7] have also highlighted the complex microbial composition found in milk from different species, emphasizing the importance of its contamination sources and its impact on cheese-making processes, especially concerning lactic acid bacteria. Some studies have even suggested that stringent hygiene measures have significantly reduced microbial populations, potentially adversely affecting the production of traditional and artisanal cheeses [8]. However, there is limited research specifically addressing the differential microbial load and microbiota of sheep milk [9,10].

Several studies have underscored the significance of milk coagulation properties (MCP) as complementary traits alongside the inherent quality of raw milk [11,12]. For that reason, there is well-documented literature on the relevance and significance of MCP, mainly in cattle [13,14] and, to a lesser extent, in sheep and goats [15,16]. Also, these studies underline the importance of improving milk MCP and explore the relationship with several factors that contribute its modification, such as breeding, feeding, farm management system, season, or milking system. Factors such as bacterial counts and the intricate interplay among microbiological characteristics and technological parameters can influence MCP and cheese production. However, research in this area, particularly regarding its impact on milk from small ruminants, remains limited [5].

Therefore, this study aims to assess the impact of composition and microorganisms present in sheep milk on its coagulation properties. Additionally, it seeks to compare the predictive performance of models incorporating microbiological variables with reference models that do not include these variables. The ultimate goal is to reach a better understanding of the microbiological factors influencing milk coagulation, providing valuable information for quality control and dairy production management.

## 2. Materials and Methods

### 2.1. *Study Design and Sample Determination*

The present study included 77 Manchega dairy sheep flocks from the region of Castilla-La Mancha (Spain), registered under Protected Designation of Origin "Queso Manchego", representing 15% of the farms listed in PDO [17]. The average flock size was 1144 ewes, ranging from 150 to 5500 ewes, with an average production of 143 kg/ewe/year, varying from 91 to 200 kg, over a mean lactation period of 121 days, ranging from 100 to 150 days. All farms employ mechanical milking systems, and ewes are milked twice daily. The milk is then filtered and stored in refrigerated tanks, maintaining a maximum temperature of 6 °C. The production system is semi-intensive, traditionally associated with grazing on natural pastures and utilizing the residues and remains of cereal crops. On average, each sheep has access to 1.44 ha, ranging from 0.5 to 1.94 ha.

The flocks were visited four times a year (once during each season) to collect bulk tank samples, which included milk from the morning milking of the visit day and the previous evening milking. Samples were stored in hermetically sealed sterilized containers and kept at 4 °C until the arrival at the Experimental Dairy Science and Cheese Laboratory at the Regional Centre for Animal Selection and Reproduction (CERSYRA-IRIAF, Valdepeñas, Ciudad Real, Spain), where they were aliquoted and sent to other facilities for further analyses.

### 2.2. *Laboratory Analysis*

Microbiological analyses were conducted at CERSYRA-IRIAF within the first 24 h following sample collection. The techniques and methods used to determine the micro-
biological counts are described in detail in a previous paper from our research group [5]. Serial dilutions were performed from each sample to inoculate 0.1 mL onto different culture media. Bacterial counts for the following groups of microorganisms were determined on PCA media (Panreac, Barcelona, Barcelon, Spain): standard plate count (SPC) was incubated in aerobic conditions at 30 °C for 72 h [18]; thermoduric bacteria (THERMO, Panreac, Barcelona, Spain) were incubated under the same conditions as SPC, after pasteurizing milk at 62.8 °C for 30 min [19]; and psychrotrophic bacteria (PSYCHRO, Panreac, Barcelona, Spain) were incubated at 6.5 °C for 10 days [20]. *Pseudomonas* spp. (PSEUDO, Panreac, Barcelona, Spain) was cultured on Cetrimide agar (Panreac, Spain) and incubated at 35 °C for 48 h [21]. Determination of *Escherichia coli* (ECOLI, Basingstoke Hampshire, UK) and other coliforms (COLI, Basingstoke Hampshire, UK) was conducted using CromoIDTM Coli medium (bioMérieux, Madrid, Spain), incubated at 37 °C for 24 h [22]. Gram-positive catalase-negative cocci count (GPCNC) was determined in modified Edwards medium with colistin and oxolinic acid supplement (Oxoid, Basingstoke Hampshire, UK), incubated at 35 °C for 48 h [9]. Lactic acid bacteria (LAB) were seeded on MRS medium (Panreac, Barcelona, Spain) acidified to pH 5.7, and incubated at 30 °C for 72 h [23]. Lactate-fermenting Clostridium spore (BAB, Basingstoke Hampshire, UK) count was performed using the most probable number (MPN) technique in Bryant and Burkey Broth (BBB, Merck, Darmstadt, Germany) [24]. For the enumeration of coagulase-positive staphylococci (CPS) and coagulase-negative staphylococci (CNS), Baird-Parker RPF agar medium (bioMérieux, Madrid, Spain) was used, and incubated at 37 °C for 24 h [4]. The microbial counts were expressed in colony-forming units per mL of milk (CFU/mL), and Clostridium spores were expressed as spores/mL. A logarithmic transformation (log₁₀) was applied to both microbial counts to normalize their distribution.

Milk composition, pH, and somatic cell count (SCC) analyses were conducted at the Interprofessional Dairy Laboratory of Castilla-La Mancha (LILCAM, Talavera de la Reina, Toledo, Spain) within 48 h from sample collection. The pH was measured using a Crison Basic 20 pH meter (Crison Instruments, Barcelona, Spain). Milk composition, including fat content (FAT), total casein (CAS), and lactose (LAC) was determined using a Milkscan 6000 FT device (Foss Electric, Hillerød, Denmark), which operates on mid-infrared spectroscopy. The total SCC was quantified by flow cytometry on a Fossomatic FC (Foss Electric, Hillerød, Denmark). The SCC was subsequently expressed as somatic cell scores (SCS) by applying the following formula [25]:

$$SCS = \log_2[\text{SCC cells/mL} \times 10^{-5}] + 3$$

Milk coagulation properties were monitored at the Small Ruminant Dairy Laboratory (University of Córdoba, Spain) within the first 24 h following sample collection. Using a Formagraph viscosimeter (Foss Electric A/S, Hillerød, Denmark), key coagulation parameters were determined based on the methodology adapted from McMahon and Brown [26]. These parameters included rennet coagulation time (RCT), curd-firming time (k₂₀), curd firmness at 60 min after rennet addition (A₆₀), and curd yield (CY). Curd yield was subsequently calculated following draining of the fresh curds [5].

## 2.3. Statistical Analysis

Data analysis was performed on SAS version 9.3 (SAS Institute Inc., Cary, NC, USA). Preliminary testing of the data was carried out to determine outliers to be discarded before further analysis. After the editing, a total of 306 data points were used (Table 1). The explanatory variables (i.e., SCS, BAB, SPC, THERMO, PSYCHRO, PSEUDO, LAB, GPCNC, ECOLI, COLI, CPS, CNS) were discretized into 5 classes based on the 20th, 40th, 60th, and 80th percentiles. This approach was taken to ensure a better assessment of the pattern of microbiological count effects and to guarantee the classes were balanced [27].
<table><caption>**Table 1.** Descriptive statistics of the variables for composition, microbiological quality, and coagulation properties of sheep milk from bulk tank (n = 308).</caption><thead><tr><th>Variable</th><th>Description</th><th>Mean</th><th>SD <sup>1</sup></th><th>Q<sub>1</sub></th><th>Q<sub>3</sub></th></tr></thead><tbody><tr><td>CY (g/100 mL) <sup>2</sup></td><td>Curd yield</td><td>3.13</td><td>0.55</td><td>2.76</td><td>3.74</td></tr><tr><td>A<sub>60</sub> (mm) <sup>3</sup></td><td>Curd firmness at 60 min.</td><td>33.04</td><td>10.76</td><td>26.30</td><td>39.98</td></tr><tr><td>RCT (min) <sup>4</sup></td><td>Rennet clotting time</td><td>31.54</td><td>11.97</td><td>22.30</td><td>39.45</td></tr><tr><td>k<sub>20</sub> (min) <sup>5</sup></td><td>Curd firming time</td><td>8.79</td><td>14.17</td><td>3.0</td><td>7.0</td></tr><tr><td>SCS <sup>6</sup></td><td>Somatic cell score</td><td>6.40</td><td>0.78</td><td>5.91</td><td>6.90</td></tr><tr><td>FAT (%)</td><td>Fat content</td><td>7.87</td><td>0.83</td><td>7.32</td><td>8.35</td></tr><tr><td>CAS (%)</td><td>Casein content</td><td>4.82</td><td>0.41</td><td>4.55</td><td>5.04</td></tr><tr><td>LAC (%)</td><td>Lactose content</td><td>4.60</td><td>0.24</td><td>4.48</td><td>4.79</td></tr><tr><td>pH (−log [H+])</td><td>pH</td><td>6.61</td><td>0.15</td><td>6.56</td><td>6.69</td></tr><tr><td>BAB <sup>7</sup></td><td>Lactate-fermenting Clostridium spores</td><td>3.42</td><td>0.63</td><td>2.96</td><td>3.80</td></tr><tr><td>SPC <sup>8</sup></td><td>Total mesophilic bacteria (standard plate count)</td><td>5.65</td><td>0.77</td><td>5.11</td><td>6.13</td></tr><tr><td>THERMO <sup>8</sup></td><td>Thermodurics</td><td>3.25</td><td>0.84</td><td>2.69</td><td>3.73</td></tr><tr><td>PSYCHRO <sup>8</sup></td><td>Psychrotrophs</td><td>5.23</td><td>1.08</td><td>4.57</td><td>5.90</td></tr><tr><td>PSEUDO <sup>8</sup></td><td>Pseudomonas spp.</td><td>3.15</td><td>0.84</td><td>2.65</td><td>3.72</td></tr><tr><td>LAB <sup>8</sup></td><td>Lactic acid bacteria</td><td>4.89</td><td>0.73</td><td>4.30</td><td>5.40</td></tr><tr><td>GPCNC <sup>8</sup></td><td>Gram-positive catalase-negative cocci</td><td>4.16</td><td>1.63</td><td>3.55</td><td>4.78</td></tr><tr><td>ECOLI <sup>8</sup></td><td>Escherichia coli</td><td>1.63</td><td>1.16</td><td>0.13</td><td>2.31</td></tr><tr><td>COLI <sup>8</sup></td><td>Coliforms other than Escherichia coli</td><td>2.91</td><td>1.19</td><td>2.35</td><td>3.67</td></tr><tr><td>CPS <sup>8</sup></td><td>Coagulase-positive staphylococci</td><td>2.52</td><td>1.35</td><td>1.90</td><td>3.46</td></tr><tr><td>CNS <sup>8</sup></td><td>Coagulase-negative staphylococci</td><td>4.30</td><td>0.54</td><td>4.06</td><td>4.50</td></tr></tbody></table>

A linear mixed model (MIXED) procedure was used to investigate the associations between microbiological quality and milk coagulation properties. The model is represented as follows:

$$Y = \sum_{i=1}^{n} \beta_{i}f_{i}+\delta$$

where Y denotes milk coagulation traits (CY, A₆₀, RCT, k₂₀), βᵢ are the unknown parameters to be estimated, fᵢ denotes the explanatory variables, and δ is the error term.

Two models were developed for each dependent variable. Initially, a baseline model (BM) was established, excluding microbiological counts, with the purpose of comparison against a second model (whole set model: WSM) that included all relevant variables. The BM model was adjusted, considering flock as a random factor (77 classes), while season (4 classes: spring, summer, autumn, winter) and SCS (5 classes: 4.44–5.72; 5.72–6.23, 6.23–6.65, 6.65–7.01, 7.01–8.90 log₂ SCC) were set up as fixed effects. Additionally, covariates such as pH, FAT (%), CAS (%), and LAC (%) were incorporated.

The WSM model incorporating microbiological counts was adjusted using the same variables as BM, completing it with microbiological counts discretized into five classes as fixed factors (Table 2).

To address challenges arising from the limited sample size (308 complete cases), a step-by-step sequence was undertaken, considering the primary effects of all fixed, random, and covariate factors. A significance level of *p* < 0.05 was used for the retention criterion. In the initial phase, all conceivable single-variable models were assessed using the Akaike Information Criterion (AIC) value. Subsequently, each remaining factor was incrementally added, one by one, to the model with the lowest AIC value, and each predictor was compared based on the AIC value. This iterative process continued until the model with the lowest AIC was identified [28]. This model was deemed the most plausible and was selected as the final model.

<footnote n="1">¹ SD: standard deviation, Q₁: 1st quartile; Q₃: 3rd quartile. ² CY: curd yield. ³ A₆₀: curd firmness at 60 min.</footnote>

<footnote n="4">⁴ RCT: rennet clotting time. ⁵ k₂₀: curd firming time. ⁶ SCS: bulk tank milk somatic cells score (log₂[SCC cells/mL × 10⁻⁵] + 3). ⁷ BAB: lactate-fermenting Clostridium spores (log₁₀[spores/mL]). ⁸ SPC: total mesophilic bacteria (standard plate count), THERMO: thermodurics, PSYCHRO: psychrotrophs, PSEUDO: Pseudomonas spp., LAB: lactic acid bacteria, GPCNC: gram-positive catalase-negative cocci, ECOLI: Escherichia coli, COLI: coliforms other than Escherichia coli, CPS: coagulase-positive staphylococci, CNS: coagulase-negative staphylococci. All microbial counts are expressed as: log₁₀(UFC/mL).</footnote>
<table><caption>**Table 2.** Classes of microbial counts used in linear mixed models as fixed factors (n = 308).</caption><thead><tr><th rowspan="2">Variable</th><th colspan="5">Classes</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr></thead><tbody><tr><td>BAB<sup>1</sup></td><td>4.04–4.85</td><td>4.85–5.26</td><td>5.26–5.38</td><td>5.38–5.89</td><td>5.89–7.00</td></tr><tr><td>SPC<sup>2</sup></td><td>2.48–3.04</td><td>3.04–3.30</td><td>3.30–3.63</td><td>3.63–4.15</td><td>4.15–5.04</td></tr><tr><td>THERMO<sup>2</sup></td><td>3.99–5.01</td><td>5.01–5.41</td><td>5.41–5.70</td><td>5.70–6.20</td><td>6.20–8.10</td></tr><tr><td>PSYCHRO<sup>2</sup></td><td>0.00–2.60</td><td>2.60–3.03</td><td>3.03–3.30</td><td>3.30–3.89</td><td>3.89–5.48</td></tr><tr><td>PSEUDO<sup>2</sup></td><td>2.00–4.43</td><td>4.43–4.98</td><td>4.98–5.36</td><td>5.36–6.14</td><td>6.14–8.00</td></tr><tr><td>LAB<sup>2</sup></td><td>0.00–2.55</td><td>2.55–2.95</td><td>2.95–3.34</td><td>3.34–3.80</td><td>3.80–4.53</td></tr><tr><td>GPCNC<sup>2</sup></td><td>2.91–4.22</td><td>4.22–4.66</td><td>4.66–5.11</td><td>5.11–5.46</td><td>5.46–6.48</td></tr><tr><td>ECOLI<sup>2</sup></td><td>0.00–2.48</td><td>2.48–3.39</td><td>3.39–4.37</td><td>4.37–4.97</td><td>4.97–5.48</td></tr><tr><td>COLI<sup>2</sup></td><td>0.00–1.00</td><td>1.00–1.65</td><td>1.65–1.96</td><td>1.96–2.48</td><td>2.48–5.48</td></tr><tr><td>CPS<sup>2</sup></td><td>0.00–2.45</td><td>2.45–2.89</td><td>2.89–3.31</td><td>3.31–3.98</td><td>3.98–5.48</td></tr><tr><td>CNS<sup>2</sup></td><td>0.00–1.65</td><td>1.65–2.50</td><td>2.50–3.06</td><td>3.06–3.59</td><td>3.59–5.48</td></tr></tbody></table>

After defining the model, an analysis of collinearity was conducted using the variance inflation factor (VIF) of the regression coefficients. A severe multicollinearity issue was present if any of the VIFs exceeded 10 [28]. The normal distribution of residuals was verified using Kolmogorov-Smirnov test, while the absence of autocorrelation of residuals was assessed through the Durbin-Watson test. Heteroscedasticity was evaluated using White's test. The model fit was assessed using the adjusted coefficient of determination (R²) and the root mean squared error (RMSE) [29].

Polynomial contrasts ($p < 0.05$) were computed to delineate the pattern of effects of microbiological counts (BAB, SPC, THERMO, PSYCHRO, PSEUDO, LAB, GPCNC, ECOLI, COLI, CPS, CNS). First-order comparisons measure linear relationships, while second-order comparisons measure quadratic relationships.

### 3. Results and Discussion

The most plausible models identified for each of the analyzed coagulation characteristics are shown in Table 3. Broadly, the models exhibited poor fit (R² < 50%), except for CY. In the case of curd yield, the coefficient of determination ranged between 54.8% for the baseline model and 57.9% for the model incorporating microbiological counts. This emphasizes the fact that the inclusion of microbiological counts not only significantly improved the model fit but also enhanced relative quality across all evaluated outcomes.

A significant association was observed between season and all coagulation traits. This association persisted even in models that incorporated microbiological counts, albeit with a slightly reduced contribution. As expected, season was a significant source of variation for all investigated traits, underscoring the importance of considering this effect in the statistical model for a better evaluation of the effect of microbiological counts. It is important to note that Manchega does not exhibit as pronounced a seasonality in milk production as other Mediterranean breeds, such as Sarda [30] or Latxa [31]. Therefore, the seasonal differences found are arguably more interesting in terms of their relevance to the cheese industry.

The impact of flock on the variability of the examined coagulation traits was modest, as indicated by previous research [11,32]. Although the AIC favored models excluding the flock effect, this does not completely rule out its influence. These results are explained by the fact that the evaluated flocks follow a similar farming system, with other variables within the scope of this study exerting a more pronounced influence on coagulation characteristics than flock variability.

<footnote n="1">¹ BAB: lactate-fermenting Clostridium spores (log₁₀[spores/mL]). ² SPC: total mesophilic bacteria (standard plate count), THERMO: thermodurics, PSYCHRO: psychrotrophs, PSEUDO: *Pseudomonas* spp., LAB: lactic acid bacteria, GPCNC: gram-positive catalase-negative cocci, ECOLI: *Escherichia coli*, COLI: coliforms other than *Escherichia coli*, CPS: coagulase-positive staphylococci, CNS: coagulase-negative staphylococci. All microbial counts are expressed as: log₁₀(UFC/mL).</footnote>
<table><caption>**Table 3.** More plausible linear mixed models (F-value and significance) based on the composition of the milk (BM), and the incorporation of microorganism counts (WSM) (n = 308).</caption><thead><tr><th rowspan="2">Variable <sup>1</sup></th><th colspan="2">CY (g/100 mL)</th><th colspan="2">A<sub>60</sub> (mm)</th><th colspan="2">RCT (min)</th><th colspan="2">k<sub>20</sub> (min)</th></tr><tr><th>BM</th><th>WSM</th><th>BM</th><th>WSM</th><th>BM</th><th>WSM</th><th>BM</th><th>WSM</th></tr></thead><tbody><tr><td>Flock</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td></tr><tr><td>Season</td><td>31.27 ***</td><td>26.47 ***</td><td>4.35 ***</td><td>2.55 *</td><td>51.99 ***</td><td>35.56 ***</td><td>5.87 ***</td><td>4.72 **</td></tr><tr><td>SCS <sup>2</sup></td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>2.27 *</td><td>ns</td><td>ns</td><td>ns</td></tr><tr><td>FAT (%)</td><td>21.71 ***</td><td>13.30 ***</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td></tr><tr><td>CAS (%)</td><td>26.07 ***</td><td>14.71 ***</td><td>7.60 ***</td><td>6.98 ***</td><td>ns</td><td>ns</td><td>5.77 **</td><td>5.10 *</td></tr><tr><td>LAC (%)</td><td>ns</td><td>ns</td><td>ns</td><td>3.75 *</td><td>70.75 ***</td><td>12.81 ***</td><td>ns</td><td>ns</td></tr><tr><td>pH (−log[H+])</td><td>2.76 *</td><td>ns</td><td>ns</td><td>ns</td><td>ns</td><td>55.69 ***</td><td>8.51 ***</td><td>5.68 *</td></tr><tr><td>BAB <sup>3</sup></td><td>-</td><td>ns</td><td>-</td><td>6.20 ***</td><td>-</td><td>4.44 ***</td><td>-</td><td>ns</td></tr><tr><td>SPC <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>3.04 **</td><td>-</td><td>6.25 ***</td><td>-</td><td>ns</td></tr><tr><td>THERMO <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>PSYCHRO <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>PSEUDO <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>LAB <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>GPCNC <sup>4</sup></td><td>-</td><td>3.28 **</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>2.97 *</td></tr><tr><td>ECOLI <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>COLI <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>CPS <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>CNS <sup>4</sup></td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td><td>-</td><td>ns</td></tr><tr><td>AIC <sup>5</sup></td><td>-592.06</td><td>-445.62</td><td>1409.53</td><td>1116.48</td><td>1356.69</td><td>1067.44</td><td>754.21</td><td>634.08</td></tr><tr><td>RMSE <sup>6</sup></td><td>0.37</td><td>0.33</td><td>10.31</td><td>9.76</td><td>9.38</td><td>8.82</td><td>3.77</td><td>3.03</td></tr><tr><td>R<sup>2</sup> <sup>7</sup></td><td>54.80</td><td>57.9</td><td>8.50</td><td>13.80</td><td>39.1</td><td>44.44</td><td>10.40</td><td>13.10</td></tr></tbody></table>

The level of somatic cells was only significantly associated with RCT. This result is in line with a previous study by Pazzola et al. [33] focused on the Sarda breed, in which higher somatic cell counts were linked to slower coagulations. Moreover, Caballero-Villalobos et al. [34] observed in individual samples from Manchega that an increase in SCC led to a rise in RCT and k<sub>20</sub> and a decrease in curd firmness measured at 30 min. This could be attributed to an increase in casein breakdown as a consequence of increased plasmin activity, the main enzyme associated with somatic cells in milk [35].

Regarding CY, the most plausible baseline model was adjusted considering season, CAS, FAT, and pH, in that order of importance. The potential inclusion of microbiological counts led to a better model, which specifically incorporated the level of gram-positive catalase-negative cocci (GPCNC), thus excluding pH. As expected, CY was positively associated with FAT and CAS, which is consistent with the correlations found for these quality parameters and CY in Manchega sheep milk [5]. In other breeds, such as Sarda, Assaf, or Churra, this relationship has also been documented [30,36,37]. CY was also negatively associated with the level of GPCNC, with the linear pattern of the relationship being confirmed across the classes of GPCNC (Figure 1).

It has been described that high microbial counts could be associated with an enzymatic load acting proteolytically and lipolytically on the milk, leading to a decrease in its fat and protein content, indirectly affecting coagulation properties [38]. In terms of season, CY was higher in autumn and winter milk than in spring and summer. The higher yields obtained in these seasons could be related to the increased levels of the major components of the milk, as has also been shown in various studies [39,40].

<footnote n="1">¹ CY: curd yield, A<sub>60</sub>: curd firmness at 60 min., RCT: rennet clotting time, k<sub>20</sub>: curd firming time. BM: baseline model, WM: whole set of variables model. ² SCS: bulk tank milk somatic cells score (log<sub>2</sub>[SCC cells/mL × 10<sup>-5</sup>] + 3).</footnote>

<footnote n="3">³ BAB: lactate-fermenting *Clostridium* spores (log<sub>10</sub>[spores/mL]). ⁴ SPC: total mesophilic bacteria (standard plate count), THERMO: thermodurics, PSYCHRO: psychrotrophs, PSEUDO: *Pseudomonas* spp., LAB: lactic acid bacteria, GPCNC: gram-positive catalase-negative cocci, ECOLI: *Escherichia coli*, COLI: coliforms other than *Escherichia coli*, CPS: coagulase-positive staphylococci, CNS: coagulase-negative staphylococci. All microbial counts are expressed as: log<sub>10</sub>(UFC/mL).</footnote>

<footnote n="5">⁵ AIC: Akaike information criterion. ⁶ RMSE: root mean squared error. ⁷ R<sup>2</sup>: Determination coefficient. *: *p* < 0.05, **: *p* < 0.01, ***: *p* < 0.001, ns: *p* > 0.05.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Least squares means and 95% CI of curd yield (CY) across significant variables (I: FAT (%); II: CAS (%); III: GPCNC, gram-positive catalase-negative cocci (log₁₀[CFU/mL]); IV: season). When the trend line is presented, it means that the linear contrast (---) was significant at *p* < 0.05. Means with different superscript letters are different at *p* < 0.05 (SNK test).</figcaption></figure>

Regarding A60, the best baseline model was adjusted with season and CAS, the latter being the most important factor. However, the model fit was the poorest of all, with a coefficient of determination of 8.5%. The possibility of including microbiological counts led to improving the model, incorporating BAB and SPC, and also including LAC. The level of BAB acquired a similar importance to CAS for A60. The association of A60 was positive with CAS and SPC and negative with LAC and BAB (Figure 2). SPC and BAB levels showed a good fit with linear and quadratic patterns. The associations found are consistent with the correlations for these quality parameters and A60 in Manchega [5] and Sarda [11,12] sheep milk, which indicate that the decrease in curd firmness is inversely proportional to the increase in lactose content, with the relationship between firmer curds and the level of CAS being noteworthy. An opposite effect of the vegetative forms of total mesophilic microorganisms (SPC) and spore-forming forms in milk (BAB) on curd hardness was observed, for which no previous references have been found, and which will be appropriate to address in future studies. A60 follows a seasonal pattern similar to CY, with increased curd firmness in autumn and winter. Todaro et al. [41] found greater hardness in spring and autumn, possibly related to the production system itself, which conditions milk quality.

The most plausible baseline model for RCT was adjusted with LAC, season, and somatic cell level, achieving a coefficient of determination of 39.1%. LAC and season were the most relevant factors for the onset of coagulation. In contrast, the model that incorporated microbiological levels was more plausible, obtaining a coefficient of determination of 44.4%. This model included the levels of SPC and BAB and also incorporated pH, this time excluding somatic cell levels. The association of RCT was positive with pH, LAC, and BAB and negative with the SPC. These results are in line with those reported by other authors [11,12], who reported an increase in coagulation times when individual sheep milk samples have a more alkaline pH, although this same study shows a negative relationship between lactose content and RCT. The mechanisms that exactly explain the association between lactose and
technological traits are still unknown, although it could be associated with a modification of the percentage and composition of minerals and proteins, in accordance with the role of lactose as an osmotic regulator of milk [11,12]. BAB followed a linear pattern, while linear and quadratic patterns were confirmed for SPC (Figure 3). It is confirmed that different types of microorganisms do not have the same effect on RCT, as was also observed for curd hardness; on one hand, active vegetative cells, such as SPC, consume lactose to assure growth and, consequently, decrease RCT, while spore-forming microorganisms like BAB are associated with a longer duration of coagulation, which would have certain negative connotations in optimizing the cheese-making process. Regarding season, RCT is significantly higher in spring and lower in autumn, similar to results obtained in Latxa sheep, with increased RCT from April to July [31], and in dairy cattle, with the lower value in autumn [42]. Conversely, other studies reported the lowest coagulation times in spring and the highest in winter for bulk tank sheep milk [41]. In any case, in our study, differences in coagulation time depending on season are much higher than those of the cited studies, which, as previously mentioned, may be related to the production system of this particular breed.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Least squares means and 95% CI of maximum curd firmness (A<sub>60</sub>) across significant variables (I: CAS (%); II: LAC (%); III: BAB, butyric acid bacteria (log<sub>10</sub>[spores/mL]); IV: SPC, total mesophilic bacteria (standard plate count) (log<sub>10</sub>[CFU/mL]); V: season). When the trend line is presented, it means that the linear (---) and/or quadratic (----) contrast was significant at *p* < 0.05. Means with different superscript letters are different at *p* < 0.05 (SNK test).</figcaption></figure>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Figure 3.** Least squares means and 95% CI of rennet clotting time (RCT) across significant variables (I: pH; II: LAC (%); III: BAB, butyric acid bacteria (log₁₀[spores/mL]); IV: SPC, total mesophilic bacteria (standard plate count) (log₁₀[CFU/mL]); V: Season). When the trend line is presented, it means that the linear (---) and/or quadratic (---) contrast was significant at *p* < 0.05. Means with different superscript letters are different at *p* < 0.05 (SNK test).

The best baseline model for k20 was adjusted with season, CAS, and pH, achieving a
poor coefficient of determination of 10.4%. The model with microbiological counts added
the GPCNC level and increased the coefficient of determination to 13.1%. The association
of k20 was positively associated with pH, similar to what was observed for RCT, and
negatively associated with CAS, similar to GPCNC, which showed a negative linear pattern
(Figure 4). The relationship between k20 and CAS follows the trend already highlighted
in previous studies, particularly when coagulation times are high (RCT > 30 min) [5].
According to some authors, high microorganism counts, as in this case, GPCNC, affect
coagulation by causing acidification of the milk and a faster firming of the curd [38]. Similar
to RCT, k20 was significantly higher in spring and lower in autumn. Others, on the contrary,
reported lower levels for k20 in spring and winter [41], differences that could be due, as
discussed for other parameters, to either the production system of Valle del Belice sheep,
which is more extensive and dependent on the availability of pastures, or to differences in
the stage of lactation [43,44].

<figure><img src="image_10.png" /><figcaption>**Figure 4.** Least squares means and 95% CI of curd-firming time ($k_{20}$) across significant variables (I: pH; II: CAS (%); III: GPCNC, gram-positive catalase-negative cocci ($\log_{10}$[CFU/mL]); IV: season). When the trend line is presented, it means that the linear contrast (— —) was significant at *p* < 0.05. Means with different superscript letters are different at *p* < 0.05 (SNK test).</figcaption></figure>

**4. Conclusions**

The applied models demonstrated limited fitting for most milk coagulation proper-
ties, except for curd yield (CY). Incorporating microbiological analysis into the models
significantly enhanced both the fit and the relative quality of the results, highlighting the
importance of microbiological factors in dairy industry quality control systems. Significant
associations were found between seasonal variations and all coagulation properties, even
when microbiological counts were considered. During autumn and winter, the coagulation
process initiates earlier (RCT) and proceeds at a slower pace ($k_{20}$), resulting in firmer curds
($A_{60}$) and yielding better curd production (CY). Both lactate-fermenting Clostridium spores
(BAB) and total mesophilic bacteria (SPC) showed significant associations with $A_{60}$ and
RCT, albeit with opposing effects. While BAB concentration reduced curd firmness and
increased coagulation onset, SPC concentration enhanced curd firmness and promoted
coagulation onset. Coagulation speed and curd yield were significantly associated with
gram-positive catalase-negative cocci (GPCNC). As GPCNC concentration increased, coagu-
lation speed accelerated, and curd yield decreased. These results emphasize the complexity
of interactions among microbiological factors and milk composition from farms and their
collective impact on coagulation characteristics and seasonal variations, a crucial aspect for
optimizing cheese production.

**Author Contributions:** Conceptualization, R.A. and J.M.P.; methodology, A.G. and J.M.P.; software, J.M.P.; validation, L.J. and J.C.-V.; formal analysis, L.J., J.M.P., B.O., N.A. and A.C.; investigation, R.A., L.J., B.O., J.C.-V. and A.G.; resources, R.A.; data curation, J.M.P., R.A. and L.J.; writing—original draft preparation, R.A.; writing—review and editing, R.A., L.J., J.C.-V., N.A. and A.C.; supervision, R.A.
and J.M.P.; project administration, R.A.; funding acquisition, R.A. and J.M.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** The present study was funded by the National Institute of Research and Agrarian and Food Technology—INIA (Government of Spain) as part of the project RTA2011-00057-C02-01.

**Institutional Review Board Statement:** Ethical approval was not required for this project, since all samples were bulk tank milk.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

**Acknowledgments:** The authors wish to acknowledge the National Manchega Breeders Association (AGRAMA) and the Regulatory Council of the PDO *Queso Manchego* for their valuable collaboration in this work.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Scintu, M.F.; Piredda, G. Typicity and biodiversity of goat and sheep milk products. *Small Rum. Res.* **2007**, *68*, 221–231. [CrossRef]

2. Regulation (EU) No 1151/2012 of the European Parliament and of the Council of 21 November 2012 on quality schemes for agricultural products and foodstuffs. *Off. J. Eur. Union* **2012**, *343*, 1–29.

3. Pirisi, A.; Lauret, A.; Dubeuf, J.P. Basic and incentive payments for goat and sheep milk in relation to quality. *Small Rum. Res.* **2007**, *68*, 167–178. [CrossRef]

4. Jayarao, B.M.; Pillai, S.R.; Sawant, A.A.; Wolfgang, D.R.; Hegde, N.V. Guidelines for monitoring bulk tank milk somatic cell and bacterial counts. *J. Dairy Sci.* **2004**, *8*, 3561–3573. [CrossRef]

5. Jiménez, L.; Caballero-Villalobos, J.; Garzón, A.; Oliete, B.; Pérez-Guzmán, M.D.; Arias, R. Exploring the relationships between coagulation, composition, and hygienic quality of bulk tank milk from Manchega sheep. *Small Rum. Res.* **2023**, *228*, 107106. [CrossRef]

6. Quigley, L.; O'Sullivan, O.; Stanton, C.; Beresford, T.P.; Ross, R.P.; Fitzgerald, G.F.; Cotter, P.D. The complex microbiota of raw milk. *FEMS Microbiol. Rev.* **2013**, *37*, 664–698. [CrossRef]

7. Toquet, M.; Gómez-Martín, A.; Bataller, E. Review of the bacterial composition of healthy milk, mastitis milk and colostrum in small ruminants. *Res. Vet. Sci.* **2023**, *140*, 1–5. [CrossRef] [PubMed]

8. Monsallier, F.; Verdier-Metz, I.; Agabriel, C.; Martin, B.; Montel, M.-C. Variability of microbial teat skin flora in relation to farming practices and individual dairy cow characteristics. *Dairy Sci. Technol.* **2012**, *92*, 265–278. [CrossRef]

9. De Garnica, M.L.; Linage, B.; Carriedo, J.A.; de La Fuente, L.F.; García-Jimeno, M.C.; Santos, J.A.; Gonzalo, C. Relationship among specific bacterial counts and total bacterial and somatic cell counts and factors influencing their variation in ovine bulk tank milk. *J. Dairy Sci.* **2013**, *96*, 1021–1029. [CrossRef] [PubMed]

10. Esteban-Blanco, C.; Gutierrez-Gil, B.; Marina, H.; Pelayo, R.; Suarez-Vega, A.; Acedo, A.; Arranz, J. The milk microbiota of the Spanish Churra sheep breed: New insights into the complexity of the Milk microbiome of dairy species. *Animals* **2020**, *10*, 1463. [CrossRef]

11. Pazzola, M.; Dettori, M.L.; Cipolat-Gotet, C.; Cecchinato, A.; Bittante, G.; Vacca, G.M. Phenotypic factors affecting coagulation properties of milk from Sarda ewes. *J. Dairy Sci.* **2014**, *97*, 7247–7257. [CrossRef] [PubMed]

12. Paschino, P.; Vacca, G.M.; Dettori, M.L.; Pazzola, M. An approach for the estimation of somatic cells' effect in Sarda sheep milk based on the analysis of milk traits and coagulation properties. *Small Rum. Res.* **2019**, *171*, 77–81. [CrossRef]

13. Nilsson, K.; Stålhammar, H.; Stenholdt Hansen, M.; Lindmark-Månsson, H.; Duchemin, S.; Fikse, F.; de Koning, D.J.; Paulsson, M.; Glantz, M. Characterisation of non-coagulating milk and effects of milk composition and physical properties on rennet-induced coagulation in Swedish Red Dairy Cattle. *Int. Dairy J.* **2019**, *95*, 50–57. [CrossRef]

14. Franzoi, M.; Costa, A.; Vigolo, V.; Penasa, M.; De Marchi, M. Effect of pasteurization on coagulation properties of bovine milk and the role of major composition traits and protein fractions. *J. Food Compos. Anal.* **2022**, *114*, 104808. [CrossRef]

15. Pazzola, M. Coagulation Traits of Sheep and Goat Milk. *Animals* **2019**, *9*, 540. [CrossRef]

16. Vacca, G.M.; Stocco, G.; Dettori, M.L.; Bittante, G.; Pazzola, M. Goat cheese yield and recovery of fat, protein, and total solids in curd are affected by milk coagulation properties. *J. Dairy Sci.* **2020**, *103*, 1352–1365. [CrossRef]

17. PDO Queso Manchego. *Protected Designation of Origin "Manchego Cheese": Activities Report of the PDO "Manchego Cheese 2023"*; Unpublished Data; Fundación Consejo Regulador de la D.O. Queso Manchego: Valdepeñas, Spain, 2023.

18. ISO 4833-1:2013; Microbiology of the Food Chain—Horizontal Method for the Enumeration of Microorganisms—Part 1: Colony Count at 30 °C by the Pour Plate Technique. International Organization for Standardization: Geneva, Switzerland, 2013.
19. ISO 4833-2:2013; Microbiology of the Food Chain—Horizontal Method for the Enumeration of Microorganisms—Part 2: Colony Count at 30 °C by the Surface Plating Technique. International Organization for Standardization: Geneva, Switzerland, 2013.

20. ISO 17410:2019; Microbiology of the Food Chain—Horizontal Method for the Enumeration of Psychrotrophic Microorganisms. International Organization for Standardization: Geneva, Switzerland, 2019.

21. Dore, S.; Fadda, A.; Fresi, S.; Denti, G.; Puggioni, G.; Cannas, E.A. Pseudomonas aeruginosa: Persistence in the milking machine: Healing experiences. In *Special Issue of the International Dairy Federation, Proceedings of the IDF 5th International Symposium on the Challenge to Sheep and Goats Milk Sectors, Alghero, Italy, 18–20 April 2007*; International Dairy Federation: Schaerbeek, Belgium, 2008; pp. 79–80.

22. ISO 16649-2:2001; Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Beta-Glucuronidase-Positive Escherichia coli—Part 2: Colony-Count Technique at 44 Degrees C Using 5-bromo-4-chloro-3-indolyl Beta-Dglucuronide. International Organization for Standardization: Geneva, Switzerland, 2001.

23. ISO 15214:1998; Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Mesophilic Lactic Acid Bacteria—Colony-Count Technique at 30 Degrees C. International Organization for Standardization: Geneva, Switzerland, 1998.

24. United States Department of Agriculture. *Most Probable Number Tables, MLG Appendix 2.02*; United States Department of Agriculture: Washington, DC, USA, 2003; pp. 1–8.

25. Ali, A.K.A.; Shook, G.E. An optimum transformation for somatic cells concentration in milk. *J. Dairy Sci.* **1980**, *63*, 487-490.
[CrossRef]

26. McMahon, D.J.; Brown, R.J. Evaluation of formagraph for comparing rennet solutions. *J. Dairy Sci.* **1982**, *65*, 1639-1642. [CrossRef]

27. Mariani, E.; Cipolat-Gotet, C.; Stefanon, B.; Zeconyi, A.; Stocco, G.; Sandri, M.; Ablondi, M.; Mountricha, M.; Summer, A. Effect of total and differential somatic cell count on yield, composition and predicted coagulation properties from individual dairy cows. *Int. Dairy J.* **2022**, *75*, 298-307. [CrossRef]

28. Quintana, A.R.; Perea, J.M.; García-Béjar, B.; Jiménez, L.; Garzón, A.; Arias, R. Dominant Yeast Community in Raw Sheep's Milk and Potential Transfers of Yeast Species in Relation to Farming Practices. *Animals* **2020**, *10*, 906. [CrossRef]

29. Alin, A. Wiley Interdisciplinary Reviews: Computational Statistics. *Multicollinearity* **2010**, *2*, 370-374.

30. Manca, M.G.; Serdino, J.; Gaspa, G.; Urgeghe, P.; Ibba, I.; Contu, M.; Fresi, P.; Macciotta, N.P.P. Derivation of multivariate indices of milk composition, coagulation properties, and individual cheese yield in dairy sheep. *J. Dairy Sci.* **2016**, *99*, 4547-4557. [CrossRef]

31. Abilleira, E.; Virto, M.; Nájera, A.I.; Salmerón, J.; Albisu, M.; Pérez-Elortondo, F.J.; Ruiz de Gordoa, J.C.; de Renobales, M.; Barron, L.J.R. Effects of seasonal changes in feeding management under part-time grazing on the evolution of the composition and coagulation properties of raw milk from ewes. *J. Dairy Sci.* **2010**, *93*, 3902-3909. [CrossRef]

32. Cecchinato, A.; Cipolat-Gotet, C.; Casellas, J.; Penasa, M.; Rossoni, A.; Bittante, G. Genetic analysis of rennet coagulation time, curd-firming rate, and curd firmness assessed over an extended testing period using mechanical and near-infrared instruments. *J. Dairy Sci.* **2013**, *96*, 50-62. [CrossRef]

33. Pazzola, M.; Cipolat-Gotet, C.; Bittante, G.; Cecchinato, A.; Dettori, M.L.; Vacca, G.M. Phenotypic and genetic relationships between indicators of the mammary gland health status and milk composition, coagulation, and curd firming in dairy sheep. *J. Dairy Sci.* **2018**, *101*, 1-12. [CrossRef] [PubMed]

34. Caballero-Villalobos, J.; Arias, R.; McSweeney, P.L.H.; Garzón, A. The effect of plasmin activity over rennet coagulation in ewe milk. In *Proceedings of the IDF International Symposium on Cheese Science and Technology & the IDF Symposium on Concentration and Drying Technologies of Dairy Products, Dublin, Ireland, 11–13 April 2016*; p. 120.

35. Pinto, G.; Cairá, S.; Nicolai, M.A.; Mauriello, R.; Cuollo, M.; Pirisi, A.; Piredda, G.; Chianese, L.; Addeo, F. Proteolysis and partial dephosphorylation of casein are affected by high somatic cell counts in sheep milk. *Food Res. Int.* **2013**, *53*, 510-521. [CrossRef]

36. Milan, M.J.; Caja, G.; Gonzalez-Gonzalez, R.; Fernandez-Perez, A.M.; Such, X. Structure and performance of Awassi and Assaf dairy sheep farms in northwestern Spain. *J. Dairy Sci.* **2011**, *94*, 771-784. [CrossRef] [PubMed]

37. Othmane, M.H.; Carriedo, J.A.; De La Fuente, L.F.; San Primitivo, F. Factors affecting test-day milk composition in dairy ewes, and relationships amongst various milk components. *J. Dairy Res.* **2022**, *69*, 53-62. [CrossRef] [PubMed]

38. Sevi, A.; Albenzio, M.; Taibi, L.; Dantone, D.; Massa, S.; Annicchiarico, G. Changes of somatic cell count through lactation and their effects on nutritional, renneting and bacteriological characteristics of ewe's milk. *Adv. Food Sci.* **1999**, *21*, 122-127.

39. Banks, J.M. Cheese yield. In *Cheese Problems Solved*; McSweeney, P.L.H., Ed.; Woodhead Publishing Series in Food Science, Technology and Nutrition; Woodhead Publishing: Sawston, UK, 2007; pp. 100-114.

40. Jaramillo, D.P.; Zamora, A.; Guamis, B.; Rodriguez, M.; Trujillo, A.J. Cheesemaking aptitude of two Spanish dairy ewe breeds: Changes during lactation and relationship between physicochemical and technological properties. *Small Rum. Res.* **2008**, *78*, 48-55. [CrossRef]

41. Todaro, M.; Scatassa, M.L.; Bonanno, A. The quality assessment of Valle del Belice sheep milk and cheeses produced in the hot summer season in Sicily. *Dairy Sci. Technol.* **2014**, *94*, 225-239. [CrossRef]

42. De Marchi, M.; Dal Zotto, R.; Cassandro, M.; Bittante, G. Milk Coagulation Ability of Five Dairy Cattle Breeds. *J. Dairy Sci.* **2007**, *90*, 3986-3992. [CrossRef]
43. Pellegrini, O.; Remeuf, F.; Rivemale, M. Evolution des caracteristiques physico-chimiques et des parametres de coagulation du lait de brebis collecte dans la region de Roquefort. *Lait* **1994**, *74*, 425-442. [CrossRef]

44. Nájera, A.I.; Barron, L.J.R.; Ribeiro, P.; Pelissier, F.; Abilleira, E.; Perez-Elortondo, F.J.; Albisu, M.; Salmeron, J.; Ruiz de Gordoa, J.C.; Virto, M.; et al. Seasonal changes in the technological and compositional quality of ewe's raw miles from commercial flocks under part-time grazing. *J. Dairy Res.* **2009**, *76*, 301-307. [CrossRef] [PubMed]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.