Effects of Mixed Adding Crude Extracts of β-Glucosidases from Three Different Non-Saccharomyces Yeast Strains on the Quality of Cabernet Sauvignon Wines

The aim of this study was to investigate the effects of crude extracts of β-glucosidase from Issatchenkia terricola SLY-4, Pichia kudriavzevii F2-24 and Metschnikowia pulcherrima HX-13 (termed as SLY-4E, F2-24E and HX-13E) on the flavor complexity and typicality of Cabernet Sauvignon wines. The grape must was fermented using Saccharomyces cerevisiae with single or mixed SLY-4E, F2-24E and HX-13E. The physicochemical characteristics, volatile aroma compounds, total anthocyanins and sensory attributes of the wines were determined. Adding SLY-4E, F2-24E and HX-13E in wines resulted in a decrease in the anthocyanin content, total acids and volatile acids in wines but an increase in the content of terpenes, benzene derivatives, higher alcohols and esters, which may enhance wine sensory qualities and result in loss of wine color. Different adding strategies of β-glucosidase led to a variety of effects on wine aroma. S/H/F-Ew significantly increased the content of benzene derivatives, higher alcohols and long-chain fatty acid esters, which enhanced the fruity and floral flavor of wines. F2-24E significantly increased the content of short- and medium-chain fatty acid esters, acetate esters and carbonyl compounds. The results indicated that the mixed addition of non-Saccharomyces crude extracts and co-fermentation with S. cerevisiae could further improve wine flavor quality.


Introduction
Wine is popular among consumers for its nutritional and healthy properties. Fermentation with Saccharomyces cerevisiae alone is a common method for commercial wine production to ensure its uniform flavor quality but usually causes poor flavor complexity and typicality of wine [1], resulting in the reduction in competitiveness. Under this context, improving the flavor complexity and typicality of wine has become the common goal pursued by scientists. It was reported that β-glucosidases could hydrolyze nonvolatile flavor precursors and generate aroma substances with volatile flavor [2,3]. Therefore, cofermenting with non-Saccharomyces and S. cerevisiae yeasts [4][5][6][7], or adding crude extracts or purified enzymes from non-Saccharomyces yeasts before fermentation [8][9][10] could improve the flavor complexity and typicality of wines. Moreover, different non-Saccharomyces yeast strains or their crude extracts of enzyme might produce different aroma compound profiles in wines [10]. Previous research indicated that Issatchenkia terricola SLY-4, Pichia kudriavzevii F2-24 and Metschnikowia pulcherrima HX-13 with β-glucosidases activities or their crude extracts of β-glucosidase (termed as SLY-4E, F2-24E and HX-13E) could improve the flavor complexity and typicality of wines and present different aroma compound profiles by co-fermenting with S. cerevisiae [10,11]. However, it is still unclear whether adding a combination of SLY-4E, F2-24E and HX-13E into must can further improve the flavor complexity and typicality of wines.

Analyzing the β-Glucosidase Activity
The β-glucosidase activity was analyzed according to the method described by Zhang et al. (2020) [10] with modifications. A total of 0.1 mL of crude extracts of β-glucosidase and 0.2 mL of 5 mmol/L p-nitrophenyl β-D-glucopyranoside (p-NPG) in phosphate buffer were mixed and incubated at 50 • C for 30 min; then, the reaction was inactivated by 2.0 mL 1.0 mol/L Na 2 CO 3 solution. The absorption value of the reaction solution was measured at 400 nm (Multiskan GO, Thermo Scientific, Waltham, MA, USA). The distilled water (200 µL) was added to replace p-NPG in control. One unit activity of β-glucosidase (U) was defined as the quantity of enzyme required to produce 1.0 µmol p-nitrophenol per minute under the reaction condition described.

Stability of β-Glucosidase in Crude Extracts
The SLY-4E, HX-13E and F2-24E extracts were added into modified P-C buffer (pH 4.0, containing 23.7% glucose and 12% ethanol) and stored at 20 • C for 14 days, respectively. The β-glucosidase activity was determined every 2 days.

Wine Fermentation
Must (237.1 g/L sugar calculated as glucose and 5.3 g/L acid calculated as tartaric acid) was incubated at 80 • C for 45 min. Then, 200 mL sterilized must containing 60 mg/L total SO 2 was macerated at 4 • C for 12 h in a 250 mL glass bottle. A total of 1U of SLY-4E, HX-13E and F2-24E extract was added to the must according to fermentation strategies (Table 1). Fermentation was carried out at 20 • C by inoculating 10 6 CFU/mL S. cerevisiae. Fermentation was considered and finished when the residual sugar was below 4 g/L.

Analysis of Wines
The content of alcohol, total acids and volatile acids was determined through methods recommended by the International Organization of Vine and Wine (OIV, 2021).
The anthocyanins content in wines was determined according to the methods described by Jungmin Lee et al. (2007) [12]. After centrifugation (8000 rpm, 20 min) of the wines (10 mL), the supernatant (1 mL) was collected and diluted appropriately with KCl buffer (0.025 mol/L, pH 1.0) and CH 3 COONa buffer (0.4 mol/L, pH 4.5), respectively. Then, the absorbance value of wines was determined at 520 nm and 700 nm after standing for 15 min in the dark at room temperature and distilled water was used as control. The anthocyanin content was expressed as cyanidin glucoside and calculated as follows: M represents the molecular weight of cyanidin glucoside (449.4 g/mol). DF represents the dilution rate of wines. ε represents the molar extinction coefficient for cyanidin glucoside (26,900 L/mol.cm). L represents the width of the cuvette (1 cm).
The volatile compounds were analyzed by headspace solid-phase micro-extraction according to Shi et al. (2019) [11]. Volatile compounds were extracted by 50/30 µm DVB/CAR/PDMS fiber (Supelco, Bellefonte, PA, USA) and analyzed by gas chromatograph (Agilent 6890N, Agilent, CA, USA) with a DB-5 capillary column (30 m × 0.32 mm × 0.25 µm) coupled to a mass spectrometer (Agilent 5975B, Agilent, CA, USA). Wines (8 mL), NaCl (2 g) and cyclohexanone (40 µg/L, 100 µL) as internal standard were added into a headspace bottle (20 mL) and bathed at 40 • C for 15 min. The fiber was pushed into the headspace for 30 min and immediately desorbed in the injector of gas chromatography at 250 • C for 5 min. Temperature increasing conditions of GC analysis were as follows: 40 • C to 130 • C at 3 • C/min, and then to 250 • C at 4 • C/min. The injector and detector temperatures were set at 250 • C and 260 • C, respectively. The mass spectrometer was operated in electron impact ionization mode at 70 eV, and ion source temperature was 250 • C. Volatile compounds were identified by comparing the MS fragmentation pattern of each compound with that in database Wiley 7.0 and NIST05. The following formula was used for the calculation of the compound content: Compound content (µg/mL) = GC peak areas of the compound × Quality of cyclohexanone (µg) GC peak area of cyclohexanone × Volume of wine sample (mL) The sensory quality of wines was analyzed according to the methods described by Shi et al. (2019) [11] with little modification.

Data Analyses
Data statistics and graphs were performed by Graphpad prism 6.0. One-way analysis of variance (ANOVA) was conducted by SPSS 19.0 (SPSS Inc., Chicago, IL, USA). Principal component analysis (PCA) was performed with SIMCA-P 14.1 (Umetrics AB, Umea, Sweden). Heatmap visualization analysis of the wines and their flavor compounds was performed with TBtools_JRE1.6 after the Z-score standardization.

Data Analyses
Data statistics and graphs were performed by Graphpad prism 6.0. One-way analysis of variance (ANOVA) was conducted by SPSS 19.0 (SPSS Inc., Chicago, IL, USA). Principal component analysis (PCA) was performed with SIMCA-P 14.1 (Umetrics AB, Umea, Sweden). Heatmap visualization analysis of the wines and their flavor compounds was performed with TBtools_JRE1.6 after the Z-score standardization.
In fact, the activity of β-glucosidase might be affected by the concentration of glucose or ethanol [14,15]. It was increased at low ethanol concentration by increasing the glycosyl transferase activity but it was decreased at high concentration of glucose (10%) or ethanol (≥10% v/v) by altering its protein structure [16][17][18][19]. These results indicated the activity of SLY-4E, HX-13E and F2-24E was affected by a high concentration of glucose and ethanol but was relatively stable during wine fermentation.

Physicochemical Characteristics of Wines
Physicochemical characteristics showed that the content of alcohol (12.07-12.56%) and total acids (5.20-5.81 g/L) were not significantly different among the wines ( Table 2). Compared with SCw (0.68 g/L), the volatile acids in wines (0.49-0.65 g/L) with added βglucosidase, except H-Ew (0.68 g/L), decreased, which had no relation with the adding In fact, the activity of β-glucosidase might be affected by the concentration of glucose or ethanol [14,15]. It was increased at low ethanol concentration by increasing the glycosyl transferase activity but it was decreased at high concentration of glucose (10%) or ethanol (≥10% v/v) by altering its protein structure [16][17][18][19]. These results indicated the activity of SLY-4E, HX-13E and F2-24E was affected by a high concentration of glucose and ethanol but was relatively stable during wine fermentation.

Physicochemical Characteristics of Wines
Physicochemical characteristics showed that the content of alcohol (12.07-12.56%) and total acids (5.20-5.81 g/L) were not significantly different among the wines ( Table 2). Compared with SCw (0.68 g/L), the volatile acids in wines (0.49-0.65 g/L) with added β-glucosidase, except H-Ew (0.68 g/L), decreased, which had no relation with the adding strategies of β-glucosidase ( Table 2). The volatile acid content of wines also decreased after being fermented by yeast strains with β-glucosidase activity [20,21]. Excessive volatile acids (>1.2 g/L) would bring undesirable flavors to wine [22]. These results indicate that adding SLY-4E, F2-24E and HX-13 into must would improve flavors by decreasing the volatile acid content. The reason for the effect of crude extracts of β-glucosidase on the volatile acid content in wines is still unclear. Further, the addition of SLY-4E, F2-24E and HX-13E in wines could also decrease the content of anthocyanin (163. .96 mg/L) ( Table 2). In particular, it shows that the anthocyanin content of H-Ew was 163.30 mg/L, which was lower than S-Ew and F-Ew. However, the anthocyanin contents of S/H-Ew (173.08 mg/L) and H/F-Ew (172.17 mg/L) and S/H/F-Ew (184.96 mg/L) were higher than H-Ew, which showed that the combinational addition of SLY-4E, F2-24E and HX-13E could attenuate the declining trend. Consistent with the result of Vernocchi et al. (2011) [21], the wines fermented by S. cerevisiae with higher β-glucosidase activity had lower anthocyanins concentrations. β-glucosidases can hydrolyze the glycosidic bond of anthocyanin to release free anthocyanidins, and free anthocyanidins easily degraded into colorless compounds, which would have a negative effect on the color of wine [23,24]. The results showed that adding crude extracts and purified β-glucosidase from yeasts could decrease the anthocyanin content of wines, and those from different yeasts would exhibit different hydrolyzing activity of anthocyanins. Note: different letters in the same column showed significant difference (p < 0.05).

Varietal Aroma Compounds of Wines
Eight varietal aroma compounds were detected in the wines. Compared with SCw (1.14 mg/L), the content of varietal aroma compounds in wines with added SLY-4E, F2-24E and HX-13E (1.26-1.86 mg/L) was significantly higher ( Table 3). The content of varietal aroma compounds in S/H/F-Ew was higher than S/H-Ew, S/F-Ew and H/F-Ew, which was higher than SLY-4E, F2-24E, and HX-13E. All the wines with added F2-24E had a higher content of varietal aroma compounds than other wines with SLY-4E or HX-13E. González-Pombo et al. (2011) [8] also reported that the purified β-glucosidase from I. terricola could significantly increase the content of varietal aroma compounds in white Muscat wine. It reflected that adding β-glucosidase could increase the varietal aroma content of wines, which usually exist as nonvolatile odorless glycosides, and β-glucosidase could hydrolyze the aroma compounds to release volatile aroma compounds and enhance the content of varietal aroma compounds [20,[25][26][27]. In addition, the content of volatile aroma compounds and the number of β-glucosidase combined added in wines are directly proportional. Therefore, it is likely that there was a synergistic effect between SLY-4E, F2-24E and HX-13E, which performed better in releasing varietal aroma compounds. Zietsman et al. (2010) [28] found a synergistic effect of both α-L-arabinofuranosidase and β-glucosidase towards diglycosidically bonded monoterpenes, which could enhance the varietal aroma compound content. However, there are no reports about the synergistic effect of β-glucosidase in increasing the content of varietal aroma compounds. As the PCA plot shows (Figure 2), 79.8% variance was explained by eight varietal aroma compounds, and PC1 and PC2 accounted for 65.2% and 14.6% variance, respectively. S/F-Ew, H/F-Ew and S/H/F-Ew grouped with 1-hexanol. S/H/F-Ew grouped with citronellol, 1-octen-3-ol and geraniol. S/F-Ew grouped with geranyl acetone, while other wines did not group with any of the varietal aroma compounds. Among the varietal aroma compounds, the odor active values (OAV) of linalool and 1-octen-3-ol were above 1.0, which would present floral and fruity flavor. These results indicate that adding mixed SLY-4E, F2-24E and HX-13E crude extract could effectively promote the releasing of all the detected varietal aroma compounds. Adding S/F-E could promote the releasing of geranyl acetone. Using various β-glucosidase to release varietal aroma compounds had a synergistic effect on releasing varietal aroma compounds, especially 1-hexanol. There are no reports about the synergistic effect of βglucosidase in releasing different aroma compounds. Systematic studies on the mechanism of the synergistic effect of β-glucosidases in increasing the varietal aroma compound should be investigated in the future. aroma compounds. As the PCA plot shows (Figure 2), 79.8% variance was explained by eight varietal aroma compounds, and PC1 and PC2 accounted for 65.2% and 14.6% variance, respectively. S/F-Ew, H/F-Ew and S/H/F-Ew grouped with 1-hexanol. S/H/F-Ew grouped with citronellol, 1-octen-3-ol and geraniol. S/F-Ew grouped with geranyl acetone, while other wines did not group with any of the varietal aroma compounds. Among the varietal aroma compounds, the odor active values (OAV) of linalool and 1-octen-3-ol were above 1.0, which would present floral and fruity flavor. These results indicate that adding mixed SLY-4E, F2-24E and HX-13E crude extract could effectively promote the releasing of all the detected varietal aroma compounds. Adding S/F-E could promote the releasing of geranyl acetone. Using various β-glucosidase to release varietal aroma compounds had a synergistic effect on releasing varietal aroma compounds, especially 1-hexanol. There are no reports about the synergistic effect of β-glucosidase in releasing different aroma compounds. Systematic studies on the mechanism of the synergistic effect of βglucosidases in increasing the varietal aroma compound should be investigated in the future.

Fermentative Aroma Compounds
The content of fermentative aroma compounds in S-Ew, H-Ew, F-Ew, S/H-Ew, S/F-Ew, H/F-Ew and S/H/F-Ew (338.21-514.14 mg/L) was significantly higher than in SCw (221.66 mg/L). S/H/F-Ew contained the highest aroma compounds (514.14 mg/L) ( Table  3). Among the fermentative aroma compounds, the content of benzene derivatives (93.70-136.26 181 mg/L), higher alcohols (188.42-336.42 mg/L) and esters (29.84-48.31 mg/L) in wines with added crude extract of β-glucosidase were higher, compared with those of SCw, but the content of fatty acids and carbonyl compounds showed no significant difference (Figure 3).
Although many studies have reported that the addition of β-glucosidase impacted the fermentation aroma, β-glucosidase from different strains affected the fermentation aroma differently [29]. It was shown that β-glucosidase from H. uvarum and Candida easanensis JK8 increased the contents of benzene derivatives in wines [30,31]. Phenyl ethanol was the most abundant compound in five benzene derivatives, which is one of the aglycons of aroma precursors in grapes and could be released by β-glucosidase [32,33]. Moreover, phenyl ethanol was the synthetic material for phenylethyl acetate and

Fermentative Aroma Compounds
The content of fermentative aroma compounds in S-Ew, H-Ew, F-Ew, S/H-Ew, S/F-Ew, H/F-Ew and S/H/F-Ew (338.21-514.14 mg/L) was significantly higher than in SCw (221.66 mg/L). S/H/F-Ew contained the highest aroma compounds (514.14 mg/L) ( Table 3). Among the fermentative aroma compounds, the content of benzene derivatives (93.70-136.26 181 mg/L), higher alcohols (188.42-336.42 mg/L) and esters (29.84-48.31 mg/L) in wines with added crude extract of β-glucosidase were higher, compared with those of SCw, but the content of fatty acids and carbonyl compounds showed no significant difference (Figure 3).
Although many studies have reported that the addition of β-glucosidase impacted the fermentation aroma, β-glucosidase from different strains affected the fermentation aroma differently [29]. It was shown that β-glucosidase from H. uvarum and Candida easanensis JK8 increased the contents of benzene derivatives in wines [30,31]. Phenyl ethanol was the most abundant compound in five benzene derivatives, which is one of the aglycons of aroma precursors in grapes and could be released by β-glucosidase [32,33]. Moreover, phenyl ethanol was the synthetic material for phenylethyl acetate and phenylacetaldehyde generated from phenylalanine by the Strecker degradation pathway, which may be affected by amino acid content and composition in musts [34]. Higher alcohols play a crucial role in wine aroma and could enhance floral and fruity flavor in wines [35]. It was reported that β-glucosidase from H. uvarum and Kloeckera apiculata increased the higher alcohol content of wine, especially isoamyl alcohol [2,36]. However, according to Zhang et al. (2020) [10], adding SLY-4E, F2-24E and HX-13E decreased the higher alcohol content. These differences could be attributed to the grape must used in this study containing a high amount of nitrogen. It was confirmed that amino acids could be transformed into higher alcohols via the Ehrlich pathway, and increased nitrogen in musts can result in higher alcohol contents in wines [37]. Esters present fruity and floral characteristics to wine, contributing to the presence of wine flavor [38]. β-glucosidase from different origins resulted in various levels of increase in the content of esters, showing high strain-dependent variations. Crude βglucosidase from Pichia fermentans and H. uvarum increased the content of ethyl octanoate, ethyl decanoate and ethyl laurate, which are medium-chain fatty acids esters [9,30]. On the other hand, β-glucosidase from Pichia anomala could increase the concentrations of ethyl hexanoate and ethyl decanoate [39]. It is possible that the substrate specificity of β-glucosidase varies from the different yeast strains and became more suitable for substrate hydrolysis in the ester formation biosynthetic pathways [40]. These results showed that adding β-glucosidases to wines could change the fermentative aroma profiles by increasing the content of benzene derivatives, higher alcohols and esters.  Wines were classified into five groups (Figure 4), which include (i) SCw, (ii) F-Ew, (iii) H/F-Ew, S/H-Ew and S/F-Ew, (iv) S/H/F-Ew, and (v) S-Ew and H-Ew. Hierarchical cluster analysis (HCA) showed that all the detected fermentative aroma compounds were classed into I, II, III and IV. SCw was abundant in Class I, which contained nonanal, 2-methylbutyric acid and isopentanoic acid and presented herb aroma. S-Ew and H-Ew were abundant in Class I and Class IV, which contained phenylacetaldehyde, ethyl heptanoate, 2-methyl-1-butanol and carbonyl compounds (2,3-pentanedione, decanal and 2,4-dimethoxybenzaldehyde). F-Ew was abundant in Class III and IV, including short-and medium-chain fatty acid esters, acetate esters and carbonyl compounds. In Class III, the OVA of isoamyl acetate, ethyl butyrate, ethyl hexanoate and ethyl octanoate was greater than 1, which improved fruity aroma in wines. H/F-Ew, S/H-Ew and S/F-Ew and S/H/F-Ew were higher in Class II, including benzene derivatives (except phenylacetaldehyde), higher alcohols (except 2-methyl-1-butanol), ethyl sebacate, long-chain fatty acid esters (ethyl laurate, ethyl tetradecanoate and hexadecanoic acid, and ethyl ester), and other esters and fatty acids. In Class II, the OVA of phenyl ethanol, phenethyl acetate, isoamyl alcohol, 3-methyl-1-pentanol, 1-nonanol, ethyl decanoate, ethyl laurate and isoamyl octanoate were greater than 1. These compounds presented fruity, fatty, floral and earthy notes and improved the complexity of aroma in H/F-Ew, S/H-Ew, S/F-Ew and S/H/F-Ew. This indicated that F2-24E was more favorable than SLY-4E and HX-13E in releasing short-and medium-chain fatty acid esters, acetate esters and carbonyl compounds when added alone. Adding mixed SLY-4E, F2-24E and HX-13E crude extract mainly improved the production of benzene derivatives, esters and higher alcohols, and S/H/F-Ew had the strongest ability to release these aroma substances. This reflects that mixed addition of SLY-4E, F2-24E and HX-13E affected the increase in acetate esters, short-and medium-chain fatty acids ester. However, there are no studies reporting the effect of mixed addition of β-glucosidase in releasing fermentative aroma compounds. These different aroma profiles of wines are probably caused by the other components in SLY-4E, F2-24E and HX-13E crude extract which affect the must composition of wines and need to further be determined after purified SLY-4E, F2-24E and HX-13E.

Sensory Analysis
Finally, a sensory analysis was carried out to evaluate the quality attributes of the wines, including flavor, aroma and color, which have effects on consumer acceptance ( Figure 5). Samples S-Ew, H-Ew, F-Ew, S/F-Ew, H/F-Ew, S/F-Ew and S/H/F-Ew had higher scores in floral (7.00, 7.11-8.22), fruity (7.22, 7.76-8.50), taste (6.00, 6.00-8.67) and acceptance (6.17, 6.33-8.50) but lower scores in appearance (9.17, 8.33-8.83) compared to SCw. Ma et al., (2017) reported that crude β-glucosidases from P. fermentans could promote the release of varietal aroma compounds and fortified fruity and floral traits [9]. Taste qualities in wine mainly contained sweetness, sourness and bitterness, contributed by sugars, organic acids and ethanol, respectively. Acetic acid is mostly responsible for the sour and vinegary smell and taste in wines [41,42]. We analyzed the volatile acid of wines, which mainly contained acetic acid. Therefore, adding SLY-4E, F2-24E and HX-13E to wines decreased the content of volatile acids and improved wine taste. In addition, there were other compounds that imparted better taste to S/H/F-Ew and this should be investigated in the future. The lower appearances score in wines with added SLY-4E, F2-24E and HX-13E was probably due to the decreased anthocyanin content [43]. S/F-Ew and S/H/F-Ew had higher appearance scores and anthocyanin content than other adding strategies, which demonstrated that S/F-E and S/H/F-E were more favorable for wine color development. Wines' quality attributes include flavor, aroma and color, which impact consumer acceptance. Results showed that adding SLY-4E, F2-24E and HX-13E could improve wine sensory qualities. scores in floral (7.00, 7.11-8.22), fruity (7.22, 7.76-8.50), taste (6.00, 6.00-8.67) and acceptance (6.17, 6.33-8.50) but lower scores in appearance (9.17, 8.33-8.83) compared to SCw. Ma et al., (2017) reported that crude β-glucosidases from P. fermentans could promote the release of varietal aroma compounds and fortified fruity and floral traits [9]. Taste qualities in wine mainly contained sweetness, sourness and bitterness, contributed by sugars, organic acids and ethanol, respectively. Acetic acid is mostly responsible for the sour and vinegary smell and taste in wines [41,42]. We analyzed the volatile acid of wines, which mainly contained acetic acid. Therefore, adding SLY-4E, F2-24E and HX-13E to wines decreased the content of volatile acids and improved wine taste. In addition, there were other compounds that imparted better taste to S/H/F-Ew and this should be investigated in the future. The lower appearances score in wines with added SLY-4E, F2-24E and HX-13E was probably due to the decreased anthocyanin content [43]. S/F-Ew and S/H/F-Ew had higher appearance scores and anthocyanin content than other adding strategies, which demonstrated that S/F-E and S/H/F-E were more favorable for wine color development. Wines' quality attributes include flavor, aroma and color, which impact consumer acceptance. Results showed that adding SLY-4E, F2-24E and HX-13E could improve wine sensory qualities.

Conclusions
In this study, three crude extracts of β-glucosidases from non-Saccharomyces yeasts SLY-4 (I. terricola), F2-24 (P. kudriavzevii) and HX-13 (M. pulcherrima) were evaluated under simulated winemaking conditions. Adding SLY-4E, F2-24E and HX-13E to wines decreased the content of total acids and volatile acids, and exhibited different hydrolyzing activity of anthocyanins. Moreover, adding single or mixed of SLY-4E, F2-24E and HX-13E to wines could increase the content of varietal aroma and fermentative aroma

Conclusions
In this study, three crude extracts of β-glucosidases from non-Saccharomyces yeasts SLY-4 (I. terricola), F2-24 (P. kudriavzevii) and HX-13 (M. pulcherrima) were evaluated under simulated winemaking conditions. Adding SLY-4E, F2-24E and HX-13E to wines decreased the content of total acids and volatile acids, and exhibited different hydrolyzing activity of anthocyanins. Moreover, adding single or mixed of SLY-4E, F2-24E and HX-13E to wines could increase the content of varietal aroma and fermentative aroma compounds, especially terpenes, benzene derivatives, higher alcohols, and esters, which enhanced the fruity and floral flavor of wines. Hence, we expect that the potential application of crude extracts from the non-Saccharomyces yeasts in winemaking can provide an approach to improve flavor complexity and characteristics of wines. However, it is still unclear what effect is possessed between multiple β-glucosidases and the effect on volatile acids and aroma compounds. In the future, investigations on the mechanism of the synergistic effect of β-glucosidases on increasing aroma compound varieties and content might be required for further studies.   Other esters Hexanoic acid, 3-methylbutyl ester 0.04 ± 0.00 e 0.05 ± 0.00 d 0.06 ± 0.00 c 0.09 ± 0.01 a 0.08 ± 0.00 b 0.08 ± 0.00 b 0.05 ± 0.00 d 0.06 ± 0.00 c 1 [35] <0.1 Apple, Pineapple [49] Hexanoic acid, 2-methylbutyl ester Nd 0.03 ± 0.00 b 0.02 ± 0.00 c 0.03 ± 0.00 b 0.03 ± 0.00 b 0.04 ± 0.00 a 0.02 ± 0.00 c 0.03 ± 0.00 b Isoamyl octanoate 0.09 ± 0.01 f 0.14 ± 0.00 e 0. 16   Note: SPSS software was used for ANOVA analysis, and the different letters in the same row show significant difference (p < 0.05). ND means no substance was detected.