fluids

Article Natural Yogurt Stabilized with Hydrocolloids from Butternut Squash ( moschata) Seeds: Effect on Physicochemical, Rheological Properties and Sensory Perception

Sergio A. Rojas-Torres, Somaris E. Quintana and Luis Alberto García-Zapateiro *

Research Group of Complex Fluid Engineering and Food Rheology, University of Cartagena, Cartagena 130015, Colombia; [email protected] (S.A.R.-T.); [email protected] (S.E.Q.) * Correspondence: [email protected]; Tel.: +57-675-2024

Abstract: Stabilizers are ingredients employed to improve the technological properties of products. The food industry and consumers have recently become interested in the development of natural ingredients. In this work, the effects of hydrocolloids from butternut squash () seeds (HBSS) as stabilizers on the physicochemical, rheological, and sensory properties of natural yogurt were examined. HBSS improved the yogurt’s physical stability and physicochemical properties, decreasing syneresis and modifying the samples’ rheological properties, improving the assessment of sensory characteristics. The samples presented shear thinning behavior characterized by a decrease in viscosity with the increase of the shear rate; nevertheless, the samples showed a two-step yield   stress. HBSS is an alternative as a natural stabilizer for the development of microstructured products.

Citation: Rojas-Torres, S.A.; Keywords: butternut squash (Cucurbita moschata); seeds; pulp; hydrocolloid; rheology; sensorial Quintana, S.E.; García-Zapateiro, L.A. Natural Yogurt Stabilized with Hydrocolloids from Butternut Squash (Cucurbita moschata) Seeds: Effect on Physicochemical, Rheological 1. Introduction Properties and Sensory Perception. Yogurt is a widely consumed dairy food that is internationally recognized for its health Fluids 2021, 6, 251. https://doi.org/ benefits, nutritional value, and digestibility [1,2]. Yogurts are produced by milk-controlled 10.3390/fluids6070251 fermentation combining symbiotic cultures of Streptococcus thermophilus and Lactobacillus delbrueckii spp. bulgaricus, resulting in a product with creamy characteristics, a typical Academic Editor: aroma, and a slightly acidic taste [3,4]. Antonio Santamaría The most important attributes for consumer acceptance of yogurts are texture and firmness, related to viscosity and stability. The different ingredients employed for yo- Received: 17 June 2021 gurt production such as nonfat dry milk, milk concentrate, and whey protein, Accepted: 7 July 2021 improve texture, mouthfeel, appearance, viscosity, and consistency and prevent whey Published: 9 July 2021 separation [5,6]. Moreover, the addition of phytochemicals, stabilizers, or other materials such as pectin, gelatin, inulin, or dietary fiber may enhance some sensory acceptance of Publisher’s Note: MDPI stays neutral yogurt and decrease syneresis [7]. Different hydrocolloids have been used as stabilizers for with regard to jurisdictional claims in yogurt to improve their acceptance and extend their shelf life [8]. published maps and institutional affil- Hydrocolloids are polysaccharides and utilized in industry to act as thickeners iations. in relatively low concentrations due to their high molecular weight and technological functionality [9]. They are generally used as stabilizers due to their hydrophilic properties such as water retention, emulsion stability, and texture modification to improve texture or enhance properties such as appearance, mouthfeel, consistency, viscosity, and prevent whey Copyright: © 2021 by the authors. separation temperature flocculation [6,10,11]. The rheological behavior of solution-added Licensee MDPI, Basel, Switzerland. hydrocolloids is affected by the hydrocolloid backbone’s structural features and its side This article is an open access article chains, molecular weight, and conformation of the hydrocolloid molecules, as well as the distributed under the terms and solvent conditions [12]. conditions of the Creative Commons Over the years, numerous papers on the natural and non-toxic stabilizer benefits Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ for human health have been published [13]; recently, pectin-rich ingredients have been 4.0/). investigated to improve the gel characteristics, microstructure, texture, and rheology of

Fluids 2021, 6, 251. https://doi.org/10.3390/fluids6070251 https://www.mdpi.com/journal/fluids Fluids 2021, 6, 251 2 of 17

yogurt through interactions with a casein network (i.e., orange fiber, apple pomace, carrot cell) [14–16]. Plant seeds have good potential as a new hydrocolloid source regarding their safety, usability, and low production costs. These include carob, flaxseed, and white mustard as potential sources of hydrocolloids. In the case of butternut squash (Cucurbita moschata) seeds, they are a source of pro- teins [17], polysaccharides [18], and bioactive components such as , provitamins, phytosterols, phospholipids, and fatty acids [19–21], also present technological properties such as solubility and emulsification capacity [22,23]. However, there are few reports on squash seed characterization and functional properties and their use in the development of food products. The present study was carried out to evaluate the effect of butternut squash seed hydrocolloids on the physicochemical, rheological, and sensory properties of yogurt.

2. Materials and Methods 2.1. Materials Butternut squash was purchased at the local food market in the city of Cartagena (Colombia). Ethanol (99.5% purity) and hexane were obtained from Panreac (Barcelona, Spain). NaOH, acetic acid, and phenolphthalein were purchased from Sigma–Aldrich (St. Louis, MO, USA). Xanthan gum was purchased from Tecnas (Medellin, Colombia). Commercial pasteurized and homogenized milk and skim milk powder were purchased from a local Colombian market. Starter culture (Lactobacillus bulgaricus and Streptococ- cus thermophilus) strains were used in yogurt production. All other reagents were of analytical grade.

2.2. Methods 2.2.1. Extraction of Hydrocolloids from Butternut Squash Seeds Squash seeds were selected by color, washed with sodium hypochlorite 100 ppm for 5 min, and washed with distilled water; after that, the seeds were freeze-dried at −50 ◦C and 0.02 Pa for 8 h employing the Labconco Freezone 1.5 Liter Benchtop Freeze Dry equipment (Kansas City, MO, USA). Dried seeds were degreased using an oil extractor (Dulong model DL-ZYJ05) to obtain an oil-free material. Next, the degreased seeds were ground employing an MF 10 basic microfine grinder drive coupled with an MF 10.1 (cutting–grinding head, IKA, Germany) to reduce the particle size. The extraction of hydrocolloids was done following the procedures described by Orgulloso-Bautista et al. [24] and Quintana et al. [25] with some modifications. The solubilization was carried out with distillate water at pH 10 for 4 h at 80 ◦C. The pH was adjusted employing acetic acid and NaOH. The mixture was separated by centrifugation for 15 min at 4000 rpm, and the supernatant was recollected. After that, the viscous solution was mixed with ethanol in a 1:1 ratio to precipitate the hydrocolloid-based extract and mixed for 2 h at 4.0 ± 0.5 ◦C. The mixture was centrifugated, and the precipitate was recollected and lyophilized during 48 h. Hydrocolloids from butternut squash (Cucurbita moschata) seeds were milled and stored at 5 ◦C until use.

2.2.2. Production of Natural Yogurt Different formulations of yogurt were prepared to employ different percentages of xanthan gum (XG) and hydrocolloids of butternut squash seed (HBSS) as stabilizers, as shown in Table1. A control sample (M1) was prepared free of stabilizers; the other samples were prepared using 1% of a blend of XG and HBSS (M3, M4, and M5). Samples employing only XG or HBSS were prepared (M2 and M6, respectively). Fluids 2021, 6, 251 3 of 17

Table 1. Experimental design of yogurt with the addition of hydrocolloids of butternut squash seed as stabilizers.

Xanthan Gum Hydrocolloids Sample Code mg/100 g mg/100 g M1 0 0 M2 1.00 0 M3 0.75 0.25 M4 0.50 0.50 M5 0.25 0.75 M6 0 1.00

Yogurt was prepared according to Ladejvardi et al. [26] and Aziznia et al. [27] with some modifications. Fresh skim milk was supplemented with skim milk powder until reaching 16% soluble solids. The blend milk was heated at 90 ◦C for 5 min and rapidly cooled to 42 ± 1 ◦C to inoculate the lactic starters. The cultured milk was divided into 6 equal parts (120 mL) in sterile glass containers to develop the samples with different percentages of XG–HBSS blends; after that, all samples were incubated at 40 ◦C until a pH value of 4.5 was achieved. Each different batch of blends was agitated and stored at 5 ◦C until analysis. The formulations were manufactured in triplicate.

2.2.3. Physicochemical Properties and Proximal Composition of Yogurt The physicochemical properties of yogurt were determinate after 24 h of processing following the method described by the Association of Official Analytical Chemistry [28]. The pH values were determined according to method 942.05/90 using a digital pH meter (Model HI 9124, Hanna Instruments, Woonsocket, RI, USA). Titratable acidity determina- tion (TTA) was analyzed according to method number 967.21; 10 g of the samples were diluted into 250 mL using distilled water. Next, an aliquot of 50 mL was blended with 0.2 mL phenolphthalein as the indicator and was titrated with 0.1 N NaOH; TTA was expressed as mg of lactic acid per 100 g of the sample. Density was measured using a gravimetric method with grease pycnometers [29]. Moisture was determined by dehydration in an oven at 105 ◦C; ash was determined by incineration at 550 ◦C until constant weight. The determination was done using hexane as a solvent from Soxhlet extraction for 4 h. The total protein content was determined by Kjeldahl and total per difference [28].

2.2.4. Syneresis Measurement The syneresis of yogurt was determined after 25 days of storage, following the proce- dures described by Säker and Rodriguez [30] with modifications. 10 g of yogurt sample were centrifuged at 5000 rpm for 20 min at 10 ◦C. After centrifugation, the clear super- natant was poured off, weighed, and used to determine the percentage (w/w) of syneresis (Equation (1)): g supernatant whey % Syneresis = × 100 (1) g sample weight

2.2.5. Color Parameters of Yogurt The color parameters of lightness (L∗), red–green color (a∗ : + : red; − : green), and yellow–blue color (b∗ : + : yellow; − : blue) of the yogurt were measured by the colorime- ter. The whiteness index (WI) and color variation (∆E∗) were calculated by Equations (2) and (3), respectively:

h i0.5 WI = 100 − (100 − L∗)2 + a∗2 + b∗2 (2)

where the subscript m stands for the sample of yogurt formulated with XG and HBSS, and c stands for the control sample (M1). Fluids 2021, 6, 251 4 of 17

2.2.6. Rheological Analysis The rheological characterization of yogurt was carried out after 48 h of preparation in a controlled-stress rheometer (Modular Advanced Rheometer System Haake Mars 60, Thermo-Scientific, Dreieich, Germany) equipped with a coaxial cylinder (inner radius 11.60 mm, outer radius 12.54 mm, cylinder length 37.6 mm) based on Quintana et al. [31]. Each sample was held at the temperature of assay for 600 s to ensure the same thermal and mechanical history for each sample. Viscous flow tests were carried out at a steady state, analyzing the variation of apparent viscosity (defined by the relationship between shear stress and shear rate) within a range of shear rates between 0.001 and 1000 s−1 at 5, 10, 15, and 25 ◦C. Small-amplitude oscillatory shear (SAOS) tests were performed to obtain viscoelastic responses. In this way, stress sweeps were carried out at a frequency of 1 Hz, applying an ascending series of stress values from 0.001 to 1000 Pa at 5 ◦C to determine the linear viscoelasticity interval. Frequency sweeps were performed to obtain the mechanical spectrum using a stress value within the linear viscoelastic range in a frequency range between 10−2 and 102 rad·s−1. The data recorded included the storage modulus (G0), which provided the elastic component; the loss modulus (G00 ), which was related to the viscous components of the material; and the loss tangent (Tan δ), which was the ratio G00 /G0 and provided the ratio of elastic to the viscous response of the material under consideration.

2.2.7. Sensorial Analysis Sensory evaluation was conducted by 30 tasters (15 males and 15 females, age 20–30 years, healthy subjects, lactose tolerant) using quantitative descriptive analysis. The panelists were instructed to evaluate each sample individually. Sensory analysis of the samples consisted of the evaluation of color, flavor, and consistency attributes employing the following hedonic scale descriptors: I dislike a lot = 1, I dislike a little = 2, I neither like nor dislike = 3, I like it = 4, and I like it a lot = 5. Equal amounts of each yogurt were prepared in 90 mL glass containers that were labeled, identified with three-digit random numbers, and equilibrated at room temperature for at least 1 h before consumption. The samples were served to panelists with a spoon and room-temperature water to cleanse the palate before presenting the samples, as described in technical guide GTC 165 [32]. The results were expressed as a percentage of acceptability.

2.2.8. Statistical Analysis The data obtained were analyzed by ANOVA (unidirectional) using Statgraphics software (Centurion Version 16.1) to determine statistically significant differences (p < 0.05) between the samples submitted to the characterizations.

3. Results and Discussion 3.1. Physicochemical Properties of Yogurt The physicochemical properties of the yogurt are shown in Table2. The pH and acidity contributed to the organoleptic properties of the yogurt. During the incubation process, the pH reached values between 4.54 and 4.57, contributing to the characteristic odor and flavor [33] due to the breakdown of lactose by lactic acid bacteria [34]. A decrease in this value favored the contraction of protein clots formed by lactic bacteria, affecting the yogurt’s sensory quality [35]. The acidity values (0.82 and 0.83% lactic acid) did not present significant variation (p > 0.05), preserving the flavor and texture of the yogurt with quality [36]. Next, the addition of XG and HBSS did not influence the density of yogurt (p > 0.05), presenting values between 1.06 and 1.08 g/mL that were within the normalized range [33]. Fluids 2021, 6, x FOR PEER REVIEW 5 of 17

Fluids 2021, 6, 251 5 of 17

Table 2. Physicochemical properties of yogurts with the addition of hydrocolloids of butternut squash seed as stabilizers. Table 2. Physicochemical properties of yogurts with the addition of hydrocolloids of butternut Acidity Density squashSample seed Code as stabilizers. pH % Acid Lactic g/mL M1 4.54 ± 0.01 a 0.85Acidity ± 0.005 a 1.06Density ± 0.00 a Sample Code pH M2 4.55 ± 0.01 a % 0.84 Acid ± 0.005 Lactic a 1.07 ±g/mL 0.00 a M3M1 4.54 4.57± ± 0.01 aa 0.85 0.83 ± 0.0050.005 aa 1.081.06 ± ±0.000.00 a a M4M2 4.55 4.56± ± 0.010.05 aa 0.84 0.83 ± 0.0050.005 aa 1.071.07 ± ±0.000.00 a a M5M3 4.57 4.57± ± 0.01 ab 0.83 0.82 ± 0.0050.005 aa 1.071.08 ± ±0.000.00 a a M6M4 4.56 4.56± ± 0.05 aa 0.83 ± 0.0050.005 aa 1.071.07 ± ±0.000.00 a a b a a Results are expressedM5 as mean ± standard4.57 ± 0.01 deviation. Different0.82 letters± 0.005 in the same row1.07 express± 0.00 sta- a a a tistically significantM6 differences (p < 4.56 0.05).± 0.05 0.83 ± 0.005 1.07 ± 0.00 Results are expressed as mean ± standard deviation. Different letters in the same row express statistically significant differences (p < 0.05). The syneresis in the product was a factor that required special attention due to the secretionThe of water syneresis on the in thegel productstructure. was Syne a factorresis can that reflect required a lower special structural attention stabiliza- due to the tion,secretion becoming of waterone of on the the worst gel structure. defects in Syneresis the final canproduct, reflect negatively a lower structural affecting stabilization, the con- sumer’sbecoming acceptance one of theand worst inferring defects the inmicrobiological the final product, problem negatively [37]. affectingThe syneresis the consumer’s of yo- gurtacceptance was evaluated and inferringduring 25 the days microbiological of storage (5 ± problem0.1 °C) (see [37 Figure]. The 1). syneresis The stabilizer-free of yogurt was sampleevaluated (M1) duringpresented 25 daysthe highest of storage syneresis (5 ± 0.1 percentages◦C) (see Figure compared1). The stabilizer-freeto the others sampleand showed(M1) an presented increase the over highest time, syneresisfollowed by percentages M2, M3, M4, compared M5, and to M6. the others and showed an increase over time, followed by M2, M3, M4, M5, and M6.

100 M1 90 M2 M3 80 M4 M5 70 M6 60 50 40 30 Syneresis (%) 20 10 0 0 5 10 15 20 25 30 Days

FigureFigure 1. Syneresis 1. Syneresis (%) (%)of yogurts of yogurts with with the the addition addition of ofhydrocollo hydrocolloidsids from from butternut butternut squash squash seeds seeds as as stabilizersstabilizers through through 25 25 days days of of storage. storage. Results Results represent represent average average values values ±± SDSD (n ( n= =3). 3).

SamplesSamples with with XG, XG, HBSS, HBSS, and and their their blends blends presented presented a significant a significant decrease decrease in their in their syneresissyneresis percentage percentage (p (

a linear increase with the addition of HBSS. The brightness of yogurt is related to fat globules and protein particle size, influencing the light reflectance of particles and scattering ability [40,41].

Table 3. Color parameters of yogurts with the addition of hydrocolloids of butternut squash seed as stabilizers. L*: luminosity; ∆E : change of color; WI: whiteness index.

Sample Code L* ∆E WI M1 62.81 ± 1.94 ª – 60.14 ± 0.36 ª M2 69.59 ± 0.40 ª 25.60 ± 2.63 ª 67.05 ± 0.36 b M3 70.11 ± 0.44 ab 58.20 ± 6.20 b 70.81 ± 0.72 ab M4 72.48 ± 1.25 b 61.13 ± 7.75 b 70.85 ± 0.93 ab M5 74.07 ± 0.68 b 79.15 ± 10.51 c 71.16 ± 2.18 b M6 76.11 ± 0.50 b 93.36 ± 5.87 d 73.88 ± 0.38 b Results are expressed as mean ± standard deviation. Different letters in the same row express statistically significant differences (p < 0.05).

∆E represents the change in the color of the yogurt formulated, concerning the refer- ence color (M1). The addition of XG and HBSS produced significant differences (p < 0.05), whereas M6 presented a difference of 75% from the control sample (M1). Similarly, the sam- ples presented a linear increase of WI (p < 0.05) with the addition of the HBSS. Changes in the WI can be related to changes in the structure, resulting in different light scattering [42].

3.3. Rheological Properties 3.3.1. Steady Shear Rate The viscous flow of the formulated yogurts containing different concentrations of XG and BHSS was related to their intrinsic characteristics such as kinetic stability and mean particle size. The samples presented a heterogeneous, three-dimensional network composed of continuously connected casein strands and whey protein aggregates [43,44]. All samples presented a nonlinear viscosity–shear rate relation (Figures A1 and A2), showing a decrease of apparent viscosity with the increase in shear rate, indicating a typical, non-Newtonian flow behavior-type shear thinning. This behavior was possibly due to the aligned, in- direction flow of stiff polymer molecules, decreasing the interaction between adjacent polymer molecule chains [45,46] and disrupting the flocculated casein micelle network while shear force was applied [47]. This viscous flow behavior can be observed in samples M2 and M6 in Figure2a. Nevertheless, the curve inflection points suggested apparent yield stress [48] in all cases. In the same way, yogurt with the addition of HBSS as a stabilizer showed decreased apparent yield stress in all batches, related to higher degrees of crosslinking of milk proteins in acidified milk [48]. The results suggested that shear (approximately 10 s−1) affected the yogurt structure. Figure2b shows the full viscosity–shear stress curves of samples M2 and M6 (AppendixA). All samples presented a two-step yielding point from the decline in vis- cosity. The extrapolation approach reported by Zhu et al. [49] was employed to obtain and compare the yield stress values of samples. The two stress values were associated with a static yield stress τo−1 (first decline) and fluidic yield stress, τ0−2 (second decline). The same behavior has been reported in the literature, e.g., in yogurts formulated with monk extract as a sweetener [50], yogurt containing selecan [51], and fat-free stirred yogurt [44], when the authors associated this behavior with the apparent yield stress [50] that could be due to breakage of the yogurt structure and the sensitive reaction of yogurts to shear stress due to its quasi-stable status [51]. However, there is no evidence of any analysis reports of these two areas regarding this behavior. Table4 shows the low shear viscosity ( η0) and yield stress (τo) parameters in the static and fluidic yield stress of all samples at different temperatures. Fluids 2021, 6, x FOR PEER REVIEW 6 of 17

luminosity values (L*) are attributed to their white color and bright. The samples pre- sented significant differences (p < 0.05) in comparison with stabilizer-free yogurt (M1), showing a linear increase with the addition of HBSS. The brightness of yogurt is related to fat globules and protein particle size, influencing the light reflectance of particles and scattering ability [40,41].

Table 3. Color parameters of yogurts with the addition of hydrocolloids of butternut squash seed as stabilizers. L*: luminosity; ∆𝐸: change of color; WI: whiteness index.

Sample Code L* ∆𝐄 WI M1 62.81 ± 1.94 ª -- 60.14 ± 0.36 ª M2 69.59 ± 0.40 ª 25.60 ± 2.63 ª 67.05 ± 0.36 b M3 70.11 ± 0.44 ab 58.20 ± 6.20 b 70.81 ± 0.72 ab M4 72.48 ± 1.25 b 61.13 ± 7.75 b 70.85 ± 0.93 ab M5 74.07 ± 0.68 b 79.15 ± 10.51 c 71.16 ± 2.18 b M6 76.11 ± 0.50 b 93.36 ± 5.87 d 73.88 ± 0.38 b Results are expressed as mean ± standard deviation. Different letters in the same row express sta- tistically significant differences (p < 0.05).

∆E represents the change in the color of the yogurt formulated, concerning the refer- ence color (M1). The addition of XG and HBSS produced significant differences (p < 0.05), whereas M6 presented a difference of 75% from the control sample (M1). Similarly, the samples presented a linear increase of WI (p < 0.05) with the addition of the HBSS. Changes in the WI can be related to changes in the structure, resulting in different light scattering [42].

3.3. Rheological Properties 3.3.1. Steady Shear Rate The viscous flow of the formulated yogurts containing different concentrations of XG and BHSS was related to their intrinsic characteristics such as kinetic stability and mean particle size. The samples presented a heterogeneous, three-dimensional network com- posed of continuously connected casein strands and whey protein aggregates [43,44]. All samples presented a nonlinear viscosity–shear rate relation (Figures A1 and A2), showing a decrease of apparent viscosity with the increase in shear rate, indicating a typical, non- Newtonian flow behavior-type shear thinning. This behavior was possibly due to the aligned, in-direction flow of stiff polymer molecules, decreasing the interaction between adjacent polymer molecule chains [45,46] and disrupting the flocculated casein micelle Fluids 2021, 6, 251 7 of 17 network while shear force was applied [47]. This viscous flow behavior can be observed in samples M2 and M6 in Figure 2a.

5 °C M2 5 °C 10 °C 10 °C 2 102 . 10 15 °C . 15 °C 20 °C 25 °c M6 ) / Pa) / s

) / Pa s 5 °C η 101 η 101 10 °C 15 °C 25 °C

100 100

10-1 10-1 Apparent viscosity ( viscosity Apparent Apparent viscosity ( viscosity Apparent Fluids 2021, 6, x FOR PEER REVIEW 7 of 17 10-2 -2 -4 -3 -2 -1 0 1 2 3 4 10 10 10 10 10 10 10 10 10 10 10-4 10-3 10-2 10-1 100 101 102 103 104 . -1 . Shear rate (γ) /s Shear rate (γ) /s-1

(a) (b)

103 M2 5 °C 5 °C 10 °C 10 °C 15 °C 102 15 °C 25 °C . . 2 10 M6 25 °C 5 °C 10 °C

) / Pa/ ) s 15 °C ) / Pa s η η 25 °C 101 101

0 10 100

-1 -1 Apparent viscosity ( viscosity Apparent 10 ( viscosity Apparent 10

-2 -1 0 1 2 10 10 10 10 10 10-2 10-1 100 101 102 Shear stress (τ) / Pa Shear stress (τ) / Pa

(c) (d)

Figure 2.Figure Flow 2.curvesFlow showing curves showing apparent apparent viscosity viscosity as functions as functions of the shear of the rate shear of samples rate of samples M1 (a) M1 (a) and M2 andand M2M6 and(b) and M6 (theb) and shear the stress shear of stress M1 ( ofc) M1and ( cM2) and and M2 M6 and (d) M6 of (yogurtsd) of yogurts with the with addition the addition of of hydrocolloidshydrocolloids from frombutternut butternut squash squash seeds seeds as stabilizers. as stabilizers.

Nevertheless, the curve inflection points suggested apparent yield stress [48] in all cases. InTable the same 4. Low way, shear yogurt viscosity with (η 0the) and addi yieldtion stress of HBSS (τo) parametersas a stabilizer of yogurts showed with decreased the addition of apparenthydrocolloids yield stress of in butternut all batches, squash related seeds to as stabilizershigher degrees at different of crosslinking temperatures. of milk pro- −1 teins in acidified milk [48]. TheTemp. results suggestedη that shearη (approximatelyτ 10 s ) affectedτ Sample Code 0 (First) 0 (Second) o−1 0−2 the yogurt structure. ◦C Pa·s Pa·s Pa Pa Figure 2bM1 shows the full 5 viscosity–shear 79.20 astress curves3.09 of * samples 0.15M2 band M6 (Ap-26.52 b pendix A). All samples presented10 a two-step 51.12 yieldinga point2.13 from+ the decline0.14 b in viscosity.30.44 b The extrapolation approach reported15 by Zhu 76.78 eta al. [49] was1.48 employed+ to0.13 obtainb and com-28.93 b a b b pare the yield stress values of25 samples. The 66.20 two stress values0.67 were associated0.13 with a static29.72 M2 5 256.3 b 2.99 * 0.66 a 55.61 b yield stress 𝜏 (first decline)10 and fluidic84.83 yieldb stress, 𝜏1.96+ (second decline).0.23 a The same22.14 b behavior has been reported in15 the literature,193.70 e.g.,b in yogurts1.68 + ·formulated0.37 witha monk fruit62.35 b b + a b extract as a sweetener [50], yogurt25 containing126.78 selecan [51],0.87 and fat-free stirred0.28 yogurt [44],40.96 M3 5 299.50 c 8.10 * 0.61 a 24.25 d when the authors associated10 this behavior 211.04 withc the apparent3.50 + yield stress0.43 [50]a that could20.07 d be due to breakage of the yogurt15 structure 265.04 andc the sensitive3.06 + ·reaction of0.52 yogurtsa to shear19.28 d stress due to its quasi-stable25 status [51]. 183.84However,c there1.30 is no evidence 0.28 ofa any analysis22.32 d d ab d reports of theseM4 two areas regarding 5 this behavior.93.93 8.15 * 0.15 23.32 10 220.07 d 4.28 + 0.48 ab 24.09 d Table 4 shows the low 15shear viscosity90.97 (𝜂d ) and yield1.92 stress+ (𝜏) parameters0.14 ab in 23.06the d static and fluidic yield stress25 of all samples66.28 at differentd temperatures.0.80 + 0.15 ab 23.03 d M5 5 37.52 e 3.04 * 0.12 b 35.82 c 10 30.75 e 2.05 + 0.12 b 27.57 b Table 4. Low shear viscosity (𝜂 ) and yield stress (𝜏) parameters of yogurts with the addition of 15 28.01 e 1.33 + 0.08 b 19.42 b hydrocolloids of butternut squash seeds as stabilizers at different temperatures. 25 29.30 e 2.06 0.12 b 26.87 b M6 5 14.40 f 8.65 * 0.05 ab 27.69 b Temp. 𝜼𝟎( 𝑭𝒊𝒓𝒔𝒕) 𝜼𝟎 (𝑺𝒆𝒄𝒐𝒏𝒅) 𝝉𝒐𝟏 𝝉𝟎𝟐 Sample Code 10 20.01 f 1.42 + 0.07 ab 19.75 c Pa·s Pa·s Pa Pa °C15 18.92 f 1.38 + 0.91 ab 12.51 c M1 255 79.207.20 f a 3.090.05 * 0.150.05 b ab 26.52 b11.07 c Different letters in the same10columns express51.12 statistically a significant2.13 + differences0.14 ( pb < 0.05) for30.44 temperature, b and different symbols express statistically15 significant76.78 differencesa (1.48p < 0.05) + for hydrocolloids.0.13 b 28.93 b 25 66.20 a 0.67 0.13 b 29.72 b M2 5 256.3 b 2.99 * 0.66 a 55.61 b 10 84.83 b 1.96 + 0.23 a 22.14 b 15 193.70 b 1.68 +· 0.37 a 62.35 b 25 126.78 b 0.87 + 0.28 a 40.96 b M3 5 299.50 c 8.10 * 0.61 a 24.25 d 10 211.04 c 3.50 + 0.43 a 20.07 d 15 265.04 c 3.06 +· 0.52 a 19.28 d 25 183.84 c 1.30 0.28 a 22.32 d M4 5 93.93 d 8.15 * 0.15 ab 23.32 d 10 220.07 d 4.28 + 0.48 ab 24.09 d 15 90.97 d 1.92 + 0.14 ab 23.06 d 25 66.28 d 0.80 + 0.15 ab 23.03 d Fluids 2021, 6, x FOR PEER REVIEW 8 of 17

M5 5 37.52 e 3.04 * 0.12 b 35.82 c 10 30.75 e 2.05 + 0.12 b 27.57 b 15 28.01 e 1.33 + 0.08 b 19.42 b 25 29.30 e 2.06 0.12 b 26.87 b M6 5 14.40 f 8.65 * 0.05 ab 27.69 b 10 20.01 f 1.42 + 0.07 ab 19.75 c 15 18.92 f 1.38 + 0.91 ab 12.51 c 25 7.20 f 0.05 0.05 ab 11.07 c Different letters in the same columns express statistically significant differences (p < 0.05) for tem- perature, and different symbols express statistically significant differences (p < 0.05) for hydrocol- Fluids 2021, 6, 251 loids. 8 of 17

From 𝜂, sample M1 (stabilizer-free) presented the lowest value. 𝜂 () did not show significant differences (p > 0.05) with the temperature, but it presented variation From η , sample M1 (stabilizer-free) presented the lowest value. η did not show with the percentage0 of gums (p < 0.05). The samples prepared with XG (M2)0 (First showed) the significant differences (p > 0.05) with the temperature, but it presented variation with the highest 𝜂 , i.e., the highest first-Newtonian region, followed by samples prepared with percentage of gums (p < 0.05). The samples prepared with XG (M2) showed the highest η , XG + BHSS, although the decrease of XG showed a reduction of 𝜂 . This high viscosity 0 i.e., the highest first-Newtonian region, followed by samples prepared with XG + BHSS, was attributed to dominating Brownian Forces at low shear rates that randomized the although the decrease of XG showed a reduction of η . This high viscosity was attributed to particles and formed many doublets and aggregates [52].0 𝜂 decreased with the dominating Brownian Forces at low shear rates that randomized () the particles and formed temperature (p < 0.05) and did not present significative differences with the percentage of many doublets and aggregates [52]. η0 (second) decreased with the temperature (p < 0.05) hydrocolloids.and did not present significative differences with the percentage of hydrocolloids. For the yield stress values, 𝜏 and 𝜏 did not present significant differences (p > For the yield stress values, τo−1 and τ0−2 did not present significant differences 0.05)(p with> 0.05) the with temperature the temperature but varied but wi variedth the with percentages the percentages of hydrocolloids of hydrocolloids (p < 0.05). (p < In 0.05). addition, 𝜏 in yogurt was linked to the breakdown of the interconnected network of In addition, τo−1 in yogurt was linked to the breakdown of the interconnected network of flocs/aggregates.flocs/aggregates. In contrast, In contrast, the thefurther further collapse collapse of flocs/aggregates of flocs/aggregates into into smaller smaller flocs flocs or or individual particles was associated with 𝜏 [53]. Based on these designs and rheological individual particles was associated withτ0−2 [53]. Based on these designs and rheological results,results, the theuse useof HBSS of HBSS as an as additive an additive for th fore theformulation formulation of microstructure of microstructure food food matrices matrices cancan be employed be employed alone alone or blended or blended with with other other hydrocolloids. hydrocolloids. TheThe temperature temperature dependence dependence of apparent of apparent viscosity viscosity at a at specific a specific temperature temperature was was evaluatedevaluated using using an Arrhenius-type an Arrhenius-type equation: equation:

𝐸 𝜂=𝜂 · E·  (3) η =η · 𝑅𝑇 a · (3) i RT −1 where 𝜂 is the viscosity at 10 s (𝜂 ), 𝜂 is a constant describing the viscosity coef- −1 . −1 ficientwhere at aη referenceis the viscosity temperature at 10 s (Pa·s),(ηγ= 10𝐸s− is1 ), theηi is activation a constant energy describing (J·mol the), viscosityR is the gas coeffi- −1 −1 −1 constantcient at(8.3144 a reference J·mol temperature·K ), and T (K) (Pa the·s), temperatureEa is the activation [54]. The energy activation (J·mol energy), R is(𝐸 the) is gas −1 −1 an importantconstant (8.3144 parameter J·mol of kinetic·K ), studies. and T (K) The the calculation temperature was [54 obtained]. The activation by plotting energy the ln ( Ea) (𝜂 is an important) against 1/T parameter (Figure 3). of kinetic studies. The calculation was obtained by plotting the . ln (ηγ=10s−1 ) against 1/T (Figure3).

2,0 M1 . M2 1,5 M3 M4

) / Pa s ) / M5 -1 1,0 M6 Arhenius equation = 10 s =

. γ η 0,5

0,0

Viscosity ln( -0,5

-1,0 0,0034 0,0035 0,0036 Temperature (1/T) / K-1

−1 . Figure 3. Viscosity at 10−1 s (η = s−1 ) versus 1/temperatures fitted to the Arrhenius equations of Figure 3. Viscosity at 10 s (𝜂 γ10) versus 1/temperatures fitted to the Arrhenius equations of yogurtyogurt samples samples with with additions additions of hydrocolloid of hydrocolloidss of butternut of butternut squash squash seeds seeds as stabilizers. as stabilizers. According to the linear fitting slope, the parameters adjusted to the Arrhenius equation are shown in Table5. The Ea of yogurts were 485.09 and 734.03 J/mol, which had a high degree of fit (R2 > 0.92). In this case, the sign inside the exponential function was positive because the viscosity decreased when the temperature increased [55]. The addition of stabilizers increased Ea values, showing a slight increase with the addition of HBSS. Therefore, the activation energy values indicated that the viscosity of the higher amounts of HBBS in the yogurt samples exhibited higher temperature dependence. Fluids 2021, 6, x FOR PEER REVIEW 9 of 17

According to the linear fitting slope, the parameters adjusted to the Arrhenius equa- tion are shown in Table 5. The Ea of yogurts were 485.09 and 734.03 J/mol, which had a high degree of fit (R2 > 0.92). In this case, the sign inside the exponential function was positive because the viscosity decreased when the temperature increased [55]. The addi- tion of stabilizers increased Ea values, showing a slight increase with the addition of HBSS. Therefore, the activation energy values indicated that the viscosity of the higher amounts Fluids 2021, 6, 251 9 of 17 of HBBS in the yogurt samples exhibited higher temperature dependence.

Table 5. Arrhenius equation parameters for yogurts with the addition of hydrocolloids of butternut squashTable 5.seedsArrhenius as stabilizers. equation parameters for yogurts with the addition of hydrocolloids of butternut squash seeds as stabilizers. 𝜼 Ea Sample Code 𝒊 R2 Pa·s J/mol η Ea SampleM1 Code i 0.063 × 10−6 a 528.55 a R2 0.98 Pa·s J/mol M2 3.728 × 10−11 a 561.05 b 0.92 −6 a a M1M3 0.063 × 10 2.314 × 10−12 a 528.55 485.09 c 0.98 0.98 × −11 a b M2M4 3.728 10 4.55 × 10−16 a 561.05 734.03 d 0.92 0.98 M3 2.314 × 10−12 a 485.09 c 0.98 M5 2.60 × 10−11 a 679.55 e 0.98 M4 4.55 × 10−16 a 734.03 d 0.98 M6 5.24 × 10−9 a 723.26 f 0.99 M5 2.60 × 10−11 a 679.55 e 0.98 Results are expressed as mean ± standard deviation. Different letters in the same column express M6 5.24 × 10−9 a 723.26 f 0.99 statistically significant differences (p < 0.05). Results are expressed as mean ± standard deviation. Different letters in the same column express statistically significant differences (p < 0.05). 3.3.2. Stress Sweep 3.3.2.The Stress stress Sweep sweep was done in all samples to determine the viscoelastic linear region (VLR),The evaluating stress sweep the behavior was done ofin the all storage samples (𝐺 to) determineand loss moduli the viscoelastic (𝐺 ). It was linear observed region that(VLR), in all evaluating cases, under the behaviorlow-stress of conditions, the storage the (G0 )modulus and loss moduli𝐺 and (𝐺G00 ).were It was linear, observed and whenthat inthe all stress cases, was under sufficiently low-stress large, conditions, there was the a modulus substantialG0 and decrease,G00 were showing linear, the and slope when tendedthe stress to zero. was sufficientlyAt this point, large, the there linearity was apr substantialofile deviated; decrease, it was showing considered the slope the critical tended stressto zero. considered At this point,as the thestatic linearity yield stress profile 𝜏. deviated; Next, a second it was decrease considered was the observed; critical stressnev- ertheless,considered the as stress the static kept yield increasing, stress τ cand. Next, it was a second difficult decrease to determine was observed; the dynamic nevertheless, yield stress.the stress Thus, kept the increasing, obtained curves and it corroborated was difficult tothe determine two-yield the point’s dynamic presence yield in stress. all cases, Thus, recognizedthe obtained by curvesthe slope’s corroborated change in the the two-yield stress sweep point’s curve presence (Figure in 4). all cases, recognized by the slope’s change in the stress sweep curve (Figure4).

103

102

101

100

G' G'' 10-1 M1 M2 M3 -2 M4

Loss modulus (G'') / Pa 10 M5 Storage modulus (G') / Pa M6

10-3 10-3 10-2 10-1 100 101 102 103 104 Stress (τ) / Pa

Figure 4. Stress sweep of yogurt samples with the addition of hydrocolloids of butternut squash seed Figure 4. Stress sweep◦ of yogurt samples with the addition of hydrocolloids of butternut squash seedas stabilizers as stabilizers at 5 atC. 5 °C. For all samples, different linear and nonlinear viscoelastic regions could be observed. For all samples, different linear and nonlinear viscoelastic regions could be observed. As shown in Table6, in the linear viscoelastic regions, G0 were higher than G00 ; remaining As shown in Table 6, in the linear viscoelastic regions, 𝐺 were higher than 𝐺; remain- constant indicated that the yogurt behaved like a viscoelastic solid. The samples prepared ing constant indicated that the yogurt behaved like a viscoelastic solid. The samples pre- without stabilizer (M1) and with XG or HBSS alone (M2 and M6) did not present significa- pared without stabilizer (M1) and with XG or HBSS alone (M2 and M6) did not present tive differences (p < 0.05) in τc, and the interaction of XG–HBSS (M3 and M4) presented higher τc, values, rationed with the crosslinking of the molecules. Fluids 2021, 6, x FOR PEER REVIEW 10 of 17

significative differences (p < 0.05) in 𝜏, and the interaction of XG–HBSS (M3 and M4) Fluids 2021, 6, 251 presented higher 𝜏, values, rationed with the crosslinking of the molecules. 10 of 17

Table 6. Critical stress 𝜏 and elastic 𝐺 and viscous 𝐺 moduli of yogurt samples with the addi- tion of hydrocolloids of butternut squash seeds as stabilizers at 5 °C. 0 00 Table 6. Critical stress τc and elastic G and viscous G moduli of yogurt samples with the addition ◦ Crossover of hydrocolloids of butternut squash𝝉 seeds as stabilizers𝑮 at 5 C. 𝑮′ Sample Code 𝒄 𝑮 = 𝑮 Pa Pa Pa 0 0 CrossoverPa τc G G 0 ” Sample Code a a a Gcd= G M1 Pa0.15 6.359Pa 3.07Pa 2.57 M2 0.15 a 45.34 b 11.60 b 4.92 bPa M1M3 0.150.25a ab 6.35953.74a c 14.103.07 c a 6.492.57 a cd M2M4 0.150.14a a 45.3415.39 bd 5.2811.60 d b 3.094.92 c b M3M5 0.250.53ab b 53.742.152 ce 1.6414.10 e c 1.476.49 a a M4M6 0.140.15a a 15.3911.65d f 4.535.28 f d 2.483.09 d c Results areM5 expressed as mean0.53 ± standardb deviatio2.152n. Differente letters 1.64in thee same column1.47 expressa statisticallyM6 significant differences 0.15 a (p < 0.05). 11.65 f 4.53 f 2.48 d Results are expressed as mean ± standard deviation. Different letters in the same column express statistically 3.3.3.significant Frequency differences Sweep (p < 0.05). 3.3.3.The Frequency structural Sweepstability of the storage (𝐺 ) and loss moduli (𝐺 ), the complex modulus (𝐺∗), and the loss tangent (tan δ) are the main rheological parameters reported for yogurts The structural stability of the storage (G0) and loss moduli (G00 ), the complex modulus [56–58]. The frequency sweeps of yogurt showed than the storage moduli (𝐺) were (G∗), and the loss tangent (tan δ) are the main rheological parameters reported for yo- higher than the loss moduli (𝐺) over the entire range of frequencies studied, indicating gurts [56–58]. The frequency sweeps of yogurt showed than the storage moduli (G0) were that the elastic components were greater than the viscous, as shown in Figure 5. higher than the loss moduli (G00 ) over the entire range of frequencies studied, indicating that the elastic components were greater than the viscous, as shown in Figure5.

102

101

100 G' G'' M1 M2 Loss modulus (G'') / Pa

Storagemodulus (G') / Pa M3 M4 -1 M5 10 M6

10-1 10 0 101 Angular frequency (ω) / rad. s-1

Figure 5. Frequency sweep of yogurts with additions of hydrocolloids of butternut squash seeds as Figure 5. Frequency sweep of yogurts with additions of hydrocolloids of butternut squash seeds as stabilizers at 5 ◦C. stabilizers at 5 °C. The solid-like properties (G0 > G00 ) were observed to be dominant in XG + HBSS yogurt.The solid-like Higher G properties0 and no crossover (𝐺 >𝐺 with) wereG00 observedwere documented to be dominant as strong in XG gel-like + HBSS behavior, yo- gurt.attributed Higher to𝐺 the and interaction no crossover of casein with with𝐺 were stabilizers documented because as electrostatic strong gel-like bond behavior, formation attributedat a low to pH the might interaction have led of casein to dense with packing stabilizers of the because protein electr gelostatic structure bond with formation hydrocol- at loidsa low [ 56pH,59 might]. Values have of ledG0 toat denseω = 1 radpacking·s−1 (Tableof the7 protein) were usedgel structure to compare with the hydrocol- effects of loidsstabilizers [56,59]. onValues the samples, of 𝐺 at presenting 𝜔 = 1 rad·s a decrease−1 (Table with7) were the additionused to compare of HBSS (thep < 0.05),effects while of stabilizersG0 was inon the the following samples, order:presenting M2 ≥ a M3decrease≥ M5 with > M6. the Nevertheless, addition of HBSS the stabilizer-free (p < 0.05), whilesample 𝐺 waspresented in the crosslinkingfollowing order: at high M2frequencies, ≥ M3 ≥ M5 showing > M6. Nevertheless, the weakest behavior,the stabilizer- which freedid sample not occur presented with the crosslinking samples that at contained high frequencies, XG and HBSS, showing and the stabilizerweakest behavior, free sample whichcould did lose not the occur desired with characteristics the samples that by the contained consumer XG when and HBSS, the set and of internal the stabilizer macromolec- free sampleular networks could lose of thethe yogurtdesired lost characteristics the ability to trapby the internal consumer fluids when [60]. Moreover,the set of theinternal sample macromolecularwith XG (M2) presented networks a of high theG yogurt0 in the lost same the way ability of the to samplestrap internal with 0.25fluids and [60]. 0.5 More- g/100 g over,of HBSSthe sample (M3 and with M4, XG respectively),(M2) presented explained a high 𝐺 by in the the ability same way of XG of tothe interact samples with with the milk protein, altering the gel network [16]. The samples prepared only with HBSS (M6) presented intermediate values of G0 and G00 with gel-like behavior. The differences seen in the rheological behavior may have been due to the creation of protein–casein complexes. Similar values were obtained by Kumar and Sasmai [61] in the development of milk gel with the addition of Cucurbita seed extract. Fluids 2021, 6, x FOR PEER REVIEW 11 of 17

0.25 and 0.5 g/100g of HBSS (M3 and M4, respectively), explained by the ability of XG to interact with the milk protein, altering the gel network [16]. The samples prepared only with HBSS (M6) presented intermediate values of 𝐺 and 𝐺 with gel-like behavior. The Fluids 2021, 6, 251 differences seen in the rheological behavior may have been due to the creation of protein– 11 of 17 casein complexes. Similar values were obtained by Kumar and Sasmai [61] in the devel- opment of milk gel with the addition of Cucurbita seed extract.

0 00 Table 7.Table Values 7. Valuesfor the forviscoelastic the viscoelastic parameters parameters of 𝐺 of(storageG (storage modulus) modulus) and and𝐺 (lossG (loss modulus) modulus) in in the ◦ −1 the frequencyfrequency sweep sweep of yogurt of yogurt samples samples at 5 at °C 5 atC 1 atrad·s 1 rad−1. ·s .

0 G𝑮 𝑮G” SampleSample Code Code Tan 𝜹Tan δ Pa PaPa M1 1.25 a 0.70 a 0.56 a a a a M2 M132.04 1.25 b 8.970.70 a 0.28 b0.56 b a b M3 M2 32.0434.35 c 9.338.97 a 0.28 b0.28 c a b M4 M3 34.353.80 d 2.089.33 a 0.55 a0.28 a a M5 M4 10.883.80 de 3.352.08 a 0.30 b0.55 M6 M5 10.882.05 f e 1.303.35 a a 0.63 a0.30 b Different letters inM6 the same column express2.05 statisticallyf significant1.30 differencesa (p < 0.05). 0.63 a Different letters in the same column express statistically significant differences (p < 0.05). The obtained values of the loss tangent (Tan 𝛿 =𝐺⁄𝐺) were between 0.27 and 0.63, = 00 0) indicating aThe wide obtained range of values behavior of the between loss tangent a weak (Tan gelδ andG a/ strongG were elastic between gel. The 0.27 ad- and 0.63, dition ofindicating XG decreased a wide the range Tan 𝛿 of associated behavior with between the formation a weak gel of anda stronger a strong gel; elastic next, the gel. The addition of XG decreased the Tan δ associated with the formation of a stronger gel; next, the addition of HBSS (sample with XG + HBSS) indicated a modification of the gel structure. addition of HBSS (sample with XG + HBSS) indicated a modification of the gel structure. Samples with XB+HBSS presented similar behavior, and then the employee of HBSS (M6) Samples with XB+HBSS presented similar behavior, and then the employee of HBSS (M6) increased the viscous behavior associated with the weaker gel structure. The addition of increased the viscous behavior associated with the weaker gel structure. The addition of HBSS is an alternative for developing food products to increase their viscosity or thickness HBSS is an alternative for developing food products to increase their viscosity or thickness and improve the water binding ability and texture of yogurt. The choice of hydrocolloid and improve the water binding ability and texture of yogurt. The choice of hydrocolloid depends on the attributes desired in the final product and the conditions of processing. depends on the attributes desired in the final product and the conditions of processing.

3.4. Sensory3.4. Sensory Analysis Analysis SensorySensory evaluation evaluation of food of products food products is a quality is a quality index of index organoleptic of organoleptic properties properties such such as texture,as texture, flavor, flavor,and appearance and appearance that is thatvital is for vital the foroverall the overall acceptability acceptability by consumers. by consumers. The concentrationThe concentration of hydrocolloids of hydrocolloids plays an plays important an important role in rolegoverning in governing yogurt’s yogurt’s stability stability and textureand texture and finds and a findscorrelation a correlation with improv with improveded quality quality and sensory and sensory perception perception [57,58]. [ 57,58]. FigureFigure 6 shows6 shows the average the average results results of the ofevaluated the evaluated parameters. parameters.

M1 5,0

4,5 M6 M2 4,0

3,5

3,0

M5 M3

M4

Figure 6. Global quality of yogurts with additions of hydrocolloids from butternut squash seed Figure 6. Global quality of yogurts with additions of hydrocolloids from butternut squash seed as as stabilizers. stabilizers. The samples prepared only with HBSS (M6) showed the highest average values for all Theassessed samples items, prepared which only corresponded with HBSS to (M6) “I like showed it a lot”, the while highest the samples average prepared values for with XG all assessedand blended items, which with hydrocolloids corresponded corresponded to “I like it a lot”, to “I while neither the like samples nor dislike”. prepared Global with quality XG andinvolved blended all with attributes. hydrocolloids The color corresponded characteristic to “I showed neither statisticallike nor dislike”. differences Global (p < 0.05) between control samples with XG + HBSS, which resulted in greater scattering of light, and therefore, the products contained less luminosity and whiteness [62]. Nevertheless, similar results were obtained with the control samples and yogurt prepared with HBSS, produced by the scattering of light by the milk constituents: fat globules, casein micelles, colloidal calcium phosphate, some pigments, and riboflavin [63]. This sensory parameter’s importance lies in the fact that it is one of the leading causes of a dairy product being Fluids 2021, 6, 251 12 of 17

bought by the consumer or rejected. However, it does not reflect the same taste or nutri- tional value [62]. Similarly, the sample with HBSS presented the highest odor acceptance (p < 0.05) (Table8).

Table 8. Sensorial analysis of yogurts with the addition of hydrocolloids of butternut squash seeds as stabilizers.

Parameters Color Odor Flavor Consistency M1 4.10 ± 0.71 ab 3.70 ± 0.70 ab 4.16 ± 0.59 a 3.36 ± 1.32 a M2 4.00 ± 0.72 a 3.76 ± 0.77 ab 3.90 ± 0.99 a 3.23 ± 1.07 a M3 3.83 ± 0.87 a 3.56 ± 0.72 a 3.86 ± 1.10 a 3.46 ± 0.93 a M4 3.86 ± 0.86 a 3.56 ± 0.81 a 3.90 ± 0.92 a 3.26 ± 0.98 a M5 4.16 ± 0.74 ab 3.96 ± 0.55 bc 3.96 ± 0.96 a 3.76 ± 0.93 a M6 4.45 ± 0.72 ab 4.26 ± 0.78 c 4.26 ± 0.86 a 3.63 ± 1.29 a Results are expressed as mean ± standard deviation. Different letters in the same column express statistically significant differences (p < 0.05).

According to the panelists, the addition of blends of hydrocolloids (XG+HBSS) did not significantly affect the consistency (p > 0.05). The sensorial parameters such as creaminess, thickness, consistency, smoothness, and taste perception were related strongly to the rheo- logical properties. The addition of hydrocolloids modified the rheology and improved the physical stability and overall mouthfeel properties to increase desirable general characteris- tics [64]. The samples with a moderate viscosity and storage and loss moduli presented the highest values of the evaluation of sensorial parameters. Considering the sensory analysis results altogether, it could be concluded that the application of HBSS as a food additive could be possible, with positive effects on the texture and sensory attributes and without the consequent appearance of unwanted effects.

3.5. Proximal Composition of a Selected Standardized Product Table9 shows the data from the proximal composition of M6, considered the best sample due to its attributes upon sensory evaluation carried out in comparison with the control sample. The proximal composition results indicated that the selected yogurt (M6) presented values of 8.18 ± 0.43%, a moisture value of 79.58 ± 0.47%, a fat value of 4.68 ± 0.11%, a fiber value of 0.89 ± 0.02%, a protein value of 6.10 ± 0.18%, and an ash value of 0.56 ± 0.05%. The data were compared with common and probiotic yogurts and presented higher carbohydrate and protein contents related to the content of HBSS due to the high content of carbohydrate and proteins in butternut squash seeds [17].

Table 9. Proximal compositions of yogurt with additions of hydrocolloids from butternut squash seeds (M6), common, and probiotic yogurt.

Parameters M6 Common Yogurt * Probiotic Yogurt * Moisture (%) 79.58 ± 0.47 88.72 88.74 Ash (%) 0.56 ± 0.05 0.72 0.71 Fat (%) 4.68 ± 0.11 2.45 2.42 Carbohydrates (%) 8.18 ± 0.43 5.24 5.20 Proteins (%) 6.10 ± 0.18 2.90 2.92 Fiber (%) 0.89 ± 0.02 – – * Source: [65]. Fluids 2021, 6, 251 13 of 17

4. Conclusions The hydrocolloids from butternut squash (Cucurbita moschata) seeds were evaluated as stabilizers in yogurt samples to preserve the physicochemical properties, flavor, and texture of samples as well as their activity on yogurt starter bacteria. Yogurts present a non-Newtonian flow behavior-type shear thinning, with inflection points in the curves (approximately 10 s−1) related to the presence of a yield point. The evaluation of stress in the function of the shear rate showed a two-step yielding point: static yield stress (first decline) and fluidic yield stress (second decline) associated with the breakdown of the interconnected network and the further collapse of flocs/aggregates into smaller flocs or individual particles; in addition, the viscoelastic properties of the yogurts presented solid- like properties in a typical strong gel-like behavior. HBSS contributed to the increase in the desirable overall characteristics that presented the highest valorization related to the highest WI and acceptance values for panelists. The addition of HBSS is an alternative for increasing the viscosity or thickness of food products and improving the water binding ability and texture of yogurt. The choice of hydrocolloid depends on the attributes desired in the final product and the conditions of processing. The employment of different hydrocolloids affects product yield, texture, mouthfeel, sensory properties, and consumer acceptance.

Author Contributions: Conceptualization, S.A.R.-T., S.E.Q. and L.A.G.-Z.; Methodology, S.A.R.-T. and L.A.G.-Z.; Software, S.E.Q.; Validation, S.E.Q. and L.A.G.-Z.; Formal Analysis, S.A.R.-T., S.E.Q. and L.A.G.-Z.; Investigation, S.A.R.-T. and L.A.G.-Z.; Resources, L.A.G.-Z.; Writing—Original Draft Preparation, S.A.R.-T., S.E.Q. and L.A.G.-Z.; Writing—Review and Editing, S.E.Q. and L.A.G.-Z.; Supervision, L.A.G.-Z.; Project Administration, L.A.G.-Z.; Funding Acquisition, S.E.Q. and L.A.G.-Z. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the Ministry of Science, Technology, and Innovation— MinCiencias (Contract 368-2019, Research Project No. 110780864755). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Conflicts of Interest: The authors declare no conflict of interest. Fluids 2021, 6, x FOR PEER REVIEW 13 of 17

4. Conclusions The hydrocolloids from butternut squash (Cucurbita moschata) seeds were evaluated as stabilizers in yogurt samples to preserve the physicochemical properties, flavor, and texture of samples as well as their activity on yogurt starter bacteria. Yogurts present a non-Newtonian flow behavior-type shear thinning, with inflection points in the curves (approximately 10 s−1) related to the presence of a yield point. The evaluation of stress in the function of the shear rate showed a two-step yielding point: static yield stress (first decline) and fluidic yield stress (second decline) associated with the breakdown of the interconnected network and the further collapse of flocs/aggregates into smaller flocs or individual particles; in addition, the viscoelastic properties of the yogurts presented solid- like properties in a typical strong gel-like behavior. HBSS contributed to the increase in the desirable overall characteristics that presented the highest valorization related to the highest WI and acceptance values for panelists. The addition of HBSS is an alternative for increasing the viscosity or thickness of food products and improving the water binding ability and texture of yogurt. The choice of hydrocolloid depends on the attributes desired in the final product and the conditions of processing. The employment of different hydro- colloids affects product yield, texture, mouthfeel, sensory properties, and consumer ac- ceptance.

Author Contributions: Conceptualization, S.A.R.-T., S.E.Q., and L.A.G.-Z.; Methodology, S.A.R.-T. and L.A.G.-Z.; Software, S.E.Q.; Validation, S.E.Q. and L.A.G.-Z.; Formal Analysis, S.A.R.-T., S.E.Q., and L.A.G.-Z.; Investigation, S.A.R.-T. and L.A.G.-Z.; Resources, L.A.G.-Z.; Writing—Original Draft Preparation, S.A.R.-T., S.E.Q. and L.A.G.-Z.; Writing—Review and Editing, S.E.Q. and L.A.G.-Z.; Supervision, L.A.G.-Z.; Project Administration, L.A.G.-Z.; Funding Acquisition, S.E.Q. and L.A.G.-Z. All authors have read and agreed to the published version of the manuscript.

Funding: This research was funded by the Ministry of Science, Technology, and Innovation—Min- Ciencias (Contract 368-2019, Research Project No. 110780864755). Institutional Review Board Statement: Not applicable. Fluids 2021, 6, 251 Informed Consent Statement: Not applicable. 14 of 17 Conflicts of Interest: The authors declare no conflict of interest.

AppendixAppendix A A

5 °C 5 °C 10 °C 10 °C 102 2 . 15 °C . 10 15 °C 20 °C 20 °C ) / Pa s ) / Pa s η 1 10 η 1 10

100 100

10-1 -1

Apparent viscosity ( 10 Apparent viscosity (

Fluids 2021, 6, x10 FOR-2 PEER REVIEW 14 of 17 10-4 10-3 10-2 10-1 100 101 102 103 104 10-2 10-4 10-3 10-2 10-1 100 101 102 103 104 γ. -1 Shear rate ( ) /s . Shear rate (γ) /s-1

(a) (b)

5 °C 5 °C 10 °C 10 °C 2 . 10 15 °C . 102 15 °C 20 °C 20 °C ) / Pa s ) / Pa s ) η η 101 101

100 100

10-1 Apparent viscosity ( Apparent viscosity ( viscosity Apparent 10-1

10-2 -4 -3 -2 -1 0 1 2 3 4 10 10 10 10 10 10 10 10 10 10-4 10-3 10-2 10-1 100 101 102 103 104 . -1 . Shear rate (γ) /s Shear rate (γ) /s-1

(c) (d)

5 °C 5 °C 10 °C 10 °C 2 2 . 10 15 °C . 10 15 °C 25 °C 25 °C ) / Pa s / Pa ) ) / Pa s / Pa ) η η 101 101

100 100

10-1 10-1 Apparent viscosity ( viscosity Apparent ( viscosity Apparent

10-2 10-2 10-4 10-3 10-2 10-1 100 101 102 103 104 10-4 10-3 10-2 10-1 100 101 102 103 104 . . Shear rate (γ) /s-1 Shear rate (γ) /s-1

(e) (f)

FigureFigure A1. A1.Flow Flow curves curves show show apparent apparent viscosity viscosity in in the the shear shear rate rate of of yogurts yogurts with with the the addition addition of of hydrocolloids hydrocolloids of of butternut butternut squashsquash seeds seeds as as stabilizers: stabilizers: (a) ( M1,a) M1, (b )(b M2,) M2, (c) ( M3,c) M3, (d )(d M4,) M4, (e) ( M5,e) M5, (f) ( M6.f) M6.

5 °C 5 °C 10 °C 10 °C 2 2 . 10 15 °C . 10 15 °C 25 °C 20 °C ) / Pa s ) / Pa s ) / η η 101 101

100 100

-1 -1

Apparent viscosity ( viscosity Apparent 10 ( viscosity Apparent 10

10-2 10-1 100 101 102 10-2 10-1 100 101 102 Shear stress (τ) / Pa Shear stress (τ) / Pa

(a) (b)

5 °C 5 °C 10 °C 10 °C 2 2 15 °C . 10 15 °C . 10 20 °C 20 °C ) / Pa s ) / ) / Pa s ) / η

η 1 101 10

100 100

-1 -1 ( viscosity Apparent 10

Apparent viscosity ( 10

10-2 10-1 100 101 102 10-2 10-1 100 101 102 Shear stress (τ) / Pa Shear stress (τ) / Pa

(c) (d) Fluids 2021, 6, x FOR PEER REVIEW 14 of 17

5 °C 5 °C 10 °C 10 °C 2 . 10 15 °C . 102 15 °C 20 °C 20 °C ) / Pa s ) / Pa s ) η η 101 101

100 100

10-1 Apparent viscosity ( Apparent viscosity ( viscosity Apparent 10-1

10-2 -4 -3 -2 -1 0 1 2 3 4 10 10 10 10 10 10 10 10 10 10-4 10-3 10-2 10-1 100 101 102 103 104 . -1 . Shear rate (γ) /s Shear rate (γ) /s-1

(c) (d)

5 °C 5 °C 10 °C 10 °C 2 2 . 10 15 °C . 10 15 °C 25 °C 25 °C ) / Pa s / Pa ) ) / Pa s / Pa ) η η 101 101

100 100

10-1 10-1 Apparent viscosity ( viscosity Apparent ( viscosity Apparent

10-2 10-2 10-4 10-3 10-2 10-1 100 101 102 103 104 10-4 10-3 10-2 10-1 100 101 102 103 104 . . Shear rate (γ) /s-1 Shear rate (γ) /s-1

(e) (f) Fluids 2021Figure, 6, 251 A1. Flow curves show apparent viscosity in the shear rate of yogurts with the addition of hydrocolloids of butternut15 of 17 squash seeds as stabilizers: (a) M1, (b) M2, (c) M3, (d) M4, (e) M5, (f) M6.

5 °C 5 °C 10 °C 10 °C 2 2 . 10 15 °C . 10 15 °C 25 °C 20 °C ) / Pa s ) / Pa s ) / η η 101 101

100 100

-1 -1

Apparent viscosity ( viscosity Apparent 10 ( viscosity Apparent 10

10-2 10-1 100 101 102 10-2 10-1 100 101 102 Shear stress (τ) / Pa Shear stress (τ) / Pa

(a) (b)

5 °C 5 °C 10 °C 10 °C 2 2 15 °C . 10 15 °C . 10 20 °C 20 °C ) / Pa s ) / ) / Pa s ) / η

η 1 101 10

100 100

-1 -1 ( viscosity Apparent 10 Fluids 2021, 6Apparent viscosity ( , 10x FOR PEER REVIEW 15 of 17

10-2 10-1 100 101 102 10-2 10-1 100 101 102 Shear stress (τ) / Pa Shear stress (τ) / Pa

(c) (d)

5 °C 5 °C 10 °C 10 °C 2 2 . 10 15 °C . 10 15 °C 20 °C 20 °C ) / Pa s Pa / ) ) ) / Pa s η η 1 101 10

0 100 10

-1

-1 Apparent viscosity ( 10 Apparent viscosity ( 10

-2 -1 0 1 2 -2 -1 0 1 2 10 10 10 10 10 10 10 10 10 10 τ Shear stress (τ) / Pa Shear stress ( ) / Pa

(e) (f) FigureFigure A2. A2.Flow Flow curves curves showing showing apparent apparent viscosity viscosity in functionsin functions of shear of shear stress stress of yogurts of yogurts with with the addition the addition of hydrocolloids of hydrocol- ofloids butternut of butternut squash squash seeds as seeds stabilizers as stabilizers (a) M1, ( (ab)) M1, M2, ( (bc) M2, M3, ( cd)) M3, M4, ( (de)) M4, M5, ( (ef) M6.M5, (f) M6.

ReferencesReferences 1.1. Desobry-Banon,Desobry-Banon, S.; S.; Vetier, Vetier, N.; N.; Hardy, Hardy, J. HealthJ. Health Benefits Benefits of of Yogurt Yogurt Consumption. Consumption. A A Review. Review.Int. Int. J. J. Food Food Prop. Prop.1999 1999, ,2 ,2, 1–12. 1–12, [CrossRefdoi:10.1080/10942919909524585.] 2.2. Shah,Shah, N.P. N.P. Health Health Benefits Benefits of of Yogurt Yogurt and and Fermented Fermented Milks. Milks. In InManufacturing Manufacturing Yogurt Yogurt and and Fermented Fermented Milks Milks; Blackwell; Blackwell Publishing: Publishing: Ames,Ames, IA, IA, USA, USA, 2007; 2007; pp. pp. 327–351. 327–351. 3.3. Andrade,Andrade, R.D.; R.D.; Arteaga, Arteaga, M.R.; M.R.; Simanca, Simanca, M.M. M.M. Efecto Efecto Del SalvadoDel Salvado de Trigo de Trigo En El En Comportamiento El Comportamiento Reológico Reológico Del Yogurt Del deYogurt Leche de deLeche Búfala. de Inf.Búfala. Tecnol Inf.ógica Tecnológica2010, 21 ,2010 117–124., 21, 117–124, [CrossRef doi:10.4067/s0718-07642010000500015.] 4.4. Mendoza,Mendoza, R.; R.; Trujillo, Trujillo, Y.; Y.; Duran, Duran, D. D. Evaluaci Evaluaciónón Del Del Almid Almidónón de de Ñame Ñame Espino Espino (Dioscorea (Dioscorea rotundata rotundata) Como) Como Estabilizante Estabilizante En En La La ElaboraciElaboraciónón de de Yogur Yogur Entero Entero Tipo Tipo Batido—Dialnet. Batido—Dialnet.Rev. Rev. Fac. Fac. De De Cienc. Cienc. B áBásicassicas 2007 2007, 5,, 5 97–105., 97–105. 5.5. Lee,Lee, W.J.; W.J.; Lucey, Lucey, J.A. J.A. Formation Formation and and Physical Physical Properties Properties of of Yogurt. Yogurt.Asian-Australas. Asian-Australas. J. Anim.J. Anim. Sci. Sci.2010 2010, 23, 23, 1127–1136., 1127–1136. [CrossRef ] 6.6. Sanchez,Sanchez, C.; C.; Zuniga-Lopez, Zuniga-Lopez, R.; Schmitt, Schmitt, C.; C.; Despond, Despond, S.; S.; Hardy, Hardy, J. Microstructure J. Microstructure of Acid-Induced of Acid-Induced Skim Skim Milk-Locust Milk-Locust Bean Bean Gum- Gum-XanthanXanthan Gels. Gels. Int. DairyInt. Dairy J. 2000 J. 2000, 10, ,199–212,10, 199–212. doi:10.1016/S0958-6946(00)00030-3. [CrossRef] 7.7. Sigdel,Sigdel, A.; A.; Ojha, Ojha, P.; P.; Karki, Karki, T.B. T.B. Phytochemicals Phytochemicals and and Syneresis Syneresis of of Osmo-Dried Osmo-Dried Mulberry Mulberry Incorporated Incorporated Yoghurt. Yoghurt.Food Food Sci. Sci. Nutr. Nutr. 20182018, 6, ,6 1045–1052., 1045–1052, [CrossRef doi:10.1002/fsn3.645.] 8. Cárdenas Casanova, A.; Alvites Gutierrez, H.; Valladares Castillo, G.; Obregón Domínguez, J.; Vásquez-Villalobos, V. Optimi- zation by Mixtures Design of Syneresis and Sensory Texture of Natural Smoothie Yogurt Using Three Types of Hydrocolloids. Agroind. Sci. 2013, 3, 35–40, doi:10.17268/agroind.science.2013.01.04. 9. Balestra, F.; Petracci, M. Technofunctional Ingredients for Meat Products. In Sustainable Meat Production and Processing; Elsevier: Amsterdam, The Netherlands, 2019; pp. 45–68. 10. Fiszman, S.M.; Salvador, A. Effect of Gelatine on the Texture of Yoghurt and of Acid-Heat-Induced Milk Gels. Eur. Food Res. Technol. 1999, 208, 100–105, doi:10.1007/s002170050383. 11. Schmidt, K.A.; Smith, D.E. Milk Reactivity of Gum and Milk Protein Solutions. J. Dairy Sci. 1992, 75, 3290–3295, doi:10.3168/jds.S0022-0302(92)78104-1. 12. Norton, I.T.; Spyropoulos, F.; Cox, P. Practical Food Rheology: An Interpretive Approach; Blackwell: Oxford, UK, 2011. 13. Kirschweng, B.; Tátraaljai, D.; Földes, E.; Pukánszky, B. Natural Antioxidants as Stabilizers for Polymers. Polym. Degrad. Stab. 2017, 145, 25–40. 14. Kieserling, K.; Vu, T.M.; Drusch, S.; Schalow, S. Impact of Pectin-Rich Orange Fibre on Gel Characteristics and Sensory Proper- ties in Lactic Acid Fermented Yoghurt. Food Hydrocoll. 2019, 94, 152–163, doi:10.1016/j.foodhyd.2019.02.051. 15. McCann, T.H.; Fabre, F.; Day, L. Microstructure, Rheology and Storage Stability of Low-Fat Yoghurt Structured by Carrot Cell Wall Particles. Food Res. Int. 2011, 44, 884–892, doi:10.1016/j.foodres.2011.01.045. 16. Wang, X.; Kristo, E.; LaPointe, G. The Effect of Apple Pomace on the Texture, Rheology and Microstructure of Set Type Yogurt. Food Hydrocoll. 2019, 91, 83–91, doi:10.1016/j.foodhyd.2019.01.004. 17. Vinayashree, S.; Vasu, P. Biochemical, Nutritional and Functional Properties of Protein Isolate and Fractions from (Cucurbita Moschata Var. Kashi Harit) Seeds. Food Chem. 2021, 340, 128177, doi:10.1016/j.foodchem.2020.128177. 18. Wang, L.; Liu, F.; Wang, A.; Yu, Z.; Xu, Y.; Yang, Y. Purification, Characterization and Bioactivity Determination of a Novel Polysaccharide from Pumpkin (Cucurbita Moschata) Seeds. Food Hydrocoll. 2017, 66, 357–364, doi:10.1016/j.foodhyd.2016.12.003. 19. Caili, F.; Huan, S.; Quanhong, L. A Review on Pharmacological Activities and Utilization Technologies of Pumpkin. Plant Foods Hum. Nutr. 2006, 61, 73–80. 20. Fruhwirth, G.O.; Hermetter, A. Seeds and Oil of the Styrian Oil Pumpkin: Components and Biological Activities. Eur. J. Lipid Sci. Technol. 2007, 109, 1128–1140, doi:10.1002/ejlt.200700105. Fluids 2021, 6, 251 16 of 17

8. Cárdenas Casanova, A.; Alvites Gutierrez, H.; Valladares Castillo, G.; Obregón Domínguez, J.; Vásquez-Villalobos, V. Optimiza- tion by Mixtures Design of Syneresis and Sensory Texture of Natural Smoothie Yogurt Using Three Types of Hydrocolloids. Agroind. Sci. 2013, 3, 35–40. [CrossRef] 9. Balestra, F.; Petracci, M. Technofunctional Ingredients for Meat Products. In Sustainable Meat Production and Processing; Elsevier: Amsterdam, The Netherlands, 2019; pp. 45–68. 10. Fiszman, S.M.; Salvador, A. Effect of Gelatine on the Texture of Yoghurt and of Acid-Heat-Induced Milk Gels. Eur. Food Res. Technol. 1999, 208, 100–105. [CrossRef] 11. Schmidt, K.A.; Smith, D.E. Milk Reactivity of Gum and Milk Protein Solutions. J. Dairy Sci. 1992, 75, 3290–3295. [CrossRef] 12. Norton, I.T.; Spyropoulos, F.; Cox, P. Practical Food Rheology: An Interpretive Approach; Blackwell: Oxford, UK, 2011. 13. Kirschweng, B.; Tátraaljai, D.; Földes, E.; Pukánszky, B. Natural Antioxidants as Stabilizers for Polymers. Polym. Degrad. Stab. 2017, 145, 25–40. [CrossRef] 14. Kieserling, K.; Vu, T.M.; Drusch, S.; Schalow, S. Impact of Pectin-Rich Orange Fibre on Gel Characteristics and Sensory Properties in Lactic Acid Fermented Yoghurt. Food Hydrocoll. 2019, 94, 152–163. [CrossRef] 15. McCann, T.H.; Fabre, F.; Day, L. Microstructure, Rheology and Storage Stability of Low-Fat Yoghurt Structured by Carrot Cell Wall Particles. Food Res. Int. 2011, 44, 884–892. [CrossRef] 16. Wang, X.; Kristo, E.; LaPointe, G. The Effect of Apple Pomace on the Texture, Rheology and Microstructure of Set Type Yogurt. Food Hydrocoll. 2019, 91, 83–91. [CrossRef] 17. Vinayashree, S.; Vasu, P. Biochemical, Nutritional and Functional Properties of Protein Isolate and Fractions from Pumpkin (Cucurbita Moschata Var. Kashi Harit) Seeds. Food Chem. 2021, 340, 128177. [CrossRef] 18. Wang, L.; Liu, F.; Wang, A.; Yu, Z.; Xu, Y.; Yang, Y. Purification, Characterization and Bioactivity Determination of a Novel Polysaccharide from Pumpkin (Cucurbita Moschata) Seeds. Food Hydrocoll. 2017, 66, 357–364. [CrossRef] 19. Caili, F.; Huan, S.; Quanhong, L. A Review on Pharmacological Activities and Utilization Technologies of Pumpkin. Plant Foods Hum. Nutr. 2006, 61, 73–80. [CrossRef][PubMed] 20. Fruhwirth, G.O.; Hermetter, A. Seeds and Oil of the Styrian Oil Pumpkin: Components and Biological Activities. Eur. J. Lipid Sci. Technol. 2007, 109, 1128–1140. [CrossRef] 21. Sener, B.; Orhan, I.; Ozcelik, B.; Kartal, M.; Aslan, S.; Ozbilen, G. Antimicrobial and Antiviral Activities of Two Seed Oil Samples of L. and Their Fatty Acid Analysis. Nat. Prod. Commun. 2007.[CrossRef] 22. Rezig, L.; Chibani, F.; Chouaibi, M.; Dalgalarrondo, M.; Hessini, K.; Guéguen, J.; Hamdi, S. Pumpkin () Seed Proteins: Sequential Extraction Processing and Fraction Characterization. J. Agric. Food Chem. 2013, 61, 7715–7721. [CrossRef] [PubMed] 23. Buˇcko,S.; Katona, J.; Popovi´c,L.; Vaštag, Ž.; Petrovi´c,L.; Vuˇcini´ce-Vasi´c,M. Investigation on Solubility, Interfacial and Emulsifying Properties of Pumpkin (Cucurbita Pepo) Seed Protein Isolate. LWT Food Sci. Technol. 2015, 64, 609–615. [CrossRef] 24. Orgulloso-Bautista, S.; Ortega-Toro, R.; Alberto, L.; Zapateiro, G. Design and Application of Hydrocolloids from Butternut Squash (Cucurbita moschata) Epidermis as a Food Additive in Mayonnaise-Type Sauces. ACS Omega 2021.[CrossRef] 25. Quintana Martínez, S.E.; Torregroza Fuentes, E.E.; García Zapateiro, L.A. Food Hydrocolloids from Butternut Squash (Cucurbita moschata) Peel: Rheological Properties and Their Use in Carica Papaya Jam. ACS Omega 2021, 6, 12114–12123. [CrossRef][PubMed] 26. Ladjevardi, Z.S.; Gharibzahedi, S.M.T.; Mousavi, M. Development of a Stable Low-Fat Yogurt Gel Using Functionality of Psyllium (Plantago Ovata Forsk) Husk Gum. Carbohydr. Polym. 2015, 125, 272–280. [CrossRef] 27. Aziznia, S.; Khosrowshahi, A.; Madadlou, A.; Rahimi, J. Whey Protein Concentrate and Gum Tragacanth as Fat Replacers in Nonfat Yogurt: Chemical, Physical, and Microstructural Properties. J. Dairy Sci. 2008, 91, 2545–2552. [CrossRef] 28. AOAC. Association of Official Analytical Chemist Official Methods of Analysis, 17th ed.; AOAC: Gaithersburg, MD, USA, 2000. 29. Santillán-Urquiza, E.; Méndez-Rojas, M.Á.; Vélez-Ruiz, J.F. Fortification of Yogurt with Nano and Micro Sized Calcium, Iron and , Effect on the Physicochemical and Rheological Properties. LWT 2017, 80, 462–469. [CrossRef] 30. Säker Vásquez, W.; Zevallos, A.R. Influencia Del Cultivo Láctico, Gelatina y Sacarosa Sobre La Viscosidad, Sinéresis y Carac- terísticas Sensoriales En Leche Fermentada. Pueblo Cont. 2016, 23, 125–136. 31. Quintana-Martinez, S.; Morales-Cano, A.; García-Zapateiro, L. Rheological Behaviour in the Interaction of Lecithin and Guar Gum for Oil-in-Water Emulsions. Czech J. Food Sci. 2018, 36, 73–80. [CrossRef] 32. ICONTEC Analisis Sensorial. Guia General-GTC 165; ICONTEC Analisis Sensorial: Bogotá, Colombia, 1999. 33. Aportela-Palacios, A.; Sosa-Morales, M.E.; Velez-Ruiz, J.F. Rheological and Physicochemical Behavior of Fortified Yogurt, With Fiber and Calcium. J. Texture Stud. 2005, 36, 333–349. [CrossRef] 34. Raju, P.N.; Pal, D. Effect of Dietary Fibers on Physico-Chemical, Sensory and Textural Properties of Misti Dahi. J. Food Sci. Technol. 2014, 51, 3124–3133. [CrossRef] 35. Marulanda, M.; Granados, C.; García-Zapateiro, L.A. Análisis Sensorial y Estimación Fisicoquímica de Vida Útil de Una Bebida Tipo Yogur a Base de Lactosuero Dulce Fermentada Con Estreptococcus Salivarius Ssp. Thermophilus y Lactobacillus Casei Ssp. Casei. Prod. Limpia 2016, 11, 94–102. [CrossRef] 36. Granata, L.A.; Morr, C.V. Improved Acid, Flavor and Volatile Compound Production in a High Protein and Fiber Soymilk Yogurt-like Product. J. Food Sci. 1996, 61, 331–336. [CrossRef] 37. Lucey, J.A. Formation and Physical Properties of Milk Protein Gels. J. Dairy Sci. 2002, 85, 281–294. [CrossRef] Fluids 2021, 6, 251 17 of 17

38. Molina Chew, I.S. Comparación de Tres Estabilizantes Comerciales Utilizados En La Elaboración de Yogurt de Leche Descremada de Vaca. Licenciatura Thesis, Universidad de San Carlos de Guatemala, de Guatemala, Guatemala, 2009. 39. Macedo Ramírez, R.C.; Vélez-Ruíz, J.F. Propiedades Fisicoquímicas y de Flujo de Un Yogur Asentado Enriquecido Con Mi- crocápsulas Que Contienen Ácidos Grasos Omega 3. Inf. Tecnol. 2015, 26, 87–96. [CrossRef] 40. Tamime, A.Y.; Robinson, R.K. Yoghurt-Science and Technology; Woodhead Publishing: Sawston, UK, 1999. 41. Ciron, C.I.E.; Gee, V.L.; Kelly, A.L.; Auty, M.A.E. Modifying the Microstructure of Low-Fat Yoghurt by Microfluidisation of Milk at Different Pressures to Enhance Rheological and Sensory Properties. Food Chem. 2012, 130, 510–519. [CrossRef] 42. Truong, T.; Palmer, M.; Bansal, N.; Bhandari, B. Effect of Milk Fat Globule Size on Functionalities and Sensory Qualities of Dairy Products; Springer: Cham, Swizerland, 2016; pp. 47–67. 43. Andoyo, R.; Guyomarc’h, F.; Burel, A.; Famelart, M.H. Spatial Arrangement of Casein Micelles and Whey Protein Aggregate Inacid Gels: Insight on Mechanisms. Food Hydrocoll. 2015, 51, 118–128. [CrossRef] 44. Laiho, S.; Williams, R.P.W.; Poelman, A.; Appelqvist, I.; Logan, A. Effect of Whey Protein Phase Volume on the Tribology, Rheology and Sensory Properties of Fat-Free Stirred Yoghurts. Food Hydrocoll. 2017, 67, 166–177. [CrossRef] 45. Hosseini-Parvar, S.H.; Matia-Merino, L.; Goh, K.K.T.; Razavi, S.M.A.; Mortazavi, S.A. Steady Shear Flow Behavior of Gum Extracted from Ocimum Basilicum L. Seed: Effect of Concentration and Temperature. J. Food Eng. 2010, 101, 236–243. [CrossRef] 46. Huang, T.; Tu, Z.; Zou, Z.; Shangguan, X.; Wang, H.; Bansal, N. Glycosylated Fish Gelatin Emulsion: Rheological, Tribological Properties and Its Application as Model Coffee Creamers. Food Hydrocoll. 2020, 102, 105552. [CrossRef] 47. Quemada, D.; Berli, C. Energy of Interaction in Colloids and Its Implications in Rheological Modeling. Adv. Colloid Interface Sci. 2002, 98, 51–85. [CrossRef] 48. Jaros, D.; Heidig, C.; Rohm, H. Enzymatic Modification Through Microbial Transglutaminase Enhances the Viscosity of Stirred Yogurt. J. Texture Stud. 2007, 38, 179–198. [CrossRef] 49. Zhu, L.; Sun, N.; Papadopoulos, K.; De Kee, D. A Slotted Plate Device for Measuring Static Yield Stress. J. Rheol. 2001, 45, 1105–1122. [CrossRef] 50. Ban, Q.; Liu, Z.; Yu, C.; Sun, X.; Jiang, Y.; Cheng, J.; Guo, M. Physiochemical, Rheological, Microstructural, and Antioxidant Properties of Yogurt Using Monk Fruit Extract as a Sweetener. J. Dairy Sci. 2020, 103, 10006–10014. [CrossRef][PubMed] 51. Fu, R.; Li, J.; Zhang, T.; Zhu, T.; Cheng, R.; Wang, S.; Zhang, J. Salecan Stabilizes the Microstructure and Improves the Rheological Performance of Yogurt. Food Hydrocoll. 2018, 81, 474–480. [CrossRef] 52. Pal, R. Effect of Droplet Size on the Rheology of Emulsions. AIChE J. 1996, 42, 3181–3190. [CrossRef] 53. Shakeel, A.; Kirichek, A.; Chassagne, C. Yield Stress Measurements of Mud Sediments Using Different Rheological Methods and Geometries: An Evidence of Two-Step Yielding. Mar. Geol. 2020, 427, 106247. [CrossRef] 54. Laity, P.R.; Holland, C. Thermo-Rheological Behaviour of Native Silk Feedstocks. Eur. Polym. J. 2017, 87, 519–534. [CrossRef] 55. Garvín, A.; Ibarz, R.; Ibarz, A. Kinetic and Thermodynamic Compensation. A Current and Practical Review for Foods. Food Res. Int. 2017, 96, 132–153. [CrossRef] 56. Crispín-Isidro, G.; Lobato-Calleros, C.; Espinosa-Andrews, H.; Alvarez-Ramirez, J.; Vernon-Carter, E.J. Effect of Inulin and Agave Fructans Addition on the Rheological, Microstructural and Sensory Properties of Reduced-Fat Stirred Yogurt. LWT Food Sci. Technol. 2015, 62, 438–444. [CrossRef] 57. Nguyen, P.T.M.; Kravchuk, O.; Bhandari, B.; Prakash, S. Effect of Different Hydrocolloids on Texture, Rheology, Tribology and Sensory Perception of Texture and Mouthfeel of Low-Fat Pot-Set Yoghurt. Food Hydrocoll. 2017, 72, 90–104. [CrossRef] 58. Xu, K.; Guo, M.; Du, J.; Zhang, Z. Okra Polysaccharide: Effect on the Texture and Microstructure of Set Yoghurt as a New Natural Stabilizer. Int. J. Biol. Macromol. 2019, 133, 117–126. [CrossRef] 59. Aguirre-Mandujano, E.; Lobato-Calleros, C.; Beristain, C.I.; Garcia, H.S.; Vernon-Carter, E.J. Microstructure and Viscoelastic Properties of Low-Fat Yoghurt Structured by Monoglyceride Gels. LWT Food Sci. Technol. 2009, 42, 938–944. [CrossRef] 60. Kavanagh, G.M.; Ross-Murphy, S.B. Rheological Characterisation of Polymer Gels. Prog. Polym. Sci. 1998, 23, 533–562. [CrossRef] 61. Kumar, A.; Sasmal, S. Rheological and Physico-Chemical Properties of Milk Gel Using Isolate of Pumpkin (Cucurbita moschata) Seeds: A New Source of Milk Clotting Peptidase. Food Hydrocoll. 2020, 106, 105866. [CrossRef] 62. Harper, W.; Hall, C. Dairy Technology and Engineering, 2th ed.; AVI. Publishers: Lubumbashi, Zaire, 1981. 63. Carpenter, R.P. Análisis Sensorial En El Desarrollo y Control de La Calidad de Alimentos; Editorial Acribia, S.A.: Zaragoza, Spain, 2002. 64. Marcotte, M.; Hoshahili, A.R.T.; Ramaswamy, H.S. Rheological Properties of Selected Hydrocolloids as a Function of Concentration and Temperature. Food Res. Int. 2001, 34, 695–703. [CrossRef] 65. Lollo, P.C.B.; De Moura, C.S.; Morato, P.N.; Cruz, A.G.; Castro, W.; de Freitas Castro, W.; Betim, C.B.; Nisishima, L.; Faria, J.; de Assis, F.; et al. Probiotic Yogurt Offers Higher Immune-Protection than Probiotic Whey Beverage. Food Res. Int. 2013, 54, 118–124. [CrossRef]