Identification of the Causal Agent of Aqueous Spot Disease of Sweet Cherries (Prunus avium L.) from the Jerte Valley (Cáceres, Spain)

The pre and postharvest disease named ‘aqueous spot’ is an emerging risk for sweet cherries growing in Jerte Valley (Cáceres, Spain). Early stages of the disease appear in the tree, but it is usually detected after harvesting, during the postharvest period. Symptoms include the appearance of skin discolouration and translucency in the shoulder areas. At the most advanced stages, a mycelium of white colour partially or completely covers the fruit. This manuscript provides a detailed description of the microbes involved in this disease, such as bacteria, yeasts, and moulds. Microbes of different cherry cultivars were studied during two consecutive seasons (2019 and 2020). The counts of bacteria and yeast in damaged tissues were higher (7.05 and 6.38 log10 CFU/g for total aerobic mesophilic microbes and yeasts, respectively) than sound tissues (6.08 and 5.19 log10 CFU/g, respectively). The Enterobacterales order dominated the bacteria population. Among yeasts, Yarrowia lipolytica, in 2019, and Metschnikowia pulcherrima and Metschnikowia viticola, in 2020, were consistently isolated from all samples. The presence of moulds was inconsistently detected at the early stage of this disease by plate counts. However, microscopic observations revealed the presence of hyphae in cherry flesh. Different pathogenic moulds were identified, although white mycelium, identified as Botrytis cinerea by molecular methods, was consistently isolated at later stages. Inoculation tests confirmed the involvement of white-mycelium B. cinerea in the development of this new postharvest disease in the Jerte Valley. Its combination with Enterobacterales enhanced the evolution of rotting, whereas the combination with yeasts decreased and delayed the symptoms. This work presents the first report of a consortia of microorganisms implicated in the development of ‘aqueous spot’, an emerging disease in sweet cherry cultivars in the Jerte Valley.


Introduction
The Jerte Valley (Extremadura, Southwest of Spain) is the most important sweet cherry (Prunus avium L.) production area in Spain. This region has an area of 7524 ha, with a potential yield of more than 41,000 tons [1]. The production is mainly exported to European countries, such as the United Kingdom, Germany, France, and Italy. Besides the autochthonous 'Picota' type cherries [2][3][4], the catalogue of cultivars in the Jerte Valley is extensive, producing from the end of April to the beginning of August. Among the early The productive cherry sector in the Jerte Valley needs to control the emerging disease called 'aqueous spot'. It is important to determine the possible origin and factors that favour its occurrence. Thus, the aim of this work was to elucidate the causal agent of this pathology. For that, an extensive analysis of the microbial population present during the development of 'aqueous spot' was performed. Finally, inoculation trials confirmed the specific identification of the microbes involved in this disease.

Sample Collection
All cherry samples, at commercial maturity stage 4 according Serradilla et al. [25], were provided by the 'Agrupación de Cooperativas del Valle del Jerte SC.L.' (https://www.ac-vallejerte.es/, accessed on 25 July 2021) during the 2019 and 2020 seasons. They were transported under refrigerated conditions (<5 °C) to the Agricultural Engineering School facilities from the University of Extremadura (Badajoz, Spain) and processed in less than 24 h.

Sample Preparation, Microbiological Counts and Isolation of Microbes
To elucidate the causal agent of 'aqueous spot' disease, 10 cherries with initial symptoms by batch were selected for microbial analysis. Initial symptoms were characterised

Sample Collection
All cherry samples, at commercial maturity stage 4 according Serradilla et al. [25], were provided by the 'Agrupación de Cooperativas del Valle del Jerte SC.L.' (https:// www.ac-vallejerte.es/, accessed on 25 July 2021) during the 2019 and 2020 seasons. They were transported under refrigerated conditions (<5 • C) to the Agricultural Engineering School facilities from the University of Extremadura (Badajoz, Spain) and processed in less than 24 h.

Sample Preparation, Microbiological Counts and Isolation of Microbes
To elucidate the causal agent of 'aqueous spot' disease, 10 cherries with initial symptoms by batch were selected for microbial analysis. Initial symptoms were characterised by small translucent areas with a loss of firmness to the touch ( Figure 1A,B). Damaged tissue and sound tissue of cherries were analysed separately. The damaged tissue area was aseptically separated with a scalpel, weighted, diluted 10 times with sterile peptone water and homogenized for 120 s in a Stomacher instrument (Lab-Blender 400, Seward, London, the UK). After homogenisation, decimal serial dilutions were performed with sterile peptone water. Finally, 0.1 mL was seeded on plate count agar (PCA, 30 • C/24 h) for total aerobic mesophilic microorganisms (TAM) counts and potato dextrose agar acidified to pH 3.5 with a sterilised solution of tartaric acid (PDA, 25 • C/5 days) for yeast and mould counts. The same procedure was performed with sound tissues. The microbial analysis was performed in triplicate by each batch.
In parallel, cherry samples were immersed in 100 mg/L sodium hypochlorite solution for 3 min to sanitise the surfaces of the cherries. Cherries were air-dried, and then microbiological analysis was carried out as above in damaged and sound cherry tissues. As a control, cherries without 'aqueous spot' symptoms were analysed. The analyses were performed in triplicate, and microbial counts were expressed as log 10 CFU/g of cherries.
The isolation of microbes, bacteria and fungi was performed based on the morphology of colonies. Three colonies of each morphology detected from the most diluted plates were isolated by streaking in the specific culture media (PCA and PDA). Once pure colonies were obtained, cultures were stored at −80 • C in 40% sterile glycerol.
An alternative procedure for the isolation of moulds was performed. Three cherries by sample with 'aqueous spot' symptoms were individually stored in sterile plastic bottles at room temperature (~20-25 • C) until mycelia development was evident ( Figure 1E,F). The mould was directly inoculated onto acidified PDA plates with a loop. After incubation at 25 • C/5 days, mould colonies were isolated as above. As controls, cherries without 'aqueous spot' symptoms were also stored in sterile plastic bottles at room temperature.

Identification of Isolated Microbes
Genomic DNA from microbial isolates, bacteria, and fungi was extracted using the NucleoSpin ® Microbial DNA Kit (Macherey-Nagel, Düren, Germany) according to the manufacturer's instructions. The DNA concentration and quality were assessed by a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). To identify the isolates at the species level, for bacteria, the 16S rRNA gene was amplified by PCR using the universal primers 27F and 1492R as described by Ruiz-Moyano et al. [26]. For fungi, the internal transcribed spacer ITS1/ITS2-5.8S rRNA region was amplified as detailed by Gallardo et al. [27] using the primer pairs ITS1 and ITS4 described by White et al. [28]. The resulting amplicons were visualised after electrophoresis on a 1% agarose gel stained with Midori Green Advance (Nippon, Tokyo, Japan). The Generuler 100 bp plus DNA ladder (Thermo-Fisher Scientific, San Jose, CA, USA) was used as the reference. The PCR products were purified using the Genejet PCR Purification Kit (Thermo-Fisher Scientific) and then sequenced at the Service of Bioscience Applied Techniques (STAB) of the University of Extremadura (Badajoz, Spain). The sequences were analysed and edited with Chromas Pro version 1.49 beta (Technelysium, City, QL, Australia) and checked by nucleotide-nucleotide comparison at EMBL/GenBank databases using the BLAST algorithm. The identities of the isolates were determined based on the highest score.
Additionally, identities of the yeast isolates identified by ITS sequencing as Metschnikowia spp. were confirmed by amplification of the D1/D2 domain of the 26S large ribosomal subunit (LSU) of rRNA using the primers NL1 and NL4 [29], following the PCR conditions described by Gallardo et al. [27] and subsequent sequencing of the PCR amplicons as above.

Microscopic Observation of Cherry Tissues
Cherries with initial 'aqueous spot' symptoms ( Figure 1A,B) from all batches were selected for tissue observation. The damaged portion was aseptically obtained and smoothly homogenised. After that, direct observations were performed. Microscopic observations were made at 40× magnification with sub-stage illumination (DMLS, Leica, Buccinasco, MI, Italy). The structures were visualised with a Leica DM 2000 LED microscope.

Inoculation Assays
The effect of microbes isolated from samples with symptoms of disease (Section 2.2) on cherries was assessed by direct inoculation. Undamaged cherries of 'Burlat', 'California' and 'Lapins' cultivars were used. After sanitisation, cherries were injured with a sterile tip (one wound per cherry,~3 mm of diameter × 3 mm of depth).
Pure cultures of bacteria were grown in nutrient broth at 30 • C for 24 h and recovered by centrifugation at 10,000× g for 3 min. Pellets were suspended in sterile water. Pure cultures of yeast were grown in acidified PDA at 25 • C for 48 h. A loop of yeast was suspended in 500 µL of sterile water. Pure mould cultures were grown in acidified PDA at 25 • C for 10 days. Conidia was recovered with a sterile 0.05% solution of tween 80. Microbe solutions were quantified by direct observation using an optical microscope in a Neubauer chamber. As inoculation suspensions, bacteria and yeasts were diluted to 10 7 cells/mL, and conidia was diluted to 10 5 spores/mL.
Ten injured cherries were inoculated with 5 µL of each microbe solution for the individual inoculations. Each cherry was stored individually in a sterile plastic bottle at 25 • C for 10 days.
A selection of bacteria, yeasts and moulds was made based on their frequency of isolation, incidence and symptoms in the assays mentioned above for combinations of double and triple inoculations. Then, combined inocula were carried out with mould × bacteria, mould × yeast and finally mould × bacteria × yeast.
In all inoculation tests, the incidence percentages and rotting symptoms were recorded daily, and pictures were taken periodically. The assays were performed in triplicate. Percentages of incidence were converted into Bliss angular values before analysis.
To know the capacity of the potential agents to infect sound sweet cherries, an additional inoculation assay was carried out by immersion in the selected moulds. To do this, 10 cherries, free of defects, were immersed for 3 min in a solution of conidia at a concentration of 10 5 spores/mL. After air-drying in an aseptic environment, cherries were placed in sterile plastic containers and incubated at 25 • C for five days. The assays were performed in triplicate.

Genetic Characterisation of Botrytis Cinerea Morphotypes
To characterise Botrytis cinerea isolates, two representative isolates from the four morphotypes detected in inoculation assays (Section 2.5) were selected. A multi-locus approach using four phylogenetic markers was used. These were partial fragments of ITS, LSU 26S rDNA, calmodulin and β-tubulin regions. The PCR amplification of the ITS region was carried out using the primers ITS1 and ITS4 [28], LSU 26S rDNA using the primers NL1 and NL4 [29], calmodulin using primers CMD5 and CMD6 [30], and β-tubulin using primers Bet2a and Bet2b [31]. PCR amplification was performed in 50 µL of reaction mixtures containing 10 ng of genomic DNA, 50 pmol of the forward primer and reverse primer, 0.2 mM of each dNTP, 0.1 vol of 10X PCR buffer, and 1.25 U Green DreamTaq DNA polymerase (Fermentas, Thermo Fisher Scientific Inc., Waltham, MA, USA). PCR reactions were run in a T100™ Thermal Cycler (BioRad, Hercules, CA, USA) using an initial denaturation at 94 • C for 4 min, followed by 35 cycles of 94 • C for 1 min, annealing at 55 • C for 1 min, elongation at 72 • C for 1 min, and a final extension at 72 • C for 10 min. PCR amplicons were visualised and sequenced as above.

Analysis of Results
Data were evaluated statistically using SPSS for Windows (version 19.0; SPSS Inc., Chicago, IL, USA). Mean values of microbial counts were compared by one-way analysis of variance (ANOVA). Multiple comparisons of the means were performed using Tukey's honestly significant difference (HSD) test (p ≤ 0.050).
The analyses of TAM and yeast counts per season showed no differences between damaged tissues in both years. Meanwhile, sound tissues from 2020, had significantly (p < 0.050) lower counts than samples from 2019. However, the analyses of total counts from both seasons showed higher counts (p < 0.001) in damaged tissues (7.05 ± 0.98 and 6.39 ± 0.92 log 10 CFU/g of cherries for TAM and yeast, respectively) compared to sound tissues (6.09 ± 1.24 and 5.19 ± 1.22 log 10 CFU/g of cherries for TAM and yeast, respectively). Sanitisation of cherry surfaces reduced TAM and yeast populations by around 1 and 2 log 10 CFU/g, respectively, with higher counts (p < 0.050) for damaged tissues than sound tissues ( Figure 2C). Figure S1 shows in detail the TAM and yeast counts after sanitisation in all cultivars collected in 2019 and 2020 seasons. type strains (CECT2850, CECT20516, CECT20518, and CECT20793). Evolutionary ana yses were conducted using MEGAX [33].
Finally, molecular detection of the transposable elements Boty and Flipper was con ducted as described in Muñoz et al. [35].

Analysis of Results
Data were evaluated statistically using SPSS for Windows (version 19.0; SPSS Inc Chicago, IL, USA). Mean values of microbial counts were compared by one-way analysi of variance (ANOVA). Multiple comparisons of the means were performed using Tukey' honestly significant difference (HSD) test (p ≤ 0.050).
The analyses of TAM and yeast counts per season showed no differences betwee damaged tissues in both years. Meanwhile, sound tissues from 2020, had significantly ( < 0.050) lower counts than samples from 2019. However, the analyses of total counts from both seasons showed higher counts (p < 0.001) in damaged tissues (7.05 ± 0.98 and 6.39   The presence of moulds was inconsistently detected on cherries. Mould colonies appeared on the least diluted plates, and media counts were always lower than the detection limit for microbial analyses (2 log 10 CFU/g). However, microscopic observations of damaged tissues from all samples showed the presence of mycelium within the flesh, as shown in Figure 3A. Additionally, the final stage of this disease in the affected cherries inside the plastic bottles included complete coverage of the fruit with white, occasionally grey, mycelium ( Figure 3B). ( Figure 2C). Figure S1 shows in detail the TAM and yeast counts after sanitisation in all cultivars collected in 2019 and 2020 seasons. The presence of moulds was inconsistently detected on cherries. Mould colonies appeared on the least diluted plates, and media counts were always lower than the detection limit for microbial analyses (2 log10 CFU/g). However, microscopic observations of damaged tissues from all samples showed the presence of mycelium within the flesh, as shown in Figure 3A. Additionally, the final stage of this disease in the affected cherries inside the plastic bottles included complete coverage of the fruit with white, occasionally grey, mycelium ( Figure 3B).

Species Identification of Bacteria, Yeast, and Moulds
The bacteria identified on the cherries were dominated by species belonging to the Enterobacterales order, isolated in all cases and on both samples and seasons. Rahnella aquatilis, Pantoea agglomerans, Tatumella tyseos, Hafnia paralvei, Rosenbergiella nectarea, Erwinia spp., and Lonsdalea quercina were the species identified (Table 1). Regarding Erwinia spp., a more specific identification was not possible because BLAST gave a similar

Species Identification of Bacteria, Yeast, and Moulds
The bacteria identified on the cherries were dominated by species belonging to the Enterobacterales order, isolated in all cases and on both samples and seasons. Rahnella aquatilis, Pantoea agglomerans, Tatumella tyseos, Hafnia paralvei, Rosenbergiella nectarea, Erwinia spp. and Lonsdalea quercina were the species identified (Table 1). Regarding Erwinia spp., a more specific identification was not possible because BLAST gave a similar percentage of identities with different accession numbers in GenBank of non-pathogenic species such as E. billingiae, E. tasmaniensis, and E. persicina. Bacillus thuringiensis and Agrococcus lahaulensis were also identified in 2019, while Leuconostoc mesenteroides and Pseudomonas spp. were identified in 2020. All Enterobacterales species were found in all samples tested, both damaged and sound tissues, sanitised or not, as well as cherries without any disease symptoms. Twelve yeast species were identified in cherry samples affected by 'aqueous spot' (Table 1). Candida oleophila, Candida railesensis, Aureobasidium pullulans, Filobasidium wieringae, Filobasidium magnum, Pichia kluyveri, and Hanseniaspora uvarum were found in both seasons, although they were not present in all samples. Yarrowia lypolitica, isolated from all samples, and Rhodotorula nothofagi were found in 2019; while Metschnikowia pulcherrima and Metschnikowia viticola, isolated from all samples, were found in 2020. The identification percentage of Metschnikowia species of ITS1-5.8S rDNA-ITS4 was below 97%, due to hairpins in the genomic ribosomal sequences [36]. The identification was confirmed by sequencing the 26S rDNA large subunit, with the percentage of identities above 99% (Table 1). Eleven out to twelve species were isolated from sound tissues. Geotrichum candidum was isolated in one 2020 sample from the lowest dilutions. The type of tissue (damaged or sound) and the hygienic conditions did not modify the identified species.
Mould presence was scarce in cherries with initial disease symptoms, as mentioned above. Eleven species were identified, although none were isolated from all samples. Penicillium crustosum, Mucor plumbeus, Alternaria alternata, Cladosporium cladosporioides, and Botrytis cinerea were isolated from samples in both seasons. Schyzophylum commune and Trichodema atroviride were found in samples from 2019, while Fusarium oxysporum, Mucor racemosus, Monilia laxa, and Alternaria infectoria were isolated from samples in 2020.
When the disease advanced, cherries were completely covered by primarily white or grey mycelia. The isolation and identification of this mould consistently confirmed the presence of different morphologies of B. cinerea in the last stages of the disease.

Inoculation of Bacteria, Yeast, and Moulds in Wounded Cherries
The results of the individual microbial inoculations of wounded cherries are shown in Table 1. When the rotting of cherries was observed, the inoculation of bacteria usually generated the browning of tissues and a slight liquefication around the wounds. Erwinia spp. produced bubbles inside of the wounds. Incidence percentages ranged from 0.0% (inoculations with B. cereus and A. lahaulensis) to 83.3% (L. quercina).
Yeast inoculation, in general, did not modify the tissue around wounds. Hanseniaspora uvarum, M. pulcherrima and M. viticola spoiled 13.3%, 20.0%, and 33.3% of cherries, respectively. The observed symptoms were mainly depression of wounds and slight browning of the tissues (Figure 4C,D). In contrast, G. candidum infected 100% of the inoculated wounds, with visual observation of yeast development ( Figure 4E). However, the symptoms were not indicative of 'aqueous spot' evolution.
Most mould species infected all inoculated wounds (100.0% of incidence) with symptoms as shows Figure 4F (F. oxisporum), Figure 4G (P. crustosum), Figure 4H (A. alternata), Figure 4I (S. commune), Figure 4J (M. laxa) and Figure 4K1,K4 (B. cinerea). The other inoculated mould species presented less infection capacity. Mucor plumbeus and T. atroviride presented 0.0% incidence, C. cladosporioides presented 36.7%, and M. racemosus and A. infectoria showed 56.7%. Interestingly, the infection evolution of the inoculated moulds demonstrated that B. cinerea infections ( Figure 4K1) showed similar symptoms to cherries in the last stage of this disease ( Figures 1F and 3B), and other morphotypes less prevalent in the 'aqueous spot' were developed ( Figure 4K2,K4). Neither of these moulds exhibited the ability to infect sound cherries by immersion in the conidia suspension.
Genetic characterisations of the different morphotypes of B. cinerea were carried out, selecting two representative isolates of each morphotype. Four reference strains were included in this work. The phylogenetic tree constructed with the concatenated sequences of ITS1-5.8S rDAN-ITS4 + LSU 26S rDNA+ calmodulin gene + β-tubulin is presented in Figure 5. The representative isolates of morphotype K1 (white mycelium) were clustered separately from the rest of the morphotypes and reference strains. The K1 sequences did not show differences in LSU 26S rDNA sequences with respect to the rest of B. cinerea strains, and only one single nucleotide polymorphism (substitution of A by T) was detected in ITS1-5.8S rDNA-ITS4. On the other hand, significant differences were found between white mycelium morphotypes and the rest in the other analysed genes, with percentages of identities 89.3% to 98.6% for calmodulin sequences and 94.1% to 94.3% for β-tubulin sequences. High genetic differences among K1 morphotypes and the rest of the strains were revealed by analysing ISSR-PCR profiles. Thus, the profiles with (CAG) 4 , (GTG) 5 , and UBC809 were unique and consistent in white mycelium B. cinerea strains.  Genetic characterisations of the different morphotypes of B. cinerea were carried out, selecting two representative isolates of each morphotype. Four reference strains were included in this work. The phylogenetic tree constructed with the concatenated sequences of ITS1-5.8S rDAN-ITS4 + LSU 26S rDNA+ calmodulin gene + β-tubulin is presented in Figure 5. The representative isolates of morphotype K1 (white mycelium) were clustered separately from the rest of the morphotypes and reference strains. The K1 sequences did not show differences in LSU 26S rDNA sequences with respect to the rest of B. cinerea strains, and only one single nucleotide polymorphism (substitution of A by T) was detected in ITS1-5.8S rDNA-ITS4. On the other hand, significant differences were found be- The results of the double and triple inoculations are shown in Table 2. Botrytis cinerea M2 was selected as a base for this assay. The combination with bacteria better represented the 'aqueous spot' symptoms than B. cinerea M2 alone, with 100% of cherries infected. Among them, the combination with L. quercina and L. mesenteroides showed a higher degree of rotting because the browning area was larger.
percentages of identities 89.3% to 98.6% for calmodulin sequences and 94.1% to 94.3% for β-tubulin sequences. High genetic differences among K1 morphotypes and the rest of the strains were revealed by analysing ISSR-PCR profiles. Thus, the profiles with (CAG)4, (GTG)5, and UBC809 were unique and consistent in white mycelium B. cinerea strains. The results of the double and triple inoculations are shown in Table 2. Botrytis cinerea M2 was selected as a base for this assay. The combination with bacteria better represented the 'aqueous spot' symptoms than B. cinerea M2 alone, with 100% of cherries infected. Among them, the combination with L. quercina and L. mesenteroides showed a higher degree of rotting because the browning area was larger.
Co-inoculation of B. cinerea M2 with yeast presented a disorder incidence between 26.7% (combined with M. pulcherrima) and 100% (Y. lipolytica, C. railenensis, H. uvarum, and M. viticola). Nevertheless, the degree of rotting was lower in all combinations, and mycelia development was delayed.
Lastly, triple combinations were tested, using L. quercina and T. terrea as bacteria, and Y. lipolytica, A. pullulans, H. uvarum, and M. viticola based on their ubiquity and symptoms. All combinations developed symptoms associated with this disorder (100% incidence). However, the degree of rotting was lower in all cases than B. cinerea M2 alone. The areas of degraded tissues were smaller, and the appearance of visual mycelia was delayed.  Co-inoculation of B. cinerea M2 with yeast presented a disorder incidence between 26.7% (combined with M. pulcherrima) and 100% (Y. lipolytica, C. railenensis, H. uvarum, and M. viticola). Nevertheless, the degree of rotting was lower in all combinations, and mycelia development was delayed.
Lastly, triple combinations were tested, using L. quercina and T. terrea as bacteria, and Y. lipolytica, A. pullulans, H. uvarum, and M. viticola based on their ubiquity and symptoms. All combinations developed symptoms associated with this disorder (100% incidence). However, the degree of rotting was lower in all cases than B. cinerea M2 alone. The areas of degraded tissues were smaller, and the appearance of visual mycelia was delayed.

Discussion
The emerging disease called 'aqueous spot' is a serious issue for cherry producers in the Jerte Valley (Cáceres, Spain). The initial symptoms begin before harvest but are often detected during the commercialisation of cherries. The events or agents that trigger this disease have remained unclear until now. A common feature in fruits deeply affected by 'aqueous spot' is the external development of a characteristic mycelium, which is primarily white and eventually covers the cherries completely. Mould counts in the early stages of 'aqueous spot' were beneath the detection thresholds, with inconsistent presence of various well-known rotting species. Monilia laxa [9], A. alternata [37], P. crustosum [11], F. oxysporum [38] and B. cinerea [13] have been associated with different diseases of sweet cherries. Inoculation trials showed a high disease incidence in most cases. White mycelial were observed in the wounded cherries inoculated with F. oxysporum ( Figure 4F), A. alternata ( Figure 4H), S. commune ( Figure 4G), and B. cinerea ( Figure 4K1). However, the evolution and morphology of white B. cinerea was uniquely compatible with 'aqueous spot' symptoms. The consistent isolation of white B. cinerea at the last stages of the disease confirmed its association with 'aqueous spot'. Mycelium presence within the flesh of cherries at the early stages of 'aqueous spot' development, even in surface-sanitised fruits, placed this mould at the origin of symptoms. The mode of infection of this morphotype of B. cinerea remains unclear; however, growth failure with immersion inoculation demonstrated that wounds or latent infections [39] may be involved in its presence and development. The analysis of the levels of expression of some virulence genes during the infection could help to understand it. The fact that early cherry cultivars, such as 'Burlat', and a few late cultivars, such as 'Lapins' and 'Sweetheart', are prone to suffer 'aqueous spot' will need to be studied in future works. Previous studies indicated that cherry cultivars are susceptible to a wide range of pathogens [40], which may be influenced by canopy structure, crop load and blooming times, as well as the susceptibility to bounds, cracking, or mechanical damage due to genotypic differences [14,41].
In this work, Botrytis spp. isolates associated with 'aqueous spot' pathology displayed a wide morphological ( Figure 4K) and genetic variability ( Figure 5). Previous works have demonstrated the extensive phenotypic and genetic diversity of B. cinerea from different origins [35,42]; however, the reports of white mycelium Botrytis are not common. It has been postulated the division of B. cinerea into two sympatric species, Vacuma and Transposa, defined by the presence of the transposable elements Boty and Flipper [43], which are related to its pathogeny and resistance to fungicides. Morphotypes K1, K2, K4, and CECT 2850 were classified as Boty, whereas the K3 morphotype was classified as Vacuma (did not contain any transposable elements) ( Figure 5). The strains of another type were classified as Transposa, which had both transposable elements.
Bacteria and yeast populations were modified by the development of symptoms of this disease. Microbial counts (TAM and yeasts) of damaged cherries increased around one log 10 CFU/g of cherries compared to sound tissues. In general, mean counts, even in sound tissues, were higher than those reported by Serradilla et al. [4]. Sanitising cherries with sodium hypochlorite reduced counts [44], but they were still higher than previously reported. This fact could be associated with increased spoilage of cherries in 2020. Wide variability in the microbial dominance among samples was detected, although the mean counts indicated that bacteria counts in damaged tissues were higher (around 0.5 log 10 CFU/g, p < 0.050). Members of the order Enterobacterales dominated the bacterial population of damaged and sound tissues in all samples from both seasons. Rahnella spp., Erwinia spp. and R. nectarea have been reported as an epiphytic population of phyllosphere from different species of fruit trees [4,45,46]. The presence of T. terrea could be associated with vectors, such as Drosophila suzukii [47]. This group of bacteria are usually not related to fruit rot or spoilage, beyond a recent report on P. agglomerans [48]. Overall, the bacteria isolated from damaged cherries had a limited capacity for spoiling tissues around the wounds in individual inoculations. Lonsdalea quercina showed the highest occurrence ( Figure 4A), producing fluid and browning of tissues around the wound. References have associated this species with canker disease in poplar trees [49]. Some bacteria, such as Erwinia spp., produced bubbles in the inoculated area, a symptom not associated with 'aqueous spot'. The combination of bacteria and B. cinerea M2 showed similar incidence and symptoms to 'aqueous spot' than M2 inoculated alone. However, the combinations with L. quercina and L. mesenterioides showed more rotting tissues around wounds. Bacteria, by itself, showed low pathogenic capacity. Therefore, the results revealed that events that originate 'aqueous spot' symptoms favoured bacteria multiplication, and the symptoms could be enhanced.
The yeast population was increased by 'aqueous spot' events, as previously mentioned. The species identified commonly inhabit fruits' surfaces. Specifically, A. pullulans, H. uvarum, P. kluyveri and M. pulcherrima have been identified in cherries [4,50]. Differences between yeast species were observed between seasons. Y. lipolytica was isolated in all samples in 2019, while M. pulcherrima and M. viticola appeared in all samples in 2020. Overall, the results of individual inoculations showed less incidence of yeasts, except for G. candidum (100% incidence), previously associated with sour rot in stone fruit [20]. The involvement of yeasts at the initial stages of 'aqueous spot' would be unlikely according to the results of the combined inoculations. The inhibition of B. cinerea M2 was detected in all combinations, decreasing the incidence in double inoculations, and delaying the appearance of visual mycelia in double and triple inoculations. The colonisation of wounds and competition for space and nutrients could explain this finding. Additionally, several strains of this yeast species have shown antagonistic activities, such as M. pulcherrima, H. uvarum [50][51][52], C. oleophila [53], and A. pullulans [54].
The present work describes for the first time the microbial population involved in the development of 'aqueous spot', an emerging disease of sweet cherries from the Jerte Valley. Genetic identifications and inoculation trials indicated that white mycelium B. cinerea were found at the early stages of this disease. This morphotype showed a wide genetic divergence with other B. cinerea strains, remarkably in calmodulin and β-tubulin gene sequences and ISSR-PCR profiles. Members of the order Enterobacterales showed a synergistic effect with B. cinerea, increasing the degree of rotting, while yeasts, at the initial stage, had an inhibitory effect on symptom development. Although this work demonstrates the involvement of white mycelium B. cinerea in this disease, further research using molecular biology techniques such as qRT-PCR or next-generation sequencing is needed to fully understand the origin and development of 'aqueous spot' in order to develop strategies for its control.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/foods10102281/s1, Figure S1: Counts of total aerobic mesophilic and yeast in damaged and sound tissue from all sweet cheery cultivars after sanitisation. Funding: This study was funded by the Spanish Ministry and Innovation (MCI), Spanish State Research Agency (AEI), and the European Regional Development Fund (ERDF, EU) under the grant RTI2018-096882-BI00 and by Junta de Extremadura Government and ERDF under grants GR18165 and GR18196.
Data Availability Statement: Data sharing does not apply to this article.