Impact of Chemical and Alternative Fungicides Applied to Grapevine cv Nebbiolo on Microbial Ecology and Chemical-Physical Grape Characteristics at Harvest

Viticulture is a cropping system in which treatment against fungal diseases (in particular powdery and downy mildews) can be extremely frequent. Accordingly, a reduction in antimicrobial treatments and the application of environmentally-friendly compounds are becoming increasingly important for a more sustainable viticulture. In addition to their effect against pathogens, the impact of these products on the quality of the grapes is very important for the oenological industries, but unfortunately at present few data are available. We evaluated the effect of the application of biocontrol products and resistance inducers in the vineyard on the mechanical properties, microbial ecology, technological and phenolic maturity of Vitis vinifera “Nebbiolo” grapes at harvest. The yield and vigor of vines were not influenced by the treatments, nor were the production of primary and secondary metabolites. However, the active ingredients influenced the mechanical properties of the skin (hardness and thickness). A significant hardening of the skin was detected when laminarin and chito-oligosaccharides were used, and sulfur induced a thickening of the skin with potential consequences for wine quality. Furthermore, the yeast community present on grape berries was influenced by the treatments. The abundance of Aureobasidium pullulans, the dominant species on the grape berry, changed in response to the compounds used. In addition, Alternaria sp. was reduced in some treatments with a potentially positive effect on the quality and the safety of the grapes. This study provides an overview of the effect of biocontrol products and resistance inducers on microbial ecology and “Nebbiolo” grape quality, contributing to the establishment of more sustainable and effective defense strategies in viticulture.

1 Department of Agricultural, Forest and Food Sciences, University of Turin, Grugliasco, Italy, 2 Agroinnova-Centre of Competence for the Innovation in the Agro-Environmental Sector, University of Turin, Grugliasco, Italy, 3 ANT-NET srl, Turin, Italy, 4 Institute for Sustainable Plant Protection, National Research Council (IPSP-CNR), Turin, Italy Viticulture is a cropping system in which treatment against fungal diseases (in particular powdery and downy mildews) can be extremely frequent. Accordingly, a reduction in antimicrobial treatments and the application of environmentally-friendly compounds are becoming increasingly important for a more sustainable viticulture. In addition to their effect against pathogens, the impact of these products on the quality of the grapes is very important for the oenological industries, but unfortunately at present few data are available. We evaluated the effect of the application of biocontrol products and resistance inducers in the vineyard on the mechanical properties, microbial ecology, technological and phenolic maturity of Vitis vinifera "Nebbiolo" grapes at harvest. The yield and vigor of vines were not influenced by the treatments, nor were the production of primary and secondary metabolites. However, the active ingredients influenced the mechanical properties of the skin (hardness and thickness). A significant hardening of the skin was detected when laminarin and chito-oligosaccharides were used, and sulfur induced a thickening of the skin with potential consequences for wine quality. Furthermore, the yeast community present on grape berries was influenced by the treatments. The abundance of Aureobasidium pullulans, the dominant species on the grape berry, changed in response to the compounds used. In addition, Alternaria sp. was reduced in some treatments with a potentially positive effect on the quality and the safety of the grapes. This study provides an overview of the effect of biocontrol products and resistance inducers on microbial ecology and "Nebbiolo" grape quality, contributing to the establishment of more sustainable and effective defense strategies in viticulture.

INTRODUCTION
Powdery and downy mildews are among the key fungal diseases of European grape-growing regions, including Italy, significantly impacting crop quality and quantity. Powdery mildew is caused by Erysiphe necator, a polycyclic fungus infecting all green tissues with gray-white mycelia on their surface. E. necator can cause a reduction in grape yield of up to 45%, as a consequence of smaller diseased berries (Calonnec et al., 2004). In addition, severely infected clusters have higher total acidity than healthy ones, a higher concentration of phenylacetic acid, acetic acid and γ-decalactone, and a lower content of total anthocyanin and 3-mercaptohexanol (Calonnec et al., 2004;Pinar et al., 2017). Downy mildew is caused by the oomycete Plasmopara viticola, affecting the leaves and berries under warm, moist and humid environmental conditions. P. viticola may reduce grape yield and affects the aroma and flavors in wine, with organoleptic impact of the pathogen infection differing according to specific grape variety and impacts on berry quality and on the ability of a wine to age (Pons et al., 2018). These diseases are mainly controlled using chemical fungicides in conventional vineyards, where problems with pesticide residues and resistance to active ingredients are increasing. Copper and sulfur are used on organic farms (Gessler et al., 2011;Gadoury et al., 2012).
Viticulture is one of the most-frequently treated cropping systems (Komarek et al., 2010;Gill and Garg, 2014;Tsakirakis et al., 2014;Provost and Pedneault, 2016). Consequently, restrictions have been applied by the European Commission on the number of pesticide treatments (Directive 2009/128/EC) and on the maximum quantity per year of copper fungicides (Regulation 2002/473/EC). Recently, copper, which is still considered the most important antifungal product in organic viticulture, has been added to the list of candidates for substitution (European Commission Implementing Regulation 2018/84) and further limited to 4 kg per hectare/year spread over 7 years from February 1st 2019 (European Commission Implementing Regulation 2018/1981 of 13 December 2018). In parallel, there is increasing demand by consumers for organic and more environmentally-friendly products and this trend is affecting the market of the viticulture and wine sectors (Schäufele and Hamm, 2017). This situation has encouraged the development of alternative crop protection strategies and of the use of more environmentally-friendly products, such as biofungicides, resistance inducers and biostimulants against powdery and downy mildews (Gadoury et al., 2012;Thind, 2017;Gutiérrez-Gamboa et al., 2019), with a view to achieving sustainable production (Flores, 2018). As a result, an increasing number of farmers have begun to apply only biocontrol products to grapevines, considering that they can be applied up to a few days before harvest and during full bloom. Biocontrol products are plant protection products containing micro-organisms (fungi, bacteria, viruses, etc.) or natural substances of mineral, plant or microbial origin, with a lower environmental impact compared to synthetic chemical products. Examples of these products are laminarin, chitosan and bicarbonate. In particular, laminarin is extracted from the brown algae Laminaria digitata, is known to induce resistance against downy and powdery mildew and gray mold in grapes (Aziz et al., 2003;Gauthier et al., 2014;Chalal et al., 2015;Romanazzi et al., 2016;Pugliese et al., 2018). According to Trouvelot et al. (2014), chitosan has been studied as plant protectant for more than 30 years and chitosan and chito-oligosaccharides have also been applied to control downy mildew on detached leaves Trotel-Aziz et al., 2006), with a variable control efficacy, and a decrease in the concentration of several amino acids has been observed on treated grapes (Garde-Cerdàn et al., 2017). Salts such as bicarbonates have been used successfully for the control of powdery and downy mildews (Deliopoulos et al., 2010;Lukas et al., 2016).
Due to a lack of information in the scientific literature, agronomists and wine growers are often not aware of the impact that these new products can have indirectly on the quality of grapes. In particular, few data are available on their influence and/or interaction in the synthesis of primary and secondary metabolites, on the composition of the skin cell wall and related mechanical properties, and on the microbiota present on the grape clusters.
The above-mentioned induced resistance can represent an attractive strategy for increasing the production of plant defense compounds in grapevine, such as phenolic compounds (Flamini et al., 2013). For red wine production, especially, the content and chemical characteristic of phenolic compounds in berry skins and seeds at harvest play an essential role in the sensory properties of the resulting product, impacting on the chromatic characteristics, astringency and bitterness (Preys et al., 2006;Paissoni et al., 2018). Moreover, the final phenolic composition of wine is related to the capacity of the grape skin and seeds to release these compounds during the maceration-fermentation process. Regarding berry skins, this phenomenon, conventionally called "extractability, " may be influenced by phytosanitary treatments affecting physical-mechanical characteristics (i.e., hardness and thickness) (Rolle et al., 2012a).
The grape surface is colonized by different microorganisms, mainly filamentous fungi, bacteria and yeasts, which in relation to the final product may have positive or negative effects (Milanović et al., 2013;Alessandria et al., 2015). Several environmental factors, such as geographical location, climatic conditions (temperature, rainfall, wind) and vineyard treatments, affect the microbial communities present on the grape surface at harvest (Barata et al., 2012;Bokulich et al., 2014;Zhang et al., 2019). Many of those factors may also interact, impacting on the microbial diversity and, consequently, wine quality (Barata et al., 2012). The importance of the grape microbiota, especially in the first stages of the spontaneous or inoculated alcoholic fermentation, is well documented (Cravero et al., 2016). In particular, non-inoculated fermentation, a practice that winemakers regard as enhancing regional typicity, is also popularly considered to enhance wine quality through the production of a wide spectrum of secondary compounds (Ciani et al., 2010). Considering the importance of the indigenous must microbiota on wine composition, there is an increasing interest in understanding the effect of vineyard plant protection products on species distribution and persistence on the grapes at harvest. To this end, studies have investigated the impact of organic and conventional agrochemicals on grape microbiota in a wide range of grape varieties (Schmid et al., 2011;Piao et al., 2015;Agarbati et al., 2019). Despite these extensive researches, little is known about the impact of a wide spectrum of biocontrol products and of resistance inducers on microbial diversity at harvest.
In the present paper, alternative strategies against downy and powdery mildews, based on the application of biocontrol products and resistance inducers, were applied in a vineyard of "Nebbiolo, " a wine grape variety commonly used in the production of renowned Barolo, Barbaresco and Roero wines. Our main purpose was to evaluate the effect of these products on (i) the parameters involved in the technological and phenolic maturity of the grapes at harvest, and (ii) the microbial ecology, as determined by culture-dependent and cultureindependent approaches.

Field Trial
The field trial was conducted on wine grape "Nebbiolo" (Vitis vinifera L.) in a vineyard located in Piobesi D'Alba (North-West Italy) (GPS: 44.731760, 7.988324, hill area) during 2017. Vines were planted in parallel and contiguous rows, vertically trained, Guyot pruned and all grafted onto the same rootstock (Kober 5BB). The distance between vines was 0.90 × 2.5 m. Vineyard management during the growing season was uniform in the experimental site and in compliance with regional best agricultural practices. The experimental scheme included four randomized blocks per treatment, each containing eight plants along a row, and the treatments were assigned to plots using a random-number generator. Best agricultural practices were considered in order to avoid drift of product from one plot to another, such as avoiding carrying out the treatments if wind speed was higher than 2 km/h and the direction was perpendicular to the rows. The experimental products being tested were applied with a hand-pulled 2-stroke engine sprayer at a pressure of 15 bar, distributing 400-600 l/ha. The plants were treated in accordance to the guidelines EPPO/OEPP PP 1/31 (3) (EPPO, 2019).
The commercial antifungal products were used according to the manufacturer's instructions. The products were compared with a non-treated control (CTR), a standard control strategy (S+Met) consisting of sulfur and metiram (until the end of flowering) and sulfur and copper hydroxide (until the end of the season) applied in conventional farms in the region, and a standard control strategy (S) applied in organic farms in the region.

Agronomic Parameters and Grape Sampling
At harvest, total fruit yield (kg) and number of clusters were determined for each treatment and block. Mean cluster weight was calculated. Plant vigor was estimated in winter as pruned wood weight.
Berry sampling was performed in the vineyard at harvest. For each treatment, 200 grape berries (50 per block) with pedicel attached were randomly selected and aseptically collected into sterile sampling bags. The samples were transported in the laboratory under refrigerated conditions and analyzed immediately.

Berry Skin Phenolic Composition
The evaluation of the berry skin phenolic composition (i.e., anthocyanin, oligomeric and polymeric flavonol contents) was conducted for each treatment and block, following the Di Stefano andCravero (1991) method with Río Segade et al. (2014) modifications. Ten berries for each sample were weighed and peeled, and the skins obtained were quickly placed in 40 mL of a hydroalcoholic buffer solution prepared as follows: 12% (v/v) ethanol, 5 g/L tartaric acid, 2 g/L sodium metabisulfite, and afterward the pH was adjusted to pH 3.2. Then, the skins were macerated for 24 h, and afterward homogenized using an Ultra-Turrax T25 homogenizer (IKA, Staufen, Germany) set at 8000 rpm for 1 min. After centrifugation (3000 × g, 20 • C, 15 min), the supernatant was taken and brought to 50 mL with the Copper hydroxide 600 7-9 From fruit set to harvest 8 S+Met, standard control strategy applied in conventional farms; S, standard control strategy applied in organic farms in the region. *Applications started on April 24th 2017 and ended on August 25th 2017. Bunch closure on June 7th 2017, fruit set on June 16th 2017.
above-mentioned buffer solution and used as grape skin extract in subsequent analysis. The phenolic indices were determined by spectrophotometry, using a Shimadzu UV-1800 machine (Shimadzu Corporation, Kyoto, Japan). The total anthocyanins and flavonoids content was determined after dilution of the grape skin extract using a 70:30:1 (v/v) ethanol:water:hydrochloric acid (concentrated) solution. Proanthocyanidins assay required acid-catalysis of the polymerized flavanols under heat. Vanillin assay was performed by binding the vanillin reagent to the free 6 and 8 positions of monomeric-oligomeric flavanols in acid medium.

Grape Texture Analysis
Ten intact berries were randomly taken at harvest from each treatment and block, thus resulting in 40 berries for each treatment, and used for the determination of their instrumental mechanical properties. The skin parameters detected were berry skin break force (F sk ) and berry skin thickness (Sp sk ), following the method proposed by Letaief et al. (2008). Briefly, for the first test each intact berry was placed alone on the platform of a TA.XTplus Texture Analyser (Stable Micro Systems, Godalming, Surrey, United Kingdom) and was subjected to penetration by a standardized stainless steel needle (P/2N probe, Stable Micro System) at a constant speed of 1 mm/s. The maximum resistance of the berry skin to the needle penetration is defined as the berry's F sk , expressed in N. For the second test, Sp sk , a portion of the skin was obtained from each berry by peeling and removing pulp residues. The berry skin portion was placed on the flat platform of the above-mentioned instrument, and a probe composed of a flat 2 mm stainless steel cylinder (P/2 probe, Stable Micro Systems) began the descent at 0.2 mm/s toward the sample. In this way, the machine acted as a precise caliper and was able to detect the thickness of the berry portion (Sp sk ), expressed in µm.

Microbiological Analysis of Grapes
For each treatment-block combination, 50 g of single grape berries were placed in sterile bottles containing 450 mL of a solution of 0.1% peptone and 0.001% Tween 80 (Sigma, Milan, Italy). The bottles were stirred for 30 min at 150 rpm to facilitate the release of cells from the berry surface of the grapes in the solution. The rinse solution was serially diluted in peptone water (0.1% peptone). Aliquots of the dilutions were plated as follows. Wallerstein Laboratory Nutrient Medium agar (WLN, Biogenetics, Milan, Italy) was used for enumeration of fungi. The medium was spread plated and incubated at 30 • C for 5 days. MRS agar (Biogenetics) supplemented with Delvocid (25 µg/mL, DSM Specialties, Heerlen, Netherlands) was used for enumeration of lactic acid bacteria and was pour plated and incubated at 30 • C for 1 week. GYC [10 g/L yeast extract (Biogenetics), 20 g/L CaCO 3 (Sigma), 20 g/L agar, 20 mL/L ethanol (99% v/v, Sigma)] was used for enumeration of acetic acid bacteria and spread plated and incubated at 30 • C for 1 week. After incubation, the colonies were counted and mean and standard deviation values were calculated. Yeast colonies formed on WLN plates were counted on the basis of the color and morphology as described previously by Urso et al. (2008). Ten yeast isolates from each colony morphotype were picked, purified by streaking on WLN medium and subjected to molecular identification to species level.

Molecular Identification of Yeast Isolates
Isolates randomly selected from the WLN plates were purified by streaking on WLN and then grown overnight in 1 mL YPD medium (2% w/v yeast extract, 1% w/v peptone, 2% w/v dextrose, all from Biogenetics) at 30 • C. The broth culture was centrifuged at 14,000 × g for 5 min and the pellet was subjected to DNA extraction according to Cocolin et al. (2000). The extracted DNA was quantified using a Nanodrop ND-1000 spectrophotometer (Celbio, Milan, Italy) and standardized at 100 ng/µL. The ITS1-5.8S rRNA-ITS2 (ITS) region was PCR amplified and subjected to restriction fragment length polymorphism (RFLP) for all isolates (Esteve-Zarzoso et al., 1999). The ITS region was amplified with primers ITS1 (5 -TCCGTAGGTGAACCTGCGG -3 ) and ITS4 (5 -TCCTCCGCTTATTGATATGC 3 ) (White et al., 1990). The reaction mix was 50 µL and contained 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl 2 , 0.2 mM of deoxynucleoside triphosphates (dNTPs), 1.25 U of Taq Polymerase (Applied Biosystems, Milan, Italy), 0.2 µM of each primer and 100 ng of template DNA. Amplification was carried out using a PTC-200 DNA Engine MJ Research Thermal Cycler (Bio-Rad, Milan, Italy), as described by Esteve-Zarzoso et al. (1999), and the PCR products were checked by electrophoresis on 1.5% (w/v) agarose gel. The PCR products were subsequently digested with endonucleases HinfI, HaeIII, CfoI (Promega, Milan, Italy), according to the supplier's instructions. The restriction fragments were separated by electrophoresis in 3% agarose gel and stained with ethidium bromide. PCR and RFLP fragment lengths were used to identify yeasts by comparing the restriction bands with those available in the literature (Guillamon et al., 1998;Esteve-Zarzoso et al., 1999;Granchi et al., 1999;Sabate et al., 2002). Identification to species level was confirmed by sequencing the D1-D2 loop of the 26S rRNA encoding gene, after amplification using primers NL1/NL4 (Kurtzman and Robnett, 1997) to obtain a PCR product, which was sequenced by a commercial facility (Eurofins, Germany).

Mycobiota Amplicon Target Sequencing
Ten mL of the grape rinse was centrifuged for 10 min at 14,000 × g and the supernatant was discarded. The pellet was used for total DNA extraction using the Master Pure Complete DNA and RNA Purification Kit (Illumina Inc, San Diego, CA, United States) following the manufacturer's instructions. DNA was quantified using the Qubit dsDNA assay kit (Thermo Fisher Scientific) and standardized at 5 ng/µL. DNA was used to assess the mycobiota by targeting the D1-D2 domain of the 26S rRNA gene using a modification of the primer LS2F (5 -GAG TCG AGT TGT TTG GGA AT-3 ) and NL4R (5 -GGT CCG TGT TTC AAG ACG G-3 ), as recently described (Mota-Gutierrez et al., 2019).
PCR mixtures (25 µL) were prepared using 12.5 µL of the 2X Kapa HiFi HotStart ReadyMix Taq (Roche, Milan, Italy), 1 µM of each primer, 2.5 µL of DNA template, and PCR-grade water. Thirty cycles of 30 s of denaturation (95 • C), 30 s of primer annealing (55 • C), and 30 s of primer elongation (72 • C), followed by a final elongation step (72 • C) of 10 min, were performed. The amplicons' size and quality were assessed using a Bio-Rad Experion Workstation (Bio-Rad, Milan, Italy). In a second PCR step, amplicons were combined with the sequencing adapters and dual indices using the Nextera XT Index Kit (Illumina, San Diego, CA, United States), forming the multiplexed pairedend libraries, and quantified using the Qubit dsDNA Assay Kit (Thermo Fisher Scientific, Milan, Italy). Individual libraries were diluted to 4 nM, denaturated with 0.2 N NaOH and spiked with 30% (v/v) PhiX. The combined pool library and PhiX were diluted to 12 pM and paired-end sequencing was performed on a MiSeq platform, using MiSeq Reagent Kit v3 2 × 250 cycle (Illumina, San Diego, CA, United States), following the standard Illumina sequencing protocol. The software used for the base calling and Illumina barcode demultiplexing processes were the Miseq Control Software version. The 26S raw read data were deposited in the Sequence Read Archive of NCBI under BioProjectID PRJNA575160.

Wine Production and Analysis
Experimental wines were produced from three treatments: CTR, potassium bicarbonate (K-Bic), and S+Met. A standard vinification protocol was applied for each treatment considered and two replicates of 100 kg of grapes were processed in the experimental cellar of University of Turin under controlled conditions. In brief (Boccacci et al., 2020): grapes harvested were crushed in a TEMA destemmer-crusher (Enoveneta, Piazzola sul Brenta, Italy). The mash was added with 25 mg/L of potassium metabisulfite. After about 6 h, selected Saccharomyces cerevisiae yeasts (Lalvin BRL97, Lallemand, Inc., Montreal, Canada) were inoculated at a dose of 20 g/hL. Two punch-downs per day were carried out in the first 3 days, then two pumping-overs per day until the end of maceration, which lasted 10 days. The end of macerations was followed by the gentle pressing of the pomace cap using a PMA 4 pneumatic press (Velo SpA, Altivole, Italy) with a maximum pressure of 1.2 bar, and a small aliquot of the press wine was joined to the free-run wine. Wine was inoculated with Oenococcus oeni Lalvin VP41 strain (Lallemand) to induce malolactic fermentation (MLF). Once MLF was completed, the wines were racked to remove lees, and free SO 2 concentration was adjusted to 70 mg/L. The alcoholic fermentation (AF) and MLF were carried out at controlled temperatures of 27 ± 2 and 20 ± 1 • C, respectively. Wines were cold-stabilized at 0 • C for 2 weeks, filtered (Seitz K300 grade filter sheets, Pall Corporation, Port Washington, NY, United States) and then bottled in glass bottles of 0.75 L with cork stoppers.
After 5 months of bottle storage, wines were analyzed. Following parameters were determined: ethanol content, residual sugars, glycerol, and organic acids by HPLC using the chromatographic condition proposed by Giordano et al. (2009), while pH and titratable acidity were evaluated using OIV (2017) methods. Wine phenolic composition was determined through spectrophotometric methods using an UV-1800 spectrophotometer (Shimazdu Corp., Kyoto, Japan): total polyphenolic index was assessed by Folin-Ciocalteu reagent assay, proanthocyanidins determined after acid hydrolysis at 100 • C using a ferrous salt (FeSO 4 ) as catalyst according to the Bate-Smith reaction and expressed as mg of cyanidin chloride/L of wine, monomeric and oligomeric forms of flavanols were evaluated as flavanols reactive to vanillin and expressed as mg of (+)-catechin/L of wine (Petrozziello et al., 2018). Total anthocyanins were determined by measuring absorbance at 536-540 nm after dilution with a hydroalcoholic solution composed of ethanol: water: 37% hydrochloric acid (70:30:1, v/v) and expressed as mg of malvidin-3-glucoside chloride/L of wine (Petrozziello et al., 2018). After the acquisition of visible spectra of undiluted samples using 1-mm optical path cuvettes, color intensity was calculated as the sum of absorbance measured at 420, 520, and 620 nm (A420 + A520+ A620 on an optical path of 10 mm) and hue was obtained as the ratio of absorbances measured at 420 and 520 nm (A420/A520) following the method OIV-MA-AS2-07B (OIV, 2017). The wine color was also evaluated by the CIE L * a * b * parameters, namely lightness (L * ), red/green color coordinate (a * ), and yellow/blue color coordinate (b * ), according to the method OIV-MA-AS2-11. The total color difference ( E * ) between two samples (for example between a fined wine and the respective control wine) was calculated using the following expression:
OTUs were clustered at 99% similarity by means of UCLUST clustering methods (Edgar, 2010) and representative sequences of each cluster were used to assign taxonomy to a reference 26S database by means of the RDPII classifier (Wang et al., 2007). An ad hoc database (Mota-Gutierrez et al., 2019) was created with the reference 26S gene sequences and the taxonomy. This database was used in the QIIME pipeline to obtain the OTU table and the downstream output. Sequences were double-checked using the BlastN search tool 1 to confirm the taxonomy assignment. In order to avoid biases due to different sequencing depths, all samples were rarefied at 16 564 reads after raw read quality filtering, and the OTU table was filtered for OTUs occurring at 0.2% of the relative abundance in at least two samples. Statistics and plotting were carried out in the R environment 2 . Alpha diversity indexes (Shannon, chao1 and n • of observed species) were calculated by the function vegan in R (Dixon, 2003) and was further analyzed using the Pairwise Kruskal-Wallis test in order to find significant differences in microbial taxon abundance according to treatment. The OTU table was used to perform the Adonis statistical test with 999 permutations through the vegan package of R. The OTU table displays the highest taxonomy resolution that was reached; when the taxonomy assignment was not able to reach species level, the genus or family name was displayed. The OTU table was used to build a principal-component analysis (PCA) as a function of treatment by using the made4 package of R.
Kruskal-Wallis tests were used to find significant differences in microbial taxa abundance according to treatment. Pairwise Spearman's correlations were used to study the relationships between the relative abundance of mycobiota with disease index as well as with the berry skin mechanical properties. The correlation plots were visualized in R using the corrplot package of R. P values were adjusted for multiple testing using the Benjamini-Hochberg Procedure, which assesses the false discovery rate (FDR).
To test for significant differences in agronomic parameters, berry mechanical properties and must chemical compositions, the analysis of variance (ANOVA) was used (SAS statistical software, version 8.2, SAS Institute, Cary, NC, United States). Significant differences were highlighted at p < 0.05.

Disease Incidence
The "Nebbiolo" vineyard used in the study was not infected with P. viticola and no symptoms caused by downy mildew were observed. Conversely, 58% of clusters were affected by E. necator, with a severity of 8% on non-treated control (CTR). Comprehensive details of the phytopathological results of this trial are reported in Figure 1. Potassium bicarbonate (K-Bic) was the most efficient in reducing the disease, similar to the standard treatments (S+Met and S). In addition, treatments using acibenzolar-S-methyl (AcS-Met), potassium phosphonate (K-Pho), laminarin (Lam), chito-oligosaccharides and oligogalacturonides (Chito) and calcium (CaO) were also effective compared to the standard treatment with sulphur (S). Fosetyl-Al (Fos-Al) and metiram (Met) had a lower efficacy and similar to Chito. Electrolyzed water (EOW) was not effective.

Agronomic Parameters and Chemical Analysis of Grape Must
The data on yield components collected at harvest, as well as the vigor estimation at the end of the vegetative season, did not produce significant differences related to the various treatments applied ( Table 2). The incidence of fungal pathogens on leaves and berries reported above did not seem to have substantial effects, compared with the CTR, on the quantitative performance of the vines in the particular season in which the trial was conducted.
In addition, the chemical parameters, which define the socalled technological maturity of the grapes, did not show significant differences among the different treatments (Table 3). Sugar content, pH and total acidity presented homogeneous values among the different treatments, with data in agreement with previous studies (Cagnasso et al., 2011). Therefore, for the primary metabolites also, the treatments used to control downy and powdery mildews did not have a substantial interference or elicitor effect. Other studies have also highlighted that grape musts from Montepulciano grapevines treated with laminarin and copper hydroxide were not significantly different from control in terms of sugar content, pH and titratable acidity (Romanazzi et al., 2016). FIGURE 1 | Efficacy of the treatments against Erysiphe necator on the incidence and severity on grape clusters, compared to non-treated control (CTR) and standard treatments for conventional (S+Met) and organic (S) farms in the region. Different letters mean significant differences (p < 0.05). Treatment descriptions are reported in Table 1.

Phenolic Composition of Berry Skin
The effects of the treatments on the grape skin phenolic composition are reported in Table 4, using the main spectrophotometric indexes used by operators of the viticultureoenology sector. Due to high variability of data among replicates, no statistical differences in anthocyanin, polymeric flavanol (proanthocyanidins) and oligomeric flavanol (reactive to vanillin) contents were found among the different treatments. However, low average values of total anthocyanin content were detected using Calcium oxide (CaO), while high total flavanols content (polymeric + oligomeric) was detected after the application of chito-oligosaccharides and oligogalacturonides (Chito). Despite these small differences, it is possible to hypothesize that in "Nebbiolo" there is no real impact from the different treatments in the biosynthetic pathways of these secondary metabolites, as well as no elicitor effect.

Mechanical Properties of Grapes
The effect of treatments on berry skin hardness, evaluated by skin break force parameter (F sk ), and skin thickness (Sp sk ) are reported in Table 3. Significant differences for F sk were not found for any treatment when compared to CTR. Lam and Chito showed the significant highest values of F sk (0.627 and 0.629 N, respectively) only with respect to K-Bic, whose grapes  Data expressed as average ± standard deviation (n = 4 for grape must parameters, n = 40 for berry skin mechanical parameters). Different letters inside each column means significant differences at p < 0.05 (Tukey-b test). Treatment descriptions are reported in Table 1. are characterized by lowest values of F sk (0.547 N). K-Bic was one the best treatments against powdery mildew (Figure 1) with efficiency similar to standard chemical strategy S+Met, and more effective than Lam and Chito. The F sk values found agree with data previously published for this variety (Rolle et al., 2012b). From a technological point of view, higher F sk produces an increase in anthocyanin extraction yield and slower extraction kinetics during simulated maceration of "Nebbiolo" grape skin in wine-like solution . The last aspect is particularly important for wine grape varieties rich in 3 -hydroxylated anthocyanins, such as "Nebbiolo, " because these pigments are extracted preferentially during the initial phase of maceration and may be easily oxidized by the enzymes present in the juice (González-Neves et al., 2008). Despite the high standard deviation of data within each treatment, significant differences were found in berry Sp sk with respect to CTR. The values of the Sp sk parameter ranged from 176 to 209 µm for treatments CTR and S, respectively. Although all treatments impacted on the value of this mechanical property when compared to data for CTR berries, only after the application of Fos-Al, EOW, K-Bic, S+Met, and S a significant increase in skin thickness was observed (+22 -33 µm). Interestingly, the three most effective treatments against powdery mildew (S, S+Met and K-Bic) have shown higher values of skin thickness and this could suggest in "Nebbiolo" a correlation between this parameter and the resistance to E. nector. However, the same higher Sp sk values were detected in EOW, the worst treatment for the control of powdery mildew without significant differences to the CTR, therefore further studies are needed. In general, even small variations in thickness can cause differences in terms of the extractability of phenolic substances, as reported by Río Segade et al. (2011) on "Galician" grapes. However, for "Nebbiolo" grapes treated with a specific inactive LSA extract in the vineyard, the berry skin thickening observed did not have a detrimental impact on anthocyanin extractability (Giacosa et al., 2019). Furthermore, skin thickness and hardness are factors that contribute to berries' response to climatic adversities, regulate on-vine and off-vine withering kinetics  and limit damage during harvest, packing, transport and storage (Kök and Çelik, 2004).

Microbiological Quality of Grapes
Three microbial groups on the grape surface were analyzed by viable count determination. Lactic Acid Bacteria (LAB, on MRS agar) were determined to be present at a concentration <10 colony forming units (CFU)/g, while Acetic Acid Bacteria (AAB, on GYC agar) were present at a concentration <100 CFU/g. The count of viable fungi is reported in Figure 2A. For all treatments, the load was between 3.5 and 4.3 Log 10 CFU/g, corresponding to values generally reported in literature for mature grapes (Prakitchaiwattana et al., 2004). All treatments had a significantly lower fungal load compared to the CTR, which was the only treatment that had a mean load above 4 Log 10 CFU/g. The lowest concentration (mean value 3.46 Log 10 CFU/g) was registered for Chito and K-Bic. Aside from CTR, the highest fungal count was registered for Fos-Al (3.98 Log 10 CFU/g) and this value was significantly higher than the loads for Chito and K-Bic.
Ten yeast isolates of each colony morphotype on WLN medium were randomly selected from different samples and subjected to PCR-RFLP of the ITS1-5.8S-ITS2 region in order to obtain an identification to species level. The species identification was confirmed for selected isolates by sequencing of the D1/D2 loop of the 26S rRNA encoding gene. The fungal identification results are shown in Figure 2B. In total, four different nonfermenting yeast species were detected on the grapes by the culture-dependent analysis. Aureobasidium pullulans was isolated across all treatments and in the CTR. These data are in accordance with previous studies that demonstrated A. pullulans is the most abundant species on grape berries at harvest time (Barata et al., 2012). Although this species is irrelevant to the winemaking industry due to an inability to tolerate ethanol, these organisms represent the resident microbiota of wine grapes (Urso et al., 2008;Alessandria et al., 2015). The presence of this microorganism on grape berries was influenced by the treatments, the frequency of isolation of A. pullulans ranging from around 20% for K-Pho to 100% for K-Bic, S+Met and CaO. The presence of higher populations in K-Bic, S+Met and CaO may influence the overall microbiota present on grape berries due to the antagonism of A. pullulans toward other yeasts and molds (Bozoudi and Tsaltas, 2018). Interestingly, an A. pullulans-based biopesticide has been developed and commercialized to take advantage of its ability as a biocontrol agent against pathogens (Chi et al., 2009). Consequently, the complete dominance of A. pullulans on grapes obtained from K-Bic, S+Met and CaO could facilitate the control of grape pathogens, as previously reported by Comitini and Ciani (2008). Interestingly, K-Bic and S+Met were two of the most effective treatments against powdery mildew ( Table 1).
The second most common species was Cryptococcus flavus, which was identified in nine out of 11 treatments and in the CTR. This species was a minor population in all treatments except for K-Pho in which it was dominant, isolated with a frequency of almost 50% of the total colonies. Cryptococcus terrestris was the third species isolated in five out of 11 treatments and in CTR, but with a frequency always below 10%. Rhodotorula graminis was isolated only in K-Pho and in CTR: these last two treatments exhibited the highest level of biodiversity with the presence of four fungal species. Overall, these data suggest that the fungal communities that colonized "Nebbiolo" grapes were influenced by antifungal compounds. Thus, differences could be explained by the selective effect of copper and sulfur-based fungicides on the yeast ecosystem, as previously demonstrated by Grangeteau et al. (2017a) for A. pullulans, Hanseniaspora guilliermondii, Hanseniaspora uvarum, Metschnikowia sp., Pichia membranifaciens, Saccharomyces cerevisiae, and Starmerella bacillaris yeast isolates.

Fungal Communities on Grapes Determined by Direct Sequencing
An appreciation of the microbial ecosystems associated with grape berry surfaces is fundamental to understanding the health status and the microbes' contribution during alcoholic fermentation, with consequent potential impact on wine quality (Barata et al., 2012). The investigation of yeast biodiversity using culture-dependent methods can lead to incomplete microbial detection, due to low abundances or the viable but non-culturable (VBNC) yeast species (Alessandria et al., 2015). The development of next-generation sequencing (NGS) technology makes it possible to obtain more detailed information associated with wine grapes and subsequent alcoholic fermentation (Bokulich et al., 2014;Grangeteau et al., 2017b).
In this study, 5 064 209 raw reads (2 × 250 bp) were obtained and 4 573 288 reads passed the filters applied through QIIME, with an average value of 103 938 ± 46 172 reads/sample, and a median sequence length of 389 bp (Supplementary Table S1). The rarefaction analysis and the estimated sample coverage indicated that there was satisfactory coverage for all the samples (ESC average 97.36%). Moreover, Shannon index (Figure 3 and Supplementary Table S1) showed a higher level of complexity (p < 0.05) when comparing Chito treatment to EOW as well as Chito compared to Met. No significant difference in the alphadiversity was observed between the treatments (Figure 3 and Supplementary Table S1). Interestingly, a significant increase in the number of OTUs detected (p < 0.05) was observed in Chito when compared to S+Met and Met (Figure 3 and Supplementary Table S1).
Through PCA, based on the relative abundance of the OTUs (Supplementary Figure S1), a stratification of the samples according to treatment was observed. Excluding Fos-Al, which clustered together with CTR, the other treatments were well separated from the CTR (Anosim statistic p < 0.001). However, no clear difference between the treatments could be observed.
At harvest time 24 fungal species were identified by NGS, while only four yeast species were identified by culturedependent methods. Regarding the most abundant OTUs shared in the entire dataset (Table 5), it was possible to observe a FIGURE 2 | Fungal load, as determined by microbiological analysis on WLN medium, expressed in Log 10 CFU/mL. The bars represent mean values and standard deviation from four replicates for each treatment. Different letters mean significant differences (p < 0.05) (A). Fungal species biodiversity: for each treatment, the frequency of isolation of the four species detected is expressed as a percentage. Isolates grown on WLN medium were grouped based on colony morphology and representatives of each colony type were identified to species level by molecular methods (B). Treatment descriptions are reported in Table 1. common mycobiota shared among treatments. The main OTUs detected were Aureobasidium (ranging from 24 to 38% of the relative abundance), Cladosporium cladosporioides (ranging from 21 to 29%), Alternaria tenuissima (ranging from 10 to 17%), Pseudeurotiaceae (ranging from 6 to 18%), Didymella negriana (ranging from 2 to 7%), Saccharomyces cerevisiae (ranging from 0.5 to 8%) and Botrytis cinerea (ranging from 0.2 to 13%). Regarding S. cerevisiae, Fleet et al. (2002) reported the low presence (populations of less than 10 CFU/g) of this species found on wine grapes using culture-dependent FIGURE 3 | Boxplots to describe α-diversity measures of grape mycobiota. Individual points and brackets represent the richness estimate and the theoretical standard error range, respectively. Treatment descriptions are reported in Table 1. techniques and proposed the use of enrichment culture instead of direct agar plating. The NGS results showed S. cerevisiae present in all the different grapes, indicating the importance of this technique to get a better view into microbial ecosystems (Bokulich et al., 2014).
The employment of different fungicides and resistance inducers did not cause a dramatic change in the yeast genera detected by culture-independent technique. According to statistical tests, only the minor OTU fractions (ranging from 0.09 to 12%) were affected by the different treatments ( Table 5). More specifically, Alternaria spp. were reduced by S and Met, compared to the CTR (FDR < 0.05). It should be pointed out that Alternaria is recognized as being among the main mycobiota populations of grape at harvest and are considered to be indicators of grape spoilage (Tournas and Katsouda, 2005). In addition, Alternaria is known to produce secondary metabolites (mycotoxins) that can contaminate first grapes and later the wine (Serra et al., 2006). Furthermore, Didymellaceae was reduced by S+Met, EOW and Met; Heterophoma novae-verbascicola by the use of S+Met and Met; Vishniacozyma carnescens and V. victoriae by the use of Chito, K-Bic, S and Met, while Filobasidium magnus was affected by Chito and S (FDR < 0.05) compared to the CTR. Nevertheless, the relevance for the microbiological quality of the grapes of the differences described above is not expected to be important due to the low abundance of these populations. Interestingly, it was possible to detect by a culture-independent approach E. necator in grape berries of the CTR and all the treatments applied. It should be noted, that E. necator was the only population that could be considered of relevance due to the high relative abundance obtained in all the samples (range 8-20%, average 15%). Regarding the effect of the various treatments on the population of E. necator, we observed a reduction when EOW was used only if compared with Chito and Lam (Figure 4).
The mycobiota composition was also tested for possible correlations with the disease index and berry skin mechanical properties obtained for each of the treatments as well as the CTR. As shown in Supplementary Figure S2, we observed a positive correlation between Cladosporium cladosporioides and berry skin break force (FDR < 0.05). C. cladosporioides was the second most abundant population in samples from all treatments and the CTR. However, no significant differences were observed in the relative abundance between the various treatments. This result implies that this species was not influenced by the treatment applied in the vineyard. The presence of Cladosporium genus, in particular Cladosporium herbarum and Cladosporium cladosporioides, on grapes is associated with the development of Cladosporium rot (Barata et al., 2012). Recently, it has been shown that the population of Cladosporium increases as grapes Only OTUs with an incidence above 0.2% of the relative abundance in at least two samples are shown. Abundance of OTUs in the biological replicates was averaged.
mature and it may develop on the surface of apparently healthy grapes (Briceño and Latorre, 2008). The positive correlation observed, i.e., increased grape skin force associated with the presence of C. cladosporioides, would require further investigation in order to confirm this correlation for other grape varieties and in subsequent harvests, and more importantly to establish the cause-effect relationship of the observations made in this study. Also the presence of Heterophoma novae-verbascicola was correlated with berry skin break force (FDR < 0.05). Furthermore, a direct association of % cluster and berries affected by powdery mildew with Filobasidium stepposum, Vishniacozyma victoriae and Vishniacozyma carnescens was observed. In the case of these last correlations revealed by statistical analysis, the OTUs involved were present at a very low abundance (<1%), are considered minor and therefore it would be important to confirm this association using RNA as a target as it is a better indicator of vitality compared to DNA, and if possible to obtain quantitative data by performing viable counts.

Wine Chemical-Physical Characteristics
The effects on wine of K-Bic (the most effective biocontrol product in reducing powdery mildew) were evaluated in comparison with wine produced from CTR and the standard control strategy (S+Met). Significant differences (p < 0.05) were found among treatments for almost all parameters tested (Supplementary Table S2). The ethanol content well reflected the average values found for soluble solids on the grape must, while the differences in acidity parameters were not always in line with grape berry measurements. All wines performed the malo-lactic fermentation, leading to no residual malic acid detected in the wines, and the tartaric acid stabilization. Among acidity values, it is worth noting the acetic acid trend, highlighting reduced contents for K-Bic and particularly for S+Met treatment, resulting in significant differences (p < 0.001) for all wines produced (Supplementary Table S2).
Despite the reduced ethanol content, the S+Met trial achieved the highest quantities of anthocyanins and total phenols (p < 0.01), but these differences were not reflected in the flavanols content, analyzed by both high-mass tannins (proanthocyanidins assay; p > 0.05) and low-mass compounds (vanillin assay; p > 0.05). The differences in the anthocyanin richness, the principal compounds responsible for the red color of young wines, had a decisive impact on the color intensity and on the CIEL * a * b * a * coordinate (p < 0.01). Both K-Bic and S+Met treatments resulted in a wine color with a E * above 3 when compared to CTR, a threshold indicating the potential ability of the human eye to recognize the color differences during a tasting (Martínez et al., 2001).  Table 1. CONCLUSION This study was focused on possible indirect effects of vineyard treatments with chemical, biocontrol products and resistance inducers against downy and powdery mildews on "Nebbiolo" grapes at harvest. Treatment with potassium bicarbonate (K-Bic) was the most effective in reducing powdery mildew on grape clusters with an efficiency similar to the standard control strategies S+Met and S applied, respectively, in conventional and organic farms.
The yield and vigor of vines, as well as the primary and secondary metabolites produced in the berries were not influenced by the different treatments. On the contrary, the active ingredients present in the treatments modified the mechanical properties of the skin (hardness and thickness). In particular, a slight hardening of the skin was detected when laminarin (Lam) and chito-oligosaccharides and oligogalacturonides (Chito) were used, while treatments with sulfur (S and S+Met) and K-Bic induced a thickening of the skin. Interestingly, these treatments were the most effective treatment against powdery mildew. In addition to mechanical properties, the total yeast community present on grape berries was also modified in response to antifungal treatments. A. pullulans was the dominant species and its relative abundance was influenced by the compounds used in this study. In particular, A. pullulans is also known as biocontrol agent against pathogens and a higher population was detected in berries treated with K-Bic and S+Met, products particularly effective against powdery mildew. Alternaria was reduced by treatments and this can positively affect the quality and the safety of the grapes. The effect of such treatments on population kinetics during spontaneous or inoculated alcoholic fermentation deserves further investigation.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the Sequence Read Archive of NCBI, BioProjectID PRJNA575160.

AUTHOR CONTRIBUTIONS
KR, SG, and MP performed most of the experiments, analyzed and compared the data. KR, VE, and IF performed the microbiological analysis and bioinformatic elaboration of data. SG and SR carried out the chemical and texture analyses.
MP and MM performed the pathological analyses. IG carried out the agronomical measurements, analyzed the corresponding data. MG supported the pathological analyses and critically revised the manuscript. GG, LR, MP, and KR conceived the project, supervised all the experiments and wrote the manuscript. All authors read and approved the manuscript.

FUNDING
This research was financially supported by the SAFEGRAPE project (Sustainable grapevine protection against fungal pathogens) funded by Fondazione Cassa di Risparmio di Cuneo.

ACKNOWLEDGMENTS
Authors gratefully acknowledge the contribution of Tenuta Carretta for vineyard management and supplying grapes. We thank "Vignaioli Piemontesi" and Dr. Danila Cuozzo for technical assistance in the collection of agronomic data.