Fungal Microbiota of Sea Buckthorn Berries at Two Ripening Stages and Volatile Profiling of Potential Biocontrol Yeasts

Sea buckthorn, Hippophae rhamnoides L., has considerable potential for landscape reclamation, food, medicinal, and cosmetics industries. In this study, we analyzed fungal microorganism populations associated with carposphere of sea buckthorn harvested in Lithuania. An amplicon metagenomic approach based on the ITS2 region of fungal rDNA was used to reveal the ripening-affected fungal community alterations on sea buckthorn berries. According to alpha and beta diversity analyses, depending on the ripening stage, sea buckthorn displayed significantly different fungal communities. Unripe berries were shown to be prevalent by Aureobasidium, Taphrina, and Cladosporium, while ripe berries were dominated by Aureobasidium and Metschnikowia. The selected yeast strains from unripe and mature berries were applied for volatile organic compounds identification by gas chromatography and mass spectrometry techniques. It was demonstrated that the patterns of volatiles of four yeast species tested were distinct from each other. The current study for the first time revealed the alterations of fungal microorganism communities colonizing the surface of sea buckthorn berries at different ripening stages. The novel information on specific volatile profiles of cultivable sea buckthorn-associated yeasts with a potential role in biocontrol is important for the development of the strategies for plant cultivation and disease management, as well as for the improvement of the quality and preservation of the postharvest berries. Management of the fungal microorganisms present on the surface of berries might be a powerful instrument for control of phytopathogenic and potentially antagonistic microorganisms affecting development and quality of the berries.


Introduction
Sea buckthorn (Hippophae rhamnoides L.) is garden and landscaping shrub widespread in Europe, Asia, and North America, distinguished for good adaptation to barren soils and harsh climates [1]. The ecological impact of sea buckthorn in preventing soil erosion, improving its properties, and reducing pollution has been demonstrated [2]. Economically, the plant species has been considered as a source of healthy nutritious food, giving rise to an increased demand for medicine and cosmetics [3,4]. developing natural tools for disease management, postharvest protection, and the quality of the sea buckthorn berries with potential in food production. The data on volatile profiling of selected yeasts could be relevant for improving plant productivity and generating new environmentally friendly formulations for pest control and plant disease management.

Ethics Statement
The collection of sea buckthorn berries was authorized by the private owner, fully acknowledged in a paper. Endangered or protected species were not involved; therefore, no specific permissions were required.

Sample Collection
Hippophae rhamnoides L. berries were aseptically sampled from the organic private plantation located in the Molėtai region of Lithuania (GPS coordinates: 55 • 15 12.2"N, 25 • 26 23.1"E). The sea buckthorn berries were collected from three different places of sea buckthorn plantation separated by more than 50 meters. The berries were randomly selected from three bushes and combined into one biological replicate at each location. The experiment was carried out at two different growth stages: at the beginning of berry ripening (UB) in July and ripe berries (RB) in September 2018. The samples were collected into sterile plastic bags, transported to the laboratory and processed within 2-3 h after harvesting.
For metagenomic analysis, 300 g of the berries were placed into flasks with 500 mL of sterile 0.05M phosphate buffer (pH 6.8) for 30 min and incubated at room temperature with shaking at 120 rpm. Outwashes were filtered through 420 µm filters and followed centrifugation at 12,000× g for 20 min. The pellets were stored at −20 • C for subsequent DNA extraction.

Metagenomic Analysis of Sea Buckthorn Fungal Microbiota
For metagenomic analysis, DNA was isolated from collected sediments (40 mg of pellet per sample) using a Genomic DNA purification kit (Thermo Fisher Scientific Baltics, Vilnius, Lithuania) and following the manufacturer's protocol. The quantity and quality of extracted DNA were measured using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific). DNA samples of fungal microorganism populations were amplified using the ITS2 region-specific primers: ITS3-KYO2 (5 -GATGAAGAACGYAGYRAA-3 ) and ITS4 (5 -TCCTCCGCTTATTGATATGC-3 ) [36]. Amplicon libraries were prepared using modified Illumina adapters (www.illumina.com), validated on an Agilent Technologies Bioanalyzer DNA 1000 and paired-end sequenced (2 × 300 bp) using Illumina MiSeq platform (BaseClear B.V. Leiden Netherlands). Complete data sets were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database (Accession number PRJNA590349).

Bioinformatics Analysis
Sequences were pre-processed, quality filtered, and analyzed using QIIME2 version 2018.4 [37]. Forward and reverse reads were trimmed with Cutadapt 2.8 [38]. Primer-free fastq files were processed using Divisive Amplicon Denoising Algorithm 2 (DADA2) plugin implemented in QIIME2. Reads with higher than 2.6 expected error rates were discarded. Chimeric sequences were eliminated using 'consensus' method. Alpha and beta-diversity analyses were performed with the q2-diversity plugin in QIIME2 at a sampling depth of 43458. Principal coordinates analysis (PCoA) was performed based on non-phylogenetic Bray-Curtis distances and visualized with the make_2d_plots.py script. To assign taxonomy of obtained Amplicon sequence variants (ASV) we used the QIIME2 q2-feature-classifier plugin and the v7.2 UNITE database as the reference at 97% identity [39].

Cultivable Yeast Identification
Yeast strains were streak-plated from single colonies onto YPD-agar medium supplemented with 50 µg mL -1 chloramphenicol and incubated for 3 days at 25 • C. The colors and the structure of the colonies were recorded. For the microscopy study, slides were prepared from 1-week-old cultures and observed under the light microscope (Leica DM750, Wetzlar, Germany), by recording micrographs with the digital camera (Leica ICC50 HD, Wetzlar, Germany).
DNA was isolated from fresh yeast culture (24 h) by using the Genomic DNA purification kit (Thermo Fisher Scientific Baltics, Vilnius, Lithuania) following the manufacturer's instructions. For identification of yeast, PCR amplification of the region between the 18S rRNA and 28S rRNA genes was performed using ITS1 (5 -TCCGTAGGTGAACCTGCGG-3 ) and ITS4 (5 -TCCTCCGCTTATTGATATGC-3 ) primers. The reaction mixtures (50 µL) contained 5 µL of DreamTaq green buffer, 1 µL of 2 mM dNTP mix, 1 µL of each primer (10 µmol L -1 ), 2.5 unit of DreamTaq DNA polymerase (all from Thermo Fisher Scientific Baltics, Vilnius, Lithuania) and 1 µL of DNA template (5 ng). The PCR conditions were as follows: initial denaturation at 94 • C for 5 min; 25 cycles of 94 • C for 1 min, 53 • C for 1 min 30 s, 72 • C for 2 min; a final elongation at 72 • C for 10 min. For preliminary yeast species identification, PCR products were digested with CfoI and Hinf I enzymes and restriction profiles were checked by 1% agarose gel electrophoresis. PCR products were purified using a GeneJet PCR purification kit (Thermo Fisher Scientific Baltics, Vilnius, Lithuania), following the manufacturer's instructions and sequenced using ITS1 and/or ITS4 primers at BaseClear (Leiden, Netherlands). The generated sequences were compared with those found in the FASTA network service of the EMBL-EBI database (http://www.ebi.ac.uk/Tools/sss/fasta/nucleotide.html) and deposited in the National Center for Biotechnology Information (NCBI) under accession numbers MN700626, MN700627, MN700628, and MN700629.

Sampling and Identification of Volatile Organic Compounds (VOCs) Produced by Yeasts
The solid-phase micro-extraction (SPME) technique was applied to collect the headspace volatiles of yeasts. Yeasts were cultivated in polystyrene Petri dishes (Ø 55 mm × 14 mm) on YPD-agar (14 mL) for 2 days at 25 • C. For sampling background volatiles, YPD-agar plates without yeast were used as control samples. SPME fiber coated with a polydimethylsiloxane-divinylbenzene polymer (DVB/PDMS, 65 mm coating layer thickness, Supelco, Pennsylvania, USA) was used. Before each collection, purification of SPME fibers was conducted at 240 • C for about 10 min in a GC injector. The needle of SPME syringe was placed above the yeast layer through a small hole made in a Petri dish; the fiber was exposed to the headspace for 60 min at room temperature. After sampling, the fiber was inserted for 2 min into the injection port of gas chromatograph for desorption of the volatiles.
The volatiles collected were analyzed using Shimadzu gas chromatograph GC-2010 coupled with Shimadzu mass selective detector MS-QP 2010 Plus (Kyoto, Japan). In the GC, Restek Stabil-Wax column (30 m × 0.25 mm × 0.25 µm, Bellefonte, PA, USA) was used. The oven temperature was programmed as follows: the initial temperature was maintained isothermally at 40 • C for 1 min; afterward, it was raised to 200 • C at a rate of 5 • C min-1, then increased to 240 • C at a rate of 10 • C min-1, and maintained isothermally for 11 min. The GC injector temperature was set at 240 • C. Helium was used as carrier gas (1.5 mL min -1 ). The relative amount of each of the compounds was counted based on the area of the chromatographic peak. The volatile compounds were identified according to their mass spectra and their retention indexes in a NIST version 2.0 mass spectral library (National Institute of Standards and Technology, USA) as well as co-chromatography of synthetic standards. Software GCMS solution version 2.71 (Shimadzu, Kyoto, Japan) was applied for data analysis.

Statistical Data Analysis
Nonparametric Mann-Whitney U test was applied to evaluate differences of volatile amounts between yeast and control samples, i.e., YPD-agar plates without yeast, using the Statistica 6.0 program package (StatSoft, Inc., Tulsa, OK, USA). To assess and visualize the associations between odor blends of four yeast species and volatile compounds, the principal component analysis (PCA) was performed using Canoco 4.5 software (Biometris Plant Research International, Wageningen, The Netherlands). Amounts of styrene, 2,5-dimethyl pyrazine, 2-ethylhexanol and methoxy-phenyl-oxime derived from yeast samples did not differ from those collected from the control; hence, they were not included in PCA analysis. Significantly smaller amounts of the compounds collected from the control samples were subtracted from the amounts released by yeasts. By Canoco 4.5 software, absolute amounts expressed as area under chromatographic peak were log-transformed, scaled dividing each value by its standard deviation and centered.

Diversity and Richness of Fungal Microorganism Communities
To characterize the diversity of fungal microorganism community distributed on sea buckthorn at different ripening stages, we harvested berries from the three distinct shrubs at private organic plantation and generated a combined sample as one biological replicate. For metagenomics analysis, three biological replicate samples of unripe berries (UB) obtained from different locations of the sea buckthorn plantation and three corresponding to ripe berries (RB) were used. Sequencing of the ITS2 rRNA partial gene region of DNA samples generated a total of 726,634 reads with an average of 121,105 sequences per sample (range 897,20-164,922) ( Table 1). Following quality filtering, 454,538 sequences were recovered. The clustering of the sequences generated a total of 1615 unique amplicon sequence variants (ASVs) (1014 [338±70, hereafter median for 3 samples ± standard deviation] for fungal ITS2 from unripe berries and 601 [200±35.6] ASV for ripe berries). The observed number of ASV reveals a significant difference in fungal community richness between UB and RB sample groups (ANOVA p=0.038). The microbial community analysis according to the nonparametric Kruskal-Wallis test showed a statistically significant decrease (p-value 0.046) in fungal microorganism diversity amongst RB samples in comparison to UB samples. In agreement with ASV data, Shannon's diversity and Simpson indexes estimate that unripe sea buckthorn berries have a higher fungal microorganism diversity than ripe berries. Rarefaction curve approached plateaus indicating that the sequencing depth was sufficient to capture the microbial diversity of each sample ( Figure S1). Principal coordinate analysis (PCoA) accomplished with the representative ASVs demonstrated a clear separation of UB and RB samples, indicating a difference in the composition of the fungal microbiota ( Figure 1). Clustering of triplicate sample sequence libraries using the Bray-Curtis dissimilarity index revealed consistency of microbial community structure among replicates.

Composition of Fungal Microbiota on Unripe and Ripe Sea Buckthorn Berries
A total of four phyla, 69 orders, and 144 families were identified in the present study (Table S1). Based on the average abundance analysis, the dominating phyla on both sea buckthorn samples were Ascomycota (85.0% for UB and 97.6% for RB, respectively) followed by the Basidiomycota (10.6% and 1.7% in UB and RB, respectively) and other fungi ( Figure 2A; Table S1). The impact of ripening on the fungal microorganism community was notable starting at the class level. Dothideomycetes were present in large numbers in both samples (54.0% for UB and 44.6% for RB, respectively) ( Figure 2B). However, on unripe berries (UB) the most abundant classes were Taphrinomycetes (18.4%) and Tremellomycetes (7.6%) when in RB samples Saccharomycetes (47.6%) prevailed. At the family level, in UB samples we found dominating Aureobasidiaceae and Taphrinaceae, while RB samples were prevalent by Aureobasidiaceae and Metschnikowiaceae ( Figure 2C). In total, 196 genera were identified in this study; the most abundant were Aureobasidium (31.4%), Taphrina (18.4%), and Cladosporium (8.1%) for UB samples, and Metschnikowia (44.5%) with Aureobasidium (40.6%) for RB samples ( Figure 2D). level. The taxonomic groups comprising less than 1% of the total composition were assigned to "Other".
The Venn diagram presents the distribution of unique fungal ASVs in UB and RB samples ( Figure 3). A total of 887 unique ASVs found in this study, 198 were shared by both UB and RB samples, while 493 ASVs were exclusive to UB and 196 to RB samples, respectively. A heatmap diagram illustrates the distribution of the most abundant fungal ASVs (Figure 4). Hierarchical cluster analysis showed that the microbial community structure displays different patterns between UB and RB groups, while samples collected within the same ripening stage feature similar community structure. The core microbiomes detected in both samples comprised of ASV1 and ASV3 assigned to Aureobasidium pullulans ( Figure 4; Table S2). At the beginning of ripening, besides the Aureobasidium pullulans the most abundant ASVs were ASV4 (7.6%) and ASV14 (2.4%) (both Taphrina carpini), ASV10 (4.7%) (Taphrina), ASV6 (6.8%) (Cladosporium), and ASV20 (1.6%) (Filobasidium wieringae). During the ripening of the berries, their prevalence diminished to less than 1%. A dramatic change was observed between the beginning and the end of the ripening period, mainly due to the increase of fungal microorganisms belonging to Metschnikowia spp. (ASVs 2, 5, 7, 9, 11, 22, and 30). The overall abundance of ASVs assigned to Metschnikowia in RB berries reached 44%, while on the unripe berries they consisted only 0.7%.

Sea-Buckthorn-Associated Cultivable Potential Biocontrol Yeasts and Their Produced Volatiles
By applying cultivation techniques, more than 20 yeast strains were isolated from the surface of sea buckthorn berries. Yeasts were subjected for primary species identification by morphological colony analysis and bright field microscopy with following fingerprinting of PCR-amplified sequences of internal transcribed spacers 1 and 2 including 5.8S ribosomal RNA gene. For further analysis of VOC production, selected strains proceeded sequencing-based species level confirmation and deposition of sequences in NCBI.
The gas chromatographic-mass spectrometric analysis revealed 29 compounds that were exclusively present in the headspace of four yeast species and occurred significantly more abundant compared to those of control samples. The yeast released volatiles were composed of 11 esters, 6 alcohols, 5 volatile fatty acids, 3 ketones, 2 aromatics, and 2 compounds bearing both ester and aromatic moieties ( Table 2). The principal component analysis revealed that the volatile blend of any of the four yeast species was clearly separated from each other. The first principal component axis accounted for 56% of the total variation and to a large extent explained separation between volatile blends of C. wieringae versus H. uvarum yeasts and A. pullulans versus M. pulcherrima yeasts ( Figure 6A). Volatile profiles of C. wieringae yeasts were characterized by five alcohols (9, 11, 14, 15, and 23), three volatile fatty acids (25, 26, and 28), and two ketones (12 and 20) which were exclusively present or occurred in the large amounts compared to those released by H. uvarum yeasts ( Figure 6B, Table 2). Ten esters were unique for the volatile blends of H. uvarum yeasts in addition to acetic acid (22), 2-hydroxy-3-butanone (18), and 2-phenylethanol (33) with amounts larger compared to those of C. wieringae yeasts ( Figure 6B, Table 2). M. pulcherrima yeasts released the larger amounts of esters (1, 3, 4, and 10), while A. pullulans yeasts released the higher amounts of two alcohols (ethanol (2), butanol (11)) ( Table 2).  Table 2.

Discussion
The current study, for the first time, uncovered the structure of fungal microorganism communities colonizing the surface of sea buckthorn berries at different ripening stages in natural habitat. For a comprehensive analysis, we have applied next generation sequencing and observed the wide diversity of eukaryotic microorganisms on tested berries. The total abundance of fungal microorganisms changed significantly during the period of fruit ripening. The microbial alpha-diversity-estimated by the Simpson's index, Shannon's diversity index, and Pielou's evenness index-was highest among the fungal community on unripe berries. The analysis of the fungal microbiota structure revealed clear differentiation between unripe and mature berries, as demonstrated by the Bray-Curtis beta-diversity analysis. According to our data, unripe sea buckthorn berries were dominated by Aureobasidium, Taphrina, Cladosporium, and Filobasidium (syn. Cryptococcus) genera microorganisms, while on ripe berries prevailed Aureobasidium and Metschnikowia (Table S1, Figure 2). It was previously demonstrated that Cryptococcus, Rhodotorula, Sporobolomyces, and Aureobasidium pullulans yeasts were associated with the early stage of maturation of fruits and berries, such as apples, grapes, nectarines, etc. [29,32,42]. In our previous work, it was demonstrated that Dothioraceae family fungal microorganisms prevailed on sea buckthorn berries, which at the lower taxonomic level were attributed to uncultured Aureobasidium [33]. Our results were partially consistent with previous findings, since most fungal microorganisms present on sea buckthorn were represented by Aureobasidium. However, the overall structure of fungal microbiota varies most likely depending on the differences in climatic conditions, berry harvesting period (previously-September 2016; in this study-July 2018), localization of sea buckthorn plantation, etc.
Various microorganisms act as biological control agents of phytopathogens and have a role in postharvest decay and toxin contamination by emitting volatile organic compounds [40,43]. Several biocontrol yeasts (e.g., Candida oleophila, Aureobasidium pullulans, Metschnikowia fructicola, Cryptococcus albidus, Saccharomyces cerevisiae) have already been used for the production of commercial plant protection products, encouraging the deeper exploration of the fascinating traits and application possibilities of other yeasts as well [16,44,45]. Even the different strains of the same yeast species are able to produce highly different amounts of the same VOCs (e.g., A. pullulans [43]), thus it is important to investigate volatiles of yeasts strains within the species that are already proven to possess biocontrol traits (like Aureobasidium pullulans [43,46], Metschnikowia pulcherrima [46], and Hanseniaspora uvarum [47]), or are less investigated, but closely related to bioactive yeasts (e.g., potential biocontrol yeast Cryptococcus wieringae which belongs to the same genus as C. albidus, C. laurentii, and C. flavus [16,[48][49][50]).
Aureobasidium pullulans is well-known yeast-like fungus hosting on plants, mostly present in a cooler climate like that of central Europe and less related to grapes grown in warm areas [51]. It was revealed that some Aureobasidium species demonstrated biological control activity against phytopathogens, such as Botrytis, Bacillus, Colletotrichum, Penicillium [43,52,53]. The antagonistic activity of A. pullulans mainly includes nutrient competition and production of volatile organic compounds or glucanases, chitinases, and extracellular proteases [40,46]. Taphrina and Cladosporium were the next numerous fungal microorganisms observed on unripe sea buckthorn berries in our study (Table S1, Figure 2). Taphrina genus fungi are known as biotrophic plant pathogens that cause plant and fruit disease symptoms, which may be ascribed to the ability to synthesize and modulate different plant hormones [54]. Cladosporium is considered as ubiquitous fungi, with some species providing numerous antifungal agents [55], while others causing plant or human diseases [56,57]. Cryptococcus encloses species producing biocontrol agents against many pathogens [16,50], while certain species can cause infections in humans and animals [58].
During the ripening of different fruits and berries, the oxidative or weakly fermentative yeast populations-such as Candida, Hanseniaspora, Metschnikowia, and Pichia-increased in frequency [42]. Hanseniaspora spp. have been frequently found on the surface of different fruits, e.g. apples, grapes, or strawberries [59,60]. Numerous yeasts belonging to the Hanseniaspora genus display antagonistic properties against the development of fruit spoilage causing molds [40,61]. Metschnikowia spp. isolated from fruit are sugar-loving yeasts and their populations are most likely stimulated by juices diffusing from the damaged surface of mature fruit [42]. This yeast genus includes several biocontrol agents, which are very effective against fruit spoilage fungi [62]. The biocontrol abilities of M. pulcherrima have been mainly assigned either to competition for nutrients or to the production of volatile organic compounds [15,46]. The competition of M. pulcherrima for iron was reported to play a substantial role in biocontrol interactions [63]. Metschnikowia spp. are linked to several groups of nectarivorous insects, suggesting that this specificity could be correlated to the services provided by the yeast to the insect or vice versa [64,65].
Yeasts produce chemo-diverse blends of volatiles comprised of aldehydes, ketones, alcohols, volatile fatty acids, esters, aromatics, heterocyclic compounds, hydrocarbons, terpenoids, and some others [66,67]. Ethanol is probably the most known volatile metabolite of yeasts derived from carbohydrate catabolism and fermentation [66]; it is one of the major volatiles collected from the headspace of H. uvarum, M. pulcherrima, and C. wieringae yeasts. Amounts of ethyl acetate, another metabolite of carbohydrate catabolism, prevailed in the VOCs produced exclusively by H. uvarum and M. pulcherrima yeasts. It is a valuable solvent for extracting compounds with antimicrobial properties [68] and possess minor inhibitory effect for Botrytis cinerea [69]. In addition to ethyl acetate, 10 esters were unique for the volatile blends of H. uvarum yeasts; these compounds originated from amino acid synthesis and degradation pathways [66]. Prevalence of esters in volatile blends of H. uvarum has been reported in previous studies [70,71]; however, ethanol being the third largest component of the blend in our experiments, has not been determined by Babcock et al. [71] or Piper et al. [70]. High production of ethanol, 3-methylbutanol, 2-phenylethanol, ethyl acetate, and 3-methylbutyl acetate were reported for M. pulcherrima and A. pullulans yeasts [15]. The data published are comparable with the results of our study in terms of large production of alcohols by both yeast species; however, we were not able to detect ethyl acetate in the headspace of A. pullulans. To date, there is no report on VOCs of C. wieringae yeast. Volatile organic compounds produced by A. pullulans that were found in our study (namely 2-phenylethanol, 2-methylpropanol, 2-methylbutanol, and 3-methylbutanol) had been previously shown to control postharvest fruit pathogens (B. cinerea, Colletotrichum acutatum, Penicillium expansum, Penicillium digitatum, and Penicillium italicum) [43]. All these alcohols were produced by A. pullulans SB-18-2, M. pulcherrima SB-18-34, C. wieringae SB-18-25, and H. uvarum SB-18-31 but in different amounts, suggesting the potential biocontrol capabilities of these strains. 2-Phenylethanol is one of the main VOCs produced by various yeast species (e.g. Cyberlindnera jadinii, Lachancea thermotolerans, Candida intermedia, Candida friedrichii, Saccharomyces cerevisiae, Wickerhamomyces anomalus) affecting numerous pathogens (e.g. Aspergillus carbonarius, Aspergillus flavus, Aspergillus ochraceus, Penicillium digitatum, Phyllosticta citricarpa) (see [40] and references therein). Among tested strains, H. uvarum SB-18-31 and M. pulcherrima SB-18-34 produced the highest amount of 2-phenylethanol, whereas A. pullulans SB-18-2 excelled at ethanol generation, while volatilome of C. wieringae SB-18-25 was dominated by 3-methylbutanol. Ethyl acetate was the dominating volatile of H. uvarum SB-18-31 and M. pulcherrima SB-18-34. 2-Phenylethyl acetate was produced in significant amounts (comparing to other tested yeast strains) by H. uvarum SB-18-31 and had been previously reported to contribute to growth inhibition of A. ochraceus growth and ochratoxin A production during processing of coffee [47].
Due to the extremely sensitive olfaction system, insects are capable to use microbial volatiles to search for suitable feeding and oviposition sites [18,72]. Emissions from microorganisms associated with different development stages of fruits and berries contribute to the information about habitat suitability [73] and it was demonstrated that insects recognized yeast communities based on their specific volatile profiles [74][75][76]. Most of the volatile compounds (2-phenyl ethanol, 2-methylbutanol, 3-methylbutanol, ethyl acetate, 2-phenylethyl acetate, ethyl propionate, and others) produced by yeast in our study had been suggested to play a role in modulating behavior of Rhagoletis batava flies [67]. Examples of microbial VOC application for pest control already demonstrated their high potential for use in integrated pest management programs [77]. Identification of attractive and repellent odors from the yeasts associated with berries at preferred and unsuitable stages for pest insect feeding and oviposition would provide a base for the development of environment-friendly pull-push pest control for Rhagoletis batava flies, a major pest of the sea buckthorn.

Conclusions
The metagenomic approach based on the sequencing of ITS2 region of fungal rDNA revealed changes in the fungal community associated with the ripening of sea buckthorn berries. Investigation of alpha and beta diversity of the fungal community, analysis of representative microbial ASVs showed a clear separation among inhabitants of unripe (UB) and ripe berries (RB). Based on results of metagenomics analysis, unripe berry samples prevalently had Aureobasidium, Taphrina, and Cladosporium genera yeasts; while ripe berries were dominated by Aureobasidium and Metschnikowia. Potential biocontrol yeasts Aureobasidium pullulans and Cryptococcus wieringae were isolated from sea buckthorn surface at the early harvesting stage, while Metschnikowia pulcherrima and Hanseniaspora uvarum were amplified from ripe berries and applied for the analysis of produced volatile organic compounds. Eleven esters, six alcohols, five volatile fatty acids, three ketones, two aromatics, and two compounds bearing both ester and aromatic moieties were isolated from the headspaces of four yeast species and identified by gas chromatography and mass spectrometry techniques. Principal component analysis revealed that volatile blends of all four yeast species were clearly separated from each other.
Supplementary Materials: The following are available online at http://www.mdpi.com/2076-2607/8/3/456/s1, Figure S1: Rarefaction curves at a genetic distance of 3% for each sample, Table S1: Fungal taxonomy abundance count for sea buckthorn samples at different ripening stages from phylum to species level, Table S2: Relative abundance of the most common ASVs.